Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Sintesi Tecnica: Metodi degli elementi virtuali funzione di corrente – pressione per il problema dell'interfaccia Stokes–Darcy

Enunciato del Problema
L'articolo affronta la simulazione numerica del problema accoppiato di interfaccia Stokes–Darcy, una sfida fondamentale nella meccanica dei fluidi che governa l'interazione tra regioni a flusso libero (Stokes) e mezzi porosi (Darcy). L'accoppiamento è ottenuto tramite la continuità delle velocità normali (conservazione della massa), l'equilibrio delle tensioni normali e la condizione di scivolamento tangenziale di Beavers–Joseph–Saffman. Sebbene esistano vari metodi a Elementi Finiti (FE), Discontinuous Galerkin (dG) e Hybrid High–Order (HHO), che spesso si basano su coppie velocità–pressione inf–sup stabili, questo lavoro si concentra su una specifica riformulazione in domini semplicemente connessi bidimensionali. In questo contesto, il vincolo di incompressibilità è soddisfatto identicamente introducendo una funzione di corrente scalare, eliminando così la variabile pressione dalle equazioni di Stokes e riducendo il sistema a un problema biarmonico del quarto ordine accoppiato con un'equazione della pressione di Darcy del secondo ordine.

Metodologia
Gli autori propongono una nuova discretizzazione mediante il Metodo degli Elementi Virtuali (VEM) per questa formulazione funzione di corrente–pressione. La metodologia sfrutta la capacità del VEM di gestire mesh poligonali generiche, il che è particolarmente vantaggioso per geometrie di interfaccia complesse e mesh non corrispondenti (non-matching).

1. Formulazione: Il problema continuo è espresso in forma primale-primale. Il dominio di Stokes utilizza una funzione di corrente $\chi \in H^2$, mentre il dominio di Darcy utilizza un potenziale di pressione $\phi \in H^1$. La formulazione debole risulta in un problema a punto di sella perturbato.
2. Discretizzazione:
   • Dominio di Stokes: Viene impiegato uno spazio VEM di ordine minimo conforme $C^1$. Questo comporta l'uso di polinomi quadratici ($P_2$) e richiede la conformità $H^2$, ottenuta attraverso spazi virtuali potenziati con gradi di libertà (DoF) definiti da valori nodali e gradienti nodali.
   • Dominio di Darcy: Viene utilizzato uno spazio VEM di ordine minimo conforme $C^0$ per la pressione, basato su polinomi lineari ($P_1$).
   • Accoppiamento: Le condizioni di interfaccia sono imposte naturalmente all'interno della formulazione debole. Il metodo accoglie mesh non corrispondenti all'interfaccia senza necessità di rimesh.
3. Implementazione: Lo schema è implementato utilizzando la libreria open-source orientata agli oggetti vem++. Una caratteristica specifica, la classe vem::vemMesh2dDofFilter, viene utilizzata per gestire i DoF attraverso i diversi sottodomini, permettendo la definizione simultanea di diversi spazi VEM su una singola mesh e il filtraggio degli stessi in base ai marcatori degli elementi.
4. Stabilizzazione: Le forme bilineari discrete includono termini di stabilizzazione ($S_S^K$ e $S_D^K$) per garantire la coercitività sul nucleo (kernel) degli operatori di proiezione, seguendo le pratiche standard VEM (ad esempio, ricetta diagonale e DOFI–DOFI).

Contributi Chiave
L'articolo delinea cinque contributi primari:

• Analisi Continua: Sviluppo di un'analisi continua per il problema dell'interfaccia Stokes–Darcy in formulazione funzione di corrente–pressione, che evita intrinsecamente la necessità di condizioni di stabilità inf–sup.
• Schema Discreto: Progettazione di una discretizzazione VEM $C^1$–$C^0$ di ordine minimo. Gli autori evidenziano l'efficienza computazionale di questo schema, notando che esso coinvolge solo $3N_S + N_D + 4N_\Sigma$ gradi di libertà (dove $N_S$, $N_D$ e $N_\Sigma$ sono i conteggi dei vertici interni nel dominio di Stokes, nel dominio di Darcy e sull'interfaccia, rispettivamente).
• Analisi dell'Errore: Derivazione di stime dell'errore a priori per la discretizzazione proposta, provando tassi di convergenza ottimali sotto ipotesi di regolarità standard.
• Implementazione Software: Integrazione di nuove funzionalità nella libreria vem++ per facilitare l'implementazione diretta di problemi di interfaccia tramite filtraggio della mesh.
• Validazione Numerica: Esecuzione di esperimenti numerici per confermare i risultati teorici.

Risultati
Gli esperimenti numerici validano i tassi di convergenza teorici e la robustezza del metodo:

• Soluzioni Costruite: Due esperimenti con soluzioni costruite (uno con valori di interfaccia nulli, uno con valori non nulli) su vari tipi di mesh (quadrilaterale, non convessa, Voronoi e Voronoi perturbata) dimostrano un decadimento asintotico lineare dell'errore ($O(h)$), coerente con le previsioni teoriche per lo schema di ordine minimo.
• Filtrazione a Fondo Cieco (Dead-End Filtration): Una simulazione di un filtro a fondo cieco a forma di quarto di anello recupera con successo i risultati di studi precedenti, dimostrando la capacità del metodo di gestire interfacce curve e condizioni al contorno miste.
• Organi Bioartificiali: Una simulazione del flusso sanguigno in un organo bioartificiale (basato su incapsulamento) con costanti di permeabilità variabili illustra la capacità del metodo di modellare flussi fisiologici complessi e l'influenza della permeabilità sul flusso attraverso l'interfaccia.

Significatività e Rivendicazioni
Il lavoro si posiziona come un ponte per colmare una lacuna nella letteratura riguardante l'estensione delle formulazioni della funzione di corrente ai problemi accoppiati Stokes–Darcy utilizzando gli Elementi Virtuali. Gli autori affermano che la struttura proposta offre una "strategia completamente priva di inf–sup", semplificando la progettazione degli spazi discreti rispetto agli approoli tradizionali velocità–pressione. La capacità di gestire geometrie di interfaccia irregolari e mesh non corrispondenti senza rimesh è evidenziata come un vantaggio chiave del framework VEM.

Gli autori rimangono modesti riguardo alle sfide future, notando esplicitamente che l'estensione di questa formulazione della funzione di corrente a domini multiplamente connessi (dove la funzione di corrente non è determinata univocamente da condizioni al contorno standard) e la scalabilità del framework a tre dimensioni (che richiederebbe potenziali vettoriali e spazi $C^1$ di ordine superiore per funzioni vettoriali) rimangono sfide aperte. Inoltre, l'integrazione di tecniche prive di stabilizzazione e la derivazione di stimatori dell'errore a posteriori robusti sono identificati come direzioni per la ricerca futura.