A Shakhov-based Bhatnagar-Gross-Krook model for polyatomic molecules and for atomic as well as polyatomic mixtures

Cet article présente l'extension du modèle cinétique SBGK dans le code PICLas pour simuler des mélanges d'atomes et de molécules polyatomiques hors équilibre, en démontrant sa supériorité par rapport à la méthode ESBGK pour la capture précise des chocs dans des écoulements hypersoniques.

L'essentiel

🚀 Le "GPS" des gaz : Une nouvelle méthode pour prédire le comportement de l'air (et du vide)

Imaginez que vous essayez de prédire comment l'air se comporte autour d'une fusée qui file à une vitesse folle, ou comment les gaz s'écoulent dans un microchip tout petit. C'est un cauchemar pour les mathématiciens !

Dans l'espace ou dans les micro-machines, l'air ne se comporte pas toujours comme un fluide lisse (comme de l'eau dans une rivière). Parfois, il se comporte comme une foule de balles de ping-pong qui rebondissent les unes sur les autres de manière chaotique.

Ce papier présente une nouvelle "recette" (un modèle mathématique) appelée SBGK, créée par des chercheurs de l'Université de Stuttgart, pour simuler ces situations complexes beaucoup plus vite et plus précisément que les méthodes actuelles.

Voici comment ça marche, image par image :

1. Le Problème : Deux mondes qui ne se parlent pas

Pour simuler les gaz, on a deux outils principaux :

• DSMC (La méthode des balles) : C'est comme si vous preniez des millions de balles de ping-pong et que vous les laissiez rebondir une par une. C'est ultra-précis, mais c'est très lent et coûteux en énergie de calcul, surtout quand il y a beaucoup de balles (régime continu).
• CFD (La méthode du fluide) : C'est comme si vous traitiez l'air comme un liquide lisse. C'est très rapide, mais ça ne fonctionne pas bien quand les balles sont trop espacées (régime raréfié).

Le défi : Comment avoir la précision de la méthode "balles" avec la vitesse de la méthode "fluide" ? C'est là qu'intervient le modèle BGK (une version simplifiée des équations de la physique).

2. La Solution : Le modèle "Shakhov-BGK" (SBGK)

Les chercheurs ont amélioré un modèle existant (appelé BGK) pour qu'il soit plus intelligent.

• L'analogie du chef d'orchestre : Imaginez que le gaz est un orchestre.
  • Les atomes (comme l'hélium) sont des musiciens simples qui ne jouent que d'une note (le mouvement vers l'avant).
  • Les molécules (comme l'oxygène ou le CO2) sont des musiciens complexes : ils bougent, ils tournent sur eux-mêmes (rotation) et ils vibrent (comme une corde de guitare).
• L'ancien modèle (ESBGK) : C'était un chef d'orchestre un peu pressé. Il disait : "Allez, tout le monde joue la même chose !" Cela marchait bien, mais quand il y avait un choc violent (comme un mur de son), il avait du mal à voir exactement où commençait le chaos. Il "floutait" un peu le début du choc.
• Le nouveau modèle (SBGK) : C'est un chef d'orchestre très précis. Il sait exactement comment chaque type de musicien (atome ou molécule complexe) doit se comporter, même quand ils vibrent et tournent.

3. Ce que le papier a fait de nouveau

Les chercheurs ont pris ce modèle SBGK et l'ont adapté pour deux choses difficiles :

1. Les mélanges : Gérer des gaz composés de différents ingrédients (ex: un mélange d'azote et d'oxygène, ou d'hélium et d'argon).
2. Les molécules complexes : Gérer les molécules qui ont de l'énergie interne (vibration et rotation), ce qui est crucial pour les gaz comme le CO2 ou l'azote à haute température.

Ils ont aussi créé une "règle de mélange" pour calculer la viscosité (l'épaisseur du gaz) et la conductivité thermique (comment il conduit la chaleur) sans avoir à tout recalculer à chaque fois.

4. Les Tests : Le "Stress Test"

Pour vérifier si leur recette fonctionne, ils l'ont testée dans deux situations extrêmes :

• Test 1 : Le Couette Supersonique (Le tapis roulant géant)
  Imaginez deux plaques de métal avec du gaz entre elles. L'une bouge vite vers la droite, l'autre vers la gauche. C'est comme un tapis roulant qui fait tourner le gaz.

  • Résultat : Le nouveau modèle a prédit la température du gaz avec une précision quasi parfaite par rapport à la méthode lente (DSMC), même pour des mélanges d'atomes très différents (comme l'hélium et l'argon).

• Test 2 : Le Cône à 70° (L'entrée dans l'atmosphère)
  Ils ont simulé un objet pointu (un cône) traversant l'air à une vitesse hypersonique (comme une fusée qui rentre dans l'atmosphère). Cela crée un choc (une onde de pression très forte).

  • Le verdict : C'est ici que le nouveau modèle brille. Là où l'ancien modèle (ESBGK) voyait le choc un peu trop tôt et de manière floue, le nouveau modèle SBGK a vu le choc exactement là où il devrait être, avec une netteté incroyable. C'est comme passer d'une photo floue à une photo 4K.

5. Pourquoi c'est important ?

Ce travail est une étape clé pour le futur de l'aérospatiale et des technologies spatiales.

• Gain de temps : Ce modèle est beaucoup plus rapide que la méthode des "balles" (DSMC).
• Précision : Il est plus précis que les autres modèles rapides pour voir les chocs et les zones de turbulence.
• Applications : Cela aidera à concevoir de meilleurs satellites, des fusées pour Mars (où l'atmosphère est du CO2), et des micro-machines.

En résumé :
Les chercheurs ont créé un "simulateur de gaz" ultra-rapide et ultra-précis qui comprend la complexité des molécules qui vibrent et tournent. C'est comme donner à un ordinateur une carte routière détaillée pour naviguer dans le chaos des gaz, là où les autres cartes étaient un peu floues. Cela permettra de mieux concevoir nos voyages dans l'espace et nos technologies du futur.

Résumé technique

1. Problématique et Contexte

Les simulations numériques de la dynamique des fluides dans les applications spatiales, les écoulements micro/nanoscopiques et les technologies du vide présentent des défis majeurs. Ces écoulements impliquent souvent de grands gradients de densité, s'étendant du régime continu au régime d'écoulement moléculaire libre.

• Limites des méthodes existantes : La méthode de Monte Carlo par simulation directe (DSMC) est très précise pour les régimes non continus, mais son coût computationnel devient prohibitif dans les régimes de transition et continus. À l'inverse, la dynamique des fluides computationnelle (CFD) est efficace en régime continu mais peine à coupler avec la DSMC en raison de différences fondamentales (bruit statistique de la DSMC vs équations déterministes de la CFD).
• Solution potentielle : Les méthodes basées sur des particules pour les régimes continus, telles que les modèles Bhatnagar-Gross-Krook (BGK), offrent une alternative prometteuse. Elles approximent l'intégrale de collision de Boltzmann par un processus de relaxation, évitant ainsi la résolution explicite du libre parcours moyen et permettant des pas de temps plus grands.
• Le défi spécifique : Bien que les modèles BGK aient été étendus aux gaz monoatomiques et diatomiques, leur application aux mélanges de gaz polyatomiques (incluant atomes et molécules) avec des degrés de liberté internes non équilibrés (rotation, vibration) reste complexe. Les modèles existants (comme le modèle ESBGK) peuvent présenter des imprécisions dans la capture des chocs forts. L'objectif est donc d'étendre le modèle Shakhov-BGK (SBGK) pour gérer ces mélanges complexes tout en conservant une efficacité computationnelle supérieure à la DSMC.

2. Méthodologie

Les auteurs ont implémenté une extension du modèle SBGK dans le code open-source PICLas, en utilisant une approche stochastique basée sur des particules.

A. Formulation du Modèle

Le modèle repose sur l'approximation de l'opérateur BGK, où la fonction de distribution $f_s$ d'une espèce $s$ relaxe vers une distribution cible $f_s^t$ :
$$\partial_t f_s + \mathbf{v} \cdot \partial_x f_s = \nu (f_s^t - f_s)$$

• Distribution Cible pour les Mélanges :
  • Pour les atomes, la distribution cible ne contient que la composante translationnelle (Maxwellienne).
  • Pour les molécules, la distribution est factorisée en composantes translationnelle, rotationnelle et vibrationnelle, tenant compte des états quantifiés de vibration.
• Correction de Shakhov : Pour obtenir le bon nombre de Prandtl ($Pr$), la distribution cible translationnelle est modifiée par un terme correctif dépendant du flux de chaleur. Cela permet de corriger le flux thermique sans affecter le tenseur des contraintes visqueuses.
• Termes de Relaxation : Le modèle utilise un seul terme de relaxation par espèce. Cela simplifie considérablement le calcul par rapport aux modèles à multiples termes de relaxation, tout en respectant le principe d'indifférenciabilité (le modèle se réduit au cas mono-espèce si toutes les espèces sont identiques).

B. Propriétés Macroscopiques et Transport

• Conservation : Le modèle garantit la conservation de la masse, de la quantité de mouvement et de l'énergie totale.
• Nombre de Prandtl : Le paramètre de Shakhov ($Pr_S$) est dérivé analytiquement pour reproduire le nombre de Prandtl physique du mélange ($Pr_{phys}$) en fonction des fractions massiques et des degrés de liberté internes ($\xi_{int}$). La relation est donnée par $Pr_{phys} = Pr_S / \alpha$.
• Viscosité et Conductivité Thermique : Pour les mélanges, la viscosité et la conductivité thermique sont calculées en utilisant la première approximation des intégrales de collision (méthode de Chapman-Enskog) basée sur le modèle de sphère dure variable (VHS). Cette approche est jugée plus précise que la règle de mélange de Wilke pour les cas étudiés.
• Températures de Relaxation : Les températures effectives de relaxation (translationnelle, rotationnelle, vibrationnelle) sont déterminées pour suivre la loi de Landau-Teller, assurant un taux de relaxation correct pour les modes internes.

C. Implémentation

L'implémentation dans PICLas utilise une technique d'échantillonnage par acceptation-rejet pour générer les nouveaux états translationnels à partir de la distribution cible de Shakhov. Les énergies internes (rotation/vibration) sont gérées selon les modèles établis précédemment pour les gaz polyatomiques.

3. Résultats et Validation

Le modèle a été validé sur deux types d'écoulements : un écoulement de Couette supersonique et un écoulement hypersonique autour d'un cône émoussé à 70°. Les résultats sont comparés à la méthode de référence DSMC et au modèle ESBGK.

A. Écoulement de Couette Supersonique

Trois cas ont été testés :

1. N2 pur : Excellente concordance entre SBGK et DSMC pour les températures translationnelle et rotationnelle.
2. Mélange Ar-He (50/50) : Bonne concordance globale, mais une légère déviation observée. Les auteurs attribuent cela au grand rapport de masse entre l'argon et l'hélium, qui rend l'approximation du nombre de Prandtl effectif plus difficile.
3. Mélange N2-N (50/50) : Accord quasi-parfait. Le rapport de masse plus faible entre l'azote atomique et moléculaire améliore la précision par rapport au cas Ar-He.

B. Écoulement Hypersonique autour d'un Cône à 70°

Trois configurations gazeuses ont été simulées (N2 pur, mélange CO2-N2, mélange ternaire N2-O2-N) avec des conditions d'entrée hypersoniques.

• Structure du Choc : Le résultat le plus significatif concerne la capture de la structure du choc.
  • Le modèle ESBGK prédit une montée de la température translationnelle trop précoce, élargissant artificiellement la région du choc.
  • Le modèle SBGK proposé capture l'origine et l'épaisseur du choc avec une précision bien supérieure, correspondant très étroitement aux résultats DSMC.
• Températures Internes : Pour les mélanges polyatomiques (CO2-N2), le SBGK reproduit correctement les températures de vibration et de rotation, même dans les zones de non-équilibre fort derrière le choc.
• Flux de Chaleur : Les flux de chaleur calculés sur la surface du cône (côté amont et arrière) montrent un excellent accord entre SBGK, ESBGK, DSMC et les données expérimentales. Cependant, la précision accrue de SBGK dans la zone de choc se traduit par une meilleure fiabilité globale pour les prédictions de charges thermiques transitoires.

4. Contributions Clés

1. Extension du modèle SBGK : Développement d'un modèle BGK de type Shakhov capable de gérer simultanément des mélanges d'atomes et de molécules polyatomiques, incluant les états vibrationnels quantifiés.
2. Approche à terme unique : Utilisation d'un seul terme de relaxation par espèce pour les mélanges, offrant un compromis optimal entre précision physique (conservation, nombre de Prandtl correct) et coût computationnel réduit par rapport aux modèles multi-termes.
3. Précision supérieure dans les chocs : Démonstration que le modèle SBGK surpasse le modèle ESBGK dans la capture de la structure fine des chocs forts, un point critique pour les simulations d'entrée atmosphérique.
4. Validation rigoureuse : Comparaison exhaustive avec la DSMC et l'ESBGK sur des géométries complexes et des mélanges gazeux réalistes (y compris des scénarios d'entrée sur Mars simulés via CO2-N2).

5. Signification et Perspectives

Ce travail représente une avancée significative pour la simulation des écoulements gazeux non continus et en transition. En offrant une alternative plus précise et potentiellement moins coûteuse que la DSMC pour les régimes de Knudsen faibles à modérés, le modèle SBGK étendu facilite le couplage hybride DSMC-BGK pour des problèmes multi-échelles complexes (par exemple, l'entrée atmosphérique de satellites ou de véhicules spatiaux).

La capacité à traiter correctement les degrés de liberté internes non équilibrés dans des mélanges complexes ouvre la voie à des simulations plus réalistes de phénomènes aérothermiques. Les travaux futurs visent à intégrer les réactions chimiques dans ce cadre, ce qui rendrait le modèle complet pour la simulation de plasmas réactifs et d'écoulements hypersoniques à haute température.