A posteriori error analysis and adaptivity of a space-time finite element method for the wave equation in second order formulation

A continuación se presenta un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Título: Análisis de error a posteriori y adaptividad de un método de elementos finitos espacio-temporal para la ecuación de onda en formulación de segundo orden.

Autores: Zhaonan Dong, Emmanuil H. Georgoulis, Lorenzo Mascotto y Zuodong Wang.

---

1. Problema Planteado

El trabajo se centra en la aproximación numérica de la ecuación de onda lineal en su formulación de segundo orden:
$$u'' - \Delta u = f$$
en un dominio espacial $D$ y un intervalo temporal $J=(0, T]$, con condiciones de contorno de Dirichlet homogéneas y condiciones iniciales dadas.

Desafíos principales:

• Las soluciones de ecuaciones hiperbólicas de segundo orden suelen presentar características espacio-temporales altamente localizadas (frentes de onda, pulsos, interacciones).
• Los métodos estándar basados en mallas cuasi-uniformes requieren un número excesivo de grados de libertad para capturar estas características en dominios grandes y tiempos largos.
• Los métodos explícitos sufren restricciones de estabilidad tipo CFL, mientras que los métodos implícitos espacio-temporales, aunque permiten mallas variables, carecían hasta ahora de análisis de error a posteriori riguroso para órdenes arbitrarios con modificación dinámica de la malla.

2. Metodología

Los autores proponen y analizan un método de elementos finitos espacio-temporal (ST-FEM) de orden arbitrario.

• Discretización Espacial: Se utilizan elementos finitos conformes $H^1$ de grado $p$ sobre mallas simpliciales o de cajas.
• Discretización Temporal:
  • Se emplean polinomios por tramos continuos en el tiempo de grado $q$.
  • Para la segunda derivada temporal, se utiliza una aproximación de tipo Galerkin Discontinuo (DG) con flujo hacia arriba (upwind).
  • Característica clave: El método permite la modificación dinámica de la malla en cada paso de tiempo. Cuando la malla cambia entre intervalos temporales consecutivos, la aproximación se vuelve discontinua en el tiempo en los nodos de tiempo, lo que requiere un tratamiento especial.
• Proyección entre mallas: Cuando la malla cambia de $\mathcal{T}^{n-1}_h$ a $\mathcal{T}^n_h$, la solución se proyecta del espacio anterior al actual mediante un operador de proyección ortogonal $P^n_p$ (puede ser interpolación lagrangiana, proyección $L^2$ o $H^1$).

3. Contribuciones Clave

A. Límites de Error A Posteriori Rigurosos

El principal aporte teórico es la demostración de límites superiores de error a posteriori completamente computables en la norma $L^\infty(L^2)$ (máximo en el tiempo de la norma $L^2$ espacial).

• Reconstrucciones No Estándar: Para manejar la discontinuidad temporal causada por el cambio de malla, los autores diseñan reconstrucciones temporales y espaciales novedosas:
  1. Reconstrucción Temporal (Tipo Hermite): Utiliza una expansión de base de Hermite para promediar los coeficientes en los nodos temporales, asegurando la continuidad de la primera derivada temporal y eliminando la discontinuidad de la solución, manteniendo el control del error.
  2. Reconstrucción Elíptica: Se utiliza una reconstrucción elíptica no estándar para manejar el error espacial, motivada por desarrollos previos en problemas parabólicos.
• Control de Constantes: A diferencia de trabajos anteriores, este análisis proporciona constantes explícitas en los límites de error, incluyendo la dependencia de los órdenes espaciales ($p$) y temporales ($q$).
• Función de Prueba de Baker: Se utiliza una función de prueba inspirada en la construcción clásica de Baker para acotar el error en la norma $L^\infty(L^2)$.

B. Estimadores de Error Desglosados

El error total se descompone en contribuciones computables:

1. Error Inicial ($\eta_{init}$): Depende de la aproximación de los datos iniciales.
2. Oscilación de la Fuente ($\eta_f$): Relacionada con la aproximación del término forzante $f$.
3. Error Temporal ($\eta_{time}$): Asociado a la discretización temporal y la proyección de la derivada temporal.
4. Error Espacial ($\eta_{space}$): Incluye estimadores de tipo residual y de salto en las facetas.
5. Error por Cambio de Malla ($\eta_{mesh}$): Un término nuevo y crucial que cuantifica el error introducido específicamente por la modificación de la malla entre pasos de tiempo.

C. Algoritmo Adaptativo Espacio-Temporal

Se propone un algoritmo adaptativo completo que utiliza los estimadores anteriores para:

• Refinar o coarsenar (acercar) la malla espacial localmente.
• Ajustar el tamaño del paso de tiempo ($\tau_n$) dinámicamente.
• Utilizar la estrategia de marcado de Dörfler para seleccionar elementos.

4. Resultados Numéricos

• Disipación de Energía: Se analizó la disipación de energía del método.
  • En mallas estáticas, la disipación es exponencialmente pequeña para órdenes temporales altos ($q$).
  • En mallas dinámicas, la proyección de la solución entre mallas introduce disipación adicional, pero se mantiene bajo control. Se observó que aumentar el orden temporal $q$ reduce significativamente la disipación.
• Eficiencia de los Estimadores:
  • Los índices de efectividad (relación entre el estimador y el error real) se mantienen acotados y no crecen desmesuradamente con el refinamiento de la malla, a diferencia de otros métodos reportados en la literatura.
  • Los estimadores capturan correctamente la convergencia óptima tanto en versiones $h$ (refinamiento de malla) como $p$ (aumento de grado polinomial).
• Efecto del Cambio de Malla:
  • Se demostró que cambios masivos de malla pueden degradar la tasa de convergencia a $O(\tau^q)$ (pérdida de un orden de convergencia) si no se controlan.
  • La proyección $H^1$ se identificó como superior a la interpolación lagrangiana o la proyección $L^2$ para transferir la solución entre mallas, minimizando el error de proyección y el término $\eta_{mesh}$.
• Algoritmo Adaptativo:
  • Las pruebas numéricas con soluciones suaves y singulares (ondas con discontinuidades) mostraron que el algoritmo adaptativo distribuye eficientemente los grados de libertad, concentrando la malla en las regiones de alta actividad de la onda.
  • El método logra una estabilidad a largo plazo superior a la de los métodos con malla uniforme, donde el error crece linealmente con el tiempo de simulación.

5. Significado e Impacto

• Rigor Teórico: Este trabajo cierra una brecha importante en la literatura al proporcionar el primer análisis de error a posteriori riguroso y con constantes explícitas para métodos de Galerkin espacio-temporales de orden arbitrario que admiten cambio dinámico de malla.
• Eficiencia Computacional: Al permitir la adaptación de la malla en tiempo y espacio, el método reduce drásticamente el costo computacional necesario para simular ondas en dominios grandes durante tiempos largos, evitando el uso de mallas uniformemente finas.
• Aplicabilidad: La metodología es extensible a otros problemas hiperbólicos de segundo orden, ecuaciones de onda con coeficientes variables, problemas amortiguados y otros operadores elípticos lineales.
• Control de Calidad: La capacidad de estimar el error introducido por el cambio de malla ($\eta_{mesh}$) es fundamental para garantizar la fiabilidad de las simulaciones adaptativas en problemas de propagación de ondas.

En resumen, el artículo establece un marco teórico y práctico sólido para la simulación eficiente y precisa de la ecuación de onda mediante métodos espacio-temporales adaptativos, superando las limitaciones de los métodos tradicionales en cuanto a la gestión de la complejidad computacional y el control de errores.