$hp$-a posteriori error estimates for hybrid high-order methods applied to biharmonic problems

Planteamiento del Problema
Este artículo aborda la estimación de error a posteriori para los métodos Híbridos de Alto Orden (HHO) aplicados al problema biarmónico ($\Delta^2 u = f$) con condiciones de contorno de empotramiento ($u = \partial_n u = 0$ en $\partial\Omega$). El problema se plantea en dominios politópicos de Lipschitz en dos y tres dimensiones con una topología general (que potencialmente contiene agujeros). El objetivo es derivar estimadores de error basados en residuos que sean robustos con respecto al tamaño de la malla ($h$), el grado polinómico ($p$ o $k$) y la topología del dominio.

Metodología
Los autores desarrollan un marco para acotar el error total descomponiéndolo en componentes conformes y no conformes.

1. Descomposición del Error: El error se divide en una parte conforme (medida en la norma de energía de una función reconstruida en $H^2_0$) y una parte no conforme (la distancia entre la solución discreta y el espacio conforme).
2. Control del Error No Conforme:
   • Para acotar el error no conforme, los autores construyen una reconstrucción conforme en $H^2_0$ de la solución discreta utilizando una partición de la unidad $C^1$.
   • Esta partición se construye utilizando espacios de elementos finitos conformes en $C^1$ (específicamente, la construcción de [40] basada en divisiones de Alfeld de la malla).
   • Crucialmente, los autores emplean descomposiciones de Helmholtz locales en parches de vértices (estrellas). A diferencia de enfoques previos que dependían de descomposiciones de Helmholtz globales (que introducen constantes que dependen del número de agujeros en el dominio), este enfoque local asegura que la constante de estabilidad dependa únicamente de la regularidad de la forma de la malla.
3. Control del Error Conforme:
   • Se derivan dos límites superiores distintos para el error conforme, cada uno de los cuales utiliza un operador de interpolación específico:
     • Configuración 1: Utiliza el operador de interpolación de elementos finitos híbridos canónicos ($C^{k+2}_T$). Este enfoque se basa en la propiedad de que la composición de este operador con la reconstrucción HHO produce una proyección elíptica $H^2$. El límite resultante involucra solo el término de estabilización y la oscilación de los datos, pero depende de un supuesto numérico respecto al error de interpolación (Supuesto 3.4).
     • Configuración 2: Utiliza el operador de interpolación $hp$ modificado de Babuška–Suri ($I^{k+2}_{mBS}$). Este enfoque no requiere supuestos adicionales, pero introduce términos de indicadores de error adicionales, específicamente residuos de volumen y saltos de flujo normal y tangencial.
4. Estimadores de Residuos: El artículo define indicadores de error locales que comprenden términos de estabilización, residuos de volumen, saltos de flujo normal y saltos tangenciales, escalados por factores $hp$ apropiados.

Contribuciones Clave

• Primer estimador $hp$ para HHO en problemas biarmónicos: Este trabajo proporciona el primer estimador de error basado en residuos $hp$ para métodos HHO aplicados al problema biarmónico en dominios generales de 2D y 3D.
• Descomposición de Helmholtz Local con Robustez Topológica: El artículo establece un nuevo estimulo de estabilidad para descomposiciones de Helmholtz locales en estrellas de vértices. Demuestra que la constante de estabilidad depende únicamente de la regularidad de la forma de la malla, evitando así la degradación de las constantes asociadas con el número de agujeros en el dominio, un problema conocido en los enfoques de descomposición global.
• Límites Superiores Duales: La derivación de dos límites superiores distintos para el error conforme ofrece flexibilidad. Un límite es "limpio" (depende solo de la estabilización y la oscilación) pero condicional a un supuesto verificado numéricamente; el otro es incondicional pero incluye términos de residuo más complejos.
• Extensión a la Topología General: La metodología está diseñada explícitamente para manejar dominios con topología arbitraria (múltiples agujeros), una característica difícil de abordar frecuentemente en la estimación $hp$ para problemas de cuarto orden.

Resultos

• Límites Teóricos: El teorema principal (Teorema 5.6) establece que el error total (norma de Hessiano roto más estabilización) está acotado por la suma de los indicadores de error propuestos. Se demuestra que los límites son $h$-óptimos y $p$-subóptimos (por un orden como máximo, o potencialmente medio orden bajo el supuesto verificado).
• Verificación Numérica:
  • Verificación del Supuesto: Experimentos numéricos en 2D y 3D confirman la validez del Supuesto 3.4 respecto al error de interpolación del operador híbrido canónico.
  • Tasas de Convergencia: Pruebas en soluciones suaves en 3D demuestran que tanto el error de energía como el estimador convergen a la tasa óptima $O(\text{DoFs}^{-(k+1)/d})$.
  • Índices de Efectividad: Los estimadores arrojan índices de efectividad moderados (que oscilan aproximadamente entre 1.2 y 2.0) a través de varios grados polinómicos y refinamientos de malla.
  • Robustez Topológica: Experimentos en dominios 2D con un número creciente de agujeros (0, 1, 4, 8) muestran que los índices de efectividad permanecen estables, confirmando la afirmación teórica de que el estimador es robusto con respecto a la topología del dominio.
  • Desempeño Adaptativo: Algoritmos de refinamiento $h$-adaptativo impulsados por el estimador recuperan con éxito las tasas de convergencia óptimas para soluciones singulares tanto en 2D (dominio en forma de L) como en 3D (singularidad puntual), superando al refinamiento uniforme.

Significancia
Los autores posicionan este trabajo como el primero en emplear una partición de la unidad $C^1$ combinada con descomposiciones de Helmholtz locales en estrellas de vértices para acotar errores no conformes en el entorno $hp$ para problemas biarmónicos. La significancia principal radica en la capacidad de controlar el error no conforme sin incurrir en constantes de estabilidad que se degradan con la complejidad topológica del dominio (número de agujeros). Además, los autores señalan que, aunque el enfoque se desarrolla para métodos HHO, la técnica para acotar errores no conformes mediante descomposiciones de Helmholtz locales es potencialmente aplicable a otros métodos no conformes, como los métodos de Penalización Discontinua de Interior (IPDG). El artículo no pretende proporcionar límites inferiores de error, señalando que tales límites seguirían de adaptaciones de argumentos existentes pero probablemente exhibirían una suboptimidad $p$ pesimista debido a las limitaciones de las técnicas $hp$ actuales para construir funciones burbuja y operadores de extensión específicos.