Distributed Roughness-Induced Transition on a Blunt Body at Mach 6: a Numerical Investigation

Cette étude présente la première simulation numérique directe de la rugosité distribuée sur un cylindre émoussé à Mach 6, révélant que l'arrangement de la rugosité dicte le mécanisme de transition en favorisant soit des modes de stries sinueuses, soit des ondes de T-S bidimensionnelles, ces dernières étant uniquement maintenues par une boucle de rétroaction acoustique interne qui élimine le besoin d'un forçage externe.

L'essentiel

Imaginez un vaisseau spatial fonçant à travers l'atmosphère à six fois la vitesse du son. Pour survivre à la chaleur intense, sa surface est recouverte d'un matériau spécial qui se consume lentement (ablation) pour protéger le vaisseau. Cependant, à mesure que ce matériau brûle, il ne laisse pas une surface parfaitement lisse ; il laisse derrière lui une texture bosselée et rugueuse, semblable à du papier de verre.

Ce document est une simulation informatique à haute vitesse qui pose une question simple mais cruciale : comment ces minuscules bosses sur la surface transforment-elles un écoulement d'air lisse et ordonné en un écoulement d'air turbulent et chaotique ?

Voici le détail de leurs découvertes en utilisant des analogies de la vie quotidienne :

1. La configuration : Le cylindre « papier de verre »

Les chercheurs ont construit un modèle numérique d'un cylindre émoussé (comme le nez d'une fusée) volant à Mach 6. Au lieu d'une surface lisse, ils l'ont recouverte de minuscules « bosses » artificielles (rugosité) pour imiter la texture sablonneuse laissée par le matériau en combustion.

Ils ont testé trois façons différentes de disposer ces bosses :

• Alignées : Comme des soldats debout en rangs et colonnes parfaits.
• Décalées : Comme un mur de briques, où les bosses d'une rangée sont décalées par rapport à la rangée située derrière elle.
• Aléatoires : Comme des cailloux éparpillés sur un trottoir sans aucun motif.

2. L'ancienne théorie vs la nouvelle découverte

Pendant longtemps, les scientifiques ont pensé que la transition vers la turbulence était causée par une « accumulation lente » d'énergie, semblable à une balançoire qui prend de la hauteur si on la pousse juste assez bien au fil du temps. C'est ce qu'on appelle la « croissance transitoire ».

La découverte du document :
La simulation a montré que cette théorie de l'« accumulation lente » n'explique pas vraiment ce qui se passe ici. Les bosses à la surface ne se contentent pas d'amplifier lentement l'énergie ; elles agissent comme des destabilisateurs. Elles prennent l'écoulement de l'air et le rendent immédiatement instable, le transformant en un type d'onde qui croît très rapidement.

Voyez cela de cette façon : l'ancienne théorie pensait que les bosses poussaient doucement un domino pour le faire tomber. La nouvelle découverte montre que les bosses sont en fait en train de donner un coup de pied au domino, le faisant s'écraser contre le suivant immédiatement.

3. Les deux types d'« ondes »

Selon la façon dont les bosses étaient disposées, l'écoulement de l'air réagissait de deux manières différentes :

• Le « Serpent » (Mode sinusoïdal) : Lorsque les bosses étaient alignées (rangées parfaites), l'écoulement de l'air commençait à onduler de gauche à droite comme un serpent. C'est un vacillement très spécifique et organisé.
• L'« Onde plate » (Ondes de Tollmien-Schlichting ou T-S) : Lorsque les bosses étaient décalées ou aléatoires, l'écoulement de l'air commençait à onduler de haut en bas selon un motif d'onde plane en 2D. C'est un type classique d'onde que l'on trouve habituellement dans des flux d'air beaucoup plus lents et à basse vitesse, ce qui est surprenant à trouver dans cet environnement à haute vitesse.

L'idée clé : La disposition des bosses dictait quelle « danse » l'air allait exécuter. La danse du « serpent » se produisait avec les bosses alignées, tandis que la danse de l'« onde plate » se produisait avec les autres.

4. Le final en « Épingle à cheveux »

Une fois que ces ondes sont devenues assez fortes, elles ont déclenché l'étape finale du crash. Les « stries » d'air constantes créées par les bosses (qui sont comme de longs rubans invisibles d'air lent) se sont soudainement tordues et ont cassé en vortex en épingle à cheveux.

Imaginez un élastique qui est tendu. Soudain, il se tord et forme une boucle qui ressemble à une épingle à cheveux. Ces boucles sont la naissance de la turbulence. Une fois que ces épingles à cheveux se forment, l'air lisse se brise complètement en un chaos total, et la chaleur sur la surface du vaisseau spatial grimpe de manière spectaculaire.

5. La surprise de la « Chambre d'écho »

L'une des découvertes les plus fascinantes fut la manière dont la turbulence a commencé pour la première fois dans les cas décalés et aléatoires.

D'habitude, les scientifiques pensent qu'il faut une « poussée » externe (comme une rafale de vent ou une vibration) pour démarrer ces ondes. Mais la simulation a montré un processus auto-entretenu :

1. La turbulence commence en un point sur le cylindre.
2. Comme l'air derrière l'onde de choc se déplace plus lentement que le son (subsonique), le bruit de cette turbulence remonte vers l'amont, comme un écho.
3. Cet « écho » frappe la partie lisse de la surface située devant la turbulence et excite l'air à cet endroit, provoquant le début d'une nouvelle turbulence.
4. Cela crée une boucle de rétroaction : la turbulence produit du bruit, le bruit remonte vers l'amont, et le bruit crée de la nouvelle turbulence.

C'est comme un micro qui capte sa propre sortie de haut-parleur et crée un cycle de larsen strident, mais dans ce cas, le « sifflement » est l'air qui devient turbulent.

Résumé

Ce document a utilisé un super-ordinateur pour observer l'écoulement de l'air sur un cylindre bosselé à haute vitesse. Ils ont découvert que :

• Le motif des bosses décide exactement de la manière dont l'air devient turbulent.
• L'ancienne idée d'une « accumulation lente d'énergie » n'est pas le principal coupable ; au contraire, les bosses déstabilisent directement des ondes spécifiques.
• Ces ondes croissent jusqu'à se tordre en formes d'« épingle à cheveux », provoquant le passage de l'air au chaos.
• Dans certains cas, la turbulence crée son propre « écho » qui remonte vers l'amont pour relancer le processus, sans avoir besoin d'une aide extérieure.

Cela aide les ingénieurs à comprendre que les minuscules bosses aléatoires laissées par la combustion des boucliers thermiques ne sont pas de simples imperfections mineures ; elles sont les architectes principaux de la manière et du moment où la surface d'un vaisseau spatial devient dangereusement chaude.

Résumé technique

Résumé Technique : Transition induite par la rugosité distribuée sur un corps émoussé à Mach 6

Énoncé du problème
La prédiction de la transition laminaire-turbulente (LTT) sur des géométries émoussées à haute vitesse demeure un défi critique pour la conception technique, particulièrement pour les véhicules d'entrée atmosphérique où l'ablation de surface crée une rugosité distribuée. Bien que des corrélations empiriques existent (ex. Hollis 2017, 2019, 2025), elles manquent d'un fondement physique clair. Les efforts théoriques antérieurs, tels que la corrélation de Reshotko & Tumin (2004), reposent sur des mécanismes de croissance transitoire mais font face à des limitations : ils supposent une mise à l'échelle linéaire des perturbations avec la hauteur de rugosité (valable uniquement pour de faibles valeurs de $k/\delta$), nécessitent l'ajustement de constantes à partir de données expérimentales, et peinent à expliquer la transition dans les cas où les hauteurs de rugosité dépassent l'épaisseur de la couche limite ($k/\delta > 1$). De plus, les investigations numériques antérieures ont souvent simulé des patchs de rugosité plutôt qu'une rugosité entièrement distribuée, ou n'ont pas observé d'amplification d'instabilité modale, laissant les mécanismes spécifiques pilotant la transition dans les écoulements de corps émoussés hypersoniques à rugosité pleinement distribuée peu clairs.

Méthodologie
Cette étude présente la première simulation numérique directe (DNS) d'un cylindre émoussé en flux transversal à Mach 6 avec des éléments de rugosité distribués sur toute la surface. Les simulations utilisent le solveur CHAMPS pour intégrer les équations de Navier-Stokes compressibles (CNSE) pour un gaz parfait. La géométrie est un cylindre de 0,1524 m de diamètre, modélisé d'après les conditions expérimentales de Hollis (2017) (Run 38, Test 6975), avec un nombre de Mach en flux libre de 5,98.

La rugosité est modélisée comme une distribution sinusoïdale avec une amplitude fixe ($A \approx 35 \%$ de l'épaisseur de la couche limite) et une longueur d'onde correspondant à une rugosité de grain de sable nominale. Trois configurations de déphasage de l'ordonnée a été étudiée :

1. Décalée (Staggered) : Les rangées sont déphasées de $180^\circ$.
2. Alignée (Aligned) : Les rangées sont en phase (décalage de $0^\circ$).
3. Aléatoire (Random) : Les rangées présentent des déphasages distribués de manière aléatoire.

Pour analyser le processus de transition, les auteurs ont employé :

• DNS : Pour capturer l'évolution non linéaire complète de l'écoulement.
• Équations de perturbation linéaire (LDE) : Résolues sur les écoulements moyens moyennés issus de la DNS pour isoler les mécanismes d'instabilité par forçage impulsionnel.
• Théorie de la stabilité linéaire (LST) : Appliquée aux profils moyennés sur l'ordonnée pour comparaison.
• Analyse spectrale : Transformées de Fourier rapides (FFT) des champs de vitesse et de pression pour identifier les fréquences et les nombres d'onde dominants.

Résultats clés

• Écoulement moyen et stries : Toutes les configurations génèrent des stries de vitesse longitudinale via un mécanisme de "lift-up" piloté par des paires de vortex contre-rotatifs (CVP). La configuration alignée produit les stries stationnaires les plus fortes en raison d'un déphasage cohérent élevant le fluide à faible flux de masse plus haut dans la couche limite. Les configurations décalée et aléatoire produisent des motifs de stries plus faibles et plus complexes. Contrairement aux prédictions de la théorie de la croissance transitoire optimale, les amplitudes de stries ont présenté une faible croissance relative (moins d'un facteur 2) et sont restées non linéaires dans tout le domaine, remettant en question l'utilité de la croissance transitoire linéaire comme outil prédictif pour ce régime.

• Mécanismes d'instabilité : Le processus de transition a été piloté par des instabilités convectives, mais le type de mode spécifique dépendait du déphasage de la rugosité :

  • Cas aligné : Dominé par un mode de strie sinusoïde fondamental. Les ondes de Tollmien-Schlichting (T-S) étaient faiblement déstabilisées.
  • Cas décalé et aléatoire : Dominés par des ondes T-S 2D (mode varieux). Ces configurations ont présenté une déstabilisation significative des ondes T-S, conduisant à une croissance exponentielle rapide.
  • Rupture (Breakdown) : Dans tous les cas, une fois que les ondes d'instabilité dominantes ont atteint une amplitude suffisante, elles ont déclenché la formation de vortex en épingle de cheveux à partir des stries stationnaires, menant à la rupture et à la turbulence.

• Localisations de transition : La DNS a prédit des localisations de transition ($\Theta_{tr}$) nettement plus en aval (environ $11^\circ$–$17^\circ$) que les corrélations empiriques de pointe (qui prédisaient $6^\circ$–$13^\circ$). Le cas aligné a été le premier à entrer en transition parmi les cas de rugosité en raison du forçage requis pour initier la LTT, tandis que les cas décalé et aléatoire ont été plus tardifs.

• Mécanisme de rétroaction : Une découverte novatrice dans les cas décalé et aléatoire est l'identification d'une boucle de rétroaction. La turbulence générée dans la région subsonique du cylindre diffuse des ondes acoustiques qui se propagent vers l'amont. Ces ondes acoustiques lentes excitent les ondes T-S sur le côté opposé du cylindre (à travers la ligne de stagnation) à leur point neutre. Ce mécanisme, qui ne nécessite aucun forçage externe, soutient le processus de transition.

Signification et affirmations
Les auteurs affirment que ce travail fournit une nouvelle compréhension physique de la transition induite par la rugosité sur les corps émoussés hypersoniques, remettant en question le recours prédominant à la théorie de la croissance transitoire. Les principales contributions sont :

1. Rejet de la croissance transitoire comme moteur primaire : L'étude démontre que pour une rugosité entièrement distribuée avec $k/\delta \approx 0,35$, le développement non linéaire des stries et les faibles rapports d'amplification observés font de la théorie de la croissance transitoire optimale un prédicteur insuffisant du processus de transition réel.
2. Déstabilisation des instabilités modales : L'article établit que la rugosité distribuée peut efficacement déstabiliser les instabilités modales (spécifiquement les ondes T-S), même sur des corps émoussés où l'écoulement sur paroi lisse est typiquement stable par mode propre. La configuration de la rugosité dicte si la transition est pilotée par des modes de strie ou par des ondes T-S.
3. Découverte d'une boucle de rétroaction acoustique : L'identification d'un mécanisme de rétroaction auto-entretenu, où les acoustiques générées par la turbulence excitent les ondes T-S de l'autre côté du corps, offre une explication potentielle à la transition dans les régions subsoniques des écoulements hypersoniques, un phénomène auparavant inexpliqué par les mécanismes linéaires.
4. Sensibilité au déphasage : Les résultats soulignent l'extrême sensibilité de la LTT au déphasage relatif des éléments de rugosité, suggérant que des topologies simples peuvent borner le comportement de modèles de rugosité plus complexes ou aléatoires.

Les auteurs concluent que, bien que le mécanisme de rétroaction acoustique proposé nécessite une preuve plus rigoureuse, l'exposition de ces mécanismes d'instabilité via la DNS haute fidélité suggère un changement fondamental dans la compréhension des perturbations responsables de la LTT sur les corps émoussés hypersoniques.