Monolithic kinetic algorithm for heterogeneous porous media systems using a continuous one-domain approach

Sintesi Tecnica: Algoritmo Cinetico Monolitico per Sistemi di Mezzi Porosi Eterogenei

Enunciato del Problema
La simulazione del flusso di fluidi in mezzi porosi eterogenei, in particolare in sistemi che accoppiano regioni di fluido libero con strutture porose, presenta sfide significative. Gli approcci tradizionali si affidano spesso a formulazioni "a due domini", che risolvono equazioni governanti distinte per le regioni di fluido libero e porose e richiedono condizioni al contorno interfacciali esplicite (ad esempio, condizioni di Beavers-Joseph) per accoppiarle. Alternativamente, gli approcci "a un dominio" utilizzano un unico insieme di equazioni, ma storicamente hanno faticato a mantenere la coerenza nei mezzi eterogenei.

Le formulazioni esistenti del Metodo Lattice Boltzmann (LBM) per i mezzi porosi soffrono spesso di una scalatura incoerente della porosità nelle equazioni macroscopiche recuperate. Nello specifico, molti modelli precedenti non riescono a scalare correttamente il gradiente di pressione e i termini di stress viscoso con la porosità quando la porosità varia spazialmente. Inoltre, molte formulazioni esprimono lo stress viscoso di Brinkman in termini di gradienti della velocità intrinseca piuttosto che della velocità superficiale (Darcy), portando a imprecisioni nei domini eterogenei. Inoltre, i modelli LBM standard spesso mancano della flessibilità necessaria per regolare indipendentemente la viscosità bulk effettiva o la velocità del suono, limitando l'efficienza computazionale o la stabilità.

Metodologia
Gli autori propongono un modello Lattice Boltzmann monolitico a dominio singolo basato su un approccio continuo a un dominio. La metodologia si basa sui seguenti pilastri:

1. Modello Macroscopico: L'idrodinamica target è derivata tramite il metodo della media volumetrica, risultando in un modello Darcy-Forchheimer-Brinkman generalizzato. Questo modello è valido per l'intero dominio computazionale, trattando la transizione tra regioni di fluido libero e porose come una variazione continua delle proprietà fisiche (porosità $\epsilon$). Le equazioni governanti includono:

   • Un'equazione di continuità per la densità mediata nel volume ($\epsilon\rho$).
   • Un'equazione del momento che incorpora il trascinamento di Darcy, il trascinamento inerziale di Forchheimer e un tensore di stress viscoso effettivo ($\sigma_{eff}$).
   • Fondamentalmente, lo stress viscoso effettivo è espresso in termini di gradienti della velocità superficiale ($v$), non della velocità intrinseca ($u = v/\epsilon$), garantendo la coerenza con la formulazione di Brinkman per porosità variabile.

2. Solutore Cinetico: Il modello utilizza un'unica equazione cinetica generalizzata su reticoli standard (ad esempio, D2Q9, D3Q19).

   • Distribuzione di Equilibrio: Viene impiegata una distribuzione di equilibrio discreta in forma di prodotto per garantire i momenti corretti e l'invarianza galileiana. Ciò include termini di correzione specifici ($\Phi_{ext}$ e $\Phi_{\eta}$) per ripristinare l'isotropia nel tensore di stress viscoso e consentire la regolazione indipendente della viscosità bulk.
   • Schema di Forzante: Viene utilizzato uno schema di forzante generalizzato coerente (metodo Relax-to-Force) per incorporare il trascinamento Darcy-Forchheimer, le forze di corpo e le necessarie correzioni per i gradienti di porosità.
   • Scalatura della Porosità: I campi localmente conservati sono definiti come densità e quantità di moto mediate nel volume, incorporando naturalmente la scalatura della porosità nei momenti cinetici.
   • Valutazione della Velocità: Viene proposto un algoritmo iterativo per risolvere la relazione implicita tra la velocità macroscopica e i termini di forzante, gestendo il trascinamento non lineare Darcy-Forchheimer e le correzioni viscose.

3. Condizioni al Contorno: Gli autori adattano le standard condizioni al contorno dei nodi bagnati (schema di equilibrio, estrapolazione non-equilibrio e schema di Tamm-Mott-Smith) alla formulazione per mezzi porosi. Una procedura specifica è derivata per calcolare la densità del fluido intrinseco ai bordi utilizzando la continuità della massa, tenendo conto della presenza del termine di forzante e delle variazioni di porosità.

Contributi Chiave

• Formulazione Coerente a Un Dominio: Il documento presenta la prima formulazione LBM che recupera un tensore di stress viscoso di Brinkman con una scalatura della porosità completamente coerente, espressa in termini di velocità superficiale. Ciò consente al modello di gestire porosità spazialmente variabili senza condizioni interfacciali artificiali.
• Correzione dei Termini di Pressione e Viscosi: Gli autori derivano e implementano specifiche correzioni di forzante per rettificare le incoerenze nei termini del gradiente di pressione e dello stress viscoso che derivano dai mezzi eterogenei, garantendo che il modello recuperi il corretto limite macroscopico.
• Regolazione Indipendente dei Parametri: Il modello permette alla velocità del suono e alla viscosità bulk effettiva di essere specificate indipendentemente dalla viscosità di taglio. Ciò migliora l'efficienza computazionale (permettendo velocità del suono più elevate) e la stabilità numerica (regolando la viscosità bulk).
• Algoritmo Monolitico: L'approccio utilizza un unico equazione e solutore cinetico per l'intero dominio, semplificando l'implementazione per geometrie complesse che coinvolgono sia fluidi liberi che corpi porosi.

Risultati e Validazione
Il modello è stato validato attraverso una serie di benchmark e studi numerici:

• Analisi di Coerenza: L'analisi spettrale temporale lineare ha confermato che il modello recupera la corretta velocità del suono effettiva e i tassi di dissipazione per i modi di taglio e acustici. L'analisi ha dimostrato che l'equilibrio in forma di prodotto con le proposte correzioni preserva l'invarianza galileiana e l'isotropia, a differenza delle precedenti formulazioni di equilibrio del secondo ordine.
• Flussi Benchmark:
  • Flusso di Couette: Le simulazioni in un canale poroso omogeneo sotto gravità hanno mostrato un eccellente accordo con le soluzioni numeriche di riferimento sia per i profili di velocità che di densità attraverso vari numeri di Reynolds e di Darcy.
  • Flusso di Poiseuille: Il modello ha catturato accuratamente il flusso in canali omogenei e, crucialmente, in canali con profili di porosità eterogenei. Ha predetto correttamente gli effetti di canalizzazione del flusso dove il fluido scorre preferenzialmente attraverso le regioni ad alta porosità.
• Dimostrazioni Numeriche:
  • Flusso Attorno a un Cilindro Permeabile: Il modello ha simulato con successo il flusso confinato attorno a un cilindro poroso, identificando cinque regimi distinti di flusso basati sul numero di Reynolds e sulla porosità del nucleo. Ha catturato fenomeni quali la stabilizzazione della scia, il distacco della ricircolazione e le instabilità della scia lontana, mostrando un buon accordo qualitativo con la letteratura.
  • Cavità Guidata da Coperchio con Ostacolo Poroso: Il modello ha simulato una cavità con un ostacolo poroso interno, dimostrando come la permeabilità influenzi la formazione di vortici, la penetrazione del flusso e i modelli di ricircolazione. I risultati hanno evidenziato la transizione dal comportamento simile a un solido al comportamento permeabile all'aumentare della porosità.

Significatività e Rivendicazioni
Gli autori affermano che questo lavoro fornisce un quadro robusto e unificato per simulare il trasporto in mezzi porosi eterogenei. Eliminando la necessità di esplicite condizioni al contorno interfacciali e garantendo una scalatura coerente di tutti i termini idrodinamici, il modello è particolarmente adatto per:

• Studi di progettazione su scala di dispositivo in fase iniziale di strutture porose ingegnerizzate con porosità variabile spazialmente.
• Simulare interazioni complesse tra fluidi liberi e corpi permeabili (ad esempio, celle a combustibile, sistemi di filtrazione, combustione porosa).
• Catturare i fenomeni di trasporto in sistemi in cui la transizione tra flusso libero e mezzo poroso è continua o dove le interfacce nette sono difficili da risolvere.

Il documento conclude che, sebbene lo studio attuale si concentri sulla validazione 2D, la formulazione è intrinsecamente multidimensionale e pronta per l'estensione a geometrie 3D. Si identifica come lavoro futuro la necessità di includere il bilancio energetico (trasferimento di calore), miscele multicomponente e flussi multifase per affrontare pienamente il trasporto accoppiato in sistemi ingegneristici realistici.