An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation

Yibao Li, Junseok Kim

Research output: Contribution to journalArticle

9 Citations (Scopus)

Abstract

In this paper, we present a high-order accurate compact scheme for the phase field crystal model in two- and three-dimensional spaces. The proposed scheme is derived by combining a fourth-order compact finite difference formula in space and a backward differentiation for the time derivative term, which is second-order accurate in time. Furthermore, a nonlinearly stabilized splitting scheme is used and thus a larger time step can be allowed. Since the equations at the implicit time level are nonlinear, we introduce a Newton-type iterative method and employ a fast and efficient nonlinear multigrid solver to solve the resulting discrete system. In particular, we implement the compact scheme in the adaptive mesh refinement framework. An adaptive time step method for the phase field crystal model is also proposed. Various numerical experiments are presented and confirm the accuracy, stability, and efficiency of our proposed method.

Original language English 194-216 23 Computer Methods in Applied Mechanics and Engineering 319 https://doi.org/10.1016/j.cma.2017.02.022 Published - 2017 Jun 1

Fingerprint

crystal field theory
Crystals
Iterative methods
Derivatives
Experiments
newton

Keywords

• Adaptive mesh refinement
• Fourth-order compact scheme
• Phase-field crystal equation

ASJC Scopus subject areas

• Computational Mechanics
• Mechanics of Materials
• Mechanical Engineering
• Physics and Astronomy(all)
• Computer Science Applications

Cite this

In: Computer Methods in Applied Mechanics and Engineering, Vol. 319, 01.06.2017, p. 194-216.

Research output: Contribution to journalArticle

@article{636e08d2fc1143d3b911b8b042b7a613,
title = "An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation",
abstract = "In this paper, we present a high-order accurate compact scheme for the phase field crystal model in two- and three-dimensional spaces. The proposed scheme is derived by combining a fourth-order compact finite difference formula in space and a backward differentiation for the time derivative term, which is second-order accurate in time. Furthermore, a nonlinearly stabilized splitting scheme is used and thus a larger time step can be allowed. Since the equations at the implicit time level are nonlinear, we introduce a Newton-type iterative method and employ a fast and efficient nonlinear multigrid solver to solve the resulting discrete system. In particular, we implement the compact scheme in the adaptive mesh refinement framework. An adaptive time step method for the phase field crystal model is also proposed. Various numerical experiments are presented and confirm the accuracy, stability, and efficiency of our proposed method.",
keywords = "Adaptive mesh refinement, Adaptive time-stepping, Fourth-order compact scheme, Phase-field crystal equation",
author = "Yibao Li and Junseok Kim",
year = "2017",
month = "6",
day = "1",
doi = "10.1016/j.cma.2017.02.022",
language = "English",
volume = "319",
pages = "194--216",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0045-7825",
publisher = "Elsevier",

}

TY - JOUR

T1 - An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation

AU - Li, Yibao

AU - Kim, Junseok

PY - 2017/6/1

Y1 - 2017/6/1

N2 - In this paper, we present a high-order accurate compact scheme for the phase field crystal model in two- and three-dimensional spaces. The proposed scheme is derived by combining a fourth-order compact finite difference formula in space and a backward differentiation for the time derivative term, which is second-order accurate in time. Furthermore, a nonlinearly stabilized splitting scheme is used and thus a larger time step can be allowed. Since the equations at the implicit time level are nonlinear, we introduce a Newton-type iterative method and employ a fast and efficient nonlinear multigrid solver to solve the resulting discrete system. In particular, we implement the compact scheme in the adaptive mesh refinement framework. An adaptive time step method for the phase field crystal model is also proposed. Various numerical experiments are presented and confirm the accuracy, stability, and efficiency of our proposed method.

AB - In this paper, we present a high-order accurate compact scheme for the phase field crystal model in two- and three-dimensional spaces. The proposed scheme is derived by combining a fourth-order compact finite difference formula in space and a backward differentiation for the time derivative term, which is second-order accurate in time. Furthermore, a nonlinearly stabilized splitting scheme is used and thus a larger time step can be allowed. Since the equations at the implicit time level are nonlinear, we introduce a Newton-type iterative method and employ a fast and efficient nonlinear multigrid solver to solve the resulting discrete system. In particular, we implement the compact scheme in the adaptive mesh refinement framework. An adaptive time step method for the phase field crystal model is also proposed. Various numerical experiments are presented and confirm the accuracy, stability, and efficiency of our proposed method.

KW - Adaptive mesh refinement

KW - Adaptive time-stepping

KW - Fourth-order compact scheme

KW - Phase-field crystal equation

UR - http://www.scopus.com/inward/record.url?scp=85015394091&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85015394091&partnerID=8YFLogxK

U2 - 10.1016/j.cma.2017.02.022

DO - 10.1016/j.cma.2017.02.022

M3 - Article

VL - 319

SP - 194

EP - 216

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0045-7825

ER -