Real-time identification of parametric sloshing-induced heat and mass transfer in a horizontally oriented cylindrical tank

Questo studio presenta la prima caratterizzazione sperimentale del trasferimento di calore e massa indotto dal moto di *sloshing* parametrico in un serbatoio cilindrico orizzontale, utilizzando un filtro di Kalman esteso (AEKF) per l'identificazione in tempo reale dei coefficienti di scambio termico durante la destratificazione del fluido.

L'essenziale

Il "Ballerino" nei Serbatoi: Perché il carburante spaziale non ama le vibrazioni

Immaginate di avere un grande contenitore cilindrico, come quello che potrebbe trovarsi sotto l'ala di un futuro aereo a idrogeno o in un razzo. Dentro, c'è un liquido molto freddo (criogenico) e, sopra, un po' di vapore.

In condizioni normali, questo sistema è molto ordinato: il liquido sta tranquillo sul fondo, il vapore sta sopra, e c'è una sorta di "stratificazione termica". È come una lasagna: strati ben distinti che non si mescolano.

Il problema: Il "Ballerino" (Lo Sloshing)
Ora, immaginate che il veicolo inizi a vibrare verticalmente, come se stesse correndo su una strada molto sconnessa o subisse delle scosse. Questa vibrazione non sposta il liquido lateralmente, ma lo fa "rimbalzare" su e giù. In fisica, questo fenomeno si chiama sloshing parametrico.

È come se cercaste di far ballare una lasagna: se la scuotete nel modo giusto, gli strati smettono di essere separati e tutto inizia a mescolarsi violentemente.

Cosa succede quando la "lasagna" si rompe?
Quando la vibrazione supera una certa soglia (una sorta di "punto di non ritorno"), il liquido non fa più solo piccole onde, ma crea dei veri e propri getti che colpiscono il soffitto del serbatoio. In quel momento succede il caos termodinamico:

1. Il Grande Mix: Il liquido freddo e il vapore caldo si scontrano e si mescolano.
2. Il Crollo della Pressione: Questo scontro improvviso fa sì che il vapore si trasformi rapidamente in liquido (condensazione). Immaginate di aspirare improvvisamente l'aria da una stanza: la pressione crolla. Questo è pericoloso perché i serbatoi devono mantenere una pressione stabile per essere sicuri.

Cosa hanno fatto gli scienziati? (L'investigatore invisibile)
Gli autori di questo studio hanno costruito un serbatoio speciale e lo hanno fatto "ballare" in laboratorio usando un fluido che si comporta in modo simile all'idrogeno.

Ma c'è una sfida: come facciamo a sapere esattamente quanto è forte lo scambio di calore mentre tutto sta scuotendo e ribollendo? Non possiamo mettere mille termometri dentro un serbatoio in movimento senza rovinare l'esperimento.

Qui entra in gioco il loro "superpotere": un algoritmo chiamato Filtro di Kalman (AEKF).
Pensate al Filtro di Kalman come a un investigatore privato molto intelligente. Lui non vede tutto direttamente, ma osserva solo pochi indizi (pressione e temperatura in pochi punti) e, usando la logica e la matematica, riesce a "indovinare" con una precisione incredibile cosa sta succedendo in ogni singolo centimetro del serbatoio, calcolando in tempo reale quanto calore si sta spostando.

Le scoperte principali (In parole povere):

• La soglia del caos: Esiste un limite preciso di vibrazione. Sotto quel limite, la "lasagna" resta ferma. Sopra quel limite, scoppia il caos.
• Il colpevole della pressione: Hanno scoperto che il calo di pressione non è causato dal calore che entra dalle pareti, ma quasi interamente dal fatto che il vapore si trasforma in liquido a causa del mix causato dalle onde.
• Il muro come "cuscinetto": Le pareti del serbatoio agiscono come una spugna termica, cercando di mitigare lo shock, ma il movimento del liquido è molto più veloce e violento.

Perché è importante?
Se vogliamo volare con aerei che usano idrogeno liquido o mandare missioni spaziali più sicure, dobbiamo sapere esattamente cosa succede quando il serbatoio trema. Questo studio fornisce la "mappa del tesoro" (o meglio, della tempesta) per progettare serbatoi che non esplodano o non perdano pressione quando iniziano a ballare.

Sommario tecnico

Riassunto Tecnico: Identificazione in tempo reale del trasferimento di calore e massa indotto da sloshing parametrico in un serbatoio cilindrico orizzontale

1. Definizione del Problema

Il lavoro affronta le sfide termodinamiche legate allo sloshing parametrico (oscillazioni del liquido indotte da accelerazioni verticali) in serbatoi cilindrici orizzontali parzialmente riempiti. Questa configurazione è critica per lo stoccaggio di propellenti criogenici (come idrogeno liquido o metano) in applicazioni aeronautiche e di trasporto terrestre.

A differenza dello sloshing laterale, l'eccitazione verticale altera periodicamente la forza di gravità efficace, innescando instabilità di tipo Faraday. In condizioni non isoterme, queste oscillazioni rompono la stratificazione termica naturale tra liquido e vapore, intensificando lo scambio di calore e massa. Ciò provoca fluttuazioni di pressione rapide e potenzialmente pericolose, ma la dinamica termodinamica in geometrie cilindriche orizzontali è finora rimasta scarsamente caratterizzata nella letteratura scientifica.

2. Metodologia Sperimentale e Modellistica

Gli autori hanno adottato un approccio metodologico duale per separare la cinematica della superficie libera dalla risposta termodinamica:

• Campagna Isoterma: Utilizzando un modulo in PMMA trasparente e tecniche di video-analisi ad alta velocità, è stata studiata la dinamica della superficie libera. È stata applicata la Proper Orthogonal Decomposition (POD) combinata con la Continuous Wavelet Transform (CWT) per identificare le frequenze naturali, i modi di oscillazione e i fenomeni di competizione tra modi.
• Campagna Non Isoterma: Utilizzando un serbatoio in acciaio AISI 316 strumentato, è stata studiata la risposta termodinamica usando il fluido surrogato 3M™ Novec™ HFE-7000. La strumentazione includeva termocoppie distribuite verticalmente (sia nel liquido che nel vapore) e trasduttori di pressione.
• Modellazione Termodinamica e Inferenza (AEKF): È stato sviluppato un modello termodinamico a parametri concentrati (lumped-parameter) che include i bilanci energetici di vapore, liquido e parete. Per stimare i coefficienti di scambio termico (Nusselt numbers) in tempo reale da dati rumorosi, è stato impiegato un Augmented Extended Kalman Filter (AEKF). Questo metodo tratta i coefficienti di scambio come stati nascosti che evolvono dinamicamente, permettendo un'inferenza adattiva durante i regimi transitori.

3. Contributi Chiave

1. Prima caratterizzazione sperimentale completa: Il lavoro fornisce la prima analisi del calo di pressione e del trasferimento di calore associato allo sloshing verticale in un cilindro orizzontale 3D.
2. Nuovo metodo di inferenza: L'uso dell'AEKF per l'analisi del trasferimento di calore e massa in sistemi non isotermi rappresenta un avanzamento rispetto ai metodi di ottimizzazione a orizzonte completo, permettendo stime online con quantificazione dell'incertezza.
3. Mappatura della stabilità parametrica: Identificazione delle regioni di instabilità (Mathieu tongues) e della sensibilità della frequenza naturale al livello di riempimento (fill level).

4. Risultati Principali

• Dinamica del Fluido: È stata identificata una soglia critica di accelerazione. Sotto tale soglia, il fluido rimane stratificato; sopra di essa, si innesca la risonanza parametrica. A un livello di riempimento del 50%, è stata osservata una competizione tra modi: l'oscillazione del modo simmetrico (2,0) (jet centrale) interagisce con il modo antisimmetrico (1,0) (onda planare), portando a fenomeni di risonanza subarmonica (period-doubling).
• Risposta Termodinamica: Lo sloshing induce una "destratificazione termica" completa. Il processo segue una sequenza a sei fasi: (1) crescita del jet, (2) burst di pressione seguito da un rapido calo dovuto alla condensazione, (3) risposta conduttiva della parete, (4) scambio calore parete-liquido, (5) miscelazione completa dei fluidi e (6) recupero della pressione.
• Meccanismi di Pressione: L'analisi della derivata della pressione ha dimostrato che il calo di pressione è guidato prevalentemente dal cambiamento di fase (condensazione), che contribuisce per circa l'80% alla variazione totale, mentre la variazione della temperatura del vapore ha un ruolo secondario e smorzante.
• Coefficienti di Nusselt: L'AEKF ha rivelato che, al momento della destratificazione, i numeri di Nusselt aumentano di due-quattro ordini di grandezza, segnando il passaggio dalla convezione naturale alla convezione forzata intensa.

5. Significato e Applicazioni

I risultati forniscono un benchmark fondamentale per la progettazione di sistemi di stoccaggio criogenico. La comprensione di come l'accelerazione verticale influenzi la stabilità termica e la pressione interna è essenziale per garantire la sicurezza e l'efficienza dei serbatoi di propellente in aerei a idrogeno, veicoli spaziali e mezzi di trasporto pesanti. La metodologia AEKF proposta offre uno strumento robusto per il monitoraggio della salute termica dei serbatoi in tempo reale.