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

Resumen Técnico: Un Marco de Elementos Finitos Paramétricos Dual-SAV Estabilizado para Flujos Geométricos Planos Constrained

Planteamiento del Problema
Las discretizaciones por elementos finitos paramétricos de flujos geométricos con restricciones enfrentan tres desafíos simultáneos: el manejo de rigidez geométrica de alto orden (que restringe los pasos de tiempo en esquemas explícitos), la gestión de la degeneración de la malla (agrupamiento o dispersión de vértices) y la imposición de restricciones globales no lineales (como la preservación de área o longitud). Los enfoques existentes a menudo luchan por combinar la regularización dinámica de la malla, la integración temporal semi-implícita estable en energía y la imposición de restricciones no lineales dentro de un único marco lineal a gran escala. Los métodos totalmente implícitos son robustos pero requieren resolver grandes sistemas no lineales, mientras que los tratamientos explícitos sufren restricciones de estabilidad severas. Además, incorporar velocidades tangenciales artificiales para la redistribución de la malla en esquemas que preservan la estructura, sin introducir fuerzas impulsoras normales espurias, sigue siendo delicado.

Metodología
El artículo propone un marco de elementos finitos paramétricos (PFEM) dual-Escalar Variable Auxiliar (SAV) estabilizado para curvas planas cerradas. La innovación central es la separación de la energía geométrica física y una energía artificial de regularización de malla en dos variables auxiliares escalares distintas ($r_g$ y $r_m$).

1. Estructura de Doble Energía:

   • Energía Geométrica Física ($E_{geom}$): Impulsa la evolución normal de la curva. Se introduce una variable auxiliar $r_g$ para linealizar las fuerzas geométricas no lineales.
   • Energía de Regularización de Malla ($E_{mesh}$): Una energía artificial (típicamente una energía de Dirichlet ponderada) utilizada exclusivamente para redistribuir los puntos de la malla tangencialmente. Se introduce una variable auxiliar separada $r_m$. Crucialmente, la fuerza impulsora de $E_{mesh}$ se proyecta únicamente sobre la dirección tangencial. Esto asegura que la redistribución de la malla mejore la parametrización sin alterar la evolución geométrica normal prescrita.

2. Discretización Semi-Implícita de Métrica Congelada:
   El esquema emplea una discretización temporal semi-implícita donde todos los términos dependientes de la geometría (métricas, fuerzas, restricciones) se evalúan sobre la curva conocida en el paso de tiempo $n$ ("métrica congelada").

   • Estabilización de Orden Cero: Las formas bilineales de estabilización (Laplaciano o biarmónico) actúan sobre la nueva velocidad en sí misma en lugar de sobre las diferencias de velocidad. Esta elección es esencial para mantener la cancelación algebraica requerida para las estimaciones de disipación de energía y evita propagar el historial de velocidades.
   • Respuesta Espacial Lineal: Una vez fijadas las variables auxiliares escalares y los multiplicadores de Lagrange, los problemas espaciales para las velocidades normal y tangencial se vuelven lineales.

3. Reducción de Bloque Algebraico para Restricciones Globales:
   El método maneja $K$ restricciones geométricas globales (por ejemplo, área, longitud) mediante una reducción de bloque algebraica.

   • Primero se resuelven los grandes problemas espaciales lineales de respuesta para generar "vectores de respuesta".
   • El sistema no lineal restante involucra únicamente la variable auxiliar geométrica ($r_g$) y los $K$ multiplicadores de Lagrange.
   • En consecuencia, la dimensión del sistema no lineal reducido es $K+1$, independiente del número de vértices de la malla $N$. Para restricciones simultáneas de área y longitud ($K=2$), esto resulta en un sistema no lineal de dimensión 3.

Contribuciones Clave
El artículo describe cinco contribuciones específicas:

1. Formulación Dual-SAV: Introducción de variables auxiliares independientes para las energías física y de malla, con la fuerza de la malla acoplada estrictamente al movimiento tangencial.
2. Esquema Semi-Implícito de Métrica Congelada: Construcción de una discretización temporal que produce problemas espaciales de respuesta lineales mientras satisface estimaciones de disipación de energía modificada discretas.
3. PFEM Totalmente Discreto: Una formulación que preserva la estructura de cancelación algebraica a nivel discreto, asegurando la disipación de energía modificada.
4. Manejo de Restricciones Independiente de la Dimensión: Derivación de una reducción de bloque algebraica que desacopla las grandes resoluciones espaciales lineales de la imposición de restricciones no lineales, reduciendo la resolución no lineal a una dimensión $K+1$.
5. Validación Numérica: Pruebas exhaustivas en flujos de segundo y cuarto orden (acortamiento de curvas, acortamiento de curvas con preservación de área, difusión de curvas y flujos tipo Helfrich).

Resultados
Los experimentos numéricos demuestran lo siguiente:

• Precisión: El flujo de acortamiento de curvas (círculo encogiéndose) exhibe una precisión temporal aproximada de primer orden, consistente con la estructura semi-implícita de Euler hacia atrás.
• Calidad de la Malla: En el acortamiento de curvas con preservación de área, la redistribución tangencial propuesta previene el agrupamiento severo de vértices observado en casos de referencia lagrangianos, manteniendo una relación de malla casi uniforme ($Q_{mesh} \approx 1$) mientras se imponen estrictamente las restricciones de área.
• Métricas No Locales: El marco incorpora exitosamente la métrica no local $H^{-1}$ para flujos de difusión de curvas sin requerir restricciones de área explícitas mediante multiplicadores de Lagrange (el área se preserva naturalmente por la propiedad de media cero de la velocidad $H^{-1}$).
• Restricciones de Alto Orden: Para flujos tipo Helfrich con restricciones simultáneas de área y longitud ($K=2$), el sistema no lineal reducido (dimensión 3) se resuelve eficientemente. Las restricciones se imponen con precisión de máquina (errores $\sim 10^{-15}$), y la energía SAV geométrica modificada disipa monótonamente.
• Escalabilidad: El análisis del costo computacional muestra que a medida que aumenta el número de vértices $N$, el número de iteraciones de Newton para el sistema no lineal reducido permanece constante (aprox. 2.2), mientras que el costo escala linealmente con el número de resoluciones de respuesta lineal.

Significado y Afirmaciones
El artículo afirma que el significado principal de este trabajo no reside en los componentes individuales (velocidad tangencial, SAV o tratamiento de restricciones), sino en su desacoplamiento algebraico dentro de un único marco de elementos finitos paramétricos. Al separar el flujo geométrico físico del mecanismo artificial de la malla y reducir la imposición de restricciones no lineales a un sistema de baja dimensión, el método logra:

• Estabilidad: Manejo robusto de la rigidez de alto orden mediante estabilización de orden cero y pasos semi-implícitos.
• Eficiencia: Evitación de grandes resoluciones de punto de silla no lineales; el costo no lineal es independiente de la resolución de la malla.
• Preservación de Estructura: Mantenimiento de la disipación de energía modificada discreta y la imposición precisa de restricciones geométricas globales.

Los autores señalan que el alcance actual se restringe a curvas planas cerradas. La extensión del marco a superficies evolutivas en $\mathbb{R}^3$ se identifica como trabajo futuro que requiere mecanismos adicionales de regularización de malla superficial y análisis de estabilidad. El artículo también reconoce que la relación entre las energías SAV modificadas y las energías geométricas originales merece una investigación adicional, particularmente para pasos de tiempo más grandes.