Augmented Lagrangian preconditioning for a simplified Ericksen--Leslie model of nematic liquid crystals

Resumo Técnico: Précondicionamento de Lagrange Aumentado para um Modelo de Ericksen–Leslie Simplificado

Definição do Problema
A solução numérica do modelo de Ericksen–Leslie simplificado para cristais líquidos nemáticos apresenta desafios significativos devido ao forte acoplamento entre o fluxo de fluido (velocidade $u$) e o campo de diretor ($n$), juntamente com a imposição simultânea de duas restrições: incompressibilidade ($\nabla \cdot u = 0$) e a condição de comprimento unitário pontual ($|n| = 1$). Embora as formulações de multiplicador de Lagrange evitem o mau condicionamento associado a pequenos parâmetros de penalidade de Ginzburg–Landau, elas resultam em sistemas de Newton com uma estrutura de ponto de sela dupla complexa. Especificamente, o sistema linearizado em cada passo de tempo é uma matriz em blocos $4 \times 4$ que acopla as variáveis velocidade–diretor com as variáveis pressão–multiplicador. Projetar solvers iterativos robustos e escaláveis para esses sistemas continua sendo um obstáculo computacional crítico.

Metodologia
Os autores propõem uma abordagem de Lagrange Aumentado (AL) totalmente discreta combinada com um précondicionador de bloco especializado.

1. Formulação: O modelo é reformulado usando um multiplicador de Lagrange escalar $q$ para impor a restrição de comprimento unitário, substituindo a não linearidade cúbica $-|\nabla n|^2 n$ por um termo bilinear $qn$. O funcional de Lagrange aumentado introduz termos de penalidade quadrática para ambas as restrições de incompressibilidade e de comprimento unitário, controlados pelos parâmetros $\gamma_u$ e $\gamma_n$. Crucialmente, as restrições são mantidas na forma de multiplicadores em vez de serem substituídas por penalidades puras, garantindo precisão enquanto melhoram as propriedades espectrais dos complementos de Schur.
2. Discretização: O sistema é discretizado no espaço usando elementos finitos Taylor–Hood ($P_2$ para velocidade e diretor, $P_1$ para pressão e multiplicador) e no tempo usando o método de Euler para trás (backward Euler). O sistema algébrico não linear resultante é linearizado via o método de Newton.
3. Estratégia de Précondicionamento: Um précondicionador de bloco é desenvolvido para o sistema de ponto de sela resultante. As incógnitas são agrupadas em um bloco primal (incrementos de velocidade e diretor) e um bloco dual (incrementos de pressão e multiplicador).
   • O bloco acoplado velocidade–diretor é aproximado por seus componentes de bloco diagonal.
   • Isso leva a aproximações desacopladas dos dois complementos de Schur de restrição:
     • Complemento de Schur da Pressão ($S_p$): Aproximado usando uma parte principal de Stokes generalizada combinada com a identidade de Sherman–Morrison–Woodbury, resultando em uma forma que envolve a matriz de massa da pressão e um termo de inverso do Laplaciano escalonado pelo parâmetro de aumento.
     • Complemento de Schur Diretor–Multiplicador ($S_q$): Aproximado ao reconhecer a natureza de ordem zero do operador de restrição. A aproximação é dominada pelo termo de penalidade, resultando em uma inversa de matriz de massa escalonada para o espaço do multiplicador.

Principais Contribuições

1. Formulação de Lagrange Aumentado: Os autores formulam um sistema AL totalmente discreto para as equações de Ericksen–Leslie simplificadas dependentes do tempo, onde ambas as restrições são impostas via multiplicadores, aumentadas por termos de penalidade para melhorar o condicionamento.
2. Derivação de Précondicionador de Bloco: Um novo précondicionador de bloco diagonal é derivado, separando a dinâmica acoplada velocidade–diretor da imposição de restrições. Os autores constroem aproximações fisicamente escalonadas para os complementos de Schur de pressão e diretor–multiplicador, contabilizando explicitamente as estruturas algébricas distintas das restrições de divergência e de comprimento unitário.
3. Avaliação Numérica Abrangente: O artigo fornece uma avaliação rigorosa da precisão do método e da robustez do solver através de soluções manufaturadas, benchmarks de energia e vários estudos de parâmetros (refinamento de malha, variação do passo de tempo, mudanças de viscosidade e geometrias complexas).

Resultados

• Precisão: Testes de solução manufaturada confirmam a precisão espacial esperada (decréscimo de ordem três em $L_2$ aproximadamente para elementos $P_2$) e convergência temporal de primeira ordem para as variáveis primárias (velocidade, diretor, pressão). O erro do multiplicador atinge um patamar de erro espacial em malhas fixas, consistente com a natureza algébrica da restrição.
• Desempenho do Solver: As contagens de iteração do FGMRES externo são quase independentes da malha e permanecem estáveis sob variações no tamanho do passo de tempo e na viscosidade do fluido.
• Sensibilidade de Parâmetros: As contagens de iteração melhoram conforme os parâmetros de aumento ($\gamma_u, \gamma_n$) aumentam. A análise espectral dos blocos de Schur pré-condicionados mostra que o número de condição do bloco do multiplicador diminui monotonicamente com o aumento de $\gamma$, enquanto o bloco da pressão permanece uniformemente limitado.
• Consistência Física: Um problema de benchmark suave demonstra o decréscimo monótono da energia total computada, consistente com a lei de dissipação de energia contínua do modelo de Ericksen–Leslie.
• Complexidade: O método lida com sucesso com domínios curvos (disco unitário) e condições iniciais complexas envolvendo defeitos topológicos (configuração de dois defeitos) sem modificação no précondicionador.

Significância e Alegações
O artigo alega que o précondicionador de bloco de Lagrange aumentado proposto aborda efetivamente o mau condicionamento associado à restrição de comprimento unitário em simulações dinâmicas de cristais líquidos. Ao tratar as restrições de incompressibilidade e de comprimento do diretor simultaneamente dentro de uma estrutura de AL, o método alcança um solver robusto que é amplamente independente do tamanho da malha, do passo de tempo e da viscosidade. Os autores enfatizam que o design do précondicionador baseia-se nas propriedades espectrais distintas das duas restrições, permitindo um desacoplamento eficiente.

O trabalho é apresentado como um passo significativo para resolver sistemas de Ericksen–Leslie monolíticos totalmente discretos. Os autores observam modestamente que seus resultados atuais são limitados ao regime incompressível e a intervalos específicos de parâmetros. Eles identificam a extensão para fluxos nemáticos compressíveis (onde as variações de densidade são significativas) como uma direção futura natural, porém desafiadora, exigindo uma reavaliação da estrutura de ponto de sela. O artigo não reivindica otimalidade universal, mas demonstra que a estratégia proposta fornece convergência estável e eficiente para o modelo simplificado testado.