A Divergence-Free Scott-Vogelius Finite Element Method for the Surface Stokes Problem

Sintesi Tecnica: Un Metodo degli Elementi Finiti Scott-Vogelius Divergenza-Free per il Problema di Stokes Superficiale

Enunciato del Problema
Il documento affronta il problema stazionario di Stokes superficiale su una superficie liscia, chiusa e orientabile bidimensionale $\gamma \subset \mathbb{R}^3$. Le equazioni governanti coinvolgono un campo di velocità tangenziale $u$ e una pressione $p$, soggetti ai vincoli che $u$ sia tangenziale alla superficie ($u \cdot \nu = 0$) e che il flusso sia incomprimibile ($\text{div}_\gamma u = 0$).

La costruzione di metodi a elementi finiti per questo problema presenta due sfide primarie:

1. Tangenzialità: L'imposizione forte del vincolo tangenziale su una geometria discreta porta spesso a fenomeni di "locking" quando si utilizzano spazi vettoriali standard conformi $H^1$, degradando le proprietà di approssimazione.
2. Incomprimibilità: Le coppie velocità-pressione standard (come Taylor-Hood) impongono solo debolmente il vincolo di divergenza nulla. Sebbene esistano metodi che impongono fortemente la tangenzialità (ad esempio, spazi BDM conformi $H(\text{div})$), essi spesso risultano in sistemi algebrici di grandi dimensioni o mancano di un'analisi di stabilità completa sulle superfici. Al contrario, i metodi che impongono fortemente la tangenzialità mantenendo la conformità $H^1$ (ad esempio, MINI o Taylor-Hood su superfici) tipicamente impongono l'incomprimibilità solo debolmente.

Gli autori mirano a costruire un metodo che imponga simultaneamente e con esattezza sia i vincoli di tangenzialità che di incomprimibilità, mantenga la dimensionalità delle standard discretizzazioni euclidee 2D ed eviti le complessità di implementazione dei costrutti a macro-elemento su superfici curve.

Metodologia
L'approccio proposto estende il framework degli elementi finiti superficiali di Kone et al. [10, 11] alla coppia Scott-Vogelius. La metodologia comprende i seguenti componenti chiave:

• Triangolazione: Il metodo utilizza una raffinazione Clough-Tocher (CT) di un'approssimazione poliedrica della superficie. Ogni triangolo della mesh iniziale viene suddiviso collegando il suo baricentro ai suoi vertici. Questo raffinamento rilassa i requisiti del grado polinomiale per la stabilità (permettendo $r \ge 2$) ed elimina i vertici "quasi singolari".
• Spazi degli Elementi Finiti:
  • Velocità ($V_h$): Definita come polinomi a pezzi continui di grado $r$ sulla triangolazione superficiale raffinata tramite CT. Per imporre la continuità debole attraverso gli spigli degli elementi (necessaria per la conformità $H(\text{div})$ sulla superficie), gli autori impiegano una trasformazione di Piola ausiliaria ai gradi di libertà (DOF) Lagrangiani. Questo unisce gli spazi locali senza richiedere termini di integrazione inter-elemento o di stabilizzazione.
  • Pressione ($Q_h$): Definita come polinomi a pezzi discontinui di grado $r-1$ sulla stessa mesh raffinata tramite CT.
• Strategia di Implementazione: A differenza dei precedenti metodi Scott-Vogelius isoparametrici basati su definizioni di macro-elemento (dove le funzioni di base sono definite su un elemento di riferimento suddiviso e mappate globalmente), questo articolo definisce gli spazi degli elementi finiti direttamente sulla triangolazione superficiale raffinata. Questo approccio evita le complesse strutture dati e le regole di quadratura specializzate richieste dalle implementazioni a macro-elemento, portando a uno schema più pratico e compatibile con i software standard per elementi finiti.
• Analisi della Stabilità: Gli autori dimostrano la condizione di inf-sup discreta (condizione di Ladyzhenskaya-Babuška-Brezzi). La prova utilizza un argomento di perturbazione:
  1. La stabilità è prima stabilita per una coppia "macro-elemento" ausiliaria definita su una triangolazione di ordine superiore leggermente diversa.
  2. La discrepanza tra lo spazio di raffinamento diretto proposto e lo spazio del macro-elemento ausiliario viene quantificata come una perturbazione $O(h)$.
  3. La stabilità viene trasferita dalla coppia ausiliaria alla coppia proposta tramite questa analisi di perturbazione.

Contributi Chiave e Risultati

1. Discretizzazione Esatta di Divergenza Nulla e Tangenzialità: La coppia Scott-Vogelius costruita su triangolazioni Clough-Tocher curve impone con esattezza sia la tangenzialità ($u \cdot \nu = 0$) che l'incomprimibilità ($\text{div}_\gamma u = 0$).
2. Implementazione Semplificata: Definendo gli spazi direttamente sulla triangolazione raffinata anziché tramite macro-elementi, gli autori forniscono uno schema significativamente più semplice da implementare pur mantenendo i benefici teorici della coppia Scott-Vogelius.
3. Stabilità e Convergenza:
   • Il metodo è dimostrato essere inf-sup stabile con una costante indipendente dalla dimensione della mesh.
   • Si deriva una convergenza di ordine ottimale nel regime isoparametrico ($k=r$). Nello specifico, la velocità converge con ordine $O(h^r)$ nella norma $H^1$ e $O(h^{r+1})$ nella norma $L^2$, mentre la pressione converge con ordine $O(h^r)$ nella norma $L^2$.
   • L'analisi include una prova di stabilità per la corrispondente coppia macro Scott-Vogelius sulle superfici e per la coppia superficiale $P_r - P_0$ come sottoprodotti.
4. Analisi della Robustezza della Pressione: Il documento investiga la robustezza della pressione (dove l'errore di velocità è indipendente dall'approssimazione della pressione). Gli autori sostengono che, per i problemi di Stokes superficiali, la robustezza della pressione non è automaticamente garantita anche per i metodi con divergenza nulla esatta, a causa degli errori di approssimazione geometrica e dell'estensione del termine sorgente. Forniscono un framework per costruire estensioni della sorgente che preservino la robustezza della pressione e dimostrano, tramite esperimenti numerici, che le estensioni standard possono distruggere tale proprietà.

Esperimenti Numerici
Gli autori validano la teoria attraverso tre esperimenti numerici sulla sfera unitaria e sul toro:

• Test 1 (Flusso a divergenza nulla): Conferma i tassi di convergenza ottimali ($O(h^{r+1})$ per la $L^2$ della velocità, $O(h^r)$ per il gradiente della velocità e la pressione) e mostra che la soluzione è priva di divergenza fino alla precisione di macchina.
• Test 2 (Nessun flusso/Velocità nulla): Dimostra che, sebbene la pressione converga in modo ottimale, l'approssimazione della velocità è non nulla a causa degli errori geometrici, evidenziando la mancanza di robustezza automatica della pressione. L'esperimento osserva anche una superconvergenza nell'errore $L^2$ della velocità ($O(h^{r+2})$), suggerendo che la stima teorica $L^2$ potrebbe non essere stringente.
• Test 3 (Toro): Conferma i tassi di convergenza ottimali su una geometria più complessa, sebbene con un certo degrado attribuito al condizionamento del sistema e agli errori di arrotondamento.

Significato
Il lavoro sostiene di colmare una lacuna nella letteratura riguardante la dinamica dei fluidi superficiali, fornendo un metodo che combina i vantaggi delle discretizzazioni con divergenza nulla esatta (leggi di conservazione, potenziale robustezza della pressione) con la praticità delle implementazioni standard degli elementi finiti. Allontanandosi dalla costruzione a macro-elemento, gli autori offrono uno schema "sostanzialmente più semplice e pratico da implementare" che mantiene il rigore teorico della coppia Scott-Vogelius su superfici curve. Il lavoro estende l'applicabilità dei metodi a divergenza nulla ai problemi superficiali, che sono critici per modellare fenomeni come i bilayer lipidici e i film sottili.