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

Deze studie combineert lineaire theorie en directe numerieke simulaties om te onthullen dat in sterk gestratificeerde stromingen zonder initiële verticale schering, horizontale scheringinstabiliteiten onvermijdelijk diapycne menging aansturen via twee afzonderlijke paden—ofwel via de directe opkomst van verticale schering of door een nietlineaire evolutie naar kolomvormige wervels—waarbij beide uiteindelijk kleine Kelvin-Helmholtz-instabiliteiten triggeren maar verschillende mengefficiënties opleveren vanwege hun excitatie van afzonderlijke verticale schalen.

De kern

Stel je de oceaan of de atmosfeer voor als een gigantische, gelaagde taart. De lagen zijn gemaakt van vloeistoffen met verschillende dichtheden (zoals verschillende smaken taart) en ze houden er niet van om gemakkelijk te mengen. Meestal bestuderen wetenschappers wat er gebeurt als je deze lagen tegelijkertijd zijwaarts én op en neer duwt (verticale afschuiving). Maar dit artikel stelt een andere vraag: wat gebeurt er als je de lagen alleen zijwaarts duwt, zonder een initiële op-en-neer beweging, in een zeer sterk gelaagde omgeving?

De onderzoekers ontdekten dat zelfs als je begint met een perfect horizontale stroming, de natuur twee verschillende "recepten" of paden heeft om chaos (turbulentie) en de lagen te mengen. Welk recept de natuur kiest, hangt volledig af van hoe je het experiment "zaait"—eigenlijk van de kleine, subtiele duw die je aan het begin aan de vloeistof geeft.

Hier is de uitsplitsing van de twee paden met eenvoudige analogieën:

De Opstelling: De Kalme Rivier

Stel je een brede, kalme rivier voor die horizontaal stroomt. Het water is gelaagd als een stapel pannenkoeken (sterke stratificatie). In het begin is de stroming glad en tweedimensionaal (het beweegt alleen naar links en rechts, niet op en neer).

Pad 1: Het "Overvolle Kamer" Effect (De Directe Route)

Hoe het begint: Je geeft de rivier een kleine, willekeurige duw overal tegelijkertijd (zoals het in de lucht gooien van een handvol confetti).
Wat er gebeurt:

1. De Rimpeling: Omdat van de lagen, begint de vloeistof niet alleen naar links en rechts te rimpelen; de vloeistof begint onmiddellijk ook op en neer te rimpelen in veel verschillende groottes tegelijkertijd. Denk hierbij aan een menigte mensen in een kamer die allemaal tegelijk proberen te bewegen, wat een chaotische, multidirectionele beweging creëert.
2. De Afschuiving: Deze rimpelingen creëren snel sterke verticale stromingen (afschuiving).
3. De Breuk: Deze verticale stromingen worden zo sterk dat ze knappen, wat kleine, gewelddadige wervelingen (eddies) creëert (zoals kleine draaikolken). Dit is de "Kelvin-Helmholtz" instabiliteit, die lijkt op de brekende golven die je ziet wanneer de wind over het water waait.
   Het Resultaat: Het mengen gebeurt efficiënt. Omdat de energie verspreid is over veel verschillende groottes van rimpelingen, is de "wrijving" (viskeuze dissipatie) lager, wat het mengproces relatief efficiënt maakt.

Pad 2: Het "Gesynchroniseerde Dans" Effect (De Indirecte Route)

Hoe het begint: Je geeft de rivier een zeer specifieke, georganiseerde duw (zoals een dirigent die met een baton zwaait om iedereen in een specifiek patroon te laten bewegen).
Wat er gebeurt:

1. De Vortex: In plaats van chaotische rimpelingen, organiseert de vloeistof zich in lange, verticale kolommen van kolkend water (zoals gigantische tornado's die in de rivier staan). Voor een lange tijd blijft de stroming perfect tweedimensionaal, bestaande uit slechts deze grote draaikolommen.
2. Het Waggelen: Uiteindelijk worden deze gigantische kolommen instabiel. Ze beginnen op een zeer specifieke, hoogfrequente manier te wankelen. De onderzoekers noemen dit een "hyperbolische instabiliteit". Stel je een tol voor die heftig begint te wankelen vlak voordat hij omvalt.
3. De Breuk: Deze heftige wankeling creëert zeer dunne, scherpe lagen van verticale afschuiving. Deze dunne lagen breken vervolgens af in kleine, gewelddadige wervelingen, net als in Pad 1.
   Het Resultaat: Het mengen vindt plaats, maar het is minder efficiënt. Waarom? Omdat dit pad extreem dunne, scherpe lagen creëert. Het kost veel energie (wrijving) om deze kleine, scherpe lagen te creëren en te breken. Het is als het proberen te snijden van een dik blok kaas met een bot mes (Pad 1) versus een scheermesje (Pad 2); het scheermesje creëert een veel scherpere, meer energie-intensieve snede.

De Belangrijkste Conclusie

Het artikel bewijst dat verticale afschuiving (op-en-neer beweging) niet aanwezig hoeft te zijn om te beginnen. Het is een onvermijdelijk bijproduct van horizontale afschuiving in sterk gelaagde vloeistoffen, mits de vloeistof dik genoeg is (hoog Reynoldsgetal).

• Als je begint met willekeurige ruis: Krijg je Pad 1 (Direct, efficiënt mengen).
• Als je begint met een specifiek patroon: Krijg je Pad 2 (Indirect, minder efficiënt mengen).

De onderzoekers gebruikten krachtige computersimulaties om aan te tonen dat deze twee paden echt en onderscheidend zijn, en dat de "receptuur" die je aan het begin kiest, bepaalt hoeveel energie verloren gaat als warmte versus hoeveel er wordt gebruikt om de lagen daadwerkelijk te mengen.

Kortom: Zelfs in een perfect kalme, gelaagde vloeistof zal een horizontale duw uiteindelijk leiden tot verticale chaos. Maar afhankelijk van hoe je de duw start, zal die chaos er anders uitzien en de lagen met verschillende niveaus van efficiëntie mengen.

Technische samenvatting

Probleemstelling
De stabiliteit van grootschalige, sterk gestratificeerde horizontale stromingen ten opzichte van driedimensionale (3D) perturbaties is een cruciaal mechanisme voor de neerwaartse energietransport en de diapycnal transport in geofysische en astrofysische systemen. Hoewel verticaal afgeschoven gestratificeerde stromingen goed bestudeerd zijn, hebben hun horizontaal afgeschoven tegenhangers (zonder initiële verticale afschering) minder aandacht gekregen. Een centrale vraag blijft: hoe gaan aanvankelijk tweedimensionale (2D), sterk gestratificeerde ($Fr \ll 1$) horizontale afscheringsstromingen over naar turbulentie en genereren zij verticale afschering? De bestaande literatuur suggereert twee verschillende paden: één waarbij de directe groei van 3D primaire modi leidt tot Kelvin-Helmholtz (KH) instabiliteiten, en een andere waarbij de nietlineaire evolutie verloopt naar kolomvormige wervels gevolgd door secundaire instabiliteiten. Het is echter onduidelijk of deze paden wederzijds uitsluitend zijn vanwege verschillen in beginvoorwaarden (bijv. hyperbooltangens versus sinusoïdaal profiel) of probleemformulering (Initial Value Problem versus body-forced), of dat ze alternatieve routes vertegenwoordigen die beschikbaar zijn binnen hetzelfde fysieke systeem. Bovendien is de impact van het gekozen pad op de efficiëntie van diapycnal menging nog niet systematisch gekwantificeerd.

Methodologie
De auteurs combineren lineaire stabiliteitstheorie met Direct Numerical Simulations (DNS) om deze transities te onderzoeken in een trippel-periodiek domein. De studie richt zich op een verticaal invariante, horizontale, plane-parallel Kolmogorov-stroom ($\bar{u}(y) = \sin(y)\mathbf{e}_x$) in een lineair gestratificeerde vloeistof met een constante buoyancy frequentie $N$. De simulaties worden uitgevoerd bij hoge Reynolds-getallen ($Re = 10.000$) en Péclet-getallen ($Pe = 10.000$), met lage Froude-getallen ($Fr = 0,15$ en $Fr = 0,1$), wat resulteert in buoyancy Reynolds-getallen ($Re_b = Fr^2 Re$) van respectievelijk 225 en 100.

Om het effect van de beginvoorwaarden op de selectie van het pad te isoleren, worden twee sets DNS uitgevoerd:

1. DNS1 (Pad P1): Geïnitialiseerd met grid-schaal witte ruis, waardoor alle onstabiele modi van de primaire horizontale afscharing kunnen groeien.
2. DNS2 (Pad P2): Geïnitialiseerd met grid-schaal ruis plus een specifieke "seed" voor de snelst groeiende 2D ($k_z=0$) eigenmode van de primaire instabiliteit.

Daarnaast voeren de auteurs een lineaire stabiliteitsanalyse uit van een "frozen-in" analytisch model van de kolomvormige stroming die in DNS2 ontstaat. Dit model benadert de nietlineair geëvolueerde 2D-stroming als een Taylor-Green wervelarray superponeerd met sinusoïdale stromingen, waardoor de berekening van eigenmodes en groeisnelheden van secundaire 3D-perturbaties mogelijk is.

Belangrijkste bijdragen en resultaten

1. Bevestiging van twee paden in een enkele opstelling:
   De studie toont aan dat beide paden naar turbulentie kunnen optreden binnen exact dezelfde modelopstelling, onderscheiden door enkel het initiële perturbatiespectrum:

• Pad P1 (Directe 3D Groei): Wanneer geïnitialiseerd met witte ruis, groeien de $k_z \neq 0$ eigenmodes van de primaire horizontale afscharing direct. Deze modi genereren verticale afschering op een breed scala aan schalen, die vervolgens onstabiel wordt voor secundaire, kleinschalige KH-instabiliteiten wanneer $Re_b$ voldoende groot is. Dit pad komt overeen met observaties in body-forced Kolmogorov-stromen.
• Pad P2 (Kolomvormige Evolutie): Wanneer de $k_z=0$ mode wordt gezaaid, domineert deze de stroming en drijft de achtergrond aan naar een langdurige, tijdsafhankelijke 2D kolomvormige (wervel-)stroming. Deze 2D-toestand is vervolgens onstabiel voor secundaire 3D-modi. De auteurs identificeren deze secundaire modi als hyperbolische instabiliteiten voortkomend uit de hyperbolische stagnatiepunten in de wervelarray. Deze instabiliteiten hebben hoge verticale golfgetallen ($k_z \sim Fr^{-1}$) en drijven de opkomst van verticale afschering aan, wat uiteindelijk tertiaire KH-instabiliteiten triggert.

2. Karakterisering van de secundaire instabiliteit in Pad P2:
   Door middel van lineaire stabiliteitsanalyse van de "frozen" kolomvormige stroming tonen de auteurs aan dat de secundaire instabiliteit een hoog-$k_z$ hyperbolische instabiliteit is. De groeisnelheid van deze instabiliteit schaalt met de lokale afscheringssnelheid van de hyperbolische punten ($\Delta_h$) en is ongeveer $\lambda \approx 0,46 \Delta_h$, grotendeels onafhankelijk van de sterkte van de stratificatie voor grote amplitudes. De instabiliteit manifesteert zich als "rimpels" in het verticale snelheidsveld gecentreerd rond hyperbolische punten, wat verschilt van de traditionele "zigzag"-instabiliteit waargenomen in geïsoleerde wervelparen.

3. Impact op mengefficiëntie:
   De studie kwantificeert de diapycnal menging en de mengefficiëntie ($\eta$) voor beide paden.

• Buoyancy Dissipatie ($\chi$): De piek van de rate van de dissipatie van de buoyancy variantie is vergelijkbaar in beide paden, aangezien deze wordt gecontroleerd door kleinschalige bewegingen ($k_z \gtrsim Fr^{-1}$) in beide gevallen.
• Viskeuze Dissipatie ($\epsilon$): Pad P2 vertoont aanzienlijk hogere viskeuze dissipatie dan Pad P1. Dit komt doordat de hyperbolische instabiliteit in P2 een smal, hoog-$k_z$ spectrum exciteert, wat leidt tot intense dissipatie in dunne verticale afscheringslagen. In contrast hiermee exciteert P1 een breed spectrum van verticale schalen, waardoor de dissipatie breder wordt verdeeld.
• Mengefficiëntie: Bijgevolg is de piek mengefficiëntie aanzienlijk hoger in Pad P1 ($\eta_{max} \approx 0,25$) dan in Pad P2 ($\eta_{max} \approx 0,165$). Dit verschil komt doordat P2 meer kinetische energie dissipeert relatief aan de hoeveelheid gemengde buoyancy.

Betekenis en Claims
Het artikel claimt dat de opkomst van verticale afschering die kleine schaal KH-instabiliteiten aandrijft, een onvermijdelijk bijproduct is van horizontale afscharing-instabiliteiten in sterk gestratificeerde stromingen bij een voldoende grote $Re_b$, mits de stroming een inflexiepunt heeft. De selectie tussen de twee paden wordt bepaald door de beginvoorwaarden: specifiek de ratio van de initiële amplitude van de $k_z=0$ mode ten opzichte van de $k_z \neq 0$ modes.

De auteurs suggereren dat de prevalentie van Pad P2 in eerdere literatuur (specifiek in studies naar hyperbooltangens afscheringslagen) mogelijk een artefact is van de initialisatiepraktijken waarbij een $k_z=0$ seed vaak wordt toegevoegd om berekeningen te versnellen, waardoor de directe groei van 3D-modi (Pad P1) wordt onderdrukt. Omgekeerd kan het gebrek aan P1 in sommige contexten te wijten zijn aan zwakkere stratificatie-limieten waarbij het bereik van onstabiele $k_z \neq 0$ modes beperkt is.

Uiteindelijk stelt dit werk vast dat horizontale afscharing-instabiliteiten sleutelspelers zijn bij het genereren van gestratificeerde turbulentie en diapycnal menging, zelfs in de afwezigheid van vooraf bestaande verticale afschering. Het benadrukt dat de specifieke route naar turbulentie (directe 3D-groei versus 2D-naar-3D transitie) het energie-dissipatiespectrum en daarmee de efficiëntie van menging fundamenteel verandert, een factor die moet worden meegewogen bij het parametriseren van turbulentie in klimaat- en geofysische modellen. De auteurs merken op dat hoewel hun resultaten robuust zijn voor geïdealiseerde stromingen, de interactie met rotatie en ambient wave velden in realistische geofysische settings een openstaande vraag blijft.