Scaling laws and local enhancements of buoyancy flux in stratified turbulent flows

Durch umfangreiche direkte numerische Simulationen geschichteter turbulenter Strömungen zeigt diese Studie auf, dass der Auftriebsfluss eine starke Intermittenz und nicht-gaußsche Statistiken aufweist, die durch großskalige Langzeitfluktuationen getrieben werden, wobei sein domänengemittelte Verhalten logarithmisch mit der Auftriebs-Reynolds-Zahl skaliert und fundamental mit konvektiven Instabilitäten verknüpft ist, welche burstsartige Energiedissipationszyklen auslösen.

Das Wesentliche

Stellen Sie sich die Atmosphäre und den Ozean nicht als glatte, fließende Flüsse vor, sondern als chaotische, brodelnde Suppen. Manchmal sind diese Suppen geschichtet wie ein Parfait – warm oben, kalt unten. Diese Schichtung wird als Stratifizierung bezeichnet. Normalerweise würde man denken, dass dies die Suppe ruhig und geordnet macht, wie einen stillen See. Doch dieses Paper enthüllt ein überraschendes Geheimnis: Selbst in diesen ruhigen, geschichteten Fluiden gibt es plötzliche, heftige „Ausbrüche“ von Chaos, die in riesigen Dimensionen auftreten, nicht nur in winzigen.

Hier ist eine Aufschlüsselung dessen, was die Forscher herausgefunden haben, unter Verwendung einfacher Analogien:

1. Die „Überraschungsparty“ in einem ruhigen Raum

In normaler, ungeschichteter Turbulenz (wie bei einem Mixer, der einen Smoothie mischt) findet Chaos hauptsächlich auf der sehr kleinen Skala statt – winzige kleine Wirbel. Aber in geschichteten Fluiden (wie dem Ozean oder der oberen Atmosphäre) fanden die Forscher heraus, dass Chaos auch groß sein kann.

Denken Sie an eine ruhige Bibliothek. Normalerweise flüstern die Leute. Aber plötzlich bricht ein massiver, unerwarteter Schrei aus einer bestimmten Ecke hervor. Dieses Paper zeigt, dass diese „Schreie“ (plötzliche Ausbrüche vertikaler Bewegungen und Temperaturveränderungen) in den Ozeanen und der Atmosphäre häufig vorkommen. Es sind nicht nur winzige Kräuselungen; es sind riesige, lokalisierte Ereignisse, die so groß wie der gesamte Fluss selbst sein können.

2. Der „Stau“ von Hitze und Bewegung

Die Wissenschaftler haben sich speziell auf den Auftriebsfluss konzentriert. Nennen wir dies den „Hitze-Bewegungs-Handschlag“. Er misst, wie viel Wärme gleichzeitig auf oder ab fließt, während Luft oder Wasser gleichzeitig auf oder ab steigt.

• Die Entdeckung: Sie fanden heraus, dass dieser „Handschlag“ unglaublich erratisch ist. Manchmal schießt heiße Luft gewaltsam nach oben, und manchmal stürzt kalte Luft gewaltsam nach unten.
• Die Analogie: Stellen Sie sich eine belebte Autobahn vor. Meistens bewegen sich die Autos mit einer stetigen Geschwindigkeit. Aber in dieser Studie fanden sie heraus, dass gelegentlich ein massiver Stau entsteht, bei dem Autos (Hitze und Bewegung) plötzlich beschleunigen oder stoppen, auf eine Weise, die völlig unvorhersehbar und extrem ist. Diese Ereignisse sind so extrem, dass die statistischen „Enden“ (Tails) der Daten riesig sind – das heißt, die „seltsamen“ Ereignisse sind viel häufiger, als die Mathematik normalerweise vorhersagt.

3. Die „Goldlöckchen-Zone“ des Chaos

Die Forscher testeten viele verschiedene Bedingungen und veränderten dabei, wie „dickflüssig“ das Fluid ist (Viskosität) und wie stark die Schichtung ist (Stratifizierung). Sie fanden eine spezifische „Goldlöckchen-Zone“ (Sweet Spot), in der diese extremen Ausbrüche am häufigsten auftreten.

• Zu viel Schichtung (Starke Stratifizierung): Das Fluid ist wie ein steifes Gel. Es vibriert nur wie eine Gitarrensaite (Wellen) und mischt sich nicht viel.
• Zu wenig Schichtung (Schwache Stratifizierung): Das Fluid ist wie ein chaotischer Mixer. Es vermischt alles, aber die Ausbrüche sind nicht so ausgeprägt.
• Genau richtig (Der Sweet Spot): Wenn die Schichtung moderat ist, wird das Fluid instabil. Es ist wie ein Stapel Jenga-Blöcke, der fast stabil ist, aber bereit, zusammenzustürzen. In dieser Zone sind die „Schreie“ (die extremen Ausbrüche) am lautesten.

4. Der Mechanismus der „Energieschuld“

Warum passieren diese Ausbrüche? Das Paper schlägt einen einfachen Mechanismus vor: Energieschuld.

Stellen Sie sich vor, Sie haben zwei Bankkonten: eines für „Auf-und-Ab-Bewegung“ (kinetische Energie) und eines für „Hitze-Potenzial“ (potenzielle Energie).

• In einer perfekten Welt wären diese Konten ausgeglichen.
• In diesen Strömungen geraten sie aus dem Gleichgewicht. Das Konto für „Auf-und-Ab-Bewegung“ wird im Vergleich zum „Hitze“-Konto zu hoch.
• Die Natur hasst Schulden. Um dieses Ungleichgewicht zu beheben, lässt das System diese zusätzliche Energie plötzlich in einem massiven, heftigen Ausbruch frei. Dies erzeugt einen Wirbel oder einen Luftzug, der die Schichten durchmischt, die Schuld begleicht, und dann beginnt der Zyklus von vorn.

5. Was das für das „Große Ganze“ bedeutet

Das Paper behauptet nicht, den Klimawandel direkt zu lösen oder Hurrikane vorherzusagen. Stattdessen liefert es ein neues Regelwerk dafür, wie sich diese Fluide verhalten.

• Die Regel: Die Intensität dieser chaotischen Ausbrüche folgt einem spezifischen mathematischen Muster (einem Potenzgesetz), basierend darauf, wie stark die Schichtung ist.
• Die Kernbotschaft: Selbst in einem stabilen, geschichteten Ozean oder Himmel können Sie nicht davon ausgehen, dass die Strömung glatt ist. Es gibt verborgene, massive „Hot Spots“ der Vermischung, die in Ausbrüchen geschehen. Wenn Sie versuchen zu modellieren, wie sich Hitze oder Schadstoffe durch den Ozean oder die Atmosphäre bewegen, müssen Sie diese plötzlichen, heftigen Spitzen berücksichtigen und nicht nur den Durchschnittsfluss.

Zusammenfassung

Dieses Paper ist wie die Entdeckung, dass ein ruhiger See gar nicht ruhig ist. Er ist voller verborgener, massiver Unterwasser-Explosionen, die in einer spezifischen „Goldlöckchen-Zone“ der Stabilität stattfinden. Diese Explosionen werden durch ein Ungleichgewicht zwischen Bewegung und Hitze angetrieben und folgen einem vorhersagbaren mathematischen Rhythmus. Das Verständnis dieses Rhythmus hilft uns zu erkennen, dass die Atmosphäre und die Ozeane weitaus „spitzer“ und unvorhersehbarer sind, als wir bisher angenommen haben.

Technische Zusammenfassung

Titel: Skalierungsgesetze und lokale Verstärkungen des Auftriebsflusses in stratifizierten turbulenten Strömungen

Problemstellung
Stratifizierte turbulente Strömungen, die in der Atmosphäre und den Ozeanen allgegenwärtig sind, weisen Intermittenz nicht nur auf kleinen Skalen, sondern auch auf großen Skalen auf, die mit der mittleren Strömung vergleichbar sind. Während voll entwickelte homogene isotrope Turbulenz (HIT) für bekannt ist, auf kleinen Skalen intermittent zu sein, gibt es Hinweise darauf, dass das Geschwindigkeitsfeld in stratifizierten Strömungen selbst nicht-gaußsche Statistiken und eine „großskalige“ Intermittenz aufweisen kann. Dieses Phänomen beinhaltet das Entstehen starker vertikaler Geschwindigkeitsentwürfe (Drafts), die als lokale Energieinjektionsmechanismen wirken und zu hoch variablen Transporteigenschaften in Richtung der Schwerkraft führen. Die spezifische Dynamik, welche diese großskaligen Ereignisse steuert, insbesondere hinsichtlich des Auftriebsflusses ($B_f$) und dessen Abhängigkeit von der Auftriebs-Reynoldszahl ($Ri_B$) und der Prandtl-Zahl ($Pr$), bedarf weiterer Klärung, um Vermischung und Dissipation in geophysikalischen Kontexten zu verstehen.

Methodik
Die Autoren führten eine umfassende parametrische Untersuchung mittels direkter numerischer Simulationen (DNS) der Boussinesq-Gleichungen durch. Die Studie nutzte den GHOST Pseudo-Spektral-Code mit hybrider Parallelisierung (MPI, OpenMP, CUDA).

• Gleichungen: Die unkomprimierbaren Navier-Stokes-Gleichungen, gekoppelt mit einer Temperaturgleichung, einschließlich Auftriebstermen, die durch die Brunt-Väisälä-Frequenz ($N$) charakterisiert sind.
• Parameter: Die Studie variierte die Froude-Zahl ($Fr$) über einen geophysikalisch relevanten Bereich ($0,01 \le Fr \le 1$), was Auftriebs-Reynolds-Zahlen ($Ri_B$) von $0,06$ bis $2300$ entspricht. Zwei Prandtl-Zahlen wurden untersucht: $Pr = 1$ und $Pr = 6$.
• Forcing-Schemata: Es wurden zwei Arten des Forcings für das Geschwindigkeitsfeld eingesetzt: ein zufälliges Pouquet-Patterson-Forcing (PP) und ein Taylor-Green-Forcing (TG).
• Umfang der Simulation: Insgesamt wurden 60 Durchläufe analysiert, einschließlich Serien mit $Pr=1$ (PP-Pr1, TG-Pr1), $Pr=6$ (PP-Pr6) sowie passive Skalar-Simulationen ($N=0$). Die Auflösungen reichten von $128^3$ bis $512^3$ Gitterpunkten.
• Analyse: Die Studie konzentrierte sich auf Hochordnungsstatistiken (Skewness und Kurtosis) der vertikalen Geschwindigkeit ($w$), der Temperaturschwankungen ($\theta$) und des Auftriebsflusses ($B_f = Nw\theta$). Die Analyse beinhaltete Volumen- und Zeitmittelungen sowie die Untersuchung von Wahrscheinlichkeitsdichtefunktionen (PDFs) und gemeinsamen PDFs.

Wesentliche Beiträge und Ergebnisse

1. Großskalige Intermittenz im Auftriebsfluss:
   Die Studie bestätigt, dass der Auftriebsfluss ($B_f$) starke nicht-Gaußsche Tails aufweist, wobei die Kurtosis-Werte Werte von $\approx 10^2$ erreichen. Dies deutet darauf hin, dass stratifizierte geophysikalische Strömungen durch hochvariable Transporteigenschaften charakterisiert werden können, selbst unter stabiler Stratifizierung. Diese Intermittenz wird mit der langzeitigen, großskaligen intermittierenden Charakteristik der vertikalen Geschwindigkeit und der Temperatur in Verbindung gebracht.

2. Skalierungsgesetze und $Ri_B$-Abhängigkeit:

• Skewness: Die Skewness des Auftriebsflusses ($S_{B_f}$) steigt mit zunehmender $Ri_B$ nach einem Potenzgesetz ($S_{B_f} \propto Ri_B^{0,28}$), wenn die Stratifizierung schwächer wird. Dieser Trend sättigt im Grenzfall des passiven Skalars ($N \to 0$), in dem die Temperatur vom Geschwindigkeitsfeld entkoppelt ist.
• Kurtosis-Maxima: Die Maxima für die Kurtosis von $w$, $\theta$ und $B_f$ treten in einem engen Bereich von $1 < Ri_B < 2$ auf. Dieser Bereich entspricht einem spezifischen Regime des Wellen-Wirbel-Gleichgewichts (wave-eddy balance).
• Mittlerer Auftriebsfluss: Der domänengemittelte Auftriebsfluss $\langle B_f \rangle$ zeigt zwei unterschiedliche Trends: Er skaliert logarithmisch mit $Ri_B$ ($\langle B_f \rangle \propto \ln(Ri_B)$) im intermediären Regime und nähert sich einem kleinen Offset an, wenn die Stratifizierung stärker wird ($Ri_B \to 0$).

3. Rolle der Prandtl-Zahl:
   Die Variation der Prandtl-Zahl ($Pr = 1$ vs. $Pr = 6$) scheint die Entwicklung extremer vertikaler Entwürfe oder das Entstehen großskaliger Intermittenz innerhalb der betrachteten Werte nicht signifikant zu beeinflussen. Die Form der Skalierungskurven für die verschiedenen $Pr$-Serien ist nahezu identisch und unterscheidet sich lediglich in der Größenordnung.

4. Mechanismus extremer Ereignisse:
   Ein einfaches Modell für die zeitliche Entwicklung von Energie und Auftriebsfluss zeigt, dass der „Energiedefekt“ zwischen vertikaler kinetischer Energie und potenzieller Energie die starken Ereignisse im Auftriebsfluss antreibt.

• Das Modell (basierend auf der Vieillefosse-Dynamik, erweitert auf stratifizierte Strömungen) legt nahe, dass ein Mangel an Äquipartition zwischen vertikaler und potenzieller Energie den Auftriebsfluss anfällig für lokale Wellendynamiken macht.
• Dieser Energie-Ungleichgewicht führt zu konvektiven Instabilitäten, der Bildung von 2D- und 3D-Wirbeln und schneller Dissipation auf der Skala einer Umsetzungszeit (turnover time), was es dem energetischen Zyklus ermöglicht, in Ausbrüchen neu zu starten.
• Gemeinsame PDFs zeigen, dass die vertikale Geschwindigkeit am Peak der Intermittenz ($Ri_B \approx 1,78$) bis zu $15\sigma_w$ erreichen kann, während die horizontale Geschwindigkeit begrenzt bleibt, was die anisotrope Natur dieser Ereignisse hervorhebt.

5. Mischeffizienz:
   Die Studie analysiert die Mischeffizienz ($\langle B_f \rangle / \langle \epsilon_V \rangle$) und die irreversible Mischeffizienz ($\hat{\Gamma}$). Beide Größen zeigen eine Abnahme, wenn die Stratifizierung stärker wird, bis sie bei niedrigen $Ri_B$ Sättigung erreichen. Die Autoren stellen fest, dass extreme Ereignisse zwar die lokalen Flüsse signifikant beeinflussen, die globale Mischeffizienz jedoch nicht verändern, wenn sie über das Strömungsvolumen und die Zeit gemittelt wird. Dennoch können die lokalen Werte von $B_f$ zwischen Regionen mit extremen Ereignissen und solchen ohne um eine Größenordnung abweichen.

Bedeutung
Die Arbeit behauptet, dass die Identifizierung dieser Skalierungsgesetze und der spezifischen Bedingungen, die extreme Ereignisse auslösen, essenziell ist, um die Modellierung und Parametrisierung geophysikalischer und klimatologischer Strömungen zu verbessern.

• Vorhersagbarkeit: Die Ergebnisse legen nahe, dass die schnelle Intensivierung von Phänomenen wie Hurrikans besser vorhergesagt werden kann, indem man Variationen der vertikalen Geschwindigkeit verfolgt, gegeben die etablierte Verbindung zwischen vertikalen Entwürfen und großskaliger Intermittenz.
• Geophysikalische Relevanz: Der identifizierte Bereich von $Ri_B$ und $Fr$, in dem diese extremen Ereignisse auftreten, ist sowohl für die Atmosphäre als auch für die Ozeane von Bedeutung.
• Modellierung: Das Verständnis der lokalen Verstärkungen des Auftriebsflusses und des Energiedefekt-Mechanismus bietet eine Grundlage für die Verfeinerung von Sub-Grid-Scale-Modellen und Parametrisierungen, die in Klima- und Ozeanmodellen verwendet werden, insbesondere bezüglich der Frage, wie Turbulenz und Wellen interagieren, um Vermischung voranzutreiben.

Die Autoren kommen zu dem Schluss, dass, obwohl die globale Mischeffizienz relativ stabil bleiben mag, die lokale, stoßartige Natur der Dissipation und des Transports, die durch diese großskaligen intermittierenden Entwürfe getrieben werden, ein kritisches Merkmal der stratifizierten Turbulenz ist, das in geophysikalischen Anwendungen berücksichtigt werden muss.