$\boldsymbol{H}(\textbf{curl})$-reconstruction of piecewise polynomial fields with application to $hp$-a posteriori nonconforming error analysis for Maxwell's equations

Título: Reconstrução H(curl) de campos polinomiais por partes com aplicação à análise de erro a posteriori não conformes hp para as equações de Maxwell

1. Problema e Motivação

O artigo aborda um desafio fundamental na análise numérica de equações de Maxwell: a estimativa de erro a posteriori para métodos de elementos finitos não conformes, especificamente no contexto de aproximações hp (refinamento de malha $h$ e aumento do grau polinomial $p$).

• Contexto: Métodos de Galerkin Descontínuo (dG) são amplamente utilizados para discretizar as equações de Maxwell devido à sua flexibilidade em lidar com descontinuidades de materiais e geometrias complexas. No entanto, como as soluções aproximadas ($E_h$) pertencem a espaços não conformes (não pertencem a $H(\text{curl})$), a análise de erro tradicional é difícil.
• O Desafio: Para realizar uma análise de erro rigorosa, é necessário comparar a solução discreta não conforme com uma solução "conforme" (que pertença a $H_0(\text{curl})$). Métodos anteriores frequentemente usavam operadores de média que resultavam em estimativas de erro subótimas em relação ao grau polinomial $p$.
• Objetivo: Desenvolver um operador de reconstrução $H(\text{curl})$-conforme para campos polinomiais por partes arbitrários, garantindo estimativas de erro ótimas em $h$ e subótimas apenas em $p$ (com perda de apenas $1/2$ ordem), e aplicá-lo à análise de erro a posteriori para o método dG simétrico de penalidade interior.

2. Metodologia

Os autores propõem uma construção explícita de um campo reconstruído $R_c(E_h) \in H_0(\text{curl}; D)$ a partir de um campo polinomial por partes $E_h$.

• Estratégia de Construção Local:

  • A reconstrução é realizada localmente sobre os "patches" de vértice (aglomerados de células que compartilham um vértice da malha).
  • Utiliza-se uma partição da unidade induzida por funções de base "hat" ($\psi_a$) para combinar as soluções locais.
  • Para cada vértice $a$, resolve-se dois problemas locais acoplados no patch $D_a$:
    1. Um problema de minimização para a parte solenoidal (divergência livre) usando decomposição de Helmholtz local.
    2. Um problema de Poisson para a parte irrotacional (gradiente de um escalar).
  • A solução global é a soma das contribuições locais: $R_c(E_h) = \sum_a (U^a + \nabla \theta^a)$.

• Ferramentas Analíticas Chave:

  • Decomposição de Helmholtz Local: Usada para separar o campo em componentes rotacionais e gradientes dentro de cada patch.
  • Desigualdade de Poincaré Quebrada (Broken Poincaré): Os autores derivam uma nova desigualdade de Poincaré que preserva o divergente e lida com saltos (jumps) no curl, essencial para controlar o erro $L^2$.
  • Regularidade Elíptica: Sob uma suposição razoável sobre a regularidade elíptica uniforme para problemas de Poisson com condições de Neumann nos patches, é possível melhorar a estimativa de erro $L^2$.

3. Principais Contribuições e Resultados

A. O Operador de Reconstrução e Estimativas de Erro

O principal teorema (Teorema 3.7) estabelece que a diferença entre o campo original $E_h$ e o campo reconstruído $R_c(E_h)$ é limitada pelos saltos (jumps) do campo original e do seu curl nas interfaces da malha.

• Norma de Curl Quebrada ($\|\nabla \times (\cdot)\|$):
  • A estimativa é ótima em $h$ e subótima em $p$ por $1/2$ ordem.
  • O erro é controlado pelos saltos tangenciais de $E_h$ e pelos saltos normais de $\nabla \times E_h$.
• Norma $L^2$:
  • Sem suposições adicionais: ótima em $h$ e subótima em $p$ por $3/2$ ordem.
  • Com a suposição de regularidade elíptica uniforme (Suposição 3.5): A suboptimalidade em $p$ é reduzida para $1/2$ ordem.
• Vantagem Topológica: Ao construir localmente nos patches de vértice (que têm topologia trivial), o método evita dependências complexas da topologia global do domínio, ao contrário de reconstruções globais que exigiriam condições de contorno específicas ou dependência de cohomologia.

B. Extensões

• Condições de Contorno Mistas: O método é adaptável para condições de Dirichlet e Neumann mistas.
• Malhas Polipodais Aglomeradas: A análise estende-se a malhas poligonais/polidrais aglomeradas, desde que existam sub-malhas simpliciais subjacentes com tamanho local comparável.
• Domínios Convexos: Em domínios convexos, uma reconstrução global (sem patches) pode atingir estimativa $L^2$ ótima em $p$, mas a constante depende da topologia global.

C. Aplicação: Análise de Erro hp para Maxwell

O operador é aplicado ao método de Galerkin Descontínuo (dG) com penalidade interior simétrica para equações de Maxwell simplificadas (operador curl-curl + termo de ordem zero).

• Resultado Principal (Teorema 6.8): É estabelecida uma estimativa superior global (reliability) para o erro a posteriori.
  • O erro total é limitado pelo indicador de erro computável $\eta$, que inclui resíduos de volume, saltos de fluxos e termos de não-conformidade.
  • Esta é a primeira análise hp completa para esse tipo de problema não conforme baseada em decomposições de Helmholtz locais.
• Eficiência Local: Também são provadas estimativas de eficiência local (o indicador não superestima o erro localmente), garantindo que o estimador seja confiável para guiar a adaptação da malha (h-adaptivity) e do grau polinomial (p-adaptivity).

4. Significado e Impacto

1. Preenchimento de Lacuna Teórica: O artigo preenche uma lacuna na literatura sobre reconstrução $H(\text{curl})$ para campos polinomiais arbitrários, garantindo propriedades de aproximação hp que métodos anteriores (baseados apenas em graus de liberdade médios) não conseguiam atingir.
2. Independência Topológica: A abordagem baseada em patches de vértice permite que as estimativas de erro sejam locais e independentes da topologia complexa do domínio global, o que é crucial para simulações em geometrias 3D complexas.
3. Aplicabilidade Prática: A análise fornece uma base teórica sólida para o uso de indicadores de erro a posteriori em códigos de elementos finitos para eletromagnetismo, permitindo a automação de refinamento de malha e grau polinomial para obter soluções precisas de forma eficiente.
4. Novas Desigualdades: A derivada de uma desigualdade de Poincaré "quebrada" que preserva o divergente é um resultado auxiliar de interesse independente para a análise de métodos numéricos em espaços de Sobolev.

Em resumo, o trabalho fornece as ferramentas matemáticas necessárias para analisar rigorosamente e controlar erros em simulações de alta ordem de equações de Maxwell usando métodos não conformes, com implicações diretas para o desenvolvimento de algoritmos adaptativos robustos em eletromagnetismo computacional.