空间分数阶扩散方程的谱及拟谱方法.pdfVIP

第38卷 第3期 应 用 数 学 学 报 Vo1．38No．3 2015年 5月 ACTAMATHEMATICAEAPPLICATAESINICA May，2015 空间分数阶扩散方程的谱及拟谱方法 郑敏玲 (湖州师范学院理学院，湖州313000) (E—mail：mlzheng@aliyun．corn) 摘 要 本文研究了空间分数阶扩散方程的数值解法．首先，利用 Galerkin谱方法对方程在空 间方向进行离散，然后，利用 一方法对时间变量进行离散．利用能量方法，得到了全离散方程 解的稳定性及误差估计．数值结果表明拟谱方法求解分数阶扩散方程数值解高精度和有效性． 关键词 空间分数阶扩散方程；稳定性；误差估计；谱方法 MR(2000)主题分类 26A33；65M06；65N12 中圉分类 O241．8 1 引言 空间分数阶扩散方程 Ou(x,t) = u( )+ ， 1 2， (1) 是与熟知的扩散过程 (亦称布朗运动)相对的反常扩散过程的一种模型，亦称超扩散方 程．这里，分数阶导数按Riemann—Liouville导数定义 Oxa = 一㈨“)击而 一dx『2．n／(_一s))～卜u(s，t)ds．· 另一种反常扩散模型称为时间分数阶扩散方程，亦称次扩散方程，可以写成 O2u( x,t) — ： — 一Ot ：O一ti一一Ox2 ++’J，I＼，r，’)，I，’ 0ll(、～』、1l，’ I(2’。，)l 关于分数阶扩散方程的数值解，一些学者对此做了不少研究．Langlands．Henry基 于分数阶导数的L1～算法研究了次扩散方程 [1】1提出了一种无条件稳定的隐式差分方 法．Yuste—Acedo基于Griinwald—Letnikov公式，提出了次扩散方程的显式方法 2【13]．关 本文 2013年 7月 19日收到．2013年 12月29日收到修改稿． 国家自然科学基金 11101140)以及湖州市自然科学基金 (2013YZ06)资助项 目 3期 郑敏玲：空间分数阶扩散方程的谱及拟谱方法 435 于空间分数阶扩散方程也有过一些研究．有趣的是，基于 Grfinwald—Letnikov公式，空 间分数阶扩散方程的隐式及显式差分方法都是不稳定的，而稳定的差分格式可以通过 移位的Griinwald—Letnikov公式得到 】．最近的文献可以参考 6『-8]．结合时间分数阶 导数的L1一算法和空间分数阶导数的移位算法， 9『，101对空间 时间分数阶方程的数值 解法进行了研究． 由于分数阶导数的非局部性质，低阶的方法用于计算分数阶导数时大大地增加了 计算的复杂度及存贮空间．因此，寻找高阶的方法是数值求解分数阶方程的必然选择． 众所周知，在偏微分方程领域里，谱方法经常被采用，因为它在许多方面要优于差分 方法 [11-13]．基于分数阶扩散方程 (1)及 (2)的变分结构，Roop和 Erin建立了分数 阶扩散方程的变分理论 i14,1~】．从此，Galerkin方法被许多学者应用于分数阶扩散方程 [16-1s】 ． Zhang—Liu-Anh运用谱方法对空间变量进行离散，然后利用隐式差分方法离散时 间变量，研究了一类空间分数阶扩散方程 [18]．据我们所知，除上面提到的文献外，有 关空间分数阶扩散方程的谱方法的出版物仍