空间四阶时间亚扩散方程有限差分方法.pdfVIP

2014年 3月 第 28卷第 1期 Mar．2014 CommunicationonAppliedMathematicsandComputation V0l1．28 NO．1 DOI 10．3969／j．issn．1006—6330．2014．01．013 空间四阶的时间亚扩散方程的有限差分方法 郭 建， 李常品， 丁恒飞 (上海大学理学院，上海 200444) 摘要 提出了两个求解空间四阶的时间亚扩散方程的数值方法，其误差阶分别为 0(7-+h) 和 0(7-+h2)．通过 Fourier方法，发现两个差分格式均为无条件稳定的．最后，通过数值 例子，验证了两个算法的有效性． 关键词 分数阶；Fourier方法；亚扩散方程；有限差分方法 2010数学分类号 65M06；65M12 中图分类号 O175．26 文献标志码 A 文章编号 1006—6330(2014)01—0096-13 Finitedifferencemethodsfortimesubdi胁 sion equationwith spacefourth-order GUO Jian， LIChang-pin， DING Heng—fei (CollegeofSciences，ShanghaiUniversity，Shanghai200444，China) Abstract TWO numericalschemesforatimesubdiffusion equationwith space fourth．orderareproposed，andtheconvergenceordersareO(T+h。)and0(~-2十2)， respectively．ByusingtheFouriermethod，itiSofundthattwo finitedifierence schemesareallunconditionallvstable．Finally．numericalexamplesaregivento testifytheemciencyofthenumericalschemes． Keywords fractionalorder；Fouriermethod；subdiffusionequation；finitediffer— encemethod 2010M athematicsSubiectClassification 65M06；65M12 ChineseLibrary Classification O175．26 0 引 言 随着科学技术的发展，经典微积分不是总能够很好地解释和模拟许多应用领域的动 力学过程，如扩散方程，所以分数阶微积分被众多学者关注并寄予厚望，在最近几十年， 通过分数阶微积分建立起的微分方程模型被广泛应用于科学工程领域 【一引． 收稿 日期 2013—06．02； 修订 日期 2013．10—13 基金项目 国家自然科学基金资助项目；上海市教育委员会科研创新重点资助项目(12ZZ084) 通信作者 李常品，研究方向为分数阶偏微分方程的数值方法、分岔混沌的应用理论和计算． E-mail：lcp@shu．edu．cn 第 1期 郭 建，等：空间四阶的时间亚扩散方程的有限差分方法 97 在有关分数阶的应用中，尤其分数阶扩散方程在水文、力学以及物理学中的应用是 近年来关注的重要研究课题．经过多位学者的努力，通过建立分数阶通量和流的本构关 系以及分数阶随机游走的广义概念后，从这两个方向又同时给出了分数阶扩散方程的一 致形式 [3]_而对求解分数阶扩散方程已有几种不同的方法 [1]_例如，Fourier变换法、 Mellin变换法、Laplace变换法等，但是只有很少类型的分数阶微分方程能够应用以上 方法求出精确解．因此，大多情况下，只能使用数值方法来进行计算． 在分数阶微分方程数值解法中，有限差分法是最常用的一类算法 [4．5