A stabilized dual-SAV parametric finite element framework for constrained planar geometric flows with mesh regularization

Resumo Técnico: Um Framework de Elementos Finitos Paramétricos Dual-SAV Estabilizado para Fluxos Geométricos Planos Confinados

Declaração do Problema
As discretizações por elementos finitos paramétricos de fluxos geométricos confinados enfrentam três desafios simultâneos: o tratamento da rigidez geométrica de alta ordem (que restringe os passos de tempo em esquemas explícitos), o gerenciamento da degeneração da malha (agrupamento ou dispersão de vértices) e a imposição de restrições globais não lineares (como preservação de área ou comprimento). Abordagens existentes frequentemente lutam para combinar regularização dinâmica da malha, passo de tempo semi-implícito estável em energia e imposição de restrições não lineares dentro de um único framework linear de grande escala. Métodos totalmente implícitos são robustos, mas exigem a resolução de grandes sistemas não lineares, enquanto tratamentos explícitos sofrem de restrições severas de estabilidade. Além disso, incorporar velocidades tangenciais artificiais para redistribuição da malha em esquemas que preservam estrutura, sem introduzir forças motrizes normais espúrias, permanece delicado.

Metodologia
O artigo propõe um framework de elementos finitos paramétricos (PFEM) dual-variável auxiliar escalar (SAV) estabilizado para curvas planas fechadas. A inovação central é a separação da energia geométrica física e de uma energia artificial de regularização de malha em duas variáveis auxiliares escalares distintas ($r_g$ e $r_m$).

1. Estrutura de Dupla Energia:

   • Energia Geométrica Física ($E_{geom}$): Impulsiona a evolução normal da curva. Uma variável auxiliar $r_g$ é introduzida para linearizar as forças geométricas não lineares.
   • Energia de Regularização de Malha ($E_{mesh}$): Uma energia artificial (tipicamente uma energia de Dirichlet ponderada) usada exclusivamente para redistribuir pontos da malha tangencialmente. Uma variável auxiliar separada $r_m$ é introduzida. Crucialmente, a força motriz de $E_{mesh}$ é projetada apenas na direção tangencial. Isso garante que a redistribuição da malha melhore a parametrização sem alterar a evolução geométrica normal prescrita.

2. Discretização Semi-Implícita de Métrica Congelada:
   O esquema emprega uma discretização de tempo semi-implícita onde todos os termos dependentes da geometria (métricas, forças, restrições) são avaliados na curva conhecida no passo de tempo $n$ ("métrica congelada").

   • Estabilização de Ordem Zero: Formas bilineares de estabilização (Laplaciano ou biharmônico) atuam sobre a própria nova velocidade em vez de diferenças de velocidade. Esta escolha é essencial para manter o cancelamento algébrico necessário para estimativas de dissipação de energia e evita propagar o histórico de velocidade.
   • Resposta Espacial Linear: Uma vez fixas as variáveis auxiliares escalares e os multiplicadores de Lagrange, os problemas espaciais para as velocidades normal e tangencial tornam-se lineares.

3. Redução de Bloco Algébrica para Restrições Globais:
   O método lida com $K$ restrições geométricas globais (por exemplo, área, comprimento) via uma redução de bloco algébrica.

   • Os grandes problemas lineares de resposta espacial são resolvidos primeiro para gerar "vetores de resposta".
   • O sistema não linear restante envolve apenas a variável auxiliar geométrica ($r_g$) e os $K$ multiplicadores de Lagrange.
   • Consequentemente, a dimensão do sistema não linear reduzido é $K+1$, independente do número de vértices da malha $N$. Para restrições simultâneas de área e comprimento ($K=2$), isso resulta em um sistema não linear de dimensão 3.

Contribuições Principais
O artigo delineia cinco contribuições específicas:

1. Formulação Dual-SAV: Introdução de variáveis auxiliares independentes para energias física e de malha, com a força da malha estritamente acoplada ao movimento tangencial.
2. Esquema Semi-Implícito de Métrica Congelada: Construção de uma discretização de tempo que produz problemas de resposta espacial lineares enquanto satisfaz estimativas de dissipação de energia modificada discreta.
3. PFEM Totalmente Discreto: Uma formulação que preserva a estrutura de cancelamento algébrico no nível discreto, garantindo dissipação de energia modificada.
4. Tratamento de Restrições Independente de Dimensão: Derivação de uma redução de bloco algébrica que desacopla as grandes resoluções lineares espaciais da imposição de restrições não lineares, reduzindo a resolução não linear para dimensão $K+1$.
5. Validação Numérica: Testes abrangentes em fluxos de segunda e quarta ordem (encurtamento de curva, encurtamento de curva preservando área, difusão de curva e fluxos do tipo Helfrich).

Resultados
Experimentos numéricos demonstram o seguinte:

• Precisão: O fluxo de encurtamento de curva (círculo encolhendo) exibe precisão temporal aproximada de primeira ordem, consistente com a estrutura de Euler para trás semi-implícita.
• Qualidade da Malha: No encurtamento de curva preservando área, a redistribuição tangencial proposta previne o agrupamento severo de vértices observado em casos de referência Lagrangiana, mantendo uma razão de malha quase uniforme ($Q_{mesh} \approx 1$) enquanto impõe estritamente restrições de área.
• Métricas Não Locais: O framework incorpora com sucesso a métrica não local $H^{-1}$ para fluxos de difusão de curva sem exigir restrições de área explícitas via multiplicadores de Lagrange (a área é naturalmente preservada pela propriedade de média zero da velocidade $H^{-1}$).
• Restrições de Alta Ordem: Para fluxos do tipo Helfrich com restrições simultâneas de área e comprimento ($K=2$), o sistema não linear reduzido (dimensão 3) é resolvido eficientemente. As restrições são impostas com precisão de máquina (erros $\sim 10^{-15}$), e a energia SAV geométrica modificada dissipa monotonicamente.
• Escalabilidade: A análise de custo computacional mostra que, à medida que o número de vértices $N$ aumenta, o número de iterações de Newton para o sistema não linear reduzido permanece constante (aprox. 2,2), enquanto o custo escala linearmente com o número de resoluções de resposta linear.

Significado e Afirmações
O artigo afirma que o significado primário deste trabalho não reside nos componentes individuais (velocidade tangencial, SAV ou tratamento de restrições), mas em seu desacoplamento algébrico dentro de um único framework de elementos finitos paramétricos. Ao separar o fluxo geométrico físico do mecanismo artificial de malha e reduzir a imposição de restrições não lineares a um sistema de baixa dimensão, o método alcança:

• Estabilidade: Tratamento robusto de rigidez de alta ordem via estabilização de ordem zero e passo semi-implícito.
• Eficiência: Evitação de grandes resoluções de ponto de sela não lineares; o custo não linear é independente da resolução da malha.
• Preservação de Estrutura: Manutenção da dissipação de energia modificada discreta e imposição precisa de restrições geométricas globais.

Os autores notam que o escopo atual está restrito a curvas planas fechadas. A extensão do framework para superfícies em evolução em $\mathbb{R}^3$ é identificada como trabalho futuro, exigindo mecanismos adicionais de regularização de malha de superfície e análise de estabilidade. O artigo também reconhece que a relação entre as energias SAV modificadas e as energias geométricas originais merece investigação adicional, particularmente para passos de tempo maiores.