Second-Order Semi-Lagrangian Exponential Time Differencing Method with Enhanced Error Estimate for the Convective Allen-Cahn Equation

被引:5
作者
Li, Jingwei [1 ]
Lan, Rihui [2 ]
Cai, Yongyong [3 ,4 ]
Ju, Lili [5 ]
Wang, Xiaoqiang [6 ]
机构
[1] Lanzhou Univ, Sch Math & Stat, Lanzhou 730000, Peoples R China
[2] Ocean Univ China, Sch Math Sci, Qingdao 266100, Shandong, Peoples R China
[3] Beijing Normal Univ, Lab Math & Complex Syst, Beijing 100875, Peoples R China
[4] Beijing Normal Univ, Sch Math Sci, Beijing 100875, Peoples R China
[5] Univ South Carolina, Dept Math, Columbia, SC 29208 USA
[6] Florida State Univ, Dept Sci Comp, Tallahassee, FL 32306 USA
基金
美国国家科学基金会; 中国国家自然科学基金;
关键词
Convective Allen-Cahn equation; Semi-Lagrangian method; Variable coefficients; Maximum bound principle; Exponential time differencing; Enhanced error estimate; FINITE-ELEMENT-METHOD; MAXIMUM PRINCIPLE; SCHEMES; CONVERGENCE; STABILITY;
D O I
10.1007/s10915-023-02316-0
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The convective Allen-Cahn (CAC) equation has been widely used for simulating multiphase flows of incompressible fluids, which contains an extra convective term but still maintains the same maximum bound principle (MBP) as the classic Allen-Cahn equation. Based on the operator splitting approach, we propose a second-order semi-Lagrangian exponential time differencing method for solving the CAC equation, that preserves the discrete MBP unconditionally. In our scheme, the AC equation part is first spatially discretized via the central finite difference scheme, then it is efficiently solved by using the exponential time differencing method with FFT-based fast implementation. The transport equation part is computed by combining the semi-Lagrangian approach with a cut-off post-processing within the finite difference framework. MBP stability and convergence analysis of our fully discretized scheme are presented. In particular, we conduct an improved error estimation for the semi-Lagrangian method with variable velocity, so that the error of our scheme is not spoiled by the reciprocal of the time step size. Extensive numerical tests in two and three dimensions are also carried out to validate the theoretical results and demonstrate the performance of our scheme.
引用
收藏
页数:29
相关论文
共 61 条
[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]  
Allievi A, 2000, INT J NUMER METH FL, V32, P439, DOI 10.1002/(SICI)1097-0363(20000229)32:4<439::AID-FLD946>3.0.CO
[3]  
2-Y
[4]  
Benque J.P., 1982, 4 INT S FINITE ELEME, P295
[5]  
Bercovier M., 1982, FINITE ELEMENT FLOW, P67
[6]   Lagrange-Galerkin methods for the incompressible Navier-Stokes equations: a review [J].
Bermejo, Rodolfo ;
Saavedra, Laura .
COMMUNICATIONS IN APPLIED AND INDUSTRIAL MATHEMATICS, 2016, 7 (03) :23-52
[7]   A new class of time discretization schemes for the solution of nonlinear PDEs [J].
Beylkin, G ;
Keiser, JM ;
Vozovoi, L .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 147 (02) :362-387
[8]  
Boukir K, 1997, INT J NUMER METH FL, V25, P1421, DOI 10.1002/(SICI)1097-0363(19971230)25:12<1421::AID-FLD334>3.0.CO
[9]  
2-A
[10]  
Cai YY, 2023, COMMUN MATH SCI, V21, P127