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

分数阶偏微分方程保高精度谱Galerkin方法以及近场动力学模型快速配置方法研究

张绪浩  
【摘要】:分数阶微积分概念历史悠久,最早源于1695年9月L' Hospital写给莱布尼茨的信件中。在分数阶微积分被提出至今300多年时间中,由于在物理和力学等学科并未获得广泛的关注与应用,而仅仅作为数学领域中的纯理论问题被诸多数学学者研究,这其中包括Euler, Lacroix, Abel, Liouville, Riemann等。随着对复杂物理现象认识程度的加深和计算机模拟能力的提高,力学与工程问题的分数阶导数建模越来越引起人们的重视。尤其对于扩散现象,研究工作者发现越来越多的扩散现象不满足Fick定律,这样的扩散过程称为反常扩散过程。在描述这些复杂系统时,由于反常扩散所具有的历史依赖与全域相关的特性恰好可以由分数阶导数来表示,因此较于整数阶动力学方程,分数阶动力学方程更能有效的描述([18,50,57,5])。无论在理论分析还是数值计算方面,新的分数阶动力学方程都为数学工作者提出了新的挑战。在数值计算方面,现在已经有很多数值求解方法,例如有限差分法([14,42]),有限体积法([91]),有限元法([60,25]),谱方法([35,94])等。这些计算方法已经被广泛的应用于反常扩散的数值模拟中。但是这些方法的误差分析均有很强的正则性假设。通过本文第二章中的反例,我们可以看到即使分数阶方程的扩散系数和右端均充分光滑,我们依然不能保证解的正则性,这是分数阶动力学方程区分于整数阶动力学方程的一个很重要的特性(在整数阶动力学方程情况下,根据方程正则性理论,方程系数和右端的光滑性可以保证方程解的正则性)。通过这一点,我们知道上述数值方法误差分析的假设条件缺少理论支撑,另外在方程解的正则性不强的情况下,即使方程的系数和右端光滑,我们采用高阶差分,高次有限元以及标准谱方法均不能达到很好的误差收敛情况。在计算效率方面,分数阶动力学方程离散得出的系数矩阵通常为满阵,如果我们假设矩阵的阶数或问题的规模为N,则系数矩阵的存储量为O(N2),如果用常用的直接方法求解分数阶动力学方程离散得出的线性系统,则计算复杂度为O(N3)。因此,一直以来,尤其是求解大规模或者多维问题时,分数阶动力学模拟是很费时间的。为解决这一问题,王宏等著名学者通过分析系数矩阵的代数结构,运用快速傅立叶变换,成功的将系数矩阵存储量降为O(N),将求解线性方程组的计算量降为每一步Krylov子空间迭代O(N log N)的计算量([79,80])。固体材料和结构的破坏问题一直是力学研究的经典问题,也是机械、航空航天、土木、水利和化工等领域关注的重点。在近场动力学理论提出之前,随着断裂力学、损伤力学等学科的发展和计算机软件硬件水平的提高,研究者提出了各种不同的力学模型和数值方法来模拟固体材料和渐进破坏的全过程。这些模型均是建立在连续介质假设之上的,他们假设介质所有的内部力均为接触力,最终控制方程绝大多数由偏微分方程所描述。传统的有限元法和有限差分法同样建立在连续介质假设的思想上,在模拟时必须明确知道断裂的位置与尺寸,这在很多现实应用中很难实现,另外随着断裂的发展,传统的有限元法或者有限差分法必须重新划分网格,有很强的网格依赖性([33])。随着不连续有限元方法的发展,在模拟固体材料断裂及发展问题上取得了一些进步,但是在模拟高维复杂断裂系统时仍有很强的局限性。为了克服连续介质力学假定与固体材料不连续这一基本矛盾,2000年,Silling基于非局部作用建模,提出了近场动力学模型,它是用积分思想表述的积分方程([68])。这一模型不在基于连续介质假设和求解微分方程来模拟破坏问题,而是将固体看成由由一些包含所有物质信息的带质量的物质点组成,点与点之间存在着相互作用,随着点与点之间距离的增加,这个作用力在减弱,因此通常人们选取一个点的δ邻域为其作用力的影响域。在该理论框架下,不连续现象自然产生,同时这一理论突破了分子动力学在计算尺度上的局限,在宏、细、微观尺度均可表现出较高的求解精度。在近场动力学提出之后,很多数值方法例如无网格方法([64,63,70])、有限元方法([15]),基于积分的有限差分方法([77])等被提出求解近场动力学模型。在有限元情况下,已经被证明数值解满足最优误差估计。然而这些方法都有一个共同的特点,特别在求解多维问题时,由于离散所得到的系数矩阵为稠密矩阵或者满阵(这取决于影响域δ的大小)。因此,类似分数阶动力学方程,系数矩阵的存储量为O(N2),求解最后线性系统的计算复杂度为O(N3)。另外如果用有限元法求解近场动力学模型,每一个系数矩阵的元素均需要计算2d次重积分,其中d是维数,但由于积分核含有奇性,则计算这一积分是很耗时的。王宏等学者同样根据矩阵的代数结构与快速傅立叶变换,成功的将矩阵存储量降为O(N),将求解最后线性系统的计算量降为每Krylov子空间迭代O(N log N)的计算量([81,82])。基于以上考虑,我们分别研究了分数阶扩散方程的保高精度谱Galerkin方法和近场动力学模型的快速配置算法.本文的安排如下:在第一章中,我们给出了在本文剩余部分用到的一些基本概念,包括分数阶导数的Riemann-Liouville导数定义和Caputo定义及其一些基本性质,另外我们给出了一些特殊矩阵的定义及其一些性质。在第二章中我们给出了一种分数阶扩散方程的保高精度谱Galerkin方法,这种方法可以保证在方程系数和右端都充分光滑的条件下,即使真解没有足够的光滑性,我们也可以保证解的高精度。并且数值解比标准谱Galerkin方法得到的数值解精度要好,因为在真解没有足够光滑性的条件下.标准谱Galerkin方法并不能达到高精度。我们同时证明了该方法的误差估计。这一章中给出的算例说明了这一保高精度谱方法的有效性。在第三章中我们提出了求解二维近场动力学模型的快速配置方法。在这一章中,我们仔细分析了由配置法离散得出的系数矩阵,经过分析我们得出系数矩阵与任何向量的乘积可以由三个block-Toeplitz-Toeplitz-block (BT-TB)矩阵与向量的乘积得到,所以系数矩阵与向量的乘法的计算复杂度为O(N log N),如果用Krylov子空间迭代法求解该线性系统,则每一步迭代的计算量可以由O(N2)降为O(NlogN).同时从本文也可以得到该矩阵的计算机存储量可以由O(N2)降为O(N).这一章中给出的算例说明了快速配置方法的有效性。在第四章中针对在近场动力学模型中积分核函数奇性大,运用Krylov子空间迭代求解由于配置法离散得出的线性系统迭代次数比较多的情况,我们提出了两种预条件矩阵。第一种预条件是block-Cireulant-Toeplitz-block (BCTB)型的,第二种预条件是block-Circulant-Circulant-block (BCCB)型的。这两种预条件对降低求解线性系统迭代次数是有效的,并且通过这一章给出的例子我们可以看到第二种预条件因为计算预条件矩阵的逆比较快速,所以计算时间会更快。在第五章中我们运用加罚的思想提出了求解一般凸区域非局部扩散模型的快速配置方法。这种方法是通过加罚将原来凸区域上的问题扩展为包含该凸区域的矩形的问题,经过配置法离散,我们可以看到系数矩阵是一个BTTB矩阵与一个对角矩阵的和,通过这种矩阵结构,我们将系数矩阵的存储由O(N2)降为O(N),将每一步Krylov子空间的计算量由O(N2)降为O(N log N).数值算例说明了这种方法的有效性。


知网文化
【相似文献】
中国期刊全文数据库 前20条
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期
11 刘勇;谢勇;;分数阶FitzHugh-Nagumo模型神经元的动力学特性及其同步[J];物理学报;2010年03期
12 刘荣花;;分数阶积分和微分函数的图像k维数研究[J];中国校外教育;2010年14期
13 曹鹤飞;张若洵;;分数阶混沌系统参数调制数字通信及电路仿真[J];河北师范大学学报(自然科学版);2011年04期
14 林浩;罗懋康;;基于分数阶微(积)分的一种降噪和白化方法[J];信息通信;2012年01期
15 刘艳芹;;一类分数阶非线性振子方程的特性研究[J];计算机工程与应用;2012年16期
16 马靖杰;;空间分数阶Edwards-Wilkinson方程的数值研究[J];山东大学学报(理学版);2012年09期
17 刘荣花;;分数阶积分和微分函数的图像k维数研究[J];中国校外教育;2010年S1期
18 辛宝贵;马军海;陈通;;分数阶双卷混沌系统的复杂性演化仿真研究[J];陕西科技大学学报(自然科学版);2011年02期
19 王妍;;多涡旋混沌同步和分数阶混沌同步[J];系统科学与数学;2011年12期
20 李宝凤;;应用Bernstein多项式求解一类分数阶微分方程[J];唐山师范学院学报;2014年02期
中国重要会议论文全文数据库 前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年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978