Stability of trapped fluid clusters in two-phase porous media flow
Abstract
A key challenge in multiphase flow through porous media is to understand and predict the conditions under which trapped fluid clusters become mobilised. Here, we investigate the stability of such clusters in two-phase flow and present a simple, quasistatic model that accurately determines the critical Bond number (that is, the critical ratio between the average pressure gradient of the flow and the surface tension) for the onset of cluster mobilisation. The model is derived by combining elementary geometrical considerations with mass conservation and a mechanical equilibrium condition, resulting in a system of coupled differential equations. Our derivation sheds new light on the mechanisms that govern cluster stability. In addition, since the number of equations equals the number of cluster openings, our model is significantly faster than direct numerical simulations of the same problem and enables efficient exploration of the system’s parameter space. Using this approach, we highlight a discrepancy with the prediction of current mean-field theories, which predict that the largest stable cluster size scales in proportion to r/Boα, where r is a typical pore size, Bo is the Bond number and α is a fixed exponent. We discuss the mechanisms that explain the breakdown of the mean field theories and we show that a scaling law of this form can only exist if α is allowed to depend on a broad set of flow and geometric parameters.
Related
- Engineered three-phase interfaces in porous electrodes for efficient in-situ H2O2 generation and synchronous removal of co-contaminants
- A thermodynamically consistent theory of multicomponent porous media: finite-deformation mechanics, compositional flow, electrostatics, reactions, and phase transformations
- Three thresholds and a computer-assisted phase diagram for the spectral stability of periodic hydroelastic Whitham waves
- Orthogonality of coastal trapped waves
- Flow-rate controls on trapped hydrogen redistribution in a heterogeneous sandstone: a synchrotron X-ray micro-CT study
BibTeX
@article{klahn2026stability,
title = {{Stability of trapped fluid clusters in two-phase porous media flow}},
author = {Mathias Klahn and Gaute Linga and Tanguy Le Borgne and Joachim Mathiesen},
journal = {Journal of Fluid Mechanics},
year = {2026},
month = oct,
doi = {10.1017/jfm.2026.11807},
url = {https://hal.science/hal-05814111v1},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Stability of trapped fluid clusters in two-phase porous media flow
AU - Mathias Klahn
AU - Gaute Linga
AU - Tanguy Le Borgne
AU - Joachim Mathiesen
PY - 2026
DA - 2026/10/19
JO - Journal of Fluid Mechanics
DO - 10.1017/jfm.2026.11807
UR - https://hal.science/hal-05814111v1
AB - A key challenge in multiphase flow through porous media is to understand and predict the conditions under which trapped fluid clusters become mobilised. Here, we investigate the stability of such clusters in two-phase flow and present a simple, quasistatic model that accurately determines the critical Bond number (that is, the critical ratio between the average pressure gradient of the flow and the surface tension) for the onset of cluster mobilisation. The model is derived by combining elementary geometrical considerations with mass conservation and a mechanical equilibrium condition, resulting in a system of coupled differential equations. Our derivation sheds new light on the mechanisms that govern cluster stability. In addition, since the number of equations equals the number of cluster openings, our model is significantly faster than direct numerical simulations of the same problem and enables efficient exploration of the system’s parameter space. Using this approach, we highlight a discrepancy with the prediction of current mean-field theories, which predict that the largest stable cluster size scales in proportion to r/Boα, where r is a typical pore size, Bo is the Bond number and α is a fixed exponent. We discuss the mechanisms that explain the breakdown of the mean field theories and we show that a scaling law of this form can only exist if α is allowed to depend on a broad set of flow and geometric parameters.
SN - 0022-1120
SN - 1469-7645
ER -