Unconstrained Scheme for Geometrically Constrained Gradient Flows

Enunciado do Problema
Este artigo aborda a aproximação numérica de fluxos de gradiente para mapas harmônicos, que são campos vetoriais que minimizam a energia de Dirichlet sujeitos a uma restrição de comprimento unitário pontual ($|u|^2=1$). Esses problemas servem como modelos para aplicações em micromagnetismo, cristais líquidos e flexão de placas não lineares. O desafio principal reside em impor a restrição geométrica durante o passo de tempo. Esquemas existentes "livres de projeção" (ex: Alouges [2], Bartels [7]) linearizam a restrição em cada passo, restringindo a velocidade discreta ao espaço tangente da variedade de restrição. Embora esses métodos ofereçam um controle robusto sobre as violações da restrição, eles exigem a resolução de sistemas de ponto de sela vetoriais degenerados a cada passo de tempo. Esses sistemas podem ser difíceis de precondicionar, particularmente em três dimensões espaciais ou para problemas de quarta ordem como a flexão de placas, onde a estabilidade inf-sup pode depender da regularidade da solução e de degenerações próximas a defeitos.

Metodologia
Os autores propõem um novo esquema "não restrito" que evita a resolução de problemas de ponto de sela, dividindo o cálculo em dois passos distintos:

1. Predição Não Restrita: Um passo de Euler retrocedido não restrito computa uma velocidade prevista $v^{k+1}$. Este passo trata a parte linear do operador de forma implícita e a não linearidade (derivada do operador de projeção $P_u = I - \tilde{u} \otimes \tilde{u}$) de forma explícita.

   • Para fluxos de gradiente $L^2$ (fluxo de calor de mapa harmônico), um termo de estabilização $\gamma(\tilde{u}^k \cdot v^{k+1}, \tilde{u}^k \cdot w)$ é adicionado ao sistema para melhorar a coercividade e controlar o erro de projeção.
   • Para fluxos gerais de espaço de Hilbert (ex: $H^1$ ou $H^2$ para flexão de placas), o esquema pode ser formulado sem esta estabilização ($\gamma=0$) utilizando uma métrica de fluxo de gradiente mais forte.

2. Projeção Pontual: A velocidade prevista é projetada pontualmente no espaço tangente discreto usando um operador de projeção nodal de baixo custo $P_{h, u^k}$.

O sistema linear resultante no primeiro passo consiste em problemas escalares simétricos e positivos definidos (SPD) desacoplados, que são significativamente mais fáceis de resolver e precondicionar do que os sistemas de ponto de sela acoplados exigidos pelos métodos anteriores.

Para garantir a estabilidade, os autores introduzem um critério a posteriori computável baseado na razão entre a norma da velocidade prevista e a norma de seu gradiente projetado. Um algoritmo de passo de tempo adaptativo (Algoritmo 1 e 2) ajusta o tamanho do passo $\tau_k$ dinamicamente para satisfazer esta condição de estabilidade, garantindo a estabilidade de energia e controlando a violação da restrição de comprimento unitário.

Principais Contribuições

1. Eficiência Algorítmica: O esquema proposto substitui a resolução de problemas de ponto de sela vetoriais por sistemas SPD desacoplados. Exemplos computacionais demonstram acelerações de aproximadamente 4–7$\times$ para o fluxo de calor de mapa harmônico e mais de 13$\times$ para problemas de flexão de placas em comparação com os esquemas padrão livres de projeção.
2. Análise de Estabilidade e Convergência: O artigo fornece uma análise completa de estabilidade e convergência para o fluxo de calor de mapa harmônico e para o cálculo de mapas harmônicos. Prova-se que o esquema é energeticamente estável sob restrições leves de passo de tempo e estabelecem-se estimativas de erro para a violação da restrição.
3. Passo de Tempo Adaptativo: Um procedimento de passo de tempo variável é proposto, regido por um critério a posteriori computável, que garante a estabilidade do esquema mesmo na presença de singularidades ou defeitos onde a constante de estabilidade do operador de projeção pode degradar.
4. Extensão para Flexão de Placas: A metodologia é estendida computacionalmente para a flexão de placas não linear (um problema de quarta ordem). Os autores constroem uma projeção tangente nodal específica para elementos de Triângulo de Kirchhoff Discreto (DKT), permitindo que o esquema não restrito lide com a restrição de isometria na flexão de placas de forma eficiente.

Resultados
Experimentos numéricos validam as descobertas teóricas:

• Fluxos Suaves: Em exemplos de fluxo de calor de mapa harmônico suaves, o esquema atinge taxas de convergência ótimas ($O(h+\tau)$ em $H^1$ e $O(h^2+\tau)$ em $L^2$).
• Fluxos Singulares: Em exemplos envolvendo blow-up de tempo finito (fluxo de calor singular), o esquema adaptativo preserva com sucesso a simetria da solução e controla as violações da restrição melhor do que os métodos de passo constante. O esquema não restrito permanece robusto onde a estrutura de ponto de sela dos métodos livres de projeção pode degenerar.
• Cristais Líquidos: Para um problema de casca de cristal líquido 3D, o esquema não restrito com passo de tempo adaptativo reduziu o tempo de computação por um fator de 48 em comparação com o esquema livre de projeção, mantendo níveis comparáveis de energia e violação de restrição.
• Flexão de Placas: Em um exemplo de flexão de placa não linear, o esquema não restrito foi encontrado 11–13 vezes mais rápido que o clássico esquema livre de projeção para passos de tempo menores, com observada convergência linear do defeito de isometria.

Significância
O artigo afirma que o esquema não restrito proposto oferece um aumento substancial na eficiência computacional ao eliminar a necessidade de solvers de ponto de sela, que são frequentemente o gargalo na simulação de fluxos de gradiente geometricamente restritos. Ao depender de sistemas SPD desacoplados e uma estratégia de passo de tempo adaptativo robusta, o método fornece uma alternativa prática e estável para aplicações complexas em micromagnetismo, cristais líquidos e flexão de placas, particularmente em cenários envolvendo altas dimensões espaciais ou operadores diferenciais de ordem superior onde o precondicionamento de sistemas de ponto de sela é desafiador.