分数阶偏微分方程交替方向有限元误差分析

摘要

分数阶偏微分方程是一类很重要的微分方程，它源于许多科学领域，例如带记忆的热传导、多孔粘弹性介质的压缩、原子反应、动力学、生命科学、药理学、病理学、牛顿流体学等.在这些实际问题中，很多数学家以及应用领域中的研究者正在尝试着用分数阶微分方程来建模.由于分数阶导数是拟微分算子，它具有保记忆性（非局部性），因此在对实际问题进行优美刻画的同时，也给数值计算带来了相当大的困难，特别是对于高维问题的数值计算.
　　交替方向有限元方法是一类计算高维偏微分方程的有效数值方法，可以把多维问题降维，分解成连续解几个一维问题.因此它既具备了交替方向存储量少、计算量低又有有限元高精度的特点.本文旨在用交替方向有限元方法研究二维的分数阶偏微分方程的数值解，设计了这类方程稳定的有效的数值格式，建立了所给格式的误差估计.本文的工作主要包括以下三个部分:
　　第一部分利用交替方向隐式有限差分法研究二维的分数阶发展方程的数值解.时间方向采用向后欧拉方法，空间方向使用二阶差商，积分项用一阶卷积求积公式逼近，得到全离散格式.然后用离散的能量方法证明了全离散格式的稳定性和收敛性.最后用数值结果和数值模拟验证了理论分析的收敛阶和所给格式的有效性.该方法能够有效地减少由全局时间依赖性所引起的对存储量的要求，从而可以计算长时间的解.
　　第二部分讨论了利用交替方向隐式Galerkin有限元方法研究二维的分数阶发展方程的数值解，空间方向采用有限元方法，时间方向使用向后欧拉格式、Grank-Nicolson格式，积分项运用相应的卷积求积公式逼近，得到全离散格式.严格证明了全离散格式的稳定性和误差估计，并用数值例子检验理论分析的正确性.该方法在空间方向上能达到任意阶精度，也能够有效地减少由全局时间依赖性所引起的对存储量的要求.
　　第三部分利用交替方向隐式Galerkin有限元方法研究二维的时间分数阶扩散-波动方程的数值解.时间方向使用Grank-Nicolson格式，空间方向采用Galerkin有限元，得到全离散格式.证明了全离散格式的稳定性和收敛性，数值例子验证其结论的正确性.

目录

摘要

1 绪论

1．1 研究背景

1．2 分数阶微积分的定义和性质

1．3 本文研究的分数阶方程

1．4 研究动机和主要工作

2 分数阶发展方程交替方向有限差分法误差分析

2．1 预备知识

2．2 交替方向有限差分格式的构造

2．3 交替方向有限差分格式稳定性分析

2．4 交替方向有限差分格式收敛性分析

2．5 数值算例

3 分数阶发展方程交替方向有限元误差分析

3．1 预备知识

3．2 交替方向有限元格式的推导

3．3 交替方向有限元格式的收敛性和稳定性

3．4 数值算例

4 分数阶发展方程Crank-Nicolson交替方向有限元误差分析

4．1 预备知识

4．2 Crank-Nicolson交替方向有限元格式的推导

4．3 Crank-Nicolson交替方向有限元格式稳定性分析

4．4 Crank-Nicolson交替方向有限元格式收敛性分析

4．5 数值算例

5 分数阶扩散-波动方程交替方向有限元误差分析

5．1 预备知识

5．2 Crank-Nicolson交替方向有限元格式的推导

5．3 交替方向Galerkin有限元格式稳定性分析

5．4 交替方向Galerkin有限元格式收敛性分析

5．5 数值算例

6 总结和未来工作展望

参考文献

攻读博士学位期间发表的学术论文

致谢

声明