一个分数阶扩散方程的定解问题的数值解法.pdf

第26卷第5期 四川理工学院学报(自然科学版) V01．26No．5 兰Q!i圭!Q旦 !!竺翌坐!!!!!!竺竺旦兰盐：翌i堡坚!!堡翌!：垒兰翌垂旦!!墨兰gl盟苎兰型!!兰呈竺兰!i垒竺翌2竺!!：：旦!i 文章编号：1673-1549(2013)05-0086-04 一个分数阶扩散方程的定解问题的数值解法 池光胜1，李功胜2，张芳1，张涛1 (1．山东凯文科技职业学院基础教学部，济南250200；2．山东理工大学理学院应用数学研究所，山东淄博255000) 摘要：研究了一个扩散系数与空间变量相关的一维空间一时间分数阶扩散方程的定解问题。 隐式差分格式，验证了这个格式是无条件稳定，并证明了它的收敛性，其收敛的阶为D(丁+h)，最 后给出了数值例子。 关键词：空间一时间分数阶扩散方程；扩散系数与空间变量相关；隐式差分格式；稳定性；收敛性 中图分类号：0242．1 文献标志码：A 通常的扩散方程为整数阶微分方程，分数阶微分方程 产生于一些反常扩散模型。在物理，工程，金融及环境问 扩散系数且D’(石)0，∥(菇)0为对流系数，S(菇，t) 题等方面得到广泛应用。最近几年在国际上掀起了一股 为源项。其中时间分数阶导数取Caputo导数： 求解分数阶微分方程研究热潮。uuF等人…通过变量变 吖7 …’ o 1 Ol 掣Ot=志FJro(￡一叼)一n曼堡铲d_，7 (一)” 咖 M 换得到分数阶对流色散方程的解；MeerschaertM等 人旧。用有限差分的方法求解分数阶偏微分方程两点边值 空间分数阶导数取Riemann—Liouville导数： 问题。数值解方面，覃平阳，张晓丹一一提出了一个具体的 ” 7 P1 1 些Ox掣=d1可旦Ox3广o㈦矿“警鸳 厂(一y) 喏’ 空间一时问分数阶对流一扩散方程的隐式差分格式，验证 了格式的稳定性和收敛陛。陈世平，刘发旺H。求解了一类 2建立差分方法 一维分数阶渗透方程，给出了隐式差分格式，并做了数值 模拟。基于上述研究，本文拟采用数值方法求解一个扩散 k=0，1，2，…，Ⅳ，其中，h和下分别是空间和时间步长， 系数与空间变量相关的空间一时间分数阶扩散方程，给出 隐式差分格式，并给出稳定性和收敛性的证明。 定方程用差商代替一阶导数： 1 问题描述 Ou(菇i，t川)Mki”一MkH+l 0x h 考虑空间一时间分数阶对流一扩散方程的初边值 对Caputo导数采用如下近似： 问题 a“u(戈。，t川) 7I-a 鼍≯=击(嘶，％≯)一 Ot“ F(2一a) 。(茹)掣+。(刈)