HYMOR: An open-source package for global modal, non-modal, and receptivity analysis in high-enthalpy hypersonic vehicles

Le papier présente HYMOR, une boîte à outils open-source capable d'effectuer des analyses modales, non-modales et de réceptivité globales pour les écoulets hypersoniques à haute enthalpie, en intégrant des modèles de gaz réel et une formulation de choc-fitting pour capturer les interactions physiques complexes inaccessibles aux méthodes locales.

L'essentiel

🚀 HYMOR : Le "Simulateur de Vol" pour les avions supersoniques de demain

Imaginez que vous voulez construire un avion capable de voler à 12 fois la vitesse du son (hypersonique) à travers l'atmosphère terrestre. À cette vitesse, l'air autour de l'avion ne se comporte plus comme un gaz normal : il devient brûlant, ionisé et chimiquement actif, comme un four à micro-ondes géant en mouvement.

Le plus grand défi ? Comprendre comment l'air passe d'un écoulement calme (laminar) à un chaos turbulent. Si cela arrive au mauvais endroit, l'avion peut surchauffer et se désintégrer.

C'est là qu'intervient HYMOR. C'est un nouveau logiciel gratuit (open-source) créé par des chercheurs du Caltech et du MIT pour prédire exactement quand et où cette turbulence va se déclencher.

Voici comment cela fonctionne, avec quelques analogies :

1. Le Problème : La "Vague" invisible

Pour prédire la turbulence, les scientifiques étudient les petites perturbations (des "vagues" invisibles) dans l'air.

• Les anciennes méthodes étaient comme regarder une goutte d'eau dans une rivière locale. Elles supposaient que la rivière était droite et calme partout. C'est utile, mais ça rate les interactions complexes entre différentes parties de la rivière.
• HYMOR, lui, regarde toute la rivière d'un coup. Il peut voir comment une vague créée au début de la rivière interagit avec un rocher à l'autre bout, même s'ils sont loin l'un de l'autre. C'est ce qu'on appelle une analyse globale.

2. La Spécificité : Le "Mur de Feu" (L'Onde de Choc)

Quand un objet va très vite, il crée une onde de choc (un mur de pression invisible) devant lui.

• L'erreur classique : La plupart des logiciels traitent ce mur comme un mur de brique épais et flou. C'est pratique pour le calcul, mais cela fausse la physique : les petites vagues qui touchent ce mur flou réagissent mal. C'est comme essayer de jouer du piano sur un coussin épais au lieu d'un clavier dur.
• La solution HYMOR : Il utilise une technique appelée "Shock-fitting". Il traite l'onde de choc comme une lame de rasoir infiniment fine et précise. Il sait exactement comment chaque petite vibration de l'air rebondit ou traverse ce mur. Cela permet d'obtenir une prédiction mathématiquement parfaite de la réponse de l'air.

3. Les Trois Super-Pouvoirs de HYMOR

HYMOR ne fait pas qu'une seule chose, il a trois modes de fonctionnement pour comprendre la transition vers la turbulence :

• Mode "Révélateur" (Analyse Modale) : Il cherche la "note de musique" qui va faire chanter l'avion. Il identifie la vibration naturelle qui va grandir et devenir dangereuse.
• Mode "Effet Papillon" (Croissance Transitoire) : Parfois, même si une vibration est stable, une petite poussée soudaine (comme un coup de vent) peut faire exploser l'énergie de l'air temporairement. HYMOR calcule quelle est la poussée la plus efficace pour créer ce chaos.
• Mode "Porte d'Entrée" (Réceptivité) : C'est le lien entre l'extérieur et l'intérieur. L'air libre contient des turbulences (vent, poussière). HYMOR calcule comment ces turbulences extérieures traversent le mur de choc (l'onde de choc) pour entrer dans la couche d'air collée à l'avion et y créer des problèmes. C'est comme voir comment une goutte de pluie extérieure traverse la vitre d'une voiture pour mouiller l'intérieur.

4. Pourquoi c'est génial ?

• C'est gratuit et ouvert : Comme le code source est public (sous licence MIT), n'importe qui peut l'utiliser, le modifier ou l'améliorer. C'est comme un Lego numérique pour les ingénieurs.
• C'est rapide : Le logiciel est écrit dans des langages modernes (MATLAB et Julia) et utilise la puissance des cartes graphiques (GPU) de votre ordinateur, comme dans les jeux vidéo, pour faire des calculs énormes en quelques secondes.
• C'est précis pour les environnements extrêmes : Il gère la chimie complexe de l'air à haute température (comme sur Mars ou lors de la rentrée atmosphérique), là où les logiciels classiques échouent souvent.

En résumé

HYMOR est une boîte à outils numérique qui permet aux ingénieurs de simuler le comportement de l'air autour de véhicules hypersoniques avec une précision chirurgicale. Au lieu de deviner où l'avion va chauffer, ils peuvent maintenant "voir" les vagues invisibles de l'air, comprendre comment elles traversent les murs de choc, et prédire exactement quand la turbulence va frapper.

C'est un pas de géant pour concevoir des avions plus sûrs, plus rapides et capables de voyager sur d'autres planètes comme Mars.

Résumé technique

Titre du papier

HYMOR : Un package open-source pour l'analyse modale, non modale et de réceptivité dans les véhicules hypersoniques à haute enthalpie.

1. Problématique

La transition laminaire-turbulente dans les écoulements hypersoniques à haute enthalpie reste l'un des problèmes majeurs non résolus en aérodynamique. La prédiction précise de l'endroit et du moment où cette transition se produit est cruciale pour dimensionner les systèmes de protection thermique et établir les marges de conception des véhicules.

Les défis spécifiques incluent :

• Complexité physique : Les effets de gaz réel (dissociation, ionisation, non-équilibre thermo-chimique) deviennent significatifs à haute enthalpie, modifiant les propriétés de l'écoulement et les caractéristiques de stabilité.
• Limitations des méthodes actuelles : Les méthodes traditionnelles (LST, PSE) sont souvent locales ou faiblement non parallèles et ne peuvent pas capturer les interactions entre des mécanismes physiques séparés spatialement. De plus, les méthodes de capture de choc (shock-capturing) introduisent une épaisseur artificielle du choc qui corrompt la réponse linéaire aux perturbations infinitésimales, faussant ainsi les analyses de réceptivité.
• Manque d'outils intégrés : Il n'existait pas auparavant de package open-source unique intégrant l'analyse de stabilité globale, la croissance transitoire et la réceptivité à l'écoulement libre pour les écoulements à haute enthalpie.

2. Méthodologie

HYMOR est une boîte à outils de calcul open-source (disponible en MATLAB et Julia) conçue pour l'analyse de stabilité linéaire globale.

A. Formulation Physique et Numérique

• Équations gouvernantes : Résolution des équations de Navier-Stokes compressibles non dimensionnalisées.
• Modèles thermo-chimiques : Le code intègre plusieurs modèles pour traiter les effets de gaz réel, allant du gaz parfait calorique (CPG) aux modèles de non-équilibre chimique (NonEq-RTVE). Les propriétés thermodynamiques effectives ($\gamma^*$, $c_v^*$) sont pré-calculées via des ajustements par fonctions de base radiale (RBF) pour accélérer les calculs (facteur $\times 1000$ par rapport à la résolution directe d'équilibre).
• Approche Shock-Fitting (Ajustement de choc) : C'est une contribution centrale. Au lieu de capturer le choc numériquement (ce qui crée une épaisseur artificielle), HYMOR traite le choc d'onde comme une discontinuité nette. Les conditions de saut de Rankine-Hugoniot sont imposées comme conditions aux limites sur une spline cubique qui évolue indépendamment de la grille. Cela élimine les artefacts numériques (comme le phénomène de carbuncle) et assure une réponse exacte à l'analyse d'interaction linéaire (LIA).
• Discrétisation : Schéma aux volumes finis d'ordre 2 sur des grilles curvilignes structurées. Intégration temporelle par Runge-Kutta d'ordre 4 (explicite) ou Euler implicite.

B. Analyses de Stabilité

Le package propose trois approches complémentaires basées sur une linéarisation automatique des opérateurs discrets autour d'un écoulement de base :

1. Analyse Modale Globale : Résolution du problème aux valeurs propres ($\lambda \mathbf{q}' = \mathbf{L}\mathbf{q}'$) pour identifier les modes instables. Pour contourner la densité mémoire des matrices, une transformation spectrale préservant la parcimonie (approximation exponentielle $\exp(\mathbf{L}\tau)$) est utilisée avec la méthode d'Arnoldi implicitement redémarrée.
2. Analyse Non Modale (Croissance Transitoire) : Calcul de la croissance maximale de l'énergie (norme de Chu) pour identifier les perturbations initiales optimales qui peuvent déclencher la transition même si l'écoulement est modalement stable.
3. Analyse de Réceptivité à l'Écoulement Libre : Modélisation de l'interaction entre les perturbations de l'écoulement libre et l'écoulement post-choc. Le choc agit comme un couplage convertissant les modes (ex: onde d'entropie en ondes acoustiques). L'analyse traite le système comme un problème entrée-sortie pour déterminer la réponse optimale du véhicule aux perturbations externes.

C. Performance

• Le code est optimisé pour les GPU (CUDA), permettant une accélération significative des solveurs de stabilité (calculs de produits matrice-vecteur dans les boucles internes).
• L'implémentation en Julia et MATLAB offre une flexibilité d'utilisation tout en maintenant des performances compétitives grâce à l'utilisation de bibliothèques BLAS/LAPACK vectorisées.

3. Résultats et Validation

Les auteurs ont validé HYMOR sur une série de cas tests comparatifs :

• Validation du Shock-Fitting : Comparaison avec des solutions de référence (Carpenter & Casper, Sinclair & Cui) pour des écoulements autour de cylindres et de sphères à différents nombres de Mach. Accord excellent pour la densité et la pression.
• Modèles de gaz réel : Validation des modèles d'équilibre et de non-équilibre chimique contre le code Eilmer (chimie à taux fini). HYMOR reproduit correctement les pics de température translationnelle-rotationnelle et les délais de relaxation chimique.
• Stabilité linéaire :
  • Écoulement de canal incompressible (convergence des valeurs propres avec les références Canuto et Nektar++).
  • Cavité entraînée (incompressible et compressible) : Détermination précise du nombre de Reynolds critique pour la bifurcation de Hopf.
• Interaction Choc-Perturbation : Validation de la conversion de modes à travers un choc 1D contre la solution analytique de l'Analyse d'Interaction Linéaire (LIA). Le code capture correctement l'amplification des perturbations et la génération de modes acoustiques à partir d'ondes d'entropie.
• Cas d'application : Application à un cône émoussé hypersonique ($M_\infty = 12$, $Re_\infty = 10^5$) avec gaz réel.
  • L'analyse modale montre que l'écoulement est stable (modes amortis).
  • L'analyse transitoire révèle une amplification d'énergie d'un facteur ~40, dominée par l'énergie entropique, via le mécanisme d'Orr.
  • L'analyse de réceptivité montre une amplification massive (~30 000) due à la combinaison de l'amplification à travers le choc (échelle $M_\infty^2$) et de la couche de cisaillement.

4. Contributions Clés

1. Intégration Unique : Premier package open-source combinant analyse modale, non modale et de réceptivité pour les écoulements hypersoniques à haute enthalpie.
2. Précision du Choc : L'utilisation systématique de la méthode shock-fitting garantit que les interactions choc-perturbation sont traitées avec une précision théorique, éliminant les erreurs liées aux méthodes de capture de choc.
3. Modélisation Avancée : Prise en compte complète des effets de gaz réel (équilibre et non-équilibre) via des modèles thermo-chimiques efficaces.
4. Efficacité Computationnelle : Utilisation de transformations spectrales pour gérer les grandes matrices creuses et exploitation des GPU pour accélérer les calculs de stabilité, rendant possible l'analyse de géométries complexes.

5. Signification

HYMOR comble un vide critique dans la communauté de l'aérodynamique hypersonique. En fournissant un outil transparent, reproductible et capable de capturer les mécanismes physiques complexes (interactions spatiales, effets de gaz réel, réceptivité), il permet aux chercheurs de mieux comprendre les mécanismes de transition vers la turbulence. Cela est essentiel pour le développement de véhicules hypersoniques plus sûrs et plus performants, en particulier pour les missions de rentrée atmosphérique où la charge thermique est critique. Le code est disponible sous licence MIT, favorisant le développement collaboratif et la reproductibilité scientifique.