Linear poroelasticity with solid incompressibility: consistent formulation and scalable numerical solution

Resumen Técnico: Poroelasticidad Lineal con Incompresibilidad del Sólido

Planteamiento del Problema
El artículo aborda el modelado numérico de medios poroelásticos donde la fase sólida es incompresible, una condición relevante para tejidos blandos e hidrogeles. Aunque las ecuaciones de poroelasticidad totalmente no lineales son generales, sus formulaciones primales suelen sufrir problemas de inconsistencia cuando la presión depende de la deformación, particularmente al imponer la incompresibilidad del sólido. Los enfoques estándar suelen basarse en la incompresibilidad total ($J=1$) o la cuasi-incompresibilidad, pero los autores buscan derivar rigurosamente un modelo donde solo la fase sólida sea incompresible, evitando las inconsistencias encontradas en formulaciones primales previas de la poroelasticidad no lineal.

Metodología
Los autores proponen una linealización consistente de las ecuaciones de poroelasticidad totalmente no lineales bajo la hipótesis de incompresibilidad del sólido. La metodología procede a través de los siguientes pasos:

1. Derivación del Modelo: Partiendo de un modelo no lineal con una porosidad Lagrangiana y un potencial de Helmholtz, los autores imponen la incompresibilidad del sólido mediante un multiplicador de Lagrange. Para evitar las inconsistencias derivadas de términos de tercer orden y la necesidad de calcular el Laplaciano de un multiplicador $L^2$, adoptan una formulación de presión mixta. Esto implica definir una variable de "presión total" ($\mu$) y linealizar el sistema alrededor de un punto de equilibrio.
2. Sistema Linealizado: El modelo resultante es un sistema de cuatro campos que involucra:
   • Desplazamiento ($d$)
   • Presión total de poros ($\mu$)
   • Porosidad ($\phi$)
   • Multiplicador de Lagrange de incompresibilidad ($\lambda$)
     El sistema se formula como un problema débil semi-discreto utilizando una discretización temporal implícita, lo que conduce a una estructura de doble punto de silla.
3. Análisis Matemático: Los autores establecen la buena adecuación (well-posedness) de los esquemas semi-discretos y totalmente discretos. Verifican la acotación de los operadores y prueban la condición inf-sup (condiciones de Ladyzhenskaya–Babuška–Brezzi). Un hallazgo clave es que la condición inf-sup discreta requiere únicamente que el espacio de prueba para el multiplicador de Lagrange ($V_{\lambda, h}$) sea un subespacio del espacio de porosidad ($V_{\phi, h}$). Esto permite el uso de pares de Elementos Finitos (FE) de orden más bajo (por ejemplo, $P_1$-continuos para el desplazamiento y $P_0$-discontinuos para la porosidad) sin requerir estabilización compleja o elementos de orden superior.
4. Equivalencia con el Modelo Biot: El artículo demuestra que el sistema propuesto de cuatro campos es matemáticamente equivalente a un modelo Biot estándar de dos campos en el límite de incompresibilidad total, con un parámetro de Lamé volumétrico modificado.
5. Estrategia de Solución Numérica: Para resolver el sistema linealizado de manera eficiente, los autores proponen un precondicionador de bloques escalable basado en el formalismo del complemento de Schur. Aprovechan el método iterativo de esfuerzo constante (fixed-stress) —interpretado como una iteración de Uzawa— para aproximar el complemento de Schur. El precondicionador utiliza una estrategia de partición de bloques:
   • El bloque de desplazamiento se gestiona mediante Multigrid Algebraico (AMG).
   • Los bloques restantes (presión, porosidad, multiplicador) utilizan un enfoque de Jacobi de bloques con AMG para el operador elíptico de presión y Jacobi para las matrices de masa.

Contribuciones Clave

• Formulación Consistente: Un modelo linealizado que impone rigurosamente la incompresibilidad del sólido sin los problemas de inconsistencia encontrados en las formulaciones primales estándar de la poroelasticidad no lineal.
• Estabilidad de Orden Bajo: La formulación admite una aproximación de FE de orden bajo estable (pares de tipo $P_1-P_0$), la cual es computacionalmente eficiente y evita la necesidad de multiplicadores de interfaz o términos de estabilización complejos.
• Solucionador Escalable: El desarrollo de un precondicionador robusto y escalable basado en la aproximación de esfuerzo constante que mantiene conteos de iteración óptimos independientemente del tamaño de la malla o la dimensión del problema (probado en 2D y 3D).
• Validación Teórica: Prueba de la buena adecuación y tasas de convergencia ($O(h^{k+1})$) para la discretización propuesta.

Resultados
Se realizaron pruebas numéricas para validar las afirmaciones teóricas:

• Prueba de Equivalencia: Una comparación del problema de Mandel en estado estacionario mostró que el modelo propuesto de cuatro campos produce resultados indistinguibles del modelo Biot incompresible estándar (diferencias del orden de la precisión de la máquina).
• Pruebas de Convergencia: Pruebas de soluciones manufacturadas en 2D y 3D confirmaron las tasas de convergencia teóricas para el desplazamiento, la presión, la porosidad y el multiplicador de Lagrange utilizando esquemas tanto de primer orden ($k=0$) como de segundo orden ($k=1$).
• Pruebas de Escalabilidad: El rendimiento del precondicionador se evaluó en un benchmark de hinchamiento (swelling) en 3D. El número promedio de iteraciones de GMRES se mantuvo constante (aproximadamente 60–70 iteraciones) a medida que los Grados de Libertad (DoF) aumentaron de ~1,000 a más de 33 millones, demostrando la escalabilidad del solucionador.

Significancia
El artículo afirma que este trabajo proporciona un marco consistente y computacionalmente eficiente para modelar fases sólidas incompresibles en poroelasticidad. Al derivar el modelo directamente de las ecuaciones no lineales y utilizar una variable de presión total, los autores cierran la brecha entre la física no lineal compleja y los solucionadores lineales estables y prácticos. La capacidad de usar elementos de orden bajo reduce el costo computacional, mientras que el precondicionador basado en el esfuerzo constante asegura que la estrategia de solución sea escalable para problemas 3D de gran escala. Los autores señalan que, aunque el modelo está actualmente linealizado, el marco está diseñado para ser consistente con la física no lineal subyacente, y el trabajo futuro explorará la aplicación de esta estrategia a sistemas totalmente no lineales y acoplados (por ejemplo, termo-poroelasticidad).