The Tempered Finite Element Method

Resumo Técnico: O Método de Elementos Finitos Temperado (TFEM)

Enunciado do Problema
O Método de Elementos Finitos (FEM) é a técnica numérica dominante na física computacional e na engenharia, renomada por sua capacidade de lidar com geometrias complexas via malhas não estruturadas. No entanto, formulações padrão de FEM exibem fragilidade ao encontrar configurações de malha específicas, particularmente aquelas contendo "capas" (triângulos com um ângulo aproximando-se de 180°) dispostas em bandas. Enquanto elementos degenerados isolados (agulhas ou capas isoladas) não impedem a convergência, bandas de capas induzem um fenômeno de "travamento" (locking). Em tais configurações, a solução numérica é forçada a ser linear ao longo da banda, independentemente da solução real, fazendo com que o tamanho efetivo da malha seja determinado pelo comprimento da banda, e não pela sua espessura. Isso resulta em uma perda de convergência para o FEM padrão. Historicamente, a "condição do ângulo máximo" (exigindo que todos os ângulos estivessem afastados de 180°) era considerada essencial para a convergência, embora trabalhos recentes tenham mostrado que ela não é estritamente necessária para elementos isolados. O desafio permanece em lidar com malhas contendo elementos de medida zero ou quase degenerados, que são relevantes para métodos de rastreamento de interface (ex: X-Mesh), onde as interfaces são rastreadas por retransmissão de malha, podendo criar bandas de capas.

Metodologia
Os autores propõem o Método de Elementos Finitos Temperado (TFEM), uma modificação do FEM clássico que permite a convergência em malhas contendo bandas de capas. A ideia central é inspirada na têmpera metalúrgica: assim como a têmpera alivia tensões internas no aço endurecido para aumentar a flexibilidade, o TFEM "suaviza" localmente as propriedades do material na banda degenerada para permitir que a solução numérica se conforme ao problema contínuo.

Matematicamente, para problemas elípticos como a equação de Poisson, a matriz de rigidez padrão $K_{ij}$ envolve o inverso do determinante Jacobiano $J$ do mapeamento do elemento. Em elementos degenerados, $J \to 0$, levando à divisão por zero ou eigenvalues infinitos. O TFEM modifica a matriz de rigidez local limitando o Jacobiano para longe de zero usando um parâmetro $J_{min}$:
$$\tilde{K}_{ij} = \int_{\hat{E}} \frac{1}{\max(J, J_{min})} [\dots] d\xi$$
Isso é equivalente a introduzir um coeficiente de difusão variável $d(x)$ na formulação fraca:
$$d(x) = \begin{cases} 1 & x \in \Omega \setminus B \\ D = \frac{J}{J_{min}} & x \in B \end{cases}$$
onde $B$ é a banda de elementos degenerados. O parâmetro $J_{min}$ é escolhido para depender do tamanho da malha $h$ (especificamente $J_{min} \sim h^{d+1}$, onde $d$ é a dimensão espacial) para equilibrar a redução do coeficiente de difusão (para evitar o travamento) com a preservação das propriedades de aproximação.

Para o caso de elementos de medida exatamente zero ($J=0$), os autores demonstram que o limite do esquema TFEM corresponde a um método de mortaring (método de argamassa) com penalização do salto da solução através da interface.

Principais Contribuições

1. Prova Teórica de Convergência: O artigo fornece uma análise de erro rigorosa para o TFEM em bandas de capas. Prova-se que, com uma escolha ótima do parâmetro de têmpera $D \sim \bar{h}/h^2$ (onde $\bar{h}$ é a largura da banda), o método alcança convergência ótima de $O(h)$ na semi-norma $H^1$, igualando a taxa do FEM padrão em malhas regulares.
2. Equivalência ao Mortaring: Para bandas de largura exatamente zero, os autores provam que o TFEM é equivalente a um método de mortar não conforme com penalização de salto, fornecendo uma ponte teórica entre o tratamento de malhas degeneradas e técnicas de acoplamento de interface.
3. Extensão para Físicas Complexas: O método é estendido com sucesso além do Laplaciano escalar para:
   • Elasticidade Linear: Demonstrando precisão comparável ao FEM padrão para análise de tensão através de interfaces.
   • Mortaring de Malhas Não Conformes: Mostrando que o TFEM pode conectar perfeitamente malhas de diferentes resoluções sem elementos de transição.
   • Elementos de Alta Ordem: Validando que o método preserva taxas de convergência ótimas para ordens polinomiais de até 4, desde que $J_{min}$ seja escalonado apropriadamente com a ordem do elemento.
   • Advecção: Aplicando o TFEM a problemas de advecção-difusão usando estabilização Streamline Upwind/Petrov-Galerkin (SUPG), observando que a têmpera é necessária apenas para o termo de estabilização.
4. Análise de Robustez: Experimentos numéricos identificam que o método é robusto em relação à constante $C$ na relação $J_{min} = C h^{d+1}$. O erro permanece próximo ao ótimo sobre uma ampla gama de valores de $C$ (abrangendo quase uma ordem de magnitude), e o condicionamento do sistema linear permanece gerenciável, especialmente com precondicionadores padrão.

Resultados

• Experimentos Numéricos: Em problemas de Poisson 2D e 3D com soluções manufaturadas, o TFEM em malhas com bandas de capas recupera as taxas de convergência do FEM padrão em malhas regulares. Especificamente, escolher $J_{min} \sim h^3$ em 2D e $J_{min} \sim h^4$ em 3D produz convergência de erro ótima em $L^2$ e $H^1$.
• Sensibilidade do Erro: O erro $H^1$ é menos sensível à escolha de $J_{min}$ do que o erro $L^2$, exibindo um patamar sobre uma ampla gama de parâmetros.
• Condicionamento: Embora o número de condição da matriz de rigidez aumente conforme a largura da banda diminui, o uso de precondicionamento de Cholesky incompleto garante que a convergência do solver iterativo permaneça independente da espessura da banda.
• Testes Aleatórios: Uma suíte de 100 testes aleatórios em 2D confirmou que, mesmo com uma escolha simples ($C=1$), os erros do TFEM são quase idênticos aos erros do FEM padrão em malhas regulares.

Significância e Alegações
O artigo afirma que o TFEM aumenta significativamente a flexibilidade da análise de elementos finitos ao relaxar os rigorosos requisitos de qualidade de malha associados a elementos degenerados. Ele permite o uso de malhas com elementos de medida zero ou quase degenerados (como as geradas por rastreamento dinâmico de interface) sem sofrer dos fenômenos de travamento que tornam o FEM padrão inoperante.

Os autores enfatizam que o método é trivial de implementar em códigos de FEM existentes, exigindo apenas uma modificação na avaliação do Jacobiano para evitar a divisão por zero. Eles posicionam o TFEM não como um substituto para o FEM padrão em malhas de alta qualidade, mas como uma extensão robusta que permite o uso de topologias de malha específicas (bandas de capas) que anteriormente eram consideradas inoperantes. O método é particularmente valioso para aplicações como o método de rastreamento de interface X-Mesh, onde manter uma malha conforme é difícil ou computacionalmente caro. O artigo conclui que, embora a escolha teórica ótima do parâmetro $C$ seja derivada, a robustez prática do método permite uma ampla gama de valores aceitáveis, tornando-o imediatamente aplicável.