Introduzione alla dinamica dei fluidi criogenici nei sistemi spaziali

I liquidi criogenici, i gas liquefatti immagazzinati a temperature inferiori a -150 °C, sono l’inondazione vitale dei moderni sistemi di propulsione spaziale e di supporto vitale. L’idrogeno liquido (LH2) e l’ossigeno liquido (LOX) alimentano i principali motori di razzi come il Raptor di SpaceX, il BE-4 di Blue Origin e il sistema di lancio spaziale della NASA.

Gli strumenti di fluidodinamica computazionale (CFD) come COMSOL Multiphysics con il modulo CFD[] offrono un quadro robusto per simulare il comportamento del fluido criogenico. Con l'accoppiamento del flusso di fluido, del trasferimento di calore, della termodinamica e del cambiamento di fase in un unico ambiente, gli ingegneri possono testare i progetti del serbatoio, le linee di alimentazione e le strategie di isolamento prima della costruzione dell'hardware.

Perché i fluidi criogenici richiedono la simulazione speciale

L'idrogeno liquido, ad esempio, ha una densità di circa 70 kg/m3 al punto di ebollizione e una conducibilità termica che cambia bruscamente vicino al punto critico. I regimi di boiling, la transizione e l'ebollizione del film possono coesistere in un singolo serbatoio a causa dell'ingresso di calore localizzato dal motore o dalla radiazione solare.

Queste complessità rendono le correlazioni empiriche non affidabili per il design. La simulazione multifisica deve contemporaneamente risolvere:

  • Trasferimento termico coniugato:[] conduzione attraverso pareti di serbatoio, convezione nel criogenico, e radiazione da fonti esterne.
  • Dinamici di cambiamento di processo:[ evaporazione, condensazione e effetti di gas reale vicino alla linea di saturazione.
  • Registri di flusso bifase:[] scorrere bubbly, stratificato, anulare e nebbia in tubi e serbatoi.
  • Interazione struttura fluida: deformazione del serbatoio sotto stress termici criogenici e accelerazioni di lancio.

L'approccio multifisico di COMSOL[] permette agli ingegneri di combinare questi fenomeni in un unico modello piuttosto che iterating tra strumenti separati.

Capacità di base di COMSOL CFD per la modellazione criogenica

Modelli di Equazioni e Turbulence

Il modulo CFD risolve le equazioni di Navier‐Stokes per massa, slancio ed energia, accoppiato con un modello di turbolenza appropriato. Per i serbatoi criogenici, il modello k‐ω SST]] si esibisce bene perché risolve flussi di prossimità, critici per prevedere il trasferimento di calore attraverso la parete del serbatoio.

Cambiamento di fase e flusso multifase

COMSOL offre due principali approcci per il cambiamento di fase:
]]] e ] I metodi di campionamento tracciano esplicitamente l'interfaccia di vapore liquido. Questi sono ideali per problemi di superficie libera come lo sloshing del serbatoio o l'impingement del getto.

Termodinamica non-Ideale

La funzione Real Gas[[[]] consente di importare equazioni fondamentali di stato come Peng‐Robinson o Span‐Wagner per previsioni di proprietà altamente accurate.Per idrogeno e ossigeno, queste equazioni sono essenziali vicino al punto critico, dove i presupposti ideali-gas danno errori superiori al 30 %.

Bilancia dei materiali personalizzati

È possibile definire tabelle di proprietà da dati NIST REFPROP[ direttamente in COMSOL, assicurando densità, calore specifico, viscosità e conducibilità termica sono funzioni di temperatura e pressione.

Processo di modellazione passo-passo in COMSOL

1. Definizione della geometria

Per un tipico serbatoio di propellante sferico, includere lo spessore della parete, porte di ingresso/uscita, baffle anti-vortex e strati di isolamento. Semplificare i filetti non critici e i fori di bullone per mantenere la mesh gestibile.

2. Assegnazione di proprietà dei materiali

Creare un materiale personalizzato per il criogenico utilizzando tavoli NIST. Ad esempio, per LH2 a 20 K a 30 K, la densità di importazione in funzione della temperatura e della pressione, e impostare la viscosità dinamica utilizzando il modello Lennard‐Jones. Per la parete del serbatoio (ad esempio, alluminio 2219 o acciaio inossidabile 316L), definire conducibilità termica, densità e coefficiente di espansione termica.

3. Interfacce di Fisica e Coupling

Aggiungere le seguenti interfacce dalla procedura guidata del modello:

  • Flusso laminare[] o Flusso turbolento (k‐ω SST)[ per il movimento fluido.
  • Trasferimento di calore in fluidi[[] e ] Trasferimento di calore in solidi[] per il trasferimento di calore coniugato.
  • Phase Field[[ (o Level Set) per il monitoraggio dell'interfaccia, se si indagano le dinamiche bifase.
  • Non-Isothermal Flow[[[]] multifisica accoppiamento per tenere conto della densità e della galleggiabilità dipendente dalla temperatura.

Per i modelli di ebollizione, aggiungere un [] Fonte di calore[] all'interfaccia pari al calore latente di vaporizzazione tempi la velocità di evaporazione. Questo può essere implementato tramite un'espressione definita dall'utente che si attiva quando la temperatura liquida raggiunge il punto di saturazione alla pressione locale.

4. Condizioni di vita

Condizioni tipiche per un serbatoio di propellante durante una fase di costa:

  • Superficie esterna della parete del serbatoio: flusso di calore da radiazione solare (ad esempio, 1,3 kW/m2 a 1 AU più infrarossi della Terra).
  • Superficie interna della parete del serbatoio: no-slip, flusso di calore coniugato.
  • Interfaccia a vapore liquido: continuità della temperatura, equilibrio di massa con velocità di cambio di fase.
  • Uscita (linea di alimentazione): limite di pressione pari alla pressione del serbatoio meno testa.
  • Ingresso (per riempimento): condizione limite di portata di massa.

Per simulazioni di microgravità, utilizzare una Forza di volume per la tensione superficiale (basata sul modello di forza superficiale continuum) e impostare l'accelerazione a zero o un piccolo valore micro-g (10−4 g).

5. Strategia di fusione

Utilizzare una maglia a strati (almeno 5-10 strati di prisma) alla parete del serbatoio per catturare gradienti di temperatura ripidi. Nella regione dell'interfaccia a due fasi, affinare gli elementi ad una dimensione di circa 1 mm (per un serbatoio di 1 m di diametro) per risolvere la curvatura dell'interfaccia.

Eseguire uno studio di convergenza delle maglie: eseguire il modello con due volte raffinatezza e confrontare i profili di temperatura e le posizioni di interfaccia.

6. Impostazioni di solvente e Passaggio del tempo

I problemi criogenici sono spesso transienti. Usare un ]Backward Differenziazione Formula (BDF)]] risolutore di ordine 1 o 2 con un passo iniziale di 0,01 s. Per la fisica accoppiata, abilitare il Completamente accoppiato] risolutore con un fattore di tempo di campio di campionamento di 0, per migliorare la convergenza.

Monitorare la convergenza controllando i residui (sotto i 10−4) e le quantità conservate (equilibrio di massa entro lo 0,1 %).

7. Procedimento e analisi post-processura

Estratti chiave per estrarre:

  • I contorni di temperatura[] su pareti di serbatoio e dominio fluido.
  • Frequenza di scarico[[ (kg/s) integrata sull'interfaccia.
  • Pressure evolution[[]] nel volume della ninnage (importante per la progettazione strutturale del serbatoio).
  • Streamlines or pathlines[[]] che mostra zone di ricircolo che possono causare stratificazione termica.
  • Posizione interfaccia[[]] in funzione del tempo, critica per la valutazione dei dispositivi di acquisizione liquida.

Gli strumenti COMSOL Valori derivati[] possono calcolare automaticamente la perdita totale di calore, la massa di vapore generata e la frazione di vuoto.

Case study: Simulare Boil‐Off in un serbatoio criogenico

Considerate un serbatoio sferico LH2 (radio 1,5 m) con pareti in alluminio spesse 1 mm e 5 cm di MLI. Il serbatoio è pieno del 50 %, a una pressione di 3 bar, con il fluido a temperatura di saturazione (circa 23 K). Il flusso termico esterno è di 150 W/m2 rispetto all’ambiente termico della sonda. L’obiettivo è quello di prevedere la velocità di ebollizione su un periodo di 24 ore.

Dopo aver impostato il modello come descritto sopra, la simulazione mostra che il MLI riduce la perdita di calore a circa 8 W/m2. La velocità di ebollizione-off si stabilizza a 0.25 g/s dopo un transito iniziale di 200 s. Le spazza parametriche rivelano che raddoppiare lo spessore MLI riduce il grado di ebollizione del 45 %, ma la penalità di peso potrebbe non essere accettabile.

Tali simulazioni informano direttamente il progetto della missione, ad esempio, possono aiutare a dimensionare un sistema di sfiato termodinamico (TVS) che raffredda attivamente il serbatoio per evitare la sovrapressione.

Applicazioni nel ciclo di vita della Spacecraft

Progettazione del sistema di alimentazione propellente

Le linee di alimentazione criogeniche devono trasferire liquido senza cavitazione o flusso di due fasi che potrebbero affamare il motore. I modelli COMSOL della linea di alimentazione, comprese le curve, le giunzioni e le valvole, possono prevedere le gocce di pressione e rilevare le regioni in cui si può verificare l'ebollizione a causa di punti caldi locali.

Convalida del sistema di protezione termica

I materiali di isolamento, inclusi gli aerogel e le coperte MLI, possono essere modellati come supporti porosi o solidi strati con efficaci conducibilità termica.

Avviare le operazioni di riempimento e di Chill‐Down

Il riempimento di un serbatoio con criogenico dalle condizioni ambientali produce violenti scosse e bollenti termiche. Le simulazioni transitorie in COMSOL possono valutare il tempo necessario per raffreddare le pareti del serbatoio a temperature sicure (spesso inferiori a 100 K) prima del riempimento principale. Questi modelli aiutano anche a ottimizzare il tasso di riempimento per evitare una formazione eccessiva del vapore che può danneggiare le valvole.

Gestione del Propellente In-Flight

La fusione in microgravità provoca il movimento liquido che può destabilizzare una sonda spaziale. Le simulazioni di struttura a fluido accoppiato in COMSOL consentono agli ingegneri di progettare baffle, diaframmi o serbatoi di tensione superficiale che mantengono il liquido insediato. Il metodo Phase Field[[]]] cattura l'interfaccia mobile durante le bruciature del motore o i propulstori di atteggiamento, fornendo forze sulle pareti del serbatoio per l'analisi strutturale.

Correlazioni di test di terra

Prima di potersi fidare di una simulazione per il volo, deve essere convalidato contro i dati di test di terra. Il modulo di stima dei parametri COMSOL può regolare i parametri incerti (ad esempio, resistenza ai contatti, fattori di visualizzazione delle radiazioni) per abbinare le letture di termocoppie da un test di sottoscala.

Risorse esterne per argomenti avanzati

Per gli ingegneri che cercano di approfondire la loro conoscenza criogenica CFD, si raccomandano le seguenti risorse:

Istruzioni future: Gemelli digitali e apprendimento automatico

La prossima frontiera nella modellazione dei fluidi criogenici è la creazione di gemelli digitali in tempo reale. I modelli di ordine ridotto (ROM) possono essere costruiti da spazza parametrici COMSOL, poi incorporati nel software di volo del veicolo per la stima dello stato. Ad esempio, una ROM che prevede la velocità di ebollizione come funzione di perdita di calore, livello di riempimento e accelerazione può funzionare in millisecondi, consentendo l'andatura propellent a bordo.

Le reti neurali addestrate sui database di simulazione COMSOL possono prevedere il collasso della pressione dopo le manovre di riorientamento più velocemente di un modello multifisico completo. Questi approcci ibridi mantengono l'accuratezza mentre il taglio dei costi computazionali da ordini di grandezza.

LiveLink for MATLAB[]] rende semplice interfacciare i dati di simulazione con tali strumenti AI, permettendo agli ingegneri di migliorare in modo iterativo le loro ROM con nuovi dati ad alta fedeltà come la missione progredisce.

Conclusioni

La modellazione della fluidodinamica criogenica in COMSOL CFD fornisce agli ingegneri spaziali un ambiente potente e unificato per affrontare le sfide multifisiche dei flussi a bassa temperatura.Dall'iniziale progettazione del concetto attraverso le operazioni di volo, queste simulazioni consentono scelte di isolamento migliori, sistemi di alimentazione più efficienti e una gestione più sicura del propellente.

Gli esempi e i flussi di lavoro presentati qui offrono un punto di partenza per la costruzione di modelli CFD criogenici affidabili. L’esplorazione spaziale spinge verso uno spazio più profondo e missioni più lunghe, la gestione accurata dei fluidi criogenici diventa ancora più critica.