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

分类号Q!幽 密级 公开 学校代码 !Q鲤2 学号至Q!QQ!!QQ!Q鱼 分数阶偏微分方程交替方向有限元误差分 析 directionfiniteelement Alternatingimplicit error for thefractional analysis partial differential equations 研 究 生 姓 名： 黎丽梅 指导教师姓名、职称： 徐大 教授 学 科 专 业： 计 算 数 学 研 究 方 向： 偏微分方程数值解 湖南师范大学学位评定委员会办公室 二零一三年十二月 摘 要 分数阶偏微分方程是一类很重要的微分方程，它源于许多科学 领域，例如带记忆的热传导、多孔粘弹性介质的压缩、原子反应、动 力学、生命科学、药理学、病理学、牛顿流体学等．在这些实际问题 中，很多数学家以及应用领域中的研究者正在尝试着用分数阶微分 方程来建模．由于分数阶导数是拟微分算子，它具有保记忆性(非局 部性)，因此在对实际问题进行优美刻画的同时，也给数值计算带来 了相当大的困难，特别是对于高维问题的数值计算． 交替方向有限元方法是一类计算高维偏微分方程的有效数值方 法，可以把多维问题降维，分解成连续解几个一维问题．因此它既 具备了交替方向存储量少、计算量低又有有限元高精度的特息本 文旨在用交替方向有限元方法研究二维的分数阶偏微分方程的数值 解，设计了这类方程稳定的有效的数值格式，建立了所给格式的误差 估计．本文的工作主要包括以下三个部分： 第一部分利用交替方向隐式有限差分法研究二维的分数阶发展 方程的数值解．时间方向采用向后欧拉方法，空间方向使用二阶差 商，积分项用一阶卷积求积公式逼近，得到全离散格式．然后用离散 的能量方法证明了全离散格式的稳定性和收敛性．最后用数值结果 和数值模拟验证了理论分析的收敛阶和所给格式的有效性．该方法 能够有效地减少由全局时间依赖性所引起的对存储量的要求，从而 可以计算长时间的解． 第二部分讨论了利用交替方向隐式Galerkin有限元方法研究二 维的分数阶发展方程的数值解，空间方向采用有限元方法，时间方 向使用向后欧拉格式、Grank-Nicolson格式，积分项运用相应的卷积 求积公式逼近，得到全离散格式．严格证明了全离散格式的稳定性 和误差估计，并用数值例子检验理论分析的正确性．该方法在空间方 向上能达到任意阶精度，也能够有效地减少由全局时间依赖性所引 起的对存储量的要求． 第三部分利用交替方向隐式Galerkin有限元方法研究二维的时 间分数阶扩散一波动方程的数值解．时间方向使用Grank-Nicolson格 式，空间方向采用Galerkin有限元，得到全离散格式．证明了全离散 格式的稳定性和收敛性，数值例子验证其结论的正确性． 关键词： 分数阶偏微分方程；有限元方法；交替方向法；有限差分法；偏 积分微分方程；卷积求积；弱奇异核；向后欧拉方法． ABSTRACT is classof Thefractionaldifferential a veryimportant partial equation differential