Multigrid Preconditioning for FEEC using Mass-Lumping and Transforming Smoothers

Sintesi Tecnica: Precondizionamento Multigrid per FEEC tramite Mass-Lumping e Smoothers Trasformanti

Problema
Il documento affronta la sfida di risolvere i sistemi lineari derivanti dalle discretizzazioni Finite Element Exterior Calculus (FEEC) di equazioni differenziali alle derivate parziali (PDE) poste sul complesso di de Rham. Nello specifico, si concentra su formulazioni miste e di punto di sella per l'operatore di Hodge-Dirac, l'Hodge-Laplaciano e sistemi magnetostatici (equazioni di Maxwell nel gauge di Coulomb). Queste discretizzazioni generano tipicamente sistemi lineari indefiniti con strutture a punto di sella. Gli standard smussatori iterativi, come Gauss-Seidel o Jacobi, sono inefficaci per questi sistemi indefiniti poiché gli operatori mancano della definitezza positiva richiesta da tali schemi di rilassamento. Sebbene lavori precedenti abbiano proposto approcci multigrid basati sui complementi di Schur, questo articolo cerca una strategia alternativa che eviti la complessità dell'inversione del complemento di Schur.

Metodologia
Il framework proposto combina il mass-lumping con gli smoothers trasformanti per costruire un precondizionatore efficace per i sistemi FEEC consistenti.

1. Mass-Lumping e Equivalenza Spettrale:
   Gli autori introducono una versione del mass-lumping degli operatori FEEC in cui le matrici di massa consistenti sono sostituite da matrici diagonali (o esplicitamente invertibili). Ciò consente l'inversione esplicita dell'isomorfismo di Riesz, un passaggio che altrimenti sarebbe proibitivamente costoso con matrici di massa consistenti dense.

   • Sotto deboli assunzioni di equivalenza di norma $h$-uniforme tra le norme consistenti e quelle mass-lumped, il documento prova una disuguaglianza di Poincaré discreta per le norme mass-lumped.
   • Questo risultato di stabilità conduce a una prova di equivalenza spettrale tra gli operatori FECC originali e quelli mass-lumped. Questa equivalenza giustifica l'uso del sistema mass-lumped come proxy per il precondizionamento del sistema originale.

2. Smoothers Trasformanti:
   Per gestire l'indefinitezza degli operatori mass-lumped, il metodo impiega smoother trasformanti. L'idea centrale è trasformare algebricamente l'operatore in una forma a blocchi in cui i blocchi diagonali sono definiti positivi, consentendo l'uso di smoother standard come Gauss-Seidel.

   • Per i problemi di Hodge-Dirac, Hodge-Laplaciano e magnetostatica, vengono costruite specifiche matrici di trasformazione sinistra e destra ($T_L, T_R$).
   • Queste trasformazioni mappano il sistema indefinito originale in un sistema a blocchi triangolare inferiore (o un sistema con blocchi diagonali definiti positivi).
   • Un rilassamento di Gauss-Seidel viene applicato a questo sistema trasformato, e il risultato viene riportato nello spazio originale.

3. Precondizionamento Multigrid:
   L'algoritmo multigrid è progettato specificamente per il sistema mass-lumped e trasformato. Questo ciclo multigrid ($M_{MG}$) funge da approssimazione dell'inversa dell'operatore mass-lumped.

   • Per risolvere il sistema FEEC originale $G_h u = g$, il precondizionatore $B$ è costruito come $B \approx J_h^{-1} M_{MG} \tilde{H}_h (J_h^*)^{-1}$, dove $J_h$ è l'isomorfismo tra gli spazi standard e mass-lumped (spesso l'identità nella pratica) e $\tilde{H}_h$ è la mappa di Riesz mass-lumped.
   • Poiché il sistema precondizionato risultante rimane indefinito, gli autori utilizzano GMRES (metodo dei residui minimi generalizzato) come solutore esterno, dato che i metodi standard come PCG o MINRES richiedono precondizionatori definiti positivi.

Contributi Chiave

• Design Algoritmico: Il documento propone un framework di precondizionamento multigrid unificato che integra il mass-lumping con gli smoother trasformanti per problemi FEEC.
• Stabilità Teorica: Stabilisce la stabilità dei sistemi mass-lumped e prova l'equivalenza spettrale con gli operatori FECC originali sotto deboli assunzioni, fornendo una base teorica per la strategia di precondizionamento.
• Gestione dell'Indefinitezza: Il lavoro dimostra come gli smoother trasformanti possano convertire operatori a blocchi indefiniti in forme accessibili alle tecniche di rilassamento classiche senza richiedere approssimazioni del complemento di Schur.
• Robustezza su Topologie Complesse: Il metodo si dimostra efficace nel gestire domini con topologie non triviali (forme armoniche), dove il deficit di rango del sistema è determinato dai numeri di Betti.

Risultati Numerici
Sono stati condotti estesi esperimenti numerici in 2D e 3D su varie mesh, incluse quelle con angoli rientranti e fori (topologia non banale).

• Convergenza: I cicli multigrid (V-cycle e W-cycle) combinati con gli smoother trasformanti hanno dimostrato una convergenza robusta. I raggi spettrali degli operatori di propagazione dell'errore sono rimasti limitati lontano dall'unità, indipendentemente dalla dimensione della mesh $h$.
• Strategie di Mass-Lumping: Sono state testate tre strategie di mass-lumping: identità scalata, somme di righe e duali baricentrici (ispirate alla Discrete Exterior Calculus). I risultati hanno indicato che queste diverse strategie di mass-lumping non hanno portato a differenze qualitative pronunciate nel comportamento della convergenza.
• Performance del Solutore: Quando utilizzato come precondizionatore per GMRES, il metodo ha prodotto conteggi di iterazioni ampiamente indipendenti dal livello di raffinamento della mesh. Il W-cycle ha generalmente fornito una convergenza più robusta rispetto al V-cycle.
• Forme Armoniche: Gli esperimenti hanno confermato che il metodo gestisce correttamente il deficit di rango associato alle forme armoniche, con il numero di autovettori non convergenti che corrisponde alla dimensione dello spazio armonico.

Significatività e Rivendicazioni
Gli autori posizionano questo lavoro principalmente come un contributo algoritmico piuttosto che come una prova formale di teoria della convergenza. Affermano esplicitamente che, sebbene abbiano provato l'equivalenza spettrale e la stabilità, una rigorosa teoria della convergenza che stabilisca limiti uniformi sui conteggi delle iterazioni GMRES rimane un problema aperto per indagini future.

La significatività dell'approccio risiede nella sua capacità di fornire un precondizionatore pratico e robusto per sistemi FEEC indefiniti utilizzando operazioni esplicite e poco costose (mass-lumping) e trasformazioni algebriche, evitando le complessità dei metodi del complemento di Schur. L'evidenza numerica suggerisce che l'approccio è una direzione valida e promettente per la risoluzione di problemi di elementi finiti misti nell'elettromagnetismo e campi correlati, giustificando ulteriori studi teorici ed estensioni ad altri problemi come il flusso di Stokes o sistemi con salti di coefficienti.