An efficient linear second order unconditionally stable direct discretization method for the phase-field crystal equation on surfaces

Abstract We develop an unconditionally stable direct discretization scheme for solving the phase-field crystal equation on surfaces. The surface is discretized by using an unstructured triangular mesh. Gradient, divergence, and Laplacian operators are defined on triangular meshes. The proposed numerical method is second-order accurate in space and time. At each time step, the proposed computational scheme results in linear elliptic equations to be solved, thus it is easy to implement the algorithm. It is proved that the proposed scheme satisfies a discrete energy-dissipation law. Therefore, it is unconditionally stable. A fast and efficient biconjugate gradients stabilized solver is used to solve the resulting discrete system. Numerical experiments are conducted to demonstrate the performance of the proposed algorithm.

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

[2]  Dong Li,et al.  On Second Order Semi-implicit Fourier Spectral Methods for 2D Cahn–Hilliard Equations , 2017, J. Sci. Comput..

[3]  A. Voigt,et al.  Capturing the complex physics behind universal grain size distributions in thin metallic films , 2013, 1304.4050.

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

[5]  Hongkai Zhao,et al.  A grid based particle method for solving partial differential equations on evolving surfaces and modeling high order geometrical motion , 2011, J. Comput. Phys..

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

[7]  Yibao Li,et al.  Computationally efficient adaptive time step method for the Cahn-Hilliard equation , 2017, Comput. Math. Appl..

[8]  Axel Voigt,et al.  The influence of membrane bound proteins on phase separation and coarsening in cell membranes. , 2012, Physical chemistry chemical physics : PCCP.

[9]  Karl Glasner,et al.  Improving the accuracy of convexity splitting methods for gradient flow equations , 2016, J. Comput. Phys..

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

[11]  M. Schanz,et al.  A high‐order FEM with exact geometry description for the Laplacian on implicitly defined surfaces , 2018 .

[12]  Martin Grant,et al.  Modeling elasticity in crystal growth. , 2001, Physical review letters.

[13]  Yibao Li,et al.  Direct Discretization Method for the Cahn–Hilliard Equation on an Evolving Surface , 2018, J. Sci. Comput..

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

[15]  Francisco Guillén-González,et al.  Second order schemes and time-step adaptivity for Allen-Cahn and Cahn-Hilliard models , 2014, Comput. Math. Appl..

[16]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

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

[18]  Charles M. Elliott,et al.  Evolving surface finite element method for the Cahn–Hilliard equation , 2013, Numerische Mathematik.

[19]  Yibao Li,et al.  An unconditionally energy-stable second-order time-accurate scheme for the Cahn-Hilliard equation on surfaces , 2017, Commun. Nonlinear Sci. Numer. Simul..

[20]  Charles M. Elliott,et al.  An Eulerian approach to transport and diffusion on evolving implicit surfaces , 2009, Comput. Vis. Sci..

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

[22]  Jaemin Shin,et al.  First and second order operator splitting methods for the phase field crystal equation , 2015, J. Comput. Phys..

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

[24]  Junseok Kim,et al.  Microphase separation patterns in diblock copolymers on curved surfaces using a nonlocal Cahn-Hilliard equation , 2015, The European physical journal. E, Soft matter.

[25]  Shimin Guo,et al.  Triply periodic minimal surface using a modified Allen-Cahn equation , 2017, Appl. Math. Comput..

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

[27]  Dynamic density functional theory of fluids , 2000 .

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

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

[30]  Badrinarayan P. Athreya,et al.  Using the phase-field crystal method in the multi-scale modeling of microstructure evolution , 2007 .