$\boldsymbol{H}(\textbf{curl})$-reconstruction of piecewise polynomial fields with application to $hp$-a posteriori nonconforming error analysis for Maxwell's equations

Aquí presento un resumen técnico detallado del artículo "H(curl)-reconstruction of piecewise polynomial fields with application to hp-a posteriori nonconforming error analysis for Maxwell's equations" (Reconstrucción H(curl) de campos polinomiales a trozos con aplicación al análisis de error a posteriori no conformes hp para las ecuaciones de Maxwell).

1. Planteamiento del Problema

El artículo aborda la necesidad de reconstruir campos vectoriales conformes en el espacio $H(\text{curl})$ a partir de campos polinomiales a trozos definidos sobre mallas simpliciales regulares. Este problema es fundamental en el análisis de error a posteriori para métodos numéricos no conformes (como el Método de Galerkin Discontinuo, dG) aplicados a las ecuaciones de Maxwell.

• Contexto: En la discretización de problemas de Maxwell (operador curl-curl), las soluciones aproximadas suelen pertenecer a espacios de elementos finitos discontinuos o no conformes. Para estimar el error de manera rigurosa, es necesario comparar la solución discreta con una función de referencia que pertenezca al espacio continuo $H_0(\text{curl})$.
• Desafío: Las reconstrucciones existentes en la literatura a menudo requieren que el campo original satisfaga ciertas propiedades de ortogonalidad (como la ortogonalidad de Galerkin) o dependen fuertemente de la topología global del dominio. Además, muchas reconstrucciones polinómicas puras no logran estimaciones óptimas en el grado polinomial $p$ (suboptimalidad en $p$).
• Objetivo: Diseñar un operador de reconstrucción $R_c$ que tome un campo polinomial a trozos arbitrario y produzca un campo en $H_0(\text{curl})$ tal que la diferencia entre ambos esté acotada por los saltos (jumps) del campo original y su rotacional, con tasas de convergencia óptimas en $h$ y subóptimas mínimas en $p$.

2. Metodología

Los autores proponen una construcción local y no polinómica basada en la descomposición de Helmholtz local y particiones de la unidad.

• Estrategia de Construcción:

  1. Partición de la Unidad: Se utiliza la partición de la unidad inducida por las funciones de base "hat" ($\psi_a$) asociadas a los vértices de la malla.
  2. Problemas Locales en Parches de Vértice: Para cada vértice $a$, se define un parche $D_a$ (la unión de las celdas que comparten el vértice). Dado que $D_a$ es un dominio Lipschitz simplemente conexo, se resuelven problemas locales bien planteados.
  3. Descomposición de Helmholtz Local: El campo reconstruido en cada parche, $E_a$, se construye como la suma de dos componentes:
     • $U^a$: Un campo en $H_0(\text{curl}; D_a) \cap H(\text{div}=0; D_a)$ que resuelve un problema variacional relacionado con el rotacional del campo original ponderado por $\psi_a$.
     • $\nabla \theta^a$: Un campo irrotacional (gradiente de un escalar) que resuelve un problema de Poisson para corregir la divergencia.
  4. Operador Global: El campo global reconstruido es la suma de las extensiones por cero de estos campos locales: $R_c(E_h) = \sum_{a \in V_h} E_a$.

• Herramientas Analíticas Clave:

  • Descomposición de Helmholtz Local: Permite separar las partes solenoidales e irrotacionales del error.
  • Desigualdad de Poincaré Rota-Discontinua: Se demuestra una nueva desigualdad de Poincaré en parches de vértice que preserva la divergencia y maneja el rotacional discontinuo.
  • Lifting de Saltos: Se utiliza un operador de levantamiento ($L_h^p$) para manejar los saltos de las condiciones de contorno y de interfaz.
  • Regularidad Elíptica: Para mejorar la estimación en norma $L^2$, se asume una propiedad de "recuperación de regularidad elíptica" uniforme para problemas de Poisson con condiciones de Neumann en los parches.

3. Contribuciones Clave

1. Operador de Reconstrucción General: Se define un operador $R_c: P_p^b(\mathcal{T}_h) \to H_0(\text{curl}; D)$ válido para cualquier campo polinomial a trozos, sin requerir que satisfaga ortogonalidades previas.
2. Independencia Topológica: Al realizar la construcción localmente en parches de vértices (que tienen topología trivial), las constantes de las estimaciones de error no dependen de la topología global del dominio (ej. agujeros, conectividad múltiple), a diferencia de las reconstrucciones globales.
3. Estimaciones hp-Óptimas:
   • Norma de Rotacional Roto (Broken Curl): La estimación es óptima en $h$ y subóptima en $p$ por un orden de $1/2$.
   • Norma $L^2$: La estimación es óptima en $h$ y subóptima en $p$ por un orden de $3/2$.
   • Mejora Condicional: Bajo la suposición de regularidad elíptica uniforme (Assumption 3.5), la suboptimalidad en $p$ para la norma $L^2$ se reduce a $1/2$.
4. Nueva Desigualdad de Poincaré: Se introduce y prueba una desigualdad de Poincaré "broken-curl" que preserva la divergencia en parches de vértices, un resultado de interés independiente.
5. Extensiones: Se discuten extensiones a:
   • Condiciones de contorno mixtas (Dirichlet/Neumann).
   • Mallas polipodales aglomeradas (agglomerated polytopal meshes).
   • Reconstrucciones globales en dominios convexos (que logran optimalidad total en $p$ para la norma $L^2$, pero con constantes dependientes de la topología).

4. Resultados Principales

• Teorema 3.7 (Acotación del Error): Establece que la diferencia entre el campo original $E_h$ y la reconstrucción $R_c(E_h)$ está acotada por los saltos tangenciales de $E_h$ y los saltos normales de su rotacional $\nabla \times E_h$ en las caras de la malla.
  $$\|R_c(E_h) - E_h\| \lesssim C \left( \sum_{K} h_K \|\llbracket E_h \rrbracket^c\|^2 + \dots \right)^{1/2}$$
  Las constantes dependen únicamente de la regularidad de la forma de la malla ($\kappa_{\mathcal{T}_h}$).
• Análisis de Error A Posteriori (Sección 6):
  • Se aplica el operador a la ecuación de Maxwell simplificada discretizada mediante un método Interior Penalty Discontinuous Galerkin (dG).
  • Se demuestra que el error no conforme puede acotarse por un estimador residual local ($\eta$).
  • Teorema 6.8 (Fiabilidad): Se prueba la fiabilidad global del estimador: $\|e\|_{\sharp} \leq C \eta$, donde $e$ es el error total y $\eta$ es la suma de indicadores locales que incluyen residuos de la ecuación, condiciones de contorno y términos de no conformidad.
  • Propiedad 6.9 (Eficiencia): Se demuestra que el estimador es eficiente localmente (acotado inferiormente por el error real más oscilación de datos).

5. Significado e Impacto

Este trabajo representa un avance significativo en el análisis numérico de ecuaciones de Maxwell por varias razones:

1. Primera Análisis hp Completo: Es la primera vez que se logra un análisis a posteriori completo de tipo $hp$ para métodos no conformes en ecuaciones de Maxwell, superando las limitaciones de análisis anteriores que solo eran óptimos en $h$.
2. Eliminación de Restricciones de Ortogonalidad: A diferencia de métodos previos que requerían que el campo discreto fuera una solución de un problema variacional específico (ortogonalidad de Galerkin), esta reconstrucción funciona para campos arbitrarios, lo que la hace más versátil para estimadores de error en esquemas adaptativos.
3. Independencia de la Topología: La capacidad de obtener constantes de error independientes de la topología del dominio es crucial para problemas en geometrías complejas o con múltiples componentes de conexión, donde las reconstrucciones globales fallan o son costosas.
4. Aplicabilidad Práctica: Aunque la reconstrucción se define teóricamente, el artículo aclara que no es necesario calcularla explícitamente en la práctica para obtener el estimador de error; el estimador se basa directamente en los saltos y residuos que la teoría utiliza para acotar el error. Esto hace que el método sea computacionalmente viable para simulaciones de ingeniería.
5. Fundamento para Métodos Avanzados: Proporciona la base teórica necesaria para el desarrollo de estrategias de refinamiento adaptativo ($h$ y $p$) eficientes en simulaciones electromagnéticas de alta precisión.

En resumen, el artículo cierra una brecha teórica importante en el análisis de métodos de elementos finitos discontinuos para Maxwell, ofreciendo herramientas robustas para la estimación de error en configuraciones generales de malla y dominio.