- Journal Home
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 7 (2010), pp. 362-382.
Published online: 2010-02
Cited by
- BibTex
- RIS
- TXT
The Cahn-Hilliard equation is often used to describe evolution of phase boundaries in phase field models for multiphase fluids. In this paper, we compare the use of the Cahn-Hilliard equation (of a constant mobility) for the phase variable with that of the singular or modified Cahn-Hilliard equation (of a variable mobility) in the context of physical derivation of the transport equation and numerical simulations of immiscible binary fluids. We show numerically that (i). both equations work fine for interfaces of small to moderate curvature in short to intermediate time scales; (ii) the Cahn-Hilliard equation renders strong dissipation in simulations of small droplets leading to dissolution of small droplets into the surrounding fluid and/or absorption of small droplets by larger droplets nearby, an artifact for immiscible binary fluids; whereas, the singular Cahn-Hilliard equation can significantly reduce the numerical dissipation around small droplets to yield physically acceptable results in intermediate time scales; (iii) the size of droplets that can be simulated by the Cahn-Hilliard equations scale inversely with the strength of the mixing free energy. Since the intermediate timescale is the time scale of interest in most transient fluid simulations, the singular Cahn-Hilliard equation proves to be the more accurate phase transporting equation for immiscible binary fluids.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.2009.09.016}, url = {http://global-sci.org/intro/article_detail/cicp/7634.html} }The Cahn-Hilliard equation is often used to describe evolution of phase boundaries in phase field models for multiphase fluids. In this paper, we compare the use of the Cahn-Hilliard equation (of a constant mobility) for the phase variable with that of the singular or modified Cahn-Hilliard equation (of a variable mobility) in the context of physical derivation of the transport equation and numerical simulations of immiscible binary fluids. We show numerically that (i). both equations work fine for interfaces of small to moderate curvature in short to intermediate time scales; (ii) the Cahn-Hilliard equation renders strong dissipation in simulations of small droplets leading to dissolution of small droplets into the surrounding fluid and/or absorption of small droplets by larger droplets nearby, an artifact for immiscible binary fluids; whereas, the singular Cahn-Hilliard equation can significantly reduce the numerical dissipation around small droplets to yield physically acceptable results in intermediate time scales; (iii) the size of droplets that can be simulated by the Cahn-Hilliard equations scale inversely with the strength of the mixing free energy. Since the intermediate timescale is the time scale of interest in most transient fluid simulations, the singular Cahn-Hilliard equation proves to be the more accurate phase transporting equation for immiscible binary fluids.