A general class of linear unconditionally energy stable schemes for the gradient flows

被引:7
作者
Tan, Zengqiang [1 ,2 ]
Tang, Huazhong [1 ,2 ,3 ]
机构
[1] Peking Univ, Ctr Appl Phys & Technol, HEDPS, Beijing 100871, Peoples R China
[2] Peking Univ, LMAM, Sch Math Sci, Beijing 100871, Peoples R China
[3] Nanchang Hangkong Univ, Nanchang 330000, Jiangxi, Peoples R China
基金
国家重点研发计划; 中国国家自然科学基金;
关键词
Gradient flows; Scalar auxiliary variable; General linear time discretization; Energy stability; Convergence; PHASE-FIELD MODELS; CONVEX SPLITTING SCHEMES; FOURIER-SPECTRAL METHODS; RUNGE-KUTTA METHODS; THIN-FILM MODEL; TUMOR-GROWTH; ALLEN-CAHN; NUMERICAL APPROXIMATIONS; EPITAXIAL-GROWTH; HIGHLY EFFICIENT;
D O I
10.1016/j.jcp.2022.111372
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
This paper studies a class of linear unconditionally energy stable schemes for the gradient flows. Such schemes are built on the SAV technique and the general linear time discretization (GLTD) as well as the linearization based on the extrapolation for the nonlinear term, and may be arbitrarily high-order accurate and very general, containing many existing SAV schemes and new SAV schemes. It is shown that the semi-discrete-in time schemes are unconditionally energy stable when the GLTD is algebraically stable, and are convergent with the order of min{q, nu} under the diagonal stability and some suitable regularity and accurate starting values, where q is the generalized stage order of the GLTD and nu denotes the number of the extrapolation points in time. The energy stability results can be easily extended to the fully discrete schemes, for example, if the Fourier spectral method is employed in space when the periodic boundary conditions are specified. Some numerical experiments on the Allen-Cahn, Cahn-Hilliard, and phase field crystal models are conducted to validate those theories as well as the effectiveness, the energy stability and the accuracy of our schemes. (C) 2022 Elsevier Inc. All rights reserved.
引用
收藏
页数:32
相关论文
共 101 条
[21]   Error Estimate of a Second Order Accurate Scalar Auxiliary Variable (SAV) Numerical Method for the Epitaxial Thin Film Equation [J].
Cheng, Qing ;
Wang, Cheng .
ADVANCES IN APPLIED MATHEMATICS AND MECHANICS, 2021, 13 (06) :1318-1354
[22]   Highly Efficient and Accurate Numerical Schemes for the Epitaxial Thin Film Growth Models by Using the SAV Approach [J].
Cheng, Qing ;
Shen, Jie ;
Yang, Xiaofeng .
JOURNAL OF SCIENTIFIC COMPUTING, 2019, 78 (03) :1467-1487
[23]   MULTIPLE SCALAR AUXILIARY VARIABLE (MSAV) APPROACH AND ITS APPLICATION TO THE PHASE-FIELD VESICLE MEMBRANE MODEL [J].
Cheng, Qing ;
Shen, Jie .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2018, 40 (06) :A3982-A4006
[24]   Search for highly stable two-step Runge-Kutta methods [J].
D'Ambrosio, R. ;
Izzo, G. ;
Jackiewicz, Z. .
APPLIED NUMERICAL MATHEMATICS, 2012, 62 (10) :1361-1379
[25]  
Dahlquist G., 1978, BIT (Nordisk Tidskrift for Informationsbehandling), V18, P384, DOI 10.1007/BF01932018
[26]  
Dahlquist G., 1976, Lect. Notes Math., P60, DOI 10.1007/BFb0080115
[27]   MAXIMUM PRINCIPLE PRESERVING EXPONENTIAL TIME DIFFERENCING SCHEMES FOR THE NONLOCAL ALLEN-CAHN EQUATION [J].
Du, Qiang ;
Ju, Lili ;
Li, Xiao ;
Qiao, Zhonghua .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2019, 57 (02) :875-898
[28]   Modeling elasticity in crystal growth [J].
Elder, KR ;
Katakowski, M ;
Haataja, M ;
Grant, M .
PHYSICAL REVIEW LETTERS, 2002, 88 (24) :2457011-2457014
[29]   THE GLOBAL DYNAMICS OF DISCRETE SEMILINEAR PARABOLIC EQUATIONS [J].
ELLIOTT, CM ;
STUART, AM .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1993, 30 (06) :1622-1663
[30]   Unconditionally gradient stable time marching the Cahn-Hilliard equation [J].
Eyre, DJ .
COMPUTATIONAL AND MATHEMATICAL MODELS OF MICROSTRUCTURAL EVOLUTION, 1998, 529 :39-46