Upscaling the Navier-Stokes-Cahn-Hilliard model for incompressible multiphase flow in inhomogeneous porous media

Questo lavoro presenta un modello macroscopico rigoroso per il flusso bifase in mezzi porosi disomogenei, derivato mediante omogeneizzazione volumetrica delle equazioni di Navier-Stokes-Cahn-Hilliard, che incorpora il comportamento di bagnatura nel potenziale chimico medio e viene validato attraverso simulazioni numeriche.

L'essenziale

Immagina di dover spiegare come l'acqua e l'olio si muovono attraverso una spugna complessa, come quella che usi in cucina o il terreno sotto i tuoi piedi. Questo è esattamente il cuore del lavoro presentato da Zhang e colleghi.

Ecco una spiegazione semplice, usando analogie quotidiane, di cosa hanno fatto e perché è importante.

1. Il Problema: La Spugna e il "Groviglio"

Immagina di voler capire come due liquidi (ad esempio acqua e petrolio) scorrono attraverso una roccia porosa o un terreno.

• La scala microscopica (Pore-scale): Se guardi attraverso un microscopio potentissimo, vedi che la roccia è piena di buchi minuscoli, curve e angoli. I liquidi devono fare i salti mortali, scivolare lungo le pareti e spingersi a vicenda. È un caos incredibile. Simulare questo per un intero campo petrolifero o una falda acquifera richiederebbe un supercomputer che impiegherebbe anni.
• La scala macroscopica (Darcy-scale): Gli ingegneri, però, non hanno bisogno di sapere cosa succede in ogni singolo buco. Vogliono sapere: "Quanto petrolio uscirà domani da questo pozzo?" o "Quanto velocemente l'inquinante si sposterà?". Hanno bisogno di una visione d'insieme, come guardare la spugna dall'alto e vedere solo il flusso generale.

Il problema è che le regole che funzionano per un singolo buco (le equazioni di Navier-Stokes) non funzionano direttamente se le applichi a tutta la spugna. Bisogna "tradurle".

2. La Soluzione: La "Fotocopia Media" (Upscaling)

Gli autori hanno usato un metodo chiamato Volume Averaging (Media Volumetrica).
Immagina di avere una foto ad altissima risoluzione della spugna piena di liquidi. Invece di guardare ogni singolo pixel, prendi un piccolo quadrato della foto (chiamato REV, o Volume Elementare Rappresentativo) e calcoli la media di tutto ciò che c'è dentro.

• L'analogia del traffico: Immagina di voler descrivere il traffico in una città.
  • Scala microscopica: Sai esattamente dove si trova ogni singola auto, quanto accelera e come cambia corsia.
  • Scala macroscopica (il loro modello): Invece, guardi un quartiere e dici: "In media, qui ci sono 50 auto che viaggiano a 30 km/h".
  • Il trucco è: come si passa dalle regole delle singole auto alle regole del traffico medio senza perdere informazioni importanti?

3. La Novità: Il "Gusto" della Spugna (Bagnabilità)

Qui arriva la parte geniale del loro lavoro.
Nelle vecchie formule, si trattava la spugna come se fosse neutra. Ma in realtà, alcune rocce "amano" l'acqua (sono idrofile) e altre "amano" l'olio (sono idrofobe). Questo si chiama bagnabilità (wettability).

• L'analogia del velcro: Immagina che le pareti dei buchi della spugna abbiano del velcro. Se l'acqua è "appiccicosa" (bagna la roccia), si attacca alle pareti e l'olio deve fare fatica a passare. Se l'olio è appiccicoso, succede il contrario.
• Il contributo degli autori: Hanno creato una formula matematica che include questo "gusto" della roccia direttamente nel calcolo della pressione. Invece di inventare formule a caso basate su esperimenti (come si faceva prima), hanno derivato matematicamente come la bagnabilità influenzi il flusso medio, partendo dalle leggi fisiche di base.

Hanno introdotto un concetto chiamato Potenziale Chimico Medio.

• Spiegazione semplice: È come se dessero un "punteggio di energia" al liquido medio. Se la roccia ama l'acqua, il punteggio dell'acqua cambia, spingendola a muoversi in modo diverso. Questo permette al modello di prevedere come i liquidi si separano o si mescolano senza dover simulare ogni singolo buco.

4. Cosa hanno fatto per dimostrarlo?

Per essere sicuri che la loro "traduzione" funzionasse, hanno usato un metodo di simulazione computerizzata chiamato Lattice Boltzmann (che è come simulare il fluido come se fosse fatto di tanti piccoli pallini che rimbalzano).

Hanno testato il modello in tre scenari:

1. Flusso in un tubo poroso: Hanno verificato che il modello calcolasse correttamente la velocità dell'acqua in una spugna uniforme. Risultato: Perfetto accordo con la teoria.
2. Il confine tra aria e spugna: Hanno simulato l'acqua che passa da un tubo vuoto a una spugna. Il modello ha gestito bene il "salto" di velocità e pressione.
3. Le dita viscose (Viscous Fingering): Questo è il test più spettacolare. Immagina di iniettare acqua in olio molto viscoso. L'acqua non spinge l'olio in modo uniforme, ma crea delle "dita" che si allungano e si dividono (come un fungo).
   • Hanno scoperto che cambiando la bagnabilità (rendendo la spugna più o meno "appiccicosa" per l'acqua), le forme di queste dita cambiavano drasticamente. Il loro modello ha previsto esattamente questo comportamento, mostrando che la "personalità" della roccia è fondamentale.

5. Perché è importante?

Prima di questo lavoro, per prevedere come si muovono i fluidi nel sottosuolo, gli ingegneri dovevano fare molte ipotesi e usare formule "empiriche" (basate su tentativi ed errori).
Ora, grazie a questo studio:

• Abbiamo un ponte matematico solido che collega il mondo microscopico (dove avvengono le cose reali) a quello macroscopico (dove prendiamo le decisioni).
• Possiamo prevedere meglio cose come: quanto petrolio possiamo estrarre, quanto velocemente l'acqua potabile viene contaminata, o come immagazzinare la CO2 nel sottosuolo.
• Il modello è unificato: funziona sia per flussi lenti che veloci, e tiene conto di quanto la roccia "ama" o "odia" i liquidi.

In sintesi:
Questi ricercatori hanno creato un "traduttore" intelligente. Prende le leggi fisiche complesse che governano ogni singolo buco di una spugna e le trasforma in una formula semplice e potente che gli ingegneri possono usare per progettare sistemi reali, tenendo conto anche di quanto la roccia sia "appiccicosa" per i liquidi. È un passo avanti fondamentale per capire il sottosuolo senza dover scavare ogni singolo centimetro.

Sommario tecnico

Titolo: Upscaling del modello Navier-Stokes-Cahn-Hilliard per flussi multifase incomprimibili in mezzi porosi non omogenei

1. Il Problema

Il flusso di fluidi multifase attraverso mezzi porosi è fondamentale in applicazioni ingegneristiche come la bonifica delle acque sotterranee, il recupero migliorato del petrolio (EOR) e lo stoccaggio geologico della CO2.

• Sfida principale: La modellazione diretta (DNS) a scala dei pori è computazionalmente proibitiva per sistemi su larga scala a causa della complessità geometrica e dell'ampio spettro di scale temporali e spaziali.
• Limiti degli approcci attuali: I modelli macroscopici tradizionali si basano sulla legge di Darcy estesa, che richiede relazioni empiriche per la permeabilità relativa e la pressione capillare. Questi modelli spesso falliscono nel catturare fenomeni complessi come l'isteresi, l'effetto dell'angolo di contatto (bagnabilità) e le dinamiche interfacciali, poiché trattano queste grandezze come funzioni empiriche della saturazione senza derivarle rigorosamente dalle equazioni microscopiche.
• Obiettivo: Sviluppare un modello macroscopico rigoroso per flussi bifase che integri le dinamiche a scala dei pori (tensione superficiale, bagnabilità) direttamente nelle equazioni di trasporto a scala di Darcy, superando la necessità di relazioni costitutive puramente empiriche.

2. Metodologia

Gli autori utilizzano il metodo di media volumetrica (Volume Averaging Method - VAM) per derivare le equazioni upscalate partendo dalle equazioni di governo a scala dei pori.

• Equazioni a scala dei pori: Il sistema è governato dalle equazioni di Navier-Stokes per la conservazione della quantità di moto e dall'equazione di Cahn-Hilliard (CH) per l'evoluzione dell'interfaccia tra due fluidi immiscibili e incomprimibili. L'approccio Cahn-Hilliard utilizza un parametro d'ordine ($\phi$) e un funzionale di energia libera che include termini di energia di volume, energia interfacciale ed energia di parete (bagnabilità).
• Procedura di Upscaling:
  1. Definizione del REV: Si definisce un Volume Elementare Rappresentativo (REV) che soddisfa vincoli di scala ($l \ll r_0 \ll L$), dove $l$ è la scala dei pori, $r_0$ la scala del REV e $L$ la scala macroscopica.
  2. Decomposizione di Gray: Le quantità fisiche vengono decomposte in una componente media intrinseca e una deviazione spaziale ($\psi = \langle\psi\rangle_\alpha + \hat{\psi}$).
  3. Teoremi di media: Si applicano i teoremi di media spaziale e temporale alle equazioni di Navier-Stokes e Cahn-Hilliard. Questo processo genera termini di chiusura non noti (derivanti dalle deviazioni spaziali e dalle integrazioni sulle interfacce solido-fluido).
  4. Problemi di chiusura: I termini non chiusi vengono modellati risolvendo problemi locali definiti all'interno del REV. Un contributo chiave è la derivazione formale del potenziale chimico medio, che incorpora direttamente le condizioni al contorno di bagnabilità.
  5. Modello risultante: Si ottiene un modello a "fluido singolo" (single-fluid) a scala di Darcy, dove la forza interfacciale e gli effetti capillari sono espressi attraverso il gradiente del potenziale chimico medio, che dipende dalla bagnabilità e dalla geometria del mezzo.

3. Contributi Chiave

• Derivazione rigorosa del potenziale chimico medio: A differenza dei modelli precedenti che sostituiscono il potenziale chimico con una funzione empirica della pressione capillare, questo lavoro deriva formalmente il potenziale chimico medio includendo esplicitamente la condizione al contorno di bagnabilità (angolo di contatto $\theta$) e l'area superficiale specifica.
• Modello unificato a fluido singolo: Il modello proposto evita la complessità dei modelli a due fluidi separati (con velocità e pressioni distinte), fornendo invece un'unica equazione di quantità di moto che recupera la legge di Darcy, l'equazione di Brinkman e le equazioni di Navier-Stokes a scala dei pori a seconda dei parametri (es. porosità, numero di Darcy).
• Inclusione della bagnabilità nella dinamica macroscopica: La bagnabilità non è un parametro aggiuntivo empirico, ma emerge naturalmente nel termine di forza interfacciale del modello macroscopico attraverso il potenziale chimico.
• Sviluppo di un metodo Lattice Boltzmann (LBM): Gli autori hanno sviluppato un algoritmo LBM specifico per risolvere le nuove equazioni upscalate, permettendo la validazione numerica del modello.

4. Risultati Numerici

Il modello è stato validato attraverso diverse simulazioni numeriche:

• Flusso di Poiseuille in mezzo poroso omogeneo: Confronto con soluzioni analitiche per flussi monofase. I risultati mostrano un eccellente accordo, confermando la capacità del modello di descrivere correttamente la resistenza viscosa e il termine di trascinamento (drag) in mezzi porosi.
• Flusso in canale fluido-poroso: Simulazione di un flusso che attraversa l'interfaccia tra una regione libera e una porosa. Il modello riproduce accuratamente il profilo di velocità e la discontinuità dello stress all'interfaccia, validando la gestione della porosità variabile.
• Problema di Buckley-Leverett: Simulazione di un'iniezione di fluido bagnante in un mezzo saturo di fluido non bagnante. Il modello cattura correttamente la formazione del fronte di shock e la sua velocità, confrontandosi bene con la soluzione analitica classica, nonostante lievi oscillazioni numeriche tipiche degli schemi LBM.
• Dita viscose (Viscous Fingering) in mezzi porosi: Simulazione dell'instabilità viscosa. I risultati dimostrano come la bagnabilità influenzi la morfologia delle dita: angoli di contatto più bassi (superfici idrofile) tendono a disperdere le dita viscose e a inibire la loro formazione, mentre angoli più alti favoriscono la ramificazione. Questo conferma la capacità del modello di catturare la fisica dell'interfaccia.

5. Significato e Conclusioni

Questo lavoro rappresenta un passo significativo nella modellazione dei flussi multifase in mezzi porosi:

• Ponte teorico: Colma il divario tra le descrizioni microscopiche (fisica dettagliata dell'interfaccia) e quelle macroscopiche (leggi di Darcy), fornendo una base teorica solida per i parametri di trasporto.
• Superamento dell'empirismo: Riduce la dipendenza da correlazioni empiriche per la pressione capillare e la permeabilità relativa, derivandole invece da principi termodinamici e condizioni al contorno fisiche.
• Versatilità: Il modello è in grado di recuperare diverse formulazioni note (Darcy, Brinkman, Stokes) come casi limite, offrendo un framework unificato per studiare fenomeni complessi in mezzi porosi eterogenei.
• Implicazioni future: Sebbene alcuni parametri di chiusura siano stati semplificati o stimati, il framework proposto offre una strada promettente per migliorare la previsione dei processi di trasporto multifase in settori come l'energia e l'ambiente, riducendo l'incertezza legata alla parametrizzazione empirica.