High-Order Invariant-Domain Preserving Continuous Finite Elements via Graph-Poisson Convex Limiting

Sintesi Tecnica: Elementi Finiti Continui ad Alto Ordine con Preservazione del Dominio Invariante tramite Limiting Convesso su Grafo-Poisson

Enunciato del Probleo
Il documento affronta l'approssimazione numerica di leggi di conservazione scalari non lineari della forma $\partial_t u + \nabla \cdot f(u) = 0$. L'obiettivo primario è la costruzione di metodi a elementi finiti (FEM) continui ad alto ordine che preservino la proprietà del dominio invariante. Tale proprietà garantisce che, se i dati iniziali giacciono all'interno di un dominio specifico $A$, la soluzione di entropia rimanga in $A$ per tutto il tempo. Sebbene gli schemi IDP (Invariant-Domain Preserving) ad alto ordine siano ben stabiliti per i metodi a Volumi Finiti (FV) e Discontinuous Galerkin (DG), costruzioni comparabili per gli elementi Lagrange continui sono meno sviluppate.

Gli approcci IDP ad alto ordine esistenti per i FEM continui affrontano due colli di bottiglia significativi:

1. Problemi di Mass Lumping: Gli schemi IDP a basso ordine richiedono una matrice di massa "lumped" con voci non negative per soddisfare il Principio del Massimo Discreto (DMP). Tuttavia, per gli elementi Lagrange ad alto ordine (ad esempio, $P_k$ con $k \ge 4$ su simplessi), le masse lumping standard possono diventare negative, impedendo l'applicazione diretta della meccanica della viscosità di grafo standard.
2. Ampiezza dello Stencil e Diffusione: L'estensione degli schemi a basso ordine basati sugli archi agli spazi ad alto ordine risulta in stencil di grafo eccessivamente ampi, portando a una diffusione numerica eccessiva.

Metodologia
Il metodo proposto scinde le discretizzazioni ad alto e basso ordine per superare questi ostacoli, collegandole tramite una nuova procedura di ricostruzione del flusso e di limiting.

1. Separazione degli Spazi:

   • Schema ad Alto Ordine: La soluzione ad alto ordine è calcolata nello spazio polinomiale desiderato ($P_k$) utilizzando un metodo di stabilizzazione a Viscosità Residua (RV), potenzialmente aumentato con stabilizzazione Continuous Interior Penalty (CIP) per polinomi di ordine pari.
   • Schema a Basso Ordine: Invece di costruire uno schema a basso ordine sullo spazio ad alto ordine, gli autori costruiscono uno schema IDP del primo ordine ottimale su una sottorete $P_1$ fine. I vertici di questa sottorete coincidono esattamente con i nodi dello spazio ad alto ordine. Ciò assicura che lo schema a basso ordine operi su un grafo $P_1$ con masse lumping strettamente positive, evitando il problema delle masse negative.

2. Ricostruzione del Flusso tramite Grafo-Poisson:
   Poiché gli aggiornamenti ad alto e basso ordine sono definiti rispetto a diverse strutture di massa (la massa consistente ad alto ordine rispetto alla massa lumping fine $P_1$), la loro differenza non può essere rappresentata direttamente come flussi conservativi sugli archi.

   • Gli autori definiscono un vettore di difetto di massa che rappresenta la differenza tra gli aggiornamenti ad alto e basso ordine.
   • Questo difetto viene ricostruito come flussi di grafo conservativi e antisimmetrici risolvendo un problema di Grafo-Poisson (un problema di Neumann discreto). Questo passaggio mappa efficacemente la differenza tra i due schemi sul grafo $P_1$ fine.

3. Limiting Convesso e Ridistribuzione della Massa:
   Il processo di limiting segue una strategia in due fasi ispirata all'Algebraic Flux Correction (AFC), ma adattata per elementi continui:

   • Fase 1 (Limiting Convesso): I flussi di grafo ricostruiti vengono limitati utilizzando coefficienti di limiting convesso. Questo passaggio impone vincoli locali (il dominio invariante) sulla soluzione, ma può introdurre un difetto di massa rispetto al target ad alto ordine poiché il limiting è eseguito sulla struttura di massa $P_1$ fine.
   • Fase 2 (Ridistribuzione della Massa basata sulla Capacità): Una fase secondaria ridistribuisce il difetto di massa rimanente. Essa utilizza la "capacità ammissibile" di ogni nodo (la distanza dai vincoli locali) per ripristinare la massa del target ad alto ordine, preservando rigorosamente i vincoli del dominio invariante. Questo passaggio si basa sull'ipotesi che la massa del target ad alto ordine ricada nell'intervallo di massa ammissibile definito dai vincoli.

Contributi Chiave

• Decoupling degli Spazi: Il documento introduce un framework in cui lo schema IDP a basso ordine è costruito su una sottorete $P_1$ fine indotta dai nodi ad alto ordine. Ciò evita il problema delle masse lumping negative inerenti alle basi Lagrange ad alto ordine e evita la diffusione eccessiva di stencil ampi.
• Ricostruzione Grafo-Poisson: Viene proposto un nuovo metodo per rappresentare la differenza tra gli aggiornamenti ad alto e basso ordine come flussi di grafo conservativi risolvendo un problema di Laplaciano di grafo. Questo sostituisce la sottrazione algebrica diretta usata nei metodi FCT/AFC tradizionali, accomodando il disallineamento tra le matrici di massa.
• Limiting a Due Fasi: La combinazione di limiting convesso sui flussi di grafo seguita da una fase di ridistribuzione della massa basata sulla capacità garantisce sia la preservazione dei domini invarianti che la conservazione globale rispetto alla massa $P_1$ fine.

Risultati Numerici
Il metodo è stato testato su diversi benchmark utilizzando la libreria FEniCSx con mesh non strutturate e gradi polinomiali fino a $P_4$:

• Trasporto Smooth: Il metodo ha dimostrato tassi di convergenza ottimali nelle norme $L_1$, $L_2$ e $L_\infty$ per problemi regolari (smooth), confermando che il limiter rimane inattivo nei regimi regolari.
• Rotazione a Tre Corpi: Il metodo ha preservato con successo la forma delle discontinuità e degli estremi regolari senza oscillazioni spurie. L'uso di vincoli rilassati (dipendenti dalla mesh) è risultato necessario per evitare il "clipping" degli estremi locali e mantenere l'accuratezza ad alto ordine.
• Equazione di Burgers: Il metodo è convergente verso la soluzione di entropia corretta per un problema contenente un punto sonico. I confronti hanno mostrato che l'uso della stabilizzazione a Viscosità Residua (RV) era cruciale; i metodi che utilizzavano solo la stabilizzazione CIP fallivano nel convergere alla soluzione di entropia per gradi polinomiali più elevati.
• Problema KPP: In un benchmark impegnativo con onde d'urto composte, il metodo RV limitato è convergente alla soluzione di entropia, mentre il metodo che utilizzava solo la stabilizzazione CIP convergevano a una soluzione debole non entropica.

Significatività e Claim
Il documento sostiene di fornire un metodo a elementi finiti continui ad alto ordine, robusto e con preservazione del dominio invariante, che supera i limiti specifici dell'applicazione delle tecniche IDP standard agli elementi Lagrange ad alto ordine. Utilizzando la separazione degli spazi e la ricostruzione Grafo-Poisson, il metodo raggiunge:

1. Robustezza: Gestisce discontinuità e punti sonici senza violare i vincoli fisici.
2. Accuratezza: Mantiene l'accuratezza ad alto ordine nelle regioni regolari e converge alla soluzione di entropia per problemi non lineari.
3. Generalità: L'approccio è applicabile ad arbitrari gradi polinomiali e mesh non strutturate, aggirando le restrizioni imposte dalle masse lumping negative nelle basi Lagrange ad alto ordine standard.

Gli autori sottolineano che la ricostruzione del flusso Grafo-Poisson è una caratteristica distintiva, che permette al metodo di funzionare come una correzione di tipo AFC dove i flussi antidiffusivi vengono ricostruiti tramite un problema basato su PDE (grafo) piuttosto che tramite una diretta sottrazione algebrica.