In this study, we propose an explicit adaptive finite difference method (FDM) for the Cahn-Hilliard (CH) equation which describes the process of phase separation. The CH equation has been successfully utilized to model and simulate diverse field applications such as complex interfacial fluid flows and materials science. To numerically solve the CH equation fast and efficiently, we use the FDM and time-adaptive narrow-band domain. For the adaptive grid, we define a narrow-band domain including the interfacial transition layer of the phase field based on an undivided finite difference and solve the numerical scheme on the narrow-band domain. The proposed numerical scheme is based on an alternating direction explicit (ADE) method. To make the scheme conservative, we apply a mass correction algorithm after each temporal iteration step. To demonstrate the superior performance of the proposed adaptive FDM for the CH equation, we present two- and three-dimensional numerical experiments and compare them with those of other previous methods.
Keywords: Adaptive finite difference scheme; Cahn–Hilliard equation; Stable numerical method.
© The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature 2022, Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.