An Extended Galerkin analysis in finite element exterior calculus

Resumen Técnico: Análisis de Galerkin Extendido en el Cálculo Exterior de Elementos Finitos (FEEC)

1. Planteamiento del Problema

El artículo aborda la discretización numérica de la ecuación de Laplace-Hodge en el marco del Cálculo Exterior de Elementos Finitos (FEEC). La ecuación modelo es:
$$(d\delta + \delta d)u = f$$
donde $u$ es una forma diferencial, $d$ es la derivada exterior y $\delta$ es la codiferente.

• Contexto: El método de Galerkin continuo (específicamente el método mixto de Arnold-Falk-Winther o AFW) es el estándar para resolver este problema, pero requiere espacios de elementos finitos conformes (continuos a través de las fronteras de los elementos).
• Desafío: Los métodos de Galerkin Discontinuo (DG) son populares por su flexibilidad y paralelización, pero su aplicación sistemática y análisis unificado en el contexto de FEEC (especialmente para el Laplaciano de Hodge mixto) ha sido limitado. Los métodos DG existentes suelen introducir solo un espacio auxiliar en las fronteras, lo que puede no capturar adecuadamente la complejidad de las condiciones de continuidad impuestas por los operadores $d$ y $\delta$.
• Objetivo: Desarrollar un marco unificado de Galerkin Extendido (XG) que permita formular métodos DG estables y consistentes para la ecuación de Laplace-Hodge, analizando su convergencia y relación con los métodos mixtos conformes clásicos.

2. Metodología: El Marco Galerkin Extendido (XG)

Los autores proponen un enfoque basado en la idea de "truco de Nitsche" extendido, introduciendo un sistema de siete campos para manejar la discontinuidad de manera controlada.

• Variables Principales y Auxiliares:
  El método aproxima tres variables principales ($u$, $du$, $\delta u$) y utiliza cuatro variables de verificación ("check variables") definidas en las fronteras de los elementos:

  1. Campos volumétricos: $\sigma_h \approx -\delta u$, $\xi_h \approx -du$, $u_h \approx u$.
  2. Campos de frontera (Check variables): $\check{\sigma}_h, \check{\xi}^*_h, \check{u}_h, \check{u}^*_h$. Estas variables actúan como multiplicadores de Lagrange débiles para imponer la continuidad inter-elemental.

• Flujos Numéricos:
  Se definen flujos numéricos ($\hat{u}, \hat{\sigma}, \dots$) que combinan promedios centrales y saltos ponderados por las variables de verificación. Por ejemplo:
  $$\hat{u}_h = \{\!\{ \text{tr } u_h \}\!\} + \check{u}_h$$
  donde $\{\!\{ \cdot \}\!\}$ es el promedio y los términos de salto se penalizan mediante parámetros de penalización ($a, b, c, d$).

• Formulación Variacional:
  Se deriva una formulación de siete campos (Ecuación 3.11) que penaliza la discontinuidad tanto respecto a $d$ como a $\delta$. Esta formulación es consistente con el problema continuo y permite eliminar las variables de verificación para obtener formulaciones reducidas (de 4 o 3 campos).

• Análisis de Estabilidad (Inf-Sup):
  El núcleo de la metodología es la demostración de una condición inf-sup (estabilidad de Babuška-Brezzi) uniforme con respecto a todos los parámetros de discretización ($h$) y parámetros de penalización.

  • Se construyen normas dependientes de la malla que incluyen términos de salto y penalización.
  • Se demuestra que el sistema resultante es un sistema de punto de silla perturbado que satisface las condiciones de estabilidad clásicas.

3. Contribuciones Clave

1. Unificación de Métodos DG: El marco XG unifica varias familias de métodos DG existentes bajo una sola estructura teórica. Al ajustar los parámetros de penalización, el método puede recuperar o aproximarse a diferentes esquemas conocidos.
2. Análisis de Siete Campos: La introducción de siete campos (tres internos y cuatro de frontera) permite una transición suave entre la conformidad completa (método AFW) y la discontinuidad total. Esto facilita un análisis de error más transparente.
3. Hibridación (Condensación): Se demuestra que el sistema de siete campos puede hibridizarse. Eliminando las variables volumétricas mediante solucionadores locales, el problema global se reduce a un sistema condensado que solo involucra dos flujos numéricos en las fronteras. Esto reduce drásticamente los grados de libertad globales.
4. Convergencia a Métodos Conformes: Se prueba teóricamente que, cuando los parámetros de penalización tienden a límites adecuados (específicamente cuando la penalización tiende a infinito o cero según el caso), la solución del método XG converge a la solución del método mixto conforme clásico (AFW) o a su dual.
5. Aplicación a Casos Específicos: El marco se aplica explícitamente a:
   • La ecuación de Poisson (casos $k=0$ y $k=3$).
   • El Laplaciano vectorial (casos $k=1$ y $k=2$ en 3D), recuperando y generalizando métodos DG para problemas de elasticidad y Maxwell.

4. Resultados Principales

• Estabilidad Uniforme: Se establece la existencia de una constante de inf-sup $C_{stab} > 0$ que es independiente de los parámetros de penalización y del tamaño de la malla $h$. Esto garantiza que el método es estable para cualquier elección razonable de parámetros.
• Estimaciones de Error Óptimas: Se derivan estimaciones de error a priori cuasi-óptimas:
  $$|||u - u_h||| \leq C h^t |u|_{H^t}$$
  donde la constante $C$ no depende de los parámetros de penalización. Esto asegura que la precisión del método no se degrade al ajustar los parámetros de estabilización.
• Reducción de Dimensionalidad: La versión hibridada del método logra un sistema global con solo dos incógnitas por cara (flujos numéricos), lo que es computacionalmente eficiente comparado con los métodos DG tradicionales que mantienen muchas variables internas.
• Consistencia con FEEC: El método respeta la estructura de la compleja de De Rham discreta, preservando propiedades topológicas y geométricas esenciales del problema continuo.

5. Significado e Impacto

Este trabajo es significativo por varias razones:

• Puente Teórico: Cierra la brecha entre los métodos de elementos finitos conformes (AFW) y los métodos DG en el contexto del cálculo exterior, proporcionando un marco teórico robusto para analizar ambos.
• Flexibilidad Computacional: Al ofrecer un marco unificado, permite a los investigadores diseñar nuevos esquemas DG para problemas complejos (como elasticidad o electromagnetismo) simplemente ajustando los parámetros dentro del marco XG, sin necesidad de desarrollar teorías de estabilidad desde cero para cada variante.
• Eficiencia: La capacidad de hibridar el método para obtener un sistema global pequeño es crucial para la implementación práctica en problemas de gran escala, reduciendo el costo de memoria y tiempo de solución de sistemas lineales.
• Generalidad: Aunque se centra en el Laplaciano de Hodge, el enfoque es aplicable a una amplia gama de problemas gobernados por complejos de Hilbert cerrados, incluyendo las ecuaciones de Maxwell y de elasticidad lineal.

En resumen, Hong, Li y Xu presentan una generalización poderosa y rigurosa de los métodos de Galerkin extendido, proporcionando tanto una herramienta teórica para el análisis de estabilidad como una metodología práctica para la construcción de esquemas numéricos eficientes y flexibles en FEEC.