变系数空间分数阶对流-扩散方程的有限差分解法.pdfVIP

第25卷第4期 沈阳大学学报( 自然科学版) V01．25，No．4 201 of 013 3年8月 Journal Science) ShenyangUniversity(Natural Aug．2 文章编号：2095—5456(2013)04—0341—04 变系数空间分数阶对流一扩散 方程的有限差分解法 马亮亮 (攀枝花学院数学与计算机学院，四川攀枝花617000) 摘 要：考虑了一个变系数空间分数阶对流一扩散方程．这个方程是将一般的对流一扩散方程中的空间二 阶导数用口(1％fl≤2)阶导数代替．提出了一个隐式差分格式，验证了这个差分格式是无条件稳定的，并证明 了它的收敛性，其收敛阶为o(r+^)，最后给出了数值例子． 关键词：对流一扩散方程；分数阶导数；隐式差分；稳定性；收敛性 中图分类号：O241．82 文献标志码：A 近年来，分数阶微分方程引起了人们的广泛 (2) 关注，应用领域也变得很广泛，包括分形动力学、 扩散和运输、混沌与湍流、随机游走、经济金融、随 1 差分格式的建立 机过程、非牛顿流体力学、材料记忆性、自相似和 多孔结构的动态过程等诸多领域[1]． 作网格剖分，令r，h分别为时间和空间步长， 在标准对流一扩散方程中，用分数阶导数 即z：一i^(声o，1，…，M；^一击)㈡一挖r 口(1p≤2)代替空间二阶导数，可得到空间分 数阶对流一扩散方程．目前，国内外学者在分数阶 (n—o，1，…，N；r一寿)．令M?为点“(z：，￡。)的近 对流一扩散方程及其衍生方程的理论分析、数值算 法、工程应用等方面做了许多有意义的研究 空间一阶导数，采用一阶差商逼近： 工作E2—4I． 堕g掣碰一丛生盘止二丛趔+D(r)，(3)‘…～、。 本文考虑如下变系数空间分数阶对流一扩散 氆 r 方程： 堕擘监立一丛堕鱼立手地立巫生+o(^)． 掣一妇∽掣+ dZ n (4) d(x,t)≮产+m∽ (Ot≤T，0zL)， 壁喽暑盟，采用G，nnwald改进型公式替代嘲： b2‘ (0zL)， u(x，O)一g(x) 学=h-一X “(o，f)一o，鱼丝笋一o． (1) 式中． 其中，1t?42，且州≥o；鼍≯Rie…n— Liouville空间分数阶导数‘5|： 叫尹一(_1)‘殴尘等笋幽． 收稿日期：2013—01—04 基金项目：国家自然科学基金资助项目；四川省科技厅资助项目(2013JY0125) 作者简介：马亮亮(1986一)，男，甘肃天水人，攀枝花学院助教． 万方数据 342 沈阳大学学报(自然科学版) 第25卷 将式