An efficient preconditioner for mixed-dimensional contact poromechanics based on the fixed stress splitting scheme

Resumen Técnico: Un Precondicionador Eficiente para la Poromecánica de Contacto de Dimensiones Mixtas Basado en el Esquema de Desacoplamiento de Esfuerzo Fijo (Fixed Stress Splitting)

Planteamiento del Problema
La simulación numérica de la poromecánica de contacto de fracturas es crítica para aplicaciones como el secuestro de CO2, la producción de energía geotérmica y el almacenamiento subterráneo de gas. Estos problemas involucran una física fuertemente acoplada: flujo de fluidos, deformación de la matriz (poromecánica de Biot) y mecánica de contacto friccional dentro de las fracturas. El modelado preciso es un desafío debido a la compleja geometría de las fracturas (a menudo representadas como objetos de menor dimensión embebidos en un dominio de mayor dimensión) y la naturaleza no suave de las condiciones de contacto (estados de adherencia, deslizamiento y apertura gobernados por la fricción de Coulomb y restricciones de no penetración).

Al resolverse mediante pasos de tiempo totalmente implícitos y el método de Newton, estos problemas generan grandes sistemas lineales totalmente acoplados. La eficiencia y robustez del solver dependen fuertemente del precondicionador. Los métodos existentes enfrentan obstáculos específicos:

1. Estructura de Punto de Silla (Saddle-Point): La formulación de Lagrange Aumentado de la mecánica de contacto introduce una estructura de punto de silla con una submatriz singular, lo que complica la aplicación de técnicas de precondicionamiento estándar.
2. Complejidad de Acoplamiento: El esquema de desacoplamiento de esfuerzo fijo (fixed stress splitting), un método bien establecido para desacoplar la poromecánica, es difícil de extender a medios fracturados con mecánica de contacto. Los solvers de punto de silla estándar suelen asumir estructuras de matriz específicas (por ejemplo, un operador tipo Laplaciano) que no se mantienen cuando el operador encapsula el problema poromecánico acoplado.
3. Escalabilidad: Los precondicionadores deben seguir siendo eficientes a medida que el tamaño de los problemas crece para incluir las extensas redes de fracturas típicas de las simulaciones de reservorios del mundo real.

Metodología
Los autores proponen un precondicionador robusto y escalable que combina aproximaciones de complemento de Schur anidados con una transformación lineal. El enfoque extiende el esquema de esfuerzo fijo tanto a los subdominios de la matriz como a los de la fractura, abordando simultáneamente la naturaleza singular de la mecánica de contacto.

La metodología procede a través de las siguientes etapas:

1. Transformación Lineal (Preprocesamiento): Para manejar la singularidad en la submatriz de la mecánica de contacto (específicamente para los estados de adherencia de la fractura), se aplica una matriz de transformación lineal a la derecha ($Q_r$) al Jacobiano. Esta transformación hace que la submatriz de la mecánica de contacto sea invertible al acoplarla con las ecuaciones de balance de fuerzas de la interfaz, inspirado en técnicas utilizadas en otros problemas de punto de silla.
2. Descomposición de Complemento de Schur Anidado: El precondicionador aproxima el factor triangular superior de la descomposición LDU de bloque del Jacobiano transformado mediante tres niveles de eliminación de complemento de Schur:
   • Nivel 1 (Mecánica de Contacto): Elimina los grados de libertad de la mecánica de contacto. La submatriz de contacto transformada se invierte exactamente (o mediante una aproximación de matriz diagonal de bloque) para desacoplar las condiciones de contacto de las fuerzas de interfaz.
   • Nivel 2 (Balance de Momento): Aplica una aproximación basada en el esfuerzo fijo para desacoplar las ecuaciones de balance de momento. Esto implica la construcción de una matriz de estabilización ($D_{55}$) con coeficientes específicos para la matriz ambiente ($l_{mat}$) y los subdominios de fractura ($l_{frac}$). Las ecuaciones acopladas de momento y de balance de fuerzas de interfaz se resuelven de forma monolítica utilizando Multigrid Algebraico (AMG).
   • Nivel 3 (Flujo de Interfaz): Elimina el flujo de fluido de la interfaz mediante la aproximación de la submatriz de flujo de la interfaz con su diagonal (se utiliza ILU(0) para el sistema de flujo de la interfaz). El complemento de Schur resultante, que representa el balance de masa de fluido, se resuelve de forma monolítica con AMG.
3. Análisis Teórico: Los autores proporcionan una justificación teórica para el esquema de desacoplamiento de esfuerzo fijo en este contexto. Al analizar un modelo linealizado simplificado, demuestran que el esquema iterativo secuencial es una contracción bajo supuestos específicos relacionados con las propiedades de los materiales, los coeficientes de fricción y la estabilidad del estado de contacto. Este análisis valida la elección de los coeficientes de estabilización utilizados en el precondicionador.

Contribuciones Clave

• Diseño de un Precondicionador Novedoso: El artículo introduce un precondicionador que integra con éxito el esquema de esfuerzo fijo con la mecánica de contacto en la poromecánica de dimensiones mixtas. Supera la singularidad de la submatriz de contacto mediante una transformación lineal y maneja el complejo acoplamiento a través de complementos de Schur anidados.
• Prueba de Convergencia Teórica: Los autores derivan garantías de convergencia para el esquema de esfuerzo fijo aplicado a la poromecánica de contacto fracturada. Demuestran que el esquema es una contracción bajo supuestos físicamente razonables, vinculando la tasa de convergencia con los coefócientes de fricción y los parámetros de elasticidad de contacto.
• Escalabilidad y Robustez: El método está diseñado para ser escalable, utilizando AMG para los subproblemas dominantes (elasticidad y balance de masa de fluido) e ILU para los flujos de interfaz, lo que lo hace adecuado para problemas a gran escala con redes de fracturas complejas.

Resultados
El precondicionador fue validado mediante una serie de experimentos numéricos utilizando el marco de trabajo PorePy:

• Sensibilidad de Parámetros: Las pruebas en un modelo de fractura única mostraron que el número de iteraciones depende principalmente del coeficiente de fricción ($F$), donde una mayor fricción requiere más iteraciones, lo cual es consistente con el análisis teórico. La rigidez normal de la fractura ($K_n$) tuvo poco efecto en la convergencia.
• Refinamiento de Malla (2D): En un modelo 2D con una sola fractura, el número de iteraciones de GMRES se mantuvo estable al usar subsolvers directos. Al reemplazar los solvers directos por AMG, el número de iteraciones aumentó ligeramente (por un factor de ~2) a medida que los grados de libertad aumentaron en tres órdenes de magnitud, demostrando una buena escalabilidad.
• Complejidad Geométrica (2D): Los experimentos con múltiples fracturas que se intersectan mostraron que el precondicionador mantuvo su eficiencia, con solo un pequeño incremento en las iteraciones de GMRES para los subsolvers AMG a medida que se refinaba la malla.
• Escalabilidad 3D: Una simulación 3D con múltiples fracturas y casi $10^6$ grados de libertad demostró que el número de iteraciones de GMRES permaneció casi constante con el refinamiento de la malla, confirmando la escalabilidad del método para problemas de gran escala. No se observaron fallos en el solver lineal en el caso 3D.

Significancia
El artículo afirma que el precondicionador propuesto ofrece una solución robusta y escalable para el desafiante problema de la poromecánica de contacto friccional en medios fracturados. Al extender el esquema de esfuerzo fijo para manejar la estructura de punto de silla de la mecánica de contacto mediante la transformación lineal y los complementos de Schur anidados, el método permite la simulación eficiente de procesos subsuperficiales complejos. Los autores enfatizan que el enfoque es aplicable a cualquier discretización de la (poro)mecánica de contacto que utilice enfoques basados en multiplicadores de Lagrange, proporcionando una vía para el modelado preciso y computacionalmente factible de la sismicidad inducida, el transporte de fluidos y la deformación en reservorios fracturados.