- Open Access
Droplet shrinkage in phase-field approaches: A comparison of the Allen-Cahn and the Cahn-Hilliard models
Phys. Rev. E 113, 055107 – Published 18 May, 2026
DOI: https://doi.org/10.1103/y6sf-3lf2
Abstract
Diffuse interface methods for multiphase flow simulations often exhibit nonphysical droplet or bubble shrinkage, particularly when based on the Cahn-Hilliard equation. This well-known artifact introduces a critical radius below which droplets vanish, thereby limiting the fidelity of simulations involving small-scale structures. In this study we examine a conservative Allen-Cahn model and demonstrate, through both analytical and numerical investigations, that it inherently avoids this shrinkage behavior. We show that the conservative Allen-Cahn model enforces an exact local balance of interfacial terms, enabling the stable preservation of droplets with radii larger than the interface thickness, regardless of the initial droplet size or domain geometry. Our analysis also provides a different theoretical treatment of the shrinkage phenomenon in the Cahn-Hilliard equation, distinguishing itself from earlier works by providing a rigorous explanation of why conservative Allen-Cahn models avoid this problem. The analysis reveals that conservative Allen-Cahn formulations achieve exact cancellation of curvature-driven terms through a local geometric force balance, while Cahn-Hilliard models exhibit systematic shrinkage due to unbalanced curvature effects that create driving forces for mass transport. Theoretical and numerical findings show that there is a shrinkage relation between interface radius and time in the form of , which leads to a power law for the disappearance time versus initial radius . Numerical experiments across a wide range of droplet sizes corroborate the theoretical predictions. These results position the conservative Allen-Cahn model as a robust and accurate phase-field approach, particularly well suited for applications demanding the resolution of fine features or long-time integration, such as multiphase transport in porous media.
Physics Subject Headings (PhySH)
Article Text
References (30)
- D. M. Anderson, G. B. McFadden, and A. A. Wheeler, Diffuse-interface methods in fluid mechanics, Annu. Rev. Fluid Mech. 30, 139 (1998).
- D. Jacqmin, Calculation of two-phase Navier-Stokes flows using phase-field modeling, J. Comput. Phys. 155, 96 (1999).
- J. W. Cahn and J. E. Hilliard, Free energy of a nonuniform system. I. Interfacial free energy, J. Chem. Phys. 28, 258 (1958).
- S. M. Allen and J. W. Cahn, A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening, Acta Metall. 27, 1085 (1979).
- P. Yue, C. Zhou, and J. J. Feng, Spontaneous shrinkage of drops and mass conservation in phase-field simulations, J. Comput. Phys. 223, 1 (2007).
- L. Zheng, T. Lee, Z. Guo, and D. Rumschitzki, Shrinkage of bubbles and drops in the lattice Boltzmann equation method for nonideal gases, Phys. Rev. E 89, 033302 (2014).
- J. E. McClure, Z. Li, M. Berrill, and T. Ramstad, The LBPM software package for simulating multiphase flow on digital images of porous rocks, Comput. Geosci. 25, 871 (2021).
- Y. Wang, C. Shu, J. Shao, J. Wu, and X. Niu, A mass-conserved diffuse interface method and its application for incompressible multiphase flows with large density ratio, J. Comput. Phys. 290, 336 (2015).
- Y. Hu, Q. He, D. Li, Y. Li, and X. Niu, On the total mass conservation and the volume preservation in the diffuse interface method, Comput. Fluids 193, 104291 (2019).
- P. Yue, C. Zhou, J. J. Feng, C. F. Ollivier-Gooch, and H. H. Hu, Phase-field simulations of interfacial dynamics in viscoelastic fluids using finite elements with adaptive meshing, J. Comput. Phys. 219, 47 (2006).
- Y. Wang, C. Shu, L. Yang, and H. Yuan, On the re-initialization of fluid interfaces in diffuse interface method, Comput. Fluids 166, 209 (2018).
- Y. Sun and C. Beckermann, Sharp interface tracking using the phase-field equation, J. Comput. Phys. 220, 626 (2007).
- P.-H. Chiu and Y.-T. Lin, A conservative phase field method for solving incompressible two-phase flows, J. Comput. Phys. 230, 185 (2011).
- M. Geier, A. Fakhari, and T. Lee, Conservative phase-field lattice Boltzmann model for interface tracking equation, Phys. Rev. E 91, 063309 (2015).
- R. Haghani, M. H. Rahimian, and A. Fakhari, Conservative phase-field lattice-Boltzmann model for ternary fluids, J. Comput. Phys. 374, 668 (2018).
- R. Haghani, A. Fakhari, and M.-H. Rahimian, Phase-change modeling based on a novel conservative phase-field method, J. Comput. Phys. 432, 110111 (2021).
- R. Haghani, H. Erfani, J. E. McClure, E. G. Flekkøy, and C. F. Berg, Color-gradient-based phase-field equation for multiphase flow, Phys. Rev. E 109, 035301 (2024).
- S. Bakhshian, S. A. Hosseini, and N. Shokri, Pore-scale characteristics of multiphase flow in heterogeneous porous media using the lattice Boltzmann method, Sci. Rep. 9, 3377 (2019).
- A. K. Gunstensen, D. H. Rothman, S. Zaleski, and G. Zanetti, Lattice Boltzmann model of immiscible fluids, Phys. Rev. A 43, 4320 (1991).
- A. Begmohammadi, R. Haghani-Hassan-Abadi, A. Fakhari, and D. Bolster, Study of phase-field lattice Boltzmann models based on the conservative Allen-Cahn equation, Phys. Rev. E 102, 023305 (2020).
- I. Steinbach, Phase-field models in materials science, Modell. Simul. Mater. Sci. Eng. 17, 073001 (2009).
- Y. Li, Q. Yu, S. Ham, S. Kwak, C. Lee, and J. Kim, A phase-field model without artificial curvature effect for the crystal growth simulation, Int. J. Heat Mass Transf. 203, 123847 (2023).
- F. Capuani, I. Pagonabarraga, and D. Frenkel, Discrete solution of the electrokinetic equations, J. Chem. Phys. 121, 973 (2004).
- R. Haghani and M.-H. Rahimian, Hybrid lattice Boltzmann finite difference model for simulation of phase change in a ternary fluid, Int. J. Heat Mass Transf. 127, 704 (2018).
- T. Lee, Effects of incompressibility on the elimination of parasitic currents in the lattice Boltzmann equation method for binary fluids, Comput. Math. Appl. 58, 987 (2009).
- R. Haghani, H. Erfani, J. E. McClure, and C. F. Berg, A note on the summation relation in phase-field equations, Phys. Fluids 35, 092112 (2023).
- J. Bao and Z. Guo, Phase-field lattice Boltzmann model with singular mobility for quasi-incompressible two-phase flows, Phys. Rev. E 109, 025302 (2024).
- J. Kim, A numerical method for the Cahn-Hilliard equation with a variable mobility, Commun. Nonlinear Sci. Numer. Simul. 12, 1560 (2007).
- T. Akai, B. Bijeljic, and M. J. Blunt, Wetting boundary condition for the color-gradient lattice Boltzmann method: Validation with analytical and experimental data, Adv. Water Resour. 116, 56 (2018).
- R. Haghani, GitHub (2026), https://github.com/haghanire/droplet_shrinkage/.