Towards a full-scale version of Yakhot's model of strong turbulence

Ce papier présente une extension complète du modèle de turbulence forte de Yakhot, offrant des expressions fermées sans paramètres libres pour les fonctions de structure d'ordre deux et trois qui décrivent avec précision la transition entre les régimes inertiels et dissipatifs sur toute l'échelle des Reynolds.

L'essentiel

🌪️ Comprendre la turbulence : Du tourbillon géant à la gouttelette fine

Imaginez que vous regardez une rivière tumultueuse ou que vous observez la fumée d'un cigare qui s'élève. Ce mouvement désordonné, ce mélange chaotique, c'est ce que les scientifiques appellent la turbulence. C'est l'un des plus grands mystères non résolus de la physique.

Ce papier, écrit par Christoph Renner, tente de résoudre une partie de ce casse-tête. Il s'agit d'une mise à jour d'un modèle théorique existant (celui de Yakhot) pour le rendre capable de décrire la turbulence à toutes les échelles, du plus grand tourbillon jusqu'à la toute petite goutte d'eau qui finit par s'arrêter.

Voici comment cela fonctionne, étape par étape :

1. Le problème : Deux mondes qui ne se parlent pas

Pour étudier la turbulence, les scientifiques utilisent des "règles" mathématiques appelées fonctions de structure. On peut les voir comme des règles de mesure pour voir à quelle vitesse l'eau bouge à différentes distances l'une de l'autre.

• Le monde des géants (Échelle inertielle) : Quand on regarde les gros tourbillons, la physique est bien comprise. Il existe une loi célèbre (la loi des 4/5 de Kolmogorov) qui fonctionne très bien ici.
• Le monde des micro-organismes (Échelle dissipative) : Quand on regarde les tout petits tourbillons, la viscosité (la "glu" du fluide) prend le dessus et freine tout. Ici, les règles changent complètement.

Le problème, c'est que le modèle de Yakhot (l'ancien modèle) était excellent pour les "géants", mais il échouait totalement pour décrire la transition vers les "micro-organismes". C'est comme avoir une carte routière parfaite pour les autoroutes, mais qui devient illisible dès qu'on entre dans une ruelle de village.

2. La découverte : Une nouvelle "règle de connexion"

L'auteur a analysé des données expérimentales très précises (mesurées dans un jet d'hélium refroidi à des températures cryogéniques, un peu comme un laboratoire de haute technologie).

Il a découvert une relation cachée, une sorte de pont invisible entre les gros mouvements et les petits.

• L'analogie : Imaginez que vous essayez de prédire le comportement d'une foule. Vous savez comment les gens bougent en groupe (les gros tourbillons). L'auteur a trouvé une formule mathématique qui explique comment, à mesure que le groupe se disperse en individus isolés (les petits tourbillons), leur mouvement change de manière très spécifique.

Cette relation, appelée équation (22) dans le texte, est une loi de puissance simple (une règle mathématique élégante) qui relie la vitesse de changement des mouvements à la taille du tourbillon.

3. La solution : Un modèle "Tout-en-un"

En utilisant cette nouvelle règle, l'auteur a pu réparer le modèle de Yakhot. Il a introduit un nouveau paramètre, une sorte de seuil de transition (noté $\rho$).

• L'analogie du thermostat : Imaginez un thermostat qui gère la température d'une maison.
  • Quand il fait très froid (petites échelles), le chauffage (la viscosité) s'active et change tout.
  • Quand il fait doux (grandes échelles), le chauffage est éteint et la nature suit son cours.
  • Le nouveau modèle définit exactement où se trouve ce "thermostat". Il dit : "Jusqu'à cette taille précise, c'est la viscosité qui règne. Au-delà, c'est la turbulence libre."

Ce seuil dépend du nombre de Reynolds, qui est simplement une mesure de la "violence" de l'écoulement. Plus l'eau coule vite et fort, plus ce seuil est petit.

4. Le résultat : Une carte complète sans trous

Le plus impressionnant de ce papier, c'est que le nouveau modèle :

1. Ne contient aucun paramètre libre : On n'a pas besoin de "deviner" des chiffres pour que ça marche. Tout est calculé à partir de la physique de base et des données réelles.
2. Est universel : Il fonctionne aussi bien pour les tout petits tourbillons (là où l'énergie est dissipée en chaleur) que pour les très gros (là où l'énergie est injectée).
3. Correspond à la réalité : Quand on compare les prédictions du modèle avec les données réelles de l'expérience, les courbes se superposent parfaitement. C'est comme si le modèle prédisait exactement comment la fumée va se disperser, du cigare jusqu'à l'air ambiant.

En résumé

Ce papier est une avancée majeure car il réussit à coudre ensemble deux morceaux de tissu qui semblaient incompatibles : la théorie des grands tourbillons et la réalité des petits tourbillons.

Grâce à une observation astucieuse dans les données expérimentales, l'auteur a trouvé la "colle" mathématique manquante. Désormais, nous avons un modèle complet, précis et sans paramètres arbitraires qui décrit la vie d'un fluide turbulent de son début à sa fin, du plus grand au plus petit. C'est un pas de géant vers la compréhension définitive de ce phénomène chaotique qui nous entoure tous les jours.

Résumé technique

1. Problématique et Contexte

La turbulence pleinement développée, homogène et isotrope reste l'un des problèmes non résolus les plus complexes de la physique des fluides. L'étude de ce phénomène repose souvent sur l'analyse des fonctions de structure $S_n(l) = \langle v(l)^n \rangle$, qui décrivent les moments statistiques des incréments de vitesse $v(l)$ sur une distance $l$.

Le modèle existant de Yakhot (1998) offre une description élégante de la plage inertielle (échelles intermédiaires) et de la transition vers les grandes échelles, en reliant les exposants de scaling aux propriétés thermodynamiques. Cependant, ce modèle présente deux lacunes majeures :

1. Il ne décrit pas correctement la convergence des fonctions de structure d'ordre pair vers leurs niveaux de grandes échelles.
2. Il ne prend pas en compte les effets dissipatifs dominants aux très petites échelles (plage de dissipation), où les lois d'échelle diffèrent radicalement de celles de la plage inertielle.

L'objectif de cet article est de combler ces lacunes en étendant le modèle de Yakhot pour couvrir l'ensemble des échelles, de la dissipation visqueuse jusqu'aux grandes échelles du système, en se basant sur des données expérimentales.

2. Méthodologie

L'approche méthodologique combine l'analyse théorique des équations de Navier-Stokes, l'extension du modèle de Yakhot et une validation rigoureuse par des données expérimentales issues d'un jet d'hélium cryogénique axisymétrique (Reynolds $Re \approx 2.3 \times 10^5$).

Étapes clés de la méthodologie :

• Extension à grande échelle : L'auteur reprend une extension récente du modèle de Yakhot qui introduit des termes supplémentaires ($D$ et $C$) pour gérer la convergence vers les grandes échelles.
• Analyse des limites à petite échelle : En utilisant la loi des quatre-cinquièmes de Kolmogorov et les conditions de symétrie, l'auteur dérive les comportements asymptotiques attendus pour $S_2(r)$ et $S_3(r)$ lorsque $r \to 0$ (où $r = l/L$ est la distance normalisée).
• Hypothèses sur le terme de dissipation :
  • Hypothèse H1 : Le terme de correction $d(r)$ est universel et indépendant de l'ordre $n$.
  • Hypothèse H2 : Le terme $d(r)$ conserve sa valeur de grande échelle (1) sur la plupart des échelles et tend vers zéro uniquement aux échelles dissipatives, n'influençant pas la dynamique à petite échelle.
• Extraction empirique de relations : En soustrayant les prédictions du modèle étendu des données expérimentales, l'auteur calcule les résidus $R_n(r)$. Une analyse de ces résidus révèle une relation de puissance surprenante entre les dérivées spatiales des fonctions de structure d'ordre pair et les moments d'ordre impair suivant.
• Clôture du système : Pour les ordres 2 et 3, la loi des quatre-cinquièmes de Kolmogorov permet de fermer le système d'équations différentielles, permettant de dériver des solutions analytiques fermées.

3. Contributions Clés

Les principales contributions de ce travail sont :

1. Découverte d'une relation empirique universelle : L'auteur établit une relation fondamentale (Éq. 22) reliant les résidus des équations de structure d'ordre pair à la fonction de structure d'ordre impair suivant :
   $$-n R_n(r) / S_{n+1}(r) \approx \tau / r^2$$
   Cette relation, valable de la plage de dissipation jusqu'à la plage inertielle, introduit un paramètre $\tau$ (ou une échelle de longueur $\rho$).

2. Introduction d'une échelle de longueur caractéristique $\rho$ : Ce paramètre marque le point de transition (crossover) entre le régime de dissipation visqueuse et le régime inertiel. Il est exprimé comme une fonction du nombre de Reynolds ($Re$) et du taux de dissipation d'énergie sans dimension ($\epsilon$).

3. Modèles fermés pour les ordres 2 et 3 : L'article dérive des expressions analytiques complètes pour les fonctions de structure du second ($S_2$) et du troisième ($S_3$) ordre. Ces modèles :

   • Sont valables sur toutes les échelles (de la dissipation aux grandes échelles).
   • Ne contiennent aucun paramètre libre (tous sont déterminés par les conditions aux limites et les paramètres physiques du flux).
   • Respectent les limites asymptotiques attendues (scaling en $r^2$ pour $S_2$ et $r^3$ pour $S_3$ aux petites échelles).

4. Correction de la symétrie : Le modèle ajuste le terme $D$ pour restaurer la symétrie $p(-u, -r) = p(u, r)$, violée par les extensions précédentes à grande échelle.

4. Résultats Principaux

• Validation expérimentale : Les modèles dérivés pour $S_2(r)$ et $S_3(r)$ montrent un accord excellent avec les données expérimentales sur toute la gamme d'échelles résolues, de la dissipation jusqu'à l'échelle du système.
• Comportement aux petites échelles : Le modèle reproduit correctement la loi de scaling $S_2(r) \propto r^2$ attendue dans la plage de dissipation, validant ainsi l'hypothèse selon laquelle le terme $d(r)$ devient négligeable à ces échelles.
• Détermination de l'échelle de transition $\rho$ :
  • La valeur expérimentale de $\rho$ est d'environ $5.8 \times 10^{-4}$.
  • La valeur théorique déduite des conditions aux limites à grande échelle est d'environ $6.1 \times 10^{-4}$.
  • Cette échelle $\rho$ est liée à l'échelle de Kolmogorov $\eta$ mais est légèrement plus grande, marquant le début de la transition vers la dissipation plutôt que la dissipation pure.
• Limites du modèle :
  • La relation empirique (22) semble moins précise pour les nombres de Reynolds plus faibles ($Re < 10^4$), où les écarts par rapport à une loi de puissance simple deviennent significatifs.
  • Le modèle n'est pas encore auto-contenu pour des ordres $n > 3$ car il dépend de l'ordre suivant ($S_{n+1}$) et aucune relation analogue n'a été trouvée pour les ordres impairs.

5. Signification et Conclusion

Ce travail représente une avancée significative vers une théorie complète de la turbulence forte. En étendant le modèle de Yakhot pour inclure la physique de la dissipation, l'auteur fournit le premier modèle fermé, sans paramètre libre et valide sur toutes les échelles pour les fonctions de structure du second et troisième ordre.

La découverte d'une relation de puissance simple reliant les résidus aux moments d'ordre supérieur suggère une structure sous-jacente profonde dans la turbulence, potentiellement liée à l'interdépendance des incréments longitudinaux et transversaux. Bien que le modèle soit actuellement limité aux ordres 2 et 3 (grâce à la loi des quatre-cinquièmes), il pose les bases théoriques pour de futures extensions vers des ordres supérieurs et offre un outil puissant pour la modélisation et la simulation numérique de la turbulence dans des régimes variés.