A Shakhov-based Bhatnagar-Gross-Krook model for polyatomic molecules and for atomic as well as polyatomic mixtures

Questo articolo presenta l'estensione del modello Shakhov-Bhatnagar-Gross-Krook (SBGK) nel codice particellare PICLas per simulare molecole poliatomiche e miscele di atomi e molecole in condizioni di non equilibrio, convalidandone l'accuratezza attraverso confronti con i metodi DSMC ed ESBGK in flussi supersonici e ipersonici.

L'essenziale

Immagina di dover prevedere come si comporta l'aria quando un razzo attraversa l'atmosfera a velocità incredibili, o come il gas si muove in un microchip minuscolo. In questi scenari, l'aria non si comporta come un fluido semplice e uniforme (come l'acqua in un fiume calmo), ma diventa un caos di singole molecole che rimbalzano, vibrano e ruotano in modo disordinato.

Questo articolo scientifico parla di un nuovo strumento matematico (un modello chiamato SBGK) creato per simulare questi comportamenti complessi in modo più veloce e preciso rispetto ai metodi attuali.

Ecco una spiegazione semplice, usando delle metafore:

1. Il Problema: Troppi Calcoli, Troppo Lento

Per capire come si muove un gas, gli scienziati usano due approcci principali:

• DSMC (Il metodo "Contare le formiche"): È come guardare ogni singola molecola di gas come se fosse una formica. È super preciso, ma se hai un miliardo di formiche, ci vuole un'eternità per calcolare dove vanno. È troppo lento per i computer attuali quando il gas è denso.
• CFD (Il metodo "Flusso d'acqua"): Tratta il gas come un fluido continuo, ignorando le singole molecole. È velocissimo, ma fallisce quando le molecole sono troppo distanti o si comportano in modo strano (come vicino a un razzo che viaggia a velocità supersonica).

L'obiettivo di questo studio è creare un ponte tra questi due mondi: un metodo che sia veloce come il fluido continuo, ma preciso come il conteggio delle singole particelle.

2. La Soluzione: Il Modello SBGK (Il "Direttore d'Orchestra")

Gli autori hanno migliorato un modello chiamato Shakhov-BGK.
Immagina il gas come un'orchestra caotica dove ogni musicista (molecola) suona a modo suo.

• I vecchi modelli (come il BGK standard) dicevano: "Ok, tutti suonate la stessa nota media". Questo funziona bene per la melodia generale, ma sbaglia i dettagli (come il calore che si sposta).
• Il modello Shakhov è come un direttore d'orchestra più intelligente. Non solo dice "suonate questa nota", ma corregge anche come l'energia (il calore) si sposta tra i musicisti. Questo permette di calcolare la viscosità (quanto il gas è "appiccicoso") e la conduttività termica (quanto bene trasmette calore) in modo perfetto.

3. La Novità: Mischia di Molecole e Atomi

Fino a poco tempo fa, questi modelli funzionavano bene solo per gas semplici (come l'ossigeno puro). Ma nello spazio o nei motori, spesso abbiamo miscele:

• Atomi: Come l'Argon o l'Elio (sono come palline lisce che rotolano).
• Molecole: Come l'Azoto o la CO2 (sono come palline con delle "braccia" che possono ruotare e vibrare).

Il problema è che le molecole hanno energie interne (ruotano e vibrano come una trottola o una molla). Quando urtano, possono scambiare energia tra il movimento in avanti (traslazione) e la rotazione/vibrazione.
Questo nuovo modello è stato aggiornato per gestire miscele complesse:

• Può gestire un gas fatto solo di atomi.
• Può gestire un gas fatto solo di molecole complesse.
• Può gestire un mix caotico di entrambi (come l'aria che entra in un'atmosfera marziana fatta di CO2 e Azoto).

4. Come l'hanno Testato? (Le Prove sul Campo)

Per vedere se il loro "direttore d'orchestra" funziona davvero, hanno simulato due scenari estremi:

• Scenario A: Il Flusso tra due pareti (Couette Flow).
  Immagina due lastre di metallo che si muovono in direzioni opposte con del gas in mezzo. Hanno provato con gas puri e miscele diverse.

  • Risultato: Il loro modello ha prodotto risultati quasi identici al metodo "Contare le formiche" (DSMC), ma molto più velocemente. C'è stata una piccola differenza solo quando hanno mischiato atomi molto pesanti con atomi molto leggeri (come Argon ed Elio), ma è un limite noto di questi calcoli.

• Scenario B: Il Cono Blunt (Il Razzo).
  Hanno simulato un flusso ipersonico (super veloce) contro un cono smussato, come la punta di un razzo che rientra nell'atmosfera. Qui si formano onde d'urto (shock), zone dove l'aria viene compressa e riscaldata all'improvviso.

  • Risultato: Questo è il punto forte. Il modello SBGK ha disegnato l'onda d'urto in modo molto più preciso rispetto al modello concorrente (ESBGK). È come se il nuovo modello vedesse il "muro" dell'aria compressa con più nitidezza, mentre il vecchio modello lo vedeva un po' sfocato.

5. Perché è Importante?

Immagina di voler progettare un satellite che deve orbitare bassissimo sopra la Terra o un veicolo che atterra su Marte.

• L'atmosfera è sottile e le molecole sono distanti (regime di "flusso libero").
• Ma vicino al veicolo, l'aria è densa e calda (regime "continuo").

I metodi attuali faticano a passare da un regime all'altro senza fare errori o impazzire nei calcoli. Questo nuovo modello SBGK è pronto per essere usato in combinazione con altri metodi per simulare questi viaggi spaziali complessi, risparmiando tempo di calcolo e fornendo dati più precisi su quanto il veicolo si scalda.

In sintesi: Gli scienziati hanno creato un "simulatore di gas" più intelligente che sa gestire miscele di atomi e molecole complesse, correggendo gli errori dei modelli precedenti specialmente quando si formano onde d'urto violente. È un passo avanti fondamentale per l'esplorazione spaziale e la tecnologia dei micro-dispositivi.

Sommario tecnico

Sintesi Tecnica: Modello Shakhov-BGK per Miscele di Molecole Poliatomiche

Titolo: A Shakhov-based Bhatnagar-Gross-Krook model for polyatomic molecules and for atomic as well as polyatomic mixtures
Autori: M. Pfeiffer e F. Tuttas (Istituto di Sistemi Aerospaziali, Università di Stoccarda)

---

1. Il Problema

Le simulazioni numeriche della fluidodinamica in applicazioni spaziali, microscopiche e nanoscopiche presentano sfide significative a causa delle grandi gradienti di densità che spaziano dal regime continuo a quello di flusso molecolare libero.

• Limiti dei metodi attuali: Il metodo Direct Simulation Monte Carlo (DSMC) è altamente accurato per flussi non in equilibrio e multi-scala, ma il suo costo computazionale diventa proibitivo nei regimi di transizione e continuo. D'altra parte, la Fluidodinamica Computazionale (CFD) classica è efficiente nel regime continuo ma fatica a gestire le interfacce con il DSMC a causa del rumore statistico e delle differenze fondamentali negli approcci.
• Necessità di modelli intermedi: I metodi basati su particelle per il regime continuo, come i modelli Bhatnagar-Gross-Krook (BGK), offrono un'alternativa promettente riducendo i costi computazionali. Tuttavia, i modelli BGK standard spesso non riescono a riprodurre correttamente il numero di Prandtl (fondamentale per la conduzione termica) e le versioni esistenti per miscele di gas (atomiche e molecolari) con gradi di libertà interni (rotazionali e vibrazionali) non in equilibrio sono complesse o computazionalmente onerose (ad esempio, modelli con multipli termini di rilassamento).

2. Metodologia

Gli autori hanno esteso l'implementazione del modello Shakhov-BGK (SBGK) nel codice open-source basato su particelle PICLas. L'approccio si basa su un singolo termine di rilassamento per specie, esteso per gestire:

• Molecole poliatomiche: Inclusione dei gradi di libertà interni (rotazionali e vibrazionali) quantizzati.
• Miscele di gas: Gestione di miscele contenenti atomi, molecole biatomiche e poliatomiche.

Aspetti chiave del modello:

• Funzione di distribuzione target: Per gli atomi, la distribuzione target è puramente traslazionale. Per le molecole, è fattorizzata in componenti traslazionali, rotazionali e vibrazionali, permettendo di modellare stati quantizzati vibrazionali (oscillatore armonico semplice).
• Correzione di Shakhov: La funzione di distribuzione target è modificata per includere un flusso di calore correttivo, garantendo la riproduzione del numero di Prandtl fisico corretto ($Pr_{phys}$).
• Proprietà di trasporto: Per determinare la viscosità e la conducibilità termica delle miscele, è stato impiegato il metodo di prima approssimazione basato sugli integrali di collisione (metodo di Chapman-Enskog), utilizzando il modello della sfera rigida variabile (VHS). Questo approccio è stato scelto per la sua maggiore accuratezza rispetto alla regola di miscelazione di Wilke per i casi studiati.
• Parametri del modello: È stato derivato un parametro specifico ($\alpha$) che dipende dalla frazione di massa, dalla densità e dai gradi di libertà interni, permettendo di calcolare il parametro di Shakhov ($Pr_S$) necessario per ottenere il numero di Prandtl fisico desiderato.
• Conservazione: Il modello soddisfa rigorosamente i principi di conservazione di massa, quantità di moto ed energia.

3. Contributi Chiave

1. Estensione del modello SBGK: Sviluppo di un modello unificato capace di simulare simultaneamente atomi e molecole (biatomiche e poliatomiche) in miscele, gestendo l'equilibrio non termico nei gradi di libertà interni.
2. Derivazione analitica del parametro di Prandtl: Formulazione esplicita per calcolare il parametro di rilassamento Shakhov in funzione delle proprietà della miscela, garantendo la correttezza fisica del trasporto di calore.
3. Efficienza computazionale: Utilizzo di un singolo termine di rilassamento per specie, che riduce drasticamente il costo computazionale rispetto ai modelli a multipli termini di rilassamento, mantenendo un numero ridotto di particelle di simulazione necessarie rispetto al DSMC.
4. Implementazione in PICLas: Integrazione completa nel codice open-source, inclusa la tecnica di campionamento acceptance-rejection per le nuove stati traslazionali.

4. Risultati

Il modello è stato validato attraverso due casi di test principali, confrontando i risultati con il DSMC (riferimento) e il modello Ellipsoidal Statistical BGK (ESBGK):

• Flusso di Couette Supersonico:

  • Testati tre casi: $N_2$ puro, miscela Ar-He, e miscela $N_2$-N.
  • Risultati: L'accordo tra SBGK e DSMC è eccellente per le temperature traslazionali e rotazionali. Per la miscela Ar-He (alto rapporto di massa), si osserva una leggera deviazione, attribuita alla difficoltà di approssimare il numero di Prandtl efficace, ma il modello rimane robusto. La miscela $N_2$-N mostra un accordo quasi perfetto.

• Flusso Ipersonico su un Cono Smussato a 70°:

  • Testati tre casi: $N_2$ puro, miscela $CO_2$-$N_2$ (simulazione ingresso su Marte), e miscela ternaria $N_2$-$O_2$-N.
  • Confronto con ESBGK: Il modello SBGK proposto supera il modello ESBGK nella cattura della struttura dell'onda d'urto. Mentre l'ESBGK tende a prevedere un inizio dell'aumento di temperatura traslazionale troppo anticipato (allargando artificialmente la regione d'urto), il SBGK riproduce con precisione lo spessore e l'inizio dell'urto, in ottimo accordo con il DSMC.
  • Flusso di calore: L'accordo sui flussi di calore sulla superficie del cono è eccellente per tutti i metodi e corrisponde anche ai dati sperimentali disponibili.

5. Significato e Impatto

Questo lavoro rappresenta un passo significativo verso la creazione di strumenti numerici ibridi efficienti per la simulazione di flussi aerospaziali complessi.

• Accuratezza nell'Urto: La capacità del modello SBGK di catturare più precisamente la regione d'urto rispetto all'ESBGK è cruciale per la progettazione di veicoli di rientro atmosferico, dove il trasferimento di calore e le forze aerodinamiche sono determinati dalla struttura dell'urto.
• Versatilità: La capacità di gestire miscele complesse con gradi di libertà interni non in equilibrio rende il modello adatto per scenari reali come l'ingresso nell'atmosfera marziana (ricco di $CO_2$) o terrestre ad alta velocità.
• Efficienza: Offrendo un compromesso ottimale tra accuratezza fisica (numero di Prandtl corretto, gestione dei gradi interni) e costo computazionale, il modello SBGK esteso facilita l'accoppiamento futuro con il DSMC per risolvere problemi multi-scala in modo più efficiente rispetto all'uso esclusivo del DSMC.

In conclusione, l'estensione del modello Shakhov-BGK a miscele poliatomiche fornisce un metodo robusto e preciso per la fluidodinamica del gas rarefatto, superando le limitazioni dei modelli BGK standard e offrendo prestazioni superiori rispetto alle varianti ellissoidali (ESBGK) in scenari di shock forte.