Multigrid Preconditioning for FEEC using Mass-Lumping and Transforming Smoothers

Resumen Técnico: Precondicionamiento de Multigrid para FEEC mediante Mass-Lumping y Smoothers de Transformación

Planteamiento del Problema
El artículo aborda el desafío de resolver sistemas lineales que surgen de las discretizaciones de Cálculo Exterior de Elementos Finitos (FEEC) de ecuaciones diferenciales parciales (PDE) planteadas en el complejo de de Rham. Específicamente, se dirige a formulaciones mixtas y de punto de silla para el operador de Hodge-Dirac, el Hodge-Laplaciano y sistemas de magnetostática (ecuaciones de Maxwell en la lazo de Coulomb). Estas discretizaciones suelen generar sistemas lineales indefinidos con estructuras de punto de silla. Los suavizadores iterativos estándar, como Gauss-Seidel o Jacobi, son ineficaces para estos sistemas indefinidos debido a que los operadores carecen de la definición positiva requerida para tales esquemas de relajación. Si bien trabajos previos han propuesto enfoques de multigrid basados en complementos de Schur, este artículo busca una estrategia alternativa que evite la complejidad de la inversión del complemento de Schur.

Metodología
El marco propuesto combina el mass-lumping (agrupación de masa) con smoothers de transformación para construir un precondicionador efectivo para los sistemas consistentes de FEEC.

1. Mass-Lumping y Equivalencia Espectral:
   Los autores introducen una versión de los operadores FEEC con mass-lumping donde las matrices de masa consistentes son reemplazadas por matrices diagonales (o explícitamente invertibles). Esto permite la inversión explícita del isomorfismo de Riesz, un paso que de otro modo sería prohibitivamente costoso con matrices de masa consistentes densas.

   • Bajo supuestos de equivalencia de norma $h$-uniformes leves entre los productos internos consistentes y los de mass-lumped, el artículo demuestra una desigualdad de Poincaré discreta para las normas de mass-lumping.
   • Este resultado de estabilidad conduce a una prueba de equivalencia espectral entre los operadores FEEC consistentes originales y los operadores de mass-lumping. Esta equivalencia justifica el uso del sistema de mass-lumping como un sustituto para precondicionar el sistema original.

2. Smoothers de Transformación:
   Para manejar la indefinición de los operadores de mass-lumping, el método emplea smoothers de transformación. La idea central es transformar algebraicamente el operador en una forma de bloque donde los bloques diagonales sean definidos positivos, permitiendo el uso de smoothers estándar como Gauss-Seidel.

   • Para los problemas de Hodge-Dirac, Hodge-Laplaciano y magnetostática, se construyen matrices de transformación específicas de izquierda y derecha ($T_L, T_R$).
   • Estas transformaciones mapean el sistema indefinido original a un sistema triangular inferior por bloques (o un sistema con bloques diagonales definidos positivos).
   • Se aplica una relajación de Gauss-Seidel a este sistema transformado y el resultado se mapea de vuelta al espacio original.

3. Precondicionamiento de Multigrid:
   El algoritmo de multigrid está diseñado específicamente para el sistema de mass-lumping transformado. Este ciclo de multigrid ($M_{MG}$) sirve como una aproximación inversa para el operador de mass-lumping.

   • Para resolver el sistema FEEC original $G_h u = g$, el precondicionador $B$ se construye como $B \approx J_h^{-1} M_{MG} \tilde{H}_h (J_h^*)^{-1}$, donde $J_h$ es el isomorfismo entre los espacios estándar y de mass-lumping (a menudo la identidad en la práctica) y $\tilde{H}_h$ es el mapa de Riesz de mass-lumping.
   • Debido a que el sistema precondicionado resultante sigue siendo indefinido, los autores utilizan GMRES (método de Residuos Mínimos Generalizados) como el solver externo, ya que los métodos estándar como PCG o MINRES requieren precondicionadores definidos positivos.

Contribuciones Clave

• Diseño Algorítmico: El artículo propone un marco de precondicionamiento de multigrid unificado que integra mass-lumping con smoothers de transformación para problemas de FEEC.
• Estabilidad Teórica: Establece la estabilidad de los sistemas de mass-lumping y demuestra la equivalencia espectral con los operadores FEOC originales bajo supuestos leves, proporcionando una base teórica para la estrategia de precondicionamiento.
• Manejo de la Indefinición: El trabajo demuestra cómo los smoothers de transformación pueden convertir operadores de bloque indefinidos en formas susceptibles a las técnicas de relajación clásicas sin requerir aproximaciones del complemento de Schur.
• Robustez en Topologías Complejas: Se muestra que el método maneja eficazmente dominios con topologías no triviales (formas armónicas), donde la deficiencia de rango del sistema está determinada por los números de Betti.

Resultados Numéricos
Se realizaron experimentos numéricos extensos en 2D y 3D en diversas mallas, incluyendo aquellas con esquinas entrantes y agujeros (topología no trivial).

• Convergencia: Los ciclos de multigrid (ciclos V y W) combinados con los smoothers de transformación demostraron una convergencia robusta. Los radios espectrales de los operadores de propagación de error se mantuvieron alejados de la unidad, independientemente del tamaño de malla $h$.
• Estrategias de Mass-Lumping: Se probaron tres estrategias de mass-lumping: identidad escalada, sumas de filas y duales baricéntricos (inspirados en el Cálculo Exterior Discreto). Los resultados indicaron que estas diferentes estrategias de mass-lumping no condujeron a diferencias cualitativas pronunciadas en el comportamiento de convergencia.
• Rendimiento del Solver: Cuando se utilizó como precondicionador para GMRES, el método produjo conteos de iteración que fueron en gran medida independientes del nivel de refinamiento de la malla. El ciclo W generalmente proporcionó una convergencia más robusta que el ciclo V.
• Formas Armónicas: Los experimentos confirmaron que el método maneja correctamente la deficiencia de rango asociada con las formas armónicas, donde el número de modos propios no convergentes coincide con la dimensión del espacio armónico.

Significancia y Afirmaciones
Los autores posicionan este trabajo principalmente como una contribución algorítmica más que como una prueba de teoría de convergencia formal. Declaran explícitamente que, si bien demuestran la equivalencia espectral y la estabilidad, establecer una teoría de convergencia rigurosa que garantice límites uniformes en los conteos de iteración de GMRES sigue siendo un problema abierto para investigaciones futuras.

La importancia del enfoque radica en su capacidad para proporcionar un precondicionador práctico y robusto para sistemas FEEC indefinidos utilizando operaciones explícitas y económicas (mass-lumping) y transformaciones algebraicas, evitando las complejidades de los métodos de complemento de Schur. La evidencia numérica sugiere que este enfoque es una dirección válida y prometedora para resolver problemas de elementos finitos mixtos en electromagnetismo y campos relacionados, lo que justifica un mayor estudio teórico y su extensión a otros problemas como el flujo de Stokes o sistemas con saltos de coeficientes.