Lattice Boltzmann model for non-ideal compressible fluid dynamics

Questo articolo introduce un nuovo modello di Boltzmann su reticolo che utilizza reticoli di primi vicini e termini di correzione quasi-equilibrio per simulare flussi di fluidi comprimibili e non ideali con coerenza termodinamica e stabilità numerica, validando con successo la sua accuratezza attraverso simulazioni quantitative delle interazioni tra onde d'urto e gocce a numeri di Mach fino a 1,47.

L'essenziale

Immagina di voler simulare il comportamento di un fluido su un computer. Per lungo tempo, i computer sono stati eccellenti nel simulare fluidi "ideali", come l'acqua che scorre dolcemente in un fiume o l'aria che si muove lentamente attorno a un'ala. Questi fluidi seguono regole semplici e prevedibili.

Ma cosa succede quando il fluido è sottoposto a pressioni e temperature estreme, comportandosi come un gas denso che è quasi un liquido, o come un liquido che è quasi un gas? Questo è il mondo dei fluidi comprimibili non ideali. Pensalo come un fluido "stressato" che si comporta in modo strano, rifiutandosi di seguire le semplici regole del mondo ideale. Questo accade in tecnologie avanzate come le turbine a CO2 supercritica e i cicli energetici organici.

Il problema è che gli strumenti informatici esistenti faticano con questi fluidi stressati. O si bloccano, o forniscono risposte errate, o richiedono una potenza di calcolo così massiccia da diventare inutili.

Questo articolo introduce un modo nuovo e più intelligente per simulare questi fluidi complicati utilizzando un metodo chiamato Metodo Lattice Boltzmann (LBM). Ecco come funziona il nuovo approccio degli autori, spiegato attraverso semplici analogie:

1. Il sistema "a due binari"

La maggior parte dei vecchi metodi di simulazione cerca di tracciare tutto (massa, velocità, energia) con un unico, complicato insieme di regole. Gli autori hanno capito che questo era come cercare di guidare un'auto mentre si contemporaneamente si fanno saltare una dozzina di palline: diventa disordinato e instabile.

Invece, hanno costruito un sistema a due binari:

• Binario A (La folla): Un insieme di regole traccia la densità e la velocità del fluido (quanti particelle ci sono e dove stanno andando).
• Binario B (L'energia): Un secondo insieme di regole, separato, traccia l'energia totale.
  Separando questi elementi, il computer non si confonde. È come avere un controllore del traffico dedicato per le auto e uno separato per il carburante, assicurandosi che nessun sistema faccia crashare l'altro.

2. L'attrattore "quasi-equilibrio"

In fisica, i fluidi vogliono naturalmente stabilizzarsi in uno stato calmo chiamato "equilibrio". Tuttavia, in queste condizioni estreme, il fluido viene costantemente spinto e tirato, quindi non riesce mai a stabilizzarsi completamente.

Gli autori hanno inventato un trucco intelligente chiamato "attrattore quasi-equilibrio".

• L'analogia: Immagina un cane che insegue una palla. La palla rappresenta lo "stato di calma perfetta". Il cane rappresenta il fluido. In una situazione normale, il cane corre dritto verso la palla.
• Il problema: In questo fluido estremo, la palla continua a spostarsi o a cambiare forma. Se il cane inseguisse la palla alla cieca, potrebbe cadere da una scogliera (la simulazione diventerebbe instabile).
• La soluzione: Gli autori hanno fornito al cane un "GPS" che prevede dove la palla sarà un millesimo di secondo dopo, basandosi su come il vento (pressione) e il terreno (densità) stanno cambiando. Questo bersaglio "spostato" permette al cane di correre senza cadere dalla scogliera. Questo garantisce che la simulazione rimanga stabile anche quando il fluido si muove molto velocemente o cambia densità rapidamente.

3. Correggere il calore "spurio"

Quando i fluidi si muovono velocemente, generano calore. Nei modelli informatici standard, il calore a volte fluisce nella direzione sbagliata o crea un falso calore "fantasma" che non esiste nella realtà.

• L'analogia: È come un termostato che pensa che la stanza sia gelida perché sta misurando la corrente d'aria da una finestra, non la temperatura effettiva della stanza.
• La correzione: Gli autori hanno aggiunto un specifico "termine di correzione" alle loro equazioni. Questo agisce come un filtro che rimuove le correnti d'aria false, assicurando che il calore fluisca esattamente come dictated dalla fisica (Legge di Fourier), anche in queste condizioni estreme e non ideali.

4. I test della "tubo a shock" e della "colonna liquida"

Per dimostrare che il loro nuovo metodo funziona, non hanno fatto solo matematica; hanno eseguito test estremi:

• Il tubo a shock: Hanno simulato un'esplosione improvvisa di pressione (un'onda d'urto) che si muove attraverso un gas denso. Nei gas normali, queste onde si comportano in un certo modo. In questi gas "non ideali", le onde possono fare qualcosa di strano: uno "shock di rarefazione" (un'onda che si espande ma agisce comunque come uno shock netto). Il loro modello ha previsto con successo questo comportamento strano, che i vecchi modelli avevano mancato.
• La colonna liquida: Hanno simulato un'onda d'urto ad alta velocità che colpisce una goccia di liquido. Questo è un test molto difficile perché lo shock rimbalza, si riflette e strappa la goccia in pezzi. Il loro modello ha gestito la collisione perfettamente, corrispondendo agli esperimenti reali in cui la goccia di liquido si appiattisce e si espande esattamente come dovrebbe.

Perché questo è importante

Gli autori affermano che il loro metodo è veloce, stabile e accurato. Utilizza una griglia semplice (come una scacchiera standard) invece di aver bisogno di una griglia super-complessa e allungata. Questo significa che gli scienziati possono ora simulare questi flussi di fluidi non ideali, estremi e ad alta velocità su computer standard con alta precisione.

In sintesi: L'articolo presenta un nuovo "manuale di guida" per le simulazioni informatiche che permette loro di gestire fluidi sotto stress estremo senza bloccarsi. Utilizzando un sistema a due binari e un intelligente "GPS" per l'energia del fluido, possono prevedere con precisione come questi fluidi complessi si comportano in scenari ad alta velocità, aprendo la strada a migliori progettazioni per sistemi energetici avanzati.

Sommario tecnico

Sintesi Tecnica: Modello Lattice Boltzmann per la Dinamica dei Fluidi Compressibili Non Ideali

Enunciato del Problema
La dinamica dei fluidi compressibili non ideali (NICFD) affronta i comportamenti dei fluidi nei regimi vicino-, trans- e super-critici, dove i fluidi si discostano significativamente dalle leggi dei gas ideali. Questi regimi sono critici per le tecnologie energetiche moderne, come i cicli Rankine organici e le turbine a CO₂ supercritica. Tuttavia, i dati sperimentali in questi stati sono scarsi e difficili da acquisire. Sebbene la Simulazione Numerica Diretta (DNS) sia necessaria per comprendere la fisica complessa (ad esempio, strutture d'urto, turbolenza in vapori densi) e generare database ingegneristici, gli strumenti numerici esistenti spesso faticano con la coerenza termodinamica, la stabilità e il recupero accurato dei limiti idrodinamici per equazioni di stato (EOS) non ideali. Nello specifico, i metodi convenzionali di Fluidodinamica Computazionale (CFD) affrontano sfide nella regione spinodale termodinamicamente instabile, dove il quadrato della velocità del suono adiabatica diventa negativo, richiedendo trattamenti speciali come la regola dell'area uguale di Maxwell. Il Metodo Lattice Boltzmann (LBM), pur essendo riuscito per flussi incomprimibili e multifase, ha storicamente trascurato la termodinamica compressibile non ideale nelle sue formulazioni cinetiche, affidandosi spesso a regolarizzazioni ad hoc o stencil estesi che compromettono l'efficienza o l'invarianza galileiana.

Metodologia
Gli autori propongono un nuovo modello cinetico e la sua realizzazione Lattice Boltzmann progettata per recuperare le equazioni di Navier–Stokes–Korteweg (NSK) comprimibili per fluidi non ideali. La metodologia si basa sui seguenti pilastri:

1. Approccio a Doppia Funzione di Distribuzione (DDF): Il modello impiega due funzioni di distribuzione: $f$ per massa e quantità di moto, e $g$ per l'energia di massa ($\rho E$). Questo disaccoppiamento consente l'uso di reticoli classici di primi vicini (specificamente D3Q27) senza richiedere stencil estesi.
2. Attrattori di Quasi-Equilibrio: L'operatore di collisione utilizza una struttura simile a BGK a doppio rilassamento che coinvolge un attrattore di equilibrio locale ($f^{eq}, g^{eq}$) e un attrattore di quasi-equilibrio spostato ($f^\star_\lambda, g^\star_\lambda$).
   • L'equilibrio locale è definito utilizzando un parametro di temperatura di riferimento $\theta = P/\rho$ (lavoro di flusso termodinamico) piuttosto che la temperatura termodinamica $T$, garantendo la corretta larghezza della distribuzione di Maxwell indipendente dall'energia interna.
   • L'equilibrio spostato introduce correzioni per la forza di Korteweg (tensione superficiale), la viscosità di massa e il flusso di calore. Impiega valori spostati per la velocità di flusso ($u^\star$), la temperatura di riferimento ($\theta^\star$) e la temperatura termodinamica ($T^\star$).
3. Coerenza Termodinamica:
   • La temperatura di riferimento $\theta$ è impostata a $P/\rho$, una formulazione che garantisce che il modello rimanga valido per EOS non ideali (ad esempio, van der Waals).
   • Viene introdotto un termine di correzione nel flusso di calore ($q^c_\lambda$) per compensare le componenti spurie non di Fourier che sorgono dai gradienti di densità nei fluidi non ideali, garantendo il recupero della legge di Fourier.
   • Lo spostamento della temperatura di riferimento consente alla viscosità di massa di essere un parametro indipendente, sintonizzabile e definito positivo, aggirando il problema della viscosità di massa negativa che può sorgere nei modelli BGK standard per pressioni non ideali.
4. Discretizzazione: Il modello utilizza uno schema esplicito di accuratezza del secondo ordine derivato dall'integrazione trapezoidale lungo le caratteristiche. Impiega un reticolo standard D3Q27 e gestisce i termini sorgente (forza di Korteweg) tramite l'equilibrio spostato, evitando errori di discretizzazione dei termini sorgente espliciti.

Contributi Chiave

• Nuovo Quadro Cinetico: Il documento introduce un modello cinetico che recupera le equazioni idrodinamiche target (bilanci di massa, quantità di moto ed energia di massa) per fluidi compressibili non ideali utilizzando uno stencil cinetico minimo (primi vicini).
• Dissipazione Definita Positiva: Costruendo il modello tramite attrattori di quasi-equilibrio, gli autori garantiscono che i tassi di dissipazione di Navier–Stokes (viscosità di taglio e di massa) siano definiti positivi e invarianti galileiani, un punto di fallimento comune nei precedenti tentativi LBM per flussi compressibili non ideali.
• Robustezza Termodinamica: Il modello gestisce naturalmente le regioni di instabilità termodinamica (spinodale) dove la velocità del suono diventa immaginaria nel limite non viscoso, poiché l'LBM eredita la propagazione lungo caratteristiche discrete fisse piuttosto che fare affidamento su solver di PDE iperboliche che richiedono trattamenti speciali in queste regioni.
• Validazione di Fenomeni Non Classici: Lo studio fornisce la prima validazione quantitativa dell'LBM per le interazioni urto-goccia a numeri di Mach fino a 1,47 in un contesto di fluido non ideale, dimostrando la capacità del modello di catturare strutture d'onda complesse.

Risultati
Il modello è stato validato attraverso una gerarchia di configurazioni sempre più complesse:

1. Stabilità Lineare e Dispersione: Le simulazioni della velocità del suono, della dissipazione dell'onda di taglio e della conducibilità termica hanno mostrato un eccellente accordo con le soluzioni analitiche per un fluido di van der Waals adattato alle proprietà dell'azoto su una gamma di numeri di Mach e temperature ridotte.
2. Coerenza Multifase: Il modello ha riprodotto accuratamente le densità di coesistenza liquido-vapore e i profili dell'interfaccia, dimostrando una convergenza del secondo ordine con il raffinamento della griglia e l'accordo con la regola dell'area uguale di Maxwell.
3. Tubi a Urto Non Classici: Sono stati simulati tre casi di tubo a urto 1-D per sondare la dinamica delle onde non classica caratterizzata dal derivato fondamentale della dinamica dei gas ($\Gamma$).
   • I casi con $\Gamma < 0$ hanno catturato correttamente la formazione di urti di rarefazione (onde di espansione che si affilano) e ventagli di compressione (onde di compressione che si allargano), fenomeni distinti dalla dinamica dei gas classica.
   • Il modello è transitato con successo tra regimi classici ($\Gamma > 0$) e non classici ($\Gamma < 0$) all'interno di una singola simulazione.
4. Interazione Urto-Colonna Liquida: È stata eseguita una simulazione 2D di un'onda d'urto planare che interagisce con una colonna di liquido saturo vicino al punto critico ($M_a = 1,47$).
   • Il modello ha catturato strutture d'onda complesse, inclusi urti trasmessi, urti riflessi, steli di Mach e punti tripli.
   • Il confronto quantitativo dell'evoluzione della larghezza della colonna liquida rispetto ai dati sperimentali e ai risultati numerici precedenti ha mostrato un buon accordo.
   • L'analisi termodinamica ha rivelato che l'interazione dell'urto spinge il vapore circostante in uno stato supercritico mentre la colonna liquida rimane marginalmente subcritica, catturando accuratamente il forte accoppiamento termodinamico.

Significato e Affermazioni
Gli autori affermano che questo lavoro estende l'applicabilità dei metodi Lattice Boltzmann ai flussi compressibili non ideali ad alta velocità con uno stencil cinetico minimo. Il significato principale risiede nella fornitura di un quadro termodinamicamente coerente e numericamente stabile che non fa affidamento su stencil estesi o regolarizzazioni ad hoc. Il modello è presentato come una teoria cinetica ridotta che mira al limite idrodinamico NSK, capace di gestire l'intera gamma di comportamenti non ideali, inclusa la regione spinodale, senza le instabilità numeriche spesso associate alla CFD convenzionale in questi regimi. La simulazione riuscita delle interazioni urto-goccia ad alti numeri di Mach funge da validazione rigorosa, dimostrando il potenziale del metodo per investigare flussi complessi nelle tecnologie supercritiche e nei processi guidati da transizioni di fase.