🔬 physics

Modelling flow-driven pore closure of weakening poroelastic media

This paper presents a mathematical model coupling large-deformation poroelasticity, solute transport, and stiffness decay to analyze how chemical weakening influences flow-driven pore closure in materials like hydrogels and rocks, identifying distinct regimes of closure behavior based on the relative timescales of mechanical relaxation and material degradation.

Original authors: Matthew V. Ghosh, Matthew G. Hennessy, Andreas Münch, Sarah L. Waters

Published 2026-08-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Matthew V. Ghosh, Matthew G. Hennessy, Andreas Münch, Sarah L. Waters

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a sponge made not of foam, but of a soft, flexible solid riddled with tiny holes, all soaked with water. This is a poroelastic material, a class of substances found everywhere from the deep rock formations beneath our feet to the artificial scaffolds doctors use to grow new human tissue. These materials have a unique dual nature: they are solid enough to hold their shape, yet porous enough to let fluids flow through them. When you squeeze such a material, the fluid inside is forced to move, and the solid structure deforms in response. This interplay between the flow of liquid and the bending of the solid skeleton is a fundamental force in geology and medicine. However, there is a third, often overlooked player in this system: chemistry. The fluid flowing through these materials often carries dissolved chemicals that can interact with the solid structure itself, making it softer or weaker over time. Understanding how this chemical weakening changes the way the material moves and holds pressure is crucial for predicting whether a rock will collapse or a medical implant will fail.

A team of researchers at the University of Oxford and the University of Bristol has built a new mathematical model to explore exactly this scenario. They focused on a specific, one-dimensional setup: a block of this soft, porous material sitting in a channel, with fluid being pushed through it from one side to the other. As the fluid flows, it carries a chemical solute that slowly eats away at the stiffness of the solid skeleton. The researchers wanted to see how this gradual weakening, combined with the pressure of the flowing fluid, would affect the material's ability to stay open or if it would eventually collapse. Their work reveals that the timing of these processes is everything. They found that the system is governed by four distinct clocks: how fast the solid bounces back after being squeezed, how fast the fluid moves, how fast the chemical spreads, and how fast the material loses its strength. By comparing the speeds of these four processes, the team could predict three very different outcomes for the material.

In some situations, the material settles into a stable state where the fluid flows steadily without the pores ever closing completely. In other cases, the material weakens so much that the pores shut down instantly at one end, causing the flow to stop abruptly. But the most significant discovery, which the researchers uncovered by adding the chemical weakening to their equations, is a middle ground. They found that for certain conditions, the material can hold its shape for a while, only to suddenly collapse and close its pores at a specific, finite moment in time. This "finite-time closure" happens because the material is slowly losing its strength while the fluid continues to push against it. As the solid gets weaker, it can no longer support the pressure, and the pores squeeze shut. This is a new behavior that would not exist if the material remained strong; it is a direct result of the chemical interaction between the fluid and the solid.

To test these ideas, the researchers simulated two real-world examples. The first was a porous polymer scaffold, similar to those used in tissue engineering, which degrades when exposed to simulated body fluid. In this scenario, the material weakens very slowly compared to how fast the fluid moves. Their simulations showed that the material first compresses quickly due to the fluid pressure, and then, over a much longer period, the chemical weakening takes over, eventually leading to pore closure if the pressure is high enough. The second example involved a hydrogel that degrades in the presence of an enzyme. Here, the timescales for weakening and fluid movement were much closer, meaning the chemical and physical changes happened simultaneously, creating a more complex dance of deformation and flow.

The researchers used powerful computer simulations to track these changes, breaking the problem down into manageable pieces to see how the porosity, the stiffness, and the chemical concentration changed over time. They confirmed that when the material weakens slowly, the system behaves in a predictable, quasi-steady way until it hits a critical point. At this point, the material can no longer sustain the pressure drop across it, and the pores collapse. The team derived a mathematical formula to predict exactly when this collapse would happen for slowly weakening materials, showing that the time to closure depends on the initial strength of the material, the pressure applied, and how fast the chemical weakens it.

This work provides a clear mechanistic understanding of why some porous materials fail while others remain stable. It shows that weakening is not just a slow, passive process; it actively changes the mechanical limits of the material. For engineers designing tissue scaffolds, this means they can predict how long a scaffold will last before it collapses under the flow of body fluids. For geologists, it offers a way to understand how chemical reactions in groundwater might cause rock formations to seal off or collapse under pressure. The study does not just describe what happens; it explains why it happens, linking the chemical decay of the material directly to the physical mechanics of fluid flow. By identifying the specific conditions that lead to sudden pore closure, the researchers have provided a tool to anticipate failure in systems where soft solids and flowing fluids meet.

Technical Summary: Modelling Flow-Driven Pore Closure of Weakening Poroelastic Media

Problem Statement
Poroelastic materials, ranging from geological rocks to biomedical polymer scaffolds and hydrogels, often undergo weakening due to interactions between their solid skeleton and chemical species within the interstitial fluid. While Biot poroelasticity effectively models pressure-driven flow and elastic deformation, it does not account for the dynamic reduction of material stiffness caused by solute interaction. This paper addresses the gap in understanding how such weakening affects the time-varying mechanics of the system, specifically focusing on the phenomenon of pore closure—where local porosity decreases to zero. The authors hypothesize that weakening the solid matrix reduces the maximal pressure drop a material can sustain, potentially inducing pore closure at flow rates where it would not otherwise occur, and may lead to the loss of a steady state.

Methodology
The authors develop a novel mathematical framework that couples three distinct physical processes:

Large-deformation poroelasticity: Based on the work of MacMinn et al. (2016), the model uses a neo-Hookean constitutive law for the solid skeleton and Darcy's law for fluid flow, accounting for finite strains.
Solute transport: An advection-diffusion equation (inspired by Fiori et al., 2025) models the transport of a chemical species through the porous medium.
Weakening dynamics: A decay equation is introduced for the Young's modulus ( $E$ ), where the rate of decay is proportional to the local solute concentration. This captures the reduction in stiffness without immediate material degradation (mass loss).

The study focuses on a canonical one-dimensional geometry: a poroelastic material of initial length $L$ confined between a moving free boundary and a rigid permeable wall, subject to a fixed fluid pressure drop. The governing equations are non-dimensionalized, revealing four characteristic timescales: poroelastic relaxation ( $t_{pe}$ ), fluid advection ( $t_{adv}$ ), solute diffusion ( $t_{diff}$ ), and weakening ( $t_{weak}$ ).

The analysis employs both numerical techniques (Finite Element Method using FEniCS) and asymptotic analytical methods. Two specific parameter regimes are investigated:

Example I: A porous polymer scaffold (PLA/PCL) in simulated body fluid, characterized by slow weakening ( $t_{weak} \gg t_{pe}$ ).
Example II: A Dex-MA hydrogel degrading in the presence of an enzyme, where timescales are more comparable.

Key Contributions and Results
The paper identifies three distinct regimes of system behavior based on the interplay between the applied pressure drop, initial porosity, and the minimum attainable stiffness ( $E_{min}$ ):

Steady State Attainment: If the initial porosity is sufficiently high and the pressure drop sufficiently low, the system reaches a steady state where porosity remains positive throughout the domain.
Finite-Time Pore Closure: A novel behavior identified through the inclusion of weakening. Even if the initial conditions allow for a steady state in a non-weakening material, the dynamic reduction of stiffness causes the maximal sustainable pressure drop to decrease over time. Consequently, the pores close at a specific finite time ( $t_{crit}$ ) at the downstream boundary.
Instantaneous Pore Closure: If the applied pressure drop exceeds the limit defined by the initial material properties (even before weakening significantly occurs), the pores close instantaneously at the boundary. This corresponds to the behavior observed in non-weakening materials under excessive load (MacMinn et al., 2016).

In the regime of slow weakening (Example I), the authors derive a leading-order approximation treating the ratio of poroelastic to weakening timescales ( $\epsilon = t_{pe}/t_{weak}$ ) as a small parameter. In this limit:

The solute concentration becomes spatially uniform to leading order.
The Young's modulus decays uniformly in space but evolves in time.
The system behaves quasi-steadily, with the porosity evolution driven by the time-dependent stiffness.
An explicit formula for the time of pore closure ( $t_{crit}$ ) is derived, showing its dependence on the initial porosity, pressure drop, and minimum stiffness.

Numerical simulations confirm that the system exhibits two distinct phases of deformation: an initial rapid compression driven by fluid pressure on the poroelastic timescale, followed by a slower, weakening-dominated deformation phase.

Significance
The paper provides fundamental mechanistic insight into how chemical weakening alters the mechanical stability of poroelastic media. By coupling solute transport with large-deformation poroelasticity, the model demonstrates that weakening is a critical factor in pore closure dynamics. Specifically, the authors show that weakening can induce pore closure at flow rates that would be stable for a non-weakening material, effectively lowering the threshold for system failure.

The identification of finite-time pore closure (Case 2) represents a new behavior arising specifically from the coupling of weakening and flow, distinct from the instantaneous closure seen in static or non-weakening scenarios. The derived analytical approximations for slow weakening provide a practical tool for predicting the time to failure in biomedical applications (such as drug-delivery hydrogels or tissue scaffolds) and geological systems, linking the timescales of chemical interaction to mechanical collapse. The work highlights that the separation of timescales between poroelastic relaxation and weakening is a key determinant of whether a system maintains a steady state or undergoes progressive pore closure.

Technical Summary: Modelling Flow-Driven Pore Closure of Weakening Poroelastic Media

More like this