Multigrid Preconditioning for FEEC using Mass-Lumping and Transforming Smoothers

Résumé Technique : Préconditionnement Multigrid pour la FEEC par Mass-Lumping et Smoothers de Transformation

Énoncé du Problème
L'article traite du défi que représente la résolution de systèmes linéaires issus de discrétisations par Calcul de l'Extérieur par Éléments Finis (FEEC) d'équations aux dérivées partielles (EDP) posées sur le complexe de de Rham. Plus précisément, il cible les formulations mixtes et de points de selle pour l'opérateur de Hodge-Dirac, le Hodge-Laplacien, et les systèmes magnétostatiques (équations de Maxwell sous la jauge de Coulomb). Ces discrétisations produisent typiquement des systèmes linéaires indéfinis avec des structures de points de selle. Les lisseurs itératifs standards, tels que Gauss-Seidel ou Jacobi, sont inefficaces pour ces systèmes indéfinis car les opérateurs ne possèdent pas la positivité définie requise pour ces schémas de relaxation. Bien que des travaux antérieurs aient proposé des approches multigrid basées sur les compléments de Schur, cet article cherche une stratégie alternative qui évite la complexité de l'inversion du complément de Schur.

Méthodologie
Le cadre proposé combine le mass-lumping (aplatissement de masse) avec des smoothers de transformation pour construire un préconditionneur efficace pour les systèmes FEEC cohérents.

1. Mass-Lumping et Équivalence Spectrale :
   Les auteurs introduisent une version par mass-lumping des opérateurs FEEC où les matrices de masse cohérentes sont remplacées par des matrices diagonales (ou explicitement inversibles). Cela permet l'inversion explicite de l'isomorphisme de Riesz, une étape autrement prohibitive avec des matrices de masse cohérentes denses.

   • Sous des hypothèses de l'ordre de la norme $h$-uniforme entre les produits intérieurs cohérents et par mass-lumping, l'article prouve une inégalité de Poincaré discrète pour les normes par mass-lumping.
   • Ce résultat de stabilité conduit à une preuve d'équivalence spectrale entre les opérateurs FEEC cohérents originaux et les opérateurs par mass-lumping. Cette équivalence justifie l'utilisation du système par mass-lumping comme substitut pour préconditionner le système original.

2. Smoothers de Transformation :
   Pour gérer l'indéfinitude des opérateurs par mass-lumping, la méthode emploie des smoothers de transformation. L'idée centrale est de transformer algébriquement l'opérateur en une forme bloc où les blocs diagonaux sont définis positifs, permettant l'utilisation de smoothers standards comme Gauss-Seidel.

   • Pour les problèmes de Hodge-Dirac, Hodge-Laplacien et de magnétostatique, des matrices de transformation gauche et droite spécifiques ($T_L, T_R$) sont construites.
   • Ces transformations projettent le système indéfini original vers un système triangulaire inférieur par blocs (ou un système avec des blocs diagonaux définis positifs).
   • Une relaxation de Gauss-Seidel est appliquée à ce système transformé, et le résultat est projeté en retour dans l'espace original.

3. Préconditionnement Multigrid :
   L'algorithme multigrid est conçu spécifiquement pour le système par mass-lumping et transformé. Ce cycle multigrid ($M_{MG}$) sert d'inverse approximatif pour l'opérateur par mass-lumping.

   • Pour résoudre le système FEEC original $G_h u = g$, le préconditionneur $B$ est construit comme $B \approx J_h^{-1} M_{MG} \tilde{H}_h (J_h^*)^{-1}$, où $J_h$ est l'isomorphisme entre les espaces standard et par mass-lumping (souvent l'identité en pratique) et $\tilde{H}_h$ est la carte de Riesz par mass-lumping.
   • Parce que le système préconditionné résultant reste indéfini, les auteurs utilisent GMRES (méthode de résidus minimaux généralisée) comme solveur externe, car les méthodes standards comme PCG ou MINRES nécessitent des préconditionneurs définis positifs.

Contributions Clés

• Conception Algorithmique : L'article propose un cadre de préconditionnement multigrid unifié qui intègre le mass-lumping avec des smoothers de transformation pour les problèmes FEEC.
• Stabilité Théorique : Il établit la stabilité des systèmes par mass-lumping et prouve l'équivalence spectrale avec les opérateurs FEEC originaux sous des hypothèses modérées, fournissant une base théorique à la stratégie de préconditionnement.
• Gestion de l'Indéfinitude : Le travail démontre comment les smoothers de transformation peuvent convertir des opérateurs par blocs indéfinis en formes adaptées aux techniques de relaxation classiques sans nécessiter d'approximations du complément de Schur.
• Robustesse sur les Topologies Complexes : La méthode est montrée comme étant capable de gérer efficacement les domaines aux topologies non triviales (formes harmoniques), où la déficience de rang du système est déterminée par les nombres de Betti.

Résultats Numériques
Des expériences numériques approfondies ont été menées en 2D et 3D sur diverses mailles, incluant des mailles avec des coins rentrants et des trous (topologie non triviale).

• Convergence : Les cycles multigrid (V-cycle et W-cycle) combinés aux smoothers de transformation ont démontré une convergence robuste. Les rayons spectraux des opérateurs de propagation d'erreur sont restés bornés loin de l'unité, indépendamment de la taille de la maille $h$.
• Stratégies de Mass-Lumping : Trois stratégies de mass-lumping ont été testées : identité scalée, sommes de lignes, et duals barycentriques (inspirés du Calcul Extérieur Discret). Les résultats ont indiqué que ces différentes stratégies n'entraînaient pas de différences qualitatives prononcées dans le comportement de convergence.
• Performance du Solveur : Utilisée comme préconditionneur pour GMRES, la méthode a produit des nombres d'itérations largement indépendants du niveau de raffinement de la maille. Le W-cycle a généralement fourni une convergence plus robuste que le V-cycle.
• Formes Harmoniques : Les expériences ont confirmé que la méthode gère correctement la déficience de rang associée aux formes harmoniques, le nombre de modes propres non convergents correspondant à la dimension de l'espace harmonique.

Signification et Revendications
Les auteurs positionnent ce travail principalement comme une contribution algorithmique plutôt qu'une preuve de théorie de convergence formelle. Ils déclarent explicitement que, bien qu'ils prouvent l'équivalence spectrale et la stabilité, une théorie de convergence rigoureuse établissant des bornes uniformes sur les comptes d'itérations GMRES reste un problème ouvert pour des investigations futures.

La signification de l'approche réside dans sa capacité à fournir un préconditionneur pratique et robuste pour les systèmes FEEC indéfinis en utilisant des opérations explicites et peu coûteuses (mass-lumping) et des transformations algébriques, évitant ainsi les complexités des méthodes de complément de Schur. Les preuves numériques suggèrent que cette approche est une direction valide et prometteuse pour résoudre les problèmes de éléments finis mixtes en électromagnétisme et domaines connexes, justifiant une étude théorique et une extension supplémentaires à d'autres problèmes comme l'écoulement de Stokes ou les systèmes avec des sauts de coefficients.