Stabilization-Free H(curl) and H(div)-Conforming Virtual Element Method

Resumo Técnico: VEM de Conformidade $H(\text{curl})$ e $H(\text{div})$ Sem Estabilização

Enunciado do Problema
O Método de Elementos Virtuais (VEM) oferece vantagens significativas sobre o Método de Elementos Finitos (FEM) ao lidar com politopos estrelados gerais, sendo ideal para geometrias complexas e refinamento de malha adaptativo. No entanto, implementações padrão de VEM para espaços $H(\text{curl})$ e $H(\text{div})$ requerem termos de estabilização para garantir a coercividade. Essa necessidade surge porque o espaço de elementos virtuais $V(K)$ é maior que o espaço polinomial $P_k(K)$, fazendo com que a projeção $L^2$ padrão seja instável (ou seja, $\|\Pi_k v\|_K \not\gtrsim \|v\|_K$). Esses termos de estabilização introduzem vários inconvenientes: impactam negativamente a estabilidade numérica e os números de condição, exigem design específico para o problema (particularmente para problemas não lineares) e exibem comportamento isotrópico, que tem dificuldade com coeficientes ou malhas anisotrópicas. Abordagens existentes sem estabilização exigem graus de liberdade (DoFs) internos excessivos, o que aumenta os custos computacionais, ou modificam o espaço de projeção de formas que diluem a flexibilidade geométrica do VEM ou falham em preservar propriedades essenciais como a comutatividade com operadores diferenciais.

Metodologia
Os autores propõem um novo VEM sem estabilização para espaços de conformidade $H(\text{curl})$ e $H(\text{div})$ de ordem geral em $\mathbb{R}^3$. O núcleo da metodologia envolve a construção de novos projetores de serendipidade e espaços de serendipidade correspondentes com um número minimizado de DoFs.

1. Projetores de Serendipidade: Diferente de trabalhos anteriores que dependem de projetores triviais ou DoFs internos, os autores definem projetores específicos ($\Pi^S_k$) para espaços de faces, arestas e nós. Esses projetores mapeiam funções virtuais para espaços polinomiais usando apenas DoFs de contorno "preservados" e um conjunto específico de momentos internos. Os projetores são definidos via condições de ortogonalidade envolvendo traços normais/tangenciais e momentos internos (por exemplo, $(v - \Pi^S_k v, x_K \times p_{kr})_K = 0$).
2. Espaços de Serendipidade: Os autores definem espaços de serendipidade ($V^S$) restringindo os espaços de elementos virtuais padrão de modo que os DoFs "reduzidos" sejam determinados pelos DoFs preservados via o projetor de serendipidade. Isso garante que a projeção $L^2$ sobre um espaço polinomial de ordem suficientemente alta $P_l$ seja estável e computável usando apenas os DoFs de contorno.
3. Complexo de De Rham: O método é projetado para lidar com a cadeia completa do complexo de De Rham em $\mathbb{R}^3$ ($H^1 \to H(\text{curl}) \to H(\text{div}) \to L^2$). A construção garante comutatividade exata com operadores diferenciais ($\nabla, \nabla\times, \nabla\cdot$), preservando a estrutura topológica do problema.
4. Análise de Estabilidade: Os autores provam que, para uma ordem de projeção $l$ suficientemente grande, a equivalência de norma $\|\Pi_l v\|_K \approx \|v\|_K$ se mantém para funções no espaço de serendipidade. Isso elimina a necessidade de termos de estabilização enquanto mantém propriedades de aproximação ótimas.

Principais Contribuições

• Tratamento Completo de De Rham: O trabalho fornece um tratamento unificado para todo o complexo de De Rham em 3D, garantindo integração contínua entre diferentes operadores diferenciais, uma característica ausente em trabalhos anteriores sem estabilização limitados a $H^1$.
• Preservação de Propriedades Matemáticas: Os espaços de serendipidade propostos mantêm rigorosamente a continuidade de contorno, conformidade de operador, comutatividade exata com operadores diferenciais e B-compatibilidade.
• Eficiência Computacional: Ao minimizar o número de DoFs e utilizar projeções de serendipidade computadas exclusivamente a partir de dados de contorno, o método reduz significamente o overhead computacional. Isso é particularmente benéfico para problemas de autovalores onde os DoFs internos não podem ser eliminados via condensação estática.
• Análise de Erro Rigorosa: O artigo fornece uma prova teórica de propriedades de aproximação ótimas e estabilidade para as projeções polinomiais de alta ordem.

Resultos
As alegações teóricas são validadas através de experimentos numéricos em problemas de autovalores de Maxwell utilizando malhas poliedrais de Voronoi.

• Taxas de Convergência: Os resultados numéricos confirmam as taxas de convergência teóricas. Para soluções infinitamente regulares, o método atinge convergência de $O(h^{m+1})$ na norma $L^2$ e de $O(h^{2m+2})$ para autovalores. Para soluções singulares (onde a regularidade é limitada por $s$), as taxas são de $O(h^s)$ e $O(h^{2s})$, respectivamente.
• Custo Computacional: Mostra-se que o tempo de montagem do VEM sem Estabilização (SFVEM) proposto é apenas moderadamente maior que o FEM convencional, verificando a eficiência do método. A ordem de projeção $l$ necessária é encontrada como próxima da ordem polinomial $k$ na prática.
• Comparação: O método elimina com sucesso as complicações relacionadas à estabilização enquanto iguala a precisão do FEM convencional e do VEM estabilizado.

Significância e Alegações
Os autores alegam que este trabalho oferece uma abordagem fundamentalmente diferente para o VEM sem estabilização que supera as limitações de estratégias existentes (DoFs excessivos ou perda de flexibilidade geométrica). Ao eliminar termos de estabilização, o método melhora a estabilidade numérica e os números de condição, tornando-o particularmente valioso para:

• Formulações mistas.
• Problemas não lineares.
• Aplicações que exigem comportamento numérico robusto através de diversas configurações de malha.

O artigo conclui que, embora a técnica possa ser teoricamente generalizada para dimensões superiores, as complexidades de implementação prática e o crescimento rápido dos DoFs de contorno atualmente limitam a discussão ao cenário 3D. O método proposto representa uma alternativa teoricamente rigorosa e computacionalmente eficiente para resolver equações diferenciais parciais, avançando a aplicabilidade do VEM na matemática computacional e engenharia.