非均匀网格上时间分数阶扩散—波动方程的 BDF2 型有限元方法

数学物理学报 ›› 2025, Vol. 45 ›› Issue (4): 1268-1290.

非均匀网格上时间分数阶扩散—波动方程的 BDF2 型有限元方法

祝鹏¹(),陈艳萍^2,^*(),徐先宇³()

¹嘉兴大学数据科学学院浙江嘉兴 314001
²南京邮电大学理学院南京 210023
³湘潭大学数学与计算科学学院湖南湘潭 411100

收稿日期:2024-09-19 修回日期:2025-01-30 出版日期:2025-08-26 发布日期:2025-08-01
通讯作者: *E-mail: ypchen@njupt.edu.cn
作者简介:E-mail: pzh@zjxu.edu.cn;|xuxianyuqdu@163.com
基金资助:
国家自然科学基金天元数学访问学者项目(12426616);浙江省自然科学基金(LY23A010005);南京邮电大学引进人才科研启动基金(NY223127)

BDF2-Type Finite Element Method for Time-Fractional Diffusion-Wave Equations on Nonuniform Grids

Zhu Peng¹(),Chen Yanping^2,^*(),Xu Xianyu³()

¹School of Data Science, Jiaxing University, Zhejiang Jiaxing 314001
²School of Science, Nanjing University of Posts and Telecommunications, Nanjing 210023
³School of Mathematics and Computational Science, Xiangtan University, Hunan Xiangtan 411100

Received:2024-09-19 Revised:2025-01-30 Online:2025-08-26 Published:2025-08-01
Supported by:
Visiting scholar program of National Natural Science Foundation of China(12426616);Zhejiang Provincial Natural Science Foundation of China(LY23A010005);Natural Science Research Start-up Foundation of Recruiting Talents of Nanjing University of Posts and Telecommunications(NY223127)

摘要/Abstract

摘要：

众所周知, 非均匀网格的研究可以有效地解决分数阶 Caputo 型导数的初值奇异现象. 在非均匀网格的理论分析中, 经常采用分数阶离散 Grönwall 不等式进行误差分析, 缺乏对误差结构的具体研究. 设计了一种非均匀网格上的误差卷积结构, 用于分析时间分数阶扩散-波动方程. 将二次插值近似应用于 Caputo 型导数, 通过使用降阶法和离散互补卷积核对 Caputo 型导数进行离散, 得到了非均匀网格上的 BDF2 型有限元方法. 离散互补卷积核在算法的收敛性分析中至关重要, 因为它简化有限元理论分析的过程, 并基于卷积核和插值估计的性质构建了全局一致性误差估计. 详细估计了非均匀网格上 BDF2 有限元格式的 $L^2$-范数误差和 $H^1$-范数误差, 并通过实验验证了所提出的有限元格式与理论收敛阶之间的一致性.

关键词: 时间分数阶扩散-波动方程, 离散卷积核, BDF2 型有限元格式, 误差卷积结构, 非均匀网格

Abstract:

As is well known, the study of nonuniform grids can effectively solve the initial value singularity phenomenon of fractional Caputo -type derivatives. In the theoretical analysis of nonuniform grids, fractional discrete Grönwall inequality is often used for error analysis, but there is a lack of specific research on error structures. An error convolution structure (ECS) was designed on nonuniform grids for analyzing the time fractional diffusion wave equation. A quadratic interpolation approximation was applied Caputo -type derivatives, and the BDF2 -type finite element method on nonuniform grids was obtained by discretizing it using a reduction method and a discrete complementary convolution kernel. The discrete complementary convolution kernel is crucial in the convergence analysis of algorithms, as it simplify the process of finite element theory analysis and construct global consistency errors based on the properties of convolution kernels and interpolation estimates. The $L^2$-norm error and $H^1$-norm error of the BDF2 finite element scheme on nonuniform grids were estimated in detail, and verifies the consistency between the proposed finite element scheme and the theoretical convergence order through experiments.

Key words: time fractional diffusion-wave equation, discrete complementary convolution kernels, BDF2 finite element format, error convolutional structure, nonuniform grids

中图分类号:

O241.8

祝鹏, 陈艳萍, 徐先宇. 非均匀网格上时间分数阶扩散—波动方程的 BDF2 型有限元方法[J]. 数学物理学报, 2025, 45(4): 1268-1290.

Zhu Peng, Chen Yanping, Xu Xianyu. BDF2-Type Finite Element Method for Time-Fractional Diffusion-Wave Equations on Nonuniform Grids[J]. Acta mathematica scientia,Series A, 2025, 45(4): 1268-1290.

图/表 4

参考文献 23

[1]	Zheng X, Wang H. Variable-order space-fractional diffusion equations and a variable-order modification of constant-order fractional problems. Applicable Analysis, 2022, 1.1(6): 1848-1870
[2]	Mukherjee M, Mondal B. An integer-order SIS epidemic model having variable population and fear effect: Comparing the stability with fractional order. Journal of the Egyptian Mathematical Society, 2022, 30(1): 1-23
[3]	Chakraverty S, Jena R M, Jena S K. Computational Fractional Dynamical Systems:Fractional Differential Equations and Applications. New York: John Wiley and Sons Inc, 2022
[4]	Pandey P, Das S, Craciun E M, et al. Two-dimensional nonlinear time fractional reaction-diffusion equation in application to sub-diffusion process of the multicomponent fluid in porous media. Meccanica, 2021, 56(1): 99-115
[5]	Acioli P S, Xavier F A, Moreira D M. Mathematical model using fractional derivatives applied to the dispersion of pollutants in the planetary boundary layer. Boundary-Layer Meteorology, 2019, 1.0(2): 285-304
[6]	徐斐, 张勇. 分数阶 Burgers 方程时间周期弱解的唯一性与渐近稳定性. 数学物理学报, 2023, 43A(6): 1710-1722
	Xu F, Zhang Y. Uniqueness and asymptotic stability of time-periodic solutions for the fractional Burgers equation. Acta Math Sci, 2023, 43A(6): 1710-1722
[7]	李仁华, 王征平. 含强制位势的分数阶薛定谔泊松方程的正规化解. 数学物理学报, 2023, 43A(6): 1723-1730
	Li R H, Wang Z P. Normalized solution of fractional Schrödinger-Poisson equations with coercive potential. Acta Math Sci, 2023, 43A(6): 1723-1730
[8]	Jannelli A. A finite difference method on quasi-uniform grids for the fractional boundary-layer Blasius flow. Mathematics and Computers in Simulation, 2024, 2.5(20): 382-398
[9]	Li D, Sun W, Wu C. A novel numerical approach to time-fractional parabolic equations with nonsmooth solutions. Numer Math Theor Meth Appl, 2021, 14(2): 355-376
[10]	Yu B, Li Y, Liu J. A Positivity-Preserving and robust fast solver for Time-Fractional Convection-Diffusion problems. Journal of Scientific Computing, 2024, 98(3): 1-26
[11]	Zheng Z Y, Wang Y M. Fast High-Order compact finite difference methods based on the averaged L1 formula for a Time-Fractional Mobile-Immobile diffusion problem. Journal of Scientific Computing, 2024, 99(2): 1-35
[12]	Hu J, Alikhanov A, Efendiev Y, et al. Partially explicit time discretization for time fractional diffusion equation. Fractional Calculus and Applied Analysis, 2022, 25(5): 1908-1924
[13]	Qiao L, Guo J, Qiu W. Fast BDF2 ADI methods for the multi-dimensional tempered fractional integrodifferential equation of parabolic type. Computers and Mathematics with Applications, 2022, 1.3(7): 89-104
[14]	Yin B, Liu Y, Li H, et al. A class of efficient time-stepping methods for multi-term time-fractional reaction-diffusion-wave equations. Applied Numerical Mathematics, 2021, 1.5(4): 56-82
[15]	Liao H L, McLean W, Zhang J. A second-order scheme with nonuniform time steps for a linear reaction-sudiffusion problem. Communications in Computational Physics, 2021, 30(2): 567-601
[16]	Yuste S B, Acedo L. An explicit finite difference method and a new von Neumann-type stability analysis for fractional diffusion equations. SIAM Journal on Numerical Analysis, 2005, 42(5): 1862-1874
[17]	Fan E, Wang J, Liu Y, et al. Numerical simulations based on shifted second-order difference/finite element algorithms for the time fractional Maxwell's system. Engineering with Computers, 2022, 38(Suppl 1): 191-205
[18]	Dehghan M, Abbaszadeh M. An efficient technique based on finite difference/finite element method for solution of two-dimensional space/multi-time fractional Bloch-Torrey equations. Applied Numerical Mathematics, 2018, 1.1(12): 190-206
[19]	Liao H, McLean W, Zhang J. A discrete Gronwall inequality with applications to numerical schemes for subdiffusion problems. SIAM Journal on Numerical Analysis, 2019, 57(1): 218-237
[20]	Fang Z, Zhao J, Li H, et al. A fast time two-mesh finite volume element algorithm for the nonlinear time-fractional coupled diffusion model. Numerical Algorithms, 2023, 93(2): 863-898
[21]	Lin J, Bai J, Reutskiy S, et al. A novel RBF-based meshless method for solving time-fractional transport equations in 2D and 3D arbitrary domains. Engineering with Computers, 2023, 39(3): 1905-1922
[22]	Li M, Shi D, Pei L. Convergence and superconvergence analysis of finite element methods for the time fractional diffusion equation. Applied Numerical Mathematics, 2020, 1.1: 141-160
[23]	Liao H, Li D, Zhang J. Sharp error estimate of the nonuniform L1 formula for linear reaction-subdiffusion equations. SIAM Journal on Numerical Analysis, 2018, 56(2): 1112-1133