Journal
JOURNAL OF COMPUTATIONAL PHYSICS
Volume 419, Issue -, Pages -Publisher
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2020.109610
Keywords
Energy stable schemes; Gradient flow models; Runge-Kutta methods; Linear high-order schemes; Pseudo-spectral methods
Funding
- Foundation of Jiangsu Key Laboratory for Numerical Simulation of Large Scale Complex Systems [202002]
- Natural Science Foundation of Jiangsu Province [BK20180413]
- National Natural Science Foundation of China [11801269]
- NVIDIA Corporation
- SC EPSCoR/IDeA Program [DE-SC0020272]
- [NSF-DMS-1816783]
- [NSF-DMS-1815921]
- [OIA-1655740]
Ask authors/readers for more resources
We present a paradigm for developing arbitrarily high order, linear, unconditionally energy stable numerical algorithms for general gradient flow models. We apply the energy quadratization (EQ) technique to reformulate the gradient flow model into an equivalent one with a quadratic free energy and a modified mobility. For any positive integer k > 0 and t(n) = n Delta t, where Delta t is the time step size, we linearize the EQ-reformulated model in (t(n), t(n+1)] through an extrapolation from numerical solutions already obtained in [0,t(n)] so that the approximation error is in the order of O(Delta t(k)). Then we employ an s-stage algebraically stable Runge-Kutta method to discretize the linearized model in (t(n), t(n+1)]. For the spatial discretization, we use the Fourier pseudo-spectral method to match the order of accuracy in time. The resulting fully discrete scheme is linear, unconditionally energy stable, uniquely solvable, and can reach arbitrarily high order. Furthermore, we present a family of linear schemes based on prediction-correction methods to complement the new linear schemes. Some benchmark numerical examples are given to demonstrate the accuracy and efficiency of the schemes. (C) 2020 Elsevier Inc. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available