An energy stable, hexagonal finite difference scheme for the 2D phase field crystal amplitude equations

Tutkimustuotos: Lehtiartikkelivertaisarvioitu


  • Zhen Guan
  • Vili Heinonen
  • John Lowengrub
  • Cheng Wang
  • Steven M. Wise


  • University of California at Santa Barbara
  • University of Massachusetts Boston
  • University of Tennessee System


In this paper we construct an energy stable finite difference scheme for the amplitude expansion equations for the two-dimensional phase field crystal (PFC) model. The equations are formulated in a periodic hexagonal domain with respect to the reciprocal lattice vectors to achieve a provably unconditionally energy stable and solvable scheme. To our knowledge, this is the first such energy stable scheme for the PFC amplitude equations. The convexity of each part in the amplitude equations is analyzed, in both the semi-discrete and fully-discrete cases. Energy stability is based on a careful convexity analysis for the energy (in both the spatially continuous and discrete cases). As a result, unique solvability and unconditional energy stability are available for the resulting scheme. Moreover, we show that the scheme is point-wise stable for any time and space step sizes. An efficient multigrid solver is devised to solve the scheme, and a few numerical experiments are presented, including grain rotation and shrinkage and grain growth studies, as examples of the strength and robustness of the proposed scheme and solver.


JulkaisuJournal of Computational Physics
TilaJulkaistu - 15 syyskuuta 2016
OKM-julkaisutyyppiA1 Julkaistu artikkeli, soviteltu

ID: 6503893