Improving the accuracy of convexity splitting methods for gradient flow equations

This paper introduces numerical time discretization methods which significantly improve the accuracy of the convexity-splitting approach of Eyre (1998) 7, while retaining the same numerical cost and stability properties.A first order method is constructed by iteration of a semi-implicit method based upon decomposing the energy into convex and concave parts. A second order method is also presented based on backwards differentiation formulas. Several extrapolation procedures for iteration initialization are proposed. We show that, under broad circumstances, these methods have an energy decreasing property, leading to good numerical stability.The new schemes are tested using two evolution equations commonly used in materials science: the Cahn-Hilliard equation and the phase field crystal equation. We find that our methods can increase accuracy by many orders of magnitude in comparison to the original convexity-splitting algorithm. In addition, the optimal methods require little or no iteration, making their computation cost similar to the original algorithm.

[1]  Xiangrong Li,et al.  Nonlinear simulations of solid tumor growth using a mixture model: invasion and branching , 2009, Journal of mathematical biology.

[2]  James A. Warren,et al.  An efficient algorithm for solving the phase field crystal model , 2008, J. Comput. Phys..

[3]  M. Grant,et al.  Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  Long-Qing Chen.  Phase-Field Models for Microstructure Evolution , 2002 .

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

[6]  J. Cahn,et al.  A microscopic theory for antiphase boundary motion and its application to antiphase domain coasening , 1979 .

[7]  K. Promislow,et al.  On the unconditionally gradient stable scheme for the Cahn-Hilliard equation and its implementation with Fourier method , 2013 .

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

[9]  Tony F. Chan,et al.  Image processing and analysis - variational, PDE, wavelet, and stochastic methods , 2005 .

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

[11]  Jie Shen,et al.  Applications of semi-implicit Fourier-spectral method to phase field equations , 1998 .

[12]  Yuan Ma,et al.  An adaptive time-stepping strategy for solving the phase field crystal model , 2013, J. Comput. Phys..

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

[14]  Victor M. Calo,et al.  An energy-stable convex splitting for the phase-field crystal equation , 2014, 1405.3488.

[15]  Keith Promislow,et al.  High accuracy solutions to energy gradient flows from material science models , 2014, J. Comput. Phys..

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

[17]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .

[18]  Xiaoming Wang,et al.  A second order in time, uniquely solvable, unconditionally stable numerical scheme for Cahn-Hilliard-Navier-Stokes equation , 2014, J. Comput. Phys..

[19]  Chert,et al.  Applications of semi-implicit Fourier-spectral method to phase field equations , 2004 .

[20]  Franck Boyer,et al.  Numerical schemes for a three component Cahn-Hilliard model , 2011 .

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

[22]  Xesús Nogueira,et al.  An unconditionally energy-stable method for the phase field crystal equation , 2012 .

[23]  Karl B Glasner,et al.  Grain boundary motion arising from the gradient flow of the Aviles–Giga functional , 2006 .

[24]  Guillermo Sapiro,et al.  Fourth order partial differential equations on general geometries , 2006, J. Comput. Phys..

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