Year: 2012
Communications in Computational Physics, Vol. 11 (2012), Iss. 4 : pp. 1261–1278
Abstract
This paper studies the numerical simulations for the Cahn-Hilliard equation which describes a phase separation phenomenon. The numerical simulation of the Cahn-Hilliard model needs very long time to reach the steady state, and therefore large time-stepping methods become useful. The main objective of this work is to construct the unconditionally energy stable finite difference scheme so that the large time steps can be used in the numerical simulations. The equation is discretized by the central difference scheme in space and fully implicit second-order scheme in time. The proposed scheme is proved to be unconditionally energy stable and mass-conservative. An error estimate for the numerical solution is also obtained with second order in both space and time. By using this energy stable scheme, an adaptive time-stepping strategy is proposed, which selects time steps adaptively based on the variation of the free energy against time. The numerical experiments are presented to demonstrate the effectiveness of the adaptive time-stepping approach.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/cicp.300810.140411s
Communications in Computational Physics, Vol. 11 (2012), Iss. 4 : pp. 1261–1278
Published online: 2012-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 18