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

Resumen Técnico: Precondicionamiento de Lagrangiano Aumentado para un Modelo de Ericksen–Leslie Simplificado

Planteamiento del Problema
La solución numérica del modelo de Ericksen–Leslie simplificado para cristales líquidos nemáticos presenta desafíos significativos debido al fuerte acoplamiento entre el flujo de fluido (velocidad $u$) y el campo del director ($n$), junto con la imposición simultánea de dos restricciones: la incompresibilidad ($\nabla \cdot u = 0$) y la condición de longitud unitaria puntual ($|n| = 1$). Si bien las formulaciones de multiplicadores de Lagrange evitan el mal condicionamiento asociado con los pequeños parámetros de penalización de Ginzburg–Landau, estas resultan en sistemas de Newton con una compleja estructura de doble punto de silla. Específicamente, el sistema linealizado en cada paso de tiempo es una matriz de bloques $4 \times 4$ que acopla las variables de velocidad–director con las variables de presión–multiplicador. El diseño de solvers iterativos robustos y escalables para estos sistemas sigue siendo un obstáculo computacional crítico.

Metodología
Los autores proponen un enfoque de Lagrangiano Aumentado (AL) totalmente discretizado combinado con un precondicionador de bloques especializado.

1. Formulación: El modelo se reformula utilizando un multiplicador de Lagrange escalar $q$ para imponer la restricción de longitud unitaria, reemplazando la no linealidad cúbica $-|\nabla n|^2 n$ por un término bilineal $qn$. El funcional de Lagrangiano aumentado introduce términos de penalización cuadrática tanto para la restricción de incompresibilidad como para la de longitud unitaria, controlados por los parámetros $\gamma_u$ y $\gamma_n$. Crucialmente, las restricciones se mantienen en forma de multiplicador en lugar de ser reemplazadas por penalizaciones puras, lo que garantiza la precisión mientras mejora las propiedades espectrales de los complementos de Schur.
2. Discretización: El sistema se discretiza en el espacio utilizando elementos finitos Taylor–Hood ($P_2$ para velocidad y director, $P_1$ para presión y multiplicador) y en el tiempo mediante el método de Euler hacia atrás. El sistema algebraico no lineal resultante se linealiza mediante el método de Newton.
3. Estrategia de Precondicionamiento: Se desarrolla un precondicionador de bloques para el sistema de punto de silla resultante. Las incógnitas se agrupan en un bloque primal (incrementos de velocidad y director) y un bloque dual (incrementos de presión y multiplicador).
   • El bloque acoplado velocidad–director se aproxima mediante sus componentes de bloque diagonal.
   • Esto conduce a aproximaciones desacopladas de los dos complementos de Schur de las restricciones:
     • Complemento de Schur de Presión ($S_p$): Se aproxima utilizando una parte principal de Stokes generalizada combinada con la identidad de Sherman–Morrison–Woodbury, produciendo una forma que involucra la matriz de masa de la presión y un término de Laplaciano inverso escalado por el parámetro de aumentación.
     • Complemento de Schur Director–Multiplicador ($S_q$): Se aproxima reconociendo la naturaleza de orden cero del operador de restricción. La aproximación está dominada por el término de penalización, resultando en una inversa de la matriz de masa escalada para el espacio del multiplicador.

Contribuciones Clave

1. Formulación de Lagrangiano Aumentado: Los autores formulan un sistema de AL totalmente discretizado para las ecuaciones de Ericksen–Leslie simplificadas dependientes del tiempo, donde ambas restricciones se imponen mediante multiplicadores, aumentadas por términos de penalización para mejorar el condicionamiento.
2. Derivación del Precondicionador de Bloques: Se deriva un novedoso precondicionador de bloque diagonal que separa la dinámica acoplada de velocidad–director de la imposición de las restricciones. Los autores construyen aproximaciones físicamente escaladas para los complementos de Schur de presión y de director–multiplicador, contabilizando explícitamente las distintas estructuras algebraicas de las restricciones de divergencia y de longitud del director.
3. Evaluación Numérica Exhaustiva: El artículo proporciona una evaluación rigurosa de la precisión del método y la robustez del solver a través de soluciones manufacturadas, pruebas de energía y diversos estudios de parámetros (refinamiento de malla, variación del paso de tiempo, cambios en la viscosidad y geometrías complejas).

Resultados

• Precisión: Las pruebas de solución manufacturada confirman la precisión espacial esperada (aproximadamente un decaimiento de orden tercero en $L_2$ para elementos $P_2$) y la convergencia temporal de primer orden para las variables primarias (velocidad, director, presión). El error del multiplicador alcanza un suelo de error espacial en mallas fijas, consistente con la naturaleza algebraica de la restricción.
• Rendimiento del Solver: Los conteos de iteración de FGMRES externo son casi independientes de la malla y permanecen estables bajo variaciones en el tamaño del paso de tiempo y la viscosidad del fluido.
• Sensibilidad de Parámetros: Los conteos de iteración mejoran a medida que los parámetros de aumentación ($\gamma_u, \gamma_n$) aumentan. El análisis espectral de los bloques de Schur precondicionados muestra que el número de condición del bloque del multiplicador disminuye monotonicamente con el aumento de $\gamma$, mientras que el bloque de presión permanece uniformemente acotado.
• Consistencia Física: Un problema de referencia suave demuestra el decaimiento monótono de la energía total computada, consistente con la ley de disipación de energía continua del modelo de Ericksen–Leslie.
• Complejidad: El método maneja con éxito dominios curvos (disco unitario) y condiciones iniciales complejas que involucran defectos topológicos (configuración de dos defectos) sin modificaciones al precondicionador.

Significancia y Reivindicaciones
El artículo afirma que el precondicionador de bloque de Lagrangiano aumentado propuesto aborda eficazmente el mal condicionamiento asociado con la restricción de longitud unitaria en simulaciones dinámicas de cristales líquidos. Al tratar las restricciones de incompresibilidad y de longitud del director simultáneamente dentro de un marco de AL, el método logra un solver robusto que es ampliamente independiente del tamaño de la malla, el paso de tiempo y la viscosidad. Los autores enfatizan que el diseño del precondicionador se basa en las distintas propiedades espectrales de las dos restricciones, lo que permite un desacoplamiento eficiente.

El trabajo se presenta como un paso significativo hacia la resolución de sistemas de Ericksen–Leslie totalmente discretizados y monolíticos. Los autores señalan modestamente que sus resultados actuales se limitan al régimen incompresible y a rangos de parámetros específicos. Identifican la extensión a flujos nemáticos compresibles (donde las variaciones de densidad son significativas) como una dirección futura natural pero desafiante, que requiere una reevaluación de la estructura de punto de silla. El artículo no pretende una optimalidad universal, sino que demuestra que la estrategia propuesta proporciona una convergencia estable y eficiente para el modelo simplificado probado.