An Approximate Bounded Cochain Projection

Resumo Técnico: Uma Projeção de Co-cadeia Limitada Aproximada

Enunciado do Problema
O artigo aborda a construção de um projetor, $\Pi^{(k)}$, que mapeia um complexo de Hilbert de formas diferenciais de dimensão infinita, $H\Lambda^{(k)}(\Omega)$, para um subcomplexo de polinômios por partes de dimensão finita, $H\Lambda^{(k)}_h(\Omega)$. O objetivo é satisfazer três propriedades específicas características de uma "projeção de co-cadeia limitada":

1. Idempotência: $\Pi^{(k)}\Pi^{(k)} = \Pi^{(k)}$.
2. Limitação Uniforme: O projetor é limitado na norma de Sobolev ($H\Lambda^{(k)}$) com uma constante independente do parâmetro de malha $h$.
3. Propriedade de Comutação: O projetor comuta com a derivada exterior ($d$), tal que $d\Pi^{(k)} = \Pi^{(k+1)}d$.

Embora as projeções de co-cadeia limitadas sejam bem estabelecidas para domínios contráteis (onde o complexo de de Rham é exato), o artigo destaca uma dificuldade em domínios não contráteis. Nesses domínios, o complexo possui cohomologia não trivial (formas harmônicas). Construções padrão frequentemente falham em manter a limitação uniforme ou a comutação exata quando o subespaço de dimensão finita não captura perfeitamente a topologia do domínio (ou seja, quando o hiato entre os espaços harmônicos de dimensão infinita e finita é não nulo).

Metodologia
Os autores propõem uma construção baseada em uma formulação mista aumentada e um complexo discreto de dois níveis.

1. Formulação Mista Aumentada:
   Em vez de depender de operadores de suavização seguidos de interpolação (que podem carecer de limitação uniforme na norma $H\Lambda$), os autores definem o projetor $\Pi^{(k)}$ via um problema variacional. Dado $\phi^{(k)} \in H\Lambda^{(k)}(\Omega)$, o projetor busca $\lambda^{(i)} \in H\Lambda^{(i)}_h(\Omega)$ e uma variável harmônica $h^{(k+1)} \in H^{(k+1)}_h$ satisfazendo um sistema de equações simétrico.

   • Este sistema introduz um multiplicador de Lagrange, $h^{(k+1)}$, que impõe a restrição de que a derivada projetada pertença ao espaço das formas exatas.
   • Os autores provam a boa colocação deste sistema utilizando condições de estabilidade, incluindo a coercividade da forma bilinear em formas fechadas e condições inf-sup para a derivada exterior e o espaço harmônico.

2. Propriedades do Sistema Aumentado:

   • Idempotência: Provada ao mostrar que, se a entrada já estiver no espaço de dimensão finita, a solução do sistema misto é a própria entrada.
   • Limitação: Provada através da decomposição da solução em componentes fechadas e ortogonais, utilizando a desigualdade de Poincaré para limitar as normas $L^2$ e de gráfico.
   • Propriedade de Comutação: O artigo demonstra que o projetor comuta com $d$ se o domínio for contrátil ou se o espaço de dimensão finita for "topologicamente conforme" (ou seja, o hiato entre os espaços harmônicos discretos e contínuos é zero). Se o espaço não for topologicamente conforme, o multiplicador harmônico $h^{(k+1)}$ é não nulo, representando a inconsistência topológica, e a propriedade de comutação é satisfeita apenas até a magnitude deste erro.

3. Projeção de Co-cadeia Limitada Aproximada (Construção de Dois Níveis):
   Para abordar domínios não contráteis onde o espaço discreto pode não ser topologicamente conforme, os autores propõem a composição de dois operadores:

   • Passo 1: Um projetor $\Pi^{(k)}$ mapeando o espaço de dimensão infinita para um espaço de dimensão finita enriquecido $\tilde{H}\Lambda^{(k)}_h(\Omega)$. Este espaço enriquecido é construído para ser topologicamente conforme (hiato = 0).
   • Passo 2: Um projetor canônico $\Pi^{(k)}_C$ mapeando o espaço enriquecido $\tilde{H}\Lambda^{(k)}_h(\Omega)$ de volta para o espaço computacional alvo $H\Lambda^{(k)}_h(\Omega)$.
     A composição $\Pi^{(k)}_C \Pi^{(k)}$ serve como o projetor de co-cadeia limitada aproximado final.

Principais Contribuições e Resultados

• Construção de um Projetor Limitado: O artigo fornece uma construção rigorosa de um projetor que é idempotente e uniformemente limitado na norma de Sobolev para domínios contráteis e não contráteis.
• Tratamento de Domínios Não Contráteis: Os autores mostram que, em domínios não contráteis, a propriedade de comutação é satisfeita exatamente se o espaço harmônico discreto for topologicamente conforme. Caso contrário, o erro na comutação é quantificado pelo multiplicador harmônico $h^{(k+1)}$, que mede a discrepância topológica.
• Estabilidade Discreta: Diferente de construções anteriores que dependem de regularizações contínuas, esta abordagem utiliza uma formulação mista que evita dificuldades analíticas associadas a inversos à direita contínuos da derivada exterior.
• Comutação Aproximada: O artigo demonstra que, ao utilizar uma abordagem de dois níveis (enriquecendo o espaço para resolver a topologia antes de projetar para a malha final), pode-se obter um projetor que é limitado e comuta com a derivada exterior até uma precisão arbitrária.

Significância e Alegações
O artigo afirma apresentar uma construção de um projetor que cumpre as três propriedades definidoras de uma projeção de co-cadeia limitada. Especificamente, ele afirma que:

• Em domínios contráteis, o projetor proposto é idempotente, uniformemente limitado e comuta exatamente com a derivada exterior.
• Em domínios não contráteis, o projetor permanece idempotente e uniformemente limitado. A propriedade de comutação é satisfeita "até uma precisão arbitrária", o que significa que o erro pode ser tornado arbitrariamente pequeno ao refinar a resolução topológica do espaço intermediário (por exemplo, via refinamento-$p$ até que o multiplicador harmônico desapareça até a precisão de máquina).

O trabalho baseia-se no framework de Cálculo Exterior de Elementos Finitos (FEEC), refinando ideias existentes sobre construções locais e condições de contorno. Ele se distingue ao oferecer uma construção que é "totalmente local", no sentido de evitar inversos à direita contínuos, baseando-se, em vez disso, em problemas discretos postos em um refinamento auxiliar para lidar com restrições topológicas. Os autores posicionam isso como um método para garantir estabilidade e consistência em métodos numéricos para formas diferenciais em geometrias complexas.