An Extended Galerkin analysis in finite element exterior calculus

Resumo Técnico: Análise Galerkin Estendida no Cálculo Exterior de Elementos Finitos

1. Problema Abordado

O artigo foca na resolução numérica da equação de Hodge-Laplace no contexto do Cálculo Exterior de Elementos Finitos (FEEC - Finite Element Exterior Calculus). A equação modelo é dada por:
$$(d\delta + \delta d)u = f$$
onde $d$ é a derivada exterior e $\delta$ é a codiferencial.

O problema central é a falta de uma análise detalhada e unificada para Métodos de Galerkin Descontínuos (DG) aplicados a este problema no framework do FEEC. Embora métodos de Galerkin Contínuos (como o método misto de Arnold-Falk-Winther - AFW) sejam bem estabelecidos, a discretização DG da Laplaciana de Hodge mista envolve variáveis adicionais, fluxos numéricos e termos de estabilização complexos, especialmente quando se considera a continuidade inter-elemento tanto em relação a $d$ quanto a $\delta$.

2. Metodologia

Os autores propõem um Framework Galerkin Estendido (XG) unificado para o FEEC. A metodologia baseia-se nos seguintes pilares:

• Formulação de Sete Campos: Para a equação de Hodge-Laplace, o método introduz um sistema com sete variáveis:

  • Três campos principais aproximando $u$, $du$e$\delta u$ (denotados como $\sigma, \xi, u$).
  • Quatro variáveis de verificação ("check variables") definidas nas fronteiras dos elementos ($\check{u}, \check{u}^*, \check{\sigma}, \check{\xi}^*$).
  • Essas variáveis de fronteira são usadas para impor fracamente as condições de continuidade inter-elemento via um truque inspirado em Nitsche.

• Fluxos Numéricos e Penalização: O método utiliza fluxos numéricos que combinam médias e saltos (jumps) das soluções e seus traços. A penalização é aplicada às descontinuidades relacionadas tanto a $d$ quanto a $\delta$, permitindo que o método opere em subespaços totalmente descontínuos.

• Redução e Hibridização:

  • O sistema de sete campos pode ser reduzido a formulações com quatro ou três campos eliminando as variáveis de verificação.
  • Em um caso extremo, o método é hibridizado, resultando em um sistema condensado com apenas dois fluxos numéricos desconhecidos na fronteira, reduzindo drasticamente o número de graus de liberdade globais.

• Análise de Estabilidade (Inf-Sup):

  • Os autores derivam uma condição inf-sup discreta que é uniforme em relação a todos os parâmetros de discretização (tamanho da malha $h$) e parâmetros de penalização.
  • A análise utiliza a decomposição de Hodge discreta e desigualdades de Poincaré discretas para construir um sistema de ponto de sela perturbado, analisado via a teoria clássica de Babuška-Brezzi.

3. Principais Contribuições

1. Unificação de Métodos DG no FEEC: O framework XG proposto unifica várias classes de métodos DG existentes e gera novas famílias de métodos para a Laplaciana de Hodge.
2. Análise Unificada de Convergência: Fornece uma prova rigorosa de estabilidade e estimativas de erro quasi-ótimas para todos os parâmetros do método, algo que estava ausente na literatura para DG no FEEC.
3. Convergência para Métodos Conformes: Demonstra-se que, quando os parâmetros de penalização tendem a limites específicos, a solução do método XG converge para a solução do método misto conforme clássico (AFW). Isso valida o método XG como uma generalização robusta.
4. Hibridização Eficiente: Apresenta uma técnica de hibridização que reduz o sistema global a um problema definido apenas nas fronteiras dos elementos, facilitando a implementação e a resolução computacional.
5. Aplicação a Casos Específicos: O framework é aplicado e detalhado para casos específicos em $\mathbb{R}^3$ ($k=0, 1, 2, 3$), recuperando métodos conhecidos para a equação de Poisson (primal e mista) e gerando novos métodos DG para a Laplaciana Vetorial (problemas de Stokes e elasticidade).

4. Resultados Chave

• Estabilidade: Foi provado que o sistema reduzido satisfaz uma condição inf-sup uniforme, garantindo a existência e unicidade da solução discreta independentemente dos parâmetros de penalização (dentro de um intervalo adequado).
• Estimativas de Erro: Foram obtidas estimativas de erro a priori ótimas. A taxa de convergência depende da regularidade da solução e do grau dos polinômios utilizados, sendo $O(h^s)$ onde $s$ é limitado pela regularidade do domínio e da solução.
• Relação com AFW: O artigo estabelece matematicamente a conexão entre o método descontínuo proposto e o método misto conforme de Arnold-Falk-Winther, mostrando que o XG é um estágio intermediário que se torna conforme quando as penalidades são infinitas (ou em limites específicos).
• Novos Métodos para Laplaciana Vetorial: Para $k=1$ e $k=3$ em 3D, o método gera novas famílias de esquemas DG para a equação de Laplaciana vetorial, que são cruciais para problemas de elasticidade e eletromagnetismo.

5. Significado e Impacto

Este trabalho preenche uma lacuna significativa na teoria numérica de EDPs, conectando o poderoso framework do FEEC (que preserva propriedades topológicas e geométricas) com a flexibilidade dos métodos de Galerkin Descontínuos.

• Flexibilidade: Permite o uso de malhas não conformes e elementos de graus variados entre vizinhos, facilitando a adaptação local (h-adaptivity).
• Robustez: A análise uniforme em relação aos parâmetros de penalização torna o método robusto para diferentes configurações de malha e problemas físicos.
• Eficiência Computacional: A capacidade de hibridização reduz a complexidade do sistema global, tornando a resolução de problemas de alta dimensão ou em domínios complexos mais viável computacionalmente.
• Generalidade: Ao tratar a Laplaciana de Hodge de forma geral, o método cobre automaticamente a equação de Poisson, Maxwell e elasticidade, oferecendo uma abordagem unificada para uma vasta gama de problemas físicos.

Em suma, o artigo fornece a base teórica e algorítmica para uma nova geração de métodos numéricos que combinam a precisão geométrica do FEEC com a flexibilidade computacional dos métodos DG.