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

Il documento presenta HYMOR, un pacchetto open-source in MATLAB e Julia che fornisce un quadro computazionale per l'analisi di stabilità lineare (modale, non-modale e di ricettività) dei fl ipersonici ad alta entalpia, integrando modelli di gas reale, un trattamento preciso dello shock tramite shock-fitting e la capacità di catturare interazioni tra meccanismi fisici spazialmente separati.

L'essenziale

🚀 HYMOR: Il "Cristallo Magico" per vedere l'invisibile nei voli supersonici

Immagina di voler progettare un'astronave che viaggia a velocità incredibili, così veloci da riscaldare l'aria intorno a essa fino a farla diventare un plasma incandescente (questo è il mondo dell'ipersonica ad alto contenuto entalpico).

Il problema principale? L'aria non si comporta come un fluido normale. Diventa chimicamente attiva, si ionizza e crea onde d'urto potenti come martelli. Se l'aria passa da uno stato "calmo" (laminare) a uno "caotico" (turbolento) sulla superficie dell'astronave, il calore aumenta drasticamente e potrebbe fondere il veicolo.

Gli ingegneri devono prevedere quando e dove avverrà questo passaggio. Per farlo, hanno bisogno di capire come le piccole vibrazioni nell'aria (i "disturbi") crescono e diventano tempeste.

Qui entra in gioco HYMOR.

1. Cos'è HYMOR? (Il "Cristallo di Visione")

HYMOR non è un'astronave, ma un software open-source (gratis e modificabile da tutti) scritto in due linguaggi di programmazione molto potenti (MATLAB e Julia). È come un "cristallo di visione" che permette agli scienziati di simulare il comportamento dell'aria intorno a un'astronave con una precisione che i vecchi metodi non avevano.

Prima, gli ingegneri guardavano l'aria "a pezzi": analizzavano una piccola fetta alla volta, come se guardassero un film fotogramma per fotogramma senza vedere la scena completa. HYMOR, invece, guarda tutto il film insieme (analisi globale). Questo permette di vedere come un'onda d'urto lontana possa influenzare l'aria vicino al muso dell'astronave, collegando punti che prima sembravano scollegati.

2. Il problema del "Fumetto" (L'Onda d'Urto)

Quando un'astronave viaggia più veloce del suono, crea un'onda d'urto (un muro di pressione invisibile).

• I vecchi metodi (Shock-capturing): Immagina di disegnare un'onda d'urto su un foglio a quadretti. Poiché il foglio è fatto di quadratini, l'onda non può essere una linea netta; deve "saltare" da un quadratino all'altro. Questo crea un effetto "sfocato" o "a gradini" (come un fumetto mal disegnato). Questo sfocato inganna il computer, facendogli credere che l'aria si comporti in modo sbagliato quando tocca l'onda.
• Il metodo HYMOR (Shock-fitting): HYMOR tratta l'onda d'urto come una linea perfetta e netta, come un rasoio che taglia l'aria in due. Non la "disegna" sui quadratini, ma la fa scorrere liberamente sopra di essi. È come se avessimo un righello magico che segue l'onda perfettamente, eliminando ogni errore di "sfocatura". Questo è fondamentale perché l'interazione tra l'aria e l'onda d'urto è delicatissima: se sbagli anche di poco qui, tutto il calcolo del calore futuro è sbagliato.

3. Le tre magie di HYMOR

Il pacchetto fa tre cose principali, che possiamo paragonare a tre diversi modi di studiare una tempesta:

• Analisi Modale (Il "Terremoto"): Cerca di capire se l'aria è intrinsecamente instabile. È come chiedere: "Se lascio una goccia d'acqua in un secchio, l'acqua inizierà a tremare da sola?" Se la risposta è sì, il sistema è instabile e diventerà turbolento.
• Crescita Non-Modale (L'"Effetto Farfalla"): A volte l'aria è stabile, ma se la colpisci nel modo esatto al momento esatto, può esplodere in turbolenza. È come spingere un'altalena: se spingi nel momento sbagliato, non succede nulla; se spingi nel momento giusto (anche con poca forza), l'altalena va altissima. HYMOR cerca il "momento perfetto" per spingere l'aria e vedere quanto può crescere.
• Ricezione (L'"Ascolto"): Immagina che l'aria esterna (il "freestream") sia piena di piccoli rumori (turbolenza, polvere, vibrazioni). Come fa l'astronave a "ascoltare" questi rumori e trasformarli in vibrazioni sulla sua pelle? HYMOR simula come i disturbi esterni attraversano l'onda d'urto e vengono amplificati, diventando il seme della turbolenza.

4. Perché è importante? (La "Cucina" dell'Aria)

In questi voli, l'aria non è solo aria: è una zuppa chimica. Le molecole si rompono, si ricombinano e cambiano forma (dissociazione e ionizzazione).
HYMOR ha una "cucina" interna molto sofisticata che può gestire diverse ricette chimiche:

• Aria normale: Come quella che respiriamo.
• Aria "congelata": Dove le reazioni chimiche sono troppo lente per accadere.
• Aria "in equilibrio": Dove le reazioni chimiche avvengono istantaneamente.
• Aria "non in equilibrio": Dove le reazioni stanno accadendo ma non sono finite (come quando accendi un fiammifero e guardi la fiamma che si espande).

Questo permette di simulare voli su Marte (dove l'aria è quasi tutta anidride carbonica) o sulla Terra a quote altissime, con una precisione che prima richiedeva supercomputer costosi e chiusi.

5. La velocità e la potenza

Il software è stato scritto per essere veloce. Sfrutta le schede grafiche (GPU) dei computer moderni, proprio come fanno i videogiochi per renderizzare immagini 3D, per calcolare milioni di punti in pochi secondi. È stato testato contro casi reali e ha dimostrato di funzionare perfettamente, confermando che le sue previsioni sono affidabili.

In sintesi

HYMOR è come un laboratorio virtuale gratuito per ingegneri aerospaziali.

• Non usa "fotogrammi sfocati" ma linee perfette per le onde d'urto.
• Guarda l'intero sistema, non solo i pezzi.
• Capisce la chimica complessa dell'aria rovente.
• Aiuta a progettare scudi termici migliori, rendendo i viaggi spaziali più sicuri ed efficienti.

È uno strumento che democratizza la scienza: invece di essere un segreto custodito da poche grandi aziende, è ora disponibile per chiunque voglia studiare come l'aria si comporta quando viene spinta al limite.

Sommario tecnico

1. Il Problema

La transizione da flusso laminare a turbolento nei flussi ipersonici ad alta entalpia rimane uno dei problemi aperti più critici nell'aerodinamica ad alta velocità. La capacità di prevedere quando e dove avviene questa transizione è fondamentale per il dimensionamento dei sistemi di protezione termica e per stabilire i margini di progettazione dei veicoli.
Le sfide principali includono:

• Complessità fisica: Gli effetti dei gas reali (dissociazione chimica, ionizzazione, ecc.) diventano significativi ad alte entalpie, alterando le proprietà del flusso e la posizione dello stallo dell'onda d'urto.
• Limitazioni dei metodi tradizionali: I metodi di stabilità locale (LST) e le equazioni di stabilità parabolizzate (PSE) non riescono a catturare le interazioni tra meccanismi fisici spazialmente separati, né le instabilità non modali (crescita transitoria).
• Interazione con l'onda d'urto: Le perturbazioni del flusso libero devono attraversare l'onda d'urto bow (d'urto d'attacco) prima di raggiungere lo strato limite. I metodi numerici convenzionali basati sulla "shock-capturing" introducono uno spessore artificiale dell'urto che corrompe la risposta lineare alle perturbazioni infinitesime, rendendo inaccurata l'analisi di ricettività.

2. Metodologia e Implementazione

Il paper presenta HYMOR (HYpersonic MOdal/non-modal, and Receptivity), un framework computazionale open-source (disponibile in MATLAB e Julia) che affronta queste sfide attraverso un approccio globale.

A. Formulazione Fisica e Numerica

• Equazioni di Governo: Risoluzione delle equazioni di Navier-Stokes comprimibili non dimensionali.
• Modelli Termodinamici: Il codice supporta diversi modelli per gli effetti dei gas reali, inclusi:
  • Gas perfetto calorico (CPG).
  • Chimica congelata con equilibrio vibrazionale-rotazionale-traslazionale (Frozen-RTV).
  • Equilibrio chimico e termico (Chemical-RTV, Chemical-RTVE).
  • Non-equilibrio chimico con equilibrio termico (NonEq-RTVE).
    I modelli utilizzano fit di funzioni (radial-basis function) per le proprietà termodinamiche ($\gamma^*$, $c_v^*$) per garantire efficienza computazionale, evitando di risolvere equazioni di equilibrio in ogni cella.
• Approccio Shock-Fitting: Invece di catturare l'urto (shock-capturing), HYMOR utilizza un metodo di shock-fitting. L'onda d'urto è trattata come una discontinuità netta definita da spline cubiche che si muovono indipendentemente dalla griglia computazionale. Questo elimina lo spessore artificiale dell'urto e garantisce che l'interazione tra perturbazioni infinitesime e l'urto riproduca esattamente la risposta prevista dalla Linear Interaction Analysis (LIA).
• Discretizzazione: Schema agli elementi finiti di secondo ordine su griglie strutturate curvilinee. Integrazione temporale con Runge-Kutta esplicito (RK4) o Euler implicito.

B. Analisi di Stabilità Lineare

HYMOR implementa tre approcci globali:

1. Analisi Modale Globale: Risoluzione del problema agli autovalori $\lambda \mathbf{q}' = \mathbf{L}\mathbf{q}'$. Per gestire la grande dimensione della matrice $\mathbf{L}$ (che cresce come $O(N^4)$), viene utilizzata una trasformazione spettrale esponenziale ($\hat{\mathbf{L}} \approx \exp(\mathbf{L}\tau)$) invece del classico shift-and-invert, preservando la sparsità e permettendo l'uso efficiente di GPU.
2. Analisi Non-Modale (Crescita Transitoria): Utilizza la norma di energia di Chu (adattata per la chimica di equilibrio) per quantificare l'amplificazione transitoria delle perturbazioni. Risolve un problema di ottimizzazione generalizzato (quoziente di Rayleigh) per trovare la perturbazione iniziale che massimizza la crescita energetica.
3. Analisi di Ricettività del Flusso Libero: Studia come le perturbazioni del flusso libero, attraversando l'urto, eccitano le modalità nello strato limite. Il sistema è modellato come un sistema input-output dove l'input è la forzante periodica del flusso libero e l'output è la risposta a valle dell'urto.

3. Contributi Chiave

• Integrazione Unica: HYMOR è il primo pacchetto open-source che integra in un'unica framework l'analisi modale, non modale e di ricettività per flussi ipersonici ad alta entalpia.
• Shock-Fitting per la Stabilità Lineare: L'implementazione rigorosa dello shock-fitting garantisce che le interazioni urto-perturbazione siano fisicamente corrette, un requisito essenziale per studi di ricettività che i metodi shock-capturing non possono soddisfare.
• Efficienza Computazionale: L'uso di trasformazioni spettrali esponenziali e l'offloading delle operazioni su GPU (tramite CUDA) permettono di eseguire analisi di stabilità su griglie più grandi rispetto ai metodi tradizionali, superando i colli di bottiglia di memoria.
• Accessibilità: Il codice è rilasciato sotto licenza MIT, disponibile in MATLAB e Julia, con modelli termodinamici validati per atmosfere terrestri e marziane.

4. Risultati e Validazione

Il codice è stato verificato e validato su una serie di casi benchmark:

• Flussi Ideali: Confronto con soluzioni analitiche e dati di riferimento per flussi su cilindri e sfere (es. Carpenter & Casper, Hamilton et al.), mostrando eccellente accordo per densità, pressione e posizione dell'urto.
• Effetti Gas Reale: Validazione dei modelli Chemical-RTVE e NonEq-RTVE contro il codice Eilmer (chimica a tassi finiti). HYMOR riproduce correttamente i picchi di temperatura traslazionale-rotazionale e il rilassamento chimico dietro l'urto.
• Condizioni di Salto Rankine-Hugoniot: Verifica delle condizioni di salto con effetti termodinamici contro la Shock and Detonation Toolbox (SD Toolbox).
• Analisi di Stabilità: Confronto con Nektar++ e letteratura per flussi incompressibili (canale, cavità guidata) e compressibili, dimostrando la convergenza degli autovalori e la corretta identificazione dei modi instabili.
• Interazione Urto-Perturbazione: Verifica dell'accuratezza della conversione modale attraverso un'onda d'urto 1D, confermando che il solver riproduce la soluzione analitica LIA (generazione di onde acustiche ed entropiche da un'onda di entropia incidente).
• Caso Applicativo: Applicazione a un cono ottuso ipersonico ($M_\infty=12$, $Re_\infty=10^5$).
  • L'analisi modale ha mostrato stabilità modale (tutti gli autovalori con parte reale negativa).
  • L'analisi di crescita transitoria ha identificato un'amplificazione energetica di circa 40, guidata da meccanismi di tipo Orr e contributi entropici.
  • L'analisi di ricettività ha mostrato un'amplificazione enorme (circa 30.000), dominata dall'effetto di amplificazione attraverso l'urto ($M_\infty^2$) e dall'interazione con lo shear nello strato limite.

5. Significato

HYMOR rappresenta un avanzamento significativo per la comunità di ricerca sulla stabilità ipersonica. Fornisce uno strumento robusto e accessibile per:

1. Superare i limiti dei metodi locali: Consentendo lo studio di fenomeni globali e interazioni spaziali complesse.
2. Migliorare la predizione della transizione: Integrando la ricettività del flusso libero e gli effetti dei gas reali, offrendo una catena di predizione più completa.
3. Promuovere la riproducibilità: Essendo open-source, facilita lo sviluppo collaborativo e la verifica indipendente dei risultati nella ricerca sulla transizione ipersonica.

In sintesi, il lavoro dimostra che l'approccio shock-fitting combinato con l'analisi di stabilità globale è essenziale per comprendere accuratamente i meccanismi di transizione in regimi ad alta entalpia, dove le approssimazioni tradizionali falliscono.