Two-Dimensional $\beta$-plane Turbulence: Dual Cascade and Zonal Jets

Résumé technique : Turbulence sur un plan $\beta$ bidimensionnel : Cascade duale et jets zonaux

Énoncé du problème
L'article étudie les propriétés statistiques des équations de Navier-Stokes bidimensionnelles, incompressibles, forcées et dissipatives, soumises à la force de Coriolis sous l'approximation d'un plan $\beta$. Ce système modélise les écoulements géophysiques (tels que les courants océaniques et les jets atmosphériques) où la rotation planétaire varie linéairement avec la coordonnée méridienne ($f(x) = f_0 + \beta x_2$). Les auteurs visent à caractériser rigoureusement deux phénomènes primordiaux dans ce régime :

1. Équilibre géostrophique : La relation statistique entre la force de Coriolis et le gradient de pression aux grandes échelles.
2. Cascade duale : Le comportement du transfert d'énergie et d'enstrophie à travers les échelles en présence de rotation, en examinant spécifiquement si la force de Coriolis modifie les lois universelles de la turbulence bidimensionnelle (cascade d'énergie inverse et cascade d'enstrophie directe) dérivées pour les fluides non rotatifs.

Le cadre mathématique implique un canal périodique $D = \mathbb{T} \times [a, b]$ avec des conditions aux limites de non-glissement dans la direction méridienne, piloté par un forçage stochastique et soumis à un amortissement d'Ekman à grande échelle ($\alpha$) et une viscosité à petite échelle ($\nu$).

Méthodologie
Les auteurs emploient un cadre probabiliste rigoureux basé sur les relations de Kármán–Howarth–Monin (KHM), adaptées aux équations de Navier-Stokes stochastiques sur un plan $\beta$. Les composantes méthodologiques clés incluent :

• Mesure stationnaire et régularité : L'analyse repose sur l'existence d'une mesure stationnaire pour le système stochastique, établie dans un travail antérieur [11]. Les auteurs supposent une régularité suffisante du bruit pour garantir que la mesure stationnaire possède une régularité $H^2$, permettant la manipulation des fonctions de corrélation à deux points.
• Relations KHM faibles : Au lieu de supposer l'homogénéité spatiale ou l'isotropie, les auteurs utilisent des formulations faibles des équations KHM impliquant des fonctions tests et des moyennes spatiales sur le domaine et la sphère unité.
• Projections duales :
  1. Projection de trace (symétrique) : Pour analyser la cascade duale, les auteurs prennent la trace de la relation tensorielle KHM. Cela permet de récupérer les équations de bilan scalaires pour l'énergie et l'enstrophie.
  2. Projection antisymétrique : Pour analyser l'équilibre géostrophique et les jets zonaux, les auteurs introduisent une nouvelle fonction test antisymétrique (une projection hors-diagonale pondérée). Cette projection spécifique isole les termes qui s'annulent dans les analyses isotropes standards, spécifiquement les contributions de Coriolis et du gradient de pression.
• Dissipation anormale faible (WAD) : Pour les résultats de cascade, les auteurs adoptent l'hypothèse WAD introduite par Bedrossian et al. [5]. Cette hypothèse postule que dans la limite non visqueuse ($\nu, \alpha \to 0$), l'énergie injectée est entièrement dissipée par la friction à grande échelle et l'enstrophie injectée par la viscosité, sans nécessendre de bornes uniformes sur les moments de vorticité plus fortes que nécessaire.

Contributions clés et résultats

1. Annulation des termes de Coriolis moyennés sphériquement (La cascade duale)
L'article prouve que les contributions moyennées sphériquement de la force de Coriolis aux flux d'énergie et d'enstrophie sont identiquement nulles à toutes les échelles.

• Théorème 3.1 : Il est démontré que la trace du tenseur de corrélation vitesse-Coriolis, $\Theta(l)$, et de la corrélation d'enstrophie de Coriolis, $Q(l)$, sont exactement nulles pour toute échelle de séparation $l > 0$.
• Théorèmes 3.2 & 3.3 (Lois de la cascade duale) : Sous l'hypothèse WAD, les auteurs dérivent le comportement asymptotique des fonctions de structure du troisième ordre moyennées.
  • Cascade d'enstrophie directe : Aux petites échelles ($l_\nu \ll l \ll l_I$), la fonction de structure mixte vitesse-vorticité suit la loi $S^{(1)}_3(l) \sim -2\eta l$, identique au cas non rotatif.
  • Cascade d'énergie inverse : Aux grandes échelles ($l_I \ll l \ll l_\alpha$), la fonction de structure de la vitesse suit la loi $S^{(3)}_3(l) \sim 2\varepsilon l$.
• Signification : Ces résultats démontrent que le taux net de transfert d'énergie et d'enstrophie dans le sens isotrope, moyenné sphériquement, est inexistant de l'influence de la force de Coriolis. La force de Coriolis n'injecte ni ne dissipe d'énergie ou d'enstrophie dans la moyenne isotrope ; elle agit purement comme une contrainte géométrique.

2. Équilibre géostrophique exact et formation de jets zonaux
L'article dérive une identité statistique exacte inédite reliant la force de Coriolis au gradient de pression via une fonction de corrélation à deux points.

• Théorème 3.4 (Équilibre KHM géostrophique) : En appliquant la projection antisymétrique à la relation KHM, les auteurs dérivent l'identité :
  $$\beta l \mathbb{E}\langle (n_2)^2 (u(x) \cdot u(x+ln)) \rangle = \mathbb{E}\langle n_2 (\nabla p(x) \cdot u^\perp(x+ln) - u^\perp(x) \cdot \nabla p(x+ln)) \rangle$$
  où $n$ est la direction du vecteur de séparation.
• Interprétation :
  • Aux petites échelles, le terme $\beta l$ est négligeable, restaurant l'isotropie locale.
  • Aux grandes échelles, le terme $\beta l$ croît de manière non bornée. Pour maintenir l'équilibre, la corrélation du gradient de pression doit croître proportionnellement.
  • Cette identité fournit une signature mathématique rigoureuse de la façon dont la force de Coriolis induit l'anisotropie. Le bilan force une réorganisation de l'écoulement où les fluctuations de vitesse méridienne sont supprimées, menant à la formation de jets zonaux (bandes d'écoulement est-ouest).
• Échelle de Rhines : Les auteurs montrent, par analyse dimensionnelle (Annexe A), que cette identité retrouve l'échelle de Rhines $l_R \sim (\varepsilon/\beta^3)^{1/5}$, l'échelle de transition où la cascade inverse est arrêtée par l'effet $\beta$.

3. Limitations sur les lois longitudinales
L'article note que bien que les lois de cascade moyennées sphériquement (isotropes) soient inchangées, les auteurs ne peuvent actuellement pas prouver les lois des fonctions de structure du troisième ordre longitudinales (par exemple, les analogues de la loi 4/5) en présence de la force de Coriolis. La contribution de Coriolis à la projection longitudinale ne s'annule pas et introduit des termes d'erreur qui divergent dans la limite non visqueuse sous les hypothèses WAD standards.

Signification et revendications
Les auteurs positionnent ce travail comme fournissant une « image claire » du rôle de la force de Coriolis dans la turbulence sur un plan $\beta$ en réconciliant deux observations physiques apparemment contradictoires :

1. Conservation des taux de cascade : La force de Coriolis ne modifie pas l'amplitude de la cascade duale (taux de transfert d'énergie et d'enstrophie) lorsqu'elle est moyennée de manière isotrope.
2. Induction d'anisotropie : La force de Coriolis modifie fondamentalement la distribution spatiale de l'énergie, pilotant la formation de structures zonales à grande échelle.

L'article prétend fournir le premier équilibre statistique exact pour la turbulence sur un plan $\beta$ qui capture l'émergence des jets zonaux sans reposer sur des approximations macroscopiques ou des fermetures statistiques. L'identité dérivée (Théorème 3.4) sert de signature mathématique rigoureuse à deux points de l'équilibre géostrophique, liant directement le paramètre de Coriolis aux corrélations pression-vitesse.

Ce travail complète les résultats précédents sur les petites échelles [10] en étendant la compréhension rigoureuse aux grandes échelles et au régime de transition, confirmant que si le mécanisme de la cascade (taux de transfert) reste robuste à la rotation, la structure de l'écoulement est profondément réorganisée par celle-ci.