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

Enunciado do Problema
Este artigo aborda a estimativa de erro a posteriori para métodos Híbridos de Alta Ordem (HHO) aplicados ao problema biarmônico ($\Delta^2 u = f$) com condições de contorno engastadas ($u = \partial_n u = 0$ em $\partial\Omega$). O problema é formulado em domínios politópicos de Lipschitz em duas e três dimensões com topologia geral (potencialmente contendo buracos). O objetivo é derivar estimadores de erro $hp$ baseados em resíduos que sejam robustos em relação ao tamanho da malha ($h$), ao grau polinomial ($p$ ou $k$) e à topologia do domínio.

Metodologia
Os autores desenvolvem um arcabouço para limitar o erro total decompondo-o em componentes conformes e não conformes.

1. Decomposição do Erro: O erro é dividido em uma parte conforme (medida na norma de energia de uma função reconstruída em $H^2_0$) e uma parte não conforme (a distância entre a solução discreta e o espaço conforme).
2. Controle do Erro Não Conforme:
   • Para limitar o erro não conforme, os autores constroem uma reconstrução conforme em $H^2_0$ da solução discreta usando uma partição da unidade $C^1$.
   • Esta partição é construída utilizando espaços de elementos finitos conformes em $C^1$ (especificamente, a construção de [40] baseada em divisões de Alfeld da malha).
   • Crucialmente, os autores empregam decomposições de Helmholtz locais em patches de vértices (estrelas). Diferente de abordagens anteriores que dependiam de decomposições de Helmholtz globais (que introduzem constantes dependentes do número de buracos no domínio), esta abordagem local garante que a constante de estabilidade dependa apenas da regularidade de forma da malha.
3. Controle do Erro Conforme:
   • Duas cotas superiores distintas para o erro conforme são derivadas, cada uma utilizando um operador de interpolação específico:
     • Configuração 1: Utiliza o operador de interpolação de elementos finitos híbridos canônico ($C^{k+2}_T$). Esta abordagem baseia-se na propriedade de que a composição deste operador com a reconstrução HHO resulta em uma projeção $H^2$-elíptica. A cota resultante envolve apenas o termo de estabilização e a oscilação de dados, mas depende de uma suposição numérica referente ao erro de interpolação (Suposição 3.4).
     • Configuração 2: Utiliza o operador de interpolação $hp$ modificado de Babuška–Suri ($I^{k+2}_{mBS}$). Esta abordagem não requer suposições adicionais, mas introduz termos de indicadores de erro adicionais, especificamente resíduos de volume e saltos de fluxo normal/tangencial.
4. Estimadores de Resíduo: O artigo define indicadores de erro locais compreendendo termos de estabilização, resíduos de volume, saltos de fluxo normal e saltos tangenciais, escalonados por fatores $hp$ apropriados.

Principais Contribuições

• Primeiro estimador $hp$ para HHO em problemas biarmônicos: Este trabalho fornece o primeiro estimador de erro a posteriori baseado em resíduos para métodos HHO aplicados ao problema biarmônico em domínios gerais de 2D e 3D.
• Decomposição de Helmholtz Local com Robustez Topológica: O artigo estabelece uma nova estimativa de estabilidade para decomposições de Helmholtz locais em estrelas de vértices. Os autores provam que a constante de estabilidade depende exclusivamente da regularidade de forma da malha, evitando assim a degradação das constantes associadas ao número de buracos no domínio, um problema conhecido em abordagens de decomposição global.
• Limites Superiores Duais: A derivação de duas cotas superiores distintas para o erro conforme oferece flexibilidade. Uma cota é "limpa" (dependendo apenas de estabilização e oscilação), mas condicional a uma suposição verificada numericamente; a outra é incondicional, mas inclui termos de resíduo mais complexos.
• Extensão para Topologia Geral: A metodologia é explicitamente desenhada para lidar com domínios de topologia arbitrária (múltiplos buracos), uma característica frequentemente difícil de abordar na estimativa $hp$ para problemas de quarta ordem.

Resultos

• Limites Teóricos: O teorema principal (Teorema 5.6) estabelece que o erro total (norma de Hessiano quebrado mais estabilização) é limitado pela soma dos indicadores de erro propostos. As cotas mostram ser $h$-ótimas e $p$-subótimas (por no máximo uma ordem, ou potencialmente $1/2$ ordem sob a suposição verificada).
• Verificação Numérica:
  • Verificação da Suposição: Experimentos numéricos em 2D e 3D confirmam a validade da Suposição 3.4 em relação ao erro de interpolação do operador híbrido canônico.
  • Taxas de Convergência: Testes em soluções suaves em 3D demonstram que tanto o erro de energia quanto o estimador convergem na taxa ótima $O(\text{DoFs}^{-(k+1)/d})$.
  • Índices de Efetividade: Os estimadores produzem índices de efetividade moderados (variando aproximadamente entre 1,2 e 2,0) através de vários graus polinomiais e refinamentos de malha.
  • Robustez Topológica: Experimentos em domínios 2D com números crescentes de buracos (0, 1, 4, 8) mostram que os índices de efetividade permanecem estáveis, confirmando a afirmação teórica de que o estimador é robusto em relação à topologia do domínio.
  • Desempenho Adaptativo: Algoritmos $h$-adaptativos impulsionados pelo estimador recuperam com sucesso taxas de convergência ótimas para soluções singulares tanto em 2D (domínio em L) quanto em 3D (singularidade de ponto), superando o refinamento uniforme.

Significância
Os autores posicionam este trabalho como o primeiro a empregar uma partição da unidade $C^1$ combinada com decomposições de Helmholtz locais em estrelas de vértices para limitar erros não conformes no cenário $hp$ para problemas biarmônicos. A principal significância reside na capacidade de controlar o erro não conforme sem incorrer em constantes de estabilidade que se degradam com a complexidade topológica do domínio (número de buracos). Além disso, os autores observam que, embora a abordagem seja desenvolvida para métodos HHO, a técnica para limitar erros não conformes via decomposições de Helmholtz locais é potencialmente aplicável a outros métodos não conformes, como os métodos de Discontinuous Galerkin de Penalidade Interior (IPDG). O artigo não pretende fornecer limites inferiores de erro, observando que tais limites seguiriam de adaptações de argumentos existentes, mas provavelmente exibiriam uma subotimalidade $p$ pessimista devido às limitações das técnicas $hp$ atuais para a construção de funções bolha e operadores de extensão específicos.