微分求积方法解时间分数阶扩散方程.docVIP

微分求积方法解时间分数阶扩散方程 曹焕 张学莹 刘荟 河海大学理学院 X 关注成功！ 加关注后您将方便地在 我的关注中得到本文献的被引频次变化的通知！ 新浪微博 腾讯微博 人人网 开心网 豆瓣网 网易微博 摘????要： 微分求积 (DQ) 法是一种基于径向基函数 (RBFs) 插值的无网格方法.本文选MultiQuadrics (MQ) 函数作为径向基函数, 并采用微分求积方法解决时间分数阶扩散方程。在离散过程中, 采用有限差分法离散时间项, 采用微分求积方法离散空间项.最后, 结合数值求解的结果做出相应的误差分析. 关键词： 径向基函数; MQ函数; 微分求积法; 时间分数阶扩散方程; 作者简介：曹焕 (1992-) , 女, 山东菏泽人, 在读硕士研究生, 研究方向:计算数学 作者简介：张学莹 (1973-) , 男, 山东郯城人, 副教授, 硕士生导师, 研究方向:计算数学, zhangxy@ 收稿日期：2017-08-27 基金：教育部留学回国人员科研启动基金项目 (20145003412) The differential quadrature method solutions for the time-fractional diffusion equation CAO Huan ZHANG Xue-ying LIU Hui School of Science, Hohai University; Abstract： The differential quadrature method is an alternative radial basis functions (RBFs) meshless method.This article selects MQ function as the radial basis function, and applys the differential quadrature method to solve the time fractional diffusion equation.In the discretization formulation, a finite difference scheme and the DQ are used repectively to discretize time fractional derivative and spatial derivative terms.Finally we make the error analysis with the results of the numerical investigation example. Keyword： RBFs; MQ function; DQ method; time fractional diffusion function; Received： 2017-08-27 0 引言 近年来, 反常扩散现象在物理学和工程学等领域得到了广泛的关注[1,2], 比如:污染物的运输、渗流、损耗以及磁等离子体等问题.与正常的扩散现象相比, 反常扩散现象具有较强的远程相互作用与历史依赖性等特点, 而标准的整数阶偏微分方程已经无法准确的描述这类反常扩散行为[3].随后, 分数阶导数被陆续的证实是一种可以准确的描述这类反常扩散现象的数学模型[4].所以, 为了解决这些问题, 引入了分数阶扩散方程. 由于分数阶偏微分方程中含有分数阶的导数, 不容易求解分数阶偏微分方程的解析解, 因此很多人研究其数值解法.目前, 有限差分方法[5]已被证实为对分数阶扩散方程进行时间和空间上离散的一种有效方法.对于分数阶扩散方程的数值解法, 也可以选择其他合适的方法对其空间上的离散, 例如:边界元法[6], 有限元法[7,8], 傅里叶法[9], 无网格径向基函数法[10-13]等.与传统的有限差分方法对空间的离散相比较, 对于处理一些大规模问题, 这些方法在一定程度上减少了计算量, 同时又提高了计算效率.2003年, Shu C提出DQ法解整数阶的偏微分方程, 且取得较理想的结果.本文将这一方法运用到解分数阶偏微分方程中, 并给出误差结果分析. 1 基本理论 1.1 时间分数阶扩散方程 边界条件: 初始条件: 式 (1) ~ (3) 中:Q (X, t) , g (X, t) , u0 (x) 是给定的函数, T是总时间, /t是关于时间t的分数阶导. 本文采用Caputo时间分数阶导数的定义, 定义如下: 这里0α1. 1.2 FDM方法对时间的离散 首先把区间[0, T]分割成K个小区间, 定义tk=kτ, k=0, 1, …, K, 这里τ=T/K为时间步长.用u