分数阶偏微分方程的几类有限元方法研究应用数学专业论文.docxVIP

分数阶偏微分方程的几类有限元方法研究摘要 分数阶偏微分方程的几类有限元方法研究 摘要 由于分数阶导数能够比整数阶导数更准确地描述具有记忆和遗传性质的 材料与物理过程，分数阶微分方程在诸多领域得到了广泛应用和深入研究．但 是分数阶微分方程的解析解通常很难求出，即便能求出解析解，大多数的解也 都含有无穷级数或一些难以计算的特殊函数，所以人们更加关注分数阶微分 方程的数值解法．目前在对分数阶微分方程所提出的各种数值方法中，有限元 方法由于区域适应性强，网格剖分灵活，并且简单通用，对解的光滑性要求不 高，因而更受关注．本文研究了分数阶偏微分方程的几类有限元方法，其中重 点研究了时间分数阶微分方程的降阶有限元方法和空间分数阶微分方程及方 程组的时空有限元方法．所做的工作可以分为如下三部分． 第一部分即本文的第三章研究了时间分数阶Cable方程的标准有限元方法． 采用传统有限元方法的思想，在时间方向上采用有限差分离散，空间方向上利 用有限元近似，得到了分数阶Cable方程的Crank—Nicolson型全离散格式．我们 对时间分数阶导数的差分离散给出了不同于L1算法的新的系数，所得到的全 离散格式关于时间层数几具有更好的继承性．文中还详细给出了全离散格式的 稳定性分析和最优阶Lz模误差估计结果，而且分别展示了空间一维和二维的 数值例子，验证了理论分析的正确性． 第二部分利用基于特征投影分解(POD)理论的降阶有限元方法分别分析 了空N--维的时间分数阶扩散方程、时间分数阶Tricomi型方程以及时间分数 阶Sob01ev方程．所采用的降阶方法是：先利用一般有限元格式计算出在前面 很短时间段内的有限元解，将其作为瞬像；然后从瞬像中找到在最小平方意义 下的最优POD基，这些基函数的个数远远小于一般有限元方法的基函数的个 数，将它们张成的空间作为降阶有限元空间；最后利用降阶有限元格式求得近 似解．在时间分数阶微分方程的有限元方法中，要求t礼时的解，需要将前面所 有满足tt。的时间层上的解存储并叠加起来．由于我们在每个时间层上降低 了解空间的自由度，整体的存储量和计算量就得到了大幅度降低，从而减轻了 万方数据 分数阶微分算子的非局部性带来的计算负担．在文中，我们给出了分数阶微分 分数阶微分算子的非局部性带来的计算负担．在文中，我们给出了分数阶微分 方程的全离散有限元格式和降阶有限元格式并分别进行了稳定性和收敛性分 析，得到了有限元解和降阶有限元解的Lz模最优阶误差估计结果．数值算例也 表明，基于POD理论的降阶有限元方法能够保证在精度上不低于传统有限元 方法的情况下极大地降低了存储量，提高了计算效率．本部分的内容安排在论 文的第四章至第六章，在第五章中我们还具体给出了降阶有限元方法的算法 流程． 本文的第三部分研究了空间分数阶微分方程(组)的时空有限元方法．在第 七章中，利用时间间断而空间连续的间断时空有限元方法研究了半线性空间 分数阶扩散方程．对近似变分问题，通过适当选取试探函数得到了可以逐时间 层推进求解的时空有限元全离散格式．在理论分析中，利用Radau积分公式将 近似解表示为基于Radau积分点的Lagrange插值多项式的线性组合的形式，从 而在不对时空网格施加任何限制条件的情况下给出了弱解的存在唯一性证明． 通过引入椭圆投影算子，利用尼采(Nitche)技巧给出了分数阶椭圆投影的Lz模 估计，进而详细导出了时空有限元解的最优阶L。。(L2)模误差估计结果．在第八 章，利用时间间断的时空有限元方法构造了非线性空间分数阶反应扩散方程组 的全离散格式并给出了格式的适定性分析和误差估计，从而将间断时空有限元 方法进一步推广到了分数阶方程组．理论分析和数值模拟结果均显示，间断时 空有限元方法能够在时间和空间两个方向上同时达到高阶精度，对空间分数 阶微分方程(组)仍适用．这为今后研究更复杂的分数阶方程(组)奠定了基础． 关键词：分数阶偏微分方程；降阶有限元方法；POD理论；间断时空有限元方法 稳定性；误差估计 万方数据 SOME SOME TYPES OF FINITE ELEMENT METHoDS FOR FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS AB STRACT Since the fractional derivatives can describe the materials with memory)inheritance characters well physical process more precisely，fractional differential equations have been applied extensively and researched deeply i