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

Loading...
Thumbnail Image

Access rights

openAccess
acceptedVersion

URL

Journal Title

Journal ISSN

Volume Title

A1 Alkuperäisartikkeli tieteellisessä aikakauslehdessä

Date

2016-09-15

Major/Subject

Mcode

Degree programme

Language

en

Pages

29

Series

Journal of Computational Physics, Volume 321, pp. 1026-1054

Abstract

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.

Description

Keywords

Amplitude equations, Energy stable scheme, Hexagonal finite differences, Multigrid, Phase field crystal

Other note

Citation

Guan, Z, Heinonen, V, Lowengrub, J, Wang, C & Wise, S M 2016, ' An energy stable, hexagonal finite difference scheme for the 2D phase field crystal amplitude equations ', Journal of Computational Physics, vol. 321, pp. 1026-1054 . https://doi.org/10.1016/j.jcp.2016.06.007