Nous nous intéresserons dans cette thèse à l'étude et la caractérisation des Textures Dynamiques (TDs), avec comme application visée l'indexation dans de grandes bases de vidéos. Ce thème de recherche étant émergent, nous proposons une définition des TDs, une taxonomie de celles-ci, ainsi qu'un état de l'art. La classe de TD la plus représentative est décrite par un modèle formel qui considère les TDs comme la superposition d'ondes porteuses et de phénomènes locaux. La construction d'outils d'analyse spatio-temporelle adaptés aux TDs est notre principale contribution. D'une part, nous montrons que la transformée en curvelets 2D+T est pertinente pour la représentation de l'onde porteuse. D'autre part, dans un objectif de décomposition des séquences vidéos, nous proposons d'utiliser l'approche par Analyse en Composantes Morphologiques. Notre contribution consiste en l'apport et l'étude de nouvelles stratégies de seuillage. Ces méthodes sont testées sur plusieurs applications: segmentation spatio-temporelle, décomposition de TDs, estimation du mouvement global d'une TD, ... Nous avons de plus montré que l'Analyse en Composantes Morphologiques et les approches multi-échelles donnent des résultats significatifs pour la recherche par le contenu et l'indexation de Textures Dynamiques de la base de données DynTex. Cette thèse constitue ainsi un premier pas vers l'indexation automatique de textures dynamiques dans des séquences d'images, et ouvre la voie à de nombreux développements sur ce sujet nouveau. Enfin, le caractère générique des approches proposées permet d'envisager leurs applications dans un cadre plus large mettant en jeu par exemple des données 3D. / This report is focused on the study and the characterization of Dynamic Textures (DTs), with the aim of video indexing in large databases. This research topic being new and emerging, we propose a taxonomy, a definition of DTs and a state of the art. The most representative DT class is described by a model that considers DTs as the superposition of several wavefronts and local oscillating phenomena. The design of spatio-temporal analysis tools adapted to DT is our main contribution. We first show that the 2D+T curvelet transform is relevant for representing wavefronts. In order to analyse and better understand the DTs, we propose in a second step to adapt the Morphological Component Analysis approach using new thresholding strategies. These methods are tested on several applications: decomposition of DTs, spatio-temporal segmentation, global motion estimation of a DT, ... We have shown that Morphological Component Analysis and multi-scale approaches enable significant results for content-based retrieval applications and dynamic texture indexing on the DynTex database. This thesis constitutes a first step towards automatic indexing of DTs in image sequences and opens the way for many new developments in this topic. Moreover, the proposed approaches are generic and could be applied in a broader context, for instance the processing of 3D data.

Motivé par l'étude de tissus biologiques, ce travail contribue aux développements d'outils numériques permettant de prédire la réponse mécanique de matériaux hétérogènes non linéaires dans lesquels les énergies d'interfaces deviennent prépondérantes. Ainsi, une méthode d'homogénéisation multi échelle combinée à une technique de réduction de modèle basée sur la décomposition orthogonale aux valeurs propres est proposée dans un cadre thermique et hyperélastique. Les énergies d'interfaces entre les différentes phases des composites sont décrites par un modèle d'interface cohérent et prises en compte numériquement par une approche liant la méthode des éléments finis étendus et la méthode level-set. Une étude de l'étalement d'une cellule vivante entre deux lamelles fixes est ensuite réalisée. Les deux modèles utilisés pour les simulations montrent que l'assemblage cortex d'actine-membrane plasmique ne joue qu'un rôle minime dans la réponse mécanique cellulaire / Motivated by the study of biological tissues, this work contributes to developing numerical tools to predict the mechanical response of nonlinear heterogeneous materials in which the energies of interfaces can no longer be ignored. First, a computational homogenization strategy combined with a model reduction technique based on the proper orthogonal decomposition is implemented in the cases of large elastic deformations and highly nonlinear conduction. The interfaces between the different phases of a composite are described by means of a coherent interface model and taken into account numerically by an extended finite element method in tandem with a level-set technique. Finally, experimental results of single cell spreading between two fixed parallel microplates are exploited through finite element modelling. Our two models show that the bilayer membrane and the actin cortex do not play a significant role in the cell mechanical response

Dans une première partie, on s'intéresse à un système d'EDP stochastiques variant selon deux échelles de temps, et plus particulièrement à l'approximation de la composante lente à l'aide d'un schéma numérique efficace. On commence par montrer un principe de moyennisation, à savoir la convergence de la composante lente du système vers la solution d'une équation dite moyennée. Ensuite on prouve qu'un schéma numérique de type Euler fournit une bonne approximation d'un coefficient inconnu apparaissant dans cette équation moyennée. Finalement, on construit et on analyse un schéma de discrétisation du système à partir des résultats précédents, selon la méthodologie dite HMM (Heterogeneous Multiscale Method). On met en évidence l'ordre de convergence par rapport au paramètre d'échelle temporelle et aux différents paramètres du schéma numérique- on étudie les convergences au sens fort (approximation des trajectoires) et au sens faible (approximation des lois). Dans une seconde partie, on étudie une méthode d'approximation de solutions d'EDP paraboliques, en combinant une approche semi-lagrangienne et une discrétisation de type Monte-Carlo. On montre d'abord dans un cas simplifié que la variance dépend des pas de discrétisation- enfin on fournit des simulations numériques de solutions, afin de mettre en avant les applications possibles d'une telle méthode. / In a first part, we are interested in the behavior of a system of Stochastic PDEs with two time-scales- more precisely, we focus on the approximation of the slow component thanks to an efficient numerical scheme. We first prove an averaging principle, which states that the slow component converges to the solution of the so-called averaged equation. We then show that a numerical scheme of Euler type provides a good approximation of an unknown coefficient appearing in the averaged equation. Finally, we build and we analyze a discretization scheme based on the previous results, according to the HMM methodology (Heterogeneous Multiscale Method). We precise the orders of convergence with respect to the time-scale parameter and to the parameters of the numerical discretization- we study the convergence in a strong sense - approximation of the trajectories - and in a weak sense - approximation of the laws. In a second part, we study a method for approximating solutions of parabolic PDEs, which combines a semi-lagrangian approach and a Monte-Carlo discretization. We first show in a simplified situation that the variance depends on the discretization steps. We then provide numerical simulations of solutions, in order to show some possible applications of such a method.

Cette thèse traite de la solution numérique de deux types de problèmes stochastiques. Premièrement, nous nous intéressons aux EDS fortement oscillantes, c'est-à-dire, les systèmes composés de variables ergodiques évoluant rapidement par rapport aux autres. Nous proposons un algorithme basé sur des résultats d'homogénéisation. Il est défini par un schéma d'Euler appliqué aux variables lentes couplé avec un estimateur à pas décroissant pour approcher la limite ergodique des variables rapides. Nous prouvons la convergence forte de l'algorithme et montrons que son erreur normalisée satisfait un résultat du type théorème limite centrale généralisé. Nous proposons également une version extrapolée de l'algorithme ayant une meilleure complexité asymptotique en satisfaisant les mêmes propriétés que la version originale. Ensuite, nous étudions la solution des EDS de type McKean-Vlasov (EDSPR-MKV) associées à la solution de certains problèmes de contrôle sous un environnement formé d'un grand nombre de particules ayant des interactions du type champ-moyen. D'abord, nous présentons un nouvel algorithme, basé sur la méthode de cubature sur l'espace de Wiener, pour approcher faiblement la solution d'une EDS du type McKean-Vlasov. Il est déterministe et peut être paramétré pour atteindre tout ordre de convergence souhaité. Puis, en utilisant ce nouvel algorithme, nous construisons deux schémas pour résoudre les EDSPR-MKV découplées et nous montrons que ces schémas ont des convergences d'ordres un et deux. Enfin, nous considérons le problème de réduction de la complexité de la méthode présentée tout en respectant la vitesse de convergence énoncée. / This Ph.D. thesis deals with the numerical solution of two types of stochastic problems. First, we investigate the numeric solution to strongly oscillating SDEs, i.e. systems in which some ergodic state variables evolve quickly with respect to the remaining ones. We propose an algorithm that uses homogenization results and consists of an Euler scheme for the slow scale variables coupled with a decreasing step estimator for the ergodic averages of the fast variables. We prove the strong convergence of the algorithm as well as a generalized central limit theorem result for the normalized error distribution. In addition, we propose an extrapolated version applicable under stronger regularity assumptions and which satisfies the same properties of the original algorithm with lower asymptotic complexity. Then, we treat the problem of solving decoupled Forward Backward Stochastic Differential equations of McKean-Vlasov type (MKV-FBSDE) which appear in some stochastic control problems in an environment of a large number of particles with mean field interactions. As a first step, we propose a new algorithm, based on the cubature method on Wiener spaces, to weakly approach the solution of a McKean-Vlasov SDE. It is deterministic and can be parametrized to obtain any given order of convergence. Using this first forward approximation algorithm, we construct two procedures to solve the decoupled MKV-FBSDE and show that they converge with orders one and two under appropriate regularity conditions. Finally, we consider the problem of reducing the complexity of the presented method while preserving the presented convergence rates.

De nos jours, les antennes sont de plus en plus complexes en raison notamment de la nécessité de réaliser une reconfigurabilité en fréquence et/ou en diagramme. Les réseaux réflecteurs et les surfaces sélectives en fréquence sont des candidats particulièrement intéressants pour couvrir les besoins actuels. Cependant, en raison de leur grande taille et de la complexité géométrique croissante de leurs cellules élémentaires, l‘analyse électromagnétique complète de ces structures rayonnantes nécessite énormément de ressources informatiques (mémoire) et exige des temps de calcul prohibitifs, notamment lorsque des éléments de commande tels que des MEMS-RF sont intégrés au sein des cellules. Les techniques numériques classiques basées sur un maillage (spatial ou spectral) systématique ne parviennent pas à simuler de manière efficace de telles structures multi-échelles et nécessitent souvent des ressources informatiques difficiles d’accès pour le concepteur d'antennes. Une technique originale baptisée « Scale Changing Technique (SCT) » tente de résoudre ce problème en segmentant le réseau en de multiples domaines imbriqués les uns dans les autres et présentant divers niveaux d'échelle. Le multi-pôle par changement d’échelle, appelé « Scale Changing Network (SCN) », modélise le couplage électromagnétique entre deux niveaux d’échelle successifs. Ce multi-pôle peut être calculé en résolvant les équations de Maxwell à partir d’une Formulation par Equations Intégrales. La mise en cascade des multi-pôles par changement d’échelle permet alors le calcul de la matrice impédance (ou admittance) de surface du réseau complet. Cette matrice peut à son tour être utilisée pour simuler la diffusion électromagnétique d’une onde incidente par le réseau. Le calcul des différents multi-pôles par changement d’échelle peut être effectué séparément de sorte que le temps de simulation du réseau complet peut être considérablement réduit en parallélisant le calcul. Par ailleurs, la modification de la géométrie de la structure à une échelle donnée, lors de la phase de conception, nécessite seulement le calcul de deux multi-pôles par changement d’échelle et ne requiert pas une nouvelle simulation de toute la structure. Cette caractéristique fait de la SCT un outil de conception modulaire. Dans le cadre de cette thèse, la SCT a permis de tenir compte de la taille finie des réseaux et de modéliser efficacement les couplages électromagnétiques entre les cellules élémentaires. Des réseaux réflecteurs uniformes et non uniformes ont été simulés par la SCT et les performances numériques de la méthode ont été analysées. / Future antenna architectures especially for space applications are becoming more and more complex due to the need of reconfigurability. This reconfigurability is needed in terms of frequency, reliability, radiation pattern and power consumption. In this context, reflectarrays and frequency selective surfaces (FSSs) are particularly the hottest domains of RF design. The accurate analysis of electromagnetic (EM) scattering from such type of complex finite-sized reflectarray antenna structures is of great practical interest. However due to their large electrical size and complex cellular patterns specially when tuning elements such as RF-MEMS are also integrated within the array elements, conventional full-wave EM analysis of such multiscale structures either fail or require enormous amount of computational resources to resolve prohibitively large number of unknowns. Moreover the characterization of large structures would normally require a second step for optimization and fine-tuning of several design parameters, as the initial design procedure assumes several approximations. Therefore a full-wave analysis of the initial design of complete structure is necessary prior to fabrication to ensure that the performance conforms to the design requirements. A modular analysis technique which is capable of incorporating geometrical changes at individual cell-level without the need to rerun the entire simulation is extremely desirable at this stage. An indigenous technique called Scale Changing Technique (SCT) addresses this problem by partitioning the cellular reflectarray geometry in numerous nested domains and subdomains defined at different scale-levels in the array plane. Multi-modal networks, called Scale Changing Networks (SCNs), are then computed to model the electromagnetic interactions between any two successive partitions by method of moments (MoM) based integral equation approach. The cascade of these networks allows the computation of the equivalent surface impedance matrix of the complete array which in turn is utilized to compute far-field radiation patterns. Full-wave analysis of both passive and active (electronically tunable by RF-MEMS) reflectarrays has successfully been performed by the SCT while utilizing very small amount of computational resources as compared to conventional full wave methods. Moreover, to speed up the SCT modeling of the reflectarrays, equivalent electrical circuit models have been extracted and applied for individual design and optimization of the reflectarray phase shifter elements.

Les travaux de cette thèse se sont intéressés à l’obtention du comportement d’un composite sergé 2x2 verre/PA66 via un outil numérique basé sur une double homogénéisation : la première concerne les torons et la seconde le volume élémentaire représentatif du composite. A partir d’une campagne de caractérisation expérimentale sur le PA66, nous avons dans un premier temps identifié l’influence à la fois de l’hygrométrie et de la température sur le comportement de la matrice en quasi-statique. Ensuite, nous avons élaboré, implémenté et validé une loi de comportement isotrope élasto-plastique endommageable. Cette loi a servi à identifier par homogénéisation le comportement des torons en supposant un collage parfait des fibres et de la matrice ainsi qu’une homogénéité de l’eau dans le toron. Grâce à ces comportements identifiés et à l’élaboration d’une loi de comportement anisotrope élasto-plastique endommageable, nous avons pu déterminer les caractéristiques élastiques dans un premier temps et les comportements longitudinaux et en cisaillement dans un second temps. Les comparaisons aux résultats expérimentaux menés sur le composite offrent des résultats satisfaisants et permettent de penser que l’outil numérique développé, permettrait, à termes, d’aider à une conception rapide incluant ce genre de matériau en diminuant le nombre d’essais expérimentaux à faire. Des perspectives enfin sont proposées notamment pour l’extension aux comportements en dynamique (crash). / The objective of this PhD thesis is to create a methodology to simulate the behavior of a 2x2 glass/PA66 twill composite using homogenization. To achieve this goal, we use two steps: first one is dedicated to represent the yarn behavior and the second one to represent the behavior of the composite RVE. An experimental characterization campaign about the PA66 has been done in order to quantify the influence of both hygrometry and temperature on the behavior of the matrix in quasi-static. Then, we formulated, implemented and validated an isotropic elasto-plastic damageable constitutive law. It was used to identify by homogenization the behavior of the yarns by assuming a perfect bonding of the fibers and the matrix as well as a water homogeneity inside the yarn. Thanks to the identified behavior of the yarns and by formulating a compatible constitutive law (anisotropic elasto-plastic damageable), we have determined for the woven composite on one hand the elastic properties and in the other hand the longitudinal/transversal and shear nonlinear behaviors. The comparisons of the numerical and experimental results offer satisfactory results. It provides good outlook in short-term in regard of structures design which include this type of woven composite materials. The main advantage of the methodology is to decrease the number of required experimental tests. Some outlook dedicated to crash studies are finally suggested.

Le présent travail de thèse se consacre au développement d'un nouveau modèle micromécanique pour les composites en thermoplastique renforcé par des fibres de verre courtes. L'objectif est notamment la modélisation du comportement visco-endommageable en fatigue du PA66-GF30. Ce matériau, particulièrement utilisé dans l'industrie automobile, est sujet à une microstructure spécifique issue du procédé de moulage par injection. L'approche multi-échelles développée consiste en une méthode de Mori-Tanaka modifiée, appliquée à des renforts avec enrobage et prenant en compte l'évolution de l'endommagement à l'échelle microscopique. La description des mécanismes d'endommagement se base sur une investigation expérimentale poussée préalablement menée au sein de l'équipe. Des scénarios d'endommagement ont été proposés et incluent trois processus locaux différents : la décohésion de l'interface, la microfissuration de la matrice et les ruptures de fibres. Ceux-ci sont spécialement affectés par la microstructure. L'approche développée intègre ces cinétiques d'endommagement ainsi que la viscoélasticité non-linéaire de la matrice et la distribution d'orientation des inclusions due au procédé de fabrication. Chaque mécanisme d'endommagement est modélisé par une loi d'évolution basée sur les contraintes locales calculées à l'échelle microscopique. La loi constitutive finale, à l'échelle du volume élémentaire représentatif, est implémentée dans une bibliothèque scientifique en C++, SMART+, et est conçue pour être compatible avec une analyse de structures par éléments finis. L'identification du modèle est réalisée par rétro-ingénierie, en tirant profit de résultats expérimentaux multi-échelles, dont notamment des tests in-situ au MEB ainsi qu'une analyse qualitative et quantitative par μCT. / The current work focuses on a new micromechanical high cycle fatigue visco-damage model for short glass fiber reinforced thermoplastic composites, namely: PA66/GF30. This material, extensively used for automotive applications, has a specific microstructure which is induced by the injection process. The multi-scale developed approach is a modified Mori-Tanaka method that includes coated reinforcements and the evolution of micro-scale damage processes. Their description is based on the experimental investigations of damage mechanisms previously performed by the team. Damage chronologies have been proposed involving three different local degradation processes: fiber-matrix interface debonding/coating degradation, matrix microcracking and fiber breakage. Their occurrence strongly depends on the microstructure. The developed model integrates these damage kinetics and accounts for the complex matrix viscoelasticity and the reinforcement orientation distributions induced by the process. Each damage mechanism is introduced through an evolution law involving local stress fields computed at the microscale. The developed constitutive law at the representative volume element scale is implemented into a C++ scientific library, SMART+, and is designed to work with Finite Element Methods. The model identification is performed via reverse engineering, taking advantage of the multiscale experimental results: in-situ SEM tests as well as quantitative and qualitative μCT investigations.

La surveillance de la qualité de l'air est actuellement effectuée avec des mesures de concentration et à partir d'outils de modélisation de la dispersion atmosphérique. Ces modèles numériques évaluent les concentrations des polluants avec une résolution spatio-temporelle plus fine que les mesures. Néanmoins, les estimations fournies par ces modèles sont moins précises que les mesures. Dans ce projet de recherche, nous avons étudié les approches de couplage d'échelles et d'assimilation de données pour améliorer les estimations fournies par le modèle de dispersion atmosphérique SIRANE, dédié à l'échelle urbaine. L'approche de couplage d'échelles consiste à déterminer les conditions aux limites d'une simulation à partir d'une autre simulation à plus grande échelle. Au cours de ce travail de thèse, nous avons analysé trois méthodes afin de coupler le modèle urbain SIRANE et le modèle à méso-échelle CHIMERE. Cette étude montre que ces méthodes permettent potentiellement d'estimer la qualité de l'air à l'échelle urbaine de manière plus satisfaisante que les modèles à méso-échelle (utilisés seuls). Cependant, elles n'améliorent pas forcément la modélisation des conditions aux limites d'une simulation à l'échelle urbaine et les estimations fournies par celles-ci. Cela est a priori lié au fait que les estimations fournies par le modèle CHIMERE ne sont pas suffisamment satisfaisantes sur notre cas d'étude. Il est néanmoins possible que ces méthodes améliorent les résultats à l'échelle urbaine en utilisant une simulation à l'échelle régionale de meilleure qualité. L'approche d'assimilation de données consiste à combiner les mesures et les données modélisées afin de déterminer la meilleure estimation de l'état d'un système. Durant cette thèse, nous avons étudié trois méthodes d'assimilation de données : la méthode de débiaisement, la méthode que nous avons nommée modulation de la contribution des sources et la méthode Best Linear Unbiased Estimator. Cette étude indique que ces méthodes permettent globalement d'améliorer les estimations fournies par le modèle SIRANE. L'étude de sensibilité vis-à-vis du nombre de mesures utilisées lors de l'assimilation de données indique qu'en général, plus ce nombre est élevé plus les résultats sont satisfaisants. Enfin, les résultats montrent que les performances statistiques associées à ces trois méthodes d'assimilation de données sont globalement comparables entre elles sur notre cas d'étude. / Air quality monitoring is currently carried out with concentration measurements and with atmospheric dispersion modeling tools. These numerical models evaluate pollutant concentrations with a finer spatio-temporal resolution than measurements. Nevertheless, the estimates provided by these models are less accurate than measurements. In this research project, we studied multiscale coupling and data assimilation approaches to improve the estimates provided by the SIRANE atmospheric dispersion model, dedicated to the urban scale. The multiscale coupling approach consists in determining the boundary conditions of a simulation from another simulation on a larger scale. In this thesis work, we analyzed three methods for coupling the SIRANE model with the CHIMERE mesoscale model. This study shows that these methods can potentially estimate the air quality at the urban scale more satisfactorily than the mesoscale models (used alone). However, they do not necessarily improve the modeling of the boundary conditions of a simulation at the urban scale and the estimates provided by them. This is a priori due to the fact that the estimates provided by the CHIMERE model are not sufficiently good on our case study. It is possible, however, that these methods improve the results at the urban scale by using a better simulation at the regional scale. The data assimilation approach consists of combining the measurements and the modelled data to determine the best estimate of the system state. During this thesis, we studied three data assimilation methods : the unbiased method, the method that we called source apportionment modulation, and the Best Linear Unbiased Estimator method. This study indicates that these methods generally improve the estimates provided by the SIRANE model. The sensitivity study on the number of measurements used during the data assimilation indicates that, in general, higher is this number, more satisfactory are the results. Finally, the results show that the statistical performances associated with these three data assimilation methods are globally comparable on our case study.

Les matériaux nanostructurés suscitent un intérêt qui va croissant en raison de leurs propriétés chimiques et physiquesexceptionnelles. A cause de la complexité et du coût des développements expérimentaux à l’échelle nano, la simulationnumérique devient une alternative de plus en plus populaire aux études expérimentales. Dans ce travail de thèse, nous avons essayé de combiner des simulations à l’échelle atomique avec de la modélisation en milieu continu pour évaluer la conductivité thermique et la réponse élastique de matériaux nanostructurés. Nous avons utilisé des simulations de dynamique moléculaire pour calculer la réponse mécanique et thermique des matériaux sur des volumes à l’échelle nano. Des méthodes de micromécanique et la méthode des éléments finis, qui utilisent la mécanique des milieux continus, ont permis d’évaluer les propriétés mécaniques des matériaux à l'échelle macroscopique. Les résultats obtenus par ces simulations numériques ont été ensuite comparés avec ceux issus de l’expérience. / Nanostructured materials are gaining an ongoing demand because of their exceptional chemical and physical properties. Due to complexities and costs of experimental studies at nanoscale, computer simulations are getting more attractive asexperimental alternatives. In this PhD work, we tried to use combination of atomistic simulations and continuum modeling for the evaluation of thermal conductivity and elastic stiffness of nanostructured materials. We used molecular dynamics simulations to probe and investigate the thermal and mechanical response of materials at nanoscale. The finite element and micromechanics methods that are on the basis of continuum mechanics theories were used to evaluate the bulk properties of materials. The predicted properties are then compared with existing experimental results.

Dans les grands systèmes de réfrigération, l’utilisation de réfrigérants comme lubrifiants, à la place des huiles, simplifie la conception du système, l’allège et réduit son impact environnemental. La très faible viscosité du réfrigérant conduit à une épaisseur de film séparant les deux surfaces comparable à leur rugosité. Néanmoins, des travaux expérimentaux avec le réfrigérant R1233zd suggèrent que la lubrification est possible dans ces conditions grâce à la formation d’une couche adsorbée sur la surface d’oxyde de fer. Les analyses expérimentales in situ dans le contact sont très difficiles à cause des conditions de fort confinement et haute pression. C’est pourquoi une approche numérique multi-échelles est développée, afin d’étudier l’impact des réactions physico-chimiques à l’interface réfrigérant-surface sur les propriétés de lubrification. La théorie de la fonctionnelle de la densité est utilisée pour quantifier au niveau quantique, l’adsorption d’une molécule de réfrigérant sur une surface d’oxyde de fer. Des énergies de liaison allant de -0.92 eV à -0.22 eV sont observées et reliées à différents cas d’adsorption. Ces résultats sont exploités pour paramétrer un champ de forces interfacial, qui prédit des structures moléculaires à l’interface, différentes de celles obtenues avec des potentiels basés sur les règles de mélange classiques. Des simulations de dynamique moléculaire utilisant ce champ de forces paramétré confirment l’existence d’une couche fortement adsorbée de R1233zd sur une surface d’oxyde de fer. Avec des surfaces atomiques lisses et seulement 2 nm d’épaisseur de film de réfrigérant, les couches adsorbées résistent à des pressions allant jusqu’à 4 GPa et des vitesses de cisaillement atteignant 100 m/s. Une valeur minimale de 5 molécules de réfrigérant par nm² est nécessaire à la formation de deux couches adsorbées à 0.5 GPa. De plus, des simulations en cisaillement avec une surface rugueuse prédisent une rupture totale du film à environ 13 GPa. / In large refrigeration systems, using the refrigerant as lubricant instead of oil can help to simplify the design, lighten the systems, and reduce their environmental impact. However, the very low viscosity of refrigerants leads to ultra-thin films separating the surfaces, with a thickness comparable to surface roughness. Nevertheless, experiments with the R1233zd refrigerant suggests that lubrication is still possible in that situation thanks to an adsorbed layer formed on iron oxide surfaces. Experimental in situ analysis area is very difficult because of high confinement and high pressure. That is why a multiscale numerical approach is developed here, to explore the impact of chemical reactions and physical processes at the refrigerant-surface interface on large-scale lubrication properties. Density functional theory is used to quantify the adsorption of a refrigerant molecule on an iron oxide surface at the quantum level. Binding energies ranging from -0.92 eV to -0.22 eV are measured and related to different adsorption cases. These results are then used to parametrize an interfacial force field, whose predictions of interfacial molecular structure differs from those obtained using potentials based on standard mixing rules. Large-scale molecular dynamics simulations involving this parametrized force field confirm the existence of a strongly adsorbed layer of R1233zd molecules on iron oxide surface. With atomically smooth surfaces, and a refrigerant film thickness as small as 2 nm, the adsorbed refrigerant layers resists pressures as high as 4 GPa and sliding velocities as high as 100 m/s. A minimum value of 5 refrigerant molecules per nm² is necessary to the formation of two adsorbed layers at 0.5 GPa. Moreover, sliding simulations with a rough surface reveal total film breakdown for ca. 13 GPa.