时间分数阶偏微分方程高阶数值解法

发布时间：2017-09-23 15:22

  本文关键词：时间分数阶偏微分方程高阶数值解法

  更多相关文章： 分数阶导数 时间分数阶反常低扩散方程 时间分数阶扩散波方程 带第一类Dirichlet边值条件的四阶分数阶扩散方程 有限差分格式 收敛性 稳定性

【摘要】：本文主要研究几类时间分数阶扩散方程的高精度数值解法并给出相应的误差估计式.首先,研究一维时间分数阶反常低扩散方程的高阶数值算法,并给出相应的先验估计式.基于带权和位移的Griinwald-Letnikov算子[48,64],选取位移(p,q,r)=(0,-1,-2),利用Riemann-Liouville分数阶导数与Caputo分数阶导数在一定的光滑条件下等价,推导出一种新的三阶逼近公式来离散Caputo分数阶导数.空间方向结合平均算子技术,对方程构造了一种新的紧差分格式(称为GL3差分格式).利用离散能量分析方法,严格证明了GL3差分格式无条件稳定和收敛.其次,讨论二维时间分数阶反常低扩散方程的高阶数值算法.遵循一维时间分数阶反常低扩散方程有限差分格式的构造技巧,对二维问题建立时间和空间均为高阶的差分格式,并分析了差分格式在L1(L∞)范数意义下是无条件稳定和收敛的.对于二维情形,由于H1半范不能嵌入到无穷模,因此,二维问题的理论分析技巧与一维问题有着本质的区别.同时,通过添加小量项,对二维时间分数阶反常低扩散方程建立紧交替方向隐式差分格式(称为紧ADI格式),通过数值实验验证了紧ADI格式的收敛阶和有效性.再次,研究时间分数阶扩散波方程的有效差分格式.在方程两端作用Riem-an-n-Liouvile分数阶积分算子,将原方程转化为积分-微分方程.利用线性插值思想,推导出一种新的二阶离散公式逼近Riemann-Liouville分数阶积分算子,并对等价方程建立高阶的数值格式.通过数值算例验证该差分格式的收敛阶和有效性.最后,研究带第一类Dirichlet边值条件的四阶分数阶扩散方程的空间高精度的数值算法.在点x0,x1,x2,x3处和点xM,xM-1,xM-2,xM-3处分别对：xx(x,t)进行线性组合,推导出第一类Dirichlet边值条件的具有四阶精度的离散公式,空间方向内部点作用平均算子,用L1离散公式逼近时间分数阶导数,对方程建立空间一致四阶的差分格式.利用数学归纳方法和离散能量分析方法证明差分格式的稳定性和收敛性.
【关键词】：分数阶导数 时间分数阶反常低扩散方程 时间分数阶扩散波方程 带第一类Dirichlet边值条件的四阶分数阶扩散方程 有限差分格式 收敛性 稳定性
【学位授予单位】：东南大学
【学位级别】：硕士
【学位授予年份】：2015
【分类号】：O241.82
【目录】：

• 摘要4-5
• Abstract5-9
• 第一章 绪论9-18
• 1.1 分数阶微分算子的定义及性质9-11
• 1.2 研究背景及现状11-16
• 1.3 本文研究动机及主要工作16-18
• 第二章 一维时间分数阶反常低扩散方程的高阶紧差分格式18-38
• 2.1 引言18-19
• 2.2 高阶紧差分格式的建立19-25
• 2.3 差分格式稳定性和收敛性分析25-33
• 2.4 数值实验33-36
• 2.5 结论36-38
• 第三章 二维时间分数阶反常低扩散方程的高阶紧差分格式38-55
• 3.1 引言38-39
• 3.2 高阶紧差分格式的建立39-43
• 3.3 差分格式稳定性和收敛性分析43-48
• 3.4 紧交替方向格式48-52
• 3.5 数值实验52-53
• 3.6 结论53-55
• 第四章 时间分数阶扩散波方程的高阶紧差分格式55-64
• 4.1 引言55
• 4.2 Riemann-Liouville积分算子_0D_t~(-α)f(t)的二阶离散公式55-57
• 4.3 高阶紧差分格式的建立57-60
• 4.4 数值实验60-63
• 4.5 结论63-64
• 第五章 带第一类Dirichlet边值条件的四阶分数阶扩散方程的高阶紧差分格式64-78
• 5.1 引言64
• 5.2 高阶紧差分格式的建立64-70
• 5.3 高阶紧差分格式稳定性和收敛性分析70-75
• 5.4 数值实验75-76
• 5.5 结论76-78
• 第六章 总结与展望78-80
• 参考文献80-86
• 附录A 攻读硕士学位期间完成的工作及获得的荣誉86-87
• 附录B 致谢87

【相似文献】

中国期刊全文数据库 前10条

1 王德金;郑永爱;;分数阶混沌系统的延迟同步[J];动力学与控制学报;2010年04期

2 杨晨航,刘发旺;分数阶Relaxation-Oscillation方程的一种分数阶预估-校正方法[J];厦门大学学报(自然科学版);2005年06期

3 王发强;刘崇新;;分数阶临界混沌系统及电路实验的研究[J];物理学报;2006年08期

4 夏源;吴吉春;;分数阶对流——弥散方程的数值求解[J];南京大学学报(自然科学版);2007年04期

5 张隆阁;;一类参数不确定混沌系统的分数阶自适应同步[J];中国科技信息;2009年15期

6 陈世平;刘发旺;;一维分数阶渗透方程的数值模拟[J];高等学校计算数学学报;2010年04期

7 辛宝贵;陈通;刘艳芹;;一类分数阶混沌金融系统的复杂性演化研究[J];物理学报;2011年04期

8 黄睿晖;;分数阶微方程的迭代方法研究[J];长春理工大学学报;2011年06期

9 蒋晓芸,徐明瑜;分形介质分数阶反常守恒扩散模型及其解析解[J];山东大学学报(理学版);2003年05期

10 陈玉霞;高金峰;;一个新的分数阶混沌系统[J];郑州大学学报(理学版);2009年04期

中国重要会议论文全文数据库 前10条

1 李西成;;经皮吸收的分数阶药物动力学模型[A];中国力学学会学术大会'2009论文摘要集[C];2009年

2 谢勇;;分数阶模型神经元的动力学行为及其同步[A];第四届全国动力学与控制青年学者研讨会论文摘要集[C];2010年

3 张硕;于永光;王亚;;带有时滞和随机扰动的不确定分数阶混沌系统准同步[A];中国力学大会——2013论文摘要集[C];2013年

4 李常品;;分数阶动力学的若干关键问题及研究进展[A];中国力学大会——2013论文摘要集[C];2013年

5 李常品;;分数阶动力学简介[A];第三届海峡两岸动力学、振动与控制学术会议论文摘要集[C];2013年

6 蒋晓芸;徐明瑜;;时间依靠分数阶Schr銉dinger方程中的可动边界问题[A];中国力学学会学术大会'2009论文摘要集[C];2009年

7 王花;;分数阶混沌系统的同步在图像加密中的应用[A];第二届全国随机动力学学术会议摘要集与会议议程[C];2013年

8 王在华;;分数阶动力系统的若干问题[A];第三届全国动力学与控制青年学者研讨会论文摘要集[C];2009年

9 张硕;于永光;王莎;;带有时滞和随机扰动的分数阶混沌系统同步[A];第十四届全国非线性振动暨第十一届全国非线性动力学和运动稳定性学术会议摘要集与会议议程[C];2013年

10 李西成;;一个具有糊状区的分数阶可动边界问题的相似解研究[A];中国力学大会——2013论文摘要集[C];2013年

中国博士学位论文全文数据库 前10条

1 陈善镇;两类空间分数阶偏微分方程模型有限差分逼近的若干研究[D];山东大学;2015年

2 任永强;油藏与二氧化碳埋存问题的数值模拟与不确定性量化分析以及分数阶微分方程的数值方法[D];山东大学;2015年

3 蒋敏;分数阶微分方程理论分析与应用问题的研究[D];电子科技大学;2015年

4 卜红霞;基于分数阶傅里叶域稀疏表征的CS-SAR成像理论与算法研究[D];北京理工大学;2015年

5 杨变霞;分数阶Laplace算子的谱理论及其在微分方程中的应用[D];兰州大学;2015年

6 邵晶;几类微分系统的定性理论及其应用[D];曲阜师范大学;2015年

7 方益;分数阶Yamabe问题的一些紧性结果[D];中国科学技术大学;2015年

8 王国涛;几类分数阶非线性微分方程解的存在理论及应用[D];西安电子科技大学;2014年

9 陈明华;分数阶微分方程的高阶算法及理论分析[D];兰州大学;2015年

10 孟伟;基于分数阶拓展算子的灰色预测模型[D];南京航空航天大学;2015年

中国硕士学位论文全文数据库 前10条

1 黄志颖;非线性时间分数阶微分方程的数值解法[D];华南理工大学;2015年

2 赵九龙;基于分数阶微积分的三维图像去噪增强算法研究[D];宁夏大学;2015年

3 楚彩虹;单载波分数阶傅里叶域均衡系统及关键技术研究[D];郑州大学;2015年

4 全晓静;非线性分数阶积分方程的Adomian解法[D];宁夏大学;2015年

5 黄洁;非线性分数阶Volterra积分微分方程的小波数值解法[D];宁夏大学;2015年

6 庄峤;复合介质中时间分数阶热传导正逆问题及其应用研究[D];山东大学;2015年

7 高素娟;分数阶延迟偏微分方程的紧致有限差分方法[D];山东大学;2015年

8 赵珊珊;时—空分数阶扩散方程的快速算法以及MT-TSCR-FDE的快速数值解法[D];山东大学;2015年

9 王珍;分数阶奇异边值问题的研究[D];山东师范大学;2015年

10 冯静;一类分数阶奇异脉冲边值问题正解的存在性研究[D];山东师范大学;2015年

，

本文编号：905994