Two-point Turbulence Closures in Physical Space

Ce papier présente un cadre de fermeture statistique prédictif à deux points pour la turbulence formulé dans l'espace physique, qui évite les transformées de Fourier en utilisant des équations d'évolution discrètes et des fermetures non linéaires inspirées de l'EDQNM pour modéliser avec précision la turbulence homogène isotrope tout en offrant des avantages pour les écoulements inhomogènes et anisotropes où les méthodes spectrales sont mal conditionnées.

L'essentiel

Imaginez que vous essayez de prédire comment une casserole d'eau bouillante se déplace, ou comment la fumée tourbillonne depuis une cheminée. C'est le monde de la turbulence. C'est chaotique, désordonné, et incroyablement difficile à prédire car chaque minuscule tourbillon d'eau ou d'air affecte tous les autres tourbillons autour de lui.

Pendant des décennies, les scientifiques ont tenté de créer des « livres de règles » (modèles mathématiques) pour prédire ce chaos. Les livres de règles les plus réussis jusqu'à présent ont été rédigés dans un langage spécial appelé l'Espace Spectral. Considérez l'Espace Spectral comme regarder une peinture complexe à travers un prisme : au lieu de voir les coups de pinceau, vous voyez les couleurs spécifiques (fréquences) qui composent l'image. C'est excellent pour les choses lisses et uniformes, mais si la peinture présente des bords nets, des fissures ou des changements soudains (comme une onde de choc dans un avion supersonique), le prisme se brise et l'image devient floue.

Cet article introduit une nouvelle façon de rédiger le livre de règles. Au lieu d'utiliser le prisme (Espace Spectral), les auteurs écrivent les règles directement dans l'Espace Physique — la vue réelle, du monde réel, où vous pouvez voir l'eau et l'air.

Voici une décomposition de leur approche utilisant des analogies simples :

1. Le Problème : L'Énigme des « Trop de Variables »

Dans la turbulence, pour prédire comment un tourbillon spécifique se déplacera ensuite, vous devez savoir comment il interagit avec tous ses voisins.

• L'Ancienne Méthode (Point Unique) : Les scientifiques regardaient auparavant une seule petite goutte d'eau et devinaient ce que faisaient ses voisins en se basant sur des règles moyennes. C'est comme essayer de prédire le trafic dans une ville en regardant seulement une voiture et en devinant le comportement de toute l'autoroute. Cela échoue souvent car cela manque la vue d'ensemble.
• La Solution à Deux Points : Les auteurs ont décidé de regarder deux points à la fois. Imaginez tendre les deux mains ; vous pouvez sentir la tension et la distance entre elles. En étudiant la relation entre deux points dans le fluide, ils peuvent capturer comment l'énergie passe d'un tourbillon à un autre avec beaucoup plus de précision.

2. L'Innovation : Marcher au lieu de Voler

La plupart des modèles de turbulence avancés (comme le célèbre modèle EDQNM) reposent sur la méthode du « prisme » (transformées de Fourier) pour faire leurs calculs. C'est rapide et élégant pour les écoulements lisses et uniformes.

• L'Astuce de l'Article : Les auteurs ont réalisé que si vous restez dans l'Espace Physique (le monde réel), vous n'avez pas besoin du prisme. Au lieu de survoler la ville pour voir toute la carte, ils ont décidé de marcher dans les rues.
• Comment ils l'ont fait : Ils ont utilisé une méthode appelée Différences Finies. Imaginez que vous voulez savoir à quel point une colline est raide. Au lieu d'utiliser une lunette magique, vous mesurez simplement la hauteur du sol à vos pieds et la hauteur du sol à quelques pas de distance. En faisant cela de manière répétée sur une grille, ils peuvent calculer comment le fluide se déplace sans jamais quitter « l'espace physique ».

3. L'« Amortissement des Tourbillons » (L'Amortisseur)

La turbulence est pleine d'énergie qui doit être dissipée (perdue sous forme de chaleur). Dans les anciens modèles, ils utilisaient un « amortisseur » (appelé amortissement des tourbillons) pour empêcher les mathématiques de devenir folles.

• Les auteurs ont dû inventer un nouveau type d'amortisseur qui fonctionne dans l'Espace Physique. Ils ont créé une « viscosité intelligente » qui agit comme une éponge, absorbant l'énergie chaotique exactement là où cela est nécessaire, en fonction des conditions locales de l'écoulement.

4. Le Problème de la Pression : La Force « Fantôme »

Dans les fluides, la pression agit instantanément partout. Si vous poussez l'eau ici, la pression change là immédiatement. C'est ce qu'on appelle un effet « non local ».

• Dans les anciens modèles de « prisme », cela était facile à gérer. Dans le nouveau modèle de « marche », c'est difficile. Les auteurs ont dû résoudre une énigme mathématique complexe impliquant des intégrales triples (imaginez calculer le poids total d'un nuage en additionnant chaque goutte de pluie individuelle dans une sphère 3D). Ils ont réussi à l'écrire dans leur nouveau langage, montrant que même si c'est lourd en calcul, c'est possible.

5. Est-ce que ça a marché ? (L'Essai Routier)

Les auteurs ont testé leur nouveau livre de règles « Espace Physique » contre deux choses :

1. L'Ancien Livre de Règles : Ils l'ont comparé aux meilleurs modèles spectraux pour la turbulence lisse et en décroissance (comme de la fumée qui s'estompe lentement). Résultat : Correspondance parfaite.
2. Données Réelles : Ils l'ont comparé à des simulations d'ordinateur ultra-puissantes (Simulations Numériques Directes) de turbulence forcée (comme un ventilateur soufflant de l'air). Résultat : Il a capturé le transfert d'énergie et la « tourbillonnance » de l'écoulement avec une grande précision.

Pourquoi est-ce important ? (Selon l'Article)

L'article affirme qu'il s'agit d'une preuve de concept. Cela prouve que vous pouvez construire un modèle de turbulence de haute précision sans utiliser le « prisme » (transformées de Fourier).

Les auteurs suggèrent que c'est une première étape cruciale pour aborder des problèmes plus difficiles où le prisme se brise, tels que :

• Écoulements compressibles : De l'air se déplaçant si vite qu'il crée des ondes de choc (comme un avion supersonique).
• Discontinuités : Des écoulements avec des sauts ou des ruptures soudains.

En Résumé :
Les auteurs ont construit une nouvelle façon robuste de prédire comment les fluides turbulents se déplacent en restant dans le « monde réel » (Espace Physique) plutôt que de traduire le problème dans un langage différent (Espace Spectral). Ils ont montré qu'en utilisant une approche basée sur une grille et des astuces mathématiques ingénieuses pour gérer la pression et la perte d'énergie, ils peuvent prédire la turbulence aussi bien que les anciennes méthodes, mais avec un cadre prêt à gérer les problèmes désordonnés et aux bords nets du monde réel.

Résumé technique

1. Énoncé du problème

La modélisation traditionnelle de la turbulence repose souvent sur des fermetures mono-points (par exemple, RANS), qui peinent à capturer avec précision les interactions non locales, l'anisotropie et le transfert d'énergie d'échelle en échelle. Bien que les fermetures bi-points (par exemple, EDQNM, DIA) répondent à ces limitations, elles ont historiquement été formulées presque exclusivement dans l'espace spectral (espace de Fourier).

La dépendance aux méthodes spectrales présente des obstacles majeurs pour les applications d'ingénierie complexes :

• Discontinuités : Les méthodes spectrales sont mal posées pour les écoulements comportant des chocs ou des discontinuités (courants dans les écoulements compressibles).
• Inhomogénéité : L'extension des fermetures spectrales aux écoulements inhomogènes est mathématiquement difficile et coûteuse en calcul.
• Conditions aux limites : Les méthodes spectrales éprouvent des difficultés avec des conditions aux limites complexes par rapport aux méthodes en espace physique.

Les auteurs visent à développer un cadre de fermeture statistique bi-points prédictif formulé entièrement en espace physique. L'objectif est de conserver la fidélité physique des statistiques bi-points (capturant la non-localité et le transfert d'énergie) tout en évitant les contraintes mathématiques des transformées de Fourier, permettant ainsi l'application à des écoulements complexes, discontinus et inhomogènes.

2. Méthodologie

Les auteurs dérivent un modèle de fermeture pour la Turbulence Isotrope Homogène Incompressible (HIT) comme preuve de concept, avec l'intention de généraliser le modèle aux écoulements anisotropes et inhomogènes.

A. Équations gouvernantes et discrétisation

• Formulation en espace physique : Au lieu de transformer les équations de Navier-Stokes en espace spectral, les auteurs travaillent directement avec les équations d'évolution des tenseurs de corrélation de vitesse bi-points ($R_{ij}$) et des moments d'ordre trois.
• Dérivées discrètes : Les dérivées spatiales sont approchées à l'aide de différences finies (spécifiquement RBF-FD, Radial Basis Function-Finite Difference, pour une haute précision près de $r=0$). Cela convertit les Équations aux Dérivées Partielles (EDP) des fonctions de corrélation en un système d'Équations Différentielles Ordinaires (EDO) sous forme matricielle-vectorielle.
  • Les termes linéaires (diffusion visqueuse) sont représentés par une matrice laplacienne radiale ($\mathbf{L}$).
  • L'évolution de la fonction de corrélation longitudinale $f(r)$ est régie par une EDO matricielle impliquant des matrices de dérivation ($\mathbf{D}$) et des matrices identité ($\mathbf{I}$).

B. Stratégie de fermeture (Adaptation EDQNM)

La fermeture suit les principes phénoménologiques du modèle Eddy-Damped Quasi-Normal Markovian (EDQNM) mais l'adapte à l'espace physique :

1. Quasi-normalité : Les moments d'ordre quatre sont approximés comme des produits de moments d'ordre deux (hypothèse gaussienne).
2. Approximation markovienne : L'historique temporel des moments d'ordre trois est approximé en évaluant les spectres au temps courant, permettant une intégration temporelle analytique.
3. Amortissement tourbillonnaire : Pour corriger le comportement non physique de l'approximation quasi-normale (manque de réversibilité et relaxation incorrecte), un terme de viscosité tourbillonnaire est introduit.
   • En espace physique, cela est mis en œuvre en modifiant la viscosité moléculaire dans la solution par exponentielle de matrice.
   • La viscosité tourbillonnaire est modélisée comme une fonction de l'incrément de vitesse local (via la fonction de structure d'ordre deux), assurant des taux d'amortissement corrects aux petites échelles.

C. Traitement de la pression non locale

Un défi unique en espace physique est l'équation de Poisson pour la pression, qui introduit une non-localité via une intégrale triple sur le domaine.

• Les auteurs dérivent la contribution de la pression à l'évolution du moment d'ordre trois.
• Ils démontrent que, bien que le terme de pression implique des intégrales triples complexes (nécessitant une quadrature numérique), la contribution nette de la pression au transfert d'énergie (évolution de l'énergie spectrale) est nulle en raison de la nature solénoïdale de l'écoulement. Cela simplifie considérablement la fermeture, permettant de se concentrer sur les termes d'advection.

D. Implémentation numérique

• Maillage : Un maillage à espacement géométrique est utilisé pour la distance de corrélation $r$ afin de résoudre les gradients raides près de $r=0$ (où les dérivées sont critiques) tout en couvrant les grandes échelles.
• Intégration temporelle : Un schéma IMEX (Implicit-Explicit) est utilisé. Les termes visqueux raides sont traités implicitement (en utilisant des exponentielles de matrice ou Euler implicite), tandis que les termes non linéaires sont traités explicitement.
• Exponentielle de matrice : La solution linéaire implique le calcul de l'exponentielle de matrice de l'opérateur laplacien, résolu efficacement à l'aide de méthodes de sous-espaces de Krylov.

3. Contributions clés

1. Premier cadre EDQNM en espace physique : L'article présente la première dérivation systématique d'une fermeture de type EDQNM entièrement en espace physique, contournant le besoin de transformées de Fourier.
2. Formulation matricielle-vectorielle : L'évolution des statistiques de turbulence est formulée comme un problème d'algèbre linéaire ($\dot{\mathbf{f}} = \mathbf{A}\mathbf{f} + \mathbf{b}$), où les dérivées spatiales sont gérées via des matrices de dérivation. Cela intègre naturellement l'information sur les échelles de longueur non locales.
3. Traitement analytique de la pression : Les auteurs fournissent une dérivation rigoureuse montrant comment le terme de pression non local peut être géré en espace physique et prouvent son annulation dans le bilan énergétique, simplifiant ainsi la fermeture.
4. Amortissement tourbillonnaire en espace physique : Une nouvelle formulation pour l'amortissement tourbillonnaire est introduite en utilisant la matrice laplacienne et les fonctions de structure, remplaçant la viscosité dépendante du nombre d'onde spectral.

4. Résultats et vérification

Le modèle a été vérifié contre deux références :

1. HIT en décroissance (Comparaison avec EDQNM spectral) :
   • Le modèle en espace physique a réussi à reproduire la décroissance des fonctions de corrélation longitudinale ($f(r)$) et latérale ($g(r)$).
   • Il a correctement capturé l'échelle inertielle de Kolmogorov ($E(\kappa) \sim \kappa^{-5/3}$) et les lois de décroissance pour l'énergie cinétique turbulente et les échelles de longueur intégrales (spectres de Saffman et Batchelor).
   • De légères divergences ont été observées aux très grandes longueurs de corrélation en raison de la troncature du maillage et des oscillations de la transformée de Hankel, mais la gamme inertielle a été résolue avec précision.
2. HIT forcée (Comparaison avec DNS) :
   • Le modèle a été testé contre des données de Simulation Numérique Directe (DNS) pour une turbulence forcée statistiquement stationnaire ($Re_\lambda \approx 433$).
   • Le modèle en espace physique a prédit avec précision le spectre d'énergie moyenné dans le temps et les fonctions de structure longitudinales.
   • Il a correctement capturé l'échelle de Kolmogorov en $-5/3$ et l'amplitude correcte du moment d'ordre trois (transfert d'énergie), validant ainsi la fermeture non linéaire.

5. Importance et perspectives futures

• Combler le fossé : Ce travail comble le fossé entre la fidélité physique des fermetures bi-points et l'applicabilité des méthodes en espace physique. Il ouvre la voie à l'application de fermetures de turbulence avancées aux écoulements compressibles avec chocs, aux écoulements réactifs et aux géométries complexes où les méthodes spectrales échouent.
• Compromis computationnels : Bien que l'approche en espace physique évite les transformées de Fourier, elle introduit des coûts computationnels associés à l'évaluation d'intégrales triples (pour la pression) et à de grandes opérations matricielles. Les auteurs proposent d'utiliser des Méthodes Multipôles Rapides (FMM) et des techniques de matrices creuses pour atténuer ces coûts pour les écoulements inhomogènes/anisotropes généraux.
• Extensibilité : Le cadre est conçu pour être extensible. Les auteurs décrivent comment la méthode peut être généralisée aux écoulements anisotropes (en utilisant la décomposition SO(3) et les harmoniques sphériques) et aux écoulements inhomogènes (en utilisant des grilles grossières locales avec des motifs bi-points intégrés).

En conclusion, cet article établit un cadre viable et prédictif pour la fermeture de turbulence bi-points en espace physique, démontrant qu'un modélisation statistique de haute fidélité peut être réalisée sans transformations spectrales, ouvrant ainsi la voie à une modélisation de turbulence plus robuste dans des écoulements d'ingénierie complexes et réels.