A theoretical one-dimensional model for variable-density Rayleigh-Taylor turbulence

Dieser Beitrag nimmt das theoretische Modell von Belen'kii & Fradkin aus dem Jahr 1965 für Rayleigh-Taylor-Turbulenz mit variabler Dichte erneut auf und erweitert es, indem er zeigt, dass die vollständige Ähnlichkeitsgleichung die Merkmale nicht-Boussinesq-Strömungen präzise erfasst und dass eine massenkorrigierte vereinfachte Lösung die komplexen Mischungsprozesse, die durch das Wettstreiten zwischen Diffusion und Massenerhaltung bestimmt werden, wirksam approximiert.

Das Wesentliche

Stellen Sie sich zwei Fluide vor, eines schwer (wie Honig) und eines leicht (wie Luft), die übereinander liegen. Die Schwerkraft möchte, dass das schwere Fluid absinkt und das leichte aufsteigt, doch sie stecken in einem chaotischen, wirbelnden Kampf an der Grenzfläche fest. Dies ist die Rayleigh-Taylor-Instabilität. Beim Mischen bilden sie eine turbulente „Suppe", in der schwere Spitzen nach unten tauchen und leichte Blasen aufsteigen.

Seit Jahrzehnten versuchen Wissenschaftler vorherzusagen, wie schnell diese Mischschicht wächst. Die meisten modernen Theorien gehen davon aus, dass die Fluide „fast" die gleiche Dichte haben, und verwenden eine einfache Faustregel. Dieser Artikel jedoch nimmt eine vergessene, 60 Jahre alte Theorie von 1965 von Belen'kii und Fradkin wieder auf, die einen anderen, genaueren Weg bietet, um dieses Chaos zu betrachten, insbesondere wenn der Dichteunterschied enorm ist.

Hier ist die Aufschlüsselung dessen, was der Artikel leistet, unter Verwendung einfacher Analogien:

1. Das vergessene Rezept

Die Autoren fanden ein altes „Rezept" (ein mathematisches Modell) dafür, wie diese Fluide sich mischen. Das ursprüngliche Rezept war auf Russisch verfasst, etwas unübersichtlich zu lesen und enthielt einige Tippfehler.

• Was sie taten: Sie reinigten das Rezept, übersetzten es und schrieben es in moderner, klarer Sprache neu.
• Die Kernidee: Anstatt das Mischen als komplexe 3D-Explosion zu betrachten, behandelten sie es wie ein eindimensionales Diffusionsproblem. Stellen Sie sich die Mischschicht nicht als chaotischen Sturm vor, sondern als einen einzelnen, sich ausbreitenden Fleck auf einem Blatt Papier. Sie modellierten das Ausbreiten dieses „Flecks" unter Verwendung des Konzepts der turbulenten Diffusivität (wie schnell sich das Chaos ausbreitet).

2. Die „Logarithmus"-Regel vs. die „Lineare" Regel

Die große Entdeckung in diesem Artikel betrifft das Wachstum der Mischschicht im Laufe der Zeit.

• Die alte Sichtweise: Die meisten Wissenschaftler glaubten, dass die Wachstumsrate von einer linearen Zahl abhängt, der sogenannten Atwood-Zahl (die den Unterschied zwischen schwerem und leichtem Fluid misst). Wenn sich der Unterschied verdoppelt, verdoppelt sich auch die Mischgeschwindigkeit.
• Die neue (alte) Sichtweise: Das Modell von 1965 legt nahe, dass das Wachstum vom natürlichen Logarithmus des Dichteverhältnisses ($\ln R$) abhängt.
  • Die Analogie: Denken Sie an die Atwood-Zahl wie an eine gerade Linie in einem Diagramm. Der Logarithmus ist wie eine Kurve, die sich abflacht. Der Artikel argumentiert, dass, wenn der Dichteunterschied enorm wird (wie beim Vergleich von Blei mit Luft), die Mischung nicht linear schneller wird; sie verlangsamt ihre Wachstumsrate und folgt dieser logarithmischen Kurve. Dies passt besser zu aktuellen Computersimulationen als die alte lineare Regel.

3. Die Asymmetrie zwischen „Schwer" und „Leicht"

Wenn sich schwere und leichte Fluide mischen, verhalten sie sich nicht gleich.

• Die Beobachtung: Das schwere Fluid bildet „Spitzen", die schnell nach unten tauchen, während das leichte Fluid „Blasen" bildet, die langsamer aufsteigen.
• Die Erkenntnis des Artikels: Das alte Modell von 1965 sagt diese Asymmetrie auf natürliche Weise voraus, ohne dass zusätzliche Anpassungen nötig sind. Es zeigt, dass, wenn der Dichteunterschied wächst, die „Spitzen" viel länger werden als die „Blasen".
• Die Geschwindigkeitsverschiebung: Der Artikel zeigt zudem, dass sich die Geschwindigkeit der Mischung zur Seite des leichten Fluids hin verschiebt.
  • Die Analogie: Stellen Sie sich ein Tauziehen vor, bei dem ein Team viel schwerer ist. Das Seil bewegt sich nicht einfach zur Mitte; der gesamte Mittelpunkt des Geschehens verschiebt sich hin zum leichteren Team. Das Modell erfasst diese „Verschiebung" perfekt.

4. Der Trick der „Massenkorrektur"

Das ursprüngliche Modell von 1965 hatte eine vereinfachte Version, die leicht zu lösen war, aber einen Fehler aufwies: Es verstieß gegen das Gesetz der Massenerhaltung.

• Das Problem: Wenn man einfach die Mathematik verwendet, ist es wie ein Ballon, der beim Aufblähen magisch Luft gewinnt oder verliert. Die Gesamtmenge an „Stoff" (Masse) addiert sich nicht korrekt.
• Die Lösung: Die Autoren erkannten, dass, wenn man das gesamte Mischprofil einfach nur leicht zur Seite des leichten Fluids verschiebt, die Mathematik plötzlich perfekt funktioniert.
  • Die Analogie: Stellen Sie sich einen perfekt symmetrischen Sandhügel vor (das vereinfachte Modell). Er sieht schön aus, aber wenn man den Sand wiegt, fehlt ein wenig. Wenn man den ganzen Hügel ein paar Zentimeter nach links schiebt, gleicht sich das Gewicht aus, und er sieht plötzlich genau wie die chaotischen Realwelt-Daten aus.
  • Diese „Verschiebung" erklärt, warum die Spitzen schneller wachsen als die Blasen: Die Diffusion des „Logarithmus der Dichte" ist symmetrisch, aber die Notwendigkeit, die Masse zu erhalten, zwingt die gesamte Struktur, sich zur leichten Seite hin zu neigen.

5. Das Fazit

Der Artikel kommt zu dem Schluss, dass dieses einfache, eindimensionale Modell von 1965 tatsächlich ein „Goldmine" ist.

• Es erfasst all die seltsamen, komplexen Verhaltensweisen der Mischung bei hohen Dichten (Asymmetrie, verschiebende Geschwindigkeiten, logarithmisches Wachstum), die moderne Wissenschaftler erst kürzlich mit Supercomputern bestätigt haben.
• Es legt nahe, dass die Physik dieser Turbulenz durch einen Wettbewerb zwischen Diffusion (Ausbreitung) und Massenerhaltung (Beibehaltung der Gesamtmenge an Fluid) bestimmt wird.

Zusammenfassend: Die Autoren gruben eine alte, staubige Theorie aus, putzten sie auf und zeigten, dass sie moderne Beobachtungen der Fluidmischung besser erklärt als viele aktuelle Theorien. Sie bewiesen, dass eine einfache „Verschiebung" in der Mathematik die Fehler des alten Modells korrigiert und perfekt beschreibt, warum schwere Fluide schneller absinken als leichte Fluide aufsteigen, wenn sie sehr unterschiedliche Dichten haben.

Technische Zusammenfassung

1. Problemstellung

Die Rayleigh–Taylor (RT)-Instabilität tritt auf, wenn eine schwere Flüssigkeit in eine leichtere Flüssigkeit beschleunigt wird, was zu turbulenter Durchmischung führt. Während das späte Wachstum der Höhe der Mischschicht ($h$) als gut etabliert gilt und mit $h \sim g t^2$ skaliert, bleibt die Abhängigkeit vom Dichteverhältnis ($R = \rho_H / \rho_L$) bei Strömungen mit variabler Dichte (nicht-Boussinesq) Gegenstand von Debatten.

• Aktueller Konsens: Die meisten modernen Studien interpretieren Ergebnisse unter Verwendung der Boussinesq-Näherung, bei der die Wachstumsrate linear mit der Atwood-Zahl ($A = (R-1)/(R+1)$) skaliert, d. h. $h \approx \alpha A g t^2$.
• Die Diskrepanz: Experimentelle und numerische Daten zeigen, dass das Blasenwachstum relativ unempfindlich gegenüber $A$ ist, während das Spike-Wachstum mit steigendem $A$ signifikant zunimmt, was zu Asymmetrien in Geschwindigkeits- und Dichteprofilen führt.
• Die Lücke: Eine wegweisende Arbeit von Belen'kii und Fradkin aus dem Jahr 1965 schlug ein Modell für turbulente Diffusivität vor, das eine Skalierung von $h \propto (\ln R) g t^2$ ergibt. Diese Arbeit wurde jedoch in Russisch veröffentlicht, enthielt typografische Fehler und wurde weitgehend übersehen. Darüber hinaus stützte sich die ursprüngliche Analyse auf eine vereinfachte gewöhnliche Differentialgleichung (ODE), ohne die Eigenschaften der vollständigen Ähnlichkeitsgleichung vollständig zu untersuchen oder sie mit modernen hochgenauen Daten zu vergleichen.

2. Methodik

Die Autoren nehmen das Modell von Belen'kii und Fradkin (1965) erneut auf, klären es und erweitern es unter Verwendung eines rigorosen eindimensionalen (1D) Rahmens.

Mathematischer Rahmen:

• Steuernde Gleichungen: Die Autoren beginnen mit den Reynolds-gemittelten Navier-Stokes-Gleichungen (RANS) für Strömungen mit variabler Dichte. Sie leiten eine mittlere Dichtegleichung ab, die analog zu einem 1D-Diffusionsproblem ist:
  $$\frac{\partial \bar{\rho}}{\partial t} = \frac{\partial}{\partial y} \left( D_t \frac{\partial \bar{\rho}}{\partial y} \right)$$
• Modell für turbulente Diffusivität: Sie leiten eine spezifische Form für die turbulente Diffusivität ($D_t$) basierend auf der Prandtl-Mischwegtheorie und Energie-Skalierungsargumenten ab. Das Modell postuliert:
  $$D_t^* = K \left( \frac{1}{\bar{\rho}} \frac{\partial \bar{\rho}}{\partial y} \right)^{1/2} g^{1/2} h_t^2$$
  wobei $h_t$ die Turbulenz-Höhenskala ist.
• Selbstähnliche Lösung: Durch Einführung von Ähnlichkeitsvariablen wird die partielle Differentialgleichung auf eine nichtlineare ODE reduziert.
  • Vollständige ODE: Eine nichtlineare ODE dritter Ordnung (reduziert auf eine ODE erster Ordnung für eine transformierte Variable $x$), die die gesamte Physik erfasst.
  • Vereinfachte ODE: Eine analytische Näherung, die durch Vernachlässigung eines Terms höherer Ordnung ($x^3$) abgeleitet wurde und für kleine Dichteverhältnisse ($R \lesssim 4$) gültig ist.

Analysestrategie:

• Numerische Lösung: Die vollständige ODE wird numerisch für einen weiten Bereich von Atwood-Zahlen gelöst ($A \in [0, 0.8]$).
• Analytische Lösung: Die vereinfachte ODE wird analytisch gelöst, um die ursprüngliche $\ln R$-Skalierung wiederherzustellen.
• Massenkorrektur: Die Autoren identifizieren, dass die vereinfachte ODE aufgrund von Symmetrieannahmen die globale Massenerhaltung verletzt. Sie schlagen eine „Massenkorrektur" (räumliches Verschieben der Lösung) vor, um das vereinfachte Modell mit der physikalischen Realität in Einklang zu bringen.
• Validierung: Theoretische Profile (Dichte, Diffusivität, Wachstumsraten) werden mit Daten aus Direkter Numerischer Simulation (DNS) und experimentellen Ergebnissen aus verschiedenen Literaturquellen verglichen.

3. Hauptbeiträge

1. Wiederbelebung und Klärung: Das Papier liefert eine saubere, moderne Herleitung des Modells von Belen'kii und Fradkin aus dem Jahr 1965, korrigiert Notationen und klärt Annahmen, die zuvor verschleiert waren.
2. Analyse der vollständigen ODE: Im Gegensatz zur ursprünglichen Arbeit löst dieses Papier die vollständige Ähnlichkeitsgleichung (nicht nur die Näherung). Es zeigt, dass das vollständige Modell inhärent komplexe nicht-Boussinesq-Eigenschaften erfasst, ohne auf ad-hoc-Anpassungen angewiesen zu sein.
3. Mechanismus für Asymmetrie: Die Autoren liefern eine mechanistische Erklärung für die beobachtete Asymmetrie zwischen Spikes und Blasen. Sie zeigen, dass die natürliche Variable für turbulente Durchmischung in Strömungen mit variabler Dichte $\ln \bar{\rho}$ ist. Die Diffusion von $\ln \bar{\rho}$ ist symmetrisch, aber die exponentielle Rückabbildung auf die Dichte ($\bar{\rho}$) kombiniert mit der Einschränkung der Massenerhaltung erzwingt eine räumliche Verschiebung zur Seite der leichten Flüssigkeit.
4. Massenkorrigiertes vereinfachtes Modell: Sie zeigen, dass die einfache analytische Lösung (parabolische Diffusivität) die vollständige Lösung genau reproduzieren kann, wenn eine globale Massenkorrektur (räumliche Verschiebung) angewendet wird. Dies bietet eine rechnerisch günstige, geschlossene Näherung für die Turbulenzmodellierung.
5. Validierung des Skalierungsgesetzes: Die Studie validiert die $\ln R$-Skalierung für die Höhe der Mischschicht und zeigt, dass sie als „Median"-Schätzung fungiert, die die Streuung vorhandener experimenteller und numerischer Daten besser beschreibt als die lineare $A$-Skalierung bei hohen Dichteverhältnissen.

4. Hauptergebnisse

• Skalierung der Mischschichthöhe: Die gesamte Höhe der Mischschicht skaliert als $h \propto (\ln R) g t^2$. Obwohl bei extremen $A$-Werten Abweichungen auftreten, bleibt diese Skalierung mit der Streuung in verfügbaren DNS- und experimentellen Daten bis zu $A \approx 0.84$ konsistent.
• Asymmetrie Spike vs. Blase:
  • Das Modell sagt voraus, dass die Spike-Höhen ($h_s$) mit steigendem $A$ schneller wachsen als die Blasen-Höhen ($h_b$).
  • Das Verhältnis $h_s/h_b$ nimmt mit $A$ zu, was mit empirischen Beobachtungen übereinstimmt.
  • Der Wachstumsparameter $\alpha_s$ nimmt mit $A$ zu, während $\alpha_b$ relativ konstant bleibt (oder weniger abweicht), was erklärt, warum die Boussinesq-Näherung für Spikes bei hohem $A$ versagt.
• Profilverschiebungen:
  • Geschwindigkeit/Diffusivität: Die Profile der turbulenten Diffusivität (und Geschwindigkeit) verschieben sich systematisch zur Seite der leichten Flüssigkeit, wenn $A$ zunimmt.
  • Dichte: Die Molenbruchprofile kollabieren bei Normalisierung gut, aber die physikalischen Grenzen (Spitzen von Spikes/Blasen) werden asymmetrisch ($\lambda_s > \lambda_b$).
• Vergleich mit Daten: Die theoretischen Profile für normalisierte Diffusivität und Molenbruch zeigen eine vernünftige Übereinstimmung mit DNS-Daten (z. B. von Goh et al., Livescu et al.), insbesondere im Kern der Mischschicht. Diskrepanzen beschränken sich hauptsächlich auf die Ränder, wo molekulare Diffusion (die im Modell vernachlässigt wird) relevant wird.

5. Bedeutung

• Theoretische Einsicht: Das Papier löst eine langjährige Unklarheit in der RT-Turbulenz, indem es zeigt, dass die konkurrierenden Dynamiken zwischen der Diffusion von $\ln \bar{\rho}$ (die symmetrisch ist) und der Massenerhaltung (die eine Verschiebung erfordert) die fundamentalen Treiber der nicht-Boussinesq-Asymmetrie sind.
• Praktischer Nutzen: Die „massenkorrigierte" vereinfachte Lösung liefert einen Satz geschlossener Ausdrücke für Dichte- und Diffusivitätsprofile. Dies ist für die Entwicklung von Reynolds-gemittelten Turbulenzmodellen (RANS) von großem Wert, bei denen Recheneffizienz entscheidend ist.
• Neubewertung der Skalierung: Die Arbeit stellt die Dominanz der linearen Atwood-Zahl-Skalierung ($A$) zugunsten der logarithmischen Dichteverhältnis-Skalierung ($\ln R$) für RT-Strömungen mit variabler Dichte in Frage und legt nahe, dass letztere eine robustere theoretische Grundlage für Anwendungen mit hohem Dichteverhältnis ist (z. B. Trägheitseinschlussfusion, Astrophysik).

Zusammenfassend zeigen Goh und Blanquart, dass ein einfaches 1D-theoretisches Rahmenwerk, das ursprünglich 1965 vorgeschlagen wurde, die wesentliche Physik enthält, um moderne Beobachtungen von RT-Turbulenz mit variabler Dichte zu erklären, sofern das Modell vollständig gelöst oder für die Massenerhaltung korrigiert wird.