抛物型方程的计算方法.docVIP

分类号：O241.82 本科生毕业论文（设计） 题目： 作 者 单 位 数学与信息科学学院 作 者 姓 名 专 业 班 级 201级班 指 导 教 师 论 文 完 成 时 间 二〇一年四月 umerical computation methods for a parabolic equation Yan qian (Class 2, Grade 2011, College of Mathematics and Information Science) Advisor: Nie hua Abstract: The common methods to solve parabolic equations include differential method, finite element 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 analyze 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, stability 1 绪 论 1.1 引 言 自然界里中热的传播，溶质在液体中弥散，多孔介质中渗流等随时间发展的现象和过程，都可以用抛物型方程来描述.因此，抛物型方程是刻画自然界的一类很重要的方程.然而，很多的方程我们并不能求出它的精解确，或者表达式过于复杂，所以需要采用数值方法去计算它们的近似解. 抛物型方程最基本的计算方法当属有限差分法[]，通过离散化便可得到计算格式，该方法构造简单，易于操作.但是在处理一些复杂的边值问题时计算会很复杂，因此我们需要探讨一些新的处理手段. 有限元计算方法起源于椭圆型方程的计算，它将求解椭圆型方程的解转换为求解其变分形式的解[1],从而极大地丰富了偏微分方程的计算手段.正式由于其在椭圆型方程计算中的巨大优势，以及抛物型方程与椭圆型方程的密切联系，所以该方法很自然的被推广到了抛物型方程初边值问题的计算上[4]. 本文系统的总结了一类抛物型方程的计算方法，包括有限差分法和有限元方法.并且通过数值算例给出了两类方法的一个比较.为此,本文需要先给出一些基本的分析知识作为研究该问题的基础[6,7]，下来就给出了抛物型方程的变分形式，这个是构造有限元计算格式的基础，在此基础上，给出了有限元计算格式并讨论了其收敛性和稳定性. 1.2 准备知识 抛物型偏微分方程是一类典型的发展方程，其一般形式如下： （1.1.1） 其中是空间自变量和时间的未知函数，是关于空间变量的线性椭圆型微分算子，即 其系数，且在方程（1.1.1）的定义域中满足椭圆性条件 （1.1.2） 当是非线性椭圆型微分算子或者是的非线性函数时，则称相应的抛物型方程为非线性的. 下面给出抛物型方程的定解条件： 初值条件，不妨设初始时刻,则 （1.1.3） 第一类边值条件： （1.1.4） 第二类边值条件：