pith. sign in

arxiv: 2606.24295 · v1 · pith:GSG6B5XLnew · submitted 2026-06-23 · ⚛️ physics.geo-ph

Suppressing Spontaneous Droplet Shrinkage in Cahn-Hilliard-Stokes Microflows

A. M. Gurin This is my paper

Pith reviewed 2026-06-25 22:02 UTC · model grok-4.3

classification ⚛️ physics.geo-ph
keywords Cahn-Hilliarddroplet shrinkagephase-field methodtwo-phase flowmicroflowsporous mediavolume conservationT-junction
0
0 comments X

The pith

A curvature-dependent correction to the double-well potential counteracts spontaneous droplet shrinkage in Cahn-Hilliard microflow simulations.

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

The paper identifies spontaneous droplet shrinkage as a non-physical artifact in Cahn-Hilliard simulations of immiscible two-phase microflows. This artifact stems from the energy functional's bias toward minimizing interfacial area, which shifts equilibrium phase fractions away from 0 and 1. The authors propose a curvature-dependent correction that shifts the double-well potential to counteract this drift. Validation on an isolated droplet in a cavity, drainage in a synthetic porous medium, and flow through a T-junction shows reduced phase fraction overshoot and undershoot along with improved droplet volume conservation. A sympathetic reader would care because accurate phase fractions and volume preservation are needed for reliable predictions of real immiscible flows in confined geometries.

Core claim

The central claim is that a curvature-dependent correction shifting the double-well potential counteracts the drift in equilibrium phase fractions away from 0 and 1 caused by the energy functional's preference for reduced interfacial area.

What carries the argument

A curvature-dependent correction that shifts the double-well potential to maintain phase fractions near their physical bounds of 0 and 1.

If this is right

  • The correction suppresses phase fraction overshoot and undershoot.
  • It improves droplet volume conservation in the tested geometries.
  • The same shift works for an isolated droplet, drainage in porous media, and flow in a T-junction.

Where Pith is reading between the lines

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

  • The correction could be applied to other phase-field formulations that share the same energy bias.
  • It may reduce the need for ad-hoc fixes when scaling simulations to complex three-dimensional geological domains.
  • Direct comparison of droplet volume histories with and without the shift in high-curvature unsteady flows would test whether the fix remains parameter-free.

Load-bearing premise

The spontaneous shrinkage artifact originates solely from the energy functional's preference for reduced interfacial area, and a curvature-dependent shift to the double-well potential will counteract it across flow geometries without new non-physical artifacts.

What would settle it

A long simulation of an isolated droplet in a cubic cavity in which the interior phase fraction remains exactly 1 when the curvature correction is active but drifts measurably when it is removed.

Figures

Figures reproduced from arXiv: 2606.24295 by A. M. Gurin.

Figure 1
Figure 1. Figure 1: (a) Computational grid, (b) initial distribution of the non-wetting [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Wetting-phase distribution in the porous medium: (a) [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Saturation evolution during drainage [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schematic of a droplet in a cubic cavity. [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Phase-fraction distribution in the computational domain at steady [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Schematic of the T-Shaped junction geometry. [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Temporal evolution of the wetting-phase saturation in the T-shaped [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Spatial distribution of the phase fraction in the T-shaped junction at [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Droplet shrinkage schematic We assume that the phase is uniformly distributed within the droplet and that the integral of the phase fraction over the entire computational domain remains constant in time. Initially, inside the droplet, the phase fraction is ϕ = 1, and outside it is ϕ = 0. The diffuse interface, where the phase fraction transitions from 1 to 0, has a characteristic width ε. We assume the ini… view at source ↗
Figure 10
Figure 10. Figure 10: Schematic illustration of the double-well potential shift. [PITH_FULL_IMAGE:figures/full_fig_p016_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Effect of the correction on equilibrium phase fractions. [PITH_FULL_IMAGE:figures/full_fig_p018_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Saturation of the wetting phase in a T-shaped junction. [PITH_FULL_IMAGE:figures/full_fig_p019_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Saturation evolution during drainage with the correction eq. (13) [PITH_FULL_IMAGE:figures/full_fig_p020_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Wetting phase distribution in the porous medium with the correction [PITH_FULL_IMAGE:figures/full_fig_p020_14.png] view at source ↗
read the original abstract

This paper addresses non-physical artifact, specifically spontaneous droplet shrinkage, in Cahn-Hilliard-based simulations of immiscible two-phase microflows. That artifact arise from the energy functional's preference for reduced interfacial area, leading to shifts in equilibrium phase fractions away from the physical bounds of 0 and 1. A curvature-dependent correction is proposed that shifts the double-well potential to counteract this drift. This method is validated on the following test cases: an isolated droplet in a cubic cavity, drainage in a synthetic porous medium, and flow through a T-shaped junction. Results demonstrate significant suppression of phase fraction overshoot/undershoot and improved droplet volume conservation.

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.

Referee Report

3 major / 1 minor

Summary. The manuscript proposes a curvature-dependent correction added to the double-well potential in the Cahn-Hilliard equation for Cahn-Hilliard-Stokes simulations of immiscible two-phase microflows. The correction is intended to counteract the non-physical drift of equilibrium phase fractions away from 0 and 1 that arises from the energy functional's preference for minimizing interfacial area, thereby suppressing spontaneous droplet shrinkage. The approach is validated on three test cases—an isolated droplet in a cubic cavity, drainage in a synthetic porous medium, and flow through a T-shaped junction—with claims of reduced phase-fraction overshoot/undershoot and improved volume conservation.

Significance. If the correction can be shown to preserve the variational structure and exact chemical-potential equilibrium when coupled to Stokes flow without introducing new artifacts or geometry-specific tuning, the method would offer a parameter-free way to improve volume conservation in phase-field models of microflows. This would be particularly useful for applications involving droplets in confined geometries or porous media, where artificial shrinkage currently limits predictive accuracy.

major comments (3)
  1. [Abstract] Abstract: the validation is described only in qualitative terms ('significant suppression of phase fraction overshoot/undershoot and improved droplet volume conservation') with no quantitative metrics, error bars, baseline comparisons to the uncorrected Cahn-Hilliard model, or tabulated volume-conservation errors across the three test cases; this leaves the central claim of effectiveness weakly supported.
  2. [Method / Validation] The skeptic concern is load-bearing: because the correction is added directly to the double-well term, any mismatch between the curvature operator and the mean-curvature contribution already present in the Cahn-Hilliard chemical potential could leave a residual driving force. The manuscript must demonstrate (via derivation or numerical isolation) that chemical-potential equilibrium remains exactly satisfied when the corrected system is coupled to Stokes flow, independent of geometry.
  3. [Validation] The weakest assumption—that the artifact originates solely from interfacial-area preference and that the curvature-dependent shift counteracts it across geometries without new non-physical effects—requires explicit verification. The three test cases do not isolate residual driving forces in a geometry-independent manner (e.g., via a static droplet test with varying curvature).
minor comments (1)
  1. [Abstract] The abstract should include at least one quantitative indicator (e.g., maximum volume error or L2 phase-fraction deviation) for each test case to allow immediate assessment of improvement.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below and have revised the manuscript accordingly to provide quantitative validation, explicit derivations, and additional isolating tests.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the validation is described only in qualitative terms ('significant suppression of phase fraction overshoot/undershoot and improved droplet volume conservation') with no quantitative metrics, error bars, baseline comparisons to the uncorrected Cahn-Hilliard model, or tabulated volume-conservation errors across the three test cases; this leaves the central claim of effectiveness weakly supported.

    Authors: We agree that the original abstract lacked quantitative support. The revised abstract now reports specific metrics, including an 80% reduction in phase-fraction overshoot and volume errors reduced from O(10^{-2}) to O(10^{-4}). We have also added Table 1, which tabulates mean volume-conservation errors (with standard deviations from five independent runs) for both the baseline and corrected models across all three test cases. revision: yes

  2. Referee: [Method / Validation] The skeptic concern is load-bearing: because the correction is added directly to the double-well term, any mismatch between the curvature operator and the mean-curvature contribution already present in the Cahn-Hilliard chemical potential could leave a residual driving force. The manuscript must demonstrate (via derivation or numerical isolation) that chemical-potential equilibrium remains exactly satisfied when the corrected system is coupled to Stokes flow, independent of geometry.

    Authors: This concern is well-taken. We have inserted a new derivation (Section 2.4) showing that the curvature-dependent shift is constructed so the modified chemical potential is identically zero at c=0 and c=1 for arbitrary curvature, preserving equilibrium even under Stokes coupling. We further added a numerical isolation experiment of a static droplet (zero velocity, hence coupled to Stokes) with radii varied over an order of magnitude; phase fractions remain at equilibrium to machine precision with no detectable residual force. revision: yes

  3. Referee: [Validation] The weakest assumption—that the artifact originates solely from interfacial-area preference and that the curvature-dependent shift counteracts it across geometries without new non-physical effects—requires explicit verification. The three test cases do not isolate residual driving forces in a geometry-independent manner (e.g., via a static droplet test with varying curvature).

    Authors: We accept that the original tests did not isolate the effect in a curvature-varied, geometry-independent setting. The revised manuscript includes a dedicated static-droplet subsection in which curvature is varied systematically while holding all other parameters fixed; results confirm that the correction eliminates the driving force without introducing new artifacts, thereby supporting the underlying assumption. revision: yes

Circularity Check

0 steps flagged

Proposed correction is an independent ansatz validated on separate test cases

full rationale

The paper identifies a known artifact arising from the standard Cahn-Hilliard energy functional's interfacial-area preference, then introduces a curvature-dependent shift to the double-well potential as an explicit modeling addition. This addition is not derived from or fitted to the target quantities; instead it is proposed and then checked for efficacy on three independent geometries (isolated droplet, porous-medium drainage, T-junction). No equation reduces to its own input by construction, no parameter is fitted on a subset and then relabeled a prediction, and no load-bearing premise rests on a self-citation chain. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based solely on the abstract, no explicit free parameters, axioms, or invented entities are described beyond the standard Cahn-Hilliard framework; the correction itself is presented as a modeling adjustment without further breakdown.

pith-pipeline@v0.9.1-grok · 5632 in / 1109 out tokens · 43353 ms · 2026-06-25T22:02:33.024264+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

11 extracted references · 7 canonical work pages

  1. [1]

    The Persistent Homology of Dual Digital Image Constructions

    Satish Balay et al. “Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries”. In:Modern Software Tools for Scientific Computing. Boston, MA: Birkh¨ auser Boston, 1997, pp. 163–202.isbn: 978-1-4612-1986-6.doi:10.1007/978- 1- 4612- 1986-6_8.url:https: //doi.org/10.1007/978-1-4612-1986-6_8

    work page doi:10.1007/978- 1997

  2. [2]

    Conservation de la positivité lors de la discrétisation des problèmes d’évolution paraboliques

    Franck Boyer and Sebastian Minjeaud. “Numerical schemes for a three component Cahn-Hilliard model”. en. In:ESAIM: Mathematical Modelling and Numerical Analysis45.4 (2011), pp. 697–738.doi:10.1051/m2an/ 2010072.url:https://www.numdam.org/articles/10.1051/m2an/ 2010072/

    work page doi:10.1051/m2an/ 2011

  3. [3]

    Study of a three component Cahn- Hilliard flow model

    Boyer, Franck and Lapuerta, C´ eline. “Study of a three component Cahn- Hilliard flow model”. In:ESAIM: M2AN40.4 (2006), pp. 653–687.doi: 10 . 1051 / m2an : 2006028.url:https : / / doi . org / 10 . 1051 / m2an : 2006028

    2006

  4. [4]

    A New Lagrange Multiplier Approach for Constructing Structure Preserving Schemes, II. Bound Preserving

    Qing Cheng and Jie Shen. “A New Lagrange Multiplier Approach for Constructing Structure Preserving Schemes, II. Bound Preserving”. In: SIAM Journal on Numerical Analysis60.3 (2022), pp. 970–998.doi:10. 1137/21M144877X. eprint:https://doi.org/10.1137/21M144877X.url: https://doi.org/10.1137/21M144877X

    work page doi:10.1137/21m144877x.url: 2022

  5. [5]

    Ostwald ripening and breakup characteristics of the advective Cahn-Hilliard equation: The role of free energy functionals

    M. F. P. ten Eikelder and Makrand A. Khanwale. “Ostwald ripening and breakup characteristics of the advective Cahn-Hilliard equation: The role of free energy functionals”. In:Center for Turbulence Research Proceedings of the Summer Program. 2024

    2024

  6. [6]

    Application of the Parallel INMOST Platform to Sub- surface Flow and Transport Modelling

    Igor Konshin et al. “Application of the Parallel INMOST Platform to Sub- surface Flow and Transport Modelling”. In:Parallel Processing and Ap- plied Mathematics. Cham: Springer International Publishing, 2016, pp. 277– 286.isbn: 978-3-319-32152-3. 21

    2016

  7. [7]

    A pressure-correction and bound-preserving discretization of the phase-field method for variable density two-phase flows

    Chen Liu et al. “A pressure-correction and bound-preserving discretization of the phase-field method for variable density two-phase flows”. In:Journal of Computational Physics449 (2022), p. 110769.issn: 0021-9991.doi: https://doi.org/10.1016/j.jcp.2021.110769.url:https://www. sciencedirect.com/science/article/pii/S0021999121006641

    work page doi:10.1016/j.jcp.2021.110769.url:https://www 2022

  8. [8]

    An efficient finite ele- ment method for simulation of droplet spreading on a topologically rough surface

    Li Luo, Xiao-Ping Wang, and Xiao-Chuan Cai. “An efficient finite ele- ment method for simulation of droplet spreading on a topologically rough surface”. In:Journal of Computational Physics349 (2017), pp. 233–252. issn: 0021-9991.doi:https : / / doi . org / 10 . 1016 / j . jcp . 2017 . 08 . 010.url:https://www.sciencedirect.com/science/article/pii/ S002199...

    2017

  9. [9]

    Algebraic Multigrid for Stokes Equations

    Yvan Notay. “Algebraic Multigrid for Stokes Equations”. In:SIAM Jour- nal on Scientific Computing39.5 (2017), S88–S111.doi:10.1137/16M1071419. eprint:https://doi.org/10.1137/16M1071419.url:https://doi. org/10.1137/16M1071419

    work page doi:10.1137/16m1071419 2017

  10. [10]

    Spontaneous shrinkage of drops and mass conservation in phase-field simulations

    Pengtao Yue, Chunfeng Zhou, and James J. Feng. “Spontaneous shrinkage of drops and mass conservation in phase-field simulations”. In:Journal of Computational Physics223.1 (2007), pp. 1–9.issn: 0021-9991.doi: https://doi.org/10.1016/j.jcp.2006.11.020.url:https://www. sciencedirect.com/science/article/pii/S0021999106005857

    work page doi:10.1016/j.jcp.2006.11.020.url:https://www 2007

  11. [11]

    Spontaneous shrinkage of droplet on a wetting surface in the phase-field model

    Chunhua Zhang and Zhaoli Guo. “Spontaneous shrinkage of droplet on a wetting surface in the phase-field model”. In:Phys. Rev. E100 (6 Dec. 2019), p. 061302.doi:10 . 1103 / PhysRevE . 100 . 061302.url:https : //link.aps.org/doi/10.1103/PhysRevE.100.061302. 22

    work page doi:10.1103/physreve.100.061302 2019