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

Questo studio dimostra l'efficacia dell'adattamento della griglia basato su metriche per simulazioni aerotermiche ipersoniche in gas reali, evidenziando come tale approccio su mesh non strutturate permetta di gestire geometrie complesse, come i sistemi di controllo di reazione, mantenendo una precisione nel riscaldamento superficiale paragonabile a quella delle simulazioni con griglie strutturate.

L'essenziale

Il Problema: Disegnare una mappa per un viaggio pericoloso

Immagina di dover guidare un'auto attraverso una tempesta di fuoco (il rientro atmosferico di un veicolo spaziale) a velocità incredibili. Per farlo in sicurezza, hai bisogno di una mappa estremamente dettagliata che ti mostri esattamente dove il fuoco è più intenso e dove l'aria diventa turbolenta.

In passato, per creare questa "mappa" (che in termini tecnici si chiama mesh o griglia computazionale), gli ingegneri dovevano disegnare a mano ogni singola strada. Era come se dovessi costruire una casa usando solo mattoni quadrati perfetti. Funzionava bene per le case semplici, ma se la casa aveva torri, camini strani o scale a chiocciola (come i razzi di controllo o le forme complesse di una capsula), costruire la mappa diventava un incubo. Inoltre, se i "mattoni" non erano allineati perfettamente con la direzione del vento, la mappa dava previsioni sbagliate su quanto il veicolo si sarebbe scaldato.

La Soluzione: Un'Intelligenza Artificiale che "Ridisegna" la mappa

Questo documento descrive un nuovo metodo chiamato MIMIC (Metric-Informed Mesh Improvement Capability). Immagina MIMIC non come un disegnatore, ma come un architetto intelligente che ha un super-potere: la capacità di ridisegnare la mappa in tempo reale mentre il viaggio avviene.

Ecco come funziona, passo dopo passo:

1. La Mappa Iniziale (Grossolana): Si parte con una mappa "abbozzata", fatta di pezzi di forma irregolare (tetraedri, come piramidi di gelato) che riempiono lo spazio. È veloce da fare, ma non è precisa.
2. Il Termometro Magico (L'Indicatore): Il computer calcola la temperatura del flusso d'aria. Dove la temperatura cambia bruscamente (come quando si attraversa un muro di fuoco improvviso, chiamato shock wave), il sistema capisce: "Qui serve più dettaglio!".
3. Il Rifacimento (Adattamento): Invece di fermarsi e ridisegnare tutto da zero, il sistema MIMIC prende la sua "mappa magica" e la modifica automaticamente:
   • Dove c'è molto calore o turbolenza, stira e rimpicciolisce i pezzi della mappa per vederli meglio (come usare un microscopio su una macchia d'inchiostro).
   • Dove l'aria è calma, ingrandisce i pezzi per risparmiare tempo di calcolo.
4. Il Risultato: Alla fine, hai una mappa che è "intelligente": è densa e precisa dove serve (sullo scudo termico e dietro la capsula) e leggera dove non serve.

Due Esperimenti Chiave

L'autore ha testato questo metodo su due scenari:

1. La Mezza Sfera (Il Test Semplice)
Ha preso una semplice forma a semisfera e ha visto cosa succede se usi "mattoni" di forme diverse vicino alla superficie.

• Il risultato: Ha scoperto che se usi dei "mattoni" piatti e allungati (esagoni o prismi) vicino alla pelle del veicolo, e poi lasci che il sistema adatti il resto dello spazio, ottieni previsioni di calore molto più precise e lisce. È come se avessi un pavimento di piastrelle perfette vicino al camino, e poi un tappeto che si adatta al resto della stanza.

2. La Capsula per Marte (Il Test Complesso)
Qui la situazione si fa seria. Hanno simulato un rientro su Marte con una capsula che ha 8 piccoli razzi di controllo (RCS) sulla schiena.

• Il problema vecchio: Con i metodi tradizionali (che usano mattoni quadrati rigidi), modellare quei 8 piccoli razzi sulla schiena era quasi impossibile senza semplificare troppo la geometria. Sarebbe come cercare di disegnare i dettagli di un'auto di Formula 1 usando solo cubi di Lego grandi.
• Il successo del nuovo metodo: Il sistema MIMIC ha gestito la capsula con i suoi 8 razzi senza problemi. Ha creato una mesh che si adattava perfettamente a ogni curva e a ogni ugello.
• La sorpresa: Il sistema ha scoperto dettagli nascosti nella "scia" dietro la capsula. Ha visto come l'aria si separa e crea vortici complessi vicino ai razzi, qualcosa che i vecchi metodi rigidi non potevano vedere perché erano troppo "rigidi" per seguire quei movimenti caotici.

Perché è Importante?

Prima, per studiare questi fenomeni, gli ingegneri dovevano:

1. Passare settimane a disegnare a mano la mappa.
2. Semplificare la geometria (togliendo dettagli importanti come i razzi di controllo) per farci stare la griglia.
3. Rischiare di sbagliare i calcoli sul calore perché la griglia non era allineata col vento.

Ora, con questo metodo:

• Automazione: Il computer fa il lavoro sporco di ridisegnare la mappa.
• Complessità: Si possono studiare veicoli con forme assurde e dettagli piccoli (come i razzi) senza impazzire.
• Precisione: Si ottengono previsioni sul calore quasi perfette, paragonabili ai metodi vecchi ma con molta più flessibilità.

In Sintesi

Immagina che prima dovessi scolpire una statua di ghiaccio usando solo un martello e uno scalpello rigido: difficile e lento. Ora, hai un 3D printer intelligente che, mentre la statua viene scolpita, aggiunge o toglie materiale esattamente dove serve per catturare ogni singola goccia d'acqua che scorre sulla superficie.

Questo documento ci dice che l'era della "mappa rigida" sta finendo. Ora possiamo simulare il rientro atmosferico su Marte (e oltre) con una precisione e una velocità che prima erano impensabili, permettendo agli ingegneri di progettare veicoli spaziali più sicuri e complessi.

Sommario tecnico

1. Problema e Contesto

Le simulazioni di flussi ipersonici per veicoli di rientro atmosferico richiedono una previsione accurata del riscaldamento aerotermico sulla superficie. Gli attuali codici CFD (Computational Fluid Dynamics) all'avanguardia, come DPLR e US3D, utilizzano tipicamente mesh strutturate a blocchi composte da elementi esagonali (hexahedra) allineati con l'onda d'urto frontale per garantire precisione. Tuttavia, questo approccio presenta limitazioni significative:

• Complessità Geometrica: La generazione di mesh strutturate diventa estremamente laboriosa e spesso impossibile per geometrie complesse che includono dettagli come i sistemi di controllo di reazione (RCS), gradini o superfici irregolari.
• Limitazioni delle Mesh Non Strutturate: Le mesh non strutturate (basate su tetraedri) sono flessibili ma tendono a produrre previsioni di riscaldamento superficiale inaccurate a causa della scarsa allineamento con l'onda d'urto, che genera entropia numerica e "inquinamento" del flusso a valle.
• Dinamica della Scia: La previsione del riscaldamento sulla parte posteriore (backshell) è incerta a causa della natura caotica della scia, che richiede una risoluzione adattiva per catturare interazioni complesse (onde d'urto, strati di taglio).

L'obiettivo del lavoro è dimostrare l'applicabilità e l'efficacia dell'adattamento della mesh basato su metriche (metric-based mesh adaptation) per problemi di gas reale in regime ipersonico, mantenendo accuratezza nel riscaldamento superficiale pur gestendo geometrie complesse.

2. Metodologia

Il lavoro combina due strumenti principali:

1. US3D: Un solver di flussi comprimibili avanzato basato sul metodo dei volumi finiti, capace di gestire miscele di gas reagenti (gas reale), turbolenza (RANS/IDDES) e flussi multi-specie.
2. MIMIC (Metric-Informed Mesh Improvement Capability): Una nuova capacità sviluppata per l'adattamento anisotropo della mesh.

Flusso di lavoro MIMIC:

• Indicatore di Errore: MIMIC calcola un campo tensoriale metrico basato sull'Hessiana (gradienti di secondo ordine) della soluzione della temperatura. Questo guida dove la mesh deve essere raffinata o coarsened.
• Adattamento: Utilizza l'API open-source ParMMG per adattare esclusivamente gli elementi tetraedrici nel volume, allineandoli alle metriche calcolate.
• Struttura Ibrida: Un contributo chiave di questo studio è la capacità di gestire una mesh di strato limite ibrida e fissa composta da prismi, esagoni (hexahedra) e piramidi vicino alla parete, mentre il resto del volume è adattato dinamicamente.
• Iterazione: Il processo è iterativo. Si parte da una mesh di base grezza, si risolve il flusso (inizialmente con schemi di primo ordine per evitare effetti "carbuncle" sull'onda d'urto), si calcola la metrica, si adatta la mesh e si ripete il processo (fino a 5 iterazioni nel caso di studio).

3. Contributi Chiave

• Estensione a Gas Reale: Applicazione della capacità MIMIC a flussi ipersonici con gas reagenti (miscele di $CO_2$-$N_2$), estendendo lavori precedenti limitati a gas perfetti.
• Supporto per Esagoni nello Strato Limite: Estensione della capacità di MIMIC per supportare superfici tessellate con elementi quadrangolari estrusi in strati di esagoni, oltre ai tradizionali triangoli estrusi in prismi.
• Gestione di Geometrie Complesse: Dimostrazione della capacità di simulare veicoli di rientro con dettagli geometrici complessi (come le bocchette RCS) che sono tipicamente fuori portata per approcci a blocchi strutturati.
• Confronto Prismi vs. Esagoni: Analisi sistematica dell'impatto della topologia dello strato limite (prismi vs. esagoni) sulla previsione del flusso di calore.

4. Risultati

Lo studio presenta due casi di test principali:

A. Flusso Supersonico su un Emisfero (Gas Perfetto)

• Confronto Topologie: Sono stati confrontati mesh con strato limite in prismi ($T_p$) e in esagoni ($T_h$), con diverse spessori (80 e 110 strati).
• Risultati:
  • Le mesh non adattate mostrano distribuzioni di calore distorte e rumorose.
  • Dopo 5 iterazioni di adattamento, la qualità della soluzione migliora drasticamente.
  • Sensibilità allo spessore: La risoluzione dello strato limite è critica; 80 strati non sono sufficienti per una previsione accurata, mentre 110 strati lo sono.
  • Prismi vs. Esagoni: Le mesh con strato limite in esagoni producono distribuzioni di flusso di calore più lisce e simmetriche rispetto a quelle in prismi. Con 110 strati di esagoni e adattamento, i risultati di US3D coincidono con i dati di riferimento DPLR (entro un margine di errore del ±7.5%).

B. Flusso Ipersonico su un Capsula a 70° con RCS (Gas Reale - Atmosfera Marziana)

• Scenario: Flusso di una miscela $CO_2$-$N_2$ su una capsula di rientro con 8 getti RCS sul retro.
• Sfida: La geometria complessa dei getti RCS e la dinamica della scia.
• Risultati:
  • Adattamento: Dopo 5 iterazioni, la mesh si raffina automaticamente attorno all'onda d'urto frontale, agli strati di taglio e nella scia, catturando scale spaziali più piccole impossibili da risolvere con mesh strutturate a singolo blocco.
  • Riscaldamento Superficiale: Le previsioni del flusso di calore sulla parte anteriore (forebody) diventano simmetriche e prive di "streaks" numerici dopo le iterazioni finali, mostrando un accordo qualitativo con le simulazioni DPLR strutturate.
  • Geometria RCS: A differenza delle simulazioni DPLR che richiedono di approssimare la geometria (es. superficie liscia), l'approccio US3D adattato ha incorporato fedelmente la geometria dei getti RCS, permettendo di studiare l'ambiente locale e le interazioni nella scia.
  • Fenomeni Fisici: L'adattamento ha rivelato dettagli nella scia, come la rottura degli strati di taglio e l'evoluzione delle bolle di separazione, che non erano visibili nella mesh iniziale.

5. Significato e Conclusioni

Il lavoro dimostra che l'adattamento della mesh basato su metriche, combinato con solver non strutturati come US3D, è una strategia fattibile e potente per l'analisi aerotermica ipersonica complessa.

• Flessibilità Geometrica: Permette agli ingegneri di simulare configurazioni con dettagli geometrici complessi (come i sistemi di controllo di reazione) senza la necessità di semplificare la geometria o generare mesh strutturate laboriose.
• Accuratezza: Dimostra che è possibile ottenere livelli di accuratezza nel riscaldamento superficiale paragonabili a quelli delle mesh strutturate a blocchi (DPLR), a condizione che si utilizzi uno strato limite sufficientemente risolto (preferibilmente esagonale) e un adattamento iterativo basato sulla soluzione.
• Futuro: Questo approccio apre la strada a simulazioni monolitiche unificate che possono gestire sia il flusso esterno che le interazioni dei getti RCS, riducendo i costi di generazione della mesh e migliorando la comprensione della fisica della scia.

In sintesi, il documento valida un nuovo flusso di lavoro che supera le limitazioni tradizionali delle mesh strutturate, offrendo un compromesso ottimale tra flessibilità geometrica e accuratezza fisica per i veicoli di rientro di prossima generazione.