Energy stable compact scheme for Cahn-Hilliard equation with periodic boundary condition

Abstract We present a compact scheme to solve the Cahn–Hilliard equation with a periodic boundary condition, which is fourth-order accurate in space. We introduce schemes for two and three dimensions, which are derived from the one-dimensional compact stencil. The energy stability is completely proven for the proposed scheme based on the application of the compact method and well-known convex splitting methods. Detailed proofs of the mass conservation and unique solvability are also established. Numerical experiments are presented to demonstrate the accuracy and stability of the proposed methods.

[1]  Seunggyu Lee,et al.  Energy-minimizing wavelengths of equilibrium states for diblock copolymers in the hex-cylinder phase , 2015 .

[2]  Q. Du,et al.  A phase field approach in the numerical study of the elastic bending energy for vesicle membranes , 2004 .

[3]  Xiaofeng Yang,et al.  Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends , 2016, J. Comput. Phys..

[4]  Andrew M. Stuart,et al.  Model Problems in Numerical Stability Theory for Initial Value Problems , 1994, SIAM Rev..

[5]  J. Warren,et al.  Prediction of dendritic growth and microsegregation patterns in a binary alloy using the phase-field method , 1995 .

[6]  A. Karma,et al.  Phase-Field Simulation of Solidification , 2002 .

[7]  W. Weston Meyer,et al.  Optimal error bounds for cubic spline interpolation , 1976 .

[8]  Weizhong Dai,et al.  Fourth‐order compact schemes of a heat conduction problem with Neumann boundary conditions , 2007 .

[9]  Cheng Wang,et al.  Stable and efficient finite-difference nonlinear-multigrid schemes for the phase field crystal equation , 2009, J. Comput. Phys..

[10]  Jaemin Shin,et al.  Unconditionally stable methods for gradient flow using Convex Splitting Runge-Kutta scheme , 2017, J. Comput. Phys..

[11]  Charles M. Elliott,et al.  Numerical analysis of a model for phase separation of a multi- component alloy , 1996 .

[12]  John W. Cahn,et al.  On Spinodal Decomposition , 1961 .

[13]  Thomas J. R. Hughes,et al.  Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models , 2011, J. Comput. Phys..

[14]  Junseok Kim,et al.  A numerical method for the Cahn–Hilliard equation with a variable mobility , 2007 .

[15]  Andrea Lamorgese,et al.  Phase Field Approach to Multiphase Flow Modeling , 2011 .

[16]  Jaemin Shin,et al.  Numerical analysis of energy-minimizing wavelengths of equilibrium states for diblock copolymers , 2014 .

[17]  D. Jacqmin.  Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[18]  D. J. Eyre.  Unconditionally Gradient Stable Time Marching the Cahn-Hilliard Equation , 1998 .

[19]  Tingchun Wang,et al.  Fourth-order compact and energy conservative difference schemes for the nonlinear Schrödinger equation in two dimensions , 2013, J. Comput. Phys..

[20]  Shuguang Li,et al.  Numerical analysis for fourth-order compact conservative difference scheme to solve the 3D Rosenau-RLW equation , 2016, Comput. Math. Appl..

[21]  Jaemin Shin,et al.  A fourth-order spatial accurate and practically stable compact scheme for the Cahn-Hilliard equation , 2014 .

[22]  Jaemin Shin,et al.  Physical, mathematical, and numerical derivations of the Cahn–Hilliard equation , 2014 .

[23]  Yibao Li,et al.  A compact fourth-order finite difference scheme for the three-dimensional Cahn-Hilliard equation , 2016, Comput. Phys. Commun..

[24]  Jaemin Shin,et al.  First and second order numerical methods based on a new convex splitting for phase-field crystal equation , 2016, J. Comput. Phys..

[25]  A. Karma,et al.  Regular Article: Modeling Melt Convection in Phase-Field Simulations of Solidification , 1999 .

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