Premium
Numerical approximations for a phase field dendritic crystal growth model based on the invariant energy quadratization approach
Author(s) -
Zhao Jia,
Wang Qi,
Yang Xiaofeng
Publication year - 2016
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.5372
Subject(s) - discretization , mathematics , nonlinear system , compact finite difference , mathematical analysis , invariant (physics) , coefficient matrix , energy functional , physics , mathematical physics , quantum mechanics , eigenvalues and eigenvectors
Summary We present two accurate and efficient numerical schemes for a phase field dendritic crystal growth model, which is derived from the variation of a free‐energy functional, consisting of a temperature dependent bulk potential and a conformational entropy with a gradient‐dependent anisotropic coefficient. We introduce a novel Invariant Energy Quadratization approach to transform the free‐energy functional into a quadratic form by introducing new variables to substitute the nonlinear transformations. Based on the reformulated equivalent governing system, we develop a first and a second order semi‐discretized scheme in time for the system, in which all nonlinear terms are treated semi‐explicitly. The resulting semi‐discretized equations consist of a linear elliptic equation system at each time step, where the coefficient matrix operator is positive definite and thus, the semi‐discrete system can be solved efficiently. We further prove that the proposed schemes are unconditionally energy stable. Convergence test together with 2D and 3D numerical simulations for dendritic crystal growth are presented after the semi‐discrete schemes are fully discretized in space using the finite difference method to demonstrate the stability and the accuracy of the proposed schemes. Copyright © 2016 John Wiley & Sons, Ltd.