Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Cette thèse est consacrée à l'étude des écoulements de gaz raréfiés à travers divers micro-conduits de type circulaire et rectangulaire dans des conditions isotherme et stationnaire. L'objectif de la thèse est de contribuer à l'étude de l'interaction gaz-surface notamment en déterminant le coefficient d'accommodation de la quantité de mouvement pour différent matériaux de surface (Or, Silice, Acier inoxydable et Sulfinert) associés à différents types de gaz (hélium, azote, argon et dioxyde-de-carbone). Afin d'atteindre cet objectif, on adopte un triple point de vue : expérimental, théorique et numérique. L'aspect expérimental est réalisé par des mesures de débit massique à travers les micro-conduits, en utilisant la méthode dite « à volume constant ». L'aspect théorique original est développé à travers une nouvelle approche basée sur la résolution de l'équation de Stokes. Cette approche a permis d'écrire une expression analytique de débit massique en régime de glissement, qui prenne en compte les effets bidimensionnels dans une section de conduit rectangulaire. Cette approche complètement explicite, est conduite au deuxième ordre. Enfin l'aspect numérique permet de calculer le débit massique, en régimes transitionnel et moléculaire libre, en résolvant numériquement l'équation cinétique BGK linéarisée. La comparaison des mesures de débit massique avec l'équation analytique, en régime de glissement, ou avec les calculs numériques, en régimes transitionnel et moléculaire libre, nous a permis de déduire des coefficients de glissement et les coefficients d'accommodation correspondant à chaque couple gaz-surface dans tous les régimes de raréfaction. / This thesis is devoted to the study of rarefied gas flows through micro-channels of various cross sections (circular and rectangular) under isothermal and stationary conditions. The objective of this thesis is to contribute to the study of gas-surface interaction by determining the tangential momentum accommodation coefficient for different surface materials (gold, silica, stainless steel and Sulfinert) and associated to various gases (helium, nitrogen, argon and carbon-dioxide). To achieve this goal three aspects are considered: experimental, theoretical and numerical. The experimental aspect is considered by measuring the mass flow rate through microchannels using the constant volume technique. The theoretical aspect is considered by the development of a new approach based on the Stokes equations. This approach yields to the analytical expression of the mass flow rate in the slip regime, which takes into account the second order effects. The last aspect, numerical, is considered by the numerical simulations of the mass flow rate in the transitional and free molecular flow regimes by solving the linearized BGK kinetic model. The comparison between the measured mass flow rates and the analytically expressions in the slip regime or with the results of numerical simulations in the transitional and free molecular regimes enabled to deduce the tangential momentum accommodation coefficients corresponding to each pair gas-surface in all flow regimes.

Dans la plupart des applications concernant un fluide s'écoulant sur une surface solide, la condition de non-glissement est largement utilisée car elle est simple et produit des résultats en accord avec les expériences. Toutefois, cette condition de limite n'est plus appropriée lorsque l'écoulement considéré est à l'échelle micro ou nano-métrique. Pour modéliser cet effet à l'échelle macroscopique, les conditions aux limites de Navier ont été introduites, avec la longueur de glissement comme paramètre. Lorsque le fluide est un gaz, cette longueur est liée au coefficient d'accommodation tangentiel (TMAC) et au libre parcours moyen, selon le modèle de Maxwell. Le but de ce travail est de traiter systématiquement ce modèle par une approche multi échelle et de l'étendre en incorporant la morphologie et l'anisotropie de la surface. La thèse est composée de cinq chapitres. Après l'Introduction, les notions de base de la théorie cinétique des gaz, l'équation de Boltzmann et les solutions associées (Navier-Stokes-Fourier, Burnett, Grad, Direct Simulation Monte Carlo…) sont rappelées dans le chapitre 1. Les modèles d'interaction gaz-paroi ainsi que les modèles de glissement introduits dans le cadre le la mécanique des fluides sont aussi rappelés. Le chapitre se termine par la description de la méthode de calcul par dynamique moléculaire (MD) utilisée dans ce travail. Le chapitre 2 est dédié au développement d'une technique simple afin de simuler les écoulements induits par la pression. Le principe est de se baser sur les formules atomistiques du tenseur des contraintes (Irving Kirkwood, Méthode de Plan, Contraintes Virielles) et de modifier les conditions périodiques, tout en maintenant la différence entre l'énergie cinétique des atomes à l'entrée et à sortie du domaine de calcul. Plusieurs types de conduite sont étudiés avec cette technique. Les résultats (température, vitesses, …) sont discutés et comparés. Le chapitre 3 concerne l'étude du potentiel d'interaction gaz-paroi par la méthode ab-initio. Le code CRYSTAL 09 est utilisé pour obtenir le potentiel entre un atome d'argon (Ar) et une surface de platine (Pt) <111> en fonction de la distance. Ensuite, le potentiel atome/surface est décomposé en potentiel binaire et approché par une fonction analytique. Cette fonction est ensuite implémentée dans un code MD pour simuler les collisions gaz-paroi et déterminer le coefficient TMAC.Dans le chapitre 4, l'effet de morphologie est étudié. Le potentiel multi-corps Quantum Sutton Chen(QSC) est utilisé pour le solide Pt <100> et un des potentiels binaires étudié dans le chapitre précédent pour le couple Ar-Pt. Le potentiel QSC est nécessaire pour reproduire l'effet de surface qui affecte le résultat final. Différentes surfaces sont traitées : surface lisse, surface nanostructurée, surface aléatoire obtenue par déposition de vapeur (CVD). Le coefficient TMAC est déterminé de façon généralisée, c.à.d en fonction de l'angle du flux d'atomes incidents sur la surface. Les anisotropies de surface et le noyau de collision sont également examinés. Dans le chapitre 5, on propose un modèle de glissement anisotrope pour les fluides en fonction du tenseur d'accommodation. Le modèle est obtenu par les calculs analytiques approchés développés dans le cadre de la théorie cinétique. On a ainsi généralisé l'équation de Maxwell en montrant que le tenseur de longueur de glissement est directement lié au tenseur d'accommodation. Le modèle est en bon accord avec les résultats de la méthode MD. Concernant la simulation MD, on a développé une technique pour reproduire l'anisotropie du tenseur d'accommodation. Le mémoire de thèse se termine par une synthèse des résultats obtenus. Des perspectives pour de futures études sont proposées / In most applications concerning a fluid flowing over a solid surface, the no-slip velocity condition was widely used because it is simple and produces the results in agreement with experiments. However, this dynamical boundary condition is not appropriate when the flow under consideration is at a micro or nano length scale.In order to model this effect at the macroscopic scale, the Navier boundary conditions have been introduced, with the slip length as a parameter. When the fluid is a gas, this length is related to the tangential momentum accommodation coefficient (TMAC) and the mean free path, according to the Maxwell model. The aim of this work is to systematically address this model using a multi-scale approach and to extend it by incorporating both the morphology and the anisotropy of a surface. The thesis consists of five chapters. In Chapter 1, the basics of the kinetic theory of gases, the Boltzmann equation and related solutions (Navier-Stokes-Fourier, Burnett, Grad, Direct Simulation Monte Carlo ...) are briefly presented. The models of gas-wall interaction and slip models introduced in the fluid mechanics are also recalled. The chapter ends with a description of the computational method used for the molecular dynamics simulations performed in this work. Chapter 2 is dedicated to the development of a simple technique to simulate the pressure driven flows. The principle is to rely on the atomistic formulas of the stress tensor (Irving Kirkwood, Method of Plane, Virial Stress) and to modify the periodic conditions by maintaining the difference between the kinetic energy of the ingoing and outgoing particles of the simulation domain. Several types of channels are studied with this technique. The results (temperature, velocity ...) are discussed and compared. Chapter 3 deals with the study of the gas-wall interaction potential by the ab-initio method. The code CRYSTAL 09 is used to obtain the potential between an atom of argon (Ar) and a surface of platinum (Pt) <111> as a function of distance. Then the gas-wall potential is decomposed into binary potential and approached by an analytic function. This function is then implemented in a MD code to simulate the gas-wall collisions and determine the TMAC coefficient. In Chapter 4, the effect of morphology is studied. The multi-body Quantum Sutton Chen (QSC) potential is used for Pt <100> solid and the binary potential proposed in the previous chapter for the Ar-Pt couple is employed. The QSC potential is needed to reproduce the surface effects that affect the final results. Different surfaces are treated : smooth, nanostructured surface and, random surface obtained by Chemical vapor deposition (CVD). The TMAC is determined using a generalized approach, i.e. depending on the angle of incident flux of gas atoms on the surface. The surface anisotropy and the scattering kernel are also examined. In Chapter 5, we propose a model of anisotropic slip for fluids based on accommodation tensor. The model is obtained by the analytical approximate calculations developed in the framework of the kinetic theory. We thus generalize Maxwell's equation by showing that the slip length tensor is directly related to the accommodation tensor. The model is in good agreement with the MD results. Thanks to our MD simulations, we develop a suitable technique for reproducing the anisotropy of the accommodation tensor. The thesis ends with a conclusion section in which we suggest some perspectives for a continuation of this work

Ce travail de thèse porte sur l'étude expérimentale, numérique et théorique des écoulements isothermes dans les microtubes et les microcanaux pour différents gaz. Le volet proprement expérimental porte sur la mesure des débits. Ce volet est complété d'abord par la mise en oeuvre d'une méthode de Monte Carlo permettant d'atteindre le profil des vitesses dans des sections choisies. Ces résultats sont comparés a ceux donnés par différentes approches théoriques et numériques : approche NS continue en régime de glissement, approches cinétiques (Boltzmann linéarisé, BGK) en régime transitionnel ou en régime moléculaire libre. A travers ces comparaisons on détermine notamment le domaine de validité du régime de glissement (premier et second ordre), des valeurs du coefficient d'accommodation de la composante tangentielle de la quantité de mouvement, les valeurs du coefficient du second ordre en régime de glissement ainsi que les grandeurs dont il dépend. On analyse aussi le comportement asymptotique de l'écoulement (débit) quand le nombre de Knudsen devient très grand.