Lattice Boltzmann model for non-ideal compressible fluid dynamics

Dieser Beitrag stellt ein neuartiges Gitter-Boltzmann-Modell vor, das auf Gittern mit ersten Nachbarn und Korrekturtermen im Quasi-Gleichgewicht basiert, um kompressible, nicht-ideale Strömungen mit thermodynamischer Konsistenz und numerischer Stabilität zu simulieren, und bestätigt dessen Genauigkeit erfolgreich durch quantitative Simulationen von Stoß-Tropfen-Wechselwirkungen bei Mach-Zahlen bis zu 1,47.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen zu simulieren, wie sich eine Flüssigkeit auf einem Computer verhält. Lange Zeit waren Computer hervorragend darin, „ideale" Flüssigkeiten zu simulieren – wie Wasser, das sanft in einem Fluss fließt, oder Luft, die langsam um einen Flügel strömt. Diese Flüssigkeiten folgen einfachen, vorhersehbaren Regeln.

Was passiert jedoch, wenn sich die Flüssigkeit unter extremem Druck und Hitze befindet und sich wie ein dichtes Gas verhält, das fast eine Flüssigkeit ist, oder wie eine Flüssigkeit, die fast ein Gas ist? Dies ist die Welt der nicht-idealen kompressiblen Flüssigkeiten. Stellen Sie sich dies als eine Flüssigkeit vor, die „gestresst" ist und sich seltsam verhält, die sich weigert, den einfachen Regeln der idealen Welt zu folgen. Dies geschieht in fortschrittlichen Technologien wie Turbinen mit überkritischem CO2 und organischen Energiezyklen.

Das Problem besteht darin, dass bestehende Computerwerkzeuge mit diesen gestressten Flüssigkeiten Schwierigkeiten haben. Sie stürzen entweder ab, liefern falsche Antworten oder erfordern so enorme Rechenleistung, dass sie unbrauchbar werden.

Dieser Artikel stellt eine neue, intelligentere Methode zur Simulation dieser kniffligen Flüssigkeiten vor, die auf einer Methode namens Gitter-Boltzmann-Methode (LBM) basiert. Hier ist, wie der neue Ansatz der Autoren funktioniert, erklärt durch einfache Analogien:

1. Das „Zwei-Spur"-System

Die meisten alten Simulationsmethoden versuchen, alles (Masse, Geschwindigkeit, Energie) mit einem einzigen, komplizierten Regelwerk zu verfolgen. Die Autoren erkannten, dass dies wie der Versuch ist, ein Auto zu fahren und gleichzeitig ein Dutzend Bälle jonglieren – es wird unübersichtlich und instabil.

Stattdessen bauten sie ein Zwei-Spur-System auf:

• Spur A (Die Menge): Ein Regelwerk verfolgt die Dichte und Geschwindigkeit der Flüssigkeit (wie viele Teilchen gibt es und wohin bewegen sie sich?).
• Spur B (Die Energie): Ein zweites, separates Regelwerk verfolgt die Gesamtenergie.
  Durch diese Trennung wird der Computer nicht verwirrt. Es ist wie ein dedizierter Verkehrsleiter für Autos und ein separater für Kraftstoff, wodurch sichergestellt wird, dass keines der Systeme das andere zum Absturz bringt.

2. Der „Quasi-Gleichgewichts"-Attraktor

In der Physik möchten Flüssigkeiten natürlich in einen ruhigen Zustand namens „Gleichgewicht" übergehen. Unter diesen extremen Bedingungen wird die Flüssigkeit jedoch ständig hin und her gezogen, sodass sie sich nie ganz beruhigt.

Die Autoren erfanden einen cleveren Trick namens „Quasi-Gleichgewichts-Attraktor".

• Die Analogie: Stellen Sie sich einen Hund vor, der einem Ball hinterherjagt. Der Ball repräsentiert den „perfekten ruhigen Zustand". Der Hund repräsentiert die Flüssigkeit. In einer normalen Situation läuft der Hund geradeaus zum Ball.
• Das Problem: Bei dieser extremen Flüssigkeit bewegt sich der Ball ständig weg oder verändert seine Form. Wenn der Hund den Ball blindlings jagt, könnte er über eine Klippe rennen (die Simulation wird instabil).
• Die Lösung: Die Autoren gaben dem Hund ein „GPS", das vorhersagt, wo der Ball in einem splitternden Moment sein wird, basierend darauf, wie sich der Wind (Druck) und das Gelände (Dichte) verändern. Dieses „verschobene" Ziel ermöglicht es dem Hund, reibungslos zu laufen, ohne von der Klippe zu fallen. Dies stellt sicher, dass die Simulation stabil bleibt, selbst wenn sich die Flüssigkeit sehr schnell bewegt oder ihre Dichte schnell ändert.

3. Behebung der „spuriösen" Wärme

Wenn Flüssigkeiten sich schnell bewegen, erzeugen sie Wärme. In Standard-Computermodellen fließt die Wärme manchmal in die falsche Richtung oder erzeugt falsche „Geister"-Wärme, die in der Realität nicht existiert.

• Die Analogie: Es ist wie ein Thermostat, der glaubt, der Raum sei eiskalt, weil er den Zug von einem Fenster misst und nicht die tatsächliche Raumtemperatur.
• Die Korrektur: Die Autoren fügten ihren Gleichungen einen spezifischen „Korrekturterm" hinzu. Dies wirkt wie ein Filter, der die falschen Zugluft entfernt und sicherstellt, dass die Wärme genau so fließt, wie es die Physik vorschreibt (Fouriersches Gesetz), selbst unter diesen extremen, nicht-idealen Bedingungen.

4. Die „Stoßrohr"- und „Flüssigkeitssäulen"-Tests

Um zu beweisen, dass ihre neue Methode funktioniert, führten sie nicht nur Mathematik durch; sie führten extreme Tests durch:

• Das Stoßrohr: Sie simulierten eine plötzliche Druckexplosion (eine Stoßwelle), die sich durch ein dichtes Gas bewegt. Bei normalen Gasen verhalten sich diese Wellen auf eine bestimmte Weise. Bei diesen „nicht-idealen" Gasen können die Wellen etwas Seltsames tun: einen „Expansionsstoß" (eine Welle, die sich ausbreitet, aber dennoch wie ein scharfer Stoß wirkt). Ihr Modell sagte dieses seltsame Verhalten erfolgreich voraus, was ältere Modelle verpasst hatten.
• Die Flüssigkeitssäule: Sie simulierten eine hochgeschwindigkeits Stoßwelle, die auf einen Flüssigkeitstropfen trifft. Dies ist ein sehr schwieriger Test, da der Stoß zurückprallt, reflektiert und den Tropfen zerreißt. Ihr Modell bewältigte die Kollision perfekt und stimmte mit realen Experimenten überein, bei denen der Flüssigkeitstropfen sich genau so abflacht und ausdehnt, wie es sollte.

Warum dies wichtig ist

Die Autoren behaupten, dass ihre Methode schnell, stabil und genau ist. Sie verwendet ein einfaches Gitter (wie ein normales Schachbrett) anstatt eines superkomplexen, gestreckten Gitters. Dies bedeutet, dass Wissenschaftler diese extremen, hochgeschwindigkeits, nicht-idealen Flüssigkeitsströmungen nun mit hoher Präzision auf Standardcomputern simulieren können.

Zusammenfassend: Der Artikel stellt ein neues „Fahrerhandbuch" für Computersimulationen vor, das es ihnen ermöglicht, Flüssigkeiten unter extremem Stress zu handhaben, ohne abzustürzen. Durch die Verwendung eines Zwei-Spur-Systems und eines intelligenten „GPS" für die Energie der Flüssigkeit können sie genau vorhersagen, wie sich diese komplexen Flüssigkeiten in Hochgeschwindigkeitsszenarien verhalten, und damit den Weg für bessere Designs fortschrittlicher Energiesysteme ebnen.

Technische Zusammenfassung

Technische Zusammenfassung: Gitter-Boltzmann-Modell für nicht-ideale kompressible Strömungsdynamik

Problemstellung
Die nicht-ideale kompressible Strömungsdynamik (NICFD) befasst sich mit Fluidverhalten in nahezu-, trans- und überkritischen Regimen, in denen Fluide signifikant von idealen Gasgesetzen abweichen. Diese Regime sind für moderne Energietechnologien entscheidend, wie etwa organische Rankine-Kreisläufe und Turbinen mit überkritischem CO₂. Experimentelle Daten in diesen Zuständen sind jedoch knapp und schwer zu beschaffen. Während direkte numerische Simulation (DNS) notwendig ist, um die komplexe Physik (z. B. Stoßstrukturen, Turbulenz in dichten Dämpfen) zu verstehen und technische Datenbanken zu generieren, haben bestehende numerische Werkzeuge oft Schwierigkeiten mit thermodynamischer Konsistenz, Stabilität und der genauen Wiederherstellung hydrodynamischer Grenzen für nicht-ideale Zustandsgleichungen (EOS). Insbesondere konventionelle Methoden der Strömungsmechanik (CFD) stoßen in der thermodynamisch instabilen Spinodalregion an Grenzen, in der das Quadrat der adiabatischen Schallgeschwindigkeit negativ wird, was spezielle Behandlungen wie die Maxwell'sche Gleichflächenregel erfordert. Die Gitter-Boltzmann-Methode (LBM), obwohl erfolgreich für inkompressible und Mehrphasenströmungen, hat historisch die nicht-ideale kompressible Thermodynamik in ihren kinetischen Formulierungen vernachlässigt und sich oft auf ad-hoc-Regularisierung oder erweiterte Stencils verlassen, die entweder die Effizienz oder die Galilei-Invarianz beeinträchtigen.

Methodik
Die Autoren schlagen ein neues kinetisches Modell und seine Lattice-Boltzmann-Realisierung vor, die entwickelt wurde, um die kompressiblen Navier-Stokes-Korteweg-Gleichungen (NSK) für nicht-ideale Fluide wiederherzustellen. Die Methodik basiert auf folgenden Säulen:

1. Ansatz mit doppelter Verteilungsfunktion (DDF): Das Modell verwendet zwei Verteilungsfunktionen: $f$ für Masse und Impuls sowie $g$ für die innere Energie ($\rho E$). Diese Entkopplung ermöglicht die Verwendung klassischer Gitter mit Nachbarn erster Ordnung (insbesondere D3Q27), ohne erweiterte Stencils zu benötigen.
2. Quasi-Gleichgewichts-Attraktoren: Der Kollisionsoperator nutzt eine doppelte Relaxations-BGK-ähnliche Struktur, die einen lokalen Gleichgewichts-Attraktor ($f^{eq}, g^{eq}$) und einen verschobenen Quasi-Gleichgewichts-Attraktor ($f^\star_\lambda, g^\star_\lambda$) umfasst.
   • Das lokale Gleichgewicht wird unter Verwendung eines Referenztemperaturparameters $\theta = P/\rho$ (thermodynamische Strömungsarbeit) anstelle der thermodynamischen Temperatur $T$ definiert, was die korrekte Breite der Maxwell-Verteilung unabhängig von der inneren Energie sicherstellt.
   • Das verschobene Gleichgewicht führt Korrekturen für die Korteweg-Kraft (Oberflächenspannung), die Volumenviskosität und den Wärmestrom ein. Es verwendet verschobene Werte für Strömungsgeschwindigkeit ($u^\star$), Referenztemperatur ($\theta^\star$) und thermodynamische Temperatur ($T^\star$).
3. Thermodynamische Konsistenz:
   • Die Referenztemperatur $\theta$ wird auf $P/\rho$ gesetzt, eine Formulierung, die sicherstellt, dass das Modell für nicht-ideale EOS (z. B. van-der-Waals) gültig bleibt.
   • Ein Korrekturterm im Wärmestrom ($q^c_\lambda$) wird eingeführt, um spurartige nicht-Fourier-Komponenten zu kompensieren, die aus Dichtegradienten in nicht-idealen Fluiden entstehen, wodurch die Wiederherstellung des Fourier-Gesetzes sichergestellt wird.
   • Die Verschiebung der Referenztemperatur ermöglicht es, die Volumenviskosität als unabhängigen, einstellbaren und positiv-definiten Parameter zu behandeln, wodurch das Problem der negativen Volumenviskosität umgangen wird, das in Standard-BGK-Modellen für nicht-ideale Drücke auftreten kann.
4. Diskretisierung: Das Modell verwendet ein explizites Schema zweiter Ordnung, das aus der Trapez-Integration entlang von Charakteristiken abgeleitet ist. Es verwendet ein Standard-D3Q27-Gitter und behandelt Quellterme (Korteweg-Kraft) über das verschobene Gleichgewicht, wodurch Diskretisierungsfehler bei expliziten Quelltermen vermieden werden.

Hauptbeiträge

• Neues kinetisches Rahmenwerk: Die Arbeit stellt ein kinetisches Modell vor, das die Zielgleichungen der Hydrodynamik (Massen-, Impuls- und Energiebilanz) für kompressible nicht-ideale Fluide unter Verwendung eines minimalen kinetischen Stencils (Nachbarn erster Ordnung) wiederherstellt.
• Positiv-definite Dissipation: Durch die Konstruktion des Modells über Quasi-Gleichgewichts-Attraktoren stellen die Autoren sicher, dass die Dissipationsraten der Navier-Stokes-Gleichungen (Scherviskosität und Volumenviskosität) positiv-definit und Galilei-invariant sind, ein häufiger Schwachpunkt früherer LBM-Versuche für kompressible nicht-ideale Strömungen.
• Thermodynamische Robustheit: Das Modell behandelt thermodynamische Instabilitätsregionen (Spinodale), in denen die Schallgeschwindigkeit im reibungsfreien Grenzfall imaginär wird, auf natürliche Weise, da die LBM die Ausbreitung entlang fester diskreter Charakteristiken erbt und nicht auf hyperbolische PDE-Löser angewiesen ist, die in diesen Regionen eine spezielle Behandlung erfordern.
• Validierung nicht-klassischer Phänomene: Die Studie liefert die erste quantitative Validierung der LBM für Stoß-Tropfen-Wechselwirkungen bei Mach-Zahlen bis zu 1,47 im Kontext nicht-idealer Fluide und demonstriert die Fähigkeit des Modells, komplexe Wellenstrukturen zu erfassen.

Ergebnisse
Das Modell wurde über eine Hierarchie zunehmend komplexer Konfigurationen validiert:

1. Lineare Stabilität und Dispersion: Simulationen der Schallgeschwindigkeit, der Dissipation von Scherwellen und der Wärmeleitfähigkeit zeigten eine hervorragende Übereinstimmung mit analytischen Lösungen für eine van-der-Waals-Flüssigkeit, angepasst an Stickstoffeigenschaften, über einen Bereich von Mach-Zahlen und reduzierten Temperaturen.
2. Konsistenz in Mehrphasensystemen: Das Modell reproduzierte genau die Dichten und Grenzflächenprofile im koexistierenden Zustand von Flüssigkeit und Dampf, zeigte eine Konvergenz zweiter Ordnung bei Gitterverfeinerung und stimmte mit der Maxwell'schen Gleichflächenregel überein.
3. Nicht-klassische Stoßrohre: Drei eindimensionale Stoßrohrfälle wurden simuliert, um nicht-klassische Wellendynamiken zu untersuchen, die durch die fundamentale Ableitung der Gasdynamik ($\Gamma$) charakterisiert sind.
   • Fälle mit $\Gamma < 0$ erfassten korrekt die Bildung von Expansionsstößen (Verstärkung von Expansionswellen) und Kompressionsfächern (Ausbreitung von Kompressionswellen), Phänomene, die sich von der klassischen Gasdynamik unterscheiden.
   • Das Modell wechselte erfolgreich zwischen klassischen ($\Gamma > 0$) und nicht-klassischen ($\Gamma < 0$) Regimen innerhalb einer einzigen Simulation.
4. Stoß-Flüssigkeitssäulen-Wechselwirkung: Eine zweidimensionale Simulation einer planaren Stoßwelle, die mit einer gesättigten Flüssigkeitssäule nahe dem kritischen Punkt wechselwirkt ($M_a = 1,47$), wurde durchgeführt.
   • Das Modell erfasste komplexe Wellenstrukturen, einschließlich transmittierter Stöße, reflektierter Stöße, Mach-Stämme und Tripelpunkte.
   • Der quantitative Vergleich der Entwicklung der Breite der Flüssigkeitssäule mit experimentellen Daten und früheren numerischen Ergebnissen zeigte eine gute Übereinstimmung.
   • Eine thermodynamische Analyse ergab, dass die Stoßwechselwirkung das umgebende Dampf in einen überkritischen Zustand treibt, während die Flüssigkeitssäule marginal unterkritisch bleibt, wodurch die starke thermodynamische Kopplung genau erfasst wird.

Bedeutung und Behauptungen
Die Autoren behaupten, dass diese Arbeit die Anwendbarkeit von Gitter-Boltzmann-Methoden auf hochgeschwindigkeitskompressible Strömungen mit nicht-idealem Verhalten unter Verwendung eines minimalen kinetischen Stencils erweitert. Die primäre Bedeutung liegt in der Bereitstellung eines thermodynamisch konsistenten und numerisch stabilen Rahmens, der nicht auf erweiterte Stencils oder ad-hoc-Regularisierung angewiesen ist. Das Modell wird als reduzierte kinetische Theorie vorgestellt, die auf den NSK-hydrodynamischen Grenzfall abzielt und in der Lage ist, das gesamte Spektrum nicht-idealen Verhaltens, einschließlich der Spinodalregion, zu behandeln, ohne die numerischen Instabilitäten, die bei konventioneller CFD in diesen Regimen häufig auftreten. Die erfolgreiche Simulation von Stoß-Tropfen-Wechselwirkungen bei hohen Mach-Zahlen dient als strenge Validierung und demonstriert das Potenzial der Methode zur Untersuchung komplexer Strömungen in überkritischen Technologien und prozessen, die durch Phasenübergänge angetrieben werden.