Towards $C^0$ finite element methods for fourth-order elliptic equation. Part I: general boundary conditions

Resumo Técnico: Rumo aos Métodos de Elementos Finitos $C^0$ para Equações Elípticas de Quarta Ordem, Parte I: Condições de Contorno Gerais

Enunciado do Problema
Este artigo aborda a aproximação numérica de equações elípticas de quarta ordem, especificamente a equação biarmônica $\Delta^2 u = f$, em domínios poligonais $\Omega \subset \mathbb{R}^2$ com condições de contorno gerais. A fronteira $\partial \Omega$ é decomposta em segmentos disjuntos de Dirichlet ($\Gamma_D$) e Neumann ($\Gamma_N$). O desafio principal reside no desenvolvimento de métodos eficazes de elementos finitos $C^0$ para este problema. Embora as abordagens padrão para a equação biarmônica com condições de contorno de Navier ($\Gamma_D = \partial \Omega$) frequentemente dependam do desacoplamento do problema em um sistema de equações de Poisson, esta formulação mista "ingênua" falha em convergir para a solução verdadeira para domínios não convexos ou configurações específicas de condições de contorno mistas. Essa falha é atribuída ao paradoxo de Sapongyan, onde o espaço de solução da formulação mista desacoplada não coincide com o espaço de solução do problema biarmônico original devido a termos de baixa regularidade originados de cantos reentrantes.

Metodologia
Os autores propõem uma formulação mista modificada que decompõe o problema biarmônico em um sistema de equações de Poisson, aumentada por problemas de Poisson intermediários adicionais para garantir que a solução permaneça dentro do espaço de Sobolev correto ($H^2(\Omega)$).

1. Análise do Espaço de Solução: O artigo analisa o alcance do operador Laplace $\Delta$ mapeando $H^2(\Omega)$ (com as devidas condições de contorno) em $L^2(\Omega)$. Identifica-se que, para domínios poligonais com grandes ângulos internos, o alcance $S = \Delta(V^2)$ é um subespaço fechado próprio de $L^2(\Omega)$ com um complemento ortogonal de dimensão finita $S^\perp$.
2. Construção de Base: Os autores constroem uma base para $S^\perp$ consistindo de funções $\xi_m$ definidas próximas ao maior ângulo interno $\omega$. Estas funções são compostas por termos singulares $s_m(r, \theta)$ (derivados dos autovalores da parte angular do Laplaciano) e termos de correção regulares $\zeta_m$ obtidos através da resolução de problemas auxiliares de Poisson. A dimensão de $S^\perp$ ($d_\perp$) depende da magnitude do maior ângulo interno e da combinação específica de condições de contorno (Dirichlet/Neumann) nas arestas adjacentes.
3. Algoritmo Modificado: O algoritmo de elementos finitos $C^0$ proposto (Algoritmo 3.1) procede da seguinte forma:
   • Resolve um problema de Poisson primário para uma função auxiliar $w$.
   • Resolve problemas de Poisson auxiliares para computar as funções de base $\xi_{m,h}$ para o complemento ortogonal.
   • Projeta $w$ sobre o subespaço $S$ removendo sua componente em $S^\perp$ (computando coeficientes via um sistema de matriz de Gram).
   • Resolve um segundo problema de Poisson usando o termo de fonte projetado $w_S$ para obter a aproximação final $u_h$.
4. Extensão para Condições de Neumann: A metodologia é estendida para o caso de Neumann puro ($\Gamma_N = \partial \Omega$), onde a condição de compatibilidade $\int_\Omega f = 0$ é necessária, e a dimensão do complemento ortogonal é mostrada como sendo 1 para cantos reentrantes.

Principais Contribuições

• Resolução do Paradoxo de Sapongyan: O artigo fornece um arcabouço teórico rigoroso demonstrando que o desacoplamento ingênuo da equação biarmônica falha para domínios poligonais gerais e condições de contorno mistas. Propõe uma formulação corrigida que garante a convergência para a solução verdadeira.
• Decomposição Generalizada: Diferente de trabalhos anteriores restritos a condições de Navier, este método lida com combinações arbitrárias de condições de Dirichlet/Neumann. O número de equações de Poisson adicionais necessárias (0, 1 ou 2) é explicitamente determinado pelo maior ângulo interno e pelos tipos locais de condições de contorno.
• Análise de Erro: Os autores estabelecem estimativas de erro rigorosas para o método de elementos finitos lineares $C^0$ proposto. Provam que a taxa de convergência da solução $u_h$ na norma $H^1$ é $O(h^{\min\{1, 2\alpha\}})$, onde $\alpha$ depende do ângulo interno e das condições de contorno. Especificamente, para condições de contorno mistas com ângulos $\omega > \pi$, a taxa é limitada pela singularidade, enquanto para os casos de Navier/Neumann, taxas ótimas são frequentemente recuperadas.
• Verificação Numérica: Experimentos numéricos extensivos em vários domínios poligonais (incluindo geometrias de canto em L e cantos reentrantes) validam as previsões teóricas. Os resultados confirmam que o método ingênuo converge para soluções incorretas em casos onde o complemento ortogonal é não trivial, enquanto o método modificado proposto converge para a solução de referência (calculada via $C^0$-IPDG).

Resultos e Significância
O artigo demonstra que a formulação mista modificada proposta é bem definida e eficaz para resolver a equação biarmônica em domínios poligonais gerais. Experimentos numéricos mostram que o método atinge as taxas de convergência previstas e resolve com sucesso as discrepâncias observadas em formulações mistas ingênuas.

A significância deste trabalho reside na sua capacidade de estender a eficiência computacional dos métodos de elementos finitos $C^0$ (que utilizam polinômios contínuos por partes padrão) para problemas de quarta ordem com condições de contorno gerais, uma tarefa anteriormente dificultada pela complexidade dos elementos $C^1$ ou pelas falhas de convergência de estratégias de desacoplamento simples. Ao caracterizar a influência específica das condições de contorno e da geometria do domínio sobre o espaço de solução, os autores fornecem um arcabouço algorítmico robusto que garante a correção matemática sem sacrificar a simplicidade de implementação. O artigo serve como a Parte I de uma série, estabelecendo a base teórica para lidar com equações de quarta ordem mais complexas com termos de ordem inferior em trabalhos subsequentes.