Lattice Boltzmann model for non-ideal compressible fluid dynamics

Ce papier présente un nouveau modèle de Boltzmann sur réseau utilisant des réseaux de premiers voisins et des termes de correction de quasi-équilibre pour simuler des écoulements de fluides compressibles non idéaux avec une cohérence thermodynamique et une stabilité numérique, validant avec succès sa précision par des simulations quantitatives d'interactions choc-goutte jusqu'à des nombres de Mach de 1,47.

L'essentiel

Imaginez que vous essayez de simuler le comportement d'un fluide sur un ordinateur. Pendant longtemps, les ordinateurs ont été excellents pour simuler des fluides « idéaux » — comme l'eau s'écoulant doucement dans une rivière ou l'air se déplaçant lentement autour d'une aile. Ces fluides suivent des règles simples et prévisibles.

Mais que se passe-t-il lorsque le fluide est soumis à une pression et une chaleur extrêmes, se comportant comme un gaz dense presque liquide, ou un liquide presque gazeux ? C'est le monde des fluides compressibles non idéaux. Imaginez un fluide « stressé » qui agit bizarrement, refusant de suivre les règles simples du monde idéal. Cela se produit dans des technologies avancées comme les turbines au CO2 supercritique et les cycles énergétiques organiques.

Le problème est que les outils informatiques existants peinent avec ces fluides stressés. Ils plantent soit, donnent des réponses erronées, soit exigent une puissance de calcul si massive qu'ils deviennent inutiles.

Cet article présente une nouvelle méthode plus intelligente pour simuler ces fluides complexes, utilisant une méthode appelée méthode de Boltzmann sur réseau (LBM). Voici comment la nouvelle approche des auteurs fonctionne, expliquée par de simples analogies :

1. Le système à « deux voies »

La plupart des anciennes méthodes de simulation tentent de suivre tout (masse, vitesse, énergie) avec un seul ensemble compliqué de règles. Les auteurs ont réalisé que c'était comme essayer de conduire une voiture tout en jonglant simultanément avec une douzaine de balles — cela devient désordonné et instable.

Au lieu de cela, ils ont construit un système à deux voies :

• Voie A (La foule) : Un ensemble de règles suit la densité et la vitesse du fluide (combien de particules il y a et où elles vont).
• Voie B (L'énergie) : Un second ensemble de règles, séparé, suit l'énergie totale.
  En séparant ces éléments, l'ordinateur ne se perd pas. C'est comme avoir un contrôleur de circulation dédié pour les voitures et un autre séparé pour le carburant, garantissant qu'aucun système ne fait planter l'autre.

2. L'attracteur « quasi-équilibre »

En physique, les fluides veulent naturellement se stabiliser dans un état calme appelé « équilibre ». Cependant, dans ces conditions extrêmes, le fluide est constamment poussé et tiré, de sorte qu'il ne se stabilise jamais vraiment.

Les auteurs ont inventé une astuce ingénieuse appelée un « attracteur quasi-équilibre ».

• L'analogie : Imaginez un chien poursuivant une balle. La balle représente l'« état de calme parfait ». Le chien représente le fluide. Dans une situation normale, le chien court directement vers la balle.
• Le problème : Dans ce fluide extrême, la balle continue de s'éloigner ou de changer de forme. Si le chien poursuit la balle aveuglément, il pourrait tomber d'une falaise (la simulation devient instable).
• La solution : Les auteurs ont donné au chien un « GPS » qui prédit où la balle sera une fraction de seconde plus tard, en fonction de la façon dont le vent (pression) et le terrain (densité) changent. Cette cible « décalée » permet au chien de courir sans heurter la falaise. Cela garantit que la simulation reste stable même lorsque le fluide se déplace très vite ou que sa densité change rapidement.

3. Réparer la chaleur « parasite »

Lorsque les fluides se déplacent rapidement, ils génèrent de la chaleur. Dans les modèles informatiques standards, la chaleur circule parfois dans la mauvaise direction ou crée une chaleur « fantôme » fictive qui n'existe pas dans la réalité.

• L'analogie : C'est comme un thermostat qui pense que la pièce est gelée parce qu'il mesure le courant d'air d'une fenêtre, et non la température réelle de la pièce.
• La réparation : Les auteurs ont ajouté un « terme de correction » spécifique à leurs équations. Cela agit comme un filtre qui élimine les courants d'air factices, garantissant que la chaleur circule exactement comme le dictent les lois de la physique (loi de Fourier), même dans ces conditions extrêmes et non idéales.

4. Les tests « tube à choc » et « colonne liquide »

Pour prouver que leur nouvelle méthode fonctionne, ils n'ont pas seulement fait des mathématiques ; ils ont réalisé des tests extrêmes :

• Le tube à choc : Ils ont simulé une explosion soudaine de pression (une onde de choc) se propageant dans un gaz dense. Dans les gaz normaux, ces ondes se comportent d'une certaine manière. Dans ces gaz « non idéaux », les ondes peuvent faire quelque chose d'étrange : un « choc de détente » (une onde qui s'étale mais agit toujours comme un choc net). Leur modèle a prédit avec succès ce comportement étrange, que les anciens modèles avaient manqué.
• La colonne liquide : Ils ont simulé une onde de choc à haute vitesse frappant une goutte de liquide. C'est un test très difficile car le choc rebondit, se réfléchit et déchire la goutte. Leur modèle a géré la collision parfaitement, correspondant aux expériences réelles où la goutte de liquide s'aplatit et se dilate exactement comme elle le devrait.

Pourquoi cela compte

Les auteurs affirment que leur méthode est rapide, stable et précise. Elle utilise une grille simple (comme un échiquier standard) plutôt que d'avoir besoin d'une grille super-complexe et étirée. Cela signifie que les scientifiques peuvent maintenant simuler ces écoulements de fluides extrêmes, à haute vitesse et non idéaux, sur des ordinateurs standards avec une grande précision.

En résumé : L'article présente un nouveau « manuel du conducteur » pour les simulations informatiques qui leur permet de gérer des fluides sous stress extrême sans planter. En utilisant un système à deux voies et un « GPS » intelligent pour l'énergie du fluide, ils peuvent prédire avec précision le comportement de ces fluides complexes dans des scénarios à haute vitesse, ouvrant la voie à de meilleures conceptions pour les systèmes énergétiques avancés.

Résumé technique

Résumé technique : Modèle de Boltzmann sur réseau pour la dynamique des fluides compressibles non idéaux

Énoncé du problème
La dynamique des fluides compressibles non idéaux (NICFD) traite des comportements des fluides dans les régimes proche-critique, transcritique et supercritique, où les fluides s'écartent considérablement des lois des gaz parfaits. Ces régimes sont critiques pour les technologies énergétiques modernes, telles que les cycles de Rankine organiques et les turbines à CO₂ supercritique. Cependant, les données expérimentales dans ces états sont rares et difficiles à acquérir. Bien que la simulation numérique directe (DNS) soit nécessaire pour comprendre la physique complexe (par exemple, les structures de choc, la turbulence dans les vapeurs denses) et générer des bases de données d'ingénierie, les outils numériques existants peinent souvent à assurer la cohérence thermodynamique, la stabilité et la récupération précise des limites hydrodynamiques pour les équations d'état (EOS) non idéales. Plus précisément, les méthodes classiques de dynamique des fluides computationnelle (CFD) rencontrent des difficultés dans la région spinodale thermodynamiquement instable où le carré de la vitesse du son adiabatique devient négatif, nécessitant des traitements spéciaux comme la règle des aires égales de Maxwell. La méthode de Boltzmann sur réseau (LBM), bien que réussie pour les écoulements incompressibles et multiphasiques, a historiquement négligé la thermodynamique compressible non idéale dans ses formulations cinétiques, s'appuyant souvent sur une régularisation ad hoc ou des maillages étendus qui compromettent l'efficacité ou l'invariance galiléenne.

Méthodologie
Les auteurs proposent un modèle cinétique novateur et sa réalisation par méthode de Boltzmann sur réseau, conçus pour retrouver les équations compressibles de Navier–Stokes–Korteweg (NSK) pour les fluides non idéaux. La méthodologie repose sur les piliers suivants :

1. Approche à double fonction de distribution (DDF) : Le modèle emploie deux fonctions de distribution : $f$ pour la masse et la quantité de mouvement, et $g$ pour l'énergie volumique ($\rho E$). Ce découplage permet l'utilisation de réseaux classiques de premiers voisins (spécifiquement D3Q27) sans nécessiter de maillages étendus.
2. Attracteurs de quasi-équilibre : L'opérateur de collision utilise une structure de type BGK à double relaxation impliquant un attracteur d'équilibre local ($f^{eq}, g^{eq}$) et un attracteur de quasi-équilibre décalé ($f^\star_\lambda, g^\star_\lambda$).
   • L'équilibre local est défini en utilisant un paramètre de température de référence $\theta = P/\rho$ (travail thermodynamique de l'écoulement) plutôt que la température thermodynamique $T$, assurant ainsi la largeur correcte de la distribution de Maxwell indépendante de l'énergie interne.
   • L'équilibre décalé introduit des corrections pour la force de Korteweg (tension superficielle), la viscosité volumique et le flux de chaleur. Il utilise des valeurs décalées pour la vitesse d'écoulement ($u^\star$), la température de référence ($\theta^\star$) et la température thermodynamique ($T^\star$).
3. Cohérence thermodynamique :
   • La température de référence $\theta$ est fixée à $P/\rho$, une formulation qui garantit que le modèle reste valable pour les EOS non idéales (par exemple, van der Waals).
   • Un terme de correction dans le flux de chaleur ($q^c_\lambda$) est introduit pour compenser les composantes non fourierennes spurieuses résultant des gradients de densité dans les fluides non idéaux, assurant ainsi la récupération de la loi de Fourier.
   • Le décalage de la température de référence permet à la viscosité volumique d'être un paramètre indépendant, réglable et défini positif, contournant ainsi le problème de la viscosité volumique négative qui peut survenir dans les modèles BGK standards pour les pressions non idéales.
4. Discrétisation : Le modèle utilise un schéma explicite d'ordre deux dérivé de l'intégration trapézoïdale le long des caractéristiques. Il emploie un réseau standard D3Q27 et traite les termes sources (force de Korteweg) via l'équilibre décalé, évitant ainsi les erreurs de discrétisation des termes sources explicites.

Contributions clés

• Nouveau cadre cinétique : L'article introduit un modèle cinétique qui retrouve les équations hydrodynamiques cibles (bilans de masse, de quantité de mouvement et d'énergie volumique) pour les fluides compressibles non idéaux en utilisant un maillage cinétique minimal (premiers voisins).
• Dissipation définie positive : En construisant le modèle via des attracteurs de quasi-équilibre, les auteurs garantissent que les taux de dissipation de Navier–Stokes (viscosité de cisaillement et viscosité volumique) sont définis positifs et invariants galiléens, un point de défaillance courant dans les tentatives précédentes de LBM pour les écoulements compressibles non idéaux.
• Robustesse thermodynamique : Le modèle gère naturellement les régions d'instabilité thermodynamique (spinodale) où la vitesse du son devient imaginaire dans la limite non visqueuse, car la LBM hérite de la propagation le long de caractéristiques discrètes fixes plutôt que de compter sur des solveurs d'EDP hyperboliques nécessitant un traitement spécial dans ces régions.
• Validation des phénomènes non classiques : L'étude fournit la première validation quantitative de la LBM pour les interactions choc-goutte à des nombres de Mach allant jusqu'à 1,47 dans un contexte de fluide non idéal, démontrant la capacité du modèle à capturer des structures d'ondes complexes.

Résultats
Le modèle a été validé à travers une hiérarchie de configurations de complexité croissante :

1. Stabilité linéaire et dispersion : Les simulations de la vitesse du son, de la dissipation des ondes de cisaillement et de la conductivité thermique ont montré un excellent accord avec les solutions analytiques pour un fluide de van der Waals ajusté aux propriétés de l'azote sur une gamme de nombres de Mach et de températures réduites.
2. Cohérence multiphasique : Le modèle a reproduit avec précision les densités de coexistence liquide-vapeur et les profils d'interface, démontrant une convergence d'ordre deux avec le raffinement du maillage et un accord avec la règle des aires égales de Maxwell.
3. Tubes à choc non classiques : Trois cas de tubes à choc 1-D ont été simulés pour sonder la dynamique des ondes non classiques caractérisée par la dérivée fondamentale de la dynamique des gaz ($\Gamma$).
   • Les cas avec $\Gamma < 0$ ont correctement capturé la formation de chocs de détente (ondes de détente se raidissant) et de ventilateurs de compression (ondes de compression s'étalant), des phénomènes distincts de la dynamique des gaz classique.
   • Le modèle a réussi à transitionner entre les régimes classique ($\Gamma > 0$) et non classique ($\Gamma < 0$) au sein d'une seule simulation.
4. Interaction choc-colonne de liquide : Une simulation 2D d'une onde de choc planaire interagissant avec une colonne de liquide saturé près du point critique ($M_a = 1,47$) a été réalisée.
   • Le modèle a capturé des structures d'ondes complexes incluant des chocs transmis, des chocs réfléchis, des tiges de Mach et des points triples.
   • La comparaison quantitative de l'évolution de la largeur de la colonne de liquide par rapport aux données expérimentales et aux résultats numériques précédents a montré un bon accord.
   • L'analyse thermodynamique a révélé que l'interaction du choc pousse la vapeur environnante vers un état supercritique tandis que la colonne de liquide reste marginalement sous-critique, capturant avec précision le fort couplage thermodynamique.

Signification et revendications
Les auteurs affirment que ce travail étend l'applicabilité des méthodes de Boltzmann sur réseau aux écoulements compressibles non idéaux à haute vitesse avec un maillage cinétique minimal. La signification principale réside dans la fourniture d'un cadre thermodynamiquement cohérent et numériquement stable qui ne dépend ni de maillages étendus ni de régularisation ad hoc. Le modèle est présenté comme une théorie cinétique réduite ciblant la limite hydrodynamique NSK, capable de gérer toute la gamme des comportements non idéaux, y compris la région spinodale, sans les instabilités numériques souvent associées à la CFD conventionnelle dans ces régimes. La simulation réussie des interactions choc-goutte à des nombres de Mach élevés sert de validation rigoureuse, démontrant le potentiel de la méthode pour investiguer des écoulements complexes dans les technologies supercritiques et les processus pilotés par les transitions de phase.