分数阶反应-子扩散方程高阶隐式差分格式及稳定性分析.pdfVIP

分数阶反应-子扩散方程的高阶隐式差分格式及其稳定性分析 梁娜 叶超 咸宁学院数学与统计学院,湖北咸宁 437100 湘潭大学数学与计算科学学院,湖南湘潭411105 摘要 针对一类带初边值条件的分数阶反应-子扩散方程，构造了一种新的高阶的隐式差分格式，其局部截断误差 为( + )，并对格式的可解性做了分析，利用Fourier方法证明了格式的无条件稳定性。最后做了数值算例，从所 给的表格数据可以得出，该格式的具有非常高的精度。 关键词 分数阶反应-子扩散方程；隐式差分格式；Riemann-Liouville分数阶导数；稳定性 中图分类号：O241. 8 文献标识码：A A higher-order implicit difference scheme and its stability analysis for the fractional reaction-subdiffusion equation Department of Mathematics, Xianning College, Hubei Xianning, 437100, China Department of Mathematics, Xiangtan University, Hunan Xiantan, 411105, China A kind of fractional reaction-subdiffusion equation with initial-boundary conditions is considered. Based on the definition of Riemann-Liouville fractional derivative, a higher-order implicit finite scheme is presented. By means of Fourier method, the stability of the scheme is provided. Finally, a numerical example is given to support the theoretic analysis. fractional reaction-subdiffusion equation; Riemann-Liouville fractional derivative; implicit finite scheme; stability 分数阶微分方程被广泛地应用到自然科学与工程学领域，引起了科研人员的重视。Ralf Met- zler和Joseph Klafter在其综述性文章[1] 中指出分数阶微分方程已经成为描述反常传播过程一种基本工具。 J. Tenreiro Machado et al.介绍了分数阶微积分理论在近三十年的发展概况[10] ，分数阶导数，分数阶微 分方程的解法与应用的一些基本理论[2]可见于I. Podlubny 的经典著作。近年来，有关分数阶微分方程的 解析解问题[3, 4]取得了很多成果，然而这些解时常包含一类特殊而又复杂的函数，如Mittag-Leffler函数， 这给实际的应用带来不便。因此有必要去研究分数阶微分方程的数值解法。 迄今为止，分数阶扩散方程的高阶数值方法很少。F. Liu et al. 讨论了一种非线性反常扩散方程的数 值解法[6] ，提出了一种时间一阶，空间二阶的隐式差分格式，并利用Fourier方法对格式做了理论分析。P. ZHUANG et al.利用积分的手段对一反常子扩散方程[7]进行离散，得到了一种时间一阶，空间二阶的隐 式差分格式，采用能量方法给出了格式的稳定性与收敛性分析，并通过改进，提高了时间误差阶。M. R. Cui对一种分数阶扩散方程首次采用了一种紧致有限差分方法[5] ，得到了一种新的空间四阶，时间一阶的 隐式差分格式。Chang-Ming Chen et al.考虑了一种变阶的反常子扩散方程[8] ，基于Riemann-Liouville导 数的Gr¨uwald-Letnikov离散以及结合四阶的紧致