抛物型方程有限差分方法的应用 - 报告

2015 年 秋 季学期研究生课程考核

(读书报告、研究报告)

考核科目: 偏微分方程数值解法 学生所在院(系): 理学院数学系学生所在学科: 数学学 生 姓 名: Hiter学 号: 1XS012000学 生 类 别: 考核结果阅卷人

抛物型方程有限差分方法的应用

摘要

抛物型偏微分方程是一类比较重要的偏微分方程。热传导方程是最简单的一种抛物型方程。热传导方程研究的是热传导过程的一个简单数学模型。根据热量守恒定律和傅里叶热传导实验定律可以导出热传导方程。在本篇论文中，将先详述抛物型偏微分方程的有限差分法的相关知识，然后给出抛物型方程的两个具体的应用实例。

关键字：抛物型方程，差分格式，应用

Abstract

Parabolic partial differential equation is a kind of important partial differential equation. The heat conduction equation is one of the simplest parabolic equations. The heat conduction equation is a simple mathematical model of the heat conduction process. Heat conduction equation is derived based on the law of conservation of heat and Friyege's law of conduction. In this thesis, we first give a detailed knowledge of the finite difference method for parabolic partial differential equations, and then give two specific examples of the application of the parabolic equation.

Keywords: parabolic equation, difference scheme, application

0 前言

抛物型方程是偏微分方程中的三大方程(另两种为双曲型方程和椭圆型方程)之一，如何去研究抛物型方程的性质在《偏微分方程数值解法》的课程中占有很大的比例。如果要研究抛物型方程，我们一般从以下几个方面来研究，分别是：定界问题、格林函数、极值原理、解的正则性、抛物方程、拟线性蜕化和反应扩散方程。但是由于所学知识有限，所以我们在此只简单的介绍抛物型方程的有限差分法并给出两个应用实例。

1 抛物型方程有限差分法

1.1 简单差分法

考虑一维模型热传导方程

， (1.1)

其中为常数。是给定的连续函数。(1.1)的定解问题分两类：

第一，初值问题(Cauchy 问题)：求足够光滑的函数，满足方程(1.1)和初始条件：

， (1.2)

第二，初边值问题(也称混合问题)：求足够光滑的函数，满足方程(1.1)和初始条件：

,

及边值条件

，

假定和在相应的区域光滑，并且于，两点满足相容条件,则上述问题有唯一的充分光滑的解。

现在考虑边值问题(1.1)，(1.3)的差分逼近。取为空间步长，为时间步长，其中，是自然数，

, ； ,

将矩形域分割成矩形网格。其中 表示网格节点；

表示网格内点(位于开矩形中的网格节点)的集合；

表示位于闭矩形中的网格节点的集合；

表示-网格边界点的集合；

表示定义在网点处的待求近似解，，。

注意到在节点处的微商和差商之间的下列关系()：

可得到以下几种最简差分格式

(一)向前差分格式

，

，==0

其中，。取为网比，则进一步有

=+++

此差分格式是按层计算：首先，令，得到

=+++

于是，利用初值和边值==0，可算出第一层的，。再由取，可利用和==0算出，。如此下去，即可逐层算出所有(，)。

由于第层值可以通过第层值直接得到，如此的格式称为显格式。并视为的近似值。

若记

，，

则显格式可写成向量形式

其中

若记

那么截断误差为

== (1.5)

其中是矩形，中某一点。

事实上，+

=+

=

==。

这里

故，从而

(二)向后差