A Jacobian-free Newton-Krylov method for cell-centred finite volume solid mechanics

Resumen Técnico: Un Método de Newton-Krylov Libre de Jacobiano para la Mecánica de Sólidos de Volúmenes Finitos Centrados en Celdas

Planteamiento del Problema
Los métodos tradicionales de volúmenes finitos (FV) para la mecánica de sólidos dependen predominantemente de enfoques implícitos segregados (iteraciones de Picard) o de esquemas totalmente explícitos. Si bien los métodos segregados ofrecen simplicidad de implementación y eficiencia de memoria, a menudo sufren de una convergencia lenta, particularmente en problemas que involucran no linealidad geométrica, grandes rotaciones o acoplamientos complejos. Por el contrario, los métodos de tipo Newton-Raphson, comunes en la comunidad de elementos finitos, ofrecen convergencia cuadrática pero requieren el ensamblaje, almacenamiento e inversión explícita de la matriz Jacobiana completa. Este requisito es computacionalmente costoso y consume mucha memoria, y adaptar los marcos de trabajo de FV segregados existentes a métodos de Newton de bloque acoplado no es trivial debido a la complejidad de derivar y almacenar el verdadero Jacobiano. En consecuencia, existe la necesidad de un solver de bloque acoplado robusto y eficiente que pueda integrarse en los marcos de trabajo de FV segregados existentes sin una reestructuración mayor del código.

Metodología
Los autores proponen un enfoque de Newton-Krylov libre de Jacobiano (JFNK) diseñado para la mecánica de sólidos de volúmenes finitos centrados en celdas. La metodología central evita la formación explícita de la matriz Jacobiana ($J$) mediante la aproximación de su acción sobre un vector ($v$) utilizando diferencias finitas:
$$Jv \approx \frac{R(u + \epsilon v) - R(u)}{\epsilon}$$
donde $R$ es la función de residuo, $u$ es el vector de solución y $\epsilon$ es un escalar de perturbación. Esta aproximación se combina con el solver de subespacio de Krylov de Mínimos Residuales Generalizados (GMRES) para resolver los sistemas lineales que surgen dentro de las iteraciones de Newton.

Para asegurar la eficiencia y la robustez, el método emplea una estrategia de precondicionamiento. En lugar de utilizar el Jacobiano exacto (que no está disponible), los autores utilizan un Jacobiano aproximado de stencil compacto derivado del término de inercia y un término de difusión (basado en la estabilización de Rhie-Chow utilizada en solvers segregados). Este precondicionador ($P$) es conceptualmente similar a la matriz utilizada en los enfoques segregados tradicionales, lo que permite que el nuevo método aproveche las rutinas de ensamblaje y almacenamiento de matrices existentes dentro de marcos como la herramienta solids4foam de OpenFOAM.

La implementación utiliza la biblioteca PETSc para el solver no lineal (SNES) y el álgebra lineal, interactuando con el solver de FV solids4foam. Las características algorítmicas clave incluyen:

• Precondicionadores: Comparación de solvers directos LU (MUMPS), ILU (Incompleta LU) y Multigrid Algebraico (HYPRE BoomerAMG).
• Globalización: Uso de búsqueda de línea y continuación transitoria (predictor-corrector) para dirigir la solución hacia el dominio de convergencia cuadrática.
• Estabilización: Se añade un término de tipo Rhie-Chow al residuo para suprimir modos de energía cero (efecto checkerboard), con un factor de escala global ($\alpha$) ajustado para la velocidad de convergencia.

Contribuciones Clave

1. Primera Aplicación a la Mecánica de Sólidos de FV: Este trabajo representa la primera propuesta de un método JFNK específicamente para la mecánica de sólidos de volúmenes finitos, cerrando la brecha entre las aplicaciones de dinámica de fluidos y la mecánica de sólidos.
2. Integración Fluida: Al utilizar un Jacobiano aproximado de stencil compacto como precondicionador, el método permite la adopción de la convergencia de tipo Newton dentro de los marcos de trabajo de FV segregados existentes con una modificación mínima del código.
3. Benchmarking Exhaustivo: El método es probado rigurosamente frente a un enfoque segregado estándar a través de una serie de casos que varían en dimensionalidad (2-D/3-D), comportamiento de material (elasticidad lineal, hiperelasticidad, elastoplasticidad) y física (estática, dinámica, interacción fluido-estructura).

Resultados
El estudio compara el método JFNK frente al enfoque segregado utilizando métricas de tiempo computacional, uso de memoria y robustez:

• Eficiencia: El método JFNK supera al enfoque segregado en todos los casos elásticos lineales y no lineales. Se observaron aceleraciones de uno a dos órdenes de magnitud en muchos casos, con una aceleración máxima de 341x en un caso de membrana de Cook de elasticidad lineal en 2-D.
• Memoria: Para problemas con menos de $10^4$ grados de libertad, el uso de memoria es comparable al del método segregado. Para problemas más grandes, el uso de memoria de JFNK aumenta, alcanzando un máximo de 3.6 veces el del método segregado (en un caso de cavidad esférica), debido principalmente a los requisitos de almacenamiento del solver GMRES y el precondicionador de multigrid.
• Robustez y Limitaciones:
  • El método JFNK demostró una robustez superior en casos hiperelásticos (grandes deformaciones/rotaciones) donde el método segregado falló en converger.
  • Sin embargo, el método JFNK divergió en casos elastoplásticos (membrana de Cook con plasticidad). Los autores atribuyen esto al precondicionador de stencil compacto, el cual, basado en la elasticidad de pequeñas deformaciones, es una mala aproximación del Jacobiano para grandes deformaciones plásticas.
• Elección del Precondicionador:
  • Los Solvers Directos LU fueron los más rápidos para problemas de tamaño pequeño a moderado y casos en 2-D.
  • El Multigrid (AMG) fue más efectivo para problemas 3-D de mayor tamaño, ofreciendo mejor escalabilidad y un menor crecimiento de memoria que LU.
  • ILU(k) fue generalmente el más lento y menos robusto.
• Estabilización: Aumentar el factor de estabilización de Rhie-Chow ($\alpha$) generalmente mejoró la velocidad y la robustez de la convergencia, aunque a costa de la precisión pre-asintótica.

Significancia y Reivindicaciones
El artículo sostiene que los métodos de Newton-Krylov libres de Jacobiano son un avance prometedor para la mecánica de sólidos de volúmenes finitos, particularmente para los marcos de trabajo segregados existentes donde modificaciones mínimas permiten su adopción. El estudio demuestra que JFNK puede lograr aceleraciones significativas y una mejor robustez para problemas elásticos sin el costo prohibitivo de memoria del ensamblaje de un Jacobiano completo.

Los autores mantienen la modestia respecto a las limitaciones actuales del método, señalando explícitamente que el enfoque aún no es robusto para la elastoplasticidad con la estrategia de precondicionamiento actual. Identifican el desarrollo de precondicionadores mejorados para la plasticidad y la inclusión de correcciones de no ortogonalidad como trabajos futuros necesarios. La implementación se ha puesto a disposición pública en la herramienta solids4foam para fomentar la exploración y extensión de la comunidad.