Turbulence Modelling of Mixing Layers under Anisotropic Strain

Diese Studie untersucht den Einfluss anisotroper Dehnungsraten auf turbulente Mischschichten unter Verwendung des K-L-RANS-Modells, zeigt auf, dass ein transversaler Dehnungsabschluss die Vorhersagegenauigkeit im Vergleich zum Standardansatz für isotrope Bedingungen verbessert, und schlägt entsprechende Modifikationen für die K-$\epsilon$- und K-$\omega$-Modelle vor.

Das Wesentliche

Stellen Sie sich eine Schüssel mit zwei Flüssigkeiten unterschiedlicher Farbe vor, wie Öl und Wasser, die auf einem Tisch steht. Wenn Sie den Tisch schütteln, wird die Grenze zwischen ihnen unordentlich und beginnt sich zu vermischen. Dies ist ähnlich wie bei einer turbulenten Mischschicht in der Physik: Zwei Fluide unterschiedlicher Dichte werden zusammengedrückt und erzeugen ein chaotisches, wirbelndes Durcheinander.

Dieser Artikel handelt davon, zu verstehen, was passiert, wenn man nicht nur den Tisch schüttelt, sondern auch den gesamten Raum, in dem die Vermischung stattfindet, dehnt oder staucht.

Hier ist eine Aufschlüsselung der Geschichte des Artikels, unter Verwendung einfacher Analogien:

1. Der Schauplatz: Den Raum dehnen

In vielen realen Szenarien – wie bei einer explodierenden Stern (Supernova) oder einem komprimierten Kernfusionsbomben – ist der Raum, in dem sich die Fluide vermischen, nicht einfach nur stillstehend. Der Raum selbst dehnt sich aus oder zieht sich zusammen.

• Die Analogie: Stellen Sie sich die Mischschicht als einen Teig vor, der geknetet wird. Normalerweise untersuchen Wissenschaftler, wie sich der Teig vermengt, wenn man ihn einfach nur hin und her schiebt. Aber in diesem Artikel fragen die Autoren: „Was passiert, wenn jemand, während Sie kneten, auch den Tisch, auf dem der Teig liegt, zieht, ihn längs dehnt oder quer staucht?"
• Das Problem: Die „Dehnung" (Scherung) ist nicht in alle Richtungen gleich. Wenn Sie ein Gummiband ziehen, wird es in einer Richtung länger, aber in den anderen dünner. Dies wird als anisotrope Scherung bezeichnet. Die meisten Computermodelle, die diese Vermischungen vorhersagen, gehen davon aus, dass die Dehnung in alle Richtungen gleich ist (wie beim Aufblasen eines perfekten Ballons), was der Realität nicht entspricht.

2. Das Werkzeug: Das „K-L"-Modell

Um vorherzusagen, wie sich die Fluide vermischen, verwenden die Autoren ein Computerprogramm namens K-L-Turbulenzmodell.

• Die Analogie: Denken Sie an dieses Modell als Kochbuch zur Vorhersage von Chaos. Es verfolgt zwei Hauptzutaten:
  1. Wie viel Energie in den Wirbeln steckt (Turbulente kinetische Energie).
  2. Wie groß die Wirbel sind (Turbulente Längenskala).
• Das Modell versucht zu erraten, wie groß die Wirbel werden, wenn sich die Fluide vermischen. Der knifflige Teil ist eine Regel im Rezept namens „Volumenkompression"-Term. Diese Regel sagt dem Modell, wie sich die Größe der Wirbel ändert, wenn der gesamte Raum gestaucht oder gedehnt wird.

3. Das Experiment: Testen von drei verschiedenen Regeln

Die Autoren führten Computersimulationen durch, um zu sehen, welche „Regel" für die Volumenkompression am besten funktionierte, wenn der Raum in bestimmten Richtungen gedehnt wurde. Sie testeten drei Versionen des Rezepts:

1. Die „Durchschnitts"-Regel: Sie geht davon aus, dass die Dehnung in alle Richtungen gleich ist (die Standardeinstellung).
2. Die „Längs"-Regel: Sie geht davon aus, dass sich die Größe der Wirbel nur danach ändert, wie stark der Raum entlang der Mischrichtung gedehnt wird.
3. Die „Seitwärts"-Regel: Sie geht davon aus, dass sich die Größe der Wirbel danach ändert, wie stark der Raum quer zur Mischrichtung (senkrecht zur Strömung) gedehnt wird.

4. Die Ergebnisse: Die „Seitwärts"-Regel gewinnt

Die Autoren verglichen ihre Computer-Vorhersagen mit hochdetaillierten, hochauflösenden Simulationen (die als „perfekter" Referenzwert dienen).

• Die Erkenntnis: Die Standard-„Durchschnitts"-Regel war okay, aber nicht großartig. Die „Längs"-Regel machte die Vorhersagen tatsächlich schlechter.
• Der Gewinner: Die „Seitwärts"-Regel (unter Verwendung der transversalen Scherung) war die genaueste.
• Warum? Die Autoren erklären, dass sich beim Dehnen einer Mischschicht die großen „Wirbel" (Eddies) je nach Richtung unterschiedlich verhalten. Es stellt sich heraus, dass die Größe dieser Wirbel empfindlicher darauf reagiert, wie sich der Raum seitwärts (transversal) verändert, als darauf, wie er sich längs verändert. Indem sie die seitliche Dehnung nutzten, um die Größe der Wirbel im Rezept anzupassen, sagte das Modell die Mischbreite und die Energie viel genauer voraus.

5. Das größere Bild: Ein neues „Drei-Teile"-Rezept

Der Artikel untersuchte auch, wie man diese komplexen Gleichungen in ein „Auftrieb-Widerstand"-Modell vereinfachen kann (eine einfachere Art, über die Vermischung nachzudenken).

• Sie erkannten, dass sich die „Breite der Mischung" und die „Größe der Wirbel" tatsächlich auf unterschiedliche Kräfte reagieren. Die Breite dehnt sich mit dem Längs-Zug, aber die Wirbelgröße reagiert auf den seitlichen Druck.
• Die Schlussfolgerung: Um die beste Vorhersage zu erhalten, benötigen Sie ein Modell, das diese beiden Dinge separat behandelt. Anstatt einer Regel für alles benötigen Sie ein Drei-Teile-Modell, das die Breite und die Wirbelgröße unabhängig voneinander entwickelt.

Zusammenfassung

Kurz gesagt geht es in diesem Artikel darum, ein Computermodell zu verbessern, das verwendet wird, um vorherzusagen, wie sich Fluide vermischen, wenn der Raum um sie herum verzerrt wird. Die Autoren entdeckten, dass die Standardmethode zur Berechnung, wie sich die „Wirbel" verkleinern oder vergrößern, unter diesen spezifischen Bedingungen falsch war. Indem sie die Regel änderten, um zu betrachten, wie sich der Raum seitwärts dehnt, anstatt sie einfach zu mitteln, machten sie das Modell viel genauer. Dies hilft Wissenschaftlern, komplexe Ereignisse wie Sternexplosionen oder Fusionsenergie-Experimente besser zu verstehen, bei denen Fluide ständig auf ungleichmäßige Weise gestaucht und gedehnt werden.

Technische Zusammenfassung

1. Problemstellung

Turbulente Mischschichten, wie sie durch Rayleigh-Taylor- (RTI) und Richtmyer-Meshkov-Instabilitäten (RMI) erzeugt werden, sind grundlegend für die Physik hoher Energiedichten, einschließlich der Trägheitseinschlussfusion (ICF) und astrophysikalischer Phänomene wie Supernovae. In diesen Szenarien erfährt die Strömung häufig anisotrope Dehnungsraten aufgrund radialer Bewegung in konvergenten Geometrien (Implosionen) oder divergenten Geometrien (Explosionen).

Aktuelle Reynolds-gemittelte Navier-Stokes (RANS)-Turbulenzmodelle, insbesondere das weit verbreitete K-L-Modell (ein Zwei-Gleichungs-Modell, das für turbulente kinetische Energie $K$ und turbulente Längenskala $L$ löst), stützen sich typischerweise auf Schließungsannahmen, die für isotrope Dehnung oder inkompressible Strömungen abgeleitet wurden. Die Standard-Schließung für den Term der Volumenkompresion in der Transportgleichung der Längenskala geht davon aus, dass sich die Längenskala basierend auf der mittleren isotropen Dehnungsrate (Divergenz der Geschwindigkeit geteilt durch 3) entwickelt. In realistischen anisotropen Szenarien (z. B. sphärische Implosionen) unterscheiden sich die Dehnungsraten in radialer (axialer) und circumferentieller (transversaler) Richtung jedoch erheblich. Das Papier adressiert das Fehlen einer theoretischen Definition und einer genauen Modellierung dafür, wie sich Volumenkompresion unter diesen anisotropen Bedingungen auf die turbulente Längenskala auswirkt.

2. Methodik

Die Autoren führten eine systematische Untersuchung mittels eindimensionaler RANS-Simulationen durch, um die Effekte der Dehnungsrichtung auf die Entwicklung von Mischschichten zu isolieren.

• Numerischer Rahmen: Die Simulationen nutzten das K-L-Turbulenzmodell, implementiert im FLAMENCO-Code. Das Modell löst Transportgleichungen für Favre-gemittelte turbulente kinetische Energie ($\tilde{K}$) und turbulente Längenskala ($L$).
• Testfälle: Die Studie basierte auf dem „Quarter-Scale $\theta$-Group"-Richtmyer-Meshkov-Instabilitätsfall (schwer-zu-leicht-Konfiguration, Atwood-Zahl $A_t=0.5$, Schock-Mach-Zahl $M_a \approx 1.84$).
• Anwendung der Dehnung: Um die Effekte der Anisotropie zu isolieren, wurden gleichmäßige normale Dehnungsraten auf das Simulationsgebiet in zwei unterschiedlichen Richtungen angewendet:
  • Axiale Dehnung ($S_A$): Dehnung in Richtung der Mischung (normal zur Grenzfläche).
  • Transversale Dehnung ($S_T$): Dehnung in den Richtungen parallel zur Grenzfläche (homogene Ebene).
  • Sowohl Profile mit konstanter Geschwindigkeit als auch mit konstanter Dehnungsrate wurden getestet.
• Variationen der Schließung: Drei verschiedene Schließungen für den Term der Volumenkompresion ($\bar{\rho} L S_\phi$) in der Längenskalengleichung wurden verglichen:
  1. Isotrope Schließung (Standard): Skaliert $L$ mit der mittleren Dehnungsrate ($S_I = \nabla \cdot \mathbf{u} / 3$).
  2. Axiale Schließung: Skaliert $L$ mit der axialen Dehnungsrate ($S_A$).
  3. Transversale Schließung: Skaliert $L$ mit der transversalen Dehnungsrate ($S_T$).
• Validierung: Die RANS-Ergebnisse wurden gegen Daten aus hochauflösenden Implicit Large Eddy Simulationen (ILES) früherer Studien (Pascoe et al.) verglichen, die als Ground Truth dienten.

3. Hauptbeiträge

• Identifizierung von Schließungsdefiziten: Das Papier zeigt, dass die Standardannahme der isotropen Schließung die Entwicklung turbulenter Mischschichten unter anisotroper Dehnung, insbesondere hinsichtlich der turbulenten kinetischen Energie (TKE) und der integralen Breite, nicht genau vorhersagen kann.
• Vorschlag einer transversalen Schließung: Die Autoren schlagen eine neue Schließungsstrategie vor und validieren diese, bei der die Entwicklung der turbulenten Längenskala unter Volumenkompresion durch die transversale Dehnungsrate ($S_T$) und nicht durch die mittlere isotrope Dehnung angetrieben wird.
• Ableitung eines Drei-Gleichungs-Auftriebs-Widerstands-Modells: Durch eine Selbstähnlichkeitsanalyse leiteten die Autoren ein vereinfachtes Auftriebs-Widerstands-Modell ab. Sie stellten fest, dass eine einzelne Längenskala nicht gleichzeitig die Skalierung der integralen Breite (die der axialen Dehnung folgt) und der turbulenten Längenskala (die der transversalen Dehnung folgt) erfassen kann. Folglich entwickelten sie ein Drei-Gleichungs-Modell, das die integrale Breite und die turbulente Längenskala separat entwickelt.
• Generalisierung auf andere Modelle: Die Studie erweitert die Erkenntnisse auf andere populäre Zwei-Gleichungs-Modelle (K-$\epsilon$ und K-$\omega$) und leitet äquivalente Modifikationen ihrer Volumenkompresions-Terme basierend auf der transversalen Dehnungsschließung ab.

4. Ergebnisse

• Integrale Breite ($W$):
  • Unter axialer Dehnung wird die Breite der Mischschicht direkt durch den Geschwindigkeitsgradienten der mittleren Strömung gedehnt oder komprimiert. Alle Schließungen erfassten diesen Trend, aber die transversale Schließung lieferte die genaueste Vorhersage der turbulenten Wachstumsrate und reduzierte den mittleren absoluten prozentualen Fehler (MAPE) im Vergleich zur isotropen Schließung um den Faktor vier.
  • Unter transversaler Dehnung zeigte die isotrope Schließung fast keine Variation der Breite (fälschlicherweise wurden Scherproduktion und Kompressionseffekte ausgeglichen), während die transversale Schließung den in ILES-Daten beobachteten leichten Anstieg der Wachstumsrate korrekt vorhersagte.
• Turbulente kinetische Energie (TKE):
  • Die transversale Schließung schnitt bei der Vorhersage der TKE deutlich besser ab als die anderen Methoden.
  • Die axiale Schließung schnitt am schlechtesten ab und sagte oft das falsche Vorzeichen der Änderung der TKE unter Expansion voraus.
  • Die isotrope Schließung überschätzte den Einfluss der Dehnung auf die TKE.
• Selbstähnlichkeitsanalyse:
  • Die Analyse ergab, dass unter anisotroper Dehnung das Verhältnis zwischen dem Peak der turbulenten Längenskala und der integralen Breite ($\beta_W$) nicht konstant ist, was die in vielen K-L-Kalibrierungen verwendete Standardannahme der Selbstähnlichkeit verletzt.
  • Dies erforderte die Entwicklung des Drei-Gleichungs-Auftriebs-Widerstands-Modells, bei dem die integrale Breite mit der axialen Dehnung ($S_A$) skaliert und die turbulente Längenskala mit der transversalen Dehnung ($S_T$) skaliert.
• Modelläquivalenz: Die abgeleiteten Modifikationen für das K-L-Modell wurden erfolgreich auf die K-$\epsilon$- und K-$\omega$-Modelle abgebildet, was darauf hindeutet, dass auch der Term der Volumenkompresion in diesen Modellen angepasst werden sollte, um die transversale Dehnung und nicht nur die mittlere Divergenz zu berücksichtigen.

5. Bedeutung

Diese Arbeit stellt einen kritischen Fortschritt in der Modellierung kompressibler turbulenter Mischung in Anwendungen hoher Energiedichte dar.

• Verbesserte Vorhersagefähigkeit: Durch die Korrektur der Schließung der Volumenkompresion unter Berücksichtigung der Anisotropie kann das K-L-Modell (und seine Äquivalente) das Wachstum von Mischschichten und die TKE in Implosions-/Explosionsszenarien, die für die Trägheitseinschlussfusion (ICF) und astrophysikalische Simulationen (z. B. Supernovae) zentral sind, nun genauer vorhersagen.
• Physikalische Einsicht: Die Studie klärt die unterschiedlichen physikalischen Rollen der axialen und transversalen Dehnung: Axiale Dehnung treibt primär die geometrische Dehnung der Mischschicht an, während transversale Dehnung die Entwicklung der turbulenten Wirbel (Längenskala) und die Dissipation bestimmt.
• Praktische Anwendung: Die vorgeschlagenen Modifikationen sind relativ einfach in bestehenden RANS-Codes zu implementieren, bieten jedoch erhebliche Verbesserungen der Genauigkeit ohne die Rechenkosten von LES oder DNS. Das abgeleitete Drei-Gleichungs-Auftriebs-Widerstands-Modell bietet ein rechen-effizientes Werkzeug für ingenieurtechnische Anwendungen, bei denen vollständige RANS-Simulationen zu teuer sind.

Zusammenfassend argumentiert das Papier, dass für anisotrope Dehnungsumgebungen die turbulente Längenskala so modelliert werden sollte, dass sie sich mit der transversalen Dehnungsrate entwickelt, eine Erkenntnis, die die langjährige Annahme der isotropen Skalierung in Standard-Turbulenzschließungen in Frage stellt.