Application of Metric-Based Mesh Adaptation to Hypersonic Aerothermal Simulations Using US3D

Ce papier démontre l'efficacité de l'adaptation de maillage basée sur des métriques pour les simulations aérothermiques hypersoniques en gaz réel, prouvant qu'elle permet de prédire avec précision le chauffage de surface et de gérer des géométries complexes, telles que celles d'un capsule d'entrée atmosphérique avec jets RCS, là où les approches structurées traditionnelles échouent.

L'essentiel

🚀 Le "Zoom Intelligent" pour les vaisseaux spatiaux : Comment éviter de brûler en rentrant dans l'atmosphère

Imaginez que vous devez dessiner une carte très précise d'un territoire montagneux pour guider un pilote à travers une tempête. Si votre carte est floue, le pilote risque de s'écraser. C'est exactement le problème que rencontrent les ingénieurs de la NASA lorsqu'ils essaient de simuler la rentrée d'un vaisseau spatial dans l'atmosphère d'une planète (comme Mars).

Ce rapport, écrit par Dirk Ekelschot, raconte comment ils ont développé un outil magique pour rendre ces cartes (les simulations informatiques) beaucoup plus précises, même pour des vaisseaux très complexes.

1. Le Problème : La "Carte Floue" du Choc Thermique

Quand un vaisseau fonce à une vitesse hypersonique (plus de 20 fois la vitesse du son), il crée une onde de choc immense devant lui, comme le sillage d'un bateau, mais faite d'air surchauffé. Cette onde de choc est le point le plus critique : c'est là que la chaleur est la plus intense.

• L'ancienne méthode (Les blocs de Lego) : Pour calculer la chaleur, les ingénieurs utilisaient des grilles rigides, comme des blocs de Lego empilés. C'est très précis si le vaisseau est une forme simple (comme une sphère parfaite). Mais si le vaisseau a des tuyaux, des antennes ou des jets de contrôle (les petits propulseurs sur le dos du vaisseau), ces blocs de Lego deviennent un cauchemar à assembler. On est obligé de simplifier la forme du vaisseau, ce qui fausse les résultats.
• Le résultat : Souvent, ces grilles rigides ne "voient" pas bien l'onde de choc. C'est comme essayer de dessiner une ligne droite très fine avec un pinceau trop gros : le résultat est flou et imprécis.

2. La Solution : Le "Zoom Intelligent" (Adaptation de Maillage)

L'équipe a utilisé un nouveau système appelé MIMIC. Imaginez que vous avez une carte numérique qui peut se transformer toute seule.

• Le principe : Au lieu de dessiner toute la carte avec la même précision (ce qui prendrait des siècles de calcul), le logiciel regarde la simulation en temps réel. Là où il y a des changements brutaux (comme l'onde de choc ou les tourbillons derrière le vaisseau), le logiciel dit : "Attends, il faut plus de détails ici !" et il zoome automatiquement en ajoutant des milliers de petits points. Là où tout est calme, il dézoome pour économiser de la puissance.
• L'astuce de génie : Le rapport montre que pour que ce zoom fonctionne parfaitement près de la paroi du vaisseau (là où la chaleur est mesurée), il faut utiliser des formes géométriques spécifiques (des "hexaèdres", qui ressemblent à des briques parfaites) plutôt que des formes plus simples (des "prismes"). C'est comme si, pour peindre un tableau réaliste, il fallait utiliser des pinceaux plats et précis plutôt que des pinceaux ronds.

3. Les Deux Expériences de la NASA

Expérience A : La demi-sphère (Le test de base)
Ils ont pris une forme simple (une demi-sphère) pour tester leur méthode.

• Résultat : Ils ont comparé deux types de "briques" pour la couche près de la paroi. Ils ont découvert que les "briques" hexagonales (hexaèdres) donnaient des résultats de chaleur beaucoup plus lisses et précis que les autres, même avec le même nombre de points. C'est comme passer d'une photo granuleuse à une photo HD.

Expérience B : Le vaisseau martien complexe (Le vrai défi)
Ensuite, ils ont appliqué cela à un vrai vaisseau d'atterrissage martien (le Sample Return Lander). Ce vaisseau est complexe : il a un nez en forme de cône, mais aussi huit petits jets de contrôle (RCS) sur son dos, comme de petites fusées pour le guider.

• Le défi : Avec les anciennes méthodes rigides, il était presque impossible de modéliser ces petits jets sans rendre le calcul trop lourd.
• La réussite : Grâce à l'adaptation automatique, le logiciel a pu "zoomer" intelligemment autour de ces petits jets et dans la traînée turbulente derrière le vaisseau.
• La comparaison : Ils ont comparé leur résultat avec une simulation de référence très lourde (DPLR). Résultat ? Leur méthode, bien que plus flexible et capable de gérer les petits détails complexes, a donné des prédictions de chaleur aussi précises que la méthode lourde et rigide.

4. Pourquoi c'est important pour nous ?

Ce travail est une révolution pour deux raisons :

1. La liberté de forme : Les ingénieurs ne sont plus obligés de simplifier la forme des vaisseaux. Ils peuvent simuler des vaisseaux avec tous leurs détails (tuyaux, antennes, jets) sans perdre de temps à créer des grilles manuelles complexes. C'est comme passer d'un dessin au trait à un modèle 3D réaliste.
2. La sécurité : En prédisant exactement où et combien le vaisseau va chauffer, on peut mieux concevoir les boucliers thermiques. Cela rend les missions vers Mars (et ailleurs) plus sûres.

En résumé :
Ce rapport nous dit que grâce à un logiciel qui "apprend" où il faut zoomer (MIMIC) et en utilisant les bons types de briques géométriques, la NASA peut maintenant simuler des rentrées atmosphériques complexes avec une précision de haute qualité, tout en simplifiant énormément le travail des ingénieurs. C'est un pas de géant vers l'exploration spatiale future !

Résumé technique

1. Problématique

Les simulations de dynamique des fluides computationnelle (CFD) pour les écoulements hypersoniques, essentielles pour la conception de véhicules d'entrée atmosphérique, reposent souvent sur des maillages structurés en blocs (éléments hexaédriques) pour prédire avec précision le chauffage de surface. Cependant, la génération de tels maillages devient extrêmement complexe, voire impossible, pour des géométries présentant des caractéristiques détaillées (comme les systèmes de contrôle de réaction ou RCS) ou des écoulements complexes (interactions onde de choc-couche limite, sillage turbulent).

Les maillages non structurés (tétraédriques) offrent une flexibilité géométrique supérieure mais souffrent traditionnellement de prédictions de flux de chaleur inexactes en régime hypersonique. Cela est dû à la génération d'entropie numérique excessive au niveau de l'onde de choc bow (détachée), qui se propage vers la surface et pollue les résultats. L'objectif de ce travail est de démontrer la viabilité de l'adaptation de maillage anisotrope basée sur des métriques pour résoudre ces problèmes, en combinant la flexibilité des maillages non structurés avec la précision requise pour l'aérothermodynamique.

2. Méthodologie

L'étude utilise une approche couplée entre deux outils principaux :

• US3D : Un solveur d'écoulement compressible avancé (méthode des volumes finis centrés) capable de gérer des mélanges de gaz réactifs (air, atmosphère martienne) et des modèles de turbulence (RANS, IDDES).
• MIMIC (Metric-Informed Mesh Improvement Capability) : Une capacité d'amélioration de maillage basée sur des métriques, développée pour adapter les maillages en fonction de l'erreur d'interpolation estimée.

Le processus d'adaptation :

1. Indicateur d'erreur : Le tenseur métrique est calculé à partir du Hessien (dérivées secondes) de la solution de température. Cela permet d'identifier les zones nécessitant un raffinement (chocs, couches de cisaillement) et d'orienter les éléments du maillage.
2. Maillage Hybride Fixe : Une couche limite structurée (prismes ou hexaèdres) est générée près de la paroi et reste fixe pendant l'adaptation. Le reste du volume est maillé avec des tétraèdres qui sont adaptés dynamiquement.
3. Adaptation itérative : Le maillage est raffiné/coarsé itérativement (jusqu'à 5 itérations) en se basant sur la solution de température. Le solveur ParMMG est utilisé pour effectuer l'adaptation anisotrope des tétraèdres en parallèle.
4. Cas d'étude :
   • Cas 1 : Écoulement supersonique sur une hémisphère (gaz parfait) pour comparer l'influence de la topologie de la couche limite (prismes vs hexaèdres).
   • Cas 2 : Écoulement hypersonique (mélange réactif $CO_2$-$N_2$) sur une capsule d'entrée sphéro-cône à $70^\circ$ avec des jets RCS (Reaction Control System) sur la coque arrière, simulant l'entrée martienne.

3. Contributions Clés

• Extension des capacités de MIMIC : Pour la première fois, MIMIC est appliqué à des écoulements réactifs hypersoniques et gère désormais des maillages de couche limite hybrides contenant des hexaèdres (extrudés depuis des quadrilatères) en plus des prismes (extrudés depuis des triangles).
• Gestion de géométries complexes : Démonstration de la capacité à intégrer des détails géométriques fins (comme les buses des jets RCS) dans une simulation aérothermique monolithique, ce qui est difficilement réalisable avec des maillages structurés en blocs.
• Stratégie de démarrage robuste : Mise en place d'une procédure d'adaptation initiale utilisant des flux inviscides du premier ordre pour éviter les effets de "carbuncle" (instabilités numériques au choc) avant de passer à des schémas du second ordre une fois le maillage suffisamment aligné avec le choc.
• Comparaison Prisme vs Hexaèdre : Analyse comparative montrant que l'utilisation d'hexaèdres dans la couche limite fixe, combinée à l'adaptation, produit des distributions de flux de chaleur plus lisses et symétriques que les prismes.

4. Résultats

• Cas de l'hémisphère :

  • L'adaptation de maillage basée sur le Hessien de la température améliore considérablement la définition du choc et la prédiction du flux de chaleur.
  • Les maillages avec une couche limite en hexaèdres (110 couches) après adaptation fournissent des résultats de flux de chaleur comparables aux données de référence DPLR (maillage structuré) et sont plus symétriques que ceux obtenus avec des prismes.
  • La sensibilité à l'épaisseur de la couche limite est confirmée : un nombre insuffisant de couches (80) conduit à une résolution insuffisante, même après adaptation.

• Cas de la capsule d'entrée (70°) :

  • Qualité de la solution : Après 5 itérations d'adaptation, les contours de température et de nombre de Mach montrent une définition nette de l'onde de choc bow et des couches de cisaillement, avec une réduction drastique du bruit numérique présent sur le maillage initial.
  • Prédiction du flux de chaleur : Les stries (artefacts numériques) présentes sur le maillage initial disparaissent presque totalement après la 5ème itération. Le flux de chaleur prédit par US3D (maillage adaptatif) est qualitativement très proche de celui obtenu avec DPLR (maillage structuré), malgré la différence de topologie.
  • Géométrie RCS : Le maillage adaptatif capture avec succès les détails des jets RCS sur la coque arrière, là où une approche structurée aurait dû simplifier la géométrie.
  • Phénomènes physiques : L'adaptation révèle des interactions instables complexes dans le sillage, notamment la rupture de la couche de cisaillement supérieure en trois couches distinctes et la croissance de la bulle de séparation derrière l'escalier arrière, des phénomènes difficiles à capturer avec des maillages structurés simples.

5. Signification et Conclusion

Ce travail démontre que l'adaptation de maillage basée sur des métriques, couplée à un solveur non structuré comme US3D, est une alternative viable et puissante aux méthodes traditionnelles basées sur des maillages structurés en blocs pour l'aérothermodynamique hypersonique.

Points majeurs de l'impact :

1. Flexibilité géométrique : Permet de simuler des véhicules d'entrée avec des détails complexes (RCS, capteurs, etc.) sans simplification excessive de la géométrie.
2. Précision comparable : Avec une adaptation suffisante et une couche limite appropriée (hexaèdres), la précision des flux de chaleur atteint le niveau des codes structurés de référence (DPLR).
3. Automatisation : Réduit la charge de travail liée à la génération manuelle de maillages structurés complexes, ouvrant la voie à des études paramétriques plus rapides et à l'optimisation de formes complexes.

L'auteur conclut que cette approche est essentielle pour les futures missions d'exploration planétaire et suggère des travaux futurs incluant l'activation des jets RCS pour étudier les interactions plume-sillage et l'intégration avec d'autres solveurs (comme HyperSolve) pour améliorer encore les schémas numériques sur maillages anisotropes.