Modified augmented Lagrangian preconditioning for mixed-dimensional beam-solid coupling

Resumen Técnico: Preacondicionamiento de Lagrangiano Aumentado Modificado para el Acoplamiento de Vigas y Sólidos de Dimensiones Mixtas

Planteamiento del Problema
El artículo aborda los desafíos numéricos asociados al acoplamiento de dimensiones mixtas de cuerpos sólidos tridimensionales (3D) y vigas de Kirchhoff–Love de un solo eje (1D) libres de torsión embebidas. Este enfoque de modelado es esencial para simular materiales reforzados con fibras, como compuestos de fibras cortas o tejidos biológicos, donde las fibras son significativamente esbeltas en comparación con la matriz sólida. El acoplamiento se impone mediante multiplicadores de Lagrange para asegurar restricciones posicionales exactas entre la línea central de la viga y el sólido circundante.

Esta formulación conduce a un sistema de ecuaciones lineales con una estructura de punto de silla indefinida. El sistema presenta dificultades específicas:

1. Indefinidez: La presencia de un bloque cero en la diagonal (asociado a los multiplicadores de Lagrange) impide el uso de resolvedores estándar de definición positiva.
2. Subproblemas Singulares: El subproblema de la viga, cuando se considera de forma aislada (como requiere el preacondicionamiento por bloques), es un problema de Neumann puro debido a la falta de condiciones de contorno de Dirichlet explícitas en las fibras embebidas. Esto hace que la matriz de rigidez de la viga sea singular.
3. Sensibilidad de Parámetros: El condicionamiento del sistema está fuertemente influenciado por los parámetros de modelado, incluyendo la relación de rigidez entre la viga y el sólido ($E = E_B/E_S$), la relación de volumen viga-sólido ($V$) y la relación de tamaño de malla.
4. Limitaciones de los Métodos Existentes: Trabajos previos de los autores utilizaron la regularización por penalización, lo que reduce el tamaño del sistema pero introduce un mal condicionamiento e inconsistencia variacional. Por el contrario, los enfoques directos de multiplicadores de Lagrange requieren estrategias de preacondicionamiento especializadas que aún no han sido plenamente exploradas para este contexto específico de elasticidad de valores vectoriales y dimensiones mixtas.

Metodología
Los autores proponen una estrategia de preacondicionamiento por bloques de Lagrangiano Aumentado Modificado (MAL) para resolver el sistema de punto de silla resultante. La metodología procede de la siguiente manera:

1. Formulación Aumentada: El sistema original de multiplicadores de Lagrange se reformula en un sistema de Lagrangiano aumentado. Esto implica añadir términos de tipo penalización ($\epsilon M^T W^{-1} M$ y $\epsilon D^T W^{-1} D$) a los bloques de rigidez del sólido y de la viga. Crucialmente, esta aumentación:

   • Preserva la imposición exacta de las restricciones de acoplamiento (a diferencia de los métodos de penalización pura).
   • Regulariza el subproblema singular de la viga, haciendo que el bloque de rigidez de la viga sea invertible incluso para parámetros de penalización pequeños.
   • Mantiene la estructura de punto de silla del sistema global.

2. Construcción del Preacondicionador de Bloques: Partiendo de un preacondicionador de Lagrangiano aumentado ideal, los autores derivan variantes triangulares por bloques prácticas. Estos preacondicionadores permiten el tratamiento independiente de los bloques del sólido, la viga y el complemento de Schur:

   • Bloque del Sólido: Se preacondiciona utilizando métodos de Multigrid Algebraico (AMG) (específicamente Smoothed Aggregation y Plain Aggregation).
   • Bloque de la Viga: Dado que la matriz aumentada de la viga es de bloque diagonal (para fibras independientes) e invertible, se resuelve utilizando métodos directos (factorización LU).
   • Complemento de Schur: El núcleo del método consiste en aproximar el complemento de Schur aumentado. Los autores introducen tres variantes distintas de aproximación:
     • Variante I: Aproxima el complemento de Schur utilizando inversos aproximados dispersos (SPAI) para el bloque de la viga y una aproximación diagonal para el bloque del sólido, asumiendo un parámetro de penalización pequeño ($\epsilon \ll E_B$).
     • Variante II: Utiliza una aproximación diagonal basada en la matriz de escala $W$ (específicamente $\hat{S} \approx -\frac{1}{\epsilon}W$), lo que requiere un parámetro de penalización cuidadosamente ajustado para agrupar los autovalores.
     • Variante III: Introduce parámetros de penalización separados para las contribuciones del sólido ($\epsilon_S$) y de la viga ($\epsilon_B$). Al establecer $\epsilon_S \to 0$ y ajustar $\epsilon_B$, esta variante busca desacoplar los problemas de escala, utilizando una aproximación de complemento de Schur diagonal.

3. Implementación Numérica: Los sistemas lineales se resuelven mediante GMRES. La implementación aprovecha los marcos de trabajo Trilinos y 4C, utilizando paralelismo de memoria distribuida (MPI).

Contribuciones Clave

• Nuevo Marco de Preacondicionamiento: El artículo presenta la primera aplicación de preacondicionadores de bloques de Lagrangiano aumentado modificados diseñados específicamente para el acoplamiento de vigas y sólidos de dimensiones mixtas con multiplicadores de Lagrange.
• Gestión de Subproblemas Singulares: La estrategia de aumentación regulariza eficazmente el subproblema singular de la viga sin recurrir a técnicas complejas de proyección del espacio nulo o métodos de deflación, permitiendo el uso de resolvedores directos estándar en el bloque de la viga.
• Variantes del Complemento de Schur: La derivación sistemática y la comparación de tres estrategias de aproximación del complemento de Schur (incluyendo la introducción de parámetros de penalización independientes para los dominios del sólido y la viga en la Variante III).
• Análisis de Escalabilidad: Evaluación exhaustiva del escalamiento fuerte y débil en clústeres de memoria distribuida para problemas con millones de grados de libertad.

Resultados
Los experimentos numéricos en elementos de volumen representativos (RVE) y una aplicación de placa sándwich compuesta híbrida arrojaron los siguientes hallazgos:

• Robustez: Los preacondicionadores propuestos demuestran robustez con respecto a la relación de rigidez ($E$), la relación de volumen ($V$) y el refinamiento de la malla. La Variante III, en particular, mostró una robustez casi total respecto a los parámetros para contrastes de rigidez moderados a altos.
• Independencia de la Malla: Los métodos lograron conteos de iteración casi independientes de la malla (típicamente entre 10 y 50 iteraciones para los RVE, y alrededor de 20 para la placa compuesta) a través de diferentes densidades de discretización.
• Escalabilidad: Los métodos exhibieron un escalamiento fuerte y débil favorable en clústeres paralelos (hasta 512 procesadores). Las Variantes II y III mostraron una eficiencia de escalamiento paralelo superior (hasta un 97% de eficiencia de escalamiento débil para la configuración) en comparación con la Variante I, debido principalmente al costo computacional de la aproximación del complemento de Schur en la Variante I.
• Aplicación de Ingeniería: En la aplicación a una placa sándwich compuesta híbrida, el preacondicionador ILU ingenuo falló en converger. En contraste, los preacondicionadores MAL propuestos (especialmente la Variante IIIc) lograron la convergencia con bajos conteos de iteración y tiempos de solución significativamente reducidos (hasta 4 veces más rápido que la variante con peor desempeño) en comparación con los enfoques estándar.

Significancia
El artículo sostiene que los preacondicionadores de bloques de Lagrangiano aumentado modificados propuestos proporcionan una alternativa robusta y eficiente tanto para las formulaciones basadas en penalización pura como para los enfoques estándar de multiplicadores de Lagrange en simulaciones a gran escala de modelos de dimensiones mixtas en mecánica de sólidos y estructural. Al permitir la imposición exacta de las restricciones mientras se regularizan los subproblemas singulares, el método facilita la simulación de estructuras reforzadas con fibras complejas con altos contrastes de rigidez y variaciones en las fracciones de volumen. Los autores concluyen que este enfoque es particularmente adecuado para aplicaciones de ingeniería que involucran fibras embebidas, ofreciendo un camino viable hacia simulaciones de alta fidelidad y escalables. Se señala como trabajo futuro la extensión de este marco a formulaciones de vigas más complejas (por ejemplo, Simo–Reissner) que incluyan efectos de torsión y cortante.