Parallel energy-stable phase field crystal simulations based on domain decomposition methods

Abstract In this paper, we present a parallel numerical algorithm for solving the phase field crystal equation. In the algorithm, a semi-implicit finite difference scheme is derived based on the discrete variational derivative method. Theoretical analysis is provided to show that the scheme is unconditionally energy stable and can achieve second-order accuracy in both space and time. An adaptive time step strategy is adopted such that the time step size can be flexibly controlled based on the dynamical evolution of the problem. At each time step, a nonlinear algebraic system is constructed from the discretization of the phase field crystal equation and solved by a domain decomposition based, parallel Newton–Krylov–Schwarz method with improved boundary conditions for subdomain problems. Numerical experiments with several two and three dimensional test cases show that the proposed algorithm is second-order accurate in both space and time, energy stable with large time steps, and highly scalable to over ten thousands processor cores on the Sunway TaihuLight supercomputer.

[1]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[2]  J. Swift,et al.  Hydrodynamic fluctuations at the convective instability , 1977 .

[3]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[4]  T. Ala‐Nissila,et al.  Extended phase diagram of the three-dimensional phase field crystal model. , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[5]  Xiao-Chuan Cai,et al.  A Restricted Additive Schwarz Preconditioner for General Sparse Linear Systems , 1999, SIAM J. Sci. Comput..

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

[7]  Peter W Voorhees,et al.  Stress-induced morphological instabilities at the nanoscale examined using the phase field crystal approach , 2009 .

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

[9]  Olof B. Widlund,et al.  Domain Decomposition Algorithms with Small Overlap , 1992, SIAM J. Sci. Comput..

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

[11]  I. Gladwell,et al.  A Survey of Numerical Methods for Partial Differential Equations , 2021, An Introduction to Numerical Methods and Analysis.

[12]  Wei Ge,et al.  The Sunway TaihuLight supercomputer: system and applications , 2016, Science China Information Sciences.

[13]  Axel Voigt,et al.  Development and Analysis of a Block-Preconditioner for the Phase-Field Crystal Equation , 2015, SIAM J. Sci. Comput..

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

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

[16]  Chao Yang,et al.  A Scalable Implicit Solver for Phase Field Crystal Simulations , 2013, 2013 IEEE International Symposium on Parallel & Distributed Processing, Workshops and Phd Forum.

[17]  William Gropp,et al.  Newton-Krylov-Schwarz Methods in CFD , 1994 .

[18]  Yan Xu,et al.  Local Discontinuous Galerkin Method and High Order Semi-Implicit Scheme for the Phase Field Crystal Equation , 2016, SIAM J. Sci. Comput..

[19]  M. Grant,et al.  Phase-field crystal modeling and classical density functional theory of freezing , 2007 .

[20]  D. Furihata,et al.  Discrete Variational Derivative Method: A Structure-Preserving Numerical Method for Partial Differential Equations , 2010 .

[21]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

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

[23]  Hartmut Löwen,et al.  Derivation of the phase-field-crystal model for colloidal solidification. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[25]  Benedikt Wirth,et al.  A simple and efficient scheme for phase field crystal simulation , 2013 .

[26]  Axel Voigt,et al.  Nucleation and growth by a phase field crystal (PFC) model , 2007 .

[27]  Nikolas Provatas,et al.  Phase-Field Methods in Materials Science and Engineering , 2010 .

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

[29]  Xiao-Chuan Cai,et al.  Restricted Additive Schwarz Preconditioners with Harmonic Overlap for Symmetric Positive Definite Linear Systems , 2003, SIAM J. Numer. Anal..

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

[31]  Marcel Bauer,et al.  Numerical Methods for Partial Differential Equations , 1994 .