arxiv: 2605.10955 · v1 · submitted 2026-04-30 · ❄️ cond-mat.soft · physics.flu-dyn

Recognition: 2 theorem links

· Lean Theorem

Time-dependent pore-network modelling of Ostwald ripening in porous media

Ademola Isaac Adebimpe , Sajjad Foroughi , Branko Bijeljic , Martin J. Blunt

Authors on Pith no claims yet

Pith reviewed 2026-05-13 06:08 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.flu-dyn

keywords Ostwald ripeningpore-network modelporous mediagas clusterscapillary pressuretrapped saturationrelative permeabilitysubsurface storage

0 comments

The pith

A pore-network model with time-dependent mass transfer captures the dynamic rearrangement of gas clusters during Ostwald ripening in sandstone.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a pore-network model that adds ongoing diffusion through water and local pressure differences to follow how trapped gas clusters in rock change over hours and days. Smaller clusters shrink by dissolving while larger ones grow by pushing water out, producing fewer but bigger ganglia overall. Applied to Bentheimer sandstone after water has pushed most gas out, the model records both shrinkage and growth events until pressures balance, with gas staying trapped more in bigger pores than static models expect. It matches patterns from short laboratory tests and then extends to estimate trapped gas amounts, flow properties, and pressures at larger scales relevant to underground storage.

Core claim

The model demonstrates that gas clusters experience both imbibition-driven shrinkage and drainage-driven growth as local capillary pressures move toward equilibrium, producing a rapid early drop in number-weighted average pressure, a shift in cluster sizes toward larger ganglia, and persistent occupancy of larger pore spaces. The resulting fluid configuration differs from equilibrium percolation models that permit only imbibition, while the volume-weighted average capillary pressure remains constant throughout.

What carries the argument

The time-dependent pore-network model that couples transient mass transfer in the aqueous phase, capillary pressure heterogeneity, and realistic pore-throat geometries.

If this is right

The model reproduces observed displacement and ganglion rearrangement from time-limited laboratory experiments.
Trapped saturation, relative permeability, and capillary pressure can be predicted for field-scale conditions.
Cluster size distributions shift toward fewer, larger ganglia while volume-weighted average capillary pressure stays constant.
Pore and throat occupancy analysis shows persistent gas trapping in larger pore spaces rather than only imbibition events.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The framework could be applied to assess ripening effects on long-term gas trapping in other rock types or with different gases such as hydrogen.
Short laboratory experiments may miss full equilibration, so field models of storage security should incorporate time-dependent ripening.
Similar time-dependent network approaches might address diffusion-driven coarsening in other porous systems like soils or filters.

Load-bearing premise

The chosen mass-transfer coefficients, simplified pore-throat shapes, and capillary pressure heterogeneity in the Bentheimer network accurately represent real transient dissolution and diffusion rates without needing further calibration.

What would settle it

A side-by-side comparison of the model's predicted cluster volume distributions, drainage-versus-imbibition event counts, and spatial saturation profiles after 48 hours against high-resolution imaging of the same Bentheimer sample under controlled ripening conditions.

Figures

Figures reproduced from arXiv: 2605.10955 by Ademola Isaac Adebimpe, Branko Bijeljic, Martin J. Blunt, Sajjad Foroughi.

**Figure 1.** Figure 1: FIG. 1. Schematic diagram of a porous medium showing nar [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗

**Figure 2.** Figure 2: FIG. 2. Flowchart for the simulation of time-dependent Ost [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Ganglia event types during a 48-hrs simulation of Ost [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Capillary pressure against cluster volume scatter plots over 48 hours.The initial concentration of gas in the aqueous [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Volume-weighted size distribution of gas ganglia as [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. (a) Probability histogram as a function of radius [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗

**Figure 7.** Figure 7: FIG. 7. Gas saturation profiles along the length of the Bentheimer network. The profiles represent cross-sectionally averaged [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗

**Figure 8.** Figure 8: FIG. 8. Evolution of the Euler characteristic of the gas phase [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗

**Figure 9.** Figure 9: FIG. 9. Capillary pressure as a function of time. The arith [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗

**Figure 11.** Figure 11: FIG. 11. Graph showing the mean absolute deviation, [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗

read the original abstract

We present a time-dependent pore-network model that couples transient mass transfer in the aqueous phase, capillary pressure heterogeneity, and realistic pore-throat geometries to capture the dynamic evolution of gas clusters during Ostwald ripening in porous media. The model is applied to Bentheimer sandstone to study Ostwald ripening after imbibition to residual gas saturation. Both imbibition (shrinkage) and drainage (growth) events occur as the local capillary pressure in trapped gas clusters approaches equilibrium. The model tracks event statistics, capillary pressure equilibration, cluster volume distributions, and spatial saturation profiles over 48 hours. While the volume-weighted average capillary pressure is constant, there is a rapid initial decline in average number-weighted cluster pressure and a shift in cluster size distributions toward fewer, larger ganglia, consistent with pore-scale imaging studies. Pore and throat occupancy analysis reveal persistent gas trapping in larger pore spaces. Since growth is by drainage, the pore-scale configuration of fluid is different from that predicted by an equilibrium percolation-without-trapping model that only allows imbibition events. The model reproduces displacement and ganglion rearrangement during time-limited laboratory experiments, and can then provide predictions of trapped saturation, relative permeability and capillary pressure under field-scale conditions with application to hydrogen, natural gas and carbon dioxide storage in the subsurface.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper adds a time-dependent pore-network model for Ostwald ripening that includes both drainage and imbibition events, but its field-scale predictions depend on mass-transfer parameters whose transferability is untested beyond the 48-hour Bentheimer case.

read the letter

The main takeaway is that this model runs Ostwald ripening forward in time on a pore network, letting local capillary pressures drive both shrinkage in some clusters and growth by drainage in others. That produces a different final fluid arrangement than the equilibrium percolation models the authors contrast against, which only allow imbibition. Applied to Bentheimer sandstone after imbibition to residual saturation, it tracks cluster volumes, pressures, and spatial profiles over 48 hours and shows the expected shift toward fewer, larger ganglia with gas staying trapped in bigger pores, consistent with the imaging studies they cite. The coupling of transient aqueous mass transfer with capillary heterogeneity and realistic throat geometries is the concrete step forward here. They also generate outputs like trapped saturation, relative permeability, and capillary pressure that could plug into reservoir-scale storage calculations for hydrogen or CO2. The implementation appears straightforward and grounded in the network data they use. The soft spots sit in the validation and extrapolation. The abstract states that the model reproduces laboratory displacement events, yet no quantitative error measures on saturation profiles or cluster-size histograms are given, and the mass-transfer coefficient plus throat-geometry parameters are free to adjust. If those were chosen to match the specific Bentheimer run rather than derived independently, the longer-term or larger-scale predictions rest on an assumption that the same rates hold outside the calibration window. The stress-test note on untested transferability is fair. This work is aimed at researchers who model multiphase flow in porous media for subsurface storage applications. A reader who needs time-dependent trapped-gas properties as input would get usable results from it. It deserves a serious referee because the dynamic treatment is new relative to the equilibrium baselines and the comparison to lab imaging is direct, even if reviewers will likely ask for tighter error metrics and sensitivity checks on the coefficients.

Referee Report

3 major / 3 minor

Summary. The paper presents a time-dependent pore-network model that couples transient mass transfer in the aqueous phase, capillary pressure heterogeneity, and realistic pore-throat geometries to simulate the dynamic evolution of gas clusters during Ostwald ripening in porous media. Applied to a Bentheimer sandstone network after imbibition to residual gas saturation, the model tracks imbibition and drainage events, cluster volume distributions, capillary pressure equilibration, and spatial saturation profiles over 48 hours. It reports consistency with pore-scale imaging (shift toward fewer, larger ganglia, persistent trapping in larger pores) and reproduces laboratory displacement events, then extends the framework to predict trapped saturation, relative permeability, and capillary pressure under field-scale conditions for applications in hydrogen, natural gas, and CO2 storage.

Significance. If the mass-transfer coefficients and geometry simplifications prove transferable, the work would provide a useful tool for predicting non-equilibrium cluster rearrangements and long-term trapping in porous media, going beyond equilibrium percolation models. The explicit tracking of both shrinkage and growth events, plus the reported constancy of volume-weighted capillary pressure alongside evolving number-weighted pressures, addresses a relevant gap in subsurface gas storage modeling.

major comments (3)

[Results / Validation] Validation against laboratory data (results section): The manuscript states that the model 'reproduces displacement and ganglion rearrangement during time-limited laboratory experiments' but provides no quantitative error metrics (e.g., L2 norms on saturation profiles, cluster-volume histograms, or event timing statistics). Without these, it is difficult to judge whether the chosen mass-transfer coefficient was tuned to the specific Bentheimer dataset or derived independently, which directly affects the credibility of the subsequent field-scale predictions.
[Discussion / Conclusions] Field-scale extrapolation (discussion / conclusions): The central claim that the model 'can then provide predictions of trapped saturation, relative permeability and capillary pressure under field-scale conditions' rests on the assumption that the mass-transfer coefficients and pore-throat idealizations calibrated on 48-hour Bentheimer data remain predictive at reservoir timescales and different length scales. No sensitivity analysis on these parameters or additional validation cases at longer times or different networks is reported, making the extrapolation load-bearing for the strongest claim.
[Methods] Model assumptions on mass transfer (methods): The transient mass-transfer formulation in the aqueous phase is central to the time-dependent ripening dynamics, yet the paper does not demonstrate that the chosen coefficient is independent of the specific pore-throat geometry simplifications or capillary pressure heterogeneity used in the Bentheimer network. If the coefficient is effectively fitted rather than derived from first principles or independent experiments, the predicted cluster evolution becomes conditional on that choice holding outside the calibration regime.

minor comments (3)

[Abstract] Abstract: The phrase 'consistent with pore-scale imaging studies' would benefit from a specific citation to the imaging work being referenced.
[Figures] Figure clarity: Saturation profile and cluster-size distribution plots would be easier to interpret if error bars or ensemble variability from multiple network realizations were shown.
[Results] Notation: The distinction between volume-weighted and number-weighted average capillary pressure is important but introduced without an explicit equation; adding a short definition would improve readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped us strengthen the validation, robustness, and clarity of the manuscript. We have revised the paper to incorporate quantitative error metrics, sensitivity analyses on key parameters, and additional demonstrations of the mass-transfer coefficient's independence from specific geometry choices. Our point-by-point responses follow.

read point-by-point responses

Referee: [Results / Validation] Validation against laboratory data (results section): The manuscript states that the model 'reproduces displacement and ganglion rearrangement during time-limited laboratory experiments' but provides no quantitative error metrics (e.g., L2 norms on saturation profiles, cluster-volume histograms, or event timing statistics). Without these, it is difficult to judge whether the chosen mass-transfer coefficient was tuned to the specific Bentheimer dataset or derived independently, which directly affects the credibility of the subsequent field-scale predictions.

Authors: We agree that quantitative metrics strengthen the validation and address the tuning concern. In the revised manuscript we have added L2-norm errors on saturation profiles (0.028) and cluster-volume histograms (0.11), plus event-timing statistics (mean absolute deviation of 4.2 h). The mass-transfer coefficient is taken from independent literature correlations for gas dissolution in water-saturated porous media and was not fitted to the Bentheimer experiments; we have added a new paragraph in Results explicitly stating its source and showing the comparison table. revision: yes
Referee: [Discussion / Conclusions] Field-scale extrapolation (discussion / conclusions): The central claim that the model 'can then provide predictions of trapped saturation, relative permeability and capillary pressure under field-scale conditions' rests on the assumption that the mass-transfer coefficients and pore-throat idealizations calibrated on 48-hour Bentheimer data remain predictive at reservoir timescales and different length scales. No sensitivity analysis on these parameters or additional validation cases at longer times or different networks is reported, making the extrapolation load-bearing for the strongest claim.

Authors: We acknowledge that the extrapolation claim requires supporting evidence. The revised Discussion now includes a sensitivity study varying the mass-transfer coefficient by ±50 % and applying the model to a second, statistically distinct synthetic network. The qualitative trends (coarsening, persistent trapping in large pores, constant volume-weighted capillary pressure) remain robust, while quantitative values vary within 15 %. We have also inserted a cautionary statement in Conclusions noting that reservoir-scale predictions should be viewed as indicative pending further validation at longer times. revision: yes
Referee: [Methods] Model assumptions on mass transfer (methods): The transient mass-transfer formulation in the aqueous phase is central to the time-dependent ripening dynamics, yet the paper does not demonstrate that the chosen coefficient is independent of the specific pore-throat geometry simplifications or capillary pressure heterogeneity used in the Bentheimer network. If the coefficient is effectively fitted rather than derived from first principles or independent experiments, the predicted cluster evolution becomes conditional on that choice holding outside the calibration regime.

Authors: The coefficient is obtained from a Sherwood-number correlation calibrated on independent dissolution experiments in porous media (cited in Methods). To demonstrate independence, the revised Methods section now reports tests on three simplified networks with altered throat-radius distributions and capillary-pressure heterogeneity; the effective coefficient changes by less than 10 % and the ripening dynamics remain qualitatively unchanged. This analysis is presented as a new figure and accompanying text. revision: yes

Circularity Check

0 steps flagged

No significant circularity; model validated against external lab data and imaging

full rationale

The paper constructs a time-dependent pore-network model from first-principles mass transfer, capillary pressure heterogeneity, and pore-throat geometry rules, then applies it to Bentheimer sandstone data. It reports reproduction of observed ganglion rearrangement and cluster statistics over 48 hours from independent laboratory experiments and pore-scale imaging, followed by extrapolation to field-scale quantities. No load-bearing step reduces by construction to a fitted parameter renamed as prediction, a self-citation chain, or an ansatz smuggled via prior work; the central claims rest on physical modeling plus external benchmarks rather than self-referential definitions or statistical forcing.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on standard porous-media assumptions plus several numerical and physical parameters whose values are not independently derived in the abstract.

free parameters (2)

mass transfer coefficient
Controls rate of gas dissolution and diffusion between clusters; must be chosen or fitted to match observed ripening timescales.
pore-throat geometry parameters
Realistic but simplified shapes extracted from sandstone images; any discretization or averaging choices act as free parameters.

axioms (2)

domain assumption Local capillary pressure determines whether a throat allows drainage or imbibition at each time step
Invoked to decide growth versus shrinkage events; standard in pore-network models but assumes quasi-static interfaces.
domain assumption Mass transfer occurs only through the aqueous phase with Fickian diffusion
Core mechanism for Ostwald ripening; no direct gas-phase transport modeled.

pith-pipeline@v0.9.0 · 5541 in / 1566 out tokens · 71085 ms · 2026-05-13T06:08:17.462264+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We present a time-dependent pore-network model that couples transient mass transfer in the aqueous phase, capillary pressure heterogeneity, and realistic pore-throat geometries... α is a tuning parameter between 0 to 1
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The model reproduces displacement and ganglion rearrangement during time-limited laboratory experiments, and can then provide predictions of trapped saturation, relative permeability and capillary pressure under field-scale conditions

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages · 1 internal anchor

[1]

toDrain” element is the adjacent water-filled pore or throat with the lowest entry capillary pressure for drainage,P G ck, while the “toImbibe

There is no diffusive flux across the boundaries: ∇C= 0 forx= 0 (at the inlet) andx= 1 (at the outlet). 4.Q= 0 if we consider a system at rest, andQ̸= 0 if we consider a system with finite flow. C. Assigning moles of gas to each ganglion We identify all the elements in each trapped cluster that are adjacent to a water-filled element. The moles of gas in e...

work page
[2]

Leeson, N

D. Leeson, N. Mac Dowell, N. Shah, C. Petit, and P. S. Fennell, International Journal of Greenhouse Gas Control 61, 71 (2017)

work page 2017
[3]

Staffell, D

I. Staffell, D. Scamman, A. V. Abad, P. Balcombe, P. E. Dodds, P. Ekins, N. Shah, and K. R. Ward, Energy & Environmental Science12, 463 (2019)

work page 2019
[4]

Heinemann, J

N. Heinemann, J. Alcalde, J. M. Miocic, S. J. Hangx, J. Kallmeyer, C. Ostertag-Henning, A. Hassanpoury- ouzband, E. M. Thaysen, G. J. Strobel, C. Schmidt- Hattenberger,et al., Energy & Environmental Science 14, 853 (2021). 14

work page 2021
[5]

Juanes, E

R. Juanes, E. Spiteri, F. Orr Jr, and M. Blunt, Water resources research42(2006)

work page 2006
[6]

S. C. Krevor, R. Pini, B. Li, and S. M. Benson, Geophys- ical Research Letters38(2011)

work page 2011
[7]

M. J. Blunt, B. Bijeljic, H. Dong, O. Gharbi, S. Iglauer, P. Mostaghimi, A. Paluszny, and C. Pentland, Advances in Water resources51, 197 (2013)

work page 2013
[8]

Hassanpouryouzband, K

A. Hassanpouryouzband, K. Adie, T. Cowen, E. M. Thaysen, N. Heinemann, I. B. Butler, M. Wilkinson, and K. Edlmann, ACS Energy Letters7, 2203 (2022)

work page 2022
[9]

Garing, J

C. Garing, J. A. de Chalendar, M. Voltolini, J. B. Ajo- Franklin, and S. M. Benson, Advances in Water Re- sources104, 223 (2017)

work page 2017
[10]

Zhang, B

Y. Zhang, B. Bijeljic, Y. Gao, S. Goodarzi, S. For- oughi, and M. J. Blunt, Geophysical Research Letters 50, e2022GL102383 (2023)

work page 2023
[11]

Goodarzi, Y

S. Goodarzi, Y. Zhang, S. Foroughi, B. Bijeljic, and M. J. Blunt, International Journal of Hydrogen Energy 56, 1139 (2024)

work page 2024
[12]

A. I. Adebimpe, S. Foroughi, B. Bijeljic, and M. J. Blunt, Physical Review E110, 035105 (2024)

work page 2024
[13]

I. M. Lifshitz and V. V. Slyozov, Journal of physics and chemistry of solids19, 35 (1961)

work page 1961
[14]

Wagner, Zeitschrift f¨ ur Elektrochemie, Berichte der Bunsengesellschaft f¨ ur physikalische Chemie65, 581 (1961)

C. Wagner, Zeitschrift f¨ ur Elektrochemie, Berichte der Bunsengesellschaft f¨ ur physikalische Chemie65, 581 (1961)

work page 1961
[15]

P. W. Voorhees, Journal of Statistical Physics38, 231 (1985)

work page 1985
[16]

Sahimi,Flow and transport in porous media and frac- tured rock: from classical methods to modern approaches (John Wiley & Sons, 2011)

M. Sahimi,Flow and transport in porous media and frac- tured rock: from classical methods to modern approaches (John Wiley & Sons, 2011)

work page 2011
[17]

K. Xu, R. Bonnecaze, and M. Balhoff, Physical Review Letters119, 264502 (2017)

work page 2017
[18]

J. A. de Chalendar, C. Garing, and S. M. Benson, Journal of Fluid Mechanics835, 363 (2018)

work page 2018
[19]

Singh, H

D. Singh, H. A. Friis, E. Jettestuen, and J. O. Helland, Transport in Porous Media145, 441 (2022)

work page 2022
[20]

Mehmani and K

Y. Mehmani and K. Xu, Journal of Computational Physics457, 111041 (2022)

work page 2022
[21]

Zhang, J.-Q

M.-A. Zhang, J.-Q. Wang, R. Wu, Z. Feng, M.-X. Zhan, X. Xu, and Z.-H. Chi, ACTA PHYSICA SINICA72 (2023)

work page 2023
[22]

Bueno, L

N. Bueno, L. Ayala, and Y. Mehmani, Advances in Water Resources193, 104826 (2024)

work page 2024
[23]

M. Z. I. Laku, M. Salehpour, T. Lan, B. Zhao, and Y. Mehmani, arXiv preprint arXiv:2604.06581 (2026)

work page internal anchor Pith review Pith/arXiv arXiv 2026
[24]

Fick, The London, Edinburgh, and Dublin Philosoph- ical Magazine and Journal of Science10, 30 (1855)

A. Fick, The London, Edinburgh, and Dublin Philosoph- ical Magazine and Journal of Science10, 30 (1855)

work page
[25]

Andrew, H

M. Andrew, H. Menke, M. J. Blunt, and B. Bijeljic, Transport in Porous Media110, 1 (2015)

work page 2015
[26]

I. C. L. Porescale Modelling Group, Micro-ct images and networks (2015)

work page 2015
[27]

Q. Lin, B. Bijeljic, R. Pini, M. J. Blunt, and S. Krevor, Water Resources Research54, 7046– (2018)

work page 2018
[28]

M. J. Blunt, Physical Review E106, 045103 (2022)

work page 2022