Leveraging unstructured grids for direct numerical simulations of wall turbulence

Ce papier présente le {\eta}-maillage, un cadre de génération de maillages non structurés pour les simulations numériques directes de la turbulence pariétale qui adapte la taille des mailles à l'échelle locale de Kolmogorov, atteignant une précision comparable à celle des maillages cartésiens classiques tout en réduisant considérablement le coût computationnel, en particulier aux nombres de Reynolds élevés et sur des géométries complexes de rainures.

L'essentiel

Imaginez essayer de simuler comment l'air s'écoule au-dessus d'une voiture ou comment l'eau se déplace le long de la coque d'un navire. Pour le faire avec précision sur un ordinateur, les scientifiques utilisent une technique appelée Simulation Numérique Directe (DNS). Considérez la DNS comme la création d'un gigantesque microscope numérique en 3D qui découpe le fluide (air ou eau) en millions de petits cubes invisibles (une grille). L'ordinateur calcule ensuite comment chaque cube individuel se déplace et interagit avec ses voisins.

Le problème est que l'écoulement du fluide près d'une surface (comme le flanc d'un navire) est incroyablement chaotique et détaillé. Pour obtenir une image claire, vous avez besoin d'un nombre massif de ces petits cubes juste à côté de la surface. Cependant, à mesure que vous vous éloignez de la surface, le chaos s'atténue, et vous n'avez plus besoin de cubes aussi petits.

L'Ancienne Méthode : Le "Mur de Briques Rigide"

Traditionnellement, les scientifiques utilisaient une grille cartésienne. Imaginez construire un mur avec des briques identiques et rigides.

• Le Problème : Pour voir les détails minuscules près de la surface, vous devez rendre les briques du bas très petites. Mais comme ces briques sont rigides et connectées en ligne droite, vous êtes forcé d'utiliser ces mêmes briques minuscules jusqu'au sommet de votre mur, même là où les détails ne sont pas importants.
• Le Résultat : Vous vous retrouvez avec un mur composé de milliards de petites briques, dont la plupart sont inutiles. Cela rend la simulation informatique incroyablement lente et coûteuse, comme essayer de compter chaque grain de sable sur une plage juste pour mesurer la marée.

La Nouvelle Solution : Le "Filet Intelligent et Extensible"

Ce papier présente une nouvelle méthode appelée grille $\eta$ (grille eta). Au lieu de briques rigides, imaginez un filet de pêche intelligent et extensible.

• Fonctionnement : Les auteurs ont conçu un système où la taille des mailles du filet change automatiquement en fonction du niveau de détail nécessaire.
  • Près de la surface (La "Couche Interne") : Le filet possède des mailles très petites et serrées pour attraper les minuscules tourbillons chaotiques du fluide.
  • Plus loin (La "Couche Externe") :** À mesure que le fluide devient plus calme, le filet s'étire automatiquement, rendant les mailles beaucoup plus grandes.
• L'Ingrédient Secret : La taille de ces mailles est basée sur quelque chose appelé l'échelle de Kolmogorov (notée $\eta$). Considérez cela comme le "plus petit tourbillon possible" qui peut exister dans le fluide à une hauteur donnée. La nouvelle grille dit simplement : "Rendez la taille de la maille juste assez grande pour attraper le plus petit tourbillon à cette hauteur spécifique, et pas plus."

Pourquoi C'est Important

Les auteurs ont testé ce "filet intelligent" sur deux types différents de codes informatiques (l'un semblable à une méthode spectrale par éléments, l'autre à une méthode des volumes finis) et l'ont comparé à l'ancienne méthode de "briques rigides".

1. Précision : Les résultats étaient presque identiques. Le "filet intelligent" a capturé la physique aussi bien que les "briques rigides", avec une différence inférieure à 1 % sur des mesures clés comme la friction et la vitesse.
2. Économies Massives : C'est là que la magie opère.
   • Pour les surfaces lisses (comme un mur plat), la nouvelle grille a réduit le nombre de "briques" requises (points de grille) d'environ 90 % à haute vitesse.
   • Pour les surfaces rugueuses (comme un mur avec de minuscules rainures appelées "riblets" conçues pour réduire la traînée), les économies ont été encore plus dramatiques — jusqu'à 97 % de points de grille en moins.

L'Analogie du "Mur Rainuré"

Pour comprendre la partie des riblets, imaginez un mur couvert de minuscules rainures parallèles (comme la texture d'une balle de golf ou la peau d'un requin).

• L'Ancienne Méthode : Pour simuler cela, la méthode des briques rigides devait garder les briques minuscules partout parce que les rainures forçaient la grille à être fine jusqu'en haut. C'était comme essayer de compter chaque fil d'un pull, même les parties loin du tissu.
• La Nouvelle Méthode : Le "filet intelligent" sait qu'une fois que vous êtes à quelques centimètres au-dessus de ces minuscules rainures, l'écoulement redevient lisse. Il étire les mailles immédiatement au-dessus des rainures, ignorant les détails minuscules qui ne comptent plus.

La Conclusion

Les auteurs ont créé un cadre qui agit comme un objectif de zoom intelligent pour les simulations de fluides. Il concentre sa puissance de calcul exactement là où elle est nécessaire (près du mur) et se relâche là où elle ne l'est pas.

• Pour les Murs Lisses : Il augmente l'effort beaucoup plus lentement à mesure que la simulation grossit.
• Pour les Murs Rugueux : Il s'adapte encore mieux, rendant les simulations de surfaces complexes réduisant la traînée réalisables sur des ordinateurs qui ne pouvaient pas les gérer auparavant.

En bref, ils ont trouvé un moyen de faire le même travail de haute qualité avec une fraction de la puissance de calcul, transformant une tâche qui pourrait prendre un superordinateur un mois en une tâche réalisable en quelques jours.

Résumé technique

Résumé technique : Exploitation de maillages non structurés pour les simulations numériques directes de la turbulence pariétale

Énoncé du problème
La Simulation Numérique Directe (SND) de la turbulence pariétale est coûteuse en calcul, principalement en raison de l'exigence de résoudre l'échelle de longueur de Kolmogorov ($\eta$) sur l'ensemble du domaine. Les solveurs sur maillages cartésiens conventionnels, largement utilisés pour les surfaces lisses et non lisses (par exemple via des méthodes de frontière immergée), imposent généralement une résolution de maillage fixe ou lentement variable dans les directions streamwise ($\Delta x^+$) et spanwise ($\Delta z^+$) sur tout le domaine. Cette approche entraîne une sur-résolution significative dans les régions externes où l'échelle de Kolmogorov augmente, créant un goulot d'étranglement computationnel. Ceci est particulièrement aigu pour les simulations impliquant des surfaces complexes comme les rainures (riblets), où le maillage doit être extrêmement fin près de la surface pour résoudre les micro-canaux, pourtant cette résolution fine est maintenue inutilement sur toute l'épaisseur de la couche limite dans les cadres cartésiens.

Méthodologie
Les auteurs proposent un nouveau cadre de génération de maillages non structurés appelé le maillage $\eta$. Le principe fondamental consiste à définir les tailles de maillage dans la direction normale à la paroi ($\Delta y^+$) et dans la direction spanwise ($\Delta z^+$) proportionnelles à l'échelle locale de Kolmogorov ($\eta^+$), permettant un raffinement du maillage dans les régions d'écoulement externe où les échelles de turbulence sont plus grandes.

1. Formulation du maillage :

   • Couche interne : Une fine couche ($\sim 50$ unités visqueuses) adjacente à la paroi utilise des tailles de maillage échelonnées visqueusement similaires à la SND conventionnelle ($0.3 \lesssim \Delta y^+ \lesssim 4$, $\Delta z^+ \simeq 5$). Pour les surfaces non lisses (par exemple, les rainures), le maillage résout la plus petite longueur d'onde de surface ($\ell^+$) avec $\ell^+/30 \lesssim \Delta y^+, \Delta z^+ \lesssim 4$.
   • Région externe : Au-dessus de la couche interne, les tailles de maillage augmentent proportionnellement à l'échelle locale de Kolmogorov ($\Delta y^+ \simeq \Delta z^+ \simeq 2\eta^+$).
   • Lois d'échelle : Les auteurs affinent des ajustements semi-empiriques pour $\eta^+$ dans les régions logarithmique et de sillage pour l'écoulement turbulent en canal et les couches limites turbulentes (TBL) à gradient de pression nul (ZPG). Ils proposent des ajustements par lois de puissance par morceaux pour $\eta^+$ qui transitionnent à $y^+ \approx \delta^+/2$.
   • Adaptation aux rainures : Pour les géométries de rainures, le cadre introduit une « sous-couche de rainures » ($\delta_r^+$) remplie de cellules carrées pour résoudre les écoulements secondaires et les rouleaux de Kelvin-Helmholtz, transitionnant de manière fluide vers l'échelle basée sur $\eta$ ci-dessus.

2. Implémentation :

   • Le cadre utilise le package open-source Gmsh pour générer des maillages non structurés multi-blocs.
   • Les maillages sont testés à l'aide de deux solveurs distincts : SOD2D (Méthode des Éléments Spectraux) et OpenFOAM (Méthode des Volumes Finis).
   • Les simulations couvrent l'écoulement turbulent en canal et les TBL ZPG sur des parois lisses et diverses géométries de rainures (rainures triangulaires) jusqu'à un nombre de Reynolds de frottement $\delta_0^+ = 1000$.

Contributions clés

• Développement du cadre : Une stratégie robuste de génération de maillages non structurés multi-blocs qui adapte dynamiquement l'espacement du maillage à l'échelle locale de Kolmogorov, applicable aux surfaces lisses et complexes rugueuses.
• Lois d'échelle affinées : Ajustements semi-empiriques mis à jour pour l'échelle de Kolmogorov dans la région de sillage des TBL et des écoulements en canal, améliorant la précision du dimensionnement du maillage dans la couche externe.
• Indépendance vis-à-vis du solveur : Démonstration que le cadre de maillage $\eta$ produit des résultats cohérents et de haute fidélité à travers différents schémas de discrétisation numérique (MES et MVF).
• Formulation spécifique aux rainures : Une stratégie de maillage adaptée aux rainures qui résout la physique de la sous-couche (écoulements secondaires, rouleaux KH) tout en permettant un grossissement rapide dans l'écoulement externe, surmontant les limitations des approches cartésiennes/IB.

Résultats

• Précision : Le maillage $\eta$ produit des résultats avec une différence inférieure à 1 % par rapport aux maillages cartésiens conventionnels et aux données de référence SND/expérimentales. Les métriques incluent les coefficients de frottement pariétal ($C_f$), les profils de vitesse moyenne ($U^+$), les contraintes turbulentes et les spectrogrammes.
• Efficacité du maillage :
  • Parois lisses : Le nombre de points de maillage ($N_\eta$) évolue comme $\propto \delta_0^{+2.5}$, comparé à $\propto \delta_0^{+3.0}$ pour les maillages cartésiens avec étirement tangente hyperbolique. À $\delta_0^+ = 6000$, le maillage $\eta$ ne nécessite que $\sim 10\%$ des points de maillage d'un maillage cartésien standard ($N_\eta/N_{Tanh} \simeq 0.1$).
  • Rainures : Le gain d'efficacité est encore plus prononcé. Pour les rainures triangulaires réduisant la traînée, $N_\eta$ évolue comme $\propto \delta_0^{+2.0}$, tandis que les maillages cartésiens évoluent comme $\propto \delta_0^{+3.0}$. À $\delta_0^+ = 6000$, le maillage $\eta$ ne nécessite que $\sim 3\%$ des points de maillage d'un maillage cartésien ($N_\eta/N_{Tanh} \simeq 0.03$).
• Physique des rainures : Le cadre capture avec succès les mécanismes d'écoulement complexes au-dessus des rainures, y compris l'origine virtuelle, les écoulements secondaires et les rouleaux de Kelvin-Helmholtz, avec une convergence du maillage obtenue en utilisant des paramètres de maillage modérés ($C_y=2.0, C_z=2.5, y_{in}^+=50$).
• Validation TBL : La SND des TBL ZPG sur des rainures a montré une bonne concordance avec les données expérimentales (Baron & Quadrio ; Choi & Orchard) concernant les pourcentages de réduction de traînée et les profils de vitesse, malgré les écarts inhérents entre SND et expérimental dans les régions de sillage.

Importance
L'article affirme que le cadre de maillage $\eta$ offre une « économie de maillage énorme », rendant la SND à haut nombre de Reynolds de la turbulence pariétale sur des surfaces complexes réalisable en calcul. En alignant la résolution du maillage sur les échelles physiques locales (échelle de Kolmogorov) plutôt que de maintenir des résolutions fixes, la méthode réduit drastiquement le coût computationnel (points de maillage et contraintes de pas de temps) sans sacrifier la précision. Les auteurs soulignent que cette approche permet la faisabilité économique de la réalisation de SND de TBL sur des surfaces complexes à des nombres de Reynolds correspondant aux conditions expérimentales, une tâche précédemment entravée par le coût prohibitif des méthodes basées sur des maillages cartésiens. Ce travail valide l'efficacité des maillages non structurés pour les simulations de turbulence de haute fidélité, dépassant les contraintes des maillages cartésiens structurés.