Modelling flow-driven pore closure of weakening poroelastic media

Resumen Técnico: Modelado del Cierre de Poros Impulsado por el Flujo en Medios Poroelásticos de Debilitamiento

Planteamiento del Problema
Los materiales poroelásticos, que van desde rocas geológicas hasta andamios poliméricos biomédicos e hidrogeles, a menudo experimentan un debilitamiento debido a las interaciones entre su esqueleto sólido y las especies químicas dentro del fluido intersticial. Si bien la poroelasticidad de Biot modela eficazmente el flujo impulsado por la presión y la deformación elástica, no tiene en cuenta la reducción dinámica de la rigidez del material causada por la interacción con solutos. Este artículo aborda la brecha en la comprensión de cómo este debilitamiento afecta la mecánica variable en el tiempo del sistema, centrándose específicamente en el fenómeno del cierre de poros —donde la porosidad local disminuye a cero—. Los autores plantean la hipótesis de que el debilitamiento de la matriz sólida reduce la caída de presión máxima que un material puede soportar, induciendo potencialmente el cierre de los poros en tasas de flujo donde de otro modo no ocurriría, y puede conducir a la pérdida de un estado estacionario.

Metodología
Los autores desarrollan un nuevo marco matemático que acopla tres procesos físicos distintos:

1. Poroelasticidad de grandes deformaciones: Basándose en el trabajo de MacMinn et al. (2016), el modelo utiliza una ley constitutiva neo-Hookeana para el esqueleto sólido y la ley de Darcy para el flujo de fluidos, teniendo en cuenta deformaciones finitas.
2. Transporte de solutos: Una ecuación de advección-difusión (inspirada en Fiori et al., 2025) modela el transporte de una especie química a través del medio poroso.
3. Dinámica de debilitamiento: Se introduce una ecuación de decaimiento para el módulo de Young ($E$), donde la tasa de decaimiento es proporcional a la concentración local de soluto. Esto captura la reducción de la rigidez sin una degradación inmediata del material (pérdida de masa).

El estudio se centra en una geometría canónica unidimensional: un material poroelástico de longitud inicial $L$ confinado entre un límite libre en movimiento y una pared permeable rígida, sujeto a una caída de presión de fluido fija. Las ecuaciones gobernantes se adimensionalizan, revelando cuatro escalas de tiempo características: relajación poroelástica ($t_{pe}$), advección de fluido ($t_{adv}$), difusión de solutos ($t_{diff}$) y debilitamiento ($t_{weak}$).

El análisis emplea tanto técnicas numéricas (Método de Elementos Finitos utilizando FEniCS) como métodos analíticos asintóticos. Se investigan dos regímenes de parámetros específicos:

• Ejemplo I: Un andamio polimérico poroso (PLA/PCL) en fluido corporal simulado, caracterizado por un debilitamiento lento ($t_{weak} \gg t_{pe}$).
• Ejemplo II: Un hidrogel de Dex-MA degradándose en presencia de una enzima, donde las escalas de tiempo son más comparables.

Contribuciones Clave y Resultados
El artículo identifica tres regímenes distintos de comportamiento del sistema basados en la interacción entre la caída de presión aplicada, la porosidad inicial y la rigidez mínima alcanzable ($E_{min}$):

1. Alcance del Estado Estacionario: Si la porosidad inicial es suficientemente alta y la caída de presión suficientemente baja, el sistema alcanza un estado estacionario donde la porosidad permanece positiva en todo el dominio.
2. Cierre de Poros en Tiempo Finito: Un comportamiento novedoso identificado mediante la inclusión del debilitamiento. Incluso si las condiciones iniciales permiten un estado estacionario en un material que no se debilita, la reducción dinámica de la rigidez hace que la caída de presión máxima sostenible disminuya con el tiempo. En consecuencia, los poros se cierran en un tiempo finito específico ($t_{crit}$) en el límite aguas abajo.
3. Cierre de Poros Instantáneo: Si la caída de presión aplicada excede el límite definido por las propiedades iniciales del material (incluso antes de que el debilitamiento ocurra significativamente), los poros se cierran instantáneamente en el límite. Esto corresponde al comportamiento observado en materiales no debilitantes bajo carga excesiva (MacMinn et al., 2016).

En el régimen de debilitamiento lento (Ejemplo I), los autores derivan una aproximación de orden principal tratando la relación entre las escalas de tiempo poroelástica y de debilitamiento ($\epsilon = t_{pe}/t_{weak}$) como un parámetro pequeño. En este límite:

• La concentración de soluto es espacialmente uniforme en el orden principal.
• El módulo de Young decae uniformemente en el espacio pero evoluciona en el tiempo.
• El sistema se comporta cuasi-estacionariamente, con la evolución de la porosidad impulsada por la rigidez dependiente del tiempo.
• Se deriva una fórmula explícita para el tiempo de cierre de poros ($t_{crit}$), mostrando su dependencia de la porosidad inicial, la caída de presión y la rigidez mínima.

Las simulaciones numéricas confirman que el sistema exhibe dos fases distintas de deformación: una compresión inicial rápida impulsada por la presión del fluido en la escala de tiempo poroelástica, seguida de una fase de deformación más lenta dominada por el debilitamiento.

Significancia
El artículo proporciona una visión mecánica fundamental de cómo el debilitamiento químico altera la estabilidad mecánica de los medios poroelásticos. Al acoplar el transporte de solutos con la poroelasticidad de grandes deformaciones, el modelo demuestra que el debilitamiento es un factor crítico en la dinámica de cierre de poros. Específicamente, los autores muestran que el debilitamiento puede inducir el cierre de poros a tasas de flujo que serían estables para un material que no se debilita, reduciendo efectivamente el umbral de falla del sistema.

La identificación del cierre de poros en tiempo finito (Caso 2) representa un nuevo comportamiento que surge específicamente del acoplamiento entre el debilitamiento y el flujo, distinto del cierre instantáneo observado en escenarios estáticos o de no debilitamiento. Las aproximaciones analíticas derivadas para el debilitamiento lento proporcionan una herramienta práctica para predecir el tiempo de falla en aplicaciones biomédicas (como hidrogeles de liberación de fármacos o andamios de tejido) y sistemas geológicos, vinculando las escalas de tiempo de la interacción química con el colapso mecánico. El trabajo destaca que la separación de escalas de tiempo entre la relajación poroelástica y el debilitamiento es un determinante clave de si un sistema mantiene un estado estacionario o experimenta un cierre de poros progresivo.