Multigrid Preconditioning for FEEC using Mass-Lumping and Transforming Smoothers

Technische Samenvatting: Multigrid-preconditionering voor FEEC met behulp van Mass-Lumping en Transformerende Smoothers

Probleemstelling
Het artikel behandelt de uitdaging van het oplossen van lineaire systemen die voortvloeien uit Finite Element Exterior Calculus (FEEC)-discretisaties van partiële differentiaalvergelijkingen (PDE's) gepositioneerd op de de Rham-complex. Specifiek richt het zich op gemengde en saddle-point formuleringen voor de Hodge-Dirac-operator, de Hodge-Laplace-operator en magnetostatische systemen (Maxwell-vergelijkingen in de Coulomb-gauge). Deze discretisaties resulteren doorgaans in ongedefinieerde lineaire systemen met saddle-point structuren. Standaard iteratieve smoothers, zoals Gauss-Seidel of Jacobi, zijn ineffectief voor deze ongedefinieerde systemen omdat de operatoren de positieve definitiviteit missen die vereist is voor dergelijke relaxatieschema's. Terwijl eerder werk voorstelde om multigrid-benaderingen te baseren op Schur-complementen, streeft dit artikel naar een alternatieve strategie die de complexiteit van Schur-complement inversie vermijdt.

Methodologie
Het voorgestelde kader combineert mass-lumping met transformerende smoothers om een effectieve preconditioner te construeren voor de consistente FEEC-systemen.

1. Mass-Lumping en Spectrale Equivalentie:
   De auteurs introduceren een mass-lumped versie van de FEEC-operatoren waarbij de consistente massamatrices worden vervangen door diagonale (of expliciet inverseerbare) matrices. Dit maakt de expliciete inversie van de Riesz-isomorfisme mogelijk, een stap die anders onbetaalbaar duur zou zijn met dichte consistente massamatrices.

   • Onder milde aannames van $h$-uniforme norm-equivalentie tussen de consistente en de mass-lumped inner products, bewijzen de auteurs een discrete Poincaré-ongelijkheid voor de mass-lumped normen.
   • Dit stabiliteitsresultaat leidt tot een bewijs van spectrale equivalentie tussen de oorspronkelijke consistente FEEC-operatoren en de mass-lumped operatoren. Deze equivalentie rechtvaardigt het gebruik van de mass-lumped systemen als proxy voor het preconditioneren van het oorspronkelijke systeem.

2. Transformerende Smoothers:
   Om de ongedefinieerdheid van de mass-lumped operatoren aan te pakken, gebruikt de methode transformerende smoothers. De kern van het idee is om de operator algebraïsch te transformeren naar een blokvorm waarbij de diagonale blokken positief definitief zijn, wat het gebruik van standaard smoothers zoals Gauss-Seidel mogelijk maakt.

   • Voor de Hodge-Dirac, Hodge-Laplace en magnetostatische problemen worden specifieke linker- en rechter-transformatiematrices ($T_L, T_R$) geconstrueerd.
   • Deze transformaties mappen het oorspronkelijke ongedefinieerde systeem naar een lager-blok-triangulair systeem (of een systeem met positief definitieve diagonale blokken).
   • Een Gauss-Seidel-relaxatie wordt toegepast op dit getransformeerde systeem, en het resultaat wordt teruggekaart naar de oorspronkelijke ruimte.

3. Multigrid-Preconditionering:
   Het multigrid-algoritme is specifiek ontworpen voor het mass-lumped, getransformeerde systeem. Deze multigrid-cyclus ($M_{MG}$) dient als een benaderde inverse voor de mass-lumped operator.

   • Om het oorspronkelijke FEEC-systeem $G_h u = g$ op te lossen, wordt de preconditioner $B$ geconstrueerd als $B \approx J_h^{-1} M_{MG} \tilde{H}_h (J_h^*)^{-1}$, waarbij $J_h$ de isomorfisme is tussen de standaard en de mass-lumped ruimtes (vaak de identiteit in de praktijk) en $\tilde{H}_h$ de mass-lumped Riesz-map is.
   • Omdat het resulterende gepreconditioneerde systeem ongedefinieerd blijft, gebruiken de auteurs GMRES (Generalized Minimal Residual method) als de externe solver, aangezien standaardmethoden zoals PCG of MINRES positief definitieve preconditioners vereisen.

Belangrijkste Bijdragen

• Algoritmisch Ontwerp: Het artikel stelt een verenigd multigrid-preconditioneringskader voor dat mass-lumping combineert met transformerende smoothers voor FEEC-problemen.
• Theoretische Stabiliteit: Het vestigt de stabiliteit van de mass-lumped systemen en bewijst de spectrale equivalentie met de oorspronkelijke FEEC-operatoren onder milde aannames, wat een theoretische basis biedt voor de preconditioning-strategie.
• Aanpak van Ongedefinieerdheid: Het werk laat zien hoe transformerende smoothers ongedefinieerde blok-operatoren kunnen omzetten in vormen die geschikt zijn voor klassieke relaxatietechnieken zonder dat hiervoor Schur-complement benaderingen nodig zijn.
• Robuustheid bij Complexe Topologieën: De methode wordt getoond effectief te zijn in het afhandelen van domeinen met niet-triviale topologieën (harmonische vormen), waarbij de rangtekortkoming van het systeem wordt bepaald door de Betti-getallen.

Numerieke Resultaten
Uitgebreide numerieke experimenten werden uitgevoerd in 2D en 3D op diverse meshes, inclusclusief meshes met reentrant hoeken en gaten (niet-triviale topologie).

• Convergentie: De multigrid-cycli (V-cyclus en W-cyclus) in combinatie met de transformerende smoothers vertoonden een robuuste convergentie. De spectrale radii van de foutvoortplantingsoperatoren bleven begrensd weg van één, onafhankelijk van de mesh-grootte $h$.
• Mass-Lumping Strategieën: Drie mass-lumping strategieën werden getest: geschaalde identiteit, rijsommen en barycentrische dualen (geïnspireerd door Discrete Exterior Calculus). De resultaten gaven aan dat deze verschillende mass-lumping strategieën niet leidden tot uitgesproken kwalitatieve verschillen in convergentiegedrag.
• Solver Prestaties: Wanneer gebruikt als preconditioner voor GMRES, leverde de methode iteratie-aantallen op die grotendeels onafhankelijk waren van het niveau van mesh-verfijning. De W-cyclus bood over het algemeen een robuustere convergentie dan de V-cyclus.
• Harmonische Vormen: De experimenten bevestigden dat de methode de rangtekortkoming die geassocieerd wordt met harmonische vormen correct afhandelt, waarbij het aantal niet-convergerende eigenmodi overeenkomt met de dimensie van de harmonische ruimte.

Betekenis en Claims
De auteurs positioneren dit werk primair als een algoritmische bijdrage in plaats van een formeel bewijs van convergentietheorie. Ze stellen expliciet dat, hoewel ze spectrale equivalentie en stabiliteit bewijzen, een rigoureuze convergentietheorie die uniforme bounds op GMRES-iteratie-aantallen vaststelt, een openstaand probleem blijft voor toekomstig onderzoek.

De betekenis van de aanpak ligt in het vermogen om een praktische, robuuste preconditioner te bieden voor ongedefinieerde FEEC-systemen met behulp van expliciete, goedkope operaties (mass-lumping) en algebraïsche transformaties, waarbij de complexiteit van Schur-complement methoden wordt vermeden. Het numerieke bewijs suggereert dat de aanpak een geldige en veelbelovende richting is voor het oplossen van gemengde eindige-elementenproblemen in elektromagnetisme en gerelateerde velden, wat verdere theoretische studie en uitbreiding naar andere problemen zoals Stokes-stroming of systemen met coëfficiënt-sprongen rechtvaardigt.