An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model

In dieser Arbeit wird ein LES-Modell vorgestellt, das durch ein thermodynamisch konsistentes Vier-Gleichungs-Modell, ein neues phasengeschränktes $\Sigma$-Sprühmodell zur Vorhersage der Subgrid-Oberflächenfläche sowie eine verbesserte Modellierung der endlichen Phasenwechselraten eine präzise und effiziente Simulation von Multiphase-Systemen ermöglicht.

Das Wesentliche

Das Problem: Der „Chaos-Sprühnebel“ in der Maschine

Stellen Sie sich vor, Sie versuchen, das Verhalten eines riesigen Feuerwerks oder eines hochmodernen Raketentriebwerks am Computer zu simulieren. Das Problem ist: In diesen Systemen passiert alles gleichzeitig und extrem schnell. Flüssiger Treibstoff wird mit enormem Druck eingespritzt, zerstäubt in Milliarden winziger Tröpfchen, die dann verdampfen, explodieren oder bei plötzlichem Druckabfall förmlich „aufplatzen“ (das nennt man Blasenbildung oder Flashing).

Bisher standen Wissenschaftler vor einem Dilemma:

1. Die „Perfektionisten“: Man versucht, jedes einzelne Tröpfchen zu berechnen. Das ist so, als wollte man jedes einzelne Sandkorn an einem Strand zählen. Es ist unglaublich genau, aber es dauert Jahre, bis der Computer fertig ist.
2. Die „Schätzer“: Man nutzt vereinfachte Modelle, die nur grob sagen: „Hier ist es nass, dort ist es gasförmig.“ Das geht schnell, ist aber oft ungenau – wie eine Wettervorhersage, die sagt „es regnet“, aber nicht, ob es ein sanfter Nieselregen oder ein heftiges Gewitter ist.

Die Lösung: Das „Schlaue Team“ (Das neue LES-Modell)

Die Forscher aus Stanford und Schweden haben nun ein Modell entwickelt, das wie ein extrem intelligentes Team arbeitet. Sie nutzen einen Trick, den man „LES“ (Large Eddy Simulation) nennt.

Stellen Sie sich das wie eine Google-Maps-Navigation vor:

• Wenn Sie auf der Autobahn sind (große Strömungen), zeigt Ihnen das Modell die großen Wege und Richtungen an. Das ist einfach und schnell zu berechnen.
• Sobald Sie aber in eine enge, verwinkelte Altstadt einfahren (die winzigen Tröpfchen und die Grenzflächen zwischen Flüssigkeit und Gas), schaltet das Modell auf einen „intelligenten Schätzer“ um. Es berechnet nicht jedes Tröpfchen einzeln, aber es hat eine sehr gute mathematische Formel, die genau vorhersagt, wie sich das Chaos in den kleinen Gassen verhalten muss.

Die drei Superkräfte des Modells

Damit dieses „Team“ funktioniert, haben die Forscher drei neue Werkzeuge eingebaut:

1. Der „Thermodynamik-Türsteher“ (Finite-Rate Phase Change)

Wenn Flüssigkeit zu Gas wird, passiert das nicht magisch und sofort. Es braucht Zeit. Frühere Modelle waren oft zu gierig: Sie haben die Flüssigkeit sofort „weggezaubert“, als wäre sie verdampft. Das neue Modell hat einen „Türsteher“: Es berechnet die Geschwindigkeit des Übergangs so, dass sie physikalisch Sinn ergibt. Es achtet darauf, dass die Energie und die Chemie im Gleichgewicht bleiben – so als würde man beim Kochen darauf achten, dass die Hitze genau so schnell in den Topf zieht, wie das Wasser verdampfen kann, ohne dass die Küche explodiert.

2. Das „unsichtbare Netz“ (Subgrid Spray)

In der Simulation sind die Tröpfchen oft viel kleiner, als die Rechenzellen des Computers groß sind. Das ist, als wollten Sie ein winziges Mückenschwarm-Verhalten mit einem riesigen Pinsel malen. Die Forscher haben eine Formel (das $\Sigma$-Modell) entwickelt, die eine „unsichtbare Oberfläche“ berechnet. Sie sagt dem Computer: „Du kannst die Tröpfchen zwar nicht einzeln sehen, aber ich verspreche dir, dass hier eine riesige Oberfläche an Flüssigkeit existiert, die bereit ist, zu verdampfen.“

3. Der „Regime-Detektor“ (Dense vs. Dilute)

Das Modell ist schlau genug zu erkennen, in welcher Situation es sich befindet.

• Ist es ein „dichter Wald“ (viele Tröpfchen dicht beieinander)? Dann berechnet es, wie sie gegeneinander prallen und zerbrechen.
• Ist es eine „leere Wüste“ (einzelne Tröpfchen im Gas)? Dann berechnet es, wie sie einsam verdampfen.

Warum ist das wichtig?

Die Forscher haben ihr Modell mit echten Experimenten (dem sogenannten „ECN Spray A“-Test) verglichen. Das Ergebnis: Die Computer-Simulation und die echten Messungen aus dem Labor sahen fast identisch aus.

Was bringt uns das im Alltag?
Wenn wir Motoren effizienter machen wollen (weniger Spritverbrauch, weniger Abgase) oder Raketen sicherer bauen wollen, müssen wir verstehen, wie der Treibstoff zerstäubt. Mit diesem Modell können Ingenieure am Computer tausende Designs testen, ohne jedes Mal eine echte, teure und gefährliche Explosion im Labor auslösen zu müssen. Es ist wie ein hochpräziser Flugsimulator für die Welt der kleinsten Tröpfchen.

Technische Zusammenfassung

Zusammenfassung: LES-Modell mit endlicher Phasenwechselrate und Subgrid-Spray auf Basis eines thermodynamisch konsistenten Vier-Gleichungs-Multiphasenmodells

1. Problemstellung

Die Simulation von großskaligen, multiphasischen und multikomponentigen Systemen (wie z. B. in Einspritzsystemen von Verbrennungsmotoren oder Raketentriebwerken) ist aufgrund der extremen Dichteverhältnisse, der komplexen Grenzflächenphysik und der notwendigen Massentransfermodelle numerisch äußerst anspruchsvoll.

Bisherige Ansätze stehen vor einem Dilemma:

• Nicht-Gleichgewichtsmodelle (z. B. Sieben-Gleichungs-Modelle) sind physikalisch präzise, aber aufgrund der Kopplung von Druck-, Temperatur- und Geschwindigkeitsfeldern zwischen den Phasen extrem rechenintensiv.
• Gleichgewichtsmodelle (Homogeneous Equilibrium Model, HEM) sparen zwar Rechenzeit, indem sie mechanisches und thermisches Gleichgewicht voraussetzen, versagen jedoch bei der Vorhersage der Zeitskalen des Massentransfers in unterauflösten Bereichen (Subgrid-Skalen), was zu unphysikalischen Massentransferraten führt.
• Zudem fehlt es bei Eulerischen Subgrid-Spray-Modellen oft an einer konsistenten Kopplung, die verhindert, dass chemische Spezies unphysikalisch über die Phasengrenzen "lecken".

2. Methodik

Die Autoren erweiterten ein bestehendes hochauflösendes numerisches Verfahren (basierend auf einem Vier-Gleichungs-Modell) um zwei wesentliche Komponenten:

• Thermodynamisch begrenztes Finite-Rate-Phasenwechselmodell: Anstatt einen unendlichen Massentransfer (HEM) anzunehmen, wird die Phasenwechselrate über eine Zeitkonstante $\tau$ gesteuert. Um sicherzustellen, dass das Modell physikalisch bleibt, wird $\tau$ so definiert, dass der Massentransfer niemals über das Niveau hinausgeht, das für das Erreichen des chemischen Gleichgewichts (basierend auf der Gibbs-Energie) notwendig wäre. Als kinetisches Modell wird der Hertz-Knudsen-Ansatz verwendet, der sowohl Verdampfung als auch Flash-Verdampfung (Flashing) abdeckt.
• Phasengeschränktes $\Sigma$-Spray-Modell: Um die für den Phasenwechsel kritische Grenzflächenoberfläche ($\Sigma$) auf Subgrid-Ebene zu modellieren, wird eine modifizierte Transportgleichung für die Oberflächenarealdichte eingeführt. Die Neuerung besteht darin, die Diffusion der Oberflächenrate so zu formulieren, dass sie auf die Gasphase beschränkt bleibt ("phase-confined"). Dies verhindert das unphysikalische Lecken von Subgrid-Flüssigkeit in die Gasphase. Das Modell unterscheidet zudem zwischen dichten Sprays (turbulente Fragmentierung) und diluten Sprays (Sekundärzerfall und Koaleszenz).

3. Zentrale Beiträge (Key Contributions)

1. Effizienz durch Gleichgewichtsannahmen: Die Nutzung des Vier-Gleichungs-Modells (mechanisches und thermisches Gleichgewicht zwischen den Phasen) ermöglicht massive Rechenersparnisse gegenüber Voll-Nicht-Gleichgewichtsmodellen, ohne die Genauigkeit der Spray-Dynamik zu opfern.
2. Thermodynamische Konsistenz: Einführung eines Mechanismus, der die endliche Phasenwechselrate an die chemische Potential-Äquilibrierung koppelt. Dies verhindert das Überschießen (Overshooting) von Massentransferraten.
3. Robustes Subgrid-Modell: Entwicklung eines $\Sigma$-Modells, das die Grenzflächenfläche in LES-Simulationen konsistent transportiert und dabei die physikalische Trennung der Phasen wahrt.
4. Regime-Indikatoren: Implementierung von Indikatoren (basierend auf der Jakob-Zahl und der Hinze-Skala), um automatisch zwischen Flashing, Verdampfung sowie dichten und diluten Spray-Regimen zu unterscheiden.

4. Ergebnisse und Validierung

Das Modell wurde anhand des ECN (Engine Combustion Network) Spray A Falls validiert:

• Nicht-verdampfender Fall: Die Simulationen der projizierten Massendichte (PMD) und der projizierten Oberflächenfläche (PSA) zeigten eine exzellente Übereinstimmung mit experimentellen Daten. Die Vorhersage des mittleren Tropfendurchmessers (SMD) lag im experimentellen Bereich, was die Korrektheit des Subgrid-Spray-Modells bestätigt.
• Verdampfender Fall: Die Vorhersage der Flüssigkeits-Penetrationstiefe und der Gaswolken-Ausdehnung stimmte sehr gut mit den experimentellen Messungen überein. Dies beweist, dass die Kopplung zwischen dem $\Sigma$-Modell (Oberfläche) und dem endlichen Phasenwechselmodell (Massentransfer) die komplexe Spray-Physik korrekt erfasst.

5. Bedeutung (Significance)

Die Arbeit liefert einen leistungsfähigen und physikalisch konsistenten Rahmen für die Large Eddy Simulation (LES) von komplexen Multiphasenströmungen. Die Kombination aus Recheneffizienz (durch das Vier-Gleichungs-Modell) und hoher physikalischer Genauigkeit (durch die thermodynamisch begrenzten Raten) macht das Modell besonders wertvoll für die industrielle Anwendung in der Verbrennungsforschung, der Raketentechnik und der Brandbekämpfung, wo sowohl schnelle Phasenwechsel als auch hochgradig turbulente Sprays auftreten.