A Least-Squares Weak Galerkin Method for the Biharmonic Cauchy Problem

Resumo Técnico: Um Método de Galerkin Fraco de Mínimos Quadrados para o Problema de Cauchy Biharmônico

Enunciado do Problema
O artigo aborda o problema de Cauchy biarmônico, uma equação diferencial parcial de quarta ordem que surge em aplicações como flexão de placas, elasticidade e processamento de imagens. Diferente dos problemas de contorno padrão, o problema de Cauchy é definido em um domínio $\Omega$ onde as condições de contorno (valor da função e derivada normal) são prescritas apenas em uma porção acessível do contorno, $\Gamma_1$, enquanto nenhuma informação está disponível na porção inacessível, $\Gamma_2$.

Este problema é severamente mal posto no sentido de Hadamard. Embora a unicidade possa ser estabelecida sob suposições específicas via princípios de continuação única, a solução não depende continuamente dos dados de Cauchy. Consequentemente, pequenas perturbações nas medições de contorno podem levar a grandes erros na solução reconstruída. A instabilidade é amplificada pela natureza de quarta ordem do operador e pela necessidade de recuperar informações de contorno ausentes a partir de dados parciais.

Metodologia
Os autores propõem um método de elementos finitos de Galerkin Fraco de Mínimos Quadrados (LS-WG). A metodologia procede através dos seguintes passos:

1. Reformulação: A equação de quarta ordem $\Delta^2 u = f$ é reformulada como um sistema acoplado de duas equações de segunda ordem, introduzindo uma variável auxiliar $w = \Delta u$:
   \beginaligned \Delta w &= f, \\ \Delta u - w &= 0, \end{aligned}
   sujeito a dados de Cauchy em $\Gamma_1$.

2. Estrutura de Galerkin Fraco: O método utiliza espaços de elementos finitos fracos onde as funções são definidas como tripletos $v = \{v_0, v_b, v_n\}$, representando a aproximação do valor interior, do valor de contorno e da derivada normal, respectivamente. Esses componentes são tratados como incógnitas independentes e não são exigidos a serem contínuos através das fronteiras dos elementos ou a coincidirem com os traços dos valores interiores.

3. Laplaciano Fraco Discreto: O Laplaciano clássico é substituído por um operador de Laplaciano fraco discreto ($\Delta_{w,r,T}$), definido via uma formulação fraca que integra contra funções de teste polinomiais. Isso permite o uso de espaços de elementos finitos descontínuos.

4. Formulação de Mínimos Quadrados: O sistema acoplado é discretizado usando uma abordagem de mínimos quadrados. O método minimiza os resíduos das equações governantes e a consistência entre os componentes interior e de contorno das funções fracas. Isso resulta em um sistema linear simétrico definido positivo.

   • Estabilização: Um termo de estabilização é introduzido para impor fracamente a consistência entre os componentes interior ($v_0$) e de contorno ($v_b, v_n$).
   • Condições de Contorno: Os dados de Cauchy em $\Gamma_1$ são impostos via projeções $L^2$ sobre os componentes de contorno das funções fracas.

5. Flexibilidade de Malha: A formulação acomoda naturalmente malhas poligonais gerais, evitando a necessidade de elementos globalmente $C^1$-conformes (como os elementos de Argyris) que são difíceis de construir em malhas gerais.

Principais Contribuições
O artigo reivindica as seguintes contribuições específicas:

• Primeiro LS-WG para Cauchy Biharmônico: Apresenta o primeiro método de elementos finitos de Galerkin fraco de mínimos quadrados especificamente projetado para o problema de Cauchy biarmônico.
• Sistemas Simétricos Definidos Positivos: A formulação proposta gera um sistema linear simétrico definido positivo. Isso elimina a necessidade de condições inf-sup discretas tipicamente exigidas por métodos de elementos finitos mistos e evita estruturas de ponto de sela.
• Eliminação da Continuidade $C^1$: O método evita a construção de espaços de elementos finitos globalmente $C^1$-conformes, baseando-se em vez disso em espaços descontínuos com operadores fracos.
• Análise Teórica: Sob a suposição de que o problema contínuo admite uma solução única, os autores estabelecem a unicidade da solução LS-WG discreta. Além disso, derivam estimativas de erro de ordem ótima em uma norma de energia discreta.
• Solvers Eficientes: Devido à natureza simétrica definida positiva do sistema resultante, solvers iterativos eficientes (por exemplo, gradiente conjugado com precondicionadores multigrid) podem ser empregados.

Resultados Numéricos
Os autores validam o método através de experimentos numéricos no domínio do quadrado unitário utilizando três tipos de malhas: grades quadradas uniformes, grades triangulares e grades poligonais não convexas.

• Convergência: Os experimentos numéricos confirmam as taxas de convergência teóricas. Para soluções suaves, o método alcança ordens de convergência ótimas na norma de energia discreta e nas normas $L^2$ para vários graus polinomiais ($k=2, 3, 4$).
• Robustez: Os resultados demonstram que o método é eficaz em diferentes tipos de malha, mostrando pequenas diferenças de desempenho entre malhas quadradas, triangulares e poligonais não convexas.
• Desafios: O artigo observa que, para soluções com grandes variações perto da fronteira livre (a parte inacessível), a solução numérica torna-se mais difícil de aproximar, atingindo a precisão de máquina em níveis de grade mais finos devido aos grandes números de condição inerentes à natureza mal posta do problema.

Significância
O artigo posiciona o método LS-WG como uma alternativa robusta para problemas mal postos de quarta ordem. Ao combinar a estabilidade das formulações de mínimos quadrados com a flexibilidade da estrutura de Galerkin fraco, o método supera dois obstáculos principais na resolução de problemas de Cauchy biarmônicos: a dificuldade de construir elementos $C^1$-conformes e a complexidade de satisfazer condições inf-sup em formulações mistas. A capacidade de lidar com malhas poligonais gerais torna a abordagem particularmente adequada para geometrias complexas onde os métodos conformes tradicionais encontram dificuldades. O trabalho preenche uma lacuna na literatura, uma vez que métodos de Galerkin fraco existentes haviam sido desenvolvidos anteriormente quase exclusivamente para problemas de valor de contorno bem postos.