Modello Lattice Boltzmann per la dinamica dei fluidi comprimibili non ideali

Informazioni Fondamentali

• ID Articolo: 2510.14712
• Titolo: Lattice Boltzmann model for non-ideal compressible fluid dynamics
• Autori: S. A. Hosseini, M. Feinberg, I. V. Karlin (ETH Zurich)
• Classificazione: physics.flu-dyn
• Rivista di Pubblicazione: Journal of Fluid Mechanics
• Link Articolo: https://arxiv.org/abs/2510.14712

Riassunto

Il presente articolo propone un metodo Lattice Boltzmann per simulare il flusso di fluidi comprimibili non ideali. Il metodo impiega un reticolo a primi vicini, introducendo una serie coerente di termini correttivi attraverso un attrattore quasi-equilibrio, garantendo un tasso di dissipazione Navier-Stokes positivo definito e invariante galileiano. Questa costruzione evita la necessità di modelli estesi o regolarizzazioni temporanee, mantenendo al contempo la stabilità numerica e la coerenza termodinamica in un ampio intervallo di stati di flusso. Il modello risultante riproduce accuratamente la fluidodinamica ai livelli euleriano e Navier-Stokes. Come verifica rigorosa, gli autori dimostrano per la prima volta nel framework Lattice Boltzmann la simulazione quantitativamente accurata dell'interazione goccia-onda d'urto con numero di Mach fino a 1,47.

Contesto e Motivazione della Ricerca

Definizione del Problema

La dinamica dei fluidi comprimibili non ideali rappresenta un ramo emergente e in rapido sviluppo della meccanica dei fluidi, principalmente dovuto all'emergere di metodologie e tecnologie che operano in condizioni prossime al punto critico, transcritiche e supercritiche. Questa classe di fluidi include:

1. Fluidi supercritici: fluidi a pressioni e temperature superiori al punto critico
2. Vapori densi: gas il cui fattore di comprimibilità si discosta dal comportamento del gas ideale
3. Fluidi BZT: fluidi a complessità molecolare elevata con derivata negativa della dinamica dei gas fondamentale

Importanza

1. Esigenze di applicazioni industriali: tecnologie energetiche come i cicli Rankine organici e le turbine CO₂ supercritiche
2. Scarsità di dati sperimentali: configurazioni sperimentali complesse e difficoltà nell'acquisizione dei dati
3. Necessità di strumenti numerici: esigenza di strumenti numerici coerenti ed efficienti per comprendere fenomeni fisici complessi

Limitazioni dei Metodi Esistenti

1. Metodi CFD tradizionali: nella regione spinodale termodinamicamente instabile, il quadrato della velocità del suono adiabatica diventa negativo, causando l'evoluzione dell'equazione da iperbolica a parabolica
2. Modelli LBM esistenti: principalmente orientati verso applicazioni di ebollizione, mancano di modelli sistematici per flussi comprimibili non ideali
3. Problemi di stabilità numerica: sfide di stabilità a numeri di Mach elevati e con forti effetti non ideali

Contributi Principali

1. Propone un nuovo modello LBM a doppia funzione di distribuzione: implementa la simulazione di flussi non ideali comprimibili utilizzando il classico reticolo a primi vicini (D3Q27)
2. Introduce il meccanismo dell'attrattore quasi-equilibrio: garantisce attraverso i termini correttivi fᵢ e gᵢ un tasso di dissipazione positivo definito e invariante galileiano
3. Realizza il controllo indipendente della viscosità di volume: risolve il problema della possibile viscosità di volume negativa nella struttura BGK
4. Realizza per la prima volta la simulazione dell'interazione goccia-onda d'urto ad alto numero di Mach: raggiunge una simulazione quantitativamente accurata con Ma=1,47 nel framework LBM
5. Mantiene la coerenza termodinamica: preserva la stabilità numerica in un ampio intervallo di stati di flusso

Dettagli Metodologici

Definizione del Compito

Simulare fluidi comprimibili non ideali che soddisfano il sistema di equazioni Navier-Stokes-Korteweg (NSK):

• Conservazione della massa: ∂ₜρ + ∇·(ρu) = 0
• Conservazione della quantità di moto: ∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
• Conservazione dell'energia: ∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

dove T^K è il tensore di tensione superficiale di Korteweg e T^NS è il tensore degli sforzi viscosi.

Architettura del Modello

1. Progettazione a Doppia Funzione di Distribuzione

Utilizza due insiemi di funzioni di distribuzione {fᵢ, gᵢ}, i=1,...,27:

• fᵢ: gestisce la conservazione della massa e della quantità di moto
• gᵢ: gestisce la conservazione dell'energia

2. Equazione Lattice Boltzmann

fᵢ(x+cᵢδt, t+δt) = fᵢ + 2β(fᵢᵉᵠ - fᵢ) + (1-β)(fᵢ* - fᵢᵉᵠ)
gᵢ(x+cᵢδt, t+δt) = gᵢ + 2β(gᵢᵉᵠ - gᵢ) + (1-β)(gᵢ* - gᵢᵉᵠ)

dove β = Pδt/(2μ + Pδt) è il parametro di rilassamento.

3. Funzioni di Distribuzione di Equilibrio

Adotta una forma fattorizzata:

fᵢᵉᵠ = ρ ∏_α Ψᵢα(uα, P/ρ + uα²)

dove Ψᵢα è la funzione di base:

Ψᵢα(ξα, ζαα) = 1 - cᵢα² + ½[(3cᵢα² - 2)ζαα + cᵢαξα]

4. Termini Correttivi Quasi-Equilibrio

Correzione della funzione di distribuzione f:

• Termine di forza: Fα = -κρ∂α∇²ρ (forza di Korteweg)
• Correzione dello sforzo: Φαα contiene la correzione del tensore degli sforzi diagonali e il controllo indipendente della viscosità di volume

Correzione della funzione di distribuzione g:

• Velocità traslata: u* = u + δt F/ρ
• Temperatura traslata: T* = T - δt²F·F/(2ρ²cᵥ)
• Correzione del flusso di calore: considera gli effetti del gradiente di entalpia per gas non ideali

Punti di Innovazione Tecnica

1. Meccanismo dell'Attrattore Quasi-Equilibrio

Introducendo fᵢ e gᵢ, il sistema non si rilassa semplicemente verso l'equilibrio locale, ma verso uno stato quasi-equilibrio modificato, il che garantisce:

• Dinamica macroscopica corretta
• Stabilità numerica
• Invarianza galileiana

2. Controllo Indipendente della Viscosità di Volume

Nel modello BGK tradizionale la viscosità di volume η' = ((D+2)/D - ρcs²/P)μ può essere negativa. Attraverso il termine correttivo Φ' si realizza il controllo indipendente:

Φ' = P((D+2)/D - ρcs²/P - η/μ)(∇·u)

3. Coerenza Termodinamica

Attraverso la progettazione accurata della temperatura traslata T* e della correzione energetica, si garantisce la coerenza termodinamica sotto equazioni di stato non lineari.

Configurazione Sperimentale

Casi di Verifica

1. Verifica Fondamentale:
   • Misurazione della velocità del suono
   • Verifica della viscosità di taglio
   • Flusso termico di Couette
   • Dissipazione del modo normale
2. Verifica del Flusso Multifase:
   • Densità di coesistenza liquido-vapore
   • Coerenza dell'interfaccia e convergenza
3. Verifica del Flusso Comprimibile:
   • Tubo d'urto unidimensionale non ideale
   • Interazione onda d'urto-colonna di liquido

Parametri Materiali

Principalmente azoto (N₂) utilizzato per la verifica:

• Pressione critica: Pc = 3,4×10⁶ Pa
• Densità critica: ρc = 241,96 kg/m³
• Temperatura critica: Tc = 126,2 K
• R/cv = 0,4

Equazione di Stato

Adotta l'equazione di stato di van der Waals:

P(ρ,T) = ρRT/(1-bρ) - aρ²

Risultati Sperimentali

Risultati Principali

1. Verifica delle Grandezze Fisiche Fondamentali

• Velocità del suono: sui rami di liquido saturo e vapore saturo, perfetto accordo con la soluzione analitica nell'intervallo Tr∈0,7,1
• Viscosità di taglio: mantiene l'invarianza galileiana a diversi numeri di Mach, errore <1%
• Flusso termico di Couette: coerenza con la soluzione analitica per Pr∈{0,6,1,2,4,9} e Ma∈{0,8,1,2,1,6}

2. Proprietà del Flusso Multifase

• Densità di coesistenza: perfetta coerenza con la costruzione di area uguale di Maxwell nell'intervallo Tr∈0,3,0,99
• Convergenza dell'interfaccia: mostra buone caratteristiche di convergenza a diverse risoluzioni di griglia

3. Tubo d'Urto Non Ideale

Simulazione riuscita di tre diverse configurazioni di tubo d'urto non ideale:

• Caso I: onda di rarefazione nella regione Γ<0
• Caso II: comportamento classico dell'onda di compressione
• Caso III: onda di rarefazione mista che attraversa dinamicamente la linea Γ=0

Tutti i casi mostrano eccellente accordo con i dati di riferimento.

4. Interazione Onda d'Urto-Goccia

• Prima realizzazione nel framework LBM della simulazione ad alto numero di Mach Ma=1,47
• Cattura accurata di complesse strutture d'onda come onda incidente, onda trasmessa, onda riflessa, barra di Mach
• Evoluzione della deformazione della goccia altamente coerente con i dati sperimentali

Esperimenti di Ablazione

Verifica della convergenza del modello attraverso test di interfaccia a diverse risoluzioni di griglia:

• Quattro livelli di risoluzione con δx da 5μm a 0,1μm
• Mostra buone caratteristiche di convergenza della griglia

Risultati Chiave

1. Stabilità numerica: mantiene la stabilità a numeri di Mach elevati e con forti effetti non ideali
2. Efficienza computazionale: utilizza il reticolo standard a primi vicini, costi computazionali controllabili
3. Fedeltà fisica: riproduce accuratamente i fenomeni fisici complessi dei fluidi non ideali

Lavori Correlati

Modelli LBM Tradizionali Multifase

1. Modello a gradiente di colore (Gunstensen et al. 1991)
2. Modello pseudo-potenziale (Shan & Chen 1993)
3. Modello ad energia libera (Swift et al. 1996)

Sviluppo LBM Comprimibile

• He & Doolen (2002): primo tentativo di flusso non ideale comprimibile
• I modelli esistenti sono principalmente orientati verso applicazioni di ebollizione
• Mancanza di modelli sistematici per flussi non ideali comprimibili

Vantaggi Relativi di Questo Articolo

1. Approccio sistematico: framework completo a doppia funzione di distribuzione
2. Fondamenti teorici: analisi rigorosa multi-scala
3. Completezza della verifica: verifica completa dalle grandezze fisiche fondamentali ai fenomeni di flusso complessi

Conclusioni e Discussione

Conclusioni Principali

1. Sviluppo riuscito di un modello LBM applicabile a fluidi comprimibili non ideali
2. Realizzazione della stabilità numerica e della coerenza fisica attraverso il meccanismo dell'attrattore quasi-equilibrio
3. Prima realizzazione nel framework LBM della simulazione dell'interazione goccia-onda d'urto ad alto numero di Mach
4. Il modello mantiene la coerenza termodinamica in un ampio intervallo di stati di flusso

Limitazioni

1. Modello a tempo di rilassamento singolo: attualmente utilizza il semplice modello di collisione BGK, potrebbe non essere sufficientemente robusto per la simulazione della turbolenza
2. Limitazioni dell'equazione di stato: la verifica principale riguarda l'equazione di van der Waals, altre equazioni di stato complesse richiedono ulteriore verifica
3. Verifica tridimensionale limitata: la maggior parte della verifica è concentrata su configurazioni unidimensionali e bidimensionali

Direzioni Future

1. Estensione a modelli di collisione più robusti: per la simulazione di configurazioni turbolente
2. Equazioni di stato più complesse: come Peng-Robinson e altre
3. Applicazioni tridimensionali su larga scala: applicazione a problemi ingegneristici reali

Valutazione Approfondita

Punti di Forza

1. Rigore teorico: analisi multi-scala completa e derivazioni matematiche rigorose
2. Innovazione metodologica: il meccanismo dell'attrattore quasi-equilibrio rappresenta un'innovazione importante
3. Completezza della verifica: verifica sistematica dalla fisica fondamentale alle applicazioni complesse
4. Valore pratico: risolve importanti problemi di applicazione ingegneristica
5. Efficienza numerica: l'utilizzo di reticoli standard evita la complessità computazionale

Carenze

1. Semplicità del modello di collisione: il modello BGK ha applicabilità limitata a flussi complessi
2. Dipendenza dai parametri: la scelta di alcuni parametri correttivi manca di guida sistematica
3. Insufficienza della verifica tridimensionale: necessità di ulteriore verifica con geometrie complesse tridimensionali

Impatto

1. Contributo accademico: apre nuove direzioni per l'applicazione di LBM nel campo dei flussi comprimibili non ideali
2. Valore ingegneristico: fornisce importanti strumenti numerici per le tecnologie a fluidi supercritici
3. Significato metodologico: il concetto di attrattore quasi-equilibrio potrebbe influenzare altri metodi numerici

Scenari Applicabili

1. Tecnologie a fluidi supercritici: turbine CO₂, cicli Rankine organici
2. Flussi multifase: rottura di gocce, processi di cambio di fase
3. Dinamica delle onde d'urto: interazioni complesse di onde in flussi ad alta velocità
4. Simulazione numerica diretta: ricerca fisica fondamentale che richiede alta precisione

Bibliografia

1. Guardone, A. et al. (2024). Nonideal compressible fluid dynamics of dense vapors and supercritical fluids. Annual Review of Fluid Mechanics, 56, 241-269.
2. Hosseini, S.A. & Karlin, I.V. (2023). Lattice Boltzmann for non-ideal fluids: Fundamentals and practice. Physics Reports, 1030, 1-137.
3. Shan, X. & Chen, H. (1993). Lattice Boltzmann model for simulating flows with multiple phases and components. Physical Review E, 47(3), 1815.