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 条
[1]   ENERGY-DECAYING EXTRAPOLATED RK-SAV METHODS FOR THE ALLEN-CAHN AND CAHN-HILLIARD EQUATIONS [J].
Akrivis, Georgios ;
Li, Buyang ;
li, Dongfang .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2019, 41 (06) :A3703-A3727
[2]   MICROSCOPIC THEORY FOR ANTIPHASE BOUNDARY MOTION AND ITS APPLICATION TO ANTIPHASE DOMAIN COARSENING [J].
ALLEN, SM ;
CAHN, JW .
ACTA METALLURGICA, 1979, 27 (06) :1085-1095
[3]   Diffuse-interface methods in fluid mechanics [J].
Anderson, DM ;
McFadden, GB ;
Wheeler, AA .
ANNUAL REVIEW OF FLUID MECHANICS, 1998, 30 :139-165
[4]   Finite element approximation of nematic liquid crystal flows using a saddle-point structure [J].
Badia, Santiago ;
Guillen-Gonzalez, Francisco ;
Vicente Gutierrez-Santacreu, Juan .
JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 230 (04) :1686-1706
[5]   CONVERGENCE ANALYSIS OF A SECOND ORDER CONVEX SPLITTING SCHEME FOR THE MODIFIED PHASE FIELD CRYSTAL EQUATION [J].
Baskaran, A. ;
Lowengrub, J. S. ;
Wang, C. ;
Wise, S. M. .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2013, 51 (05) :2851-2873
[6]   Diffusive atomistic dynamics of edge dislocations in two dimensions [J].
Berry, J ;
Grant, M ;
Elder, KR .
PHYSICAL REVIEW E, 2006, 73 (03) :1-12
[7]   Phase-field simulation of solidification [J].
Boettinger, WJ ;
Warren, JA ;
Beckermann, C ;
Karma, A .
ANNUAL REVIEW OF MATERIALS RESEARCH, 2002, 32 :163-194
[8]   Adaptive phase-field computations of dendritic crystal growth [J].
Braun, RJ ;
Murray, BT .
JOURNAL OF CRYSTAL GROWTH, 1997, 174 (1-4) :41-53
[9]  
Burrage K., 1980, BIT (Nordisk Tidskrift for Informationsbehandling), V20, P185, DOI 10.1007/BF01933191