Convergence analysis and numerical implementation of a second order numerical scheme for the three-dimensional phase field crystal equation

https://doi.org/10.1016/j.camwa.2017.07.012Get rights and content
Under an Elsevier user license
open archive

Abstract

In this paper we analyze and implement a second-order-in-time numerical scheme for the three-dimensional phase field crystal (PFC) equation. The numerical scheme was proposed in Hu et al. (2009), with the unique solvability and unconditional energy stability established. However, its convergence analysis remains open. We present a detailed convergence analysis in this article, in which the maximum norm estimate of the numerical solution over grid points plays an essential role. Moreover, we outline the detailed multigrid method to solve the highly nonlinear numerical scheme over a cubic domain, and various three-dimensional numerical results are presented, including the numerical convergence test, complexity test of the multigrid solver and the polycrystal growth simulation.

Keywords

Three-dimensional phase field crystal
Finite difference
Energy stability
Second order numerical scheme
Convergence analysis
Nonlinear multigrid solver

Cited by (0)