Stabilization-Free H(curl) and H(div)-Conforming Virtual Element Method

Resumen Técnico: VEM de Conformidad $H(\text{curl})$ y $H(\text{div})$ Libre de Estabilización

Planteamiento del Problema
El Método de Elementos Virtuales (VEM) ofrece ventajas significativas sobre el Método de Elementos Finitos (FEM) al manejar politopos estrellados generales, lo que lo hace ideal para geometrías complejas y el refinamiento de malla adaptativo. Sin embargo, las implementaciones estándar de VEM para espacios $H(\text{curl})$ y $H(\text{div})$ requieren términos de estabilización para asegurar la coercitividad. Esta necesidad surge porque el espacio de elementos virtuales $V(K)$ es más grande que el espacio polinomial $P_k(K)$, lo que causa que la proyección $L^2$ estándar sea inestable (es decir, $\|\Pi_k v\|_K \not\gtrsim \|v\|_K$). Estos términos de estabilización introducen varios inconvenientes: impactan negativamente la estabilidad numérica y los números de condición, requieren un diseño específico para cada problema (particularmente para problemas no lineales) y exhiben un comportamiento isotrópico que tiene dificultades con coeficientes o mallas anisotrópicas. Los enfoques existentes libres de estabilización requieren grados de libertad (DoF) internos excesivos que aumentan los costos computacionales o modifican el espacio de proyección de formas que diluyen la flexibilidad geométrica del VEM o fallan en preservar propiedades esenciales como la conmutatividad con los operadores diferenciales.

Metodología
Los autores proponen un nuevo VEM libre de estabilización para espacios de conformidad $H(\text{curl})$ y $H(\text{div})$ de orden general en $\mathbb{R}^3$. El núcleo de la metodología consiste en la construcción de nuevos proyectores de serendipity y los correspondientes espacios de serendipity con un número minimizado de DoF.

1. Proyectores de Serendipity: A diferencia de trabajos previos que dependen de proyectores triviales o DoF internos, los autores definen proyectores específicos ($\Pi^S_k$) para los espacios de caras, aristas y nodos. Estos proyectores mapean funciones virtuales a espacios polinomiales utilizando únicamente los DoF de frontera "preservados" y un conjunto específico de momentos internos. Los proyectores se definen mediante condiciones de ortogonalidad que involucran trazas normales/tangenciales y momentos internos (por ejemplo, $(v - \Pi^S_k v, x_K \times p_{kr})_K = 0$).
2. Espacios de Serendipity: Los autores definen los espacios de serendipity ($V^S$) restringiendo los espacios de elementos virtuales estándar de tal manera que los DoF "reducidos" se determinan mediante los DoF preservados a través del proyector de serendipity. Esto asegura que la proyección $L^2$ sobre un espacio polinomial de orden suficientemente alto $P_l$ sea estable y computable utilizando únicamente los DoF de frontera.
3. Complejo de De Rham: El método está diseñado para manejar la cadena completa del complejo de De Rham en $\mathbb{R}^3$ ($H^1 \to H(\text{curl}) \to H(\text{div}) \to L^2$). La construcción asegura una conmutatividad exacta con los operadores diferenciales ($\nabla, \nabla\times, \nabla\cdot$), preservando la estructura topológica del problema.
4. Análisis de Estabilidad: Los autores demuestran que para un orden de proyección $l$ suficientemente grande, la equivalencia de norma $\|\Pi_l v\|_K \approx \|v\|_K$ se cumple para funciones en el espacio de serendipity. Esto elimina la necesidad de términos de estabilización manteniendo propiedades de aproximación óptimas.

Contribuciones Clave

• Tratamiento Completo de De Rham: El trabajo proporciona un tratamiento unificado para todo el complejo de De Rham en 3D, asegurando una integración fluida entre los diferentes operadores diferenciales, una característica ausente en trabajos previos libres de estabilización limitados a $H^1$.
• Preservación de Propiedades Matemáticas: Los espacios de serendipity propuestos mantienen rigurosamente la continuidad de frontera, la conformidad de operador, la conmutatividad exacta con los operadores diferenciales y la compatibilidad B.
• Eficiencia Computacional: Al minimizar el número de DoF y utilizar proyecciones de serendipity computadas únicamente desde los datos de frontera, el método reduce significativamente la carga computacional. Esto es particularmente beneficioso para problemas de autovalores donde los DoF internos no pueden eliminarse mediante condensación estática.
• Análisis de Error Riguroso: El artículo proporciona una prueba teórica de las propiedades de aproximación óptimas y la estabilidad para las proyecciones polinomiales de alto orden.

Resultados
Las afirmaciones teóricas son validadas mediante experimentos numéricos en problemas de autovalores de Maxwell utilizando mallas poliedrales de Voronoi.

• Tasas de Convergencia: Los resultados numéricos confirman las tasas de convergencia teóricas. Para soluciones infinitamente regulares, el método logra una convergencia de $O(h^{m+1})$ en la norma $L^2$ y de $O(h^{2m+2})$ para los autovalores. Para soluciones singulares (donde la regularidad está limitada por $s$), las tasas son $O(h^s)$ y $O(h^{2s})$, respectivamente.
• Costo Computacional: Se muestra que el tiempo de ensamblaje del VEM libre de estabilización (SFVEM) propuesto es solo moderadamente mayor que el del FEM convencional, verificando la eficiencia del método. Se encuentra que el orden de proyección $l$ requerido está cerca del orden polinomial $k$ en la práctica.
• Comparación: El método elimina con éxito las complicaciones relacionadas con la estabilización mientras iguala la precisión del FEM convencional y del VEM estabilizado.

Significancia y Reivindicaciones
Los autores afirman que este trabajo ofrece un enfoque fundamentalmente diferente de la VEM libre de estabilización que supera las limitaciones de las estrategias existentes (exceso de DoF o pérdida de flexibilidad geométrica). Al eliminar los términos de estabilización, el método mejora la estabilidad numérica y los números de condición, lo que lo hace particularmente valioso para:

• Formulaciones mixtas.
• Problemas no lineales.
• Aplicaciones que requieren un comportamiento numérico robusto a través de diversas configuraciones de malla.

El artículo concluye que, si bien la técnica puede generalizarse teóricamente a dimensiones superiores, las complejidades de implementación práctica y el rápido crecimiento de los DoF de frontera limitan actualmente la discusión al entorno 3D. El método propuesto representa una alternativa teóricamente rigurosa y computacionalmente eficiente para resolver ecuaciones en derivadas parciales, avanzando la aplicabilidad de VEM en las matemáticas computacionales y la ingeniería.