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

Dit paper introduceert HYMOR, een open-source pakket in MATLAB en Julia dat globale modale, niet-modale en receptiviteitsanalyses mogelijk maakt voor hypersonische stromingen bij hoge enthalpie, waarbij real-gas effecten en een schok-fittingformulering worden gebruikt om interacties tussen ruimtelijk gescheiden fysieke mechanismen nauwkeurig te modelleren.

De kern

🚀 HYMOR: De "Proefballon" voor Hypersonische Ruimteschepen

Stel je voor dat je een ruimteschip ontwerpt dat zo snel vliegt dat het de lucht voor zich in vuur verandert. Dit is een hypersonisch vliegtuig (sneller dan 5x de geluidssnelheid). Het grootste probleem bij zo'n snelheid is niet dat het uit elkaar valt, maar dat de hitte het smelt. Om dit te voorkomen, moeten we precies weten wanneer en waar de stroming rondom het schip van rustig (laminair) naar chaotisch (turbulent) overgaat.

Deze overgang is als het verschil tussen een rustig stromende rivier en een woelige stroomversnelling. Zodra het turbulent wordt, stijgt de wrijving en de hitte enorm.

HYMOR is een nieuwe, gratis softwaretool (een "open-source pakket") die wetenschappers helpt om deze overgang te voorspellen. Het is ontwikkeld door onderzoekers van Caltech en MIT.

Hier is hoe het werkt, vertaald in alledaagse termen:

1. De "Scherpe Mes" vs. De "Wazige Foto" (Shock-fitting)

Wanneer een schip zo snel vliegt, ontstaat er een stootgolf (een schokgolf) ervoor. Dit is een onzichtbare muur van samengeperste lucht.

• De oude manier (Shock-capturing): De meeste computersimulaties behandelen deze schokgolf als een wazige foto. Ze proberen de overgang te "vangen" door hem over een paar pixels te verspreiden. Dit introduceert kleine foutjes, alsof je een scherpe mesrand probeert te tekenen met een kwastje. Voor de precieze fysica van de overgang is dit niet goed genoeg.
• De HYMOR manier (Shock-fitting): HYMOR behandelt de schokgolf als een scherp mes. Het is een perfecte, dunne lijn waar de luchtplaatselementen direct van veranderen. Dit zorgt ervoor dat de interactie tussen de lucht en de schokgolf 100% accuraat wordt berekend, zonder die "wazige" foutjes.

2. Drie Manieren om te Kijken (Modal, Non-modal, Receptivity)

HYMOR kijkt naar het probleem op drie verschillende manieren, alsof je een auto bekijkt om te zien waarom hij trilt:

• A. De "Stabiele Rijdende Auto" (Modal Analysis):
  Kijkt naar de natuurlijke trillingen van het systeem. Zou het schip trillen als je het een klein duwtje gaf? HYMOR zoekt de "slechtste" trillingen die kunnen ontstaan. Als deze trillingen groter worden, is het schip onstabiel.

  • Vergelijking: Het is alsof je een gitaarsnaar plukt en kijkt of de toon zuiver blijft of uitloopt.

• B. De "Korte Schok" (Non-modal / Transient Growth):
  Soms is het systeem stabiel op lange termijn, maar kan een korte, krachtige duw (een tijdelijke verstoring) enorme schade aanrichten voordat het weer rustig wordt.

  • Vergelijking: Stel je een toren van kaarten voor. Hij staat stabiel (geen trillingen), maar als je er snel langs loopt (een tijdelijke verstoring), kan hij instorten. HYMOR berekent hoe groot die "duw" mag zijn voordat de toren omvalt.

• C. De "Deur naar de Storm" (Freestream Receptivity):
  Dit is misschien wel het belangrijkste deel. De lucht voor het schip zit vol met kleine onrustjes (turbulentie, geluid). Hoe komen die onrustjes het schip binnen? Ze moeten eerst door de stootgolf heen.

  • Vergelijking: De stootgolf is als een deur. Als je een rimpeling in de lucht hebt (een steen die in een meer valt), verandert de stootgolf die rimpeling volledig terwijl hij erdoorheen gaat. HYMOR berekent precies hoe die rimpeling eruitziet na het passeren van de deur, zodat we weten of hij het schip gaat verstoren.

3. De "Magische Lucht" (Hoge Entalpie & Real-gas)

Bij zulke hoge snelheden is de lucht niet meer gewoon lucht. De moleculen beginnen te trillen, te breken (dissociëren) en zelfs elektrisch geladen te worden (ioniseren). Het gedraagt zich als een "magische vloeistof" met andere eigenschappen.

• HYMOR heeft ingebouwde "vertalers" (thermochimische modellen) die begrijpen hoe deze hete, rare lucht zich gedraagt, of je nu op Aarde vliegt of op Mars (waar de lucht anders is).

4. Waarom is dit zo speciaal?

Vroeger moesten wetenschappers verschillende softwareprogramma's gebruiken voor elk deel van het probleem (één voor de basisstroom, één voor de trillingen, één voor de schokgolf).

• HYMOR is de "Zwitsers zakmes": Het doet alles in één pakket. Het berekent de basisstroom, lineariseert het automatisch (maakt het simpel voor de berekening), en voert alle drie de analyses uit.
• Het is gratis en open: Iedereen kan het gebruiken, bekijken en verbeteren.
• Het is snel: Het maakt slim gebruik van moderne grafische kaarten (GPUs) om de zware rekenkracht te leveren, net zoals videospelletjes dat doen.

Conclusie

HYMOR is als een superkrachtige, digitale windtunnel die niet alleen kijkt naar hoe het schip vliegt, maar ook precies voorspelt waar de hitte en trillingen gaan ontstaan. Door de schokgolf als een perfect scherp object te behandelen en rekening te houden met de extreme hitte, helpt deze tool ingenieurs om veiligere en efficiëntere hypersonische voertuigen te bouwen. Het is een grote stap voorwaarts in het begrijpen van de meest extreme vluchten die de mensheid ooit zal ondernemen.

Technische samenvatting

1. Probleemstelling

De overgang van laminaire naar turbulente stroming in hypersonische flows blijft een van de meest uitdagende open problemen in de snelheidsaerodynamica. Het voorspellen van wanneer en waar deze overgang optreedt is cruciaal voor het dimensioneren van thermische beschermingssystemen en het bepalen van veiligheidsmarges voor voertuigontwerpen.

Bestaande methoden voor stabiliteitsanalyse hebben beperkingen:

• Lokale methoden (LST, PSE): Deze veronderstellen vaak evenwijdige stroming of zwakke niet-parallelle stroming en kunnen interacties tussen ruimtelijk gescheiden fysische mechanismen niet vastleggen.
• Shock-capturing methoden: De meest gebruikte numerieke aanpak voor schokgolven introduceert een kunstmatige schokdikte. Dit vervormt de respons van infinitesimale verstoringen op de schok (Linear Interaction Analysis of LIA), wat essentieel is voor receptiviteitsanalyses.
• Gebrek aan geïntegreerde open-source tools: Bestaande softwarepakketten (zoals LASTRAC, STABL) zijn vaak niet open-source, beperkt tot lokale methoden, of bieden geen geïntegreerde framework voor globale modale, niet-modale en receptiviteitsanalyses in high-enthalpy regimes met reële-gas effecten.

2. Methodologie

Het paper introduceert HYMOR (HYpersonic MOdal/non-modal, and Receptivity), een open-source rekenframework (geschreven in MATLAB en Julia) dat de volgende kernmethodes combineert:

• Governing Equations: De solver lost de compressibele Navier-Stokes-vergelijkingen op, gekoppeld aan een ideaal-gas toestandsvergelijking met effectieve thermodynamische eigenschappen. Het ondersteunt diverse thermochemische modellen voor high-enthalpy flows, variërend van calorisch perfecte gassen tot chemisch niet-evenwicht (NonEq-RTVE), inclusief vibratie- en elektronische excitatie.
• Shock-fitting aanpak: In plaats van schokken op te vangen (shock-capturing), wordt de boogschock behandeld als een scherpe discontinuïteit met exacte Rankine-Hugoniot-sprongvoorwaarden. De schok wordt gemodelleerd met een kubische spline die onafhankelijk van het rekenrooster beweegt. Dit elimineert numerieke artefacten (zoals het "carbuncle"-fenomeen) en zorgt voor een exacte weergave van de interactie tussen verstoringen en de schok, essentieel voor LIA.
• Lineaire Stabiliteitsanalyse:
  • Globale Modale Analyse: Oplossen van het eigenwaardeprobleem voor het volledige ruimtelijke domein. Om de hoge rekenkosten van grote matrices te omzeilen, wordt een spectrale transformatie (exponentiële mapping) gebruikt in plaats van de traditionele "shift-and-invert" methode, wat geheugenbesparend is en GPU-versnelling mogelijk maakt.
  • Globale Transient Growth (Niet-modale): Analyse van de maximale energieversterking van verstoringen over een tijdsinterval, gebruikmakend van de Chu-energynorm (aangepast voor chemisch evenwicht). Dit identificeert mechanismen die niet door modale instabiliteit worden voorspeld.
  • Freestream Receptivity: Analyse van hoe vrijstroomverstoringen (via de schok) de grenslaag excitatie. Het model koppelt de vrijstroomverstoringen (gebaseerd op de "frozen turbulence" hypothese) aan de post-schock respons via een lineair gekoppeld systeem.
• Numerieke Implementatie: Gebruik van een tweede-orde eindige-volume schema op een curvilineair rooster. De code biedt een solver voor niet-lineaire basisstroomberekeningen met automatische linearisatie van discrete operatoren.

3. Belangrijkste Bijdragen

1. Eerste Open-Source Integratie: HYMOR is het eerste open-source pakket dat modale, niet-modale en receptiviteitsanalyses voor high-enthalpy hypersonische flows in één framework combineert.
2. Exacte Shock-Behandeling: De implementatie van shock-fitting zorgt voor een nauwkeurige lineaire interactieanalyse (LIA) zonder de vervorming die typisch is bij shock-capturing methoden.
3. Efficiënte Rekenstrategie: Door het gebruik van spectrale transformaties en GPU-versnelling (via CUDA) kunnen grote roosters worden geanalyseerd zonder dat het geheugenprobleem van de "fill-in" bij matrix-factorisatie (zoals bij shift-and-invert) optreedt.
4. Uitgebreide Thermochemie: Ondersteuning voor complexe real-gas effecten (chemisch evenwicht en niet-evenwicht) die essentieel zijn voor flows op grote hoogtes of bij Mars-missies.
5. Dual Language Support: Beschikbaar in zowel MATLAB als Julia met identieke functionaliteit, wat toegankelijkheid en prestaties optimaliseert.

4. Resultaten en Validatie

De auteurs hebben HYMOR uitgebreid gevalideerd tegen benchmarkcases:

• Shock-fitting Validatie: Vergelijkingen met analytische oplossingen en bestaande data voor flows over cilinders en bollen tonen uitstekende overeenkomst in dichtheid, druk en schokstand.
• Real-Gas Modellen: Validatie tegen de Eilmer-code voor chemisch niet-evenwicht flows (bijv. over een wig) toont dat HYMOR de juiste temperatuurspieken en relaxatiegedrag nabij de schok reproduceert.
• Stabiliteitsvalidatie: Resultaten voor incompressibele kanaalstroming en koppelingsstromingen (lid-driven cavity) komen overeen met referentiewaarden uit de literatuur (Canuto, Nektar++).
• LIA Validatie: De code reproduceert exact de mode-conversie (van entropie naar akoestische golven) over een schok, zoals voorspeld door de lineaire interactietheorie van McKenzie en Westphal.
• Toepassing: Een case-studie op een hypersonische stompe kegel (Mach 12, Re=100.000) toont:
  • Modale stabiliteit (de stroming is modaal stabiel).
  • Significante transient growth (energieversterking factor ~40) door Orr-mechanismen in de grenslaag.
  • Extreme receptiviteit (energieversterking factor ~30.000) wanneer vrijstroomverstoringen de schok passeren, waarbij de kinetische energie de dominante bijdrage levert.

5. Betekenis en Conclusie

HYMOR vult een kritieke leemte in het landschap van hypersonische stabiliteitssoftware. Door de combinatie van globale analyse, shock-fitting en real-gas modellen, biedt het onderzoekers een krachtig hulpmiddel om de fysica van overgang in extreme omstandigheden beter te begrijpen.

De open-source aard (MIT-licentie) en de beschikbaarheid in twee populaire programmeertalen bevorderen reproduceerbaarheid en samenwerking binnen de gemeenschap. De code stelt onderzoekers in staat om complexe interacties te bestuderen die met traditionele lokale methoden ontoegankelijk zijn, wat essentieel is voor het ontwerp van toekomstige hypersonische voertuigen en thermische beschermingssystemen. De prestatie-analyses tonen aan dat de code, ondanks het gebruik van hogere programmeertalen, concurrerende snelheden bereikt door efficiënte vectorisatie en GPU-gebruik.