Multicontinuum Generalized Multiscale Finite Element Method (MC-GMsFEM). Theory and applications to upscaling of two-phase flow

Sintesi Tecnica: Metodo degli Elementi Finiti Multiscala Generalizzato Multicontinuo (MC-GMsFEM)

Problema
I problemi multiscala in mezzi eterogenei vengono tipicamente affrontati costruendo modelli a scala grossolana che incorporano informazioni a scala fine tramite funzioni di base multiscala. Sebbene i metodi esistenti (come il GMsFEM) approssimino efficientemente le soluzioni a scala fine su griglie grossolane, essi generalmente mancano di un quadro sistematico per derivare equazioni differenziali macroscopiche fisicamente coerenti. La rappresentazione a scala grossolana in questi metodi coinvolge spesso gradi di libertà che mancano di continuità spaziale, oscurando la connessione con i modelli di continuo sottostanti. Inoltre, le tecniche di omogeneizzazione multicontinuo tradizionali faticano ad affrontare sistemi complessi in cui il numero di continui interagenti varia spazialmente e dinamicamente, richiedendo spesso la risoluzione di problemi di cella speciali che sono difficili da implementare.

Questo articolo affronta tali limitazioni sviluppando un quadro unificato capace di costruire funzioni di base multiscala conformi e, contemporaneamente, consentendo la derivazione sistematica delle equazioni governanti macroscopiche. Un'applicazione primaria considerata è l'upscaling del flusso bifase immiscibile in mezzi porosi eterogenei, dove sorgono molteplici continui interagenti (inclusi fasi mobili e intrappolate).

Metodologia: Framework MC-GMsFEM
Il metodo proposto (Multicontinuum Generalized Multiscale Finite Element Method - MC-GMsFEM) integra le costruzioni multiscala degli elementi finiti con le rappresentazioni multicontinuo. La metodologia centrale procede come segue:

1. Identificazione del Continuo: Il metodo assume l'identificazione dei continui globalmente (ad esempio, tramite funzioni caratteristiche $\psi_i$). Nel contesto del flusso bifase, questi continui corrispondono a distinti fluidi o regioni specifiche come i gangli intrappolati.
2. Costruzione delle Funzioni Ausiliarie:
   • Funzioni di Continuo ($N_i$): Le funzioni locali che rappresentano l' $i$-esimo continuo sono costruite risolvendo problemi locali con vincoli integrali. Nello specifico, $N_i$ è definito in modo tale da avere una media unitaria sul continuo $i$ e medie nulle su tutti gli altri continui all'interno di un blocco grossolano $K$: $\int_K N_i \psi_j = \delta_{ij} \int_K \psi_j$.
   • Funzioni di Correzione del Gradiente ($M_i$): Per tenere conto delle caratteristiche a scala fine e dei gradienti, vengono costruite funzioni ausiliarie $M_i$ (o i loro corrispettivi discreti $z_i = M_i \cdot \nabla \phi_0$). Queste soddisfano problemi locali con condizioni al contorno nulle sulle frontiere dei blocchi grossolani e vincoli di ortogonalità rispetto ai continui.
3. Rappresentazione della Funzione di Base: Le funzioni di base multiscala sono espresse come:
   $$\phi^j_i = N_i \phi^j_0 + M_i \cdot \nabla \phi^j_0$$
   dove $\phi^j_0$ sono funzioni di partizione dell'unità standard. Questa rappresentazione separa le variabili grossolane e i loro gradienti.
4. Derivazione delle Equazioni Macroscopiche: La soluzione a scala grossolana è decomposta come $u_H = \sum_i (N_i U_i + M_i \cdot \nabla U_i)$, dove $U_i$ sono le variabili del continuo. Sostituendo questa decomposizione nella formulazione variazionale delle equazioni governanti (ad esempio, equazioni ellittiche o flusso di Stokes), gli autori derivano le equazioni macroscopiche. Queste equazioni includono naturalmente termini che coinvolgono i gradienti delle variabili del continuo (termini di tipo Brinkman), che sono cruciali per catturare gli effetti nelle regioni intrappolate.
5. Applicazione al Flusso Bifase: Il framework è applicato a equazioni accoppiate di flusso di Stokes e trasporto. Il metodo tratta le lunghe regioni di fluido intrappolato come continui distinti. Le equazioni di incompressibilità e di quantità di moto sono derivate sia per formulazioni di pressione a singola pressione che a multicontinuo, risultando in modelli macroscopici che tengono conto dei termini di trasferimento tra i continui e dei vincoli di incompressibilità.

Risultati Chiave
L'articolo presenta esperimenti numerici per equazioni ellittiche in domini eterogenei e perforati, nonché simulazioni di flusso bifase in mezzi perforati.

• Problemi Ellittici: L'MC-GMsFEM è stato testato su microstrutture eterogenee e perforate. Il metodo ha dimostrato un'elevata accuratezza, con errori relativi per le soluzioni upscaled e downscaled che generalmente rimanevano al di sotto del 20% in vari casi di test (Tabella 1). Le visualizzazioni hanno confermato che le soluzioni multiscala downscaled corrispondevano strettamente alle soluzioni di riferimento a griglia fine.
• Flusso Bifase:
  • Rilevamento dei Gangli: Un algoritmo di ricerca in ampiezza (BFS) è stato utilizzato per identificare le regioni di fluido intrappolato (gangli) e trattarle come continui separati.
  • Accuratezza: Nei test riguardanti i Target domini 1 e 2 con varie realizzazioni di "finger" (Test A–E), l'approccio MC-GMsFEM (sia nelle varianti globali che locali) ha prodotto soluzioni a scala grossolana con piccoli errori relativi rispetto alle medie di riferimento a griglia fine (Tabelle 2–16).
  • Impatto dei Termini di Gradiente: Un esperimento critico (Test C) ha confrontato i modelli con e senza i termini di gradiente di velocità ($M_i \cdot \nabla U_i$). I risultati (Tabella 17) hanno mostrato che l'omissione di questi termini (riducendo il modello a un modello standard di Darcy multicontinuo) aumentava significativamente gli errori, in particolare per le velocità dei gangli intrappolati (l'errore è passato da circa il 16% al 96%). Ciò indica che i termini di gradiente sono essenziali per modellare accuratamente le regioni intrappolate.
  • Separazione del Continuo: Lo studio ha mostrato che trattare i gangli come continui distinti migliora la coerenza fisica del modello, sebbene l'impatto numerico sugli errori di concentrazione sia stato trovato trascurabile nei casi di test specifici, mentre la rappresentazione della velocità è stata significativamente raffinata.

Significato e Rivendicazioni
L'articolo sostiene che l'MC-GMsFEM offra un approccio robusto e generale per la modellazione di sistemi multiscala complessi. La sua importanza primaria risiede in:

1. Quadro Unificato: Fornisce un percorso sistematico per derivare le equazioni governanti macroscopiche direttamente dalle discretizzazioni multiscala, colmando il divario tra approssimazione numerica e omogeneizzazione del continuo.
2. Flessibilità: A differenza dei precedenti approcci multicontinuo, questo framework gestisce mezzi eterogenei in cui il numero di continui interagenti varia spazialmente e dinamicamente senza richiedere la soluzione di problemi di cella speciali per ogni configurazione.
3. Coerenza Fisica: Integrando naturalmente il metodo in un contesto standard di elementi finiti, esso permette la derivazione di equazioni macroscopiche che includono i necessari termini di gradiente, i quali si sono dimostrati critici per catturare la dinamica delle fasi intrappolate nel flusso bifase.
4. Efficienza Computazionale: Il metodo consente simulazioni efficienti su griglia grossolana preservando al contempo le informazioni a scala fine, offrendo un'alternativa computazionalmente efficiente alla modellazione completa a scala fine.

Gli autori notano con modestia che, sebbene il metodo offra vantaggi rispetto alle tecniche di omogeneizzazione esistenti, non definisce esplicitamente i coefficienti del termine di gradiente nelle equazioni macroscopiche in una forma chiusa che possa essere facilmente pre-calcolata senza risolvere problemi locali; tuttavia, l'adattabilità del metodo a continui variabili e generiche eterogeneità è evidenziata come il suo principale punto di forza.