A collocation scheme that is equivalent to discontinuous Galerkin discretizations

Resumen Técnico: Un Esquema de Colocación Equivalente a Discretizaciones de Galerkin Discontinuo

Enunciado del Problema
El artículo aborda la construcción de discretizaciones estables y de alto orden para ecuaciones diferenciales parciales (EDP) en mallas no estructuradas. Si bien los operadores de Suma por Partes (SBP) proporcionan un marco robusto para la estabilidad, los operadores SBP multidimensionales en elementos simplex (como triángulos) pueden ser complejos de construir y a menudo carecen de expresiones en forma cerrada. Un trabajo anterior de Chan introdujo un operador SBP basado en proyección (denominado aquí como operador de Colocación Modal o MC) como un paso intermedio en discretizaciones de elementos finitos. Sin embargo, las propiedades específicas de utilizar este operador directamente en un esquema de colocación, particularmente su relación con los métodos de Galerkin Discontinuo (DG), no habían sido exploradas completamente. La pregunta central es si una semidiscretización basada en este operador MC específico produce soluciones equivalentes a una discretización modal DG que utilice la misma cuadratura subyacente.

Metodología
El autor define el operador MC utilizando bases polinómicas ortonormales y reglas de cuadratura con puntos interiores positivos. La metodología central incluye:

1. Definición del Operador: Definir el operador MC de primera derivada $D_{xd} = V_{xd}V^T W$, donde $V$ es la matriz de Vandermonde de las funciones base ortonormales, $V_{xd}$ contiene sus derivadas y $W$ es la matriz diagonal de los pesos de cuadratura.
2. Demostración de Equivalencia Teórica: El artículo construye una demostración que muestra que la semidiscretización de un problema de valor inicial y frontera (IBVP) hiperbólico no lineal escalar utilizando el operador MC produce un vector de solución $u(t)$ que es exactamente la proyección del vector de solución modal DG $\tilde{u}(t)$ sobre el espacio nodal ($u(t) = V\tilde{u}(t)$). Esta demostración se basa en descomponer el espacio de soluciones en el subespacio polinómico y el espacio nulo del operador.
3. Análisis Espectral: El autor analiza los autovalores de las matrices del sistema resultante. Se demuestra que el espectro del operador MC consiste en los autovalores del operador DG correspondiente más un conjunto de autovalores nulos asociados al espacio nulo.
4. Manejo de la Inconsistencia del Espacio Nulo: El artículo reconoce que los operadores MC son generalmente inconsistentes en el espacio nulo (su espacio nulo contiene más que solo vectores constantes). Para problemas no estacionarios, estos modos extraños se inicializan en cero y permanecen en cero. Para problemas estacionarios, el autor propone la Estabilización por Proyección Local (LPS) para eliminar la singularidad del espacio nulo sin afectar la exactitud polinómica.
5. Verificación Numérica: La teoría se verifica mediante experimentos numéricos en mallas triangulares utilizando:
   • La ecuación de advección de coeficientes constantes.
   • La ecuación no lineal de Burgers.
   • Comparaciones entre discretizaciones MC "Estándar" y formulaciones "Conservadoras de Entropía" (antisimétricas).
   • Pruebas que involucran reglas de cuadratura con pesos negativos.

Contribuciones Clave

• Teorema de Equivalencia: La contribución principal es la demostración de que una discretización de colocación modal que utiliza el operador MC es matemáticamente equivalente a una discretización modal DG, siempre que compartan la misma regla de cuadratura. Esto se mantiene independientemente del número de grados de libertad (puntos de cuadratura) en el esquema de colocación.
• Propiedades Espectrales: El artículo establece que el radio espectral del operador MC está determinado efectivamente por el método DG subyacente, con la adición de autovalores nulos. Esto implica que la agrupación de nodos (a menudo problemática para los métodos de elementos espectrales en triángulos) no afecta negativamente el radio espectral de los operadores MC de la misma manera que lo hace para los operadores SBP nodales tradicionales.
• Gestión del Espacio Nulo: El trabajo aclara que, aunque los operadores MC son inconsistentes en el espacio nulo, esto no es perjudicial para simulaciones dependientes del tiempo si las condiciones iniciales se proyectan correctamente. Para problemas de estado estacionario, el artículo demuestra cómo recuperar la consistencia mediante disipación numérica basada en proyección (LPS).
• Matiz de Conservación de Entropía: El estudio revela que las discretizaciones MC conservadoras de entropía (basadas en descomposiciones antisimétricas) no son inherentemente equivalentes a DG porque el residuo puede no residir dentro del espacio polinómico. Sin embargo, la equivalencia puede recuperarse proyectando el residuo sobre la base polinómica.

Resultados

• Advección de Coeficientes Constantes: Los resultados numéricos muestran que la diferencia $L_2$ entre las soluciones MC y modal DG es de precisión de máquina ($\approx 10^{-15}$) para diversos grados polinómicos y niveles de exactitud de cuadratura.
• Radios Espectrales: Las distribuciones de autovalores confirman que los operadores MC exhiben radios espectrales favorables en comparación con los operadores SBP de norma mínima, los cuales sufren un crecimiento catastrófico debido a la agrupación de nodos en triángulos.
• Pesos Negativos: La equivalencia se mantiene incluso cuando la regla de cuadratura incluye pesos negativos, un escenario donde la teoría convencional de estabilidad energética de SBP podría fallar en predecir la estabilidad. La estabilidad está garantizada en su lugar por la equivalencia con la solución DG estable energética.
• Ecuación de Burgers: Para la ecuación no lineal de Burgers, la discretización MC "Estándar" es equivalente a DG. La variante "Conservadora de Entropía" muestra inicialmente una discrepancia porque el residuo no está generado por la base polinómica. Sin embargo, una vez que el residuo se proyecta sobre la base, la equivalencia se restablece a precisión de máquina.

Significado y Afirmaciones
El artículo afirma que el operador de Colocación Modal introducido por Chan no es meramente una herramienta intermedia, sino un esquema de discretización válido que es efectivamente un método modal DG. El significado radica en la simplificación de la construcción de esquemas estables y de alto orden en mallas simplex no estructuradas. Al utilizar el operador MC, se puede aprovechar la expresión en forma cerrada del operador mientras se conservan la estabilidad y las propiedades espectrales de los métodos DG.

El autor concluye modestamente que, aunque ambos métodos son equivalentes, la elección entre ellos depende de la aplicación. Para problemas lineales donde las matrices de masa y rigidez se precalculan, DG puede ser preferible. Sin embargo, para problemas no lineales resueltos con esquemas de marcha temporal explícitos o esquemas implícitos libres de matrices, el operador MC puede ofrecer ventajas computacionales. El artículo aboga firmemente por el uso de la colocación modal como el método preferido para construir operadores SBP multidimensionales.