Linear, Second Order and Unconditionally Energy Stable Schemes for the Viscous Cahn-Hilliard Equation with Hyperbolic Relaxation using the Invariant Energy Quadratization Method
In this paper, we consider numerical approximations for the viscous Cahn-Hilliard equation with hyperbolic relaxation. This type of equations processes energy-dissipative structure. The main challenge in solving such a diffusive system numerically is how to develop high order temporal discretization for the hyperbolic and nonlinear terms, allowing large time-marching step, while preserving the energy stability, i.e. the energy dissipative structure at the time-discrete level. We resolve this issue by developing two second-order time-marching schemes using the recently developed "Invariant Energy Quadratization" approach where all nonlinear terms are discretized semi-explicitly. In each time step, one only needs to solve a symmetric positive definite (SPD) linear system. All the proposed schemes are rigorously proven to be unconditionally energy stable, and the second-order convergence in time has been verified by time step refinement tests numerically. Various 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy, and efficiency of the proposed schemes.
X. Yang et al., "Linear, Second Order and Unconditionally Energy Stable Schemes for the Viscous Cahn-Hilliard Equation with Hyperbolic Relaxation using the Invariant Energy Quadratization Method," Journal of Computational and Applied Mathematics, vol. 343, pp. 80 - 97, Elsevier B.V., Dec 2018.
The definitive version is available at https://doi.org/10.1016/j.cam.2018.04.027
Mathematics and Statistics
Center for High Performance Computing Research
Keywords and Phrases
Cahn-Hilliard; Flory-Huggins; Linear; Phase-field; Stability; Variable mobility
International Standard Serial Number (ISSN)
Article - Journal
© 2018 Elsevier B.V., All rights reserved.
01 Dec 2018