Communications on Applied Mathematics and Computation ›› 2025, Vol. 7 ›› Issue (4): 1378-1397.doi: 10.1007/s42967-023-00326-1

• ORIGINAL PAPERS • 上一篇    下一篇

The L2-1σ/LDG Method for the Caputo Diffusion Equation with a Variable Coefficient

Qiaoqiao Dai, Dongxia Li   

  1. Department of Mathematics, Shanghai University, Shanghai, 200444, China
  • 收稿日期:2023-08-03 修回日期:2023-08-22 接受日期:2023-09-19 出版日期:2023-12-06 发布日期:2023-12-06
  • 通讯作者: Dongxia Li,E-mail:lidongxia96@shu.edu.cn E-mail:lidongxia96@shu.edu.cn

The L2-1σ/LDG Method for the Caputo Diffusion Equation with a Variable Coefficient

Qiaoqiao Dai, Dongxia Li   

  1. Department of Mathematics, Shanghai University, Shanghai, 200444, China
  • Received:2023-08-03 Revised:2023-08-22 Accepted:2023-09-19 Online:2023-12-06 Published:2023-12-06

摘要: In this paper, an efficient method is proposed to solve the Caputo diffusion equation with a variable coefficient. Since the solution of such an equation in general has a typical weak singularity near the initial time t=0, the time-fractional derivative with order in (0, 1) is discretized by L2-1σ formula on nonuniform meshes. For the spatial derivative, the local discontinuous Galerkin (LDG) method is employed. A complete theoretical analysis of the numerical stability and convergence of the derived scheme is given using a discrete fractional Gronwall inequality. Numerical experiments demonstrate the validity of the established scheme and the accuracy of the theoretical analysis results.

关键词: Local discontinuous Galerkin (LDG) method, Nonuniform meshes, L2-1σ method, Stability analysis, Error estimate, Variable coefficient

Abstract: In this paper, an efficient method is proposed to solve the Caputo diffusion equation with a variable coefficient. Since the solution of such an equation in general has a typical weak singularity near the initial time t=0, the time-fractional derivative with order in (0, 1) is discretized by L2-1σ formula on nonuniform meshes. For the spatial derivative, the local discontinuous Galerkin (LDG) method is employed. A complete theoretical analysis of the numerical stability and convergence of the derived scheme is given using a discrete fractional Gronwall inequality. Numerical experiments demonstrate the validity of the established scheme and the accuracy of the theoretical analysis results.

Key words: Local discontinuous Galerkin (LDG) method, Nonuniform meshes, L2-1σ method, Stability analysis, Error estimate, Variable coefficient

中图分类号: