收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

分数阶扩散方程的几种数值解法

苏丽娟  
【摘要】: 尽管分数阶微积分的历史几乎和整数阶的一样长,但是由于缺少相关的实际应用背景,分数阶微积分在其初期发展十分缓慢.众所周知,对于解释和模拟许多应用科学领域的动力学过程,经典微积分都是一个强有力的工具.但是,越来越多的实验和现实告诉我们,在自然界的反常动力学中有许多复杂系统,不能用经典的导数模型来描述.因此,在最近的十几年里,分数阶微积分已经被应用于几乎所有科学、工程和数学的领域中去. 物理中,反常扩散或许是一种最常研究的复杂问题.我们利用分数阶导数,可以将经典的整数阶扩散与波的偏微分方程,推广到时间和空间的分数阶上去.进而再扩展到各类非线性方程并给出其初边值问题的解,是近几年来分数阶微积分应用的一个主要领域.一般来讲,这些问题大都具非常重要的实际应用背景,如在分形和多孔介质中的弥散、半导体物理、湍流及凝聚态物理等. 本文主要研究一些分数阶扩散方程及其数值解法,共由四个彼此相关而又相互独立的章节构成.第一章简要介绍了分数阶微积分的历史、理论及其应用,以及文中将用到的一些基本知识,和相关数值解的现有研究成果;第二章和第三章研究的都是双边空间分数阶对流扩散方程,在这两章中我们分别给出了此类方程的几类不同有限差分法,主要有分数阶权平均法、改进型权平均法和特征有限差分法等;而在最后一章中,我们则是给出了时间分数阶扩散方程的一种高精度隐式数值解法. 第一章为序言.首先了介绍分数阶微积分的历史及其发展情况,并给出了几种常用的分数阶算子定义以及它们的一些基本性质,例如:Riemann-Liouville分数阶算子,Caputo分数阶算子和Grunwald-Letnikov分数阶算子等,同时还列出了几个相关的运算性质.然后,在§1.3中,我们对Mittag-Leffler型特殊函数和它的基本性质也进行了一定的叙述,这类特殊函数主要包括单参数的Mittag-Leffler函数和两参数的广义Mittag-Leffler函数.这类特殊函数常常是很多分数阶微分方程的基本解,其它具有类似性质的特殊函数还有Wright函数和H-fox函数,等等. 此外,在本章中的§1.4,我们还归纳叙述了目前为止,几类常见的分数阶微分方程的一些数值解法.例如,有限差分法,有限元法,微分变换法,Adomian区域分解法,变分迭代法,同伦摄动法,等等.同时对每种方法分别列举出了一些相关研究成果.最后,在本章的最后一节中,我们较详细的介绍了分数阶微积分在当前非线性物理复杂系统的各个领域中的应用. 在接下来的章节中,我们将研究两种不同的反常扩散模型.在第二章中,主要研究1维空间分数阶对流扩散方程.我们根据移位Grunwald公式离散Riemann-Liouville分数阶导数,从而提出了方程的分数阶权平均法.通过理论研究和算例分析,可以得知以前出现过的一些相关数值算法,它们大都是此方法的某些特例.在§2.3中,我们利用圆盘定理和矩阵法证明了分数阶权平均法的稳定性,具体理论结果由定理2.1详细给出. 然后,在§2.4中,我们又讨论了分数阶权平均法的一种新的改进格式,并再次给出了相关稳定性分析.最后,则是用数值例子来验证理论的正确性,同时又计算了分数阶权平均法的特例,分数阶Crank-Nicholson (FCN)法.显然,无条件稳定又拥有2阶时间精确度的FCN法更好一些.本章部分内容已经公开发表在Physics Letters A. 在第三章中,关于双边空间分数阶对流扩散方程,据我们所知,目前为止,它的数值解法全都是Eulerian法.结果,这些方法都具有和2阶对流扩散方程相同的数值局限性.在本章中,结合移位Griinwald-Letnikov有限差分过程以及Lagrangian法,我们在§3.3中首次提出了一种分数阶特征有限差分法(CFDM).此法保留了2阶对流扩散方程特征法和分数阶对流扩散方程有限差分法的所有数值优点.在§3.4中,我们证明了这种方法是无条件稳定、相容和收敛的,并且给出了本方法误差估计的最大值. 在§3.5中,我们给出了一个实际算例的数值模拟,并把分数阶特征有限差分法和其它的分数阶标准差分法相比较.算例结果表明,这种分数阶新CFDM在精度和稳定性上都大大优于其它已知方法,例如显式迎风差分法和隐式迎风差分法等.并且,此法对于对流占优问题,显得尤为高效、优越.本章内容已投到Journal of Computational Physics. 在第四章中,我们主要讨论的是一类时间分数阶扩散方程.反常次扩散运动是复杂系统中一个特别重要的内容,如在一些有机和无序材料中,它的运动路径被一些几何或能量因子约束着.对于反常次扩散随机游走过程的数学模型,一般扩散方程则会被Riemann-Liouville分数阶时间扩散方程所替代.分析表明,这些分数阶模型显然比经典的整数阶模型更加符合实际背景. 在本章中,首先我们我们利用移位Grunwald公式来逼近时间分数阶导数,并且使用中心差分格式去逼近1阶时间导数和2阶空间导数,从而提出了此类扩散方程的一种新的隐式差分法.它是一种三层差分格式,其中第一时间层的数值解可以由全隐式格式或Crank-Nicholson格式给出,这两种格式都是无条件稳定的.接着,我们利用广义化Fourier-Von Neumann分析法,证明了这种新方法的无条件稳定性,并导出此法关于时间的2阶精确度.最后,则是用数值试验和对比法,来验证和观察本算法的性质和特征.本章内容已投到Applied Mathematics and Computation.


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 杨光俊;分数阶积分方程[J];云南大学学报(自然科学版);1997年03期
2 张淑琴;有限区间上的分数阶扩散波方程的解[J];西北师范大学学报(自然科学版);2005年02期
3 冯康,王烈衡;非协调元的分数阶Sobolev空间[J];计算数学;1995年03期
4 卢旋珠;时间分数阶对流-扩散方程的有限差分方法[J];福州大学学报(自然科学版);2004年04期
5 刘甲国,徐明瑜;人颅骨粘弹性的分数阶模型研究[J];中国生物医学工程学报;2005年01期
6 刘慈群;非牛顿流体的广义Maxwell模型及其解[J];力学与实践;1995年02期
7 许晓军,陆启生;分数阶傅里叶变换在球面谐振腔稳定性分析中的应用[J];激光与光电子学进展;1999年03期
8 陈喆,王宏禹,邱天爽;基于分数阶傅立叶变换的模糊函数的研究[J];信号处理;2003年06期
9 胡孝博;球面波照射下实现非精确分数阶傅里叶变换[J];南京工程学院学报(自然科学版);2005年01期
10 白亿同,罗玉芳,胡永旭;Weierstrass类型函数的光滑度与分形维数之间的联系初探[J];工科数学;1995年01期
11 刘莉;周朴;;基于分数傅立叶变换的图像隐藏技术[J];国防科技大学学报;2005年06期
12 同登科,石丽娜;多孔介质中粘弹性液体广义流动分析[J];水动力学研究与进展A辑;2004年06期
13 周玉鼎;斯仁道尔吉;;同伦摄动法与几类分数阶积分微分方程的解[J];内蒙古师范大学学报(自然科学汉文版);2010年01期
14 闵富红;王执铨;;分数阶混沌系统同结构与异结构广义同步[J];控制与决策;2008年09期
15 包革军,凌怡;一类缺项的负系数解析函数[J];哈尔滨工业大学学报;1995年03期
16 陈明杰;分数阶傅里叶变换的数值实现[J];重庆大学学报(自然科学版);2003年05期
17 赵义正,杨景曙;基于分数阶Fourier变换的瞬时频率估计方法[J];安徽大学学报(自然科学版);2005年01期
18 唐建国,陆磊,邓冠铁,赵尔敦;一类函数图像的填充维数与分数阶导数[J];华中理工大学学报;1997年09期
19 王振滨,曹广益,朱新坚;基于状态空间模型分解的分数阶系统辨识算法[J];系统工程与电子技术;2004年12期
20 谢庭藩;分数阶连续模与Fourier级数的绝对收敛[J];数学研究与评论;1981年S1期
中国重要会议论文全文数据库 前10条
1 张云;赵昕;姜义成;;舰船目标的分数阶自相关瞬时成像算法研究[A];全国第二届信号处理与应用学术会议专刊[C];2008年
2 寇春海;严烨;刘健;;分数阶微分方程稳定性分析及其在HIV-1感染模型中的应用[A];中国力学学会学术大会'2009论文摘要集[C];2009年
3 苏丽娟;;双边空间分数阶对流—扩散方程的一种有限差分解法[A];中国力学学会学术大会'2009论文摘要集[C];2009年
4 陈文;;软物质和分数阶导数建模[A];中国工程物理研究院科技年报(2003)[C];2003年
5 李惠光;王平顺;李国友;;采用分数阶导数的高斯滤波微分算子边缘检测方法[A];中国自动化学会全国第九届自动化新技术学术交流会论文集[C];2004年
6 王在华;;分数阶动力系统的若干问题[A];第三届全国动力学与控制青年学者研讨会论文摘要集[C];2009年
7 尹晨;徐明瑜;谭文长;;药物控释系统中一类分数阶双可动边界问题的近似解[A];第九届全国生物力学学术会议论文汇编[C];2009年
8 曹军义;薛士明;;分数阶碰摩转子系统的非线性动力学分析[A];第十二届全国设备故障诊断学术会议论文集[C];2010年
9 魏延玉;官玉彬;王文祥;;分数阶复宗量贝塞尔函数的数值计算[A];中国电子学会真空电子学分会第十一届学术年会论文集[C];1997年
10 陈文剑;孙辉;张明辉;朱建军;朱广平;;掩埋目标探测混响抑制方法研究[A];2009年度全国物理声学会议论文集[C];2009年
中国博士学位论文全文数据库 前10条
1 杨水平;几类非线性分数阶微分方程初值问题的数值求解及动力性质[D];湘潭大学;2011年
2 杨雪花;正交样条与拟小波配置法在分数阶偏微分方程数值解中的应用[D];湖南师范大学;2014年
3 李亚宁;两类分数阶微分方程定解问题解的研究[D];兰州大学;2014年
4 赵潇;分数阶粘弹性本构参数识别[D];大连理工大学;2012年
5 陈喆;自适应与分数阶非平稳信号处理的研究[D];大连理工大学;2003年
6 李静;矩阵方程的正定解和分数阶微分方程的谱问题[D];山东大学;2014年
7 陈韵;分数阶Fourier变换在水声通信中的应用研究[D];哈尔滨工程大学;2012年
8 郭霄怡;分数阶微积分在量子力学和非牛顿流体力学研究中的某些应用[D];山东大学;2007年
9 王俊刚;时间分数阶扩散方程的两类不适定问题研究[D];兰州大学;2013年
10 尉宇;线性调频和非线性调频信号的检测与参数估计[D];华中科技大学;2005年
中国硕士学位论文全文数据库 前10条
1 汪辉;分数阶小波变换理论及其应用研究[D];西南交通大学;2014年
2 鲁卫善;生物组织传热的分数阶模型及其某些应用[D];山东大学;2014年
3 汪向艳;空间分数阶扩散方程的几种数值解法[D];宁夏大学;2014年
4 邱亮亮;交点间断有限元求解分数阶微分方程[D];兰州大学;2014年
5 贾贺芝;分数阶Cable方程的高精度数值解法[D];湘潭大学;2011年
6 管玉洁;变阶分数阶偏微分方程的非一致时间步长有限差分算法[D];山东大学;2014年
7 王佃萍;快速傅里叶变换在分数阶扩散方程中的并行实现[D];山东大学;2014年
8 叶超;分数阶偏微分方程的高阶数值格式与理论分析[D];湘潭大学;2012年
9 肖玲玲;非线性分数阶反应子扩散方程数值解法[D];湘潭大学;2010年
10 乔晓梅;分数阶脉冲微分方程解的存在性[D];山东师范大学;2014年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978