Modèle de Boltzmann sur réseau pour la dynamique des fluides compressibles non-idéaux

Informations fondamentales

• ID de l'article: 2510.14712
• Titre: Lattice Boltzmann model for non-ideal compressible fluid dynamics
• Auteurs: S. A. Hosseini, M. Feinberg, I. V. Karlin (ETH Zurich)
• Classification: physics.flu-dyn
• Journal de publication: Journal of Fluid Mechanics
• Lien de l'article: https://arxiv.org/abs/2510.14712

Résumé

Cet article propose une méthode de Boltzmann sur réseau pour simuler l'écoulement de fluides compressibles non-idéaux. La méthode utilise un réseau de premiers voisins et introduit un ensemble cohérent de termes correctifs par l'intermédiaire d'un attracteur quasi-équilibre, garantissant un taux de dissipation de Navier-Stokes défini positif et invariant galiléen. Cette construction évite le recours à des gabarits étendus ou à une régularisation temporaire, tout en maintenant la stabilité numérique et la cohérence thermodynamique sur une large gamme d'états d'écoulement. Le modèle résultant reproduit fidèlement la mécanique des fluides aux niveaux Euler et Navier-Stokes. À titre de validation rigoureuse, les auteurs démontrent pour la première fois dans le cadre de Boltzmann sur réseau une simulation quantitativement précise de l'interaction gouttelette-onde de choc jusqu'à un nombre de Mach de 1,47.

Contexte et motivation de la recherche

Définition du problème

La dynamique des fluides compressibles non-idéaux est une branche émergente et en rapide développement de la mécanique des fluides, principalement due à l'émergence de méthodes et de technologies fonctionnant à des états proches du point critique, transcritiques et supercritiques. Cette classe de fluides comprend:

1. Fluides supercritiques: fluides aux conditions de pression et température au-dessus du point critique
2. Vapeurs denses: gaz dont le facteur de compressibilité s'écarte du comportement du gaz idéal
3. Fluides BZT: fluides complexes de haut poids moléculaire présentant une dérivée négative de la dynamique gazeuse fondamentale

Importance

1. Besoins des applications industrielles: cycles de Rankine organique, turbines CO₂ supercritiques et autres technologies énergétiques
2. Rareté des données expérimentales: configurations expérimentales complexes et difficultés d'acquisition de données
3. Besoin d'outils numériques: nécessité d'outils numériques cohérents et efficaces pour comprendre les phénomènes physiques complexes

Limitations des méthodes existantes

1. Méthodes CFD traditionnelles: dans la région de la spinodale thermodynamiquement instable, le carré de la vitesse du son adiabatique devient négatif, causant la transition des équations d'évolution du type hyperbolique au type parabolique
2. Modèles LBM existants: principalement orientés vers les applications d'ébullition, manquant de modèles systématiques pour les écoulements compressibles non-idéaux
3. Problèmes de stabilité numérique: défis de stabilité à haut nombre de Mach et sous des effets non-idéaux forts

Contributions principales

1. Proposition d'un nouveau modèle LBM à double fonction de distribution: utilisant le réseau classique de premiers voisins (D3Q27) pour simuler les écoulements compressibles non-idéaux
2. Introduction du mécanisme d'attracteur quasi-équilibre: assurant via les termes correctifs fᵢ et gᵢ un taux de dissipation défini positif et invariant galiléen
3. Réalisation d'un contrôle indépendant de la viscosité de volume: résolvant le problème de viscosité de volume potentiellement négative dans la structure BGK
4. Première simulation d'interaction gouttelette-onde de choc à haut nombre de Mach: réalisant une simulation quantitativement précise à Ma=1,47 dans le cadre LBM
5. Maintien de la cohérence thermodynamique: préservant la stabilité numérique sur une large gamme d'états d'écoulement

Détails de la méthode

Définition de la tâche

Simulation de fluides compressibles non-idéaux satisfaisant le système d'équations de Navier-Stokes-Korteweg (NSK):

• Conservation de la masse: ∂ₜρ + ∇·(ρu) = 0
• Conservation de la quantité de mouvement: ∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
• Conservation de l'énergie: ∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

où T^K est le tenseur de tension superficielle de Korteweg et T^NS est le tenseur de contrainte visqueuse.

Architecture du modèle

1. Conception à double fonction de distribution

Utilisation de deux ensembles de fonctions de distribution {fᵢ, gᵢ}, i=1,...,27:

• fᵢ: traite la conservation de la masse et de la quantité de mouvement
• gᵢ: traite la conservation de l'énergie

2. Équation de Boltzmann sur réseau

fᵢ(x+cᵢδt, t+δt) = fᵢ + 2β(fᵢᵉᵠ - fᵢ) + (1-β)(fᵢ* - fᵢᵉᵠ)
gᵢ(x+cᵢδt, t+δt) = gᵢ + 2β(gᵢᵉᵠ - gᵢ) + (1-β)(gᵢ* - gᵢᵉᵠ)

où β = Pδt/(2μ + Pδt) est le paramètre de relaxation.

3. Fonctions de distribution d'équilibre

Adoption d'une forme produit:

fᵢᵉᵠ = ρ ∏_α Ψᵢα(uα, P/ρ + uα²)

où Ψᵢα est la fonction de base:

Ψᵢα(ξα, ζαα) = 1 - cᵢα² + ½[(3cᵢα² - 2)ζαα + cᵢαξα]

4. Termes correctifs quasi-équilibre

Correction de la fonction de distribution f:

• Terme de force: Fα = -κρ∂α∇²ρ (force de Korteweg)
• Correction de contrainte: Φαα incluant la correction du tenseur de contrainte diagonale et le contrôle indépendant de la viscosité de volume

Correction de la fonction de distribution g:

• Vitesse décalée: u* = u + δt F/ρ
• Température décalée: T* = T - δt²F·F/(2ρ²cᵥ)
• Correction du flux thermique: considérant les effets de gradient d'enthalpie du gaz non-idéal

Points d'innovation technique

1. Mécanisme d'attracteur quasi-équilibre

Par l'introduction de fᵢ et gᵢ, le système ne se détend pas simplement vers l'équilibre local, mais vers un état quasi-équilibre corrigé, ce qui garantit:

• Une dynamique macroscopique correcte
• La stabilité numérique
• L'invariance galiléenne

2. Contrôle indépendant de la viscosité de volume

Dans le modèle BGK traditionnel, la viscosité de volume η' = ((D+2)/D - ρcs²/P)μ peut être négative. Un contrôle indépendant est réalisé via le terme correctif Φ':

Φ' = P((D+2)/D - ρcs²/P - η/μ)(∇·u)

3. Cohérence thermodynamique

Par une conception soignée de la température décalée T* et des corrections d'énergie, la cohérence thermodynamique est assurée sous des équations d'état non-linéaires.

Configuration expérimentale

Cas de validation

1. Validation fondamentale:
   • Mesure de la vitesse du son
   • Vérification de la viscosité de cisaillement
   • Écoulement de Couette thermique
   • Amortissement du mode normal
2. Validation des écoulements multiphasiques:
   • Densités de coexistence liquide-vapeur
   • Cohérence et convergence de l'interface
3. Validation des écoulements compressibles:
   • Tube de choc unidimensionnel non-idéal
   • Interaction onde de choc-colonne liquide

Paramètres de substance

Principalement l'azote (N₂) pour la validation:

• Pression critique: Pc = 3,4×10⁶ Pa
• Densité critique: ρc = 241,96 kg/m³
• Température critique: Tc = 126,2 K
• R/cv = 0,4

Équation d'état

Adoption de l'équation d'état de van der Waals:

P(ρ,T) = ρRT/(1-bρ) - aρ²

Résultats expérimentaux

Résultats principaux

1. Vérification des grandeurs physiques fondamentales

• Vitesse du son: accord parfait avec les solutions analytiques sur les branches liquide et vapeur saturées dans la plage Tr∈0,7,1
• Viscosité de cisaillement: maintien de l'invariance galiléenne à différents nombres de Mach, erreur <1%
• Écoulement de Couette thermique: accord avec les solutions analytiques pour Pr∈{0,6,1,2,4,9} et Ma∈{0,8,1,2,1,6}

2. Propriétés des écoulements multiphasiques

• Densités de coexistence: accord complet avec la construction d'égale surface de Maxwell dans la plage Tr∈0,3,0,99
• Convergence de l'interface: démonstration de bonnes caractéristiques de convergence à différentes résolutions de maillage

3. Tube de choc non-idéal

Simulation réussie de trois configurations différentes de tube de choc non-idéal:

• Cas I: onde de raréfaction dans la région Γ<0
• Cas II: comportement d'onde de compression classique
• Cas III: onde de raréfaction mixte traversant dynamiquement la ligne Γ=0

Tous les cas montrent un excellent accord avec les données de référence.

4. Interaction onde de choc-gouttelette

• Première simulation à haut nombre de Mach (Ma=1,47) réalisée dans le cadre LBM
• Capture précise des structures d'onde complexes: onde de choc incidente, onde transmise, onde de choc réfléchie, tige de Mach
• Évolution de la déformation de la gouttelette en accord étroit avec les données expérimentales

Expériences d'ablation

Vérification de la convergence du modèle par des tests d'interface à différentes résolutions de maillage:

• Quatre niveaux de résolution avec δx de 5 μm à 0,1 μm
• Démonstration de bonnes caractéristiques de convergence de maillage

Découvertes clés

1. Stabilité numérique: maintien de la stabilité à haut nombre de Mach et sous effets non-idéaux forts
2. Efficacité computationnelle: utilisation d'un réseau standard de premiers voisins, coût de calcul maîtrisé
3. Fidélité physique: reproduction précise des phénomènes physiques complexes des fluides non-idéaux

Travaux connexes

Modèles LBM multiphasiques traditionnels

1. Modèle de gradient de couleur (Gunstensen et al. 1991)
2. Modèle de pseudo-potentiel (Shan & Chen 1993)
3. Modèle d'énergie libre (Swift et al. 1996)

Développement du LBM compressible

• He & Doolen (2002): première tentative d'écoulement compressible non-idéal
• Les modèles existants sont principalement orientés vers les applications d'ébullition
• Absence de modèles systématiques pour les écoulements compressibles non-idéaux

Avantages relatifs de cet article

1. Approche systématique: cadre complet à double fonction de distribution
2. Fondements théoriques: analyse multi-échelle rigoureuse
3. Complétude de la validation: validation complète des grandeurs physiques fondamentales aux phénomènes d'écoulement complexes

Conclusions et discussion

Conclusions principales

1. Développement réussi d'un modèle LBM applicable aux fluides compressibles non-idéaux
2. Réalisation de la stabilité numérique et de la cohérence physique via le mécanisme d'attracteur quasi-équilibre
3. Première simulation d'interaction gouttelette-onde de choc à haut nombre de Mach dans le cadre LBM
4. Le modèle maintient la cohérence thermodynamique sur une large gamme d'états d'écoulement

Limitations

1. Modèle à temps de relaxation unique: utilisation actuelle du simple modèle de collision BGK, potentiellement insuffisant pour la simulation de turbulence
2. Restrictions sur l'équation d'état: validation principalement sur l'équation de van der Waals, nécessitant une vérification supplémentaire pour d'autres équations d'état complexes
3. Validation tridimensionnelle limitée: la plupart des validations concentrées sur des configurations unidimensionnelles et bidimensionnelles

Directions futures

1. Extension à des modèles de collision plus robustes: pour la simulation de configurations turbulentes
2. Équations d'état plus complexes: telles que Peng-Robinson et autres
3. Applications tridimensionnelles à grande échelle: application à des problèmes d'ingénierie réels

Évaluation approfondie

Avantages

1. Rigueur théorique: analyse multi-échelle complète et dérivations mathématiques
2. Innovativité de la méthode: le mécanisme d'attracteur quasi-équilibre constitue une innovation importante
3. Complétude de la validation: validation systématique de la physique fondamentale aux applications complexes
4. Valeur pratique: résolution de problèmes d'applications d'ingénierie importants
5. Efficacité numérique: utilisation d'un réseau standard évitant la complexité computationnelle

Insuffisances

1. Modèle de collision simple: applicabilité limitée du modèle BGK aux écoulements complexes
2. Dépendance paramétrique: absence de directives systématiques pour le choix de certains paramètres correctifs
3. Validation tridimensionnelle insuffisante: nécessité de validations supplémentaires pour des géométries complexes tridimensionnelles

Impact

1. Contribution académique: ouverture d'une nouvelle direction pour le LBM dans le domaine des écoulements compressibles non-idéaux
2. Valeur d'ingénierie: fourniture d'un outil numérique important pour les technologies de fluides supercritiques
3. Signification méthodologique: le concept d'attracteur quasi-équilibre peut influencer d'autres méthodes numériques

Scénarios d'application

1. Technologie des fluides supercritiques: turbines CO₂, cycles de Rankine organique
2. Écoulements multiphasiques: rupture de gouttelettes, processus de changement de phase
3. Dynamique des ondes de choc: interactions d'ondes complexes dans les écoulements à haute vitesse
4. Simulation numérique directe: recherche physique fondamentale nécessitant une haute précision

Références

1. Guardone, A. et al. (2024). Nonideal compressible fluid dynamics of dense vapors and supercritical fluids. Annual Review of Fluid Mechanics, 56, 241-269.
2. Hosseini, S.A. & Karlin, I.V. (2023). Lattice Boltzmann for non-ideal fluids: Fundamentals and practice. Physics Reports, 1030, 1-137.
3. Shan, X. & Chen, H. (1993). Lattice Boltzmann model for simulating flows with multiple phases and components. Physical Review E, 47(3), 1815.