Instability and breaking of internal waves in a horizontal shear layer

Dit onderzoek combineert straalvolgingstheorie en directe numerieke simulaties om de instabiliteit en breking van interne golven in een horizontale schuiflaag te analyseren, waarbij een dimensieloze parameter $F$ wordt geïntroduceerd om te voorspellen of breking wordt veroorzaakt door golven die steil worden door refractie of door advektie van impuls, wat leidt tot significante turbulente dissipatie en menging.

De kern

Stel je voor dat de oceaan of de atmosfeer niet stil is, maar vol zit met onzichtbare golven. Dit zijn geen golven op het wateroppervlak die je kunt zien, maar interne golven die zich verplaatsen binnenin het water of de lucht, vaak tussen lagen van verschillende temperaturen of zoutgehalte. Deze golven dragen enorme hoeveelheden energie.

De vraag die deze onderzoekers zich stellen, is: Wat gebeurt er als deze interne golven in botsing komen met een stroming die langs de kant stroomt?

Hier is een uitleg van hun onderzoek, vertaald naar alledaagse taal en met een paar creatieve vergelijkingen.

1. Het Probleem: De Golf en de Stroom

Stel je een rivier voor. In het midden stroomt het water snel, aan de randen langzaam. Dit noemen we een schuiflaag (shear layer). Nu gooi je een grote, onzichtbare golf in deze rivier.

De onderzoekers wilden weten:

• Wordt de golf er sterker van door de stroming?
• Breekt de golf?
• En als hij breekt, hoe ziet dat eruit en hoeveel energie wordt er vrijgemaakt?

2. De Twee Manieren waarop een Golf "Breekt"

Het onderzoek toont aan dat er twee heel verschillende manieren zijn waarop deze golven kunnen instorten (breken), afhankelijk van hoe ze de stroming benaderen.

Manier A: De "Schaal" (Refraction)
Stel je voor dat je een auto rijdt over een weg die steeds steiler wordt. Je moet harder trappen om dezelfde snelheid te houden, en je banden worden steeds smaller en strakker.

• Wat er gebeurt: De zijwaartse stroming duwt de golf zo hard dat hij "op zijn kop" wordt geduwd. De golf wordt steeds steiler, totdat hij als een brekende oceaan golf over de top slaat.
• Het resultaat: De golf stort in door zwaartekracht. De lagen water draaien om elkaar heen (zoals een omgekeerde saladekom). Dit is een "convectieve" instabiliteit.

Manier B: De "Trekkracht" (Advection)
Stel je nu voor dat je op een skateboard staat en iemand duwt je van achteren, maar niet recht, maar schuin. Je snelheid neemt toe, maar je balans raakt verstoord op een heel andere manier.

• Wat er gebeurt: De golf "pakt" de stroming mee. Door de beweging van de golf wordt de achtergrondstroom zelf verstoord. Dit creëert een enorme verticale schuifkracht (alsof je twee lagen karton over elkaar schuift).
• Het resultaat: De golf breekt niet door omvallen, maar door wrijving en schuiven. Het lijkt meer op het ontstaan van turbulentie in een snelstromende rivier.

3. De "Voorspeller": De F-Waarde

De onderzoekers hebben een simpele formule bedacht (de F-waarde) om te voorspellen welke van de twee manieren gaat gebeuren.

• Kleine F: De golf wordt steiler en breekt door omvallen (Manier A).
• Grote F: De golf wordt verstoord door de stroming en breekt door schuiven (Manier B).

Het mooie is: hun simpele wiskundige theorie (die eigenlijk alleen voor kleine golven bedoeld was) bleek verrassend goed te werken, zelfs voor enorme, chaotische golven in hun computersimulaties. Het was alsof ze met een simpele weersvoorspelling de vorming van een orkaan konden voorspellen.

4. Het Grote Gebeuren: Turbulentie en Menging

Wanneer de golf breekt, gebeurt er iets fascinerends: Turbulentie.

• Energie-explosie: In sommige gevallen wordt er zoveel energie vrijgemaakt dat het vele malen meer is dan de energie die de golf zelf had. Het is alsof een kleine vonk een enorme brand veroorzaakt door de brandstof van de stroming zelf aan te wakkeren.
• Menging: Dit is cruciaal voor de aarde. De oceaan moet warmte en zuurstof van boven naar beneden mengen. Als golven breken, mengen ze het water.
  • Als de golf "omvalt" (Manier A), is de menging heel efficiënt.
  • Als de golf "schuift" (Manier B), is de menging minder efficiënt, maar wordt er wel meer momentum (beweging) uit de stroming gehaald, waardoor de stroming zelf vertraagt.

5. Waarom is dit belangrijk?

Dit klinkt als pure natuurkunde, maar het heeft enorme gevolgen voor ons klimaat.

• Klimaatmodellen: Om te voorspellen hoe de oceaan warmte opslaat of hoe CO2 wordt opgenomen, moeten we weten hoe goed het water gemengd wordt.
• Onzekerheid: Tot nu toe dachten wetenschappers dat deze menging altijd ongeveer hetzelfde was. Dit onderzoek toont aan dat het enorm kan variëren, afhankelijk van hoe de golf precies breekt. Soms is de menging 5 keer zo sterk als anders!

Samenvatting in één zin

De onderzoekers hebben ontdekt dat interne golven in de oceaan op twee heel verschillende manieren kunnen breken (door omvallen of door schuiven), en dat dit bepaalt hoeveel energie er vrijkomt en hoe goed het water wordt gemengd – een proces dat cruciaal is voor het klimaat van onze planeet.

De kernboodschap: Het is niet alleen belangrijk dat golven breken, maar vooral hoe ze breken. En dat "hoe" hangt af van de interactie met de stroming eromheen.

Technische samenvatting

Titel: Instabiliteit en breken van interne golven in een horizontale schuiflaag

Auteurs: Sam F. Lewin, Alexis K. Kaminski, Arun Balakrishna en Miles M.P. Couchman.
Doelwit: Journal of Fluid Mechanics (in behandeling).

---

1. Probleemstelling en Context

Interne zwaartekrachtgolven spelen een cruciale rol in het transport van energie en menging in de oceanen en de atmosfeer. Het mechanisme waarmee deze golven breken en hun energie dissiperen via kleine schaal-turbulentie is een fundamenteel probleem. Hoewel de interactie tussen interne golven en verticale schuifstromen (vertical shear) goed bestudeerd is (bijv. via kritieke niveaus waar de golfsnelheid gelijk is aan de stroomsnelheid), is de interactie met horizontale schuifstromen (waar de schuifrichting loodrecht staat op de zwaartekracht) minder onderzocht.

Deze studie onderzoekt de evolutie en instabiliteit van lineaire monochromatische interne golven die zich voortplanten in een achtergrondstroom met een horizontale schuif, beschreven door:
$$U^*(y^*) = \Delta u^* \tanh(y^*/h^*) \mathbf{e}_x$$
De focus ligt op hoe twee specifieke mechanismen voor energiegroei van de golf leiden tot breken en turbulentie:

1. Refractie: Verandering in golfsteilheid door de schuif.
2. Advectie van impuls: Verandering in stroomrichtings-velocity perturbaties door de golf die de achtergrondstroom meeneemt (vergelijkbaar met het "lift-up" mechanisme).

2. Methodologie

De auteurs combineren twee benaderingen:

A. Ray-tracing theorie (WKBJ-benadering)
Er wordt een lineaire analyse uitgevoerd onder de Boussinesq-approximatie en de veronderstelling van zwakke rotatie ($f \ll N$).

• Stralingsvergelijkingen: De golfpaden worden berekend langs stralen die tangentieel zijn aan de groepssnelheid.
• Nieuwe parameter $F$: Een dimensieloze verhouding van verstoring-energie wordt geconstrueerd om de relatieve belangrijkheid van de twee energiegroeimechanismen te kwantificeren:
  $$F^2 = \frac{k_y^2 k_z^2}{\Omega^4 |\mathbf{k}|^4} \left(\frac{\partial U}{\partial y}\right)^2$$
  • Een kleine $F$ impliceert dat energiegroei voornamelijk door refractie (toename van golfsteilheid $s$) komt, wat leidt tot convectieve instabiliteit (omkering van isopycnalen).
  • Een grote $F$ impliceert dat energiegroei voornamelijk door advectie van impuls komt, wat leidt tot een toename van de verticale schuif en mogelijk zuivere schuifinstabiliteit (Kelvin-Helmholtz).
• Voorspelling: De theorie voorspelt de lokale verdeling van kinetische en potentiële energie en de lokale Richardson-getal ($Ri_g$) om het type instabiliteit te voorspellen.

B. Directe Numerieke Simulaties (DNS)
Om de theorie te testen, worden volledig niet-lineaire DNS-simulaties uitgevoerd met de solver Diablo.

• Opstelling: Een golfpakket (Gaussian envelop) wordt gelanceerd in de horizontale schuiflaag.
• Parameterbereik: De simulaties variëren in initiële golfsteilheid ($s_0$), schuifsterkte ($\Delta u$), en golfvectoren ($k_x, k_y, k_z$) om scenario's te dekken van golfopsluiting (trapping, $k_x < 0$) tot golfreflectie ($k_x > 0$) en doorgang ($k_x = 0$).
• Doel: Het observeren van de overgang naar turbulentie, de energiebudgetten en de mengingsefficiëntie.

3. Belangrijkste Resultaten

A. Validatie van de Ray-tracing theorie
Hoewel de WKBJ-benadering strikt genomen alleen geldig is voor kleine amplitude-golven, toont de vergelijking met DNS aan dat de theorie kwalitatief uitstekende voorspellingen doet over het type instabiliteit:

• Golfopsluiting ($k_x < 0$): De golf wordt gevangen in een kritiek niveau. De steilheid $s$ neemt sterk toe.
  • Bij lage $F$: Instabiliteit wordt gedreven door convectie (overturning van dichtheidslagen).
  • Bij hoge $F$: De instabiliteit wordt gedreven door versterkte verticale schuif, zelfs als de golfsteilheid niet extreem groot wordt.
• Golfreflectie/Doorgang ($k_x \geq 0$): De golf wordt niet gevangen.
  • Bij hoge $F$: De advectie van impuls creëert sterke verticale schuif die leidt tot Kelvin-Helmholtz-billen en turbulentie, zonder dat de golf eerst "breekt" door convectie.

B. Energie- en Impulsoverdracht

• Dissipatie: In veel gevallen overschrijdt de totale turbulente dissipatie de initiële golfenergie met een factor 5 tot 7,5. Dit suggereert dat de achtergrondstroom (die lineair stabiel is) een significante energiebron is voor de turbulentie.
• Momentum: De interactie leidt tot de vorming van smalle stralen (jets) in de stroomrichting en een verandering in de dikte van de schuiflaag.
• Mixing Efficiency ($\eta$): De mengingsefficiëntie (de verhouding tussen potentiële energie die wordt gemengd en totale dissipatie) varieert sterk afhankelijk van het brekingsmechanisme:
  • Convectief breken (kleine $F$): $\eta \approx 0.42$ (zeer hoog, veel menging).
  • Schuif-gedreven breken (grote $F$): $\eta \approx 0.08$ (laag, weinig menging).
    Dit staat in contrast met de vaak aangenomen constante waarde van $\approx 0.2$ in oceanografische modellen.

C. Rol van Interne Golf-emissie
Tijdens het breken worden in sommige scenario's nieuwe interne golven gegenereerd. De energie van deze uitgezonden golven kan oplopen tot 25% van de totale gedissipeerde energie, wat aanzienlijk is en vaak wordt verwaarloosd in klassieke verticale schuif-modellen.

4. Bijdragen en Significantie

1. Unificatie van Mechanismen: De studie identificeert en kwantificeert twee distincte routes naar golfbreken in horizontale schuifstromen: convectief (door steilheid) en schuif-gedreven (door impulsadvectie). De parameter $F$ fungeert als een robuuste voorspeller voor welke route dominant is.
2. Voorspellende Kracht van Lineaire Theorie: Het is verrassend dat een lineaire ray-tracing theorie, gebaseerd op schaal-scheiding, de complexe niet-lineaire dynamiek van golfbreken en turbulentie-overgang kwalitatief correct voorspelt, zelfs wanneer de onderliggende aannames (kleine amplitude) worden geschonden.
3. Impact op Menging: De resultaten tonen aan dat de "geschiedenis" van de golf (hoe hij brak) de mengingsefficiëntie van de resulterende turbulentie bepaalt. Dit ondersteunt het concept van "history matters" in stratificatie-turbulentie.
4. Oceanografische Relevantie: De bevindingen suggereren dat horizontale schuifstromen (zoals die voorkomend in de thermocline of bij fronten) een belangrijke, maar vaak onderschatte bron kunnen zijn voor energiedissipatie en menging in de oceanen, met variabele efficiëntie afhankelijk van de golfparameters.

Conclusie

De auteurs tonen aan dat interne golven in horizontale schuifstromen via meerdere mechanismen kunnen breken, wat leidt tot een breed scala aan turbulente toestanden. De verhouding tussen kinetische en potentiële energie (gekwantificeerd door $F$) bepaalt of de instabiliteit convectief of schuif-gedreven is, wat op zijn beurt de mengingsefficiëntie en de totale energiedissipatie bepaalt. Deze inzichten zijn essentieerd voor het verbeteren van parameterisaties van menging in grootschalige oceanische en atmosferische modellen.