抛物型方程的计算方法

《抛物型方程的计算方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

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

2、抛物型方程的有限元计算方法及理论分析.关键词：差分方法，有限元方法，收敛性，稳定性NumericalcomputationmethodsforaparabolicequationYanqian(Class2,Grade2011,CollegeofMathematicsandInformationScience)Advisor:NiehuaAbstract:Thecommonmethodstosolveparabolicequationsincludedifferentialmethod,finiteelementmethodetc.Them

3、ainideaofdifferentialmethodistoconstructdifferentialschemesbydiscretizingdifferentialequationsdirectly.Finiteelementschemeisbasedonthevariationalmethodofparabolicequations.Inthisarticle,wegivesomedifferentialschemesforaparabolicequationandanalyzetheirconvergenceandstabilit

4、y.Moreover,thefiniteelementmethodandthecorrespondingtheoreticalanalysisforparabolicequationareestablished.Keywords:differentialmethod,finiteelementmethod,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、间网格节点集上的离散函数;然后，用适当的差分格式将微分算子替换成差分算子，这一过程称为半离散化.对由半离散化得到的常微分方程初值问题，再进一步对时间离散化，选用适当的求解常微分方