Efficient modified stabilized invariant energy quadratization approaches for phase-field crystal equation

The phase-field crystal equation is a sixth-order nonlinear parabolic equation and can be applied to simulate various phenomena such as epitaxial growth, material hardness, and phase transition. We propose a series of efficient modified stabilized invariant energy quadratization approaches with unconditional energy stability for the phase-field crystal model. Firstly, we propose a more suitable positive preserving function strictly in square root and consider a modified invariant energy quadratization (MIEQ) approach. Secondly, a series of efficient and suitable functionals H(ϕ) in square root are considered and the modified stabilized invariant energy quadratization (MSIEQ) approaches are developed. We prove the unconditional energy stability and optimal error estimates for the semi-discrete schemes carefully and rigorously. A comparative study of classical IEQ, MIEQ, and MSIEQ approaches is considered to show the accuracy and efficiency. Finally, we present various 2D numerical simulations to demonstrate the stability and accuracy.

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

[2]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[4]  P. Bates,et al.  International Journal of C 2009 Institute for Scientific Numerical Analysis and Modeling Computing and Information Numerical Analysis for a Nonlocal Allen-cahn Equation , 2022 .

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

[6]  Cheng Wang,et al.  An Energy-Stable and Convergent Finite-Difference Scheme for the Phase Field Crystal Equation , 2009, SIAM J. Numer. Anal..

[7]  Xiaofeng Yang,et al.  Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , 2010 .

[8]  Steven M. Wise,et al.  An Energy Stable and Convergent Finite-Difference Scheme for the Modified Phase Field Crystal Equation , 2011, SIAM J. Numer. Anal..

[9]  Jie Shen,et al.  Second-order Convex Splitting Schemes for Gradient Flows with Ehrlich-Schwoebel Type Energy: Application to Thin Film Epitaxy , 2012, SIAM J. Numer. Anal..

[10]  Zhengru Zhang,et al.  The stability and convergence of two linearized finite difference schemes for the nonlinear epitaxial growth model , 2012 .

[11]  Cheng Wang,et al.  A Linear Energy Stable Scheme for a Thin Film Model Without Slope Selection , 2012, J. Sci. Comput..

[12]  Sébastian Minjeaud.  An unconditionally stable uncoupled scheme for a triphasic Cahn–Hilliard/Navier–Stokes model , 2013 .

[13]  Steven M. Wise,et al.  Convergence Analysis of a Second Order Convex Splitting Scheme for the Modified Phase Field Crystal Equation , 2012, SIAM J. Numer. Anal..

[14]  Fengxin Chen,et al.  Convergence analysis for H1-Galerkin mixed finite element approximation of one nonlinear integro-differential model , 2013, Appl. Math. Comput..

[15]  Cheng Wang,et al.  A Linear Iteration Algorithm for a Second-Order Energy Stable Scheme for a Thin Film Model Without Slope Selection , 2014, J. Sci. Comput..

[16]  Zhaojie Zhou,et al.  Local discontinuous galerkin approximation of convection‐dominated diffusion optimal control problems with control constraints , 2014 .

[17]  A. Voigt,et al.  A Navier-Stokes phase-field crystal model for colloidal suspensions. , 2013, The Journal of chemical physics.

[18]  Fei Liu,et al.  Stabilized semi‐implicit spectral deferred correction methods for Allen–Cahn and Cahn–Hilliard equations , 2015 .

[19]  Ying Chen,et al.  Efficient, adaptive energy stable schemes for the incompressible Cahn-Hilliard Navier-Stokes phase-field models , 2016, J. Comput. Phys..

[20]  Junseok Kim,et al.  A simple and efficient finite difference method for the phase-field crystal equation on curved surfaces , 2016 .

[21]  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..

[22]  M. Grasselli,et al.  Energy stable and convergent finite element schemes for the modified phase field crystal equation , 2016 .

[23]  M. Dehghan,et al.  The numerical simulation of the phase field crystal (PFC) and modified phase field crystal (MPFC) models via global and local meshless methods , 2016 .

[24]  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..

[25]  Jia Zhao,et al.  Numerical approximations to a new phase field model for two phase flows of complex fluids , 2016 .

[26]  Zhaojie Zhou,et al.  Time-stepping discontinuous Galerkin approximation of optimal control problem governed by time fractional diffusion equation , 2018, Numerical Algorithms.

[27]  Zhifeng Weng,et al.  A Fourier spectral method for fractional-in-space Cahn–Hilliard equation , 2017 .

[28]  Yibao Li,et al.  An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation , 2017 .

[29]  Xiaofeng Yang,et al.  Numerical approximations of the Cahn-Hilliard and Allen-Cahn Equations with general nonlinear potential using the Invariant Energy Quadratization approach , 2017, 1712.02760.

[30]  Daozhi Han,et al.  Linearly first- and second-order, unconditionally energy stable schemes for the phase field crystal model , 2017, J. Comput. Phys..

[31]  Cheng Wang,et al.  A second‐order energy stable backward differentiation formula method for the epitaxial thin film equation with slope selection , 2017, 1706.01943.

[32]  Xiao Li,et al.  Stabilized linear semi-implicit schemes for the nonlocal Cahn-Hilliard equation , 2018, J. Comput. Phys..

[33]  Lili Ju,et al.  Unconditionally Energy Stable Linear Schemes for the Diffuse Interface Model with Peng–Robinson Equation of State , 2018, J. Sci. Comput..

[34]  Cheng Wang,et al.  A Second-Order Energy Stable BDF Numerical Scheme for the Cahn-Hilliard Equation , 2018 .

[35]  Ying Chen,et al.  A Uniquely Solvable, Energy Stable Numerical Scheme for the Functionalized Cahn–Hilliard Equation and Its Convergence Analysis , 2018, J. Sci. Comput..

[36]  Jiang Yang,et al.  The scalar auxiliary variable (SAV) approach for gradient flows , 2018, J. Comput. Phys..

[37]  Xiaoli Li,et al.  Efficient modified techniques of invariant energy quadratization approach for gradient flows , 2019, Appl. Math. Lett..

[38]  Wenqiang Feng,et al.  An energy stable fourth order finite difference scheme for the Cahn-Hilliard equation , 2017, J. Comput. Appl. Math..

[39]  Liquan Mei,et al.  Efficient numerical schemes with unconditional energy stabilities for the modified phase field crystal equation , 2019, Adv. Comput. Math..

[40]  Zhaojie Zhou,et al.  Spectral Galerkin approximation of optimal control problem governed by Riesz fractional differential equation , 2019, Applied Numerical Mathematics.