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

著录项

作者
张绪浩;
展开▼
作者单位

山东大学;

展开▼
授予单位山东大学;
学科计算数学
授予学位博士
导师姓名羊丹平,王宏;
年度 2016
页码
总页数
原文格式 PDF
正文语种中文
中图分类偏微分方程的数值解法;
关键词
分数阶偏微分方程; 数值计算; 保高精度谱; 近场动力学;

相似文献

中文文献
外文文献
专利

1. 时间空间分数阶对流扩散方程的谱Galerkin方法 [J] . 辛佳 ,黄健飞 ,赵维加 . 青岛大学学报（自然科学版） . 2013,第003期
2. 谱-Galerkin 方法求解分数阶偏积分微分方程 [J] . 李物兰 ,白宝钢 ,李胜军 . 湖南师范大学自然科学学报 . 2013,第006期
3. 一类空间分数阶偏微分方程的快速迭代方法探究 [J] . 张旻 . 兰州教育学院学报 . 2018,第010期
4. 一类带变系数的空间分数阶偏微分方程的Chebyshev拟谱分法 [J] . 杨银 . 工程数学学报 . 2014,第005期
5. 求解非线性时间分数阶Klein-Gordon方程的谱配置方法 [J] . 周琴 ,杨银 . 吉林大学学报（理学版） . 2018,第002期
6. 用快速自适应组合网格方法求解二阶椭圆型偏微分方程 [C] . 彭志健 . 第三届全国有限元会议 . 1992
7. 几类分数阶偏微分方程及近场动力学模型的高效数值算法及应用研究 [A] . 刘欢 . 2019

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

目录

摘要

著录项

相似文献

相关主题

期刊订阅