Méthode de Chapman-Enskog

À la fin du XIXe siècle on connaît l'équation de Boltzmann qui régit la dynamique du milieu gazeux à l'échelle microscopique et les équations d'Euler et de Navier-Stokes pour le niveau macroscopique. Le passage d'une échelle à l'autre constitue une partie du sixième problème de Hilbert. David Hilbert, auteur des énoncés des problèmes jugés majeurs à la fin du XIXe siècle pose les bases d'une méthode sous forme d'un développement qui porte son nom (1912). Il faudra attendre quelques années pour que Sydney Chapman et David Enskog proposent simultanément et indépendamment en 1916 et 1917 une solution à ce problème[1],[2],[3]. Plus récemment cette méthode a été étendue au cas d'un gaz en déséquilibre thermodynamique[4], ce dernier aspect étant aujourd'hui encore un domaine de recherche très actif.

La méthode de Chapman-Enskog est une méthode de perturbation consistant à définir la solution sous forme de série de fonctions de distribution en fonction d'un « petit paramètre » assimilable au nombre de Knudsen. À l'ordre zéro on retrouve la distribution de Maxwell-Boltzmann et les équations d'Euler. L'ordre un permet de connaître l'expression des flux de chaleur et de quantité de mouvement et celle des coefficients de transport (les coefficients de diffusion par gradients de concentration, de pression et de température, les viscosités dynamique et volumique, la conductivité) à partir des potentiels d'interaction moléculaires. Cette approche permet de retrouver les équations de Navier-Stokes et de justifier la diffusion par gradient thermique, inconnue à l'époque où sont publiés les travaux de Chapman et d'Enskog. Cette méthode permettra par la suite de calculer tous ces coefficients à partir de la connaissance de l'un d'entre eux en reconstituant à partir d'une mesure (généralement la viscosité) un potentiel d'interaction tel que le potentiel de Lennard-Jones[2].

Harold Grad a proposé une approche alternative consistant à chercher la solution par la méthode des moments de la fonction de distribution (1949). L'équation de Boltzmann est multipliée par ( est la vitesse microscopique de l'équation de Boltzmann et le produit tensoriel) et intégrée en vitesse. Dans ce type de méthode l'équation portant sur le ne moment fait apparaître le (n+1)e. Il faut donc faire une hypothèse pour « fermer » le système. Grad suppose la solution exprimée par une série tronquée de polynômes d'Hermite[5],[3]. David Levermore a plus récemment (1996) proposé une fermeture qui fait appel à une propriété générale : la solution maximise l'entropie du système de fermions que sont les particules du milieu[6]. Des codes de calcul[7] basés sur ces méthodes sont restés dans le domaine du laboratoire car n'apportant pas un gain notable en termes de domaine de validité (en termes de nombre de Knudsen) par rapport aux codes standard résolvant les équations de Navier-Stokes, lesquels ont fait l'objet de développements considérables.

La méthode de Chapman-Enskog a été étendue à l'équation de Boltzmann-Chernikov en relativité générale pour les applications en cosmologie[8].

Équations d'évolution microscopique et macroscopique

modifier

Niveau microscopique

modifier

On note   la fonction de distribution statistique de la vitesse   à l'instant   au point   pour la particule (atome ou molécule) appartenant à l'espèce  . Le nombre probable de particules dans le volume  , de vitesses   à cet instant est  . La distribution statistique   se mesure donc en s3 m−6.

  Remarque : certains auteurs utilisent une distribution de la quantité de mouvement   plutôt qu'une distribution de la vitesse.

L'équation de Boltzmann s'écrit (en l'absence de force extérieure)  , l'opérateur (ou noyau) de collision, est un opérateur intégral quadratique décrit ci-dessous, donnant l'effet des collisions que l'on supposera élastiques pour simplifier le problème : pas d'échange entre degrés de liberté internes et translation, pas de réaction chimique. On exclut donc la viscosité volumique qui résulte de ce type d'échange.

Il y a autant de distributions que d'espèces présentes dans le milieu. À chacune correspond une équation de Boltzmann couplée aux autres par les seconds membres qui représentent des collisions homogènes ( ) ou hétérogènes ( ).

La collision élastique

modifier
 
Schéma d'une interaction moléculaire élastique dans le système lié au barycentre.

Les vitesses avant interaction sont   et   dans un référentiel galiléen. Les indices représentent indifféremment une même espèce ou deux espèces différentes. Ces vitesses valent   et   après interaction. On se place dans un système centré sur le barycentre qui a une vitesse constante du fait de la conservation de la quantité de mouvement. Dans ce système qui est donc galiléen la vitesse initiale de la particule   est la vitesse relative  . Par symétrie on peut affirmer que la trajectoire sera contenue dans le plan contenant l'origine et  . On choisit un repère tel que   (voir figure). Dans ce repère la déviation est  , fonction du paramètre d'impact  , de la vitesse relative   et du potentiel d'interaction que l'on suppose ne dépendant que de la distance en les deux particules en interaction. Si cette hypothèse est rigoureuse pour l'interaction entre deux atomes, on peut la considérer utilisable pour deux molécules : le potentiel est alors un potentiel moyen statistique.

La direction de sortie d'interaction est définie par  . On peut calculer les vitesses finales à partir des considérations suivantes :

  • la conservation de la quantité de mouvement dans l'interaction implique que 
  • la vitesse relative   a un module constant du fait de la conservation de l'énergie, donc   ou   Les vitesses après l'interaction sont donc : 

De plus la conservation du moment cinétique au cours de l'interaction conduit à  . Le système décrivant la collision est réversible. Le théorème de Liouville permet donc d'écrire :  

Le noyau de collision

modifier

Le nombre probable de particules qui traversent l'aire   par unité de temps est  . Elles interagissent avec le nombre probable de particules dans le volume élémentaire  . Le nombre de particules qui disparaissent de la statistique par unité de temps est   avec : On compte de la même façon la quantité de particules qui apparaissent, soit : Compte tenu des relations données ci-dessus pour la collision l'opérateur de collision s'écrit : Cette équation est nommée équation de Wang Chang et Uhlenbeck. On peut donner une formulation équivalente en introduisant la section efficace différentielle   définie par[9]:  d'où 

Niveau macroscopique

modifier

Les variables

modifier

L'équation de Boltzmann décrit au niveau microscopique l'évolution de particules. Pour décrire chacune des espèces   au niveau macroscopique et l'ensemble de celles-ci on définit :

  • la densité particulaire   ;
  • la masse volumique   ;
  • la vitesse moyenne   ;
  • l'énergie interne   ;

On peut alors définir des valeurs pour l'ensemble des espèces, soit :

  • la vitesse d'ensemble   (vitesse barycentrique de l'ensemble des particules) ;
  • la vitesse de diffusion   (quantité cinématique) ;

Quelques variables auxiliaires (où   est le nombre d'Avogadro)

  • la fraction volumique   ;
  • la fraction massique   ;
  • la masse molaire  , de valeur moyenne   ;

Les flux

modifier

Le flux de la quantité   est par définition la quantité On définit ainsi, en notant   le produit dyadique :

  • le flux de masse   ;
  • le tenseur de pression   , qui représente le flux de quantité de mouvement. Il est symétrique par construction.
  • la pression partielle   , définie à partir de la trace du tenseur de pression.
  • le flux de chaleur   . Il représente le flux d'énergie interne.

Les flux globaux sont obtenus simplement en sommant sur  , ainsi que la pression  On peut alors définir une température à partir de l'équation d'état  

Équations d'évolution

modifier

En multipliant chacune des équations de Boltzmann successivement par chacun des invariants collisionnels   et en intégrant sur les vitesses et, s'il y a lieu, sur les espèces, on obtient les équations d'évolution macroscopiques appelées équations d'Enskog. On note   le produit contracté. Tous les seconds membres sont nuls du fait des lois de conservation : On obtient ainsi un système d'évolution pour  ,   et   dans lequel les flux  ,   et   restent à expliciter.

Équation de Boltzmann adimensionnelle

modifier

On se place dans l'hypothèse d'un milieu homogène (une seule espèce présente).

Afin d'estimer l'apport de chaque terme dans l'équation de Boltzmann il est nécessaire d'adimensionner celle-ci. Pour cela on définit les quantités de référence suivantes :

  • la température   d'où l'on déduit une vitesse de référence qui est la vitesse du son pour un gaz parfait   ;
  • la longueur de référence macroscopique   liée au problème avec laquelle on définit un temps de référence   lié à la propagation des ondes acoustiques ;
  • la pression de référence   avec lesquelles on définit une masse volumique de référence   et la quantité   ;
  • une vitesse de référence microscopique correspondant à la vitesse moyenne d'une particule   ;
  • une autre vitesse de référence microscopique correspondant à la vitesse relative moyenne entre deux particules   ;
  • une longueur de référence microscopique correspondante au libre parcours moyen    est la section efficace différentielle.

Si l'on définit à présent les variables réduites  ,  ,  ,  ,   et  , l'équation de Boltzmann s'écrit : 

  •   est le nombre de Strouhal qui pondère la variation temporelle du système,
  •   est le nombre de Knudsen dont l'inverse pondère le terme de collision, donc la tendance au retour vers l'équilibre. En pratique, pour que les équations du continu soient valides, il faut que   , typiquement.

Formulation de la solution

modifier

On écrit la solution comme une série en utilisant un paramètre   du même ordre de grandeur que le nombre de Knudsen :    respectant les lois de conservation, chacun des termes du développement doit également les respecter. D'où les contraintes sur la solution :  On porte cette approximation dans l'équation de Boltzmann et on sépare les termes correspondant à chaque puissance de  .

À l'ordre zéro

modifier

On obtient simplement  Cette équation est vérifiée si tous les termes qui la composent sont nuls, donc en particulier  ce qui implique que   ou que  Ainsi,   est donc un invariant collisionnel. Il s'écrit donc comme une combinaison linéaire des invariants collisionnels canoniques:  

En introduisant cette expression dans les équations définissant les variables macroscopiques on identifie les paramètres de ce développement et on retrouve la loi de distribution des vitesses de Maxwell, soit :  

avec  . Les flux de diffusion   sont nuls ainsi que le flux de chaleur  . Le tenseur de pression est réduit à sa trace    est le tenseur unitaire. Les équations macroscopiques correspondantes sont les équations d'Euler.

À l'ordre un

modifier

L'ordre un fait apparaître une équation intégrale de Fredholm pour l'inconnue  , soit :  La difficile résolution de cette équation[1],[2] permet de donner les diverses quantités inconnues des équations d'Enskog que l'on peut alors assimiler aux équations de Navier-Stokes.

Le flux de diffusion

modifier

Il s'obtient sous forme d'un système linéaire appelé système de Stefan-Maxwell :  où l'on voit apparaître le coefficient de diffusion binaire   et le « coefficient de diffusion thermique multicomposant »   (en fait un nombre sans dimension) lié aux coefficients ordinaires   (qui ne sont pas des coefficients de diffusion[N 1] et qui peuvent être négatifs) par :  Pour un milieu comportant   espèces, le rang de ce système est   puisque  . Sa solution formelle est : On retrouve les termes de diffusion par gradient de concentration, de pression et de température (effet Soret).   est le coefficient de diffusion multicomposant, solution d'un système linéaire faisant intervenir les coefficients binaires. Ce système étant également de rang   la solution n'est pas unique et comporte   termes indépendants. On choisit généralement   par souci de symétrie. Ce choix est arbitraire.

Il existe diverses solutions approchées du système de Stefan-Maxwell permettant d'obtenir une expression explicite du flux de diffusion sous une forme voisine de la loi de Fick, laquelle n'est exacte que pour un mélange binaire[10],[11].

Le tenseur de pression

modifier

Le tenseur de pression a une forme classique

 

  est le tenseur unité et   le tenseur des contraintes visqueuses

 

Un terme supplémentaire   apparaît dans la pression lorsque l'on prend en compte les interactions inélastiques. Son influence est faible voire totalement négligeable pour les gaz peu denses[4].

On notera que l'hypothèse de Stokes se justifie naturellement par cette approche.

Le flux de chaleur

modifier

Il est donné par

 

  est la conductivité thermique. Le dernier terme de l'équation est le corollaire de l'effet Soret et est nommé effet Dufour.

Les coefficients de transport

modifier

Les coefficients de transport s'expriment sous forme de systèmes linéaires qui font intervenir des quantités du type   qui sont développées en polynômes de Sonine-Laguerre. Les coefficients du développement s'expriment en fonction des intégrales de collision[1]. En pratique on se contente du premier ordre pour le développement et les intégrales de collision sont des fonctions de la température tabulées par divers auteurs[2]. De plus il existe des solutions approchées des systèmes linéaires donnant les divers coefficients sous forme explicite[10],[11].

Certains coefficients de transport sont analytiques lorsque l'on utilise un potentiel sphère dures. Ceux-ci sont[11] :

  • le coefficient de diffusion binaire variant en   ;
  • la viscosité variant en  , indépendante de la pression.

Les variations sont analogues lorsque l'on utilise un potentiel plus réaliste, par exemple un potentiel de Lennard-Jones : seul l'exposant de la température change un peu.

La conductivité a une forme analogue à celle de la viscosité pour un gaz monoatomique.

La fonction de distribution

modifier

La fonction de distribution est

 

  est donné dans le flux de diffusion et

 
 

Cette fonction de distribution est nécessaire pour le calcul de la couche de Knudsen qui donne les conditions à la paroi pour les équations de Navier-Stokes.

À l'ordre deux

modifier

On obtient là encore une intégrale de Fredholm pour l'inconnue  

 

David S. Burnett a proposé en 1936 une solution de cette équation[12]. Celle-ci a l'inconvénient de ne pas respecter le théorème H[13]. Il semble que la montée en ordre constitue une impasse, toutes les variantes proposées jusqu'à ce jour ne résolvant pas ce problème[14].

Notes et références

modifier
  1. Il serait plus logique de définir un « vrai » coefficient de diffusion thermique  .

Références

modifier
  1. a b et c (en) Sydney Chapman et Thomas George Cowling, The Mathematical Theory of Non-uniform Gases : an account of the kinetic theory of viscosity, thermal conduction, and diffusion in gases, Cambridge/New York/Port Chester etc., Cambridge University Press, , 422 p. (ISBN 0-521-40844-X)
  2. a b c et d (en) Joseph Oakland Hirschfelder, Charles Francis Curtiss et Robert Byron Bird, Molecular Theory of Gases and Liquids, John Wiley and Sons, (ISBN 978-0-471-40065-3)
  3. a et b (en) Gilberto Medeiros Kremer, « The Methods of Chapman-Enskog and Grad and Applications », RTO-EN-AVT 194,‎ 2011 [1]
  4. a et b Raymond Brun, Introduction à la dynamique des gaz réactifs, Toulouse, Cépaduès, , 402 p. (ISBN 978-2-36493-190-9)
  5. (en) Harold Grad, « On the Kinetic Theory of Rarefied Gases », Communications on Pure and Applied Mathematics, no 4,‎ , p. 331-407
  6. (en) Charles David Levermore, « Moment Closure Hierarchies for Kinetic Theories », Journal of Statistical Physics, vol. 83,‎ (lire en ligne)
  7. Patrick Le Tallec et Jean-Philippe Perlat, « Numerical Analysis of Levermore’s Moment System », sur HAL,
  8. (en) G. M. Kremer, « The Boltzmann Equation in Special and General Relativity », AIP Conference Proceedings, vol. 160, no 1501,‎ (lire en ligne)
  9. Bases de la mécanique physique
  10. a et b (en) Joel Henry Ferziger et Hans Gerard Kaper, Mathematical Theory of Transport Processes in Gases, North Holland Publishing, (ISBN 0-7204-2046-6)
  11. a b et c (en) Georges Duffa, Ablative Thermal Protection Systems Modeling, Reston, VA, AIAA Educational Series, , 431 p. (ISBN 978-1-62410-171-7)
  12. (en) D. Burnett, « The Distribution of Molecular Velocities and the Mean Motion in a Non-Uniform Gas », Proceedings of the London Mathematical Society, vol. s2-40, no 1,‎ , p. 382-435
  13. (en) Byung Chan Eu, Nonequilibrium Statistical Mechanics. Ensemble Method., Springer, (ISBN 978-90-481-5007-6)
  14. Cédric Villani, Limites hydrodynamiques de l'équation de Boltzmann, coll. « Séminaire Bourbaki » (no 893 [2]),

Voir aussi

modifier

Articles connexes

modifier