Scaling laws and local enhancements of buoyancy flux in stratified turbulent flows

Grâce à des simulations numériques directes approfondies de flux turbulents stratifiés, cette étude révèle que le flux de flottabilité présente une forte intermittence et des statistiques non gaussiennes pilotées par des fluctuations de grande échelle et de longue durée, son comportement moyen sur le domaine suivant une loi logarithmique par rapport au nombre de Reynolds de flottabilité et étant fondamentalement lié aux instabilités convectives qui déclenchent des cycles de dissipation d'énergie de type rafales.

L'essentiel

Imaginez l'atmosphère et l'océan non pas comme des rivières fluides et lisses, mais comme des soupes chaotiques et bouillonnantes. Parfois, ces soupes sont stratifiées comme un parfait — chaud sur le dessus, froid sur le dessous. Cette stratification est appelée stratification. Habituellement, on pourrait penser que cela rend la soupe calme et ordonnée, comme un lac paisible. Mais cet article révèle un secret surprenant : même dans ces fluides stratifiés et calmes, il existe des « poussées » de chaos soudaines et violentes qui se produisent à de très grandes échelles, et pas seulement à de petites échelles.

Voici une décomposition de ce que les chercheurs ont découvert, en utilisant des analogies simples :

1. La « fête surprise » dans une pièce calme

Dans une turbulence normale et non stratifiée (comme un mixeur mélangeant un smoothie), le chaos se produit principalement à très petite échelle — de minuscules petits tourbillons. Mais dans les fluides stratifiés (comme l'océan ou la haute atmosphère), les chercheurs ont découvert que le chaos peut être grand.

Imaginez une bibliothèque calme. Habituellement, les gens chuchotent. Mais soudain, un cri massif et inattendu éclate d'un coin spécifique. Cet article montre que dans l'océan et le ciel, ces « cris » (poussées soudaines de mouvement vertical et de changement de température) se produisent fréquemment. Ce ne sont pas de simples ondulations, ce sont des événements localisés et massifs qui peuvent être aussi grands que le flux lui-même.

2. Le « embouteillage » de chaleur et de mouvement

Les scientifiques étudiaient spécifiquement le flux de flottabilité. Appelons cela la « poignée de main chaleur-mouvement ». Cela mesure la quantité de chaleur qui monte ou descend exactement au même moment que l'air ou l'eau qui monte ou descend.

• La découverte : Ils ont découvert que cette « poignée de main » est incroyablement erratique. Parfois, l'air chaud monte violemment, et parfois l'air froid descend violemment.
• L'analogie : Imaginez une autoroute très fréquentée. La plupart du temps, les voitures circulent à une vitesse constante. Mais dans cette étude, ils ont découvert qu'occasionnellement, un énorme carambolage se produit où les voitures (la chaleur et le mouvement) accélèrent ou s'arrêtent soudainement de manière complètement imprévisible et extrême. Ces événements sont si extrêmes que les « queues » statistiques des données sont énormes — ce qui signifie que les événements « bizarres » sont beaucoup plus fréquents que ce que les mathématiques prédisent habituellement.

3. La zone « Goldilocks » du chaos

Les chercheurs ont testé de nombreuses conditions différentes, modifiant la façon dont le fluide est « épais » (viscosité) et la force de la stratification (stratification). Ils ont trouvé une « zone de réglage parfait » ou zone Goldilocks où ces poussées extrêmes se produisent le plus souvent.

• Trop de stratification (Stratification forte) : Le fluide est comme un gel rigide. Il se contente de vibrer comme une corde de guitare (ondes) et ne se mélange pas beaucoup.
• Trop peu de stratification (Faible stratification) : Le fluide est comme un mixeur chaotique. Il mélange tout, mais les poussées ne sont pas aussi distinctes.
• Juste ce qu'il faut (Le point idéal) : Lorsque la stratification est modérée, le fluide devient instable. C'est comme une pile de blocs Jenga qui est presque stable mais prête à s'effondrer. Dans cette zone, les « cris » (les poussées extrêmes) sont les plus forts.

4. Le mécanisme de la « Dette Énergétique »

Pourquoi ces poussées se produisent-elles ? L'article propose un mécanisme simple : la Dette Énergétique.

Imaginez que vous avez deux comptes bancaires : un pour le « Mouvement de haut en bas » (Énergie Cinétique) et un pour le « Potentiel de Chaleur » (Énergie Potentielle).

• Dans un monde parfait, ces comptes seraient équilibrés.
• Dans ces flux, ils se déséquilibrent. Le compte « Haut et Bas » devient trop élevé par rapport au compte « Chaleur ».
• La nature déteste la dette. Pour corriger ce déséquilibre, le système déverse soudainement toute cette énergie supplémentaire dans une poussée massive et violente. Cela crée un tourbillon ou un courant qui mélange les couches, rembourse la dette, puis le cycle recommence de zéro.

5. Ce que cela signifie pour la « vue d'ensemble »

L'article ne prétend pas résoudre le changement climatique ou prédire les ouragans directement. Au lieu de cela, il fournit un nouveau livre de règles pour comprendre comment ces fluides se comportent.

• La règle : L'intensité de ces poussées chaotiques suit un modèle mathématique spécifique (une loi de puissance) basé sur la force de la stratification.
• La conclusion : Même dans un océan ou un ciel stratifié, vous ne pouvez pas supposer que le flux est lisse. Il existe des « points chauds » cachés et massifs de mélange qui se produisent par poussées. Si vous essayez de modéliser la façon dont la chaleur ou la pollution se déplace dans l'océan ou l'atmosphère, vous devez tenir compte de ces pics soudains et violents, et pas seulement du flux moyen.

Résumé

Cet article est comme la découverte qu'un lac calme n'est pas réellement calme. Il est rempli d'explosions sous-marines massives et cachées qui se produisent dans une « zone Goldilocks » spécifique de stabilité. Ces explosions sont entraînées par un déséquilibre entre le mouvement et la chaleur, et elles suivent un rythme mathématique prévisible. Comprendre ce rythme nous aide à réaliser que l'atmosphère et les océans sont bien plus « irréguliers » et imprévisibles que nous ne le pensions auparavant.

Résumé technique

Titre : Lois d'échelle et renforcements locaux du flux de flottabilité dans les écoulements turbulents stratifiés

Énoncé du problème
Les écoulements turbulents stratifiés, prévalents dans l'atmosphère et les océans, présentent une intermittence non seulement aux petites échelles, mais aussi aux grandes échelles comparables à l'écoulement moyen. Alors que la turbulence homogène isotrope (HIT) pleinement développée est connue pour être intermittente aux petites échelles, des preuves suggèrent que dans les écoulements stratifiés, le champ de vitesse peut lui-même présenter des statistiques non gaussiennes et une intermittence à « grande échelle ». Ce phénomène implique l'émergence de puissants courants de vitesse verticale qui agissent comme des mécanismes d'injection d'énergie locale, conduisant à des propriétés de transport hautement variables le long de la direction de la gravité. La dynamique spécifique régissant ces événements à grande échelle, en particulier le flux de flottabilité ($B_f$), et sa dépendance vis-à-vis du nombre de Reynolds de flottabilité ($Ri_B$) et du nombre de Prandtl ($Pr$), nécessite une élucidation plus approfondie pour comprendre le mélange et la dissipation dans les contextes géophysiques.

Méthodologie
Les auteurs ont réalisé une large exploration paramétrique à l'aide de simulations numériques directes (DNS) des équations de Boussinesq. L'étude a utilisé le code pseudo-spectral GHOST avec une parallélisation hybride (MPI, OpenMP, CUDA).

• Équations directrices : Les équations de Navier-Stokes incompressibles couplées à une équation de température, incluant des termes de flottabilité caractérisés par la fréquence de Brunt-Väisälä ($N$).
• Paramètres : L'étude a fait varier le nombre de Froude ($Fr$) sur une plage géophysiquement pertinente ($0,01 \le Fr \le 1$), correspondant à des nombres de Reynolds de flottabilité ($Ri_B$) allant de $0,06$à$2300$. Deux nombres de Prandtl ont été examinés : $Pr = 1$ et $Pr = 6$.
• Schémas de forçage : Deux types de forçage ont été employés pour le champ de vitesse : un forçage de Pouquet-Patterson (PP) aléatoire et un forçage de Taylor-Green (TG).
• Portée de la simulation : Un total de 60 calculs a été analysé, incluant des séries avec $Pr=1$ (PP-Pr1, TG-Pr1), $Pr=6$ (PP-Pr6), et des simulations de scalaire passif ($N=0$). Les résolutions allaient de $128^3$ à $512^3$ points de grille.
• Analyse : L'étude s'est concentrée sur les statistiques d'ordre élevé (asymétrie/skewness et kurtosis) de la vitesse verticale ($w$), des fluctuations de température ($\theta$) et du flux de flottabilité ($B_f = Nw\theta$). L'analyse a inclus des moyennes volumiques et temporelles, ainsi que l'examen des fonctions de densité de probabilité (PDF) et des PDF conjointes.

Contributions clés et résultats

1. Intermittence à grande échelle du flux de flottabilité :
   L'étude confirme que le flux de flottabilité ($B_f$) présente des queues fortement non gaussiennes, avec des valeurs de kurtosis atteignant environ $10^2$. Cela indique que les écoulements géophysiques stratifiés peuvent être caractérisés par des propriétés de transport hautement variables, même sous stratification stable. Cette intermittence est associée à un comportement intermittent de longue durée et de grande échelle de la vitesse verticale et de la température.

2. Lois d'échelle et dépendance à $Ri_B$ :

• Asymétrie (Skewness) : L'asymétrie du flux de flottabilité ($S_{B_f}$) augmente avec $Ri_B$ selon une loi de puissance ($S_{B_f} \propto Ri_B^{0,28}$) à mesure que la stratification s'affaiblit. Cette tendance sature dans la limite du scalaire passif ($N \to 0$), où la température se découple du champ de vitesse.
• Pics de Kurtosis : Les maxima pour le kurtosis de $w$, $\theta$ et $B_f$ se produisent dans une plage étroite de $1 < Ri_B < 2$. Cette plage correspond à un régime spécifique de l'équilibre onde-tourbillon.
• Flux de flottabilité moyen : Le flux de flottabilité moyen du domaine $\langle B_f \rangle$ présente deux tendances distinctes : il suit une loi logarithmique avec $Ri_B$ ($\langle B_f \rangle \propto \ln(Ri_B)$) dans le régime intermédiaire, et approche un petit décalage lorsque la stratification se renforce ($Ri_B \to 0$).

3. Rôle du nombre de Prandtl :
   La variation du nombre de Prandtl ($Pr = 1$ vs. $Pr = 6$) ne semble pas affecter de manière significative le développement de courants verticaux extrêmes ou l'émergence de l'intermittence à grande échelle, du moins dans les valeurs considérées. La forme des courbes d'échelle pour les différentes séries de $Pr$ est presque identique, ne différant que par leur magnitude.

4. Mécanisme des événements extrêmes :
   Un modèle simple pour l'évolution temporelle de l'énergie et du flux de flottabilité indique que le « déficit d'énergie » entre l'énergie cinétique verticale et l'énergie potentielle pilote les événements forts de la flottabilité.

• Le modèle (basé sur la dynamique de Vieillefosse étendue aux écoulements stratifiés) suggère qu'un manque d'équipartition entre l'énergie verticale et l'énergie potentielle rend le flux de flottabilité sensible à la dynamique locale des ondes.
• Ce déséquilibre énergétique conduit à des instabilités convectives, à la formation de tourbillons 2D et 3D, et à une dissipation rapide sur l'échelle d'un temps de rotation, permettant au cycle énergétique de redémarrer par bouffées.
• Les PDF conjointes révèlent qu'au pic de l'intermittence ($Ri_B \approx 1,78$), la vitesse verticale peut atteindre jusqu'à $15\sigma_w$, tandis que la vitesse horizontale reste limitée, soulignant la nature anisotrope de ces événements.

5. Efficacité du mélange :
   L'étude analyse l'efficacité du mélange ($\langle B_f \rangle / \langle \epsilon_V \rangle$) et l'efficacité du mélange irréversible ($\hat{\Gamma}$). Ces deux quantités diminuent à mesure que la stratification s'accentue jusqu'à ce qu'elles saturent à bas $Ri_B$. Les auteurs notent que bien que les événements extrêmes affectent considérablement les flux locaux, ils n'altèrent pas globalement l'efficacité du mélange lorsqu'ils sont moyennés sur le volume et le temps de l'écoulement. Cependant, les valeurs locales de $B_f$ peuvent différer d'un ordre de grandeur entre les régions avec des événements extrêmes et celles sans eux.

Signification
L'article affirme que l'identification de ces lois d'échelle et des conditions spécifiques déclenchant les événements extrêmes est essentielle pour améliorer la modélisation et la paramétrisation des écoulements géophysiques et climatologiques.

• Prédictibilité : Les résultats suggèrent que l'intensification rapide de phénomènes tels que les ouragans pourrait être mieux prédite en suivant les variations de la vitesse verticale, étant donné le lien établi entre les courants verticaux et l'intermittence à grande échelle.
• Pertinence géophysique : La plage de $Ri_B$ et de $Fr$ où ces événements extrêmes se produisent est pertinente tant pour l'atmosphère que pour les océans.
• Modélisation : Comprendre les renforcements locaux du flux de flottabilité et le mécanisme de déficit d'énergie fournit une base pour affiner les modèles de sous-maille et les paramétrisations utilisés dans les modèles climatiques et océaniques, notamment concernant la manière dont la turbulence et les ondes interagissent pour piloter le mélange.

Les auteurs concluent que bien que l'efficacité globale du mélange puisse rester relativement stable, la nature locale et par bouffées de la dissipation et du transport pilotés par ces courants intermittents à grande échelle est une caractéristique critique de la turbulence stratifiée qui doit être prise en compte dans les applications géophysiques.