Evolution of passive scalar mixing layers in stratified and unstratified homogeneous turbulence

Hochauflösende Large-Eddy-Simulationen zeigen, dass sich das passive skalare Mischen in geschichteter homogener Turbulenz in transversaler Richtung ähnlich wie in ungeschichteten Fällen verhält, die Schichtung jedoch die vertikale Durchmischung durch Unterdrückung großskaliger Durchmischungsprozesse stark hemmt, was zu ausgeprägten Wachstumsgrenzen und Fluktuationsintensitäten führt, die spezifische Modellierungsansätze für skalare Flüsse begründen.

Das Wesentliche

Stellen Sie sich vor, Sie beobachten, wie ein Tintentropfen sich in einem Glas Wasser ausbreitet. Wenn Sie das Wasser sanft umrühren, verteilt sich die Tinte gleichmäßig. So denken Wissenschaftler gewöhnlich über das Mischen in Fluiden nach: Turbulenz wirkt wie ein riesiger Löffel, der alles zusammenrührt, bis es einheitlich ist.

Was passiert jedoch, wenn dieses Wasser geschichtet ist? Stellen Sie sich vor, das Wasser am Boden ist schwer und salzig, während das Wasser oben leicht und frisch ist. Dies nennt man Schichtung. In der realen Welt geschieht dies im Ozean (wo das tiefe Wasser dichter ist) und in der Atmosphäre (wo die Luft mit zunehmender Höhe dünner und leichter wird).

Diese Arbeit ist ein hochtechnisiertes Computerexperiment, das folgende Frage stellt: Wie verändert diese Schichtung die Art und Weise, wie ein „Fleck" (wie Verschmutzung oder Rauch) sich durch ein turbulentes Fluid ausbreitet?

Hier ist die Geschichte dessen, was sie herausfanden, aufgeschlüsselt in einfache Konzepte.

Das Setup: Zwei Arten von „Flecken"

Die Forscher erzeugten ein virtuelles Fluid, das anfangs chaotisch wirbelte (turbulent). Dann brachten sie zwei verschiedene „Flecken" (passive Skalare) ein, um zu sehen, wie sie sich verhielten:

1. Der horizontale Fleck: Eine Schicht Tinte, die sich seitlich ausbreitete, wie ein flaches Blatt Papier, das im Wasser schwebt.
2. Der vertikale Fleck: Eine Schicht Tinte, die sich auf und ab ausbreitete, wie eine vertikale Wand aus Farbe.

Sie führten zwei Simulationen durch: eine in normalem Wasser (ungeschichtet) und eine in „geschichtetem" Wasser (geschichtet).

Die große Entdeckung: Der Unterschied zwischen „Oben-Unten" und „Seitlich-Seitlich"

1. Der seitliche Fleck (transversale Schicht)

Was passierte: Als der Fleck seitlich ausgebreitet wurde, hinderte die Schichtung ihn nicht daran, sich auszubreiten. Tatsächlich breitete er sich im geschichteten Wasser sogar etwas schneller aus als im normalen Wasser.
Die Analogie: Stellen Sie sich eine Menschenmenge vor, die in einem Flur rennt. Wenn der Boden völlig flach ist (ungeschichtet), rennen sie in alle Richtungen. Wenn der Boden eine leichte, unsichtbare Neigung hat (geschichtet), rennen sie immer noch genauso gut seitwärts, vielleicht sogar etwas energischer. Der „Fleck" breitet sich in beiden Fällen breit aus.
Der Haken: Obwohl die Gesamtausbreitung ähnlich war, war der „Fleck" im geschichteten Wasser „spitzer". Statt eines sanften Gradienten hatte er schärfere, gezacktere Ränder. Er war „intermittierender", was bedeutet, dass es Taschen aus reiner Tinte und Taschen aus reinem Wasser gab, mit weniger sanftem Übergang dazwischen.

2. Der vertikale Fleck (vertikale Schicht)

Was passierte: Hier geschah das Magische (und die Einschränkung). Im normalen Wasser breitete sich der vertikale Fleck genauso leicht auf und ab aus wie der seitliche. Aber im geschichteten Wasser hörte die Ausbreitung fast vollständig auf.
Die Analogie: Stellen Sie sich vor, Sie versuchen, einen dicken Milchshake mit einem Löffel umzurühren. Wenn Sie versuchen, den Löffel auf und ab zu bewegen, widerstehen die Schichten des Shakes Ihnen. Die schweren Bestandteile wollen unten bleiben, und die leichten Bestandteile wollen oben bleiben. Die „Rühr"-Bewegung wird zusammengedrückt.
Das Ergebnis: Der vertikale Fleck wuchs am Anfang ein winziges bisschen, stieß dann aber an eine „Decke". Er konnte nicht weiter wachsen, weil die stabilen Schichten des Fluids wie ein Deckel wirkten und verhinderten, dass die Turbulenz die Dinge vertikal mischte. Das Fluid konnte immer noch seitwärts wirbeln, aber es konnte nicht auf und ab mischen.

Warum ist das wichtig? (Das „Warum" hinter dem „Was")

Die Forscher stellten fest, dass sich das Fluid in vertikaler Richtung wie eine Feder verhält. Sobald die Turbulenz versucht, eine schwere Schicht nach oben zu drücken oder eine leichte nach unten, zieht die Schwerkraft sie zurück. Dies stoppt die „Rühr"-Bewegung.

Das Fluid kann jedoch immer noch seitwärts wirbeln. Daher wird die „vertikale Länge" der Turbulenz auf eine bestimmte Größe festgelegt (bestimmt durch die Stärke der Schwerkraft und der Fluidschichten), und der Fleck kann diese Größe nicht überschreiten.

Das „Rezept" zur Vorhersage

Die Arbeit versuchte auch, ein einfaches mathematisches „Rezept" zu erstellen, um vorherzusagen, wie sich diese Flecken ausbreiten würden, ohne einen Supercomputer zu benötigen.

• Wenn Sie die Form des Flecks kennen: Sie können eine einfache Ein-Zahlen-Formel verwenden, um vorherzusagen, wie schnell er sich seitwärts ausbreitet. Das funktioniert sehr gut.
• Wenn Sie die Form nicht kennen: Sie müssen die Form raten (unter der Annahme, dass sie wie eine glatte Kurve aussieht). Wenn Sie dies tun, benötigen Sie eine Zwei-Zahlen-Formel. Dies funktioniert hervorragend, nachdem der Fleck Zeit hatte, sich mit dem wirbelnden Fluid in einen Rhythmus zu finden.

Das Fazit

• Seitliches Mischen: Stabile Schichten (wie im tiefen Ozean oder in der oberen Atmosphäre) verhindern das seitliche Mischen nicht; sie machen es möglicherweise sogar etwas intensiver und gezackter.
• Vertikales Mischen: Stabile Schichten wirken wie eine Bremse. Sie verhindern, dass sich das Fluid fast vollständig auf und ab mischt.
• Der Unterschied zwischen „Rühren" und „Mischen": Das Fluid kann immer noch seitwärts „rühren" (sich bewegen), aber es kann nicht vertikal „mischen" (vermischen), weil die Schichten einem Austausch widerstehen.

Die Autoren weisen darauf hin, dass ihr Experiment eine spezifische Art von Fluideigenschaft verwendete (Prandtl-Zahl von 0,7). Sie warnen, dass sich die Ergebnisse ändern könnten, wenn das Fluid „dicker" wäre oder andere Eigenschaften hätte (Prandtl-Zahl > 1), aufgrund eines „umgekehrten" Effekts, bei dem das Mischen seinen eigenen Auftrieb erzeugt. Aber für die Bedingungen, die sie testeten, gilt die Regel: „Seitwärts ist frei, auf-und-ab ist blockiert".

Technische Zusammenfassung

Technische Zusammenfassung: Evolution passiver Skalar-Mischschichten in geschichteter und ungeschichteter homogener Turbulenz

Problemstellung
Die Dispersion von Partikeln und Skalaren in stabil geschichteten Umgebungen, wie der oberen Troposphäre, der Stratosphäre und ozeanischen Mischschichten, stellt eine kritische Herausforderung in der atmosphärischen und ozeanischen Wissenschaft dar. Während die Auswirkungen stabiler Schichtung auf das Dichtefeld und den Turbulenzzerfall gut dokumentiert sind, bleibt die Dispersion passiver Skalare (z. B. Aerosole, Schadstoffe, Mikroplastik) weniger verstanden. Insbesondere ist unklar, wie sich eine passive Skalar-Mischschicht (PSML) entwickelt, wenn sie in homogene isotrope Turbulenz (HIT) eingebracht wird, die anschließend unter dem Einfluss eines stabilisierenden Dichtegradienten zerfällt. Die zentrale Frage lautet, wie die Einführung von Auftriebslängen und Anisotropie die selbstähnliche Entwicklung von Skalargestatistiken verändert, die in ungeschichteter HIT beobachtet werden.

Methodik
Die Autoren setzten hochauflösende Large-Eddy-Simulationen (LES) ein, um den Zerfall sowohl ungeschichteter als auch stabil geschichteter homogener Turbulenz zu untersuchen. Die Simulationen nutzten die nicht-divergenten Navier-Stokes-Gleichungen mit der nicht-hydrostatischen Boussinesq-Näherung.

• Strömungskonfiguration: Zwei Simulationen wurden durchgeführt. Die erste war ungeschichtete HIT. Die zweite war in den Anfangsbedingungen identisch, enthielt jedoch einen gleichmäßigen stabilisierenden Dichtegradienten ($g \neq 0$), wodurch die Strömung nach etwa einer Auftriebsperiode in stark geschichtete Turbulenz überging.
• Passive Skalare: Zwei passive Skalar-Mischschichten wurden initialisiert: eine mit einem mittleren Gradienten in vertikaler Richtung ($\xi_z$) und eine in transversaler Richtung ($\xi_y$). Diese wurden als Mischungsbrüche ($0 \le \xi \le 1$) modelliert.
• Numerische Details: Um diffusive Effekte zu minimieren, die Abweichungen von der Selbstähnlichkeit verursachen, verwendeten die Simulationen Hyperviskosität und Hyperdiffusivität mit gleichen Koeffizienten. Die molekulare Prandtl-Zahl wurde auf $Pr = 0.7$ gesetzt. Die Reynolds-Zahlen waren hoch ($Re_h > 15.000$ und $Re_\lambda > 450$ für den ungeschichteten Fall), was die Existenz von inertial-konvektiven und inertial-diffusiven Bereichen sicherstellte. Die Auflösung war ausreichend, um die Auftriebslängenskala aufzulösen ($L_b/\Delta \gg 1$).
• Analyse: Die Studie verfolgte die Entwicklung von Skalarenprofilbreiten, Varianz, höheren Momenten (Schiefe, Kurtosis), Flüssen und Mischungs-Effizienz. Zudem wurden Eddy-Diffusivitätsmodelle für den transversalen Skalarfluss unter zwei Bedingungen bewertet: wenn das mittlere Skalarenprofil bekannt ist und wenn es angenommen werden muss.

Hauptergebnisse

1. Transversale Mischung ($\xi_y$): In der Richtung senkrecht zur Schwerkraft ist die Entwicklung des passiven Skalars im geschichteten Fall im Großen und Ganzen ähnlich wie im ungeschichteten Fall, breitet sich jedoch etwas schneller aus.

   • Die Skalarfluktuationen sind im geschichteten Fall intensiver, mit einer Spitzenvarianz von etwa 15 % höher.
   • Die turbulente/nicht-turbulente Grenzfläche ist im geschichteten Fall intermittenter.
   • Sobald das Skalarfeld ein Quasi-Gleichgewicht mit dem Geschwindigkeitsfeld erreicht, kollabieren die normalisierten Flussprofile für beide Fälle (geschichtet und ungeschichtet), was auf ähnliche zugrunde liegende Transportmechanismen trotz der Schichtung hinweist.

2. Vertikale Mischung ($\xi_z$): Die Schichtung hemmt die vertikale Mischung drastisch.

   • Die vertikale Schichtbreite ($\delta_z$) wächst zunächst, hört aber nach etwa einer Auftriebsperiode ($t^* \approx 3$) schnell auf, sich auszudehnen.
   • Dieses Aussetzen des Wachstums wird der Unterdrückung großskaliger vertikaler Durchmischung durch Auftriebskräfte zugeschrieben.
   • Die endgültige Breite der Schicht ist proportional zur vertikalen Integrallängenskala der horizontalen Geschwindigkeit ($L_v$), die so begrenzt ist, dass eine vertikale Froude-Zahl der Größenordnung eins aufrechterhalten wird ($Fr_v \sim O(1)$).
   • Während das mittlere Profil aufgrund molekularer Diffusion ein monotonisches Wachstum zeigt, offenbart die Flussanalyse ein Phänomen des „Rückwärtsflusses", bei dem Durchmischung gelegentlich Skalar-Fasern zurück auf ihre ursprüngliche Seite fegt, was mit wellenartigem Verhalten im Inneren konsistent ist.

3. Mischungs-Effizienz: Die kumulative Mischung (Dissipation der Skalarenvarianz) im geschichteten Fall übersteigt bei späten Zeiten schließlich die des ungeschichteten Falls. Dies wird der Tatsache zugeschrieben, dass der Mischung dem Durchmischen folgt, und die horizontale Durchmischungsrate in der geschichteten Strömung proportional zur horizontalen Längenskala ist, die größer wird als im ungeschichteten Fall.

4. Modellierung des Flusses:

   • Bekanntes Profil: Wenn das mittlere Skalarenprofil bekannt ist, sagt ein Eddy-Diffusivitätsmodell mit einer einzigen Konstante ($F_y = -c_1 \delta_y v' \partial \xi / \partial y$) den Fluss für beide Fälle (geschichtet und ungeschichtet) effektiv voraus, mit $c_1 \approx 0.4$.
   • Angenommenes Profil: Wenn die Profilform (angenommen als Fehlerfunktion) und die Skalarlängenskala aus Geschwindigkeitsstatistiken geschätzt werden müssen, ist ein Zwei-Konstanten-Modell erforderlich. Dieses Modell funktioniert bei späten Zeiten gut, wenn das Skalarfeld im Quasi-Gleichgewicht mit dem Geschwindigkeitsfeld steht, wodurch die Skalarlängenskala aus der Dissipationsrate der kinetischen Energie skaliert werden kann.

Bedeutung und Behauptungen
Die Arbeit behauptet, eine detaillierte Charakterisierung der passiven Skalar-Mischung in zerfallender geschichteter Turbulenz zu liefern und die Lücke zwischen theoretischen Vorhersagen anisotroper Diffusion und numerischen Beobachtungen zu schließen.

• Anisotropie: Die Studie bestätigt, dass Schichtung eine scharfe Dichotomie in der Mischung schafft: Die horizontale Mischung bleibt aktiv und wird sogar verstärkt, während die vertikale Mischung effektiv unterdrückt wird, sobald die Strömung stark geschichtet wird.
• Modellierungsnutzen: Die Autoren zeigen, dass Standard-Eddy-Diffusivitätsansätze für geschichtete Strömungen weiterhin tragfähig sind, sofern das Skalar im Quasi-Gleichgewicht mit dem Geschwindigkeitsfeld steht. Sie betonen, dass die Annahme eines bekannten Profils die Modellierungsanforderungen im Vergleich zur Annahme einer Profilform erheblich vereinfacht.
• Einschränkungen und zukünftige Arbeit: Die Autoren vermerken bescheiden, dass ihre Ergebnisse spezifisch für eine Prandtl-Zahl von 0,7 sind. Sie stellen ausdrücklich fest, dass neuere Literatur darauf hindeutet, dass höhere Prandtl-Zahlen ($Pr > 1$) bei kleinen Skalen umgekehrte Auftriebsflüsse induzieren, was die Strömungszerfallraten und folglich die Raten der passiven Skalar-Mischung verändern könnte. Daher sind die aktuellen Erkenntnisse ohne weitere Untersuchungen nicht direkt auf Szenarien mit hoher Prandtl-Zahl anwendbar.

Die Studie kommt zu dem Schluss, dass die Schichtung zwar die vertikale Dynamik einer Mischschicht grundlegend verändert, die transversale Mischung jedoch viele Merkmale ungeschichteter Turbulenz beibehält, wenn auch mit erhöhter Intermittenz und Fluktuationsintensität.