Two pathways to diapycnal mixing in strongly stratified flows with no initial vertical shear

Diese Studie kombiniert lineare Theorie und direkte numerische Simulationen, um aufzuzeigen, dass in stark geschichteten Strömungen ohne anfängliche vertikale Scherung horizontale Scherinstabilitäten unweigerlich die diapyknische Vermischung über zwei unterschiedliche Pfade vorantreiben – entweder durch das direkte Entstehen vertikaler Scherung oder durch eine nichtlineare Entwicklung hin zu säulenförmigen Wirbeln –, wobei beide letztlich kleine Kelvin-Helmholtz-Instabilitäten auslösen, jedoch aufgrund der Anregung unterschiedlicher vertikaler Skalen verschiedene Vermischungseffizienzen aufweisen.

Das Wesentliche

Stellen Sie sich den Ozean oder die Atmosphäre als einen riesigen, geschichteten Kuchen vor. Die Schichten bestehen aus Fluiden mit unterschiedlichen Dichten (wie verschiedene Geschmacksrichtungen eines Kuchens), und sie mögen es nicht, sich leicht zu vermischen. Normalerweise untersuchen Wissenschaftler, was passiert, wenn man diese Schichten gleichzeitig seitlich und auf und ab (vertikale Scherung) drückt. Aber diese Arbeit stellt eine andere Frage: Was passiert, wenn man die Schichten nur seitlich drückt, ohne anfängliche Auf-und-Ab-Bewegungen, in einer sehr stark geschichteten Umgebung?

Die Forscher fanden heraus, dass die Natur selbst dann, wenn man mit einem vollkommen horizontalen Fluss beginnt, zwei verschiedene „Rezepte“ oder Wege hat, um Chaos (Turbulenz) und die Vermischung der Schichten zu erzeugen. Welches Rezept die Natur wählt, hängt ganz davon ab, wie man das Experiment „impft“ – also, welchen winzigen kleinen Anstoß man dem Fluid zu Beginn gibt.

Hier ist die Aufschlüsselung der beiden Wege unter Verwendung einfacher Analogien:

Das Setup: Der ruhige Fluss

Stellen Sie sich einen breiten, ruhigen Fluss vor, der horizontal fließt. Das Wasser ist geschichtet wie ein Stapel Pfannkuchen (starke Schichtung). Zuerst ist die Strömung glatt und zweidimensional (sie bewegt sich nur nach links und rechts, nicht auf und ab).

Weg 1: Der „Überfülltes Zimmer“-Effekt (Der direkte Weg)

Wie es beginnt: Man gibt dem Fluss einen winzigen, zufälligen Anstoß überall gleichzeitig (wie das Werfen einer Handvoll Konfetti in die Luft).
Was passiert:

1. Die Welle: Aufgrund der Schichten kräuselt sich das Fluid nicht einfach nur nach links und rechts; es beginnt sofort, in vielen verschiedenen Größen gleichzeitig auf und ab zu kräuseln. Denken Sie an eine Menschenmenge in einem Raum, in der alle gleichzeitig versuchen sich zu bewegen, was ein chaotisches, multidirektionales Drängeln erzeugt.
2. Die Scherung: Diese Wellen erzeugen sehr schnell starke vertikale Strömungen (Scherung).
3. Der Aufbruch: Diese vertikalen Strömungen werden so stark, dass sie aufbrechen und kleine, heftige Wirbel (wie winzige Wasserwirbel) erzeugen. Dies ist die „Kelvin-Helmholtz-Instabilität“, die aussieht wie die brechenden Wellen, die man sieht, wenn Wind über das Wasser weht.
   Das Ergebnis: Die Vermischung geschieht effizient. Da die Energie über viele verschiedene Größen von Wellen verteilt ist, ist die „Reibung“ (viskose Dissipation) geringer, was den Mischprozess relativ effizient macht.

Weg 2: Der „Synchronisierter Tanz“-Effeffekt (Der indirekte Weg)

Wie es beginnt: Man gibt dem Fluss einen sehr spezifischen, organisierten Anstoß (wie ein Dirigent, der einen Taktstock schwingt, um alle zu einem bestimmten Muster zu bewegen).
Wass passiert:

1. Der Wirbel: Anstatt chaotischer Wellen organisiert sich das Fluid in langen, vertikalen Säulen aus wirbelndem Wasser (wie riesige Tornados, die im Fluss stehen). Für eine lange Zeit bleibt die Strömung perfekt zweidimensional, bestehend aus diesen großen, wirbelnden Säulen.
2. Das Wackeln: Schließlich werden diese riesigen Säulen instabil. Sie beginnen auf eine sehr spezifische, hochfrequente Weise zu wackeln. Die Forscher nennen dies eine „hyperbolische Instabilität“. Stellen Sie sich einen Kreisel vor, der heftig wackelt, kurz bevor er umkippt.
3. Der Aufbruch: Dieses heftige Wackeln erzeugt sehr dünne, scharfe Schichten vertikaler Scherung. Diese dünnen Schichten brechen dann in kleine, heftige Wirbel auf, genau wie in Weg 1.
   Das Ergebnis: Die Vermischung findet statt, aber sie ist weniger effizient. Warum? Weil dieser Weg extrem dünne, scharfe Schichten erzeugt. Es braucht viel Energie (Reibung), um diese winzigen, scharfen Schichten zu erzeugen und zu brechen. Es ist, als würde man versuchen, einen dicken Block Käse mit einem stumpfen Messer zu schneiden (Weg 1) im Vergleich zu einer Rasierklinge (Weg 2); die Rasierklinge erzeugt einen viel schärferen, energieintensiveren Schnitt.

Das Wichtigste in Kürze

Die Arbeit beweist, dass vertikale Scherung (Auf-und-Ab-Bewegung) nicht von vornherein vorhanden sein muss, um zu starten. Sie ist ein unvermeidliches Nebenprodukt horizontaler Scherung in stark geschichteten Fluiden, vorausgesetzt, das Fluid ist dick genug (hohe Reynolds-Zahl).

• Wenn man mit zufälligem Rauschen beginnt: Erhält man Weg 1 (Direkter, effizienter Mischprozess).
• Wenn man mit einem spezifischen Muster beginnt: Erhält man Weg 2 (Indirekter, weniger effizienter Mischprozess).

Die Forscher nutzten leistungsstarke Computersimulationen, um zu zeigen, dass diese beiden Wege real, unterscheidbar sind und dass das „Rezept“, das man zu Beginn wählt, bestimmt, wie viel Energie als Wärme verschwendet wird im Vergleich dazu, wie viel für die eigentliche Vermischung der Schichten genutzt wird.

Kurz gesagt: Selbst in einem vollkommen ruhigen, geschichteten Fluid wird ein horizontaler Stoß schließlich vertikales Chaos erzeugen. Aber je nachdem, wie man den Stoß initiiert, wird dieses Chaos anders aussehen und die Schichten mit unterschiedlicher Effizienz vermischen.

Technische Zusammenfassung

Problemstellung
Die Stabilität großskaliger, stark geschichteter horizontaler Strömungen gegenüber dreidimensionalen (3D) Störungen ist ein entscheidender Mechanismus für den Downscale-Energietransfer und den diapyknen Transport in geophysikalischen und astrophysikalischen Systemen. Während vertikal gescherte geschichtete Strömungen gut untersucht sind, haben deren horizontal gescherte Gegenstücke (ohne initiale vertikale Scherung) weniger Aufmerksamkeit erhalten. Eine zentrale Frage bleibt: Wie gehen anfänglich zweidimensionale (2D), stark geschichtete ($Fr \ll 1$) horizontale Scherströmungen in Turbulenz über und erzeugen dabei vertikale Scherung? Die bisherige Literatur deutet auf zwei verschiedene Wege hin: einen, der durch das direkte Wachstum von 3D-Primärmoden zu Kelvin-Helmholtz (KH)-Instabilitäten führt, und einen anderen, der über die nichtlineare Entwicklung in Säulenwirbel (columnar vortices) gefolgt von Sekundärinstabilitäten verläuft. Es bleibt jedoch unklar, ob diese Pfade aufgrund unterschiedlicher Anfangsbedingungen (z. B. hyperbolischer Tangens gegenüber Sinusprofilen) oder Problemformulierungen (Anfangswertproblem vs. körpergetriebene Strömung) gegenseitig ausschließend sind oder ob sie alternative Routen darstellen, die innerhalb desselben physikalischen Systems verfügbar sind. Zudem wurde der Einfluss des gewählten Pfades auf die Effizienz der diapyknen Mischung bisher nicht systematisch quantifiziert.

Methodik
Die Autoren kombinieren lineare Stabilitätstheorie mit direkten numerischen Simulationen (DNS), um diese Übergänge in einem trippelperiodischen Gebiet zu untersuchen. Die Studie konzentriert sich auf eine vertikal invariante, horizontale, plan-parallele Kolmogorov-Strömung ($\bar{u}(y) = \sin(y)\mathbf{e}_x$) in einer linear geschichteten Flüssigkeit mit konstanter Auftriebsfrequenz $N$. Die Simulationen werden bei hohen Reynolds- ($Re = 10.000$) und Péclet-Zahlen ($Pe = 10.000$) durchgeführt, bei niedrigen Froude-Zahlen ($Fr = 0,15$ und $Fr = 0,1$), was zu Buoyancy-Reynolds-Zahlen ($Re_b = Fr^2 Re$) von 225 bzw. 100 führt.

Um den Effekt der Anfangsbedingungen auf die Pfadwahl zu isolieren, werden zwei Sätze von DNS durchgeführt:

1. DNS1 (Pfad P1): Initialisiert mit gitterskaligem weißem Rauschen, wodurch alle instabilen Moden der primären horizontalen Scherinstabilität wachsen können.
2. DNS2 (Pfad P2): Initialisiert mit gitterskaligem Rauschen plus einem spezifischen „Seed“ für den am schnellsten wachsenden 2D-Eigenmodus ($k_z=0$) der primären Instabilität.

Zusätzlich führen die Autoren eine lineare Stabilitätsanalyse eines „eingefrorenen“ analytischen Modells der in DNS2 entstehenden Säulenströmung durch. Dieses Modell approximiert die nichtlinear entwickelte 2D-Strömung als ein Taylor-Green-Wirbel-Array, das mit Sinusströmungen überlagert ist, wodurch die Berechnung von Eigenmoden und Wachstumsraten sekundärer 3D-Störungen ermöglicht wird.

Wesentliche Beiträge und Ergebnisse

1. Bestätigung zweier Pfade in einem einzigen Setup:
   Die Studie zeigt, dass beide Wege zur Turbulenz innerhalb desselben Modell-Setups auftreten können, wobei sie allein durch das initiale Störungspektrum unterschieden werden:

• Pfad P1 (Direktes 3D-Wachstum): Wenn mit weißem Rauschen initialisiert, wachsen $k_z \neq 0$-Eigenmoden der primären horizontalen Scherinstabilität direkt. Diese Moden erzeugen vertikale Scherung über einen breiten Skalenbereich, die anschließend bei ausreichend großem $Re_b$ instabil gegenüber sekundären, kleinskaligen KH-Instabilitäten wird. Dieser Pfad steht im Einklang mit Beobachtungen in körpergetriebenen Kolmogorov-Strömungen.
• Pfad P2 (Säulenartige Entwicklung): Wenn der $k_z=0$-Modus angeregt wird, dominiert dieser die Strömung und treibt den Hintergrund in eine langlebige, zeitabhängige 2D-Säulenströmung (vortical flow). Dieser 2D-Zustand ist anschließend instabil gegenüber sekundären 3D-Moden. Die Autoren identifizieren diese sekundären Moden als hyperbolische Instabilitäten, die aus den hyperbolischen Stagnationspunkten des Wirbel-Arrays resultieren. Diese Instabilitäten weisen hohe vertikale Wellenzahlen auf ($k_z \sim Fr^{-1}$) und treiben die Entstehung vertikaler Scherung voran, was schließlich tertiäre KH-Instabilitäten auslöst.

2. Charakterisierung der Sekundärinstabilität in Pfad P2:
   Durch die lineare Stabilitätsanalyse der eingefrorenen Säulenströmung zeigen die Autoren, dass die Sekundärinstabilität eine hochfrequente ($high-k_z$) hyperbolische Instabilität ist. Die Wachstumsrate dieser Instabilität skaliert mit der lokalen Dehnungsrate der hyperbolischen Punkte ($\Delta_h$) und ist näherungsweise $\lambda \approx 0,46 \Delta_h$, wobei sie für große Amplituden weitgehend unabhängig von der Schichtung ist. Die Instabilität manifestiert sich als „Rippel“ im vertikalen Geschwindigkeitsfeld, die um die hyperbolischen Punkte zentriert sind, und unterscheidet sich damit von der traditionellen „Zickzack“-Instabilität, die bei isolierten Wirbelpaaren beobachtet wird.

3. Einfluss auf die Mischungseffizienz:
   Die Studie quantifiziert die diapykne Mischung und die Mischungseffizienz ($\eta$) für beide Pfade.

• Buoyancy-Dissipation ($\chi$): Die Spitzenrate der Auftriebsvarianzen-Dissipation ist in beiden Pfaden vergleichbar, da sie in beiden Fällen durch kleinskalige Bewegungen ($k_z \gtrsim Fr^{-1}$) kontrolliert wird.
• Viskose Dissipation ($\epsilon$): Pfad P2 weist eine signifikant höhere viskose Dissipation auf als Pfad P1. Dies liegt daran, dass die hyperbolische Instabilität in P2 ein schmales, hochfrequentes ($high-k_z$) Spektrum anregt, was zu intensiver Dissipation in dünnen vertikalen Scherlagen führt. Im Gegensatz dazu regt P1 ein breites Spektrum vertikaler Skalen an, wodurch die Dissipation breiter verteilt wird.
• Mischungseffizienz: Folglich ist die maximale Mischungseffizienz in Pfad P1 ($\eta_{max} \approx 0,25$) signifikant höher als in Pfad P2 ($\eta_{max} \approx 0,165$). Der Unterschied ergibt sich daraus, dass P2 mehr kinetische Energie im Verhältnis zum gemischten Auftrieb dissipiert.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass die Entstehung vertikaler Scherung, die kleine KH-Instabilitäten antreibt, eine unvermeidliche Nebenwirkung von horizontalen Scherinstabilitäten in stark geschichteten Strömungen bei ausreichend großem $Re_b$ ist, sofern die Strömung einen Inflektionspunkt besitzt. Die Wahl zwischen den beiden Pfaden wird durch die Anfangsbedingungen bestimmt: speziell durch das Verhältnis der initialen Amplitude des $k_z=0$-Modus relativ zu den $k_z \neq 0$-Moden.

Die Autoren legen nahe, dass die Verbreitung von Pfad P2 in der bisherigen Literatur (speziell in Studien zu hyperbolischen Tangens-Scherlagen) ein Artefakt der Initialisierungspraktiken sein könnte, bei denen oft ein $k_z=0$-Seed hinzugefügt wird, um die Berechnungen zu beschleunigen, wodurch das direkte Wachstum von 3D-Moden (Pfad P1) unterdrückt wird. Umgekehrt kann das Fehlen von P1 in einigen Kontexten auf schwächere Schichtungslimits zurückzuführen sein, in denen der Bereich der instabilen $k_z \neq 0$-Moden begrenzt ist.

Letztlich etabliert die Arbeit, dass horizontale Scherinstabilitäten Schlüsselakteure bei der Erzeugung von geschichteter Turbulenz und diapykner Mischung sind, selbst wenn keine präexistente vertikale Scherung vorhanden ist. Sie hebt hervor, dass der spezifische Weg zur Turbulenz (direktes 3D-Wachstum vs. 2D-zu-3D-Übergang) das Energiedissipationsspektrum und damit die Effizienz der Mischung grundlegend verändert – ein Faktor, der bei der Parametrisierung von Turbulenz in Klima- und geophysikalischen Modellen berücksichtigt werden muss. Die Autoren merken an, dass während ihre Ergebnisse für idealisierte Strömungen robust sind, das Zusammenspiel mit Rotation und dem umgebenden Wellenfeld in realistischen geophysikalischen Settings eine offene Frage bleibt.