An Interface Green's Function Framework for Complete Discrete $W^{1,\infty}$ Analysis of Discontinuous Galerkin Methods

Resumen Técnico: Un Marco de Funciones de Green de Interfaz para un Análisis Completo de $W^{1,\infty}$ Discreto de Métodos de Galerkin Descontinuo

Planteamiento del Problema
El análisis de error puntual de los métodos de Galerkin Descontinuo (DG) ha sido un área de investigación significativa, donde la literatura existente establece estimaciones de error óptimas en la norma $L^\infty$ y en la seminorma $W^{1,\infty}$ rota. Sin embargo, persiste una brecha teórica crítica: la semínorma $W^{1,\infty}$ rota, definida como $\|\nabla(u - u_h)\|_{L^\infty_h(\Omega)}$, no logra controlar las discontinuidades a través de las interfaces de la malla. En consecuencia, los análisis existentes no proporcionan un control directo sobre los saltos de interfaz, los cuales son intrínsecos a las aproximaciones discontinuas y vitales para la estabilidad, la conservación y el posprocesamiento. Además, en dominios poliédricos convexos, donde los principios máximos discretos débiles no están disponibles actualmente para los métodos DG, derivar estimaciones óptimas sin factores logarítmicos subóptimos (por ejemplo, $|\ln h|^\beta$) ha sido un desafío persistente, particularmente para grados polinómicos superiores a uno.

Metodología
Los autores proponen un novedoso marco de funciones de Green de interfaz para lograr un análisis completo de $W^{1,\infty}$ discreto. La innovación metodológica central se aparta de los análisis clásicos de funciones de Green, que típicamente utilizan una única función delta discreta para estimar la solución o su gradiente. En su lugar, este marco introduce funciones de Green de interfaz generadas por la diferencia entre funciones delta discretas adyacentes.

1. Definición de Funciones de Green de Interfaz:

   • Para interfaces de malla interiores, el marco define una función de Green $\Gamma_{if}$ que satisface $-\Delta \Gamma_{if} = \delta_{h, K_+^*} - \delta_{h, K_-^*}$, donde $\delta_{h, K}$ representa una función delta discreta suave soportada en el elemento $K$. Este término fuente representa naturalmente el salto de las aproximaciones discontinuas a través de la interfaz.
   • Para caras de frontera, se define una función de Green $\Gamma_{bf}$ similar utilizando una única función delta discreta.

2. Fórmulas de Representación:
   Los saltos de las soluciones numéricas (y de sus gradientes) se representan utilizando estas funciones de Green de interfaz y sus aproximaciones DG. Esto reduce la estimación de las cantidades de interfaz al error de aproximación de las propias funciones de Green de interfaz.

3. Ingrediente Analítico Clave: Estimaciones de Energía Local:
   La demostración se basa en una nueva estimación de energía local para estas funciones de Green de interfaz, establecida mediante la descomposición diedica del dominio.

   • Los autores emplean un argumento de expansión de Taylor con respecto a la variable fuente.
   • Crucialmente, la diferencia entre las funciones delta discretas vecinas ($\delta_{h, K_+^*} - \delta_{h, K_-^*}$) crea un efecto de cancelación a nivel de la fuente singular.
   • Esta cancelación elimina la contribución singular dominante, produciendo estimaciones locales más agudas (específicamente, un factor adicional de $h$) en comparación con las estimaciones clásicas.
   • Esta mejora permite la derivación de estimaciones de interfaz sin introducir factores logarítmicos que típicamente surgen de la aplicación directa de desigualdades de traza discretas combinadas con las estimaciones de norma máxima existentes.

Contribuciones Clave y Resultados
El artículo establece una teoría de $W^{1,\infty}$ discreta "completa", definida por la norma:
$$|||v|||_{W^{1,\infty}_h(\Omega)} = \max_{K \in \mathcal{T}_h} \|\nabla v\|_{L^\infty(K)} + \max_{F \in \mathcal{E}_h} h_F^{-1} \|[v]\|_{L^\infty(F)}$$
Esta norma mide simultáneamente el gradiente roto y los saltos de interfaz escalados.

El marco se aplica a los siguientes métodos en dominios poliédricos convexos para la ecuación de Poisson:

• Método de Penalización Interior Simétrica (IPDG).
• Método Híbrido Mixto (Hybrid Mixed DG).
• Dos métodos representativos de DG Hibridizable (HDG-1 y HDG-2).

Principales Resultados Teóricos:

• Estimaciones de Salto Óptimas: Los autores demuestran estimaciones óptimas en norma máxima para los saltos de interfaz $[u - u_h]$ y los flujos numéricos, eliminando el factor logarítmico $\beta$ que aparecía anteriormente en las estimaciones en dominios poliédricos convexos.
• Límites de Error Completos: Derivan estimaciones de error óptimas para los gradientes rotos y los saltos de interfaz simultáneamente, completando el análisis puntual para estos métodos.
• Extensión a Problemas de Cuarto Orden: El marco se extiende al método de penalización interior $C^0$ para la ecuación biarmónica. Los autores derivan nuevas estimaciones de norma máxima para los saltos de la derivada normal, una cantidad que no había sido tratada previamente en el análisis puntual para este método.

Significancia y Reivindicaciones
El artículo afirma que su principal contribución es metodológica. Al introducir funciones de Green de interfaz basadas en la diferencia de funciones delta discretas, los autores superan la barrera fundamental de la falta de principios máximos discretos débiles para los métodos DG en dominios poliédricos convexos.

• Eliminación de Factores Logarítmicos: El marco elimina con éxito los factores logarítmicos subóptimos heredados de los enfoques clásicos de desigualdad de traza, proporcionando estimaciones óptimas para grados polinómicos mayores que uno.
• Generalidad: La construcción depende de la estructura discontinua local de la aproximación más que de la formulación numérica específica. Los autores afirman que la metodología es suficientemente general como para ser aplicable a una amplia clase de métodos de elementos finitos discontinuos y no conformes.
• Completitud: El trabajo proporciona el primer análisis de $W^{1,\infty}$ discreto completo que controla simultáneamente los gradientes rotos y los saltos de interfaz para los métodos DG especificados, llenando un vacío crítico en la comprensión teórica de las aproximaciones puntuales para métodos discontinuos.