Couche de Knudsen
La couche de Knudsen décrit la région proche d'une interface solide-gaz ou liquide-gaz pour laquelle le milieu situé à l'extérieur peut être décrit par les équations de Navier-Stokes. Dans cette région épaisse de quelques libres parcours moyens on sait résoudre analytiquement l'équation de Boltzmann. On établit ainsi une équation de saut qui est une relation entre une variable et sa dérivée normale à la surface. Ces équations deviennent des conditions à la limite pour la région décrite par les équations du continu[1].
Position du problème
Au voisinage de toute interface solide-gaz ou liquide-gaz une région de quelques libres parcours moyen raccorde la paroi et une distribution des vitesses et des énergies imposée par celle-ci et le milieu externe décrit par une distribution du type Chapman-Enskog, proche de la statistique de Maxwell-Boltzmann, caractéristique des écoulements décrits par les équations de Navier-Stokes. On retrouve ce problème pour une onde de choc qui n'est pas une réelle discontinuité à l'échelle du libre parcours moyen.
Toutefois l'étendue de cette région de l'espace est faible devant le rayon de courbure de la surface et on la suppose également faible devant toute dimension caractéristique portée par l'éventuelle ligne de courant qui jouxte cette surface. Ceci autorise une approche unidimensionnelle du problème[N 1].
Interaction pariétale
L'interaction d'un atome ou d'une molécule avec une paroi solide ou liquide est un problème difficile accessible seulement à la dynamique moléculaire. Dans la pratique on se contente de décrire deux cas extrêmes pour les vitesses :
- une réflexion spéculaire, équivalente à une condition de symétrie, et conduisant à un glissement (vitesse parallèle à la paroi non nulle) s'il existe un écoulement,
- une réflexion diffuse dans laquelle la particule est renvoyée aléatoirement et de manière isotrope avec une distribution maxwellienne des vitesses à la température de paroi qui est une donnée du problème. Cette composante est la cause du ralentissement de la vitesse.
On attribue un poids à la réflexion diffuse : le coefficient d'accomodation, un nombre compris entre 0 et 1. La part de la réflexion spéculaire est le complément à 1. Cette quantité est strictement positive, une réflexion totalement spéculaire étant dénuée de sens physique.
Pour les énergies internes on définit de la même façon un second coefficient d'accomodation donnant la fraction de particules réémises conformément à la statistique de Boltzmann à la température de paroi. Ce second coefficient est généralement confondu avec le premier, sauf dans les cas où l'on dispose de données précises issues d'expériences utilisant des faisceaux de particules à faibles vitesses représentatives du problème[2].
Résolution du problème
L'équation de Boltzmann est résolue dans le domaine compris entre la paroi et une hauteur arbitraire h à laquelle on suppose valide les équations du continu. On connait les conditions à la limite supérieure : il s'agit de paramètres du problème. Par ailleurs le bas du domaine obéit aux relations décrites ci-dessus.
La résolution est obtenue :
- soit en écrivant des relations de conservation pour les flux de particules montantes et descendantes prises séparément (masse, quantité de mouvement et énergie),
- soit en utilisant la méthode de Grad[3].
Les résultats sont des relations de saut liant la variable à sa dérivée en haut du domaine. Pour des particules sans énergie interne on obtient :
- pour la vitesse
- pour la température
avec
- axe normal à la paroi,
- vitesse parallèle à la paroi,
- coefficient d'accomodation,
- pression,
- conductivité.
Ces relations deviennent des conditions aux limites pour les équations de Navier-Stokes. Elles sont appliquées à la paroi : la couche de Knudsen est "écrasée". Cela suppose bien entendu qu'elle soit de très faible épaisseur.
Numériquement les sauts deviennent très faibles dès que la pression s'approche de la pression normale. C'est la justification de l'hypothèse traditionnelle pour les équations du continu : vitesse nulle à la paroi, température du gaz égale à celle de la paroi.
Cette approche peut être généralisée au cas d'un gaz comportant plusieurs espèces et des énergies internes[4].
Notes
- Le cas bidimensionnel est traité dans la référence 1.
Références
- (en) Mikhail N. Kogan, Rarefied Gas Dynamics, Springer, (ISBN 978-1-4899-6189-1)
- (en) S. C. Saxena et R. K. Joshi, Thermal Accommodation and Adsorption Coefficients of Gases, Hemisphere Publishing Corporation, (ISBN 0891168702)
- (en) Harold Grad, « On the Kinetic Theory of Rarefied Gases », Communications on Pure and Applied Mathematics, vol. 2, no 4, , p. 331-407
- (en) Roop N. Gupta, Carl D. Scott et James N. Moss, « Surface Slip equations for Multicomponent Nonequilibrium Air Flow », NASA Technical report CR-181252, (lire en ligne)