Lattice-Boltzmann-Modell für nicht-ideale kompressible Fluiddynamik

Grundinformationen

• Paper-ID: 2510.14712
• Titel: Lattice Boltzmann model for non-ideal compressible fluid dynamics
• Autoren: S. A. Hosseini, M. Feinberg, I. V. Karlin (ETH Zürich)
• Klassifizierung: physics.flu-dyn
• Veröffentlichtes Journal: Journal of Fluid Mechanics
• Paper-Link: https://arxiv.org/abs/2510.14712

Zusammenfassung

In diesem Artikel wird eine Lattice-Boltzmann-Methode zur Simulation von kompressiblen nicht-idealen Fluidströmungen vorgestellt. Das Verfahren nutzt ein Gitter mit ersten Nachbarn und führt durch einen quasi-Gleichgewichts-Attraktor einen konsistenten Satz von Korrekturgliedern ein, um eine positive Definitheit und Galilei-Invarianz der Navier-Stokes-Dissipationsrate zu gewährleisten. Diese Konstruktion vermeidet die Notwendigkeit erweiterter Schablonen oder temporärer Regularisierung und erhält gleichzeitig numerische Stabilität und thermodynamische Konsistenz über einen breiten Bereich von Strömungszuständen. Das resultierende Modell reproduziert die Fluiddynamik auf der Euler- und Navier-Stokes-Ebene präzise. Als strenge Validierung präsentieren die Autoren erstmals quantitativ genaue Simulationen von Tropfen-Stoßwellen-Wechselwirkungen mit Machzahlen bis zu 1,47 im Lattice-Boltzmann-Rahmen.

Forschungshintergrund und Motivation

Problemdefinition

Die Dynamik nicht-idealer kompressibler Fluide ist ein aufstrebendes und sich schnell entwickelndes Gebiet der Fluiddynamik, hauptsächlich aufgrund des Auftretens von Methoden und Technologien, die unter nahe-kritischen, transkritischen und überkritischen Bedingungen arbeiten. Diese Fluidklassen umfassen:

1. Überkritische Fluide: Fluide unter Druck- und Temperaturbedingungen oberhalb des kritischen Punktes
2. Dichte Dämpfe: Gase, deren Kompressibilitätsfaktor vom Verhalten idealer Gase abweicht
3. BZT-Fluide: Hochmolekulare komplexe Fluide mit negativer fundamentaler gasdynamischer Ableitung

Bedeutung

1. Anforderungen industrieller Anwendungen: Organische Rankine-Zyklen, überkritische CO₂-Turbinen und andere Energietechnologien
2. Mangel an experimentellen Daten: Experimentelle Aufbauten sind komplex und die Datenbeschaffung schwierig
3. Bedarf an numerischen Werkzeugen: Konsistente und effiziente numerische Werkzeuge sind erforderlich, um komplexe physikalische Phänomene zu verstehen

Einschränkungen bestehender Methoden

1. Traditionelle CFD-Methoden: In thermodynamisch instabilen Spinodalregionen wird das Quadrat der adiabatischen Schallgeschwindigkeit negativ, was dazu führt, dass die Evolutionsgleichungen von hyperbolisch zu parabolisch werden
2. Bestehende LBM-Modelle: Konzentrieren sich hauptsächlich auf Siedungsanwendungen und ermangeln systematischer Modelle für kompressible nicht-ideale Strömungen
3. Numerische Stabilitätsprobleme: Stabilitätschallenges bei hohen Machzahlen und starken nicht-idealen Effekten

Kernbeiträge

1. Vorschlag eines neuen Dual-Verteilungs-LBM-Modells: Verwendung eines klassischen Gitters mit ersten Nachbarn (D3Q27) zur Simulation nicht-idealer kompressibler Strömungen
2. Einführung eines quasi-Gleichgewichts-Attraktor-Mechanismus: Durch Korrekturkenngrößen fᵢ und gᵢ wird positive Definitheit und Galilei-invariante Dissipationsrate gewährleistet
3. Realisierung unabhängiger Volumenviskositätskontrolle: Lösung des Problems der möglicherweise negativen Volumenviskosität in der BGK-Struktur
4. Erstmalige Realisierung von Tropfen-Stoßwellen-Wechselwirkungssimulationen bei hoher Machzahl: Quantitativ genaue Simulation bei Ma=1,47 im LBM-Rahmen
5. Erhaltung thermodynamischer Konsistenz: Aufrechterhaltung numerischer Stabilität über einen breiten Bereich von Strömungszuständen

Methodische Details

Aufgabendefinition

Simulation von kompressiblen nicht-idealen Fluiden, die die Navier-Stokes-Korteweg (NSK)-Gleichungen erfüllen:

• Massenerhaltung: ∂ₜρ + ∇·(ρu) = 0
• Impulserhaltung: ∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
• Energieerhaltung: ∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

wobei T^K der Korteweg-Oberflächenspannungstensor und T^NS der viskose Spannungstensor ist.

Modellarchitektur

1. Dual-Verteilungs-Design

Verwendung zweier Verteilungsfunktionssätze {fᵢ, gᵢ}, i=1,...,27:

• fᵢ: Behandlung von Massen- und Impulserhaltung
• gᵢ: Behandlung von Energieerhaltung

2. Lattice-Boltzmann-Gleichung

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ᵢᵉᵠ)

wobei β = Pδt/(2μ + Pδt) der Relaxationsparameter ist.

3. Gleichgewichtsverteilungsfunktionen

Verwendung einer Produktform:

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

wobei Ψᵢα die Basisfunktion ist:

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

4. Quasi-Gleichgewichts-Korrekturkenngrößen

Korrektur der f-Verteilungsfunktion:

• Kraftterm: Fα = -κρ∂α∇²ρ (Korteweg-Kraft)
• Spannungskorrektur: Φαα enthält Korrektionen des diagonalen Spannungstensors und unabhängige Volumenviskositätskontrolle

Korrektur der g-Verteilungsfunktion:

• Verschobene Geschwindigkeit: u* = u + δt F/ρ
• Verschobene Temperatur: T* = T - δt²F·F/(2ρ²cᵥ)
• Wärmestromkorrektur: Berücksichtigung von Enthalpie-Gradienten nicht-idealer Gase

Technische Innovationen

1. Quasi-Gleichgewichts-Attraktor-Mechanismus

Durch Einführung von fᵢ und gᵢ relaxiert das System nicht einfach zum lokalen Gleichgewicht, sondern zu einem modifizierten quasi-Gleichgewichtszustand, was gewährleistet:

• Korrekte makroskopische Dynamik
• Numerische Stabilität
• Galilei-Invarianz

2. Unabhängige Volumenviskositätskontrolle

In traditionellen BGK-Modellen kann die Volumenviskosität η' = ((D+2)/D - ρcs²/P)μ negativ sein. Durch den Korrekturterm Φ' wird unabhängige Kontrolle realisiert:

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

3. Thermodynamische Konsistenz

Durch sorgfältig gestaltete verschobene Temperatur T* und Energiekorrektionen wird thermodynamische Konsistenz unter nichtlinearen Zustandsgleichungen gewährleistet.

Experimentelle Einrichtung

Validierungsfälle

1. Grundlegende Validierung:
   • Schallgeschwindigkeitsmessung
   • Verifikation der Scherenviskosität
   • Thermisches Couette-Strömung
   • Dämpfung des Normalenmodus
2. Validierung mehrphasiger Strömungen:
   • Flüssig-Dampf-Koexistenzdichte
   • Oberflächenkonsistenz und Konvergenz
3. Validierung kompressibler Strömungen:
   • Eindimensionale nicht-ideale Stoßrohre
   • Stoßwellen-Flüssigkeitssäulen-Wechselwirkung

Stoffparameter

Hauptsächlich Stickstoff (N₂) zur Validierung:

• Kritischer Druck: Pc = 3,4×10⁶ Pa
• Kritische Dichte: ρc = 241,96 kg/m³
• Kritische Temperatur: Tc = 126,2 K
• R/cv = 0,4

Zustandsgleichung

Verwendung der van-der-Waals-Zustandsgleichung:

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

Experimentelle Ergebnisse

Hauptergebnisse

1. Validierung grundlegender physikalischer Größen

• Schallgeschwindigkeit: Perfekte Übereinstimmung mit analytischen Lösungen auf gesättigten Flüssigkeits- und Dampfzweigen im Bereich Tr∈0,7,1
• Scherenviskosität: Erhaltung der Galilei-Invarianz bei verschiedenen Machzahlen mit Fehler <1%
• Thermisches Couette-Strömung: Übereinstimmung mit analytischen Lösungen bei Pr∈{0,6,1,2,4,9} und Ma∈{0,8,1,2,1,6}

2. Mehrphasige Strömungseigenschaften

• Koexistenzdichte: Vollständige Übereinstimmung mit Maxwell-Flächenkonstruktion im Bereich Tr∈0,3,0,99
• Oberflächenkonvergenz: Gute Konvergenzeigenschaften bei verschiedenen Gitterauflösungen

3. Nicht-ideale Stoßrohre

Erfolgreiche Simulation von drei verschiedenen nicht-idealen Stoßrohr-Konfigurationen:

• Fall I: Verdünnungsstoßwelle in der Γ<0-Region
• Fall II: Klassisches Kompressionswellenverhalten
• Fall III: Gemischte Verdünnungswelle, die dynamisch die Γ=0-Linie überquert

Alle Fälle zeigen ausgezeichnete Übereinstimmung mit Referenzdaten.

4. Tropfen-Stoßwellen-Wechselwirkung

• Erstmalige Realisierung von Hochmachzahl-Simulation (Ma=1,47) im LBM-Rahmen
• Genaue Erfassung komplexer Wellenstrukturen: einfallende Stoßwelle, transmittierte Welle, reflektierte Stoßwelle, Mach-Stab
• Tropfendeformationsentwicklung zeigt hohe Übereinstimmung mit experimentellen Daten

Ablationsexperimente

Verifikation der Modellkonvergenz durch Oberflächentests mit verschiedenen Gitterauflösungen:

• Vier Auflösungsstufen von δx = 5 μm bis 0,1 μm
• Demonstration guter Gitterkonvergenzeigenschaften

Schlüsselfunde

1. Numerische Stabilität: Erhaltung der Stabilität bei hohen Machzahlen und starken nicht-idealen Effekten
2. Rechnerische Effizienz: Verwendung von Standard-Gittern mit ersten Nachbarn ermöglicht kontrollierbare Rechenkosten
3. Physikalische Genauigkeit: Präzise Reproduktion komplexer physikalischer Phänomene nicht-idealer Fluide

Verwandte Arbeiten

Traditionelle LBM-Mehrphasen-Modelle

1. Farbgradienten-Modell (Gunstensen et al. 1991)
2. Pseudopotential-Modell (Shan & Chen 1993)
3. Freie-Energie-Modell (Swift et al. 1996)

Entwicklung kompressibler LBM

• He & Doolen (2002): Erste Versuche bei kompressiblen nicht-idealen Strömungen
• Bestehende Modelle konzentrieren sich hauptsächlich auf Siedungsanwendungen
• Mangel an systematischen Modellen für kompressible nicht-ideale Strömungen

Relative Vorteile dieses Artikels

1. Systematischer Ansatz: Vollständiger Dual-Verteilungs-Rahmen
2. Theoretische Grundlagen: Strenge Multiskalen-Analyse
3. Vollständige Validierung: Umfassende Validierung von grundlegenden physikalischen Größen bis zu komplexen Strömungsphänomenen

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Erfolgreiche Entwicklung eines LBM-Modells für kompressible nicht-ideale Fluide
2. Realisierung numerischer Stabilität und physikalischer Konsistenz durch quasi-Gleichgewichts-Attraktor-Mechanismus
3. Erstmalige Realisierung von Tropfen-Stoßwellen-Wechselwirkungssimulationen bei hoher Machzahl im LBM-Rahmen
4. Modell erhält thermodynamische Konsistenz über einen breiten Bereich von Strömungszuständen

Einschränkungen

1. Einfaches Einzelrelaxationszeit-Modell: Gegenwärtig wird das einfache BGK-Kollisionsmodell verwendet, das für Turbulenz-Simulationen möglicherweise nicht robust genug ist
2. Beschränkung der Zustandsgleichung: Hauptsächlich Validierung der van-der-Waals-Gleichung; andere komplexe Zustandsgleichungen erfordern weitere Verifikation
3. Begrenzte dreidimensionale Validierung: Die meisten Validierungen konzentrieren sich auf eindimensionale und zweidimensionale Konfigurationen

Zukünftige Richtungen

1. Erweiterung auf robustere Kollisionsmodelle: Für Turbulenz-Simulationen
2. Komplexere Zustandsgleichungen: Wie Peng-Robinson und andere
3. Großskalige dreidimensionale Anwendungen: Anwendung auf praktische Ingenieurprobleme

Tiefgreifende Bewertung

Stärken

1. Theoretische Strenge: Vollständige Multiskalen-Analyse und mathematische Herleitung
2. Methodische Innovation: Der quasi-Gleichgewichts-Attraktor-Mechanismus ist eine wichtige Innovation
3. Umfassende Validierung: Systematische Validierung von grundlegenden Physiken bis zu komplexen Anwendungen
4. Praktischer Wert: Lösung wichtiger Ingenieuranwendungsprobleme
5. Numerische Effizienz: Verwendung von Standard-Gittern vermeidet Rechenkomplexität

Schwächen

1. Einfaches Kollisionsmodell: BGK-Modell hat begrenzte Anwendbarkeit auf komplexe Strömungen
2. Parameterabhängigkeit: Systematische Anleitung zur Auswahl bestimmter Korrekturparameter fehlt
3. Unzureichende dreidimensionale Validierung: Weitere Validierung mit dreidimensionalen komplexen Geometrien erforderlich

Einfluss

1. Akademischer Beitrag: Eröffnet neue Richtungen für LBM in kompressiblen nicht-idealen Strömungen
2. Ingenieurwert: Bereitstellung wichtiger numerischer Werkzeuge für überkritische Fluidtechnologie
3. Methodologische Bedeutung: Das Konzept des quasi-Gleichgewichts-Attraktors könnte andere numerische Methoden beeinflussen

Anwendungsszenarien

1. Überkritische Fluidtechnologie: CO₂-Turbinen, organische Rankine-Zyklen
2. Mehrphasenströmungen: Tropfenzerfall, Phasenwechselprozesse
3. Stoßwellendynamik: Komplexe Wellenwechselwirkungen in Hochgeschwindigkeitsströmungen
4. Direkte numerische Simulation: Grundlagenforschung mit hoher Präzision

Literaturverzeichnis

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.