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

基于多项式插值逼近的分数阶偏微分方程高精度差分方法

孙红  
【摘要】:近几十年来,由于分数阶导数具有非局部性质,比整数阶导数更适合描述具有记忆和遗传性质的材料和过程.因此,分数阶微分方程更能准确地刻画许多自然界的现象,得到了越来越多的学者的关注.关于分数阶偏微分方程的数值解法无论对工程技术领域还是对数学本身都具有重要的价值.本文主要是在超收敛点处对时间分数阶波方程、多项时间分数阶波方程、时空分数阶Bloch-Torrey方程以及非线性时间分数阶四阶反应-扩散方程等初边值问题构造数值解法,并给出相应的理论分析.本文首先研究的是时间分数阶波方程的数值解法.对一维和二维时间分数阶波方程利用降阶法得到等价的方程组,然后利用L2-1σ公式(Alikhanov,J.Comput.Phys.280(2015),424-438)对等价方程组建立时间方向二阶精度,空间方向分别为二阶和四阶精度的有限差分格式.利用离散能量法,严格证明了格式在H1范数下的无条件稳定性和收敛性.同时还给出了三维时间分数阶波方程的差分格式.数值算例验证了格式的计算精度和有效性.其次,对多项时间分数阶波方程建立时间二阶精度的有限差分格式.利用降阶法得到等价方程组,再对方程组中的多项时间分数阶导数在其超收敛点处离散,从而对多项时间分数阶波方程分别建立时间和空间方向都为二阶精度的差分格式和时间二阶、空间四阶精度的差分格式.我们证明了两个格式是唯一可解的,且在最大模下是无条件稳定的和收敛的,收敛阶分别为O(τ2 + h2)和O(τ2 + h4).数值实验表明格式的有效性,验证了差分格式的理论分析精度.随后,讨论了一维和二维时空分数阶Bloch-Torrey方程的差分方法.利用L2-1σ公式来离散时间分数阶Caputo导数,分别应用分数阶二阶中心差分格式(C.Celik,M.Duman,J.Comput.Phys.231(2012),1743-1750.)和四阶紧算子(X.Zhao,Z.Z.Sun,Z.P.Zhao,SIAM J.Sci.Comput.36(2014),A2865-A2886.)对空间分数阶Riesz导数进行逼近,从而对一维和二维Bloch-Torrey方程构造有限差分格式.同时我们给出了分数阶二阶中心差分算子的权系数和的下界的一个估计式.利用离散能量法以及权系数的下界估计式,我们对格式的稳定性和收敛性给出了严格的理论证明.对二维问题,我们还给出了两个ADI格式来求解方程.数值算例验证了差分格式的有效性.最后一部分考虑了非线性时间分数阶四阶反应-扩散方程的数值逼近.首先利用降阶法,得到一个等价方程组,运用L2-1σ公式对时间Caputo分数阶导数进行离散,对空间整数阶导数采用二阶格式离散,进而构造一个三层线性化的有限差分格式.利用离散能量法,我们给出了格式在L2模下的无条件稳定性和收敛性的严格的理论证明,收敛阶为O(τ2 +h12+h22).对于差分格式的收敛性证明是理论分析的一个难点,我们主要应用二维网格函数空间的一个嵌入定理给出了格式的收敛性分析.数值实验验证了格式的理论分析的精度.


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 赵羽,蔡平,周敏东;分数阶Fourier变换的数值计算[J];哈尔滨工程大学学报;2002年06期
2 王德金;郑永爱;;分数阶混沌系统的延迟同步[J];动力学与控制学报;2010年04期
3 杨晨航,刘发旺;分数阶Relaxation-Oscillation方程的一种分数阶预估-校正方法[J];厦门大学学报(自然科学版);2005年06期
4 王发强;刘崇新;;分数阶临界混沌系统及电路实验的研究[J];物理学报;2006年08期
5 夏源;吴吉春;;分数阶对流——弥散方程的数值求解[J];南京大学学报(自然科学版);2007年04期
6 张隆阁;;一类参数不确定混沌系统的分数阶自适应同步[J];中国科技信息;2009年15期
7 陈世平;刘发旺;;一维分数阶渗透方程的数值模拟[J];高等学校计算数学学报;2010年04期
8 辛宝贵;陈通;刘艳芹;;一类分数阶混沌金融系统的复杂性演化研究[J];物理学报;2011年04期
9 黄睿晖;;分数阶微方程的迭代方法研究[J];长春理工大学学报;2011年06期
10 蒋晓芸,徐明瑜;分形介质分数阶反常守恒扩散模型及其解析解[J];山东大学学报(理学版);2003年05期
11 陈玉霞;高金峰;;一个新的分数阶混沌系统[J];郑州大学学报(理学版);2009年04期
12 刘勇;谢勇;;分数阶FitzHugh-Nagumo模型神经元的动力学特性及其同步[J];物理学报;2010年03期
13 刘荣花;;分数阶积分和微分函数的图像k维数研究[J];中国校外教育;2010年14期
14 曹鹤飞;张若洵;;分数阶混沌系统参数调制数字通信及电路仿真[J];河北师范大学学报(自然科学版);2011年04期
15 林浩;罗懋康;;基于分数阶微(积)分的一种降噪和白化方法[J];信息通信;2012年01期
16 刘艳芹;;一类分数阶非线性振子方程的特性研究[J];计算机工程与应用;2012年16期
17 马靖杰;;空间分数阶Edwards-Wilkinson方程的数值研究[J];山东大学学报(理学版);2012年09期
18 刘荣花;;分数阶积分和微分函数的图像k维数研究[J];中国校外教育;2010年S1期
19 辛宝贵;马军海;陈通;;分数阶双卷混沌系统的复杂性演化仿真研究[J];陕西科技大学学报(自然科学版);2011年02期
20 王妍;;多涡旋混沌同步和分数阶混沌同步[J];系统科学与数学;2011年12期
中国重要会议论文全文数据库 前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 尹学辉;基于分数阶PDE的图像结构保持型去噪算法研究[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年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978