Document Type
Article
Journal/Book Title/Conference
SIAM Journal on Scientific Computing
Volume
42
Issue
1
Publisher
Society for Industrial and Applied Mathematics
Publication Date
1-21-2020
Award Number
NSF, Division of Mathematical Sciences 1816783
Funder
NSF, Division of Mathematical Sciences
First Page
B135
Last Page
B156
Abstract
We present a systematic approach to developing arbitrarily high-order, unconditionally energy stable numerical schemes for thermodynamically consistent gradient flow models that satisfy energy dissipation laws. Utilizing the energy quadratization method, we formulate the gradient flow model into an equivalent form with a corresponding quadratic free energy functional. Based on the equivalent form with a quadratic energy, we propose two classes of energy stable numerical approximations. In the first approach, we use a prediction-correction strategy to improve the accuracy of linear numerical schemes. In the second approach, we adopt the Gaussian collocation method to discretize the equivalent form with a quadratic energy, arriving at an arbitrarily high-order scheme for gradient flow models. Schemes derived using both approaches are proved rigorously to be unconditionally energy stable. The proposed schemes are then implemented in four gradient flow models numerically to demonstrate their accuracy and effectiveness. Detailed numerical comparisons among these schemes are carried out as well. These numerical strategies are rather general so that they can be readily generalized to solve any thermodynamically consistent PDE models.
Recommended Citation
Gong, Yuezheng, et al. “Arbitrarily High-Order Unconditionally Energy Stable Schemes for Thermodynamically Consistent Gradient Flow Models.” SIAM Journal on Scientific Computing, vol. 42, no. 1, Jan. 2020, pp. B135–56. DOI.org (Crossref), doi:10.1137/18M1213579.
