Equivalence of mixed and nonconforming methods on general polytopal partitions. Part I: Multiscale and projection methods

Resumo Técnico: Equivalência de Métodos Mistos e Não Conformes em Partições Polipodais Gerais

1. Problema e Contexto

O artigo aborda o problema de difusão com coeficiente variável em um domínio $\Omega \subset \mathbb{R}^d$ ($d=2,3$). O objetivo é estabelecer uma equivalência rigorosa entre duas formulações clássicas de elementos finitos:

• Formulação Primal: Baseada na energia potencial, buscando o potencial $u$.
• Formulação Mista: Baseada na energia complementar, buscando simultaneamente o fluxo $\sigma$ e o potencial $u$.

O foco central é a equivalência primal/mista em redes (malhas) gerais de polítopos (células que podem ser triângulos, quadriláteros, polígonos arbitrários, etc.), considerando métodos de multiescala e métodos de projeção. O trabalho investiga como métodos conformes (mistos) e não conformes (primais) podem ser reescritos um no outro, mesmo em contextos complexos de malhas polipodais e tensores de mobilidade variáveis.

2. Metodologia

O autor adota uma perspectiva de alto nível baseada em espaços virtuais e utiliza o conceito de hibridização (introduzido originalmente por Fraeijs de Veubeke e generalizado por Arnold e Brezzi). A metodologia segue os seguintes passos:

1. Hibridização Contínua:

   • O autor define duas formulações "hibridizadas" (uma primal e uma mista) que relaxam as condições de conformidade global ($H^1_0$ e $H(\text{div})$) para o interior de cada célula, introduzindo variáveis de contorno (traços) nas faces da malha.
   • Demonstra-se que as formulações primal e mista originais são equivalentes às suas versões hibridizadas, e que existe um caminho de equivalência entre as formulações hibridizadas primal e mista.

2. Construção Virtual de Dimensão Finita (Seção 3):

   • Define-se espaços de funções virtuais locais ($U_{\ell,k}(K)$ e $\Sigma_{\ell,k}(K)$) que preservam exatamente a relação constitutiva do modelo ($\sigma = b \nabla u$).
   • Estes espaços são definidos implicitamente através de soluções de EDPs locais (problemas de Dirichlet/Neumann), garantindo que a equivalência primal/mista seja mantida no nível discreto.
   • Estabelece-se que, para esses espaços virtuais, a equivalência entre os métodos primais e mistos é válida, preservando a conformidade $H(\text{div})$ para o fluxo e uma conformidade fraca (momentos de face) para o potencial.

3. Métodos de Projeção (Seção 4):

   • Reconhece-se que os espaços virtuais definidos acima são frequentemente não computáveis na prática (contendo funções não polinomiais).
   • Para contornar isso, o autor introduz Métodos de Projeção. A ideia é substituir o espaço virtual computável por um espaço polinomial (ou analiticamente conhecido) e projetar o gradiente do potencial sobre o espaço de fluxo misto.
   • Deriva-se um critério simples de bem-postura e equivalência: a sobrejetividade do mapa de graus de liberdade (DoF) do espaço de fluxo. Se este mapa for sobrejetivo, o método misto é equivalente a um método primal de projeção.

3. Principais Contribuições

• Generalização da Equivalência Multiescala: Estabelece a primeira equivalência de "primeiro nível" entre quatro abordagens diferentes de métodos de multiescala (livres de oversampling) em malhas polipodais gerais. Isso amplia resultados anteriores (como [10]) para o contexto de métodos HHO (Hybrid High-Order) e MHM (Multiscale Hybrid-Mixed).
• Critério Prático de Bem-Postura: Fornece um critério verificável na prática para garantir a bem-postura (uniforme) e a equivalência primal/mista em malhas polipodais gerais. Este critério é uma reformulação simplificada de resultados teóricos anteriores, aplicável a métodos híbridos auto-estabilizados.
• Reinterpretação de Métodos Híbridos: Demonstra que métodos híbridos de alta ordem auto-estabilizados (como HHO e Weak Galerkin) podem ser reinterpretados como esquemas de projeção. Isso fornece equivalentes de elementos finitos mistos para esses métodos, algo que não era trivialmente conhecido.
• Validade em Malhas Arbitrárias:
  • Prova que a equivalência entre o método de Raviart-Thomas de ordem mais baixa (mistos) e o método de Crouzeix-Raviart enriquecido com "bolhas" (primal) é válida em simples arbitrários (sem restrições geométricas de regularidade), corrigindo uma lacuna na literatura que assumia malhas regulares.
  • Estende resultados análogos para hiper-retângulos e o espaço de Rannacher-Turek.
• Distinção entre Espaços Computáveis e Virtuais: Clarifica a distinção entre métodos onde a relação constitutiva é preservada exatamente (espaços virtuais) e métodos onde ela é preservada via projeção (espaços polinomiais), mostrando que ambos levam a formulações equivalentes sob certas condições.

4. Resultados Chave

• Diagrama de Equivalência: O artigo valida um diagrama de equivalência (Figura 1) que conecta:
  • Primal $\leftrightarrow$ Hibridizado Primal
  • Misto $\leftrightarrow$ Hibridizado Misto
  • Primal $\leftrightarrow$ Misto (via os dois caminhos intermediários).
• Método de Raviart-Thomas (RT) vs. Crouzeix-Raviart (CR):
  • No caso de tensores constantes por célula e malhas simpliciais, o método misto RT de ordem mais baixa é equivalente a uma versão enriquecida do método primal CR (com uma modificação no lado direito do sistema).
  • A solução do sistema misto pode ser obtida resolvendo um sistema simétrico positivo definido (SPD) condensado (primal) e reconstruindo o fluxo localmente.
• Espaços Virtuais vs. Projeção:
  • Para espaços de fluxo polinomiais de ordem superior (como BDM ou RT de ordem alta), a relação constitutiva exata $\sigma = b \nabla u$ não pode ser mantida dentro do espaço polinomial.
  • Nesses casos, a equivalência é restaurada através de um método de projeção, onde o fluxo é a projeção ortogonal do gradiente do potencial sobre o espaço de fluxo misto.

5. Significado e Impacto

Este trabalho é fundamental para a teoria de métodos numéricos em malhas polipodais (como VEM, HHO, DG, etc.) por várias razões:

1. Unificação Teórica: Unifica a compreensão de métodos de multiescala e métodos de projeção, mostrando que muitos métodos modernos são, na essência, diferentes faces da mesma moeda (equivalência primal/mista).
2. Guia para Implementação: O critério de bem-postura derivado oferece uma ferramenta prática para desenvolvedores de códigos verificarem se suas formulações híbridas são estáveis e equivalentes a métodos mistos bem compreendidos.
3. Flexibilidade Geométrica: Ao remover restrições geométricas (como a necessidade de malhas regulares para a equivalência RT/CR), o trabalho valida o uso de métodos de alta ordem em malhas complexas e distorcidas, comuns em aplicações de engenharia e física.
4. Base para Parte II: O artigo serve como a base teórica para a segunda parte da série, que abordará a implementação prática em métodos de elementos polipodais gerais, onde a conformidade $H(\text{div})$ global pode ser relaxada.

Em suma, o artigo fornece a fundação matemática rigorosa que justifica o uso de esquemas de projeção e hibridização em malhas gerais, garantindo que as propriedades de conservação e precisão dos métodos mistos sejam herdadas por métodos primais não conformes modernos.