抛物型方程的计算方法

1、 分类号：O241.82本科生毕业论文（设计） 题目： 一类抛物型方程的计算方法 作 者 单 位 数学与信息科学学院 作 者 姓 名 专 业 班 级 2011级数学与应用数学创新2班 指 导 教 师 论 文 完 成 时 间 二一五年四月 一类抛物型方程的数值计算方法 (数学与信息科学学院数学与应用数学专业2011级创新2班) 指导教师 摘 要: 抛物型方程数值求解常用方法有差分方法、有限元方法等。差分方法是一种对方程直接进行离散化后得到的差分计算格式，有限元方法是基于抛物型方程的变分形式给出的数值计算格式.本文首先给出抛物型方程的差分计算方法，并分析了相应差分格式的收敛性、稳定性等基本理论问题

2、.然后，给出抛物型方程的有限元计算方法及理论分析.关键词：差分方法，有限元方法，收敛性，稳定性 Numerical computation methods for a parabolic equation Yan qian (Class 2, Grade 2011, College of Mathematics and Information Science) Advisor: Nie huaAbstract: The common methods to solve parabolic equations include differential method, finite element

3、method etc. The main idea of differential method is to construct differential schemes by discretizing differential equations directly. Finite element scheme is based on the variational method of parabolic equations. In this article, we give some differential schemes for a parabolic equation and anal

4、yze their convergence and stability. Moreover, the finite element method and the corresponding theoretical analysis for parabolic equation are established. Key words: differential method, finite element method, convergence, stability1 绪 论1.1 引 言自然界里中热的传播，溶质在液体中弥散，多孔介质中渗流等随时间发展的现象和过程，都可以用抛物型方程来描述.因此，

5、抛物型方程是刻画自然界的一类很重要的方程.然而，很多的方程我们并不能求出它的精解确，或者表达式过于复杂，所以需要采用数值方法去计算它们的近似解.抛物型方程最基本的计算方法当属有限差分法，通过离散化便可得到计算格式，该方法构造简单，易于操作.但是在处理一些复杂的边值问题时计算会很复杂，因此我们需要探讨一些新的处理手段. 有限元计算方法起源于椭圆型方程的计算，它将求解椭圆型方程的解转换为求解其变分形式的解1,从而极大地丰富了偏微分方程的计算手段.正式由于其在椭圆型方程计算中的巨大优势，以及抛物型方程与椭圆型方程的密切联系，所以该方法很自然的被推广到了抛物型方程初边值问题的计算上4. 本文系统的总结

6、了一类抛物型方程的计算方法，包括有限差分法和有限元方法.并且通过数值算例给出了两类方法的一个比较.为此,本文需要先给出一些基本的分析知识作为研究该问题的基础6,7，下来就给出了抛物型方程的变分形式，这个是构造有限元计算格式的基础，在此基础上，给出了有限元计算格式并讨论了其收敛性和稳定性.1.2 准备知识 抛物型偏微分方程是一类典型的发展方程，其一般形式如下： （1.1.1） 其中是空间自变量和时间的未知函数，是关于空间变量的线性椭圆型微分算子，即 其系数，且在方程（1.1.1）的定义域中满足椭圆性条件 （1.1.2）当是非线性椭圆型微分算子或者是的非线性函数时，则称相应的抛物型方程为非线性的.

7、 下面给出抛物型方程的定解条件：初值条件，不妨设初始时刻,则 （1.1.3）第一类边值条件： （1.1.4）第二类边值条件： （1.1.5）第三类边值条件： （1.1.6）其中，为.2，有限差分法 本章将给出抛物型方程最基本的计算方法有限差分法。我们以一维热传导方程为例，给出其差分格式并讨论其收敛性，稳定性等基本问题.本章内容主要引用文献1.用差分法计算抛物型方程的初边值问题时，可以先考虑在区域上引入空间网格，例如在直角坐标系中采用平行于坐标轴的等距离直线族形成的矩形网格，其次，将定义在上的函数替换成定义在空间网格节点集上的离散函数;然后，用适当的差分格式将微分算子替换成差分算子，这一过程称为

8、半离散化.对由半离散化得到的常微分方程初值问题，再进一步对时间离散化，选用适当的求解常微分方程初值问题的数值方法，就得到求解抛物型方程的初边值问题的全离散化格式.接下来，将按照这一处理思路对热传导方程的差分计算格式进行探讨.2.1 差分格式 考虑一维热传导方程： ，0< （2.1.1） 其中是正常数，连续。下面给出两类定解条件：第一，初值问题：求可微函数，满足（1.1.1）和初始条件： （2.1.2）第二，初边值问题：求可微函数（1.1.1）和初始条件： （2.1.3） 以及边值条件 （2.1.4） 现在考虑边值问题（1.2.1），（1.2.3），（1.2.4）的差分格式.取步长空间和时

9、间步长，其中都是自然数.用两族平行直线和将矩形域分割成矩形网格，网格节点为.以表示网格内点集合，即位于开矩形的网点集合；表示所有位于闭矩形的网点集合；是网格界点集合. 其次，用表示定义在网点上的函数，用适当的差商代替方程（1.2.1）中相应的偏微商，便得到以下几种最简单的差分格式.2.1.1 向前差分格式 考虑 （2.1.5）其中 以表示网比，将（2.2.5）整理成易于计算的形式，使得第层值，即上标为在等式右边，第层值在等式左边，则可得到 （2.1.6）这样的话，又（2.1.6）取，利用初值条件和边值条件可计算出.再将的值带入计算，从而就可逐次迭代计算出所以的，并且视其为精确解的近似，由于第（

10、）层的值通过第层值明显表示，无需求解线性代数方程组，如此差分格式称为显示格式. 下来给出这种计算格式的误差分析：记 显然截断误差 （2.1.7）2.1.2向后差分格式 考虑 (2.1.8） 将上式改写为 (2.1.9)显然，第（）层的值不能用第层值明显表示，而是由线性代数方程组（2.1.9）确定，这样的差分格式称为隐格式.令 则截断误差为 （2.1.10） 此外，还有六点差分格式以及Richardson格式，具体可以参见文献1，都是简单的抛物型方程差分格式.2.2 差分格式的稳定性与收敛性 差分格式的稳定性概念见文献1，此处本文只给出相关的稳定性定理及实例分析.2.2.1判别稳定性的直接估计法

11、（矩阵法）命题11（必要条件）以表示矩阵的谱半径，则差分格式稳定性的必要条件是存在与无关的常数使 （2.2.1）命题21（充分条件）若是正规矩阵，及和它的共轭转置成绩可交换：，则（2.2.1）也是差分格式稳定的充分条件.推论1 若是对称矩阵，是矩阵的实系数有理函数：,则差分格式稳定的充要条件是，其中是的特征值。（只需注意是实数和矩阵的四则运算） 下面引出两个例子，来具体分析有限差分计算的稳定性判定：例11 对向前差分格式（以下设（2.2.4）中的），则，为使或当且仅当，从而.所以向前差分格式当时稳定，当时不稳定.2.2.2收敛性与敛速估计 如最简差分格式，考虑热传导方程的初边值问题： (2.2

12、.2）相应的差分格式为 （2.2.3）其向量形式如，其中为增长矩阵.那么差分逼近的截断误差 (2.2.4) 是上的任一充分光滑函数,称差分算子是边值问题（2.2.2）的相容逼近，如果相容条件 （2.2.5）成立，其中是分量为的向量，是中的范数. 先对差分解作出某种估计:将（2.2.3）的解分解为 其中满足零初值和非齐右端方程： 而满足非零初值和齐右端方程： 其中，依次为以，为分量的向量。若差分格式按初值稳定，则亦按右端稳定，于是有常数，使 从而 （2.2.6） 现在估计差分解的误差:设是热传导方程（2.2.2）的解，是差分方程（2.2.3）的解，误差 那么 ，即误差满足差分方程： ，其向量形式

13、为 ，这里依次为以为分量的向量.由估计式（2.2.5）得 （2.2.7）若相容条件（2.2.4）成立，则 ，其中表示以为分量的向量. 证明了如下：定理11：若差分方程满足相容条件，且按初值稳定，则差分解收敛到热传导方程的解，且有误差估计式（2.2.7）.推论11：当网比时，向前差分格式的解有收敛阶。对任何网比，向后差分格式的解有收敛阶，六点对称格式的解有收敛阶.3，有限元计算 有限元计算方法产生于椭圆型方程的计算8，其优越的计算性能使得很多学者开始探索将其用于发展方程的计算之中，文献3给出了这方面的具体研究.本章给出热传导方程的有限元计算格式,首先给出有限元计算的基本理论，之后建立热传导方程的

14、变分形式，从而在此基础上给出有限元计算格式.3.1 基本理论 本节先给出有限元计算的基本数学理论，包括索伯列夫空间初步和初边值条件下解的存在性与正则性1,6,7,9.3.1.1 Sobolev空间构造如下 Sobolev空间： 赋范数 （3.1.1） (3.1.2)定理3.17 上述赋范数的Sobolev空间是Banach空间. 特别的，当记为，引入下述内积 （3.1.3）定理3.27 是Hilbert空间. 另外，用记号表示映射族其中任一关于按空间的度量是连续的.类似的记号还有.3.1.2 解的存在性，正则性 先给出如下定理9，假定 为整数，则初边值问题（2.2.1)-(2.2.4)存在唯一

15、的解满足 和估计式 . （3.1.4 ）3.2 有限元计算格式 本节讨论初边值问题（1.1.1)-(1.1.4)的有限元近似.首先，给出其变分形式.记设用函数与（1.1.1)-(1.1.4)的两端做内积，有利用格林公式， 因，故上式右端第二项为零，引进双线性泛函 那么可得 ， （3.1.5）则称（3.1.5）为问题（1.1.1)-(1.1.4)的变分形式. 有限元方法的第一步，就是将求解区域刨分成有限个互不重叠的子区域，称其为单元.用表示刨分中单元的最大直径，记相应的刨分为 ，其代表了一个刨分族.以和表示中单元的外接圆和内切圆的直径.如果存在不依赖于的常数，使得 ，则称刨分族是正则的. 第二步

16、是构造的一族有限维子空间，要求它具有如下逼近性质：对于某一整数，有 对任意. (3.1.6)一般，是通过在刨分上作分片多项式差值的方法去构造的.由索伯列夫空间插值理论7，当刨分族为正则时，由上所有属于的，次数的分片多项式组成的子空间满足条件. 对于给定有限元空间后，初边值问题（1.1.1)-(1.1.4)的有限元近似定义为：求映射它满足 （3.1.7）其中是函数的某个近似.这里，可以看出（3.1.7）是变分问题（3.1.5）的一个近似. 设是空间的一个基底，则近似问题（3.1.5）又可以表示为：求解函数表达式 即确定其中系数，使得 ， （3.1.8 ） ，其中的系数. 所以，（3.1.8 ）是

17、以为未知函数的一个一阶常微分方程组.由于此处时间t仍然是一个连续变量，所以说（3.1.8 ）是问题（1.1.1)-(1.1.4)的一个半离散格式. 引入一些记号： 其中是一个Gram矩阵，它是非退化的，从而可以将（3.1.8）改写为矩阵形式，即 （3.1.9）由常微分方程基本理论可知，初值问题（3.1.9）对于任意和存在唯一的解，从而近似问题存在唯一解.3.3 收敛性分析和误差估计 本节介绍抛物问题有限元的理论分析，将通过能量估计证明有限元法的收敛性和近似解的误差估计. 首先给出方程 （3.1.10） 利用上一节的有限元空间，可以得到（3.1.10）的一个有限元近似:求，使得 （3.1.11）

18、其中是上节中所定义的的双线性泛函.因为正定，通过Lax-Milgram定理6可以得到，有限元方程（3.1.10）存在唯一的解. 下面，本文将不加证明的给出关于抛物型方程有限元计算方法的收敛性分析和误差估计，具体证明过程可以参考文献3,4. 定理1 假定空间具有逼近性质(3.1.6)，边值问题（3.1.10）的解，则由（3.1.11）所定义的近似解是收敛的，并且满足 （3.1.12） 定理2 （Gronwall不等式)设于连续且满足 （3.1.13）其中则 （3.1.14）对于半离散非齐次方程（3.1.7），此处将给出其解所满足的一个先验估计式. 定理3 半离散方程（3.1.7）的解满足 （3.

19、1.15）估计式（3.1.15）给出了半离散问题（3.1.7）的适定性.特别的，当时，由上面定理可知：即半离散齐次方程的解在范数意义下是稳定的. 下面给出非齐次半离散问题（3.1.7）解的误差估计：定理4 假定空间具有逼近性质(3.1.6)，并且近似初值满足 （3.1.16）则半离散问题（3.1.7）的解满足 （3.1.17）至此，本节就给出了抛物型问题有限元法的收敛性和误差估计，也就意味着本论文的目标，抛物型方程（热传导）的有限元计算方法完整的给出.4，总 结 将有限元方法推广到抛物型方程，可以丰富这类方程的计算手段，同时也可以有效的解决传统的差分格式存在的缺陷.近年来，自适应有限元计算方法

20、成为这一领域的研究前沿9-12,它在传统的有限元方法的基础上增加了“智能”的因素，使得该方法可以更好的刻画客观世界. 本文详细的给出了抛物型方程的简单差分计算格式和有限元计算格式.通过构造计算格式，得出数值计算的矩阵方程，以及给出了收敛性，稳定性，误差估计等一些基本的理论，比较全面的介绍了抛物型方程的一些基本计算方法.更为详尽的资料可以参考文献4. 当然，本文只给出了抛物型方程计算的半离散格式，并没有对时间变量t进行离散化，如何给出抛物型方程的全离散格式，以及探讨抛物型方程的其他计算方法，包括偏微分方程的应用，例如现在有学者通过凸几何的相关内容如Brunn-Minkowski 不等式来研究偏微

21、分方程的凸性等12,14，这些都是自己今后学习的方向. 参考文献1 李荣华.微分方程数值解法M.北京：高等教育出版社,20022 李志平.偏微分方程数值解讲义M.北京：北京大学出版社,20093 黄明游.发展方程的数值解法M.北京:科学出版社,20044 黄明游等译.抛物问题的Galerkin有限元法M.吉林大学出版社，1986.5 黄明游.发展方程的有限元方法M.上海:上海科学技术出版社，1988.6 张恭庆.泛函分析讲义M.北京：北京大学出版社，2006.7 王明新.索伯列夫空间M.北京：高等教育出版社，2010.8 P.G.Ciarlet.The Finite Element Method for Elliptic ProblemsM.Amsterdam:North-Holland Public Company,1978.9 Peano A G , Pasini A. Adaptive approximati