Early onset of secondary shear instability in Kelvin-Helmholtz braids at high Reynolds number

Dit artikel toont aan dat bij hoge Reynolds- en Richardson-getallen de secundaire schuifinstabiliteit in de braid-regio's van Kelvin-Helmholtz-billows eerder optreedt dan de verzadiging van de primaire billows, waardoor deze instabiliteit de overgang naar turbulentie en diapycnaal mengen kan beheersen.

De kern

De Stille Moordenaar in de Golf: Waarom de 'Vlecht' van een Golf sneller breekt dan de Golf zelf

Stel je voor dat je naar de oceaan kijkt op een dag met veel wind. Je ziet enorme, rollende golven. In de natuurkunde noemen we deze rollende structuren Kelvin-Helmholtz-billows (of simpelweg: de grote draaikolken).

Meestal denken we dat deze grote rollen eerst volledig moeten uitgroeien, tot ze hun maximale grootte bereiken, en dan pas uit elkaar vallen tot kleine, chaotische turbulentie. Maar dit nieuwe onderzoek laat zien dat er een geheimzinnige, snellere kracht werkt die dit proces verstoort.

1. De Grote Rol en de Dunne Vlecht

Stel je de grote draaikolken voor als twee enorme, draaiende bollen boter die tegen elkaar aan rollen. Tussen deze twee bollen zit een heel dunne, schuine strook boter. In de wetenschap noemen we dit de vlecht (of braid).

• De oude theorie: Wetenschappers dachten dat deze vlecht stabiel blijft totdat de grote bollen boter volledig zijn uitgegroeid en "vol" zijn. Pas dan zou de vlecht dun genoeg worden om te breken.
• De nieuwe ontdekking: De onderzoekers (Bouckley, Lewin en Lefauve) hebben ontdekt dat de vlecht soms veel eerder breekt, nog voordat de grote bollen boter klaar zijn.

2. De Vergelijking: Een Strakke Lijn vs. Een Knijpende Hand

Hoe kan dat? Het heeft te maken met twee tegenstrijdige krachten die op de vlecht werken:

1. De Knijpende Hand (Straling): De grote draaikolken die eromheen draaien, trekken de vlecht uit. Het is alsof je een stuk deeg uitrekt. Dit maakt de vlecht dunner en strakker. Normaal gesproken maakt dit de vlecht stabieler (zoals een strakke snaar die niet snel breekt).
2. De Verwarmde Boter (Baroclinische Kracht): Omdat de vlecht schuin staat in het water (waar het water van verschillende dichtheden is), werkt er een andere kracht. De zwaartekracht en de dichtheidsverschillen zorgen ervoor dat de vlecht van binnenuit "opwarmt" en instabiel wordt.

Het geheim: Bij een hoge dichtheid (stratificatie) werkt de "knijpende hand" langzamer (de grote bollen groeien trager), maar de "verwarmde boter" werkt juist sneller. De vlecht wordt zo dun en zo instabiel dat hij breekt, terwijl de grote bollen er nog rustig aan doen om te groeien.

3. De Analogie: De Slingerende Lijn

Stel je een slingerende lijn voor tussen twee mensen die aan het touwtrekken zijn.

• Als de mensen langzaam trekken, blijft de lijn strak en stabiel.
• Maar als de mensen plotseling heel hard in een specifieke richting trekken (de "verwarmde boter"), kan de lijn in het midden knappen, zelfs als de mensen zelf nog niet hun maximale kracht hebben bereikt.

In dit onderzoek laten ze zien dat bij een bepaalde snelheid en dichtheid (in de oceaan: hoge Reynolds- en Richardson-getallen), deze "knijpkracht" in het midden van de lijn zo sterk wordt dat de lijn voordat de trekkracht volledig is opgebouwd, al uit elkaar valt.

4. Waarom is dit belangrijk voor de oceaan?

Dit klinkt misschien als pure theorie, maar het heeft grote gevolgen voor hoe we de wereld begrijpen:

• Menging van water: In de oceaan moet zout water en zoet water, of warm en koud water, met elkaar mengen. Dit is cruciaal voor het klimaat en het leven in de zee.
• De verrassing: Vroeger dachten we dat de grote draaikolken de belangrijkste plek waren waar dit menging plaatsvond. Dit onderzoek toont aan dat de vlecht (het dunne stukje ertussen) vaak de eerste en belangrijkste plek is waar het water volledig wordt doorelkaar geschud.
• Veldwaarnemingen: Wetenschappers hebben in de echte oceaan al gemerkt dat het water vaak het meest turbulent is in die dunne strepen tussen de golven, en niet in de golven zelf. Dit artikel geeft de wiskundige verklaring waarom dat zo is.

5. De Simulatie: Een Digitale Oceaan

De onderzoekers hebben supercomputers gebruikt om deze situatie na te bootsen. Ze hebben gekeken naar waterstromen met een snelheid en dichtheid die voorkomen in de echte oceaan (van de kust tot diep in de zee).

Hun conclusie is duidelijk:

• Als het water "dik" genoeg is (hoge stratificatie), breekt de vlecht vroegtijdig.
• Dit gebeurt zelfs voordat de grote draaikolken hun maximale grootte bereiken.
• Dit betekent dat de "chaos" (turbulentie) eerder begint dan we dachten, en dat de vlecht de regisseur is van hoe het water zich mengt.

Samenvatting in één zin

Dit onderzoek laat zien dat in de oceaan de dunne, schuine strepen tussen de grote draaikolken vaak eerder uit elkaar vallen en turbulent worden dan de grote kolken zelf, omdat de krachten die ze dunner maken, sneller werken dan de krachten die ze stabiliseren.

Het is alsof de "naad" van een jas eerder scheurt dan dat de jas zelf volledig is opgemaakt.

Technische samenvatting

Probleemstelling

De overgang van laminair naar turbulent gedrag in stabiel gestratificeerde, parallelle schuifstromingen via Kelvin-Helmholtz (KH) instabiliteit is een langdurig open probleem in de fluïdamechanica. Traditioneel wordt aangenomen dat deze overgang plaatsvindt nadat primaire "billows" (wervelstructuren) tot een eindige amplitude zijn gegroeid en vervolgens instabiel worden. Echter, veldwaarnemingen in de oceanen tonen vaak aan dat intensieve turbulente menging geconcentreerd is in de "vlechtstukken" (braids) die de billows verbinden, en niet in de billows zelf.

De vraag is wanneer en waarom secundaire instabiliteiten in deze vlechtstukken optreden. Er zijn drie klassen van secundaire instabiliteiten:

1. Wervelpaaringsinstabiliteiten (vortex pairing).
2. Secundaire convectieve instabiliteiten in de kern van de billow.
3. Secundaire schuifinstabiliteit (SSI) in de vlechtstukken.

Eerdere studies (zoals Corcos & Sherman, 1976; Staquet, 1995) namen vaak aan dat de achtergrondstroom "bevroren" is na verzadiging van de primaire billow. Dit artikel betoogt dat deze aanname onjuist is voor geofysisch relevante Reynolds- en Richardson-getallen, waarbij SSI veel eerder kan optreden dan de verzadiging van de billow, gedreven door de tijdsafhankelijke verdikking en versterking van de schuif in het vlechtstuk.

Methodologie

De auteurs combineren een semi-analytisch model met hoogwaardige numerieke simulaties:

1. Theoretisch Model (Onviskeus):

   • De auteurs passen de klassieke analyse van Corcos & Sherman aan door een coördinatenstelsel te gebruiken dat is uitgelijnd met het vlechtstuk (geïclineerd onder een hoek $\phi$).
   • Ze modelleren de tijdsafhankelijke evolutie van het vlechtstuk door een uniform compressief rekveld ($\gamma$) en barocliene schuifversterking te beschouwen.
   • Ze leiden een stabiliteitscriterium af dat niet alleen gebaseerd is op de Richardson-getal ($Ri_g < 1/4$), maar vooral op de verhouding tussen rek en schuif: $\gamma/\Omega < \Gamma^*$.
   • Het model voorspelt dat bij hoge stratificatie (hoge $Ri$) de billow-groei vertraagt, maar de barocliene schuifproductie in het vlechtstuk versnelt, waardoor SSI eerder kan ontstaan.

2. Numerieke Simulaties (DNS):

   • Er zijn 16 directe numerieke simulaties (DNS) uitgevoerd met de solver Oceananigans.jl op GPU's.
   • Het parameterbereik omvat Richardson-getallen $Ri \in [0.05, 0.20]$ en Reynolds-getallen $Re \in [10^4, 10^7]$.
   • De simulaties zijn beperkt tot twee dimensies en één primaire golflengte om vortex-pairing en convectieve instabiliteiten uit te sluiten en zich puur op SSI te focussen.
   • De simulaties worden gestopt bij de verzadiging van de primaire billow ($t_{2D}$) of bij het begin van SSI ($t_S$), afhankelijk van welke eerst optreedt.

Belangrijkste Bijdragen

• Tijdsafhankelijk Model: De auteurs ontwikkelen een onviskeus, tijdsafhankelijk model voor de evolutie van het vlechtstuk, in plaats van de gebruikelijke aanname van een statische achtergrondstroom na billow-saturatie.
• Nieuw Stabiliteitscriterium: Ze introduceren en valideren een criterium voor SSI-ontstaan gebaseerd op de verhouding $\gamma/\Omega$ (rek/schuif), met een kritieke drempelwaarde $\Gamma^* \approx 0.02$.
• Mechanisme voor Vroege Ontwikkeling: Het paper biedt een mechanistische verklaring voor veldobservaties waarbij menging in vlechtstukken domineert, door aan te tonen dat SSI de primaire billow-verzadiging kan voorafgaan bij hoge $Ri$ en $Re$.

Resultaten

• Vroege Ontwikkeling van SSI: Bij voldoende hoge Richardson-getallen ($Ri \gtrsim 0.10$) en hoge Reynolds-getallen ($Re \ge 10^5$) ontwikkelt SSI zich significant vóór de verzadiging van de primaire billow ($t_S < t_{2D}$).
• Rol van Reynolds-getal:
  • Bij zeer hoge $Re$($10^6 - 10^7$) gedraagt het vlechtstuk zich als onviskeus; de simulaties komen perfect overeen met het theoretische model.
  • Bij lagere $Re$($10^4$) vertraagt viscositeit het dunner worden van het vlechtstuk, wat de barocliene schuifversterking remt en het ontstaan van SSI kan voorkomen of uitstellen, zelfs bij $Ri$ waar het model SSI voorspelt.
• Rol van Richardson-getal: Een hogere stratificatie ($Ri$) vertraagt de groei van de billow (waardoor $t_{2D}$ toeneemt), maar versnelt de schuifversterking in het vlechtstuk. Dit creëert een groer tijdsvenster ($t_{2D} - t^*$) waarin SSI kan groeien tot een eindige amplitude.
• Kritieke Drempel: De simulaties bevestigen dat SSI optreedt wanneer de verhouding $\gamma/\Omega$ de kritieke waarde $\Gamma^* \approx 0.02$ bereikt. Zonder extra ruis (in marginale gevallen) is deze drempel noodzakelijk om de instabiliteit te activeren.
• Vergelijking met andere instabiliteiten: Omdat SSI bij hoge $Ri$ en $Re$ eerder optreedt dan wervelpairing (dat sterk afhankelijk is van $Ri$ en langzamer groeit) en dan convectieve instabiliteiten (die pas na saturatie significant worden), kan SSI de overgang naar 3D-turbulentie en de daaropvolgende diapycnale menging controleren.

Betekenis en Conclusie

Dit onderzoek biedt een fundamentele verklaring voor waarnemingen in de oceanen waarbij menging voornamelijk plaatsvindt in de vlechtstukken van KH-instabiliteiten, zelfs voordat de billows volledig zijn gevormd.

• Geofysische Relevantie: De bevindingen suggereren dat bij geofysisch relevante omstandigheden ($Re \sim 10^5 - 10^8$, $Ri \sim 0.1 - 0.2$), SSI de dominante mechanisme is voor de initiële overgang naar turbulentie. Dit heeft grote gevolgen voor het modelleren van menging in oceanen en atmosfeer.
• Numerieke Implicaties: De studie benadrukt dat bij het simuleren van KH-instabiliteiten op hoge Reynolds-getallen, de oplosbaarheid van de vlechtstukken cruciaal is. Als de interface niet fijn genoeg is opgelost, kunnen spuriële instabiliteiten ontstaan of echte SSI worden gemist. Het paper biedt een criterium voor de minimale grid-resolutie nodig om SSI correct te vangen.
• Toekomstperspectief: Hoewel de huidige studie 2D is, suggereert het dat 3D-simulaties nodig zijn om de volledige turbulentie en menging te begrijpen die door deze vroege SSI wordt gegenereerd. Ook de invloed van een hoge Prandtl-getal (typisch voor oceanische zout- en temperatuurstratificatie) blijft een open vraag voor toekomstig onderzoek.

Kortom, het paper toont aan dat de dynamiek van het vlechtstuk niet statisch is, maar een actief, tijdsafhankelijk proces is dat bij hoge stratificatie en Reynolds-getallen de overgang naar turbulentie kan initiëren voordat de primaire structuur verzadigt.