Virtual element approximation of eigenvalue problems: is the stabilization of the right hand side necessary?

1. El Problema

El trabajo aborda la aproximación numérica de problemas de valores propios asociados a operadores elípticos autoadjuntos (como el problema de Laplace) utilizando el Método de Elementos Virtuales (VEM).

En la formulación variacional estándar, se busca encontrar pares $(\lambda, u)$ tales que:
$$a(u, v) = \lambda b(u, v) \quad \forall v \in V$$
donde $a(\cdot, \cdot)$ representa la forma bilineal del operador diferencial (ej. gradiente) y $b(\cdot, \cdot)$ representa la forma bilineal del producto escalar (lado derecho o masa).

El desafío central:
En el VEM, las formas bilineales continuas $a(\cdot, \cdot)$ y $b(\cdot, \cdot)$ no se pueden calcular directamente porque las funciones en el espacio virtual no son polinomios explícitos. Por lo tanto, se utilizan proyecciones polinomiales para construir las formas discretas $a_h$ y $b_h$.

• Para garantizar la estabilidad y la positividad definida de las matrices resultantes ($A$ y $B$), tradicionalmente se añade un término de estabilización a ambas formas.
• Sin embargo, la elección de los parámetros de estabilización es crítica. Una estabilización inadecuada de la matriz $B$ (lado derecho) puede introducir valores propios espurios en la parte baja del espectro, contaminando los modos fundamentales que suelen ser de interés.
• La pregunta de investigación es: ¿Es estrictamente necesario estabilizar la matriz $B$ (el lado derecho) cuando se utilizan espacios VEM estándar de bajo orden?

2. Metodología

Los autores emplean un enfoque teórico basado en la teoría de Descloux–Nassif–Rappaz para la convergencia de operadores no compactos, adaptada a problemas de valores propios.

Marco Abstracto:

1. Definición de Operadores: Se definen el operador solución continuo $T$ y el operador solución discreto $T_h$.
2. Análisis del Núcleo: Reconocen que si la forma $b_h$ no está estabilizada, puede tener un núcleo no trivial ($K_{b,h} \neq \{0\}$), lo que llevaría a valores propios infinitos. Para evitar esto, restringen el análisis al subespacio ortogonal $K_{b,h}^\perp$ con respecto a la forma $a_h$.
3. Condiciones de Convergencia: Utilizan dos propiedades clave para demostrar la convergencia de los valores propios discretos a los continuos:
   • P1: Convergencia del operador solución discreto al continuo ($\|(T - T_h)f\| \to 0$).
   • P2: Propiedad de aproximación del espacio discreto (el espacio $V_h$ debe ser capaz de aproximar las soluciones del problema fuente).

Estrategia de VEM:

• Se considera el problema de Poisson en dominios poligonales.
• Se utiliza la estabilización estándar (tipo "doF-to-doF") únicamente para la forma $a_h$ (lado izquierdo).
• Hipótesis clave: Se elimina la estabilización de la forma $b_h$ (lado derecho), definiendo $b_h$ únicamente mediante la proyección $L^2$ ($\Pi^0_k$).
• Se demuestra teóricamente que, para grados polinomiales bajos ($k=1, 2$), la propiedad P2 se cumple incluso sin estabilizar $B$, debido a una coincidencia especial entre los proyectores de gradiente y $L^2$ en estos órdenes.

3. Contribuciones Clave

1. Prueba Teórica para Bajas Órdenes: El resultado principal es la demostración rigurosa de que, para espacios VEM de orden $k=1$ y $k=2$, la estabilización de la matriz de masa ($B$) no es necesaria. Los valores propios discretos convergen a los valores propios continuos incluso si la matriz $B$ tiene un núcleo no trivial (es decir, si $b_h$ no es estrictamente positiva definida).
2. Análisis del Núcleo: Se caracteriza el núcleo de la matriz $B$ no estabilizada y se demuestra que, al restringir el problema al subespacio ortogonal $K_{b,h}^\perp$, el operador discreto sigue siendo compacto y autoadjunto, preservando la estructura espectral necesaria.
3. Evidencia Numérica para Órdenes Superiores: Aunque la teoría se limita a $k=1, 2$, los autores proporcionan evidencia numérica extensa que sugiere que el método también funciona (sin valores propios espurios y con tasas de convergencia óptimas) para órdenes más altos ($k=3, 4$) en diversas secuencias de mallas.

4. Resultados Numéricos

Los autores realizaron pruebas computacionales utilizando el código VEM desarrollado en la Universidad de Milán-Bicocca. Se probaron cuatro tipos de mallas:

• Triángulos (T)
• Cuadrados (S)
• Voronoi (V)
• Hexágonos (H)
• Mallas Diádicas (D)

Hallazgos principales:

• Dimensiones del Núcleo: Se confirmó que para $k=1$ y $k=2$, el núcleo de la matriz $B$ es trivial (dimensión 0) en mallas triangulares y cuadradas, y a menudo trivial o muy pequeño en mallas Voronoi y hexagonales. Para $k \ge 3$, el núcleo es no trivial en la mayoría de los casos.
• Convergencia Espectral:
  • En todos los casos ($k=1, 2, 3, 4$), los primeros 10 valores propios calculados convergen a los valores exactos con tasas de convergencia óptimas (orden $O(h^{2k})$ para el error en norma $L^2$ y $O(h^k)$ en norma $H^1$, aunque el artículo reporta tasas de error en los valores propios consistentes con la teoría).
  • Ausencia de Modos Espurios: No se observaron valores propios espurios en el espectro, incluso cuando la matriz $B$ tenía un núcleo no trivial (casos con $k \ge 3$).
  • Estabilidad: El método demostró ser robusto sin necesidad de ajustar parámetros de estabilización para el lado derecho.

5. Significado e Impacto

• Simplificación de la Implementación: Eliminar la estabilización del lado derecho simplifica la implementación del VEM para problemas de valores propios, reduciendo la carga computacional y la complejidad del código.
• Eliminación de la Sintonización de Parámetros: Resuelve el problema de la "sintonización" (tuning) de los parámetros de estabilización en $B$, que es una tarea delicada y propensa a errores. Si $B$ no necesita estabilización, se evita el riesgo de contaminar el espectro con valores espurios debido a una mala elección de parámetros.
• Validación de Enfoques "Stabilization-Free": El trabajo apoya la tendencia hacia métodos VEM que no requieren estabilización o que son auto-estabilizantes, demostrando que para problemas de valores propios elípticos, la estabilización de la masa no es un requisito teórico para la convergencia en los órdenes estándar.
• Generalidad: Los resultados se mantienen válidos en mallas poligonales generales (Voronoi, hexagonales), lo cual es una ventaja distintiva del VEM frente a otros métodos de elementos finitos.

En conclusión, el artículo establece que, para la aproximación de problemas de valores propios elípticos con VEM estándar, la estabilización de la matriz de masa (lado derecho) es innecesaria, garantizando la convergencia correcta y la ausencia de modos espurios, especialmente para los órdenes bajos más comunes en la práctica.