Stability analysis of Arbitrary-Lagrangian-Eulerian ADER-DG methods on classical and degenerate spacetime geometries

Resumen Técnico: Análisis de Estabilidad de Métodos ADER-DG de tipo Arbitrario-Lagrangiano-Euleriano

Planteamiento del Problema
El marco de trabajo Arbitrario-Lagrangiano-Euleriano (ALE) se utiliza ampliamente para resolver ecuaciones diferenciales parciales (EDP) hiperbólicas en mallas móviles, particularmente cuando ocurren cambios de topología. En este contexto, los métodos ALE directos conectan mallas en diferentes pasos de tiempo mediante volúmenes de control espacio-temporales. Desarrollos recientes de Gaburro et al. introdujeron elementos espacio-temporales "degenerados" (elementos de tipo sliver) para manejar cambios de topología donde las celdas pueden tener un tamaño espacial nulo al inicio o al final de un paso de tiempo, pero poseen un volumen espacio-temporal no nulo. Si bien los métodos ADER-DG (Arbitrary high-order DERivatives Discontinuous Galerkin) son efectivos para estos problemas, ha faltado un análisis de estabilidad teórica riguroso para estos esquemas —específicamente en geometrías espacio-temporales degeneradas—. La literatura existente suele depender de límites empíricos de CFL (Courant-Friedrichs-Lewy) o estudios de estabilidad limitados para órdenes polinómicos bajos, dejando un vacío en la comprensión de las restricciones de estabilidad para métodos ALE ADER-DG explícitos e implícitos de alto orden tanto en geometrías clásicas como degeneradas.

Metodología
Los autores emplean un análisis de estabilidad de von Neumann para investigar la estabilidad de los métodos ALE ADER-DG explícitos e implícitos. El estudio se centra en la ecuación de advección lineal como un sustituto (proxy) para sistemas hiperbólicos.

1. Formulación: El artículo establece primero el marco matemático para ambos esquemas ALE ADER-DG, explícito e implícito.
   • Esquema Explícito: Utiliza un enfoque de predictor-corrector. El paso del predictor construye una aproximación polinómica local en el espacio-tiempo dentro de cada volumen de control (o elemento sliver) utilizando una iteración de punto fijo. El paso del corrector actualiza la solución utilizando funciones de base móviles y flujos numéricos (de tipo Rusanov) en las interfaces.
   • Esquema Implícito: Resuelve el polinomio espacio-temporal global simultáneamente en todos los volúmenes de control mediante una iteración Newton-GMRES, integrando directamente la EDP sobre los volúmenes de control espacio-temporales.
2. Modelado de Geometría Degenerada: Para analizar geometrías degeneradas, los autores introducen un entorno subrogado (1D) donde las interfaces estándar entre volúmenes de control son reemplazadas por "elementos sliver". Estos elementos imitan los slivers de tipo hueco utilizados en cambios de topología 2D/3D, los cuales tienen un ancho espacial cero en $t_n$ y $t_{n+1}$, pero un volumen espacio-temporal no nulo.
3. Análisis de Estabilidad:
   • Geometrías Clásicas: Los autores derivan matrices de amplificación para los esquemas explícito e implícito. Computan el radio espectral $\rho$ de estas matrices sobre un rango de números CFL y ángulos de fase.
   • Geometrías Degeneradas: El análisis se extiende para incluir elementos sliver. El dominio se modela como un bloque periódico que contiene un volumen de control estándar y un elemento sliver. Se construye una matriz de amplificación dependiente del número CFL, el parámetro de ancho del sliver $\delta$ y el ángulo de fase.
   • Verificación Numérica: Los límites de estabilidad teóricos se validan computando factores de amplificación discretos sobre mallas finas de valores de CFL y grados polinómicos ($N=1$ a $9$). Los órdenes de consistencia también se verifican numéricamente.

Contribuciones Clave y Resultados

• Refinamiento de los Límites de Estabilidad Explícitos: Para geometrías clásicas, el estudio confirma que para grados polinómicos bajos ($N \le 3$), los límites de estabilidad coinciden con los límites de CFL empíricos ampliamente utilizados. Sin embargo, para órdenes superiores ($N \ge 4$), el artículo demuestra que los valores de CFL empíricos comúnmente usados en realidad violan las condiciones estrictas de estabilidad de von Neumann (donde $\rho > 1$). Los autores identifican límites de CFL rigurosos significativamente más bajos requeridos para una estabilidad estricta en esquemas explícitos de alto orden.
• Restricciones de Velocidad de la Malla: El análisis proporciona una caracterización teórica del rango admisible de velocidades de la malla para un CFL objetivo fijo, mostrando cómo el movimiento de la malla afecta las restricciones de estabilidad.
• Estabilidad Incondicional Implícita: El análisis de von Neumann y una prueba teórica (Teorema 1) confirman que el método ALE ADER-DG implícito es incondicionalmente estable para la ecuación de advección lineal en geometrías clásicas, independientemente del tamaño del paso de tiempo.
• Estabilidad en Geometrías Degeneradas:
  • Caso Explícito: La introducción de elementos sliver no degrada la estabilidad. Se encuentra que el factor de amplificación discreto para el entorno degenerado es menor o igual al del entorno clásico. Por consiguiente, los mismos límites de CFL aplicables a las geometrías clásicas son válidos para las degeneradas. Los autores señalan que el tratamiento implícito del paso del predictor sliver dentro del esquema global explícito puede aumentar ligeramente el límite de CFL aceptable.
  • Caso Implícito: El método implícito permanece incondicionalmente estable incluso en presencia de elementos sliver, siempre que el parámetro de ancho del sliver $\delta$ satisfaga una restricción geométrica específica relativa al paso de tiempo y al tamaño de la malla (Teorema 2).
• Consistencia: Experimentos numéricos confirman que tanto los esquemas explícito como el implícito mantienen su orden de consistencia esperado ($N+1$) tanto en geometrías clásicas como degeneradas.

Significancia y Reivindicaciones
El artículo afirma llenar un vacío teórico crítico al proporcionar el primer análisis de estabilidad de von Neumann riguroso para métodos ALE ADER-DG en geometrías espacio-temporales degeneradas. La significancia principal radica en validar el uso de elementos degenerados (slivers) para conectar mallas móviles con cambios de topología. Los resultados demuestran que el uso de estos elementos de tamaño cero no impone restricciones de estabilidad adicionales ni reduce el paso de tiempo permisible en comparación con las geometrías clásicas. Este hallazgo proporciona una base teórica para la aplicación práctica de los métodos ALE directos en escenarios complejos que involucran cambios de topología y allana el camino para el desarrollo de nuevos métodos basados en celdas de corte (cut cell) espacio-temporales. Los autores enfatizan que, si bien en la práctica se suelen utilizar valores de CFL empíricos (mitigados por viscosidad y limitadores), los límites rigurosos derivados son esenciales para una comprensión teórica completa del comportamiento del método.