A HIGHLY EFFICIENT AND ACCURATE NEW SCALAR AUXILIARY VARIABLE APPROACH FOR GRADIENT FLOWS

被引:105
作者
Huang, Fukeng [1 ]
Shen, Jie [1 ]
Yang, Zhiguo [1 ]
机构
[1] Purdue Univ, Dept Math, W Lafayette, IN 47907 USA
关键词
gradient flow; SAV approach; adaptive time-stepping; high-order methods; stability; ENERGY-STABLE SCHEMES; FINITE-ELEMENT APPROXIMATION; MULTICOMPONENT ALLOY; NONUNIFORM SYSTEM; NUMERICAL SCHEMES; PHASE-SEPARATION; SAV APPROACH; CAHN; 2ND-ORDER; MODEL;
D O I
10.1137/19M1298627
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We present several essential improvements to the powerful scalar auxiliary variable (SAV) approach. Firstly, by using the introduced scalar variable to control both the nonlinear and the explicit linear terms, we are able to reduce the number of linear equations with constant coefficients to be solved at each time step from two to one, so the computational cost of the new SAV approach is essentially half of the original SAV approach while keeping all its other advantages. This technique is also extended to the multiple SAV approach. Secondly, instead of discretizing the dynamical equation for the auxiliary variable, we use a first-order approximation of the energy balance equation, which allows us to construct high-order unconditionally energy-stable SAV schemes with uniform and, more importantly, variable time step sizes, enabling us to construct, for the first time, high-order unconditionally stable adaptive time-stepping backward differentiation formula schemes. Representative numerical examples are provided to demonstrate the improved efficiency and accuracy of the proposed method.
引用
收藏
页码:A2514 / A2536
页数:23
相关论文
共 34 条
[1]   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
[2]   OBSERVATION OF BOSE-EINSTEIN CONDENSATION IN A DILUTE ATOMIC VAPOR [J].
ANDERSON, MH ;
ENSHER, JR ;
MATTHEWS, MR ;
WIEMAN, CE ;
CORNELL, EA .
SCIENCE, 1995, 269 (5221) :198-201
[3]  
[Anonymous], 2011, Spectral methods, DOI [DOI 10.1007/978-3-540-71041-7, 10.1007/978-3-540-71041-7]
[4]   Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow [J].
Bao, WZ ;
Du, Q .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2004, 25 (05) :1674-1697
[5]   Finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy [J].
Barrett, JW ;
Blowey, JF .
NUMERISCHE MATHEMATIK, 1997, 77 (01) :1-34
[6]   Finite element approximation of a model for phase separation of a multi-component alloy with a concentration-dependent mobility matrix [J].
Barrett, JW ;
Blowey, JF .
IMA JOURNAL OF NUMERICAL ANALYSIS, 1998, 18 (02) :287-328
[7]   FREE ENERGY OF A NONUNIFORM SYSTEM .1. INTERFACIAL FREE ENERGY [J].
CAHN, JW ;
HILLIARD, JE .
JOURNAL OF CHEMICAL PHYSICS, 1958, 28 (02) :258-267
[8]   FREE ENERGY OF A NONUNIFORM SYSTEM .3. NUCLEATION IN A 2-COMPONENT INCOMPRESSIBLE FLUID [J].
CAHN, JW ;
HILLIARD, JE .
JOURNAL OF CHEMICAL PHYSICS, 1959, 31 (03) :688-699
[9]   Preserving energy resp dissipation in numerical PDEs using the "Average Vector Field" method [J].
Celledoni, E. ;
Grimm, V. ;
McLachlan, R. I. ;
McLaren, D. I. ;
O'Neale, D. ;
Owren, B. ;
Quispel, G. R. W. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 231 (20) :6770-6789
[10]   A SECOND ORDER BDF NUMERICAL SCHEME WITH VARIABLE STEPS FOR THE CAHN-HILLIARD EQUATION [J].
Chen, Wenbin ;
Wang, Xiaoming ;
Yan, Yue ;
Zhang, Zhuying .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2019, 57 (01) :495-525