Invariant domain preserving limiting of time explicit and time implicit discretizations for systems of conservation laws

Resumen Técnico: Limitación que Preserva Dominios Invariantes para Discretizaciones Temporales Explícitas e Implícitas

Enunciado del Problema
El artículo aborda el desafío de diseñar aproximaciones numéricas de alto orden (HO) para sistemas hiperbólicos no lineales de leyes de conservación que preserven dominios invariantes (ID). Si bien los métodos HO capturan con precisión las soluciones suaves, a menudo no satisfacen las restricciones físicas a nivel discreto, como la positividad de la densidad y la energía interna, o el cumplimiento de las desigualdades de entropía. Por el contrario, los esquemas de bajo orden (LO) (por ejemplo, Godunov, Lax-Friedrichs) son conocidos por preservar estos dominios invariantes, pero sufren de una difusión numérica excesiva. El objetivo es construir una técnica de limitación que restrinja la solución HO para que permanezca dentro del conjunto convexo de estados admisibles $D$ (por ejemplo, densidad positiva, energía interna positiva, límites de entropía) sin sacrificar la conservación ni la precisión de alto orden en regiones suaves.

Metodología
Los autores proponen un marco general de limitación aplicable a cualquier discretización espacial conservativa (Volumen Finito, Galerkin Discontinuo) y a una amplia gama de esquemas de integración temporal, incluidos los métodos de Runge-Kutta (RK) explícitos e implícitos, métodos multietapa y la integración temporal Galerkin Discontinua (DG).

1. Concepto Central: El método limita la solución HO alrededor de una solución LO que está garantizada como preservadora de ID. Generaliza el marco de Transporte Corregido por Flujo (FCT) limitando los flujos antidifusivos (la diferencia entre los residuos HO y LO), pero define los coeficientes de limitación para expresar la solución final como una combinación convexa de cantidades que preservan ID, similar al marco de limitación convexa.
2. Coeficientes de Limitación: A diferencia del FCT estándar, que impone límites a los componentes de la solución, este enfoque define coeficientes $l_{\kappa,\kappa'}^{i,j}$ tales que la solución limitada $u_h^{n+1}$ es una combinación convexa de la solución LO y la actualización HO. Los coeficientes se determinan encontrando los mayores factores de escala $k$ tales que los estados intermedios permanezcan dentro del dominio invariante $D$.
3. Procedimiento Iterativo: El limitador se integra en un algoritmo iterativo. Partiendo de la solución LO, el método refina iterativamente la solución limitada re-aplicando el limitador a la diferencia entre la solución HO y la solución limitada actual. Este proceso continúa hasta cumplir una tolerancia o alcanzar un número máximo de iteraciones, con el objetivo de recuperar la solución HO lo más cerca posible mientras se mantienen las propiedades de preservación de ID.
4. Aceleración: Se introduce un factor de aceleración heurístico $\beta \geq 1$ para escalar los coeficientes de limitación, acelerando la convergencia del proceso iterativo hacia la solución HO.
5. Adaptación a la Integración Temporal: El marco se adapta a diversos esquemas temporales definiendo un esquema LO de un solo paso (Euler hacia adelante para explícito, Euler hacia atrás para implícito) para calcular los residuos necesarios y los flujos antidifusivos, evitando la necesidad de recalcular soluciones LO en cada etapa de métodos multietapa o multivariados.
6. Discretización Espacial: Se proporcionan derivaciones específicas para flujos antidifusivos para esquemas de Volumen Finito (FV), esquemas DG generales y el Método de Elementos Espectrales Galerkin Discontinuo (DGSEM). Para DGSEM, el método combina la limitación basada en promedios de celda con un limitador de escala lineal para garantizar que los valores nodales también satisfagan las restricciones del dominio invariante.

Contribuciones Clave

• Generalización a Sistemas y Esquemas Implícitos: El método extiende los conceptos de limitación convexa a sistemas de leyes de conservación y, crucialmente, se aplica a discretizaciones temporales implícitas, un dominio donde los enfoques previos de limitación convexa estaban restringidos a esquemas explícitos.
• Limitación Basada en Flujos sin Estados Bar: El enfoque define coeficientes de limitación para formar combinaciones convexas de soluciones LO y HO sin depender de "estados bar" (estados intermedios utilizados en algunas formulaciones de limitación convexa). Esto permite el uso de cualquier flujo numérico consistente y conservativo en el esquema HO.
• Refinamiento Iterativo: La introducción de un procedimiento iterativo permite que la solución limitada converja hacia la solución HO, recuperando precisión en regiones suaves mientras preserva estrictamente los dominios invariantes.
• Heurística de Aceleración: Se propone un simple factor de escala $\beta$ para reducir significativamente el número de iteraciones requeridas para la convergencia.
• Amplia Aplicabilidad: La derivación se presenta para discretizaciones conservativas generales, con implementaciones específicas detalladas para FV, DG y DGSEM, acopladas con integradores temporales RK explícitos/implícitos, multietapa y DG.

Resultados Numéricos
El artículo valida el método mediante experimentos numéricos en problemas 1D y 2D que involucran ecuaciones hiperbólicas escalares y las ecuaciones de Euler compresibles:

• Problemas de Riemann (Ondas de Choque Fuertes): Para las ecuaciones de Euler compresibles, el esquema HO sin limitación falla debido a presiones/densidades negativas. El limitador propuesto resuelve con éxito las ondas de choque fuertes de manera nítida, manteniendo la positividad y los límites de entropía tanto para la integración temporal explícita (ERK3) como implícita (DIRK33), incluso con números CFL grandes (hasta 10).
• Inestabilidad de Kelvin-Helmholtz: En simulaciones 2D con grandes pasos de tiempo, el esquema LO amortigua excesivamente las escalas pequeñas. El esquema limitado recupera la alta resolución del esquema HO, preservando estructuras de escala fina.
• Difracción de Choques y Flujo sobre Escalones: El limitador demuestra robustez en geometrías 2D complejas, manteniendo la admisibilidad física donde el esquema HO ilimitado falla. Logra una resolución comparable a la de esquemas explícitos con pasos de tiempo pequeños, pero con la estabilidad de los métodos implícitos.
• Ecuaciones Escalares: Las pruebas en transporte lineal, la ecuación de Burgers y el problema KPP (flujo no convexo) muestran que el limitador elimina eficazmente las oscilaciones espurias y hace cumplir el principio del máximo y las desigualdades de entropía.
• Convergencia: Se demuestra que el limitador iterativo con el factor de aceleración $\beta=2$ converge a la solución HO significativamente más rápido que el enfoque estándar ($\beta=1$).

Significado y Afirmaciones
Los autores afirman que este trabajo proporciona un limitador bastante general capaz de manejar cualquier ley de conservación hiperbólica en mallas generales con una amplia gama de esquemas de integración temporal. El significado principal radica en cerrar la brecha entre la robustez de los esquemas implícitos de bajo orden y la precisión de los métodos de alto orden. Al garantizar que el esquema numérico preserve los dominios invariantes a nivel discreto, el método permite el uso de pasos de tiempo grandes en simulaciones implícitas sin el riesgo de producir soluciones no físicas (por ejemplo, densidad negativa). El artículo enfatiza que el limitador mantiene la conservación y puede aplicarse iterativamente para refinar la solución, ofreciendo una herramienta práctica para simulaciones robustas de alto orden de problemas complejos de dinámica de fluidos. El trabajo concluye señalando que las extensiones futuras considerarán la integración temporal Euler hacia atrás linealizada para problemas de estado estacionario.