Instability and breaking of internal waves in a horizontal shear layer

Cette étude combine la théorie de la trajectoire des ondes et des simulations numériques directes pour analyser la rupture des ondes internes dans un cisaillement horizontal, identifiant deux mécanismes distincts de croissance d'énergie et démontrant que la dynamique de rupture, qui peut générer une dissipation turbulente dépassant l'énergie initiale, dépend fortement du rapport d'énergie de perturbation initial.

L'essentiel

🌊 Quand les vagues intérieures rencontrent un courant : Le grand bal des océans

Imaginez l'océan (ou l'atmosphère) non pas comme une masse d'eau calme, mais comme un gâteau à plusieurs étages. Les couches supérieures sont plus chaudes et légères, les couches inférieures plus froides et lourdes. C'est ce qu'on appelle la stratification.

Dans ce "gâteau", des vagues invisibles se déplacent à l'intérieur, entre les couches. On les appelle des vagues internes. Elles sont comme des secousses qui voyagent à travers le gâteau sans le briser immédiatement.

Mais parfois, ces vagues rencontrent un courant horizontal (un vent ou un courant marin qui coule d'un côté à l'autre). C'est là que l'histoire devient passionnante, et c'est ce que l'équipe de chercheurs a étudié.

1. Le Problème : Deux façons de faire "craquer" le gâteau

Lorsqu'une vague interne traverse ce courant, elle peut devenir instable et se briser, créant de la turbulence (comme du café que l'on remue). Les chercheurs ont découvert qu'il existe deux mécanismes principaux pour faire craquer ces vagues, un peu comme deux façons différentes de faire tomber une tour de cartes :

• Le mécanisme de la "Pente Raide" (Réfraction) :
  Imaginez que vous conduisez une voiture sur une route qui monte de plus en plus vite. Plus vous montez, plus vous devez ralentir, et votre voiture semble se "tasser". Pour la vague, le courant agit comme une pente qui la force à se comprimer verticalement. Elle devient de plus en plus raide, comme une vague à la surface de l'océan prête à déferler.

  • Résultat : La vague se brise parce qu'elle est trop pentue. C'est une instabilité convective (comme une bulle d'air chaud qui monte et éclate).

• Le mécanisme du "Tapis Roulant" (Advection) :
  Imaginez maintenant que vous êtes sur un tapis roulant qui accélère. Si vous tenez une corde, le mouvement du tapis va étirer la corde dans le sens du courant. Ici, le courant horizontal "pousse" la vague, étirant le mouvement de l'eau dans la direction du flux. Cela crée des cisaillements (des frottements) violents entre les couches d'eau.

  • Résultat : La vague se brise parce qu'elle est "tirée" et étirée par le courant. C'est une instabilité de cisaillement (comme des couches de beurre qui glissent les unes sur les autres).

2. La Recette Magique : Le Paramètre "F"

Les chercheurs ont créé une sorte de "recette" mathématique (un nombre sans unité appelé F) pour prédire quelle des deux méthodes va gagner.

• Si F est petit : C'est le mécanisme de la "Pente Raide" qui domine. La vague se brise en faisant des culbutes (comme une vague de surf). L'énergie est partagée équitablement entre le mouvement de l'eau et la chaleur (la flottabilité).
• Si F est grand : C'est le mécanisme du "Tapis Roulant" qui gagne. Le courant étire la vague, créant des courants violents. L'énergie est presque tout entière dans le mouvement de l'eau (cinétique).

C'est comme si vous aviez un bouton de réglage sur votre machine à vagues : en tournant le bouton (en changeant les conditions initiales), vous choisissez si la vague va se briser en "faisant des culbutes" ou en "se faisant étirer".

3. L'Expérience : La Cuisine Numérique

Pour vérifier cette théorie, les chercheurs n'ont pas utilisé de vrais océans (trop grands !), mais des supercalculateurs pour simuler des millions de litres d'eau virtuels. Ils ont lancé des vagues dans différents courants et ont regardé ce qui se passait.

Ce qu'ils ont découvert :

• La théorie fonctionne ! Même si les mathématiques utilisées pour la prédire sont simplifiées (comme une carte routière simplifiée), elles prédisent très bien le comportement réel de la vague, même quand celle-ci devient très complexe et turbulente.
• L'effet de surprise : Dans certains cas, la turbulence générée par la rupture de la vague libère beaucoup plus d'énergie que celle contenue dans la vague elle-même au départ. C'est comme si une petite étincelle (la vague) avait déclenché un incendie géant en allumant le carburant caché dans le courant (l'énergie du fond).
• Le mélange : La façon dont la vague se brise change la qualité du "mélange".
  • Si elle se brise en faisant des culbutes, elle mélange très bien les couches (comme bien fouetter un œuf).
  • Si elle se brise en étant étirée, elle mélange moins bien, mais crée des courants très forts.

4. Pourquoi est-ce important ?

Vous vous demandez peut-être : "À quoi ça sert de savoir comment une petite vague se brise dans un courant ?"

C'est crucial pour comprendre le climat et les océans !

• Le climat : Les océans absorbent la chaleur et le carbone. Pour que cela fonctionne, l'eau de surface doit se mélanger avec l'eau profonde. Ce mélange est souvent causé par la rupture de ces vagues internes.
• La précision : Si nos modèles climatiques ne comprennent pas comment ces vagues se brisent (culbutes ou étirement), ils ne peuvent pas prédire correctement la température de la Terre ou le niveau des mers.

En résumé

Cette étude nous dit que les vagues invisibles de l'océan ne sont pas de simples passagères. Lorsqu'elles rencontrent un courant, elles peuvent se transformer en deux types de "monstres" différents selon les conditions :

1. Des vagues qui se tordent (instabilité convective).
2. Des vagues qui sont étirées (instabilité de cisaillement).

Comprendre lequel des deux va dominer permet aux scientifiques de mieux prédire comment l'océan mélange la chaleur, l'oxygène et les nutriments, un peu comme un chef qui doit savoir exactement comment mélanger ses ingrédients pour obtenir le plat parfait.

Résumé technique

1. Problématique et Contexte

L'étude se concentre sur le comportement des ondes de gravité internes se propageant dans un écoulement de cisaillement horizontal, où la direction du cisaillement est orthogonale à la gravité. Ce scénario est pertinent pour la dynamique océanique et atmosphérique, où les ondes internes transportent une énergie significative et où leur dissipation via la turbulence à petite échelle est un mécanisme clé pour le mélange et le transport d'énergie.

Contrairement au cas classique du cisaillement vertical (où les niveaux critiques sont bien étudiés), le cisaillement horizontal présente des dynamiques distinctes. L'article vise à comprendre comment les interactions entre une onde monochromatique et un cisaillement horizontal (modélisé par un profil $U(y) = \Delta u \tanh(y/h)$) peuvent mener à l'instabilité et au déferlement (breaking) de l'onde. L'objectif est d'identifier les mécanismes de croissance de l'énergie des perturbations et de prédire les conditions menant à la turbulence.

2. Méthodologie

Les auteurs combinent une approche théorique linéaire et des simulations numériques directes (DNS) non linéaires complètes :

• Théorie de ray-tracing (WKBJ) :

  • Une analyse linéaire est effectuée sous l'approximation de WKB (Wentzel-Kramers-Brillouin-Jeffreys) pour décrire la propagation d'un paquet d'ondes dans le cisaillement.
  • Les équations de ray-tracing sont résolues pour suivre l'évolution du vecteur d'onde $\mathbf{k}$, de la fréquence intrinsèque $\Omega$ et de l'amplitude de l'onde (raideur $s$).
  • Deux mécanismes distincts de croissance d'énergie sont identifiés :
    1. Réfraction : Augmentation de la raideur de l'onde due à la modification du vecteur d'onde par le cisaillement (conduisant potentiellement à une instabilité convective).
    2. Advection de la quantité de mouvement : Augmentation des perturbations de vitesse dans la direction du cisaillement ($u'$) due à l'advection du flux moyen par l'onde (mécanisme de type "lift-up").
  • Un paramètre adimensionnel clé, le rapport d'énergie de perturbation $F$, est construit pour quantifier l'importance relative de l'advection de la quantité de mouvement par rapport à la raideur de l'onde. $F$ est lié à l'inverse du nombre de Richardson gradient local ($Ri_g$).

• Simulations Numériques Directes (DNS) :

  • Des simulations non linéaires sont réalisées avec le solveur pseudo-spectral Diablo.
  • Une gamme de paramètres est explorée pour couvrir différents régimes : onde piégée ($k_x < 0$), onde réfléchie ($k_x > 0$) et onde non affectée par la réfraction ($k_x = 0$).
  • Les simulations varient l'amplitude initiale de l'onde ($s_0$), l'intensité du cisaillement ($\Delta u$) et le nombre de Reynolds.
  • Un cadre de référence incliné est utilisé pour permettre l'étude de vecteurs d'onde avec de très petites composantes horizontales ($k_x \ll 1$).

3. Contributions Clés et Résultats

A. Prédiction des mécanismes d'instabilité via le paramètre $F$

Le paramètre $F$ permet de prédire le type d'instabilité qui se développera :

• Cas $F$ faible (et $s$ grand) : La croissance de l'énergie est dominée par l'augmentation de la raideur de l'onde (réfraction). L'énergie cinétique et potentielle restent approximativement équiparties. L'instabilité résultante est un mélange d'instabilité de cisaillement et d'instabilité convective (dûe au retournement des isopycnes).
• Cas $F$ grand : L'advection de la quantité de mouvement domine, entraînant une augmentation significative des perturbations de vitesse horizontale ($u'$) et donc du cisaillement vertical. Cela peut conduire à une instabilité de cisaillement pure (type Kelvin-Helmholtz) même si la raideur de l'onde reste faible ($s < 1$).

B. Dynamiques observées dans les DNS

Les simulations confirment qualitativement les prédictions théoriques malgré les écarts avec les hypothèses linéaires :

• Piégeage d'onde ($k_x < 0$) : L'onde est réfractée vers un niveau de piégeage où sa raideur augmente. Pour de grandes amplitudes, cela génère des retournements convectifs. Pour de faibles amplitudes, l'accélération du flux moyen (mécanisme de "self-acceleration") joue un rôle crucial dans le déclenchement de l'instabilité tridimensionnelle.
• Réflexion et passage ($k_x \ge 0$) : Lorsque $F$ est élevé, l'advection de la quantité de mouvement crée un cisaillement vertical intense, déclenchant des instabilités de type Kelvin-Helmholtz (billows) et une transition rapide vers la turbulence, même sans retournement convectif initial.

C. Bilans d'énergie et de quantité de mouvement

• Dissipation turbulente : La dissipation d'énergie peut dépasser largement l'énergie initiale de l'onde (jusqu'à 7,5 fois dans certains cas). Cela suggère que le déferlement de l'onde agit comme un déclencheur pour une transition sous-critique de l'écoulement de cisaillement de fond, qui est linéairement stable.
• Transfert de quantité de mouvement : Le déferlement entraîne un freinage significatif de l'écoulement de cisaillement moyen et une modification de l'épaisseur de la couche de cisaillement.
• Efficacité de mélange : L'efficacité de mélange (rapport entre la dissipation d'énergie potentielle disponible et la dissipation totale) varie considérablement selon le mécanisme de déferlement :
  • Déferlement convectif ($F$ faible) : Efficacité élevée ($\eta \approx 0.42$).
  • Déferlement par cisaillement ($F$ élevé) : Efficacité faible ($\eta \approx 0.08$).
• Émission d'ondes : La turbulence générée peut émettre de nouvelles ondes internes, dont l'énergie totale peut atteindre 25 % de l'énergie dissipée, un phénomène moins observé dans les cas de cisaillement vertical.

4. Signification et Implications

Cette étude apporte plusieurs avancées majeures pour la compréhension de la dynamique des fluides géophysiques :

1. Validation d'un modèle prédictif simple : La théorie de ray-tracing linéaire, bien que basée sur des hypothèses simplificatrices, parvient à prédire avec succès les régimes d'instabilité et la partition énergie cinétique/énergie potentielle dans des écoulements fortement non linéaires. Le paramètre $F$ s'avère être un indicateur robuste.
2. Rôle du cisaillement horizontal : L'article démontre que le cisaillement horizontal peut être un mécanisme de déferlement aussi important que le cisaillement vertical, en particulier via le mécanisme d'advection de la quantité de mouvement (lift-up), souvent négligé dans les modèles de paramétrisation océanique.
3. Variabilité du mélange : La découverte que l'efficacité de mélange dépend fortement du mécanisme de déferlement (convectif vs cisaillement) remet en question l'hypothèse d'une efficacité de mélange constante dans les modèles océaniques globaux.
4. Source d'énergie pour la turbulence : Le fait que la dissipation puisse dépasser l'énergie de l'onde incidente indique que les ondes internes peuvent servir de "déclencheurs" pour extraire l'énergie stockée dans les courants de fond stables, un processus crucial pour le budget énergétique de l'océan.

En conclusion, ce travail établit un lien clair entre les conditions initiales d'une onde interne dans un cisaillement horizontal et la nature de la turbulence qui en résulte, offrant des perspectives pour améliorer les paramétrisations de mélange dans les modèles climatiques et océanographiques.