The jump filter in the discontinuous Galerkin method for hyperbolic conservation laws

Resumen Técnico: El Filtro de Salto (Jump Filter) en el Método de Galerkin Descontinuo para Leyes de Conservación Hiperbólicas

Planteamiento del Problema
Los esquemas numéricos lineales de alto orden, particularmente el método de Galerkin Descontinuo (DG), se utilizan ampliamente para resolver leyes de conservación hiperbólicas debido a su alta precisión y flexibilidad con mallas no estructuradas. Sin embargo, al simular soluciones que contienen discontinuidades (como ondas de choque), estos esquemas suelen generar oscilaciones espurias indeseables (el fenómeno de Gibbs). Estas oscilaciones pueden provocar la divergencia del resolvedor o la creación de estructuras artificiales. Aunque existen diversas técnicas de captura de choques —como los limitadores de Variación Total Acotada (TVB), los limitadores WENO y los métodos de viscosidad artificial—, muchas sufren de desventajas que incluyen un alto costo computacional, dificultad de implementación en mallas no estructuradas o curvas, la necesidad de descomposición característica o una disipación numérica excesiva que degrada la precisión en regiones suaves.

Metodología
Los autores proponen un novedoso filtro de salto (jump filter) diseñado para operar dentro del marco de trabajo de DG con el fin de mitigar estas oscilaciones preservando la alta precisión y la compacidad del método.

• Formulación de Viscosidad Local: El núcleo del método es la introducción de un término de viscosidad local basado en la información de salto en las interfaces de las celdas. A diferencia de la viscosidad espectral global, este enfoque es localizado. El coeficiente de viscosidad, $\sigma_j$, se determina de forma adaptativa mediante la magnitud de los saltos en la solución y sus derivadas a través de los límites de las celdas.
  • En regiones suaves donde los saltos son mínimos, la viscosidad añadida es insignificante, preservando la alta precisión.
  • En regiones con discontinuidades fuertes, la viscosidad aumenta para amortiguar los modos polinómicos de alto orden responsables de las oscilaciones.
• Filtrado Espectral mediante División Temporal (Time-Splitting): Para evitar la rigidez y el costo computacional asociados con la resolución implícita del término de difusión (lo que restringiría el paso de tiempo), los autores emplean una estrategia de división temporal. Después de cada etapa del esquema DG de Runge-Kutta explícito estándar, el filtro de salto se aplica como un paso de post-procesamiento. Esto implica resolver una ecuación diferencial ordinaria (ODE) simple para los coeficientes modales de la solución polinómica:
  $$\frac{d}{dt} u^l_j + \sigma^l_j u^l_j = 0$$
  La solución se actualiza multiplicando los coeficientes por un factor de decaimiento exponencial, $e^{-\sigma^l_j \Delta t}$. Esta operación es explícita, computacionalmente barata y no requiere descomposición característica.
• Integración de Limitador Híbrido: El filtro de salto también se integra en un marco de limitador híbrido. En este enfoque, se identifican las celdas "problemáticas" utilizando un indicador de salto. Para estas celdas, el filtro de salto se utiliza para calcular una solución polinómica de menor orden, que luego se combina con la solución DG de alto orden. Esto tiene como objetivo reducir aún más la disipación numérica en comparación con el uso del filtro por sí solo.
• Aplicabilidad: El método está formulado tanto para leyes de conservación escalares como para sistemas (específicamente las ecuaciones de Euler) en una y dos dimensiones. Opera directamente en el espacio de las variables físicas (conservativas), evitando las complejidades de la descomposición característica.

Contribuciones Clave

1. Novedoso Filtro de Salto: Introducción de un filtro de viscosidad local que utiliza la información de salto en las interfaces de las celdas para amortiguar adaptativamente los modos de alto orden.
2. Preservación de las Propiedades de DG: El análisis teórico confirma que el esquema propuesto mantiene la conservación, la estabilidad $L^2$ y las tasas de convergencia óptimas (precisión de alto orden) para soluciones suaves.
3. Eficiencia Computacional: El método está diseñado para un bajo costo computacional. No requiere descomposición característica, opera enteramente en el espacio físico y mantiene la compacidad del esquema DG (comunicación de datos local), lo que lo hace altamente adecuado para entornos de computación paralela y GPU.
4. Integración Híbrida: Demostración de la integración fluida del filtro en un marco de limitador híbrido, suprimiendo eficazmente las oscilaciones en choques fuertes mientras minimiza la disipación numérica en regiones suaves.

Resultados Numéricos
El artículo valida el método mediante extensos experimentos numéricos en las ecuaciones de Euler en 1D y 2D, cubriendo:

• Flujo Suave: Las pruebas de precisión (por ejemplo, la ecuación de Burgers, vórtice isentrópico) confirman que el filtro no degrada el orden de convergencia del esquema DG estándar en mallas finas.
• Choques Fuertes: Las pruebas que involucran el problema de Lax, ondas de explosión (blast waves) y reflexión de Mach doble muestran que el filtro controla eficazmente las oscilaciones espurias y captura perfiles de choque nítidos sin los sobreimpulsos típicos de los esquemas de alto orden sin filtrar.
• Interacciones Complejas: El método funciona bien en escenarios de interacción choque-vórtice, interacción choque-burbuja y la inestabilidad de Kelvin-Helmholtz, donde los esquemas de alto orden con el filtro resuelven estructuras de flujo intrincadas (como líneas de deslizamiento y capas de mezcla turbulenta) mejor que sus contrapartes de menor orden.
• Robustez: El esquema demuestra robustez en simulaciones de chorros astrofísicos de alto número de Mach y mantiene la positividad cuando se combina con un limitador de preservación de positividad.

Significancia y Reivindicaciones
Los autores afirman que el filtro de salto ofrece una alternativa robusta, precisa y eficiente a los limitadores tradicionales y a los métodos de viscosidad artificial. Su principal significancia radica en su simplicidad y eficiencia computacional:

• Evita la costosa descomposición característica requerida por muchos otros limitadores de alto orden.
• Preserva la compacidad del método DG, facilitando una paralelización eficiente.
• Proporciona una forma fluida de introducir la disipación necesaria solo donde es necesaria, equilibrando la capacidad de captura de choques con la preservación de la alta precisión en regiones suaves.

El artículo concluye que este enfoque aborda eficazmente el desafío de simular leyes de conservación hiperbólicas con soluciones discontinuas, ofreciendo una herramienta práctica para la captura de choques de alta resolución. Se señala como trabajo futuro la exploración de su implementación en mallas no estructuradas y problemas de estado estacionario.