An Adaptive Time-Stepping Strategy for the Cahn-Hilliard Equation

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.

[1]  C. M. Elliott,et al.  Numerical Studies of the Cahn-Hilliard Equation for Phase Separation , 1987 .

[2]  M. Gurtin,et al.  Structured phase transitions on a finite interval , 1984 .

[3]  R. Nicolaides,et al.  Numerical analysis of a continuum model of phase transition , 1991 .

[4]  Tao Tang,et al.  An Adaptive Time-Stepping Strategy for the Molecular Beam Epitaxy Models , 2011, SIAM J. Sci. Comput..

[5]  Charles M. Elliott,et al.  Error estimates with smooth and nonsmooth data for a finite element method for the Cahn-Hilliard equation , 1992 .

[6]  Jie Shen,et al.  A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method , 2003 .

[7]  Gustaf Söderlind,et al.  Automatic Control and Adaptive Time-Stepping , 2002, Numerical Algorithms.

[8]  Susan E. Minkoff,et al.  A comparison of adaptive time stepping methods for coupled flow and deformation modeling , 2006 .

[9]  J. Warren,et al.  Controlling the accuracy of unconditionally stable algorithms in the Cahn-Hilliard equation. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Yunqing Huang,et al.  Moving mesh methods with locally varying time steps , 2004 .

[11]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[12]  S. M. Choo,et al.  Conservative nonlinear difference scheme for the Cahn-Hilliard equation—II , 1998 .

[13]  C. M. Elliott,et al.  A nonconforming finite-element method for the two-dimensional Cahn-Hilliard equation , 1989 .

[14]  L. Segel,et al.  Nonlinear aspects of the Cahn-Hilliard equation , 1984 .

[15]  Ming Wang,et al.  A nonconforming finite element method for the Cahn-Hilliard equation , 2010, J. Comput. Phys..

[16]  Stig Larsson,et al.  THE CAHN-HILLIARD EQUATION , 2007 .

[17]  Tao Tang,et al.  Stability Analysis of Large Time-Stepping Methods for Epitaxial Growth Models , 2006, SIAM J. Numer. Anal..

[18]  Zhi-zhong Sun,et al.  A second-order accurate linearized difference scheme for the two-dimensional Cahn-Hilliard equation , 1995 .

[19]  Daisuke Furihata,et al.  A stable and conservative finite difference scheme for the Cahn-Hilliard equation , 2001, Numerische Mathematik.

[20]  M ChooS,et al.  Cahn‐Hilliad方程式に関する保存型非線形差分スキーム‐II , 2000 .

[21]  Yu-Lin Chou.  Applications of Discrete Functional Analysis to the Finite Difference Method , 1991 .

[22]  S. M. Wise,et al.  Unconditionally Stable Finite Difference, Nonlinear Multigrid Simulation of the Cahn-Hilliard-Hele-Shaw System of Equations , 2010, J. Sci. Comput..

[23]  Yinnian He,et al.  On large time-stepping methods for the Cahn--Hilliard equation , 2007 .