Cette thèse a pour objet l'étude du comportement asymptotique des solutions des équations des fluides de grades 2 et 3. Dans le premier chapitre, on étudie les profils asymptotiques au premier ordre des solutions des équations des fluides de grade 2 en dimension 3. On démontre que les solutions des équations des fluides de grade 2 convergent vers des solutions particulières et explicites des équations de la chaleur, lorsque le temps tend vers l'infini. Ce résultat montre en particulier que les fluides de grade 2 se comportent asymptotiquement comme les fluides newtoniens régis par les équations de Navier-Stokes. Pour cette étude, on utilise les variables d'échelles (ou variables autosimilaires), et on effectue des estimations d'énergies dans divers espaces fonctionnels, en particulier dans des espaces de Sobolev à poids polynomiaux. La description des profils asymptotiques est obtenue sous des conditions de petitesse sur les données initiales de l'équation.Le second chapitre de cette thèse traite des profils asymptotiques à l'ordre 1 des solutions des équations des fluides de grade 3 en dimension 2. À l'instar des résultats du premier chapitre, on obtient ici aussi la convergence des solutions de ces équations vers des solutions explicites des équations de la chaleur. Les outils utilisés pour cette étude sont semblables à ceux utilisés pour les fluides de grade 2 en dimension 3, à savoir les variables autosimilaires et des estimations d'énergies. Dans ce cas aussi, on conclut que les fluides de grade 3 se comportent asymptotiquement comme les fluides newtoniens.Dans le dernier chapitre, on étudie l'existence d'un attracteur pour les équations des fluides de grade 3 en dimension 2 avec des conditions périodiques. On considère donc les solutions faibles de ces équations à données initiales dans l'espace de Sobolev H¹. Ces solutions faibles définissent un semi-groupe généralisé. Ensuite, on montre que les solutions à données initiales dans H² possèdent un attracteur global pour la topologie H¹. Pour ce travail, on utilise un schéma de Galerkin, des estimations a priori et une méthode de monotonie. Les principales difficultés que l'on rencontre sont liées au peu de régularité des données initiales et au fait que l'on ne sait par si les solutions des équations des fluides de grade 3 à données H¹ sont uniques. / This thesis is devoted to the study of the asymptotic behaviour of the solutions of the second and third grades fluids equations. In the first chapter, we study the asymptotic profiles to the first order of the solutions of the second grade fluids equations in dimension 3. We show that these solutions behave asymptotically (when the time goes to infinity) like explicit solutions to the heat equations. This result shows in particular that the asymptotic behaviour of the fluids of grade 2 is the same as the one of the Newtonian fluids, modelized by the classical Navier-Stokes equations. For this study, we use scaled variables (also called self-similar variables), and we perform energy estimates in several functions spaces, including weighted Sobolev spaces. Notice that the first order asymptotic expansion that we obtain holds under smallness assumptions on the initial data.In the second chapter of this thesis, we study the asymptotic profiles to the first order of the solutions of the third grade fluids equations in dimension 2. As in the previous chapter, we establish the convergence of these solutions to explicit solutions to the heat equations. The methods that we use are very similar to the ones used in the case of the second grade fluids equations on in dimension 3, namely scaled variables and energy estimates. We also conclude that the fluids of grade 3 behave asymptotically in time like Newtonian fluids.The last chapter is devoted to the study of the existence of an attractor for the third grade fluids equations in dimension 2 with periodic boundary conditions. We consider the weak solutions of these equations with initial data in H¹. These weak solutions define a generalized semiflow on H¹. Then, we show that the solutions with initial data in H² admit a global attractor for the H¹-topology. To this end, we use a Galerkin method, a priori estimates and a monotonicity method. The main difficulties come from the lack of regularity on the solutions and from the fact that these solutions are not known to be unique.

Le prototypage virtuel des circuits pneumatiques de puissance, par exemple les circuits de freinage des véhicules industriels, constitue un enjeu important en raison de la complexité des écoulements en régime transsonique et des couplages entre les échelles locales et macroscopiques. Ces problèmes sont rencontrés lors de la conception, de la synthèse des commandes et de l'analyse des performances statiques et dynamiques de ces circuits et l'analyse. La mise au point des modèles numériques de ces systèmes induit des coûts et des temps importants par rapport à d'autres systèmes. La démarche proposée dans cette thèse repose sur la construction numérique de bases de données permettant de caractériser le comportement local et macroscopique d'un composant de circuit en fonction de la variation de certains paramètres physiques ou géométriques par rapport à un point de fonctionnement de référence. Les bases de données résultent de l'extrapolation de la solution des équations de Navier Stokes moyennées (RANS) pour le point de référence considéré obtenu à l'aide d’un logiciel de paramétrisation en mécanique des fluides (Turb’Opty). La contribution de cette thèse repose pour l'essentiel dans un travail d'analyse des solutions issues de la paramétrisation dans deux contextes différents: la tuyère De Laval et un élément "coude", des composants élémentaires de circuit. Nous avons montré que ces exemples "simples" conduisent déjà à des difficultés importantes en termes de paramétrisation du problème et du calcul des dérivées des champs aérodynamiques en raison de la taille du problème. Pour pallier cette difficulté, nous avons proposé de déraffiner le maillage et nous avons alors montré que cette démarche conduit parfois à déplacer ou à atténuer certains phénomènes (chocs). La deuxième contribution de ce travail repose sur l'évaluation de la qualité des solutions extrapolées, de leur domaine de validité et la construction des liens entre grandeurs locales et macroscopiques. Nous avons enfin proposé une démarche permettant de reconstruire la caractéristique en débit d'un composant à partir de la détermination de la solution extrapolée pour un nombre limité de points de référence. / Virtual prototyping pneumatic circuits for power transmission, for example braking circuits of trucks, is still a difficulty because of the complexity of the flow behavior in transonic conditions and of the coupling between local and macroscopic scales. These problems are met during system design, control synthesis and for static and dynamic performance analysis. Tuning accurate numerical models requires important costs and time when compared to other systems. The methodology proposed in this PhD thesis relies on numerically determining a data base that characterizes the local and macroscopic behavior of a circuit component according the variation from a reference point of some physical or geometrical parameters. The data bases are obtained from the extrapolation of the Mean Navier Stokes solution (RANS) for a given reference point with the help of a parametrization software dedicated to fluid mechanics (Turb'Flow). The main contribution of this thesis relies io the analysis of the solution obtained from the parametrization in two different cases: the De Laval nozzle and un "elbow" connecting element, which are elementary component in a circuit. We have shown that these two "simple" cases lead already to important difficulties in term of problem parametrization and calculation of the derivatives of the aerodynamic fields because of the problem dimension. In order to tackle this, we proposed to reduce the spatial discretization (mesh derefining) and we showed that this approach could sometimes lead to damp or move some phenomena (shocks). The second contribution of this work relies on evaluating the quality of the extrapolated solution and their validity domain, and on building links between local and macroscopic behavior. Finally, we proposed a method that allows the mass flow rate characteristic of a component to be determined from the calculation of the extrapolated solution issued from a limited number of reference points.

Cette travail de thèse est lié au besoin de modélisation des écoulements diphasiques en générateurs de vapeur (entrée liquide et sortie vapeur). La démarche proposée consiste à faire le choix d'une modélisation hybride de l'écoulement, en scindant la phase gaz en deux champs, modélisés de manières différentes. Ainsi, les petites bulles sphériques sont modélisées avec une approche dispersée classique avec le modèle eulérien à deux fluides, et les bulles déformées sont simulées à l'aide d'une méthode de localisation d'interface. Le travail effectué porte sur la mise en place, la vérification et la validation du modèle dédié aux larges bulles déformées, ainsi que le couplage entre les deux approches pour le gaz gaz, permettant des premiers calculs de démonstration utilisant l'approche hybride complète.

L'objectif de ce travail consiste à développer la méthode de résolution PGD (Proper Generalized Decomposition), qui est une méthode de réduction de modèle où la solution est recherchée sous forme séparée, à la résolution des équations de Navier-Stokes. Dans un premier temps, cette méthode est appliquée à la résolution d'équations modèles disposant d'une solution analytique. L' équation de diffusion stationnaire 2D et 3D, l'équation de diffusion instationnaire 2D et les équations de Burgers et Stokes sont traitées. Nous montrons que dans tous ces cas la méthode PGD permet de retrouver les solutions analytiques avec une précision équivalente au modèle standard. Nous mettons également en évidence la supériorité de la PGD par rapport au modèle standard en terme de temps de calcul. En effet, dans tous ces cas, laPGD se montre beaucoup plus rapide que le solveur standard (plusieurs dizaine de fois). La résolution des équations de Navier-Stokes isothermes et anisothermes est ensuite effectuée par une discrétisation volumes finis sur un maillage décalé où le couplage vitesse-pression a été géré à l'aide d'un schéma de prédiction-correction. Dans ce cas une décomposition PGD sur les variables d'espaces uniquement a été choisie. Pour les écoulements incompressibles 2D stationnaire ou instationnaire, de type cavité entrainée et/ou différentiellement chauffé, les résultats obtenus par résolution PGD sont similaires à ceux du solveur standard avec un gain de temps significatif (la PGD est une dizaine de fois plus rapide que le solveur standard). Enfin ce travail introduit une première approche de la résolution des équations de transferts par méthode PGD en formulation spectrale. Sur les différents problèmes traités, à savoir l'équation de diffusion stationnaire, l'équation de Darcy et les équations de Navier-Sokes, la PGD a montré une précision aussi bonne que le solveur standard. Un gain de temps a été observé pour le cas de l'équation de Poisson, par contre, concernant le problème de Darcy ou les équations de Navier-Stokes les performances de la PGD en terme de temps de calcul peuvent encore être améliorées.

Dans l'industrie d'extraction pétrolière, l'efficacité des séparateurs eau-huile pour la production offshore est cruciale. L'objet de ce travail est de mettre en place les outils numériques nécessaires à la modélisation du fonctionnement de ces systèmes. Les phénomènes physiques entrant en jeu sont principalement : la présence d'interfaces entre des fluides non miscibles, la viscosité de ces fluides, et les effets de tension superficielle. Les modèles physiques et numériques correspondants ont été implémentés dans le cadre de la méthode numérique SPH (Smoothed Particle Hydrodynamics) développée au L.M.F.. Cette méthode numérique appartient à la classe des méthodes particulaires (sans maillage), suivant une approche d'écoulement compressible et avec une résolution explicite. Pour modéliser au mieux les écoulements bifluides, la formulation historique de la SPH a été enrichie par deux approches différentes, développées simultanément. Chacune d'entre elles a été validée séparément. La physique supplémentaire a été rajoutée par des modèles communs qui ont été validés sur différents cas tests tels que l'écoulement de Poiseuille, les instabilités de Rayleigh-Taylor, des cas d'envahissement ou l'évolution de bulles dans un liquide. Ce dernier cas a permis la comparaison aux outils de conception utilisés dans le procédé d'ingénierie de SAIPEM S.A., par l'intermédiaire d'une validation sur la loi de Stokes. Finalement, les capacités de la méthode sont illustrées sur la séparation eau-huile dans un séparateur de géométrie simplifiée.

Cette thèse porte sur l'étude expérimentale et numérique des forces fluides s'exerçant sur un cylindre ou un faisceau de cylindres oscillant latéralement en écoulement axial. Les paramètres de ce système sont l'amplitude, la fréquence d'oscillation, la vitesse de l'écoulement axial, le confinement et le rapport entre la longueur et le diamètre du cylindre. L'objectif de la thèse est de déterminer l'amortissement fluide, c'est-à-dire la force dissipative, créé par l'écoulement axial. Les résultats des travaux de thèse ont pour application industrielle la détermination de l'amortissement fluide des assemblages combustibles dans le cœur des centrales nucléaire lors d'un séisme. On s'intéresse dans le cadre de la thèse uniquement aux configurations pour lesquelles la vitesse d'oscillation latérale est faible devant la vitesse d'écoulement axial. Dans une première partie, on étudie le cas d'un cylindre non confiné oscillant en écoulement axial. Deux méthodes différentes sont utilisées : une dynamique et une quasi-statique. En dynamique, l'amortissement fluide est mesuré par la décroissance des oscillations pour un cylindre en oscillations libres. Dans l'approche quasi-statique, l'amortissement se déduit de la force normale mesurée sur un cylindre incliné. La gamme des faibles ratios entre la vitesse latérale et la vitesse axiale correspond, en statique, à une gamme de faibles inclinaisons où le cylindre est en écoulement quasi-axial. Le cas de l'approche quasi-statique a été traité expérimentalement dans des expériences en soufflerie et numériquement avec des calculs CFD. L'analyse des forces fluides montre que pour des inclinaisons inférieures à 5°, une portance linéaire avec l'angle contribue majoritairement à l'amortissement. L'origine de la force de portance est discutée à partir des mesures de pression à la paroi et de vitesses fluides. Les résultats de l'approche quasi-statique sont comparés aux résultats donnés par les expériences dynamiques. Dans une seconde partie, une étude expérimentale d'un faisceau confiné de 40 cylindres oscillant dans un écoulement axial d'eau est réalisée. La force et le déplacement sont mesurés simultanément. Les coefficients de masse ajoutée et d'amortissement sont évalués et leur variation en fonction de l'amplitude, la fréquence d'oscillation et la vitesse de l'écoulement axial est caractérisée. On montre ainsi que le coefficient d'amortissement dépend d'un unique paramètre : le rapport entre l'amplitude de la vitesse d'oscillation et la vitesse de l'écoulement axial.

La réduction de charge de fluides frigorigènes dans les systèmes de production de froid est un enjeu important s'inscrivant dans les politiques environnementales sur la contribution de ces fluides à l'effet de serre. Une des voies menant à la réduction de charge est l'augmentation de la compacité des échangeurs de chaleur, conduisant inévitablement à la conception d'évaporateurs plus complexes. Néanmoins, si de nombreuses études ont été publiées sur l'analyse thermo-hydraulique de l'ébullition convective de frigorigènes dans les tubes horizontaux, très peu ont été menées sur des écoulements diphasiques dans des géométries non conventionnelles. Cette thèse vise à étudier les caractéristiques des écoulements diphasiques avec changement de phase des frigorigènes de troisième (HFC) et de quatrième génération tels que le HFO-1234yf dans des singularités (coudes de retour, contractions). Ainsi, un banc d'essais expérimental a été spécifiquement conçu et construit pour mener des expériences sur l'ebullition convective. Ce banc d'essais a permis de visualiser les régimes d'écoulement et de mesurer les chutes de pression dans les singularités. Des visualisations de régimes d'écoulement ont été réalisées pour les frigorigènes HFO-1234yf, R-134a et R-410A dans des tubes droits et dans des singularités. Les résultats ont été confrontés avec diverses méthodes de prédiction de la litérature. Les perturbations sur l'écoulement causées par ces singularités et leurs effets sur le comportement hydrodynamique des fluides frigorigènes ont également été etudiés. Une serie d'éxperience a été menée pour déterminer le longueur de perturbation en amont et en aval des singularités. Finalement, les valeurs expérimentales des chutes de pression ont été confrontées à différentes corrélations de la littérature. Des nouvelles méthodes de prédiction de la chute de pression pour des écoulements diphasiques dans des singularités sont présentées. / The refrigerant charge reduction in HVAC\R systems is an important issue because it falls within environmental policies regarding refrigerants contributions to the greenhouse effect. A way to move toward charge reduction is to increase the compactness of heat exchangers, which means more complex designs of the evaporators. Nevertheless, while a large amount of studies have has been published on the thermal and hydraulic analysis of flow boiling of refrigerants in horizontal tubes, very little attention has been given to flow boiling in geometries different from straight tubes. This PhD thesis aims at studying the flow boiling characteristics of third generation (HFC) and fourth generation fluids such as HFO-1234yf in geometries which modify the fluid dynamics and two-phase flow with respect to horizontal straight tubes. To achieve this goal, an experimental test facility was specifically designed and built to conduct refrigerant evaporation experiments. This test facility allowed to perform flow regimes visualizations and pressure drop measurements in singularities (such as sudden contractions and return bends). First, two-phase flow regimes visualizations have been carried out using HFO-1234yf, R-134a and R-410A either in straight tubes or in singularities. A qualitative analysis of the flow behavior and also several comparisons to flow pattern prediction methods from the literature were conducted. The second objective of this work was to characterize the flow disturbances caused by singularities such as sudden contractions and return bends, and to study their effects on the hydrodynamic performance (e.g. pressure drop) of refrigerants.Lastly, pressure drop databases for third and fourth generation refrigerants are presented.

Cette thèse porte sur l'écoulement à l'intérieur d'une turbine de type Hélice. Deux méthodes de mesures ont été utilisées afin de remplir les objectifs de l'étude. La précision des instruments et l'incertitude des mesures effectuées à l'aide d'outils de vélocimétrie laser tels que la vélocimétrie Doppler laser (LDV) et la vélocimétrie par images de particules (PIV) ont été analysées. Trois phénomènes principaux sont étudiés à l'aide du LDV et du PIV. Le premier porte sur le débalancement circonférentiel de l'écoulement dans toute la machine. Le second porte sur les interactions rotor-stator se produisant entre les directrices et les aubes de la roue. Le troisième englobe les écoulements à l'entrée de l'aspirateur, le vortex à charge partielle (torche) et à forte charge ainsi que le sillage des directrices. Une attention particulière est portée sur les fréquences caractéristiques trouvées à l'aide des mesures LDV à cet endroit. Parallèlement aux mesures expérimentales, des simulations numériques ont été effectuées et comparées avec les données expérimentales disponibles. Cette comparaison permet d'améliorer les simulations numériques en posant différentes conditions limites et modèles de turbulence. Une meilleure utilisation et une amélioration des modèles RANS, URANS et LES avec l'utilisation de la base de données expérimentales est discutée dans l'étude. En conclusion à la recherche, le débalancement circonférentiel de l'écoulement est identifié avec les méthodes expérimentales et numériques. Le débalancement est de l'ordre de 10% de la vitesse de référence crête à crête aux avants directrices et est amorti à 5% à l'interface. Le débalancement maximal est de 12% à l'intérieur du cône divergent de l'aspirateur. L'intensité turbulente dans les canaux d'admission a été estimée à 8%. Les résultats expérimentaux ont démontré que les interactions rotor-stator sont faibles à l'intérieur de ce type de machine car il y a une dissipation rapide du sillage des directrices avant la roue. La torche à charge partielle et les recirculations se produisant dans l'aspirateur sont présentées avec leurs effets sur l'efficacité hydraulique et la configuration de l'écoulement. Plusieurs structures et vortex sont identifiés tels que le vortex se produisant dans l'espace en porte-à-faux des directrices, le vortex de bout d'aubes, la torche et les vortex contrarotatifs des canaux d'admission. Une base de donnée a été crée pour des recherches futures portant sur l'écoulement d'une turbine Hélice.

La thèse s'appuie sur une cartographie et un échantillonnage structural et pétrographique systématique dans une zone de l'Himalaya Central du Népal, entre les transversales de la Kali Gandaki à l'Ouest (E83°30') de la Trisuli Gandaki à l'Est (E85°20'), sur une surface d'environ 9000 km2. Cette zone correspond essentiellement aux formations métamorphiques du Moyen Himalaya et du Cristallin du Haut Himalaya, séparées par un chevauchement majeur et omniprésent en Himalaya, le Grand Chevauchement Central (ou MCT). La thèse est découpée en quatre parties, traitant respectivement de la lithologie des grands ensembles structuraux, du métamorphisme puis de la déformation dans la zone du MCT, et enfin d'une synthèse sur la nature et le fonctionnement du MCT. La première partie commence par une description des formations qui composent le Cristallin du Haut Himalaya, que l'on peut ici subdiviser en trois formations, métapélites alumineuses à la base, gneiss à minéraux calciques au dessus, orthogneiss (anciens granites à env. 500 Ma) au sommet. On montre qu'il y a d'importantes variations d'épaisseur latérale (épaisseur totale augmentant de 5 à 12 km d'Ouest en Est). Les formations du Moyen Himalaya (Formations du Moyen Pays Népalais), en dessous du MCT, sont séparées en deux ensembles : ensemble inférieur, série de grès épaisse de 2000 à 4000m, se terminant par des niveaux quartzitiques, et dessinant un grand anticlinal; ensemble supérieur, comprenant de manière caractéristique des niveaux carbonatés et carbonés (schistes noirs), préservé sur les flancs Nord et Sud de l'anticlinal. Au Sud, le métamorphisme et la déformation sont peu intenses, et l'auteur propose que ces formations soient la suite stratigraphique normale (la couverture) de l'ensemble gréso-quartzitique inférieur. Des restes de fossiles sont mentionnés dans les schistes noirs, mais non datés. Au Nord, juste sous le MCT, ces formations sont métamorphiques et très déformées (elles sont reprises dans la déformation associée au MCT), mais se parallélisent néanmoins bien avec les formations méridionales. La seconde partie est consacrée à l'étude du métamorphisme dans la zone du chevauchement central, métamorphisme caractérisé par une disposition "inverse" (i.e. métamorphisme décroissant lorsque l'on descend dans la pile structurale) qui a intrigué de longue date les géologues himalayens. Un inventaire systématique des minéraux et paragénèses minérales est effectué, une attention toute particulière étant donnée au calage des minéraux de métamorphisme par rapport à la déformation. Ceci permet de définir cartographiquement différentes zones métamorphiques et les isogrades qui les séparent. L'amplitude du relief permet de déterminer la géométrie 3D de ces isogrades, qui sont en première approximation parallèles au MCT. Les conditions P-Tp associées sont définies par l'étude des paragénèses des roches alumineuses dans le système KFMASH. Les équilibres chloritoïde-grenat-disthène sont particulièrement discutés. Un trajet P-Tp à travers la zone du MCT est établi, trajet qui ne montre pas de discontinuité au passage du MCT, et un fort gradient thermique en dessous du MCT. Les données métamorphiques (thermiques, structurales, à échelle du minéral et à échelle de la carte) permettent ainsi d'établir le lien entre métamorphisme inverse et fonctionnement du MCT. Les données thermobarométriques tardi-métamorphiques sont précisées par une étude systématique des inclusions fluides (à H2O +/- CO2 et sels) contenues dans les lentilles de quartz d'exsudation syn-schisteuses. L'auteur montre que la structure barométrique dans la zone du MCT ne peut s'expliquer sans impliquer une forte érosion syn-fonctionnement du MCT. La troisième partie détaille la déformation de part et d'autre du MCT, de l'échelle de l'affleurement jusqu'à l'échelle du réseau : l'importance des structures mésoscopiques à symétrie monoclinique est soulignée, en particulier les amandes métriques à millimétriques associant plan de cisaillement et plan d'aplatissement. Une étude systématique de l'orientation préférentielle des axes C du quartz dans les niveaux de quartzites a été réalisée (en collaboration avec J.L. Bouchez) : comme à échelle de l'affleurement ou de la lame-mince, la fabrique de réseau a une symétrie monoclinique prédominante, de plus en plus forte quand on se rapproche du MCT. Cette partie décrit donc un certain nombre de structures de cisaillement, pour certaines inédites à l'époque, et les utilise pour montrer que le MCT n'est pas un simple chevauchement, mais une zone de cisaillement hétérogène d'épaisseur plurikilométrique. La quatrième partie récapitule les principaux éléments nouveaux apportés à l'histoire thermique et structurale du MCT. Toutes les nouvelles données du travail confirment la relation causale entre MCT et métamorphisme inverse, et ont des implications sur le caractère thermique chaud de la croûte indienne lors du chevauchement, et sur l'importance du rôle de l'érosion. Le modèle du "fer à repasser" (P. Le Fort) proposé à l'époque pour expliquer ce métamorphisme inverse est discuté et ses limites soulignées.

Dans cette thèse, nous adressons trois projets avec l'application de microfluidique Avec le Vélocimétrie de Réflexion Totale Interne, nous avons réalisés le nanoPTV des fluides à 800 nm près de parois du solide. Nous arrivons à une précision sans précédent, par la détermination précise de la position du parois, et par la simulation de Langevin, en tenant compte des nombreux sources de biais physique, comme le mouvement Brownien, effet du cisaillement, la répulsion électrostatique entres les particules et le parois, et la défocalisation de la lentille. Nous obtenons ±5 nm and ± 10 nm de précision sur la longueur de glissement pour la solution de sucrose et de l'eau. La condition de non-glissement sur la surface hydrophile est confirmée, et un glissement sur la surface hydrophobe est observé. Nous collaborons avec A. Leshansky pour étudier la formation des gouttes sur une intersection entre un canal confiné et un réservoir profond. Cette phénomène est appelé le "step emulsificaiton". La dynamique de la formation des gouttes est étudiée expérimentalement de façon approfondie. La théorie est basée sur la dynamique des fluides dans un canal Hele-Shaw, avec les effets de forces capillaires. Nous arrivons à expliquer le mécanique du fluides derrière la formation des gouttes, inclus les taille des gouttes. Nous collaborons avec un groupe des entreprises pétrolières (AEC), pour étudier le mouvement des nano particules dans un micro model de milieux poreux. Ces particules sont supposé de faire transition une fois en contact avec l'huile ou expériencer un changement de la température. L'injection des particules dans les réservoirs de l'huile et de gaz permet de collecter l'information sur la distribution et la quantité de l'huile et de gaz. Avant l'application en mass dans l'industrie, c'est favorable de les tester dans un micro model, qui possèdes une structures similaire aux pores des roches. Nous avons testé les nano particules synthétisés par les autres membres de l'AEC, et confirmé que l'idée du micro model est une méthode efficace de prédire la performance des particules sous sol.