Stabilization-Free H(curl) and H(div)-Conforming Virtual Element Method

Sintesi Tecnica: Metodo degli Elementi Virtuali per spazi $H(\text{curl})$ e $H(\text{div})$ privi di stabilizzazione

Definizione del Problema
Il Metodo degli Elementi Virtuali (VEM) offre significativi vantaggi rispetto al Metodo degli Elementi Finiti (FEM), consentendo la gestione di politopi stellati generici, il che lo rende ideale per geometrie complesse e raffinamenti adattivi della mesh. Tuttavia, le implementazioni standard del VEM per gli spazi $H(\text{curl})$ e $H(\text{div})$ richiedono termini di stabilizzazione per garantire la coercitività. Questa necessità deriva dal fatto che lo spazio degli elementi virtuali $V(K)$ è più ampio dello spazio polinomiale $P_k(K)$, rendendo instabile la proiezione $L^2$ standard (ovvero $\|\Pi_k v\|_K \not\gtrsim \|v\|_K$). Tali termini di stabilizzazione introducono diversi svantaggi: impattano negativamente sulla stabilità numerica e sui numeri di condizionamento, richiedono una progettazione specifica per il problema (particolarmente per problemi non lineari) ed esibiscono un comportamento isotropo che fatica a gestire coefficienti o mesh anisotropi. Gli approcci esistenti privi di stabilizzazione richiedono o un numero eccessivo di gradi di libertà (DoF) interni, che aumentano i costi computazionali, o modificano lo spazio di proiezione in modi che diluiscono la flessibilità geometrica del VEM o non preservano proprietà essenziali come la commutatività con gli operatori differenziali.

Metodologia
Gli autori propongono un nuovo VEM privo di stabilizzazione per spazi conformi $H(\text{curl})$ e $H(\text{div})$ di ordine generale in $\mathbb{R}^3$. Il nucleo della metodologia consiste nella costruzione di nuovi proiettori serendipity e dei corrispondenti spazi serendipity con un numero minimizzato di DoF.

1. Proiettori Serendipity: A differenza dei lavori precedenti che si affidano a proiettori banali o DoF interni, gli autori definiscono proiettori specifici ($\Pi^S_k$) per gli spazi di facce, spigoli e nodi. Questi proiettori mappano le funzioni virtuali in spazi polinomiali utilizzando solo i DoF di frontiera "preservati" e un insieme specifico di momenti interni. I proiettori sono definiti tramite condizioni di ortogonalità che coinvolgono le tracce normali/tangenziali e i momenti interni (ad esempio, $(v - \Pi^S_k v, x_K \times p_{kr})_K = 0$).
2. Spazi Serendipity: Gli autori definiscono gli spazi serendipity ($V^S$) vincolando gli spazi degli elementi virtuali standard in modo che i DoF "ridotti" siano determinati dai DoF preservati tramite il proiettore serendipity. Ciò garantisce che la proiezione $L^2$ su uno spazio polinomiale di ordine sufficientemente elevato $P_l$ sia stabile e computabile utilizzando solo i DoF di frontiera.
3. Complesso di De Rham: Il metodo è progettato per gestire l'intera catena del complesso di De Rham in $\mathbb{R}^3$ ($H^1 \to H(\text{curl}) \to H(\text{div}) \to L^2$). La costruzione assicura la commutatività esatta con gli operatori differenziali ($\nabla, \nabla\times, \nabla\cdot$), preservando la struttura topologica del problema.
4. Analisi della Stabilità: Gli autori dimostrano che per un ordine di proiezione $l$ sufficientemente grande, l'equivalenza di norma $\|\Pi_l v\|_K \approx \|v\|_K$ è soddisfatta per le funzioni nello spazio serendipity. Ciò elimina la necessità di termini di stabilizzazione mantenendo al contempo proprietà di approssimazione ottimali.

Contributi Chiave

• Trattamento Completo di De Rham: Il lavoro fornisce un trattamento unificato per l'intero complesso di De Rham in 3D, garantendo un'integrazione fluida tra i diversi operatori differenziali, una caratteristica assente nei precedenti lavori privi di stabilizzazione limitati a $H^1$.
• Preservazione delle Proprietà Matematiche: Gli spazi serendipity proposti mantengono rigorosamente la continuità di frontiera, la conformità dell'operatore, la commutatività esatta con gli operatori differenziali e la B-compatibilità.
• Efficienza Computazionale: Minimizzando il numero di DoF e utilizzando proiezioni serendipity calcolate esclusivamente dai dati di frontiera, il metodo riduce significativamente il carico computazionale. Ciò è particolarmente vantaggioso per i problemi di autovalori dove i DoF interni non possono essere eliminati tramite condensazione statica.
• Analisi Rigorosa dell'Errore: Il documento fornisce una prova teorica delle proprietà di approssimazione ottimali e della stabilità per le proiezioni polinomiali di alto ordine.

Risultati
Le affermazioni teoriche sono validate attraverso esperimenti numerici su problemi di autovalori di Maxwell utilizzando mesh poliedriche di Voronoi.

• Tassi di Convergenza: I risultati numerici confermano i tassi di convergenza teorici. Per soluzioni infinitamente regolari, il metodo raggiunge una convergenza $O(h^{m+1})$ nella norma $L^2$ e $O(h^{2m+2})$ per gli autovalori. Per soluzioni singolari (dove la regolarità è limitata da $s$), i tassi sono rispettivamente $O(h^s)$ e $O(h^{2s})$.
• Costo Computazionale: Si dimostra che il tempo di assemblaggio del proposto SFVEM (Stabilization-Free VEM) è solo moderatamente superiore al FEM convenzionale, verificando l'efficienza del metodo. L'ordine di proiezione $l$ richiesto si rivela essere vicino all'ordine polinomiale $k$ nella pratica.
• Confronto: Il metodo elimina con successo le complicazioni legate alla stabilizzazione eguagliando l'accuratezza del FEM convenzionale e del VEM stabilizzato.

Significatività e Rivendicazioni
Gli autori affermano che questo lavoro offre un approccio fondamentalmente diverso al VEM privo di stabilizzazione che supera i limiti delle strategie esistenti (eccessivi DoF o perdita di flessibilità geometrica). Eliminando i termini di stabilizzazione, il metodo migliora la stabilità numerica e i numeri di condizionamento, rendendolo particolarmente prezioso per:

• Formulazioni miste.
• Problemi non lineari.
• Applicazioni che richiedono un comportamento numerico robusto attraverso diverse configurazioni di mesh.

Il documento conclude che, sebbene la tecnica possa essere teoricamente generalizzata a dimensioni superiori, le complessità di implementazione pratica e la rapida crescita dei DoF di frontiera limitano attualmente la discussione all'ambiente 3D. Il metodo proposto rappresenta un'alternativa teoricamente rigorosa e computazionalmente efficiente per la risoluzione di equazioni alle derivate parziali, facendo avanzare l'applicabilità del VEM nella matematica computazionale e nell'ingegneria.