第十章偏微分方程数值解法

1、第十章 偏微分方程数值解法偏微分方程问题，其求解十分困难。除少数特殊情况外，绝大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。§1 差分方法的基本概念1.1 几类偏微分方程的定解问题椭圆型方程：其最典型、最简单的形式是泊松（Poisson）方程特别地，当时，即为拉普拉斯（Laplace）方程，又称为调和方程 Poisson方程的第一边值问题为 其中为以为边界的有界区域，为分段光滑曲线，称为定解区域，分别为，上的已知连续函数。第二类和第三类边界条件可统一表示为 其中为边界的外法线方向。当时为第二类边界条件，时为第三类边界条件

2、。 抛物型方程：其最简单的形式为一维热传导方程方程可以有两种不同类型的定解问题：初值问题初边值问题 其中，为已知函数，且满足连接条件边界条件称为第一类边界条件。第二类和第三类边界条件为 其中，。当时，为第二类边界条件，否则称为第三类边界条件。双曲型方程:最简单形式为一阶双曲型方程物理中常见的一维振动与波动问题可用二阶波动方程描述，它是双曲型方程的典型形式。方程的初值问题为边界条件一般也有三类，最简单的初边值问题为 1.2 差分方法的基本概念差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有

3、限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：（1）选取网格； （2）对微分方程及定解条件选择差分近似，列出差分格式； （3）求解差分格式； （4）讨论差分格式解对于微分方程解的收敛性及误差估计。下面，用一个简单的例子来说明用差分方法求解偏微分方程问题的一般过程及差

4、分方法的基本概念。设有一阶双曲型方程初值问题。 （1） 选取网格： 2h -h 0 h 2h 3h首先对定解区域作网格剖分，最简单常用一种网格是用两族分别平行于轴与轴的等距直线，将分成许多小矩形区域。这些直线称为网格线，其交点称为网格点，也称为节点，和分别称作方向和方向的步长。这种网格称为矩形网格。（2） 对微分方程及定解条件选择差分近似，列出差分格式： 如果用向前差商表示一阶偏导数，即 其中。方程 在节点处可表示为 其中。由于当足够小时，在式中略去，就得到一个与方程相近似的差分方程 此处，可看作是问题的解在节点处的近似值。同初值条件 结合，就得到求问题的数值解的差分格式。式 称为差分方程的截

5、断误差。如果一个差分方程的截断误差为，则称差分方程对是阶精度，对是阶精度的。显然，截断误差的阶数越大，差分方程对微分方程的逼近越好。若网格步长趋于0时，差分方程的截断误差也趋于0，则称差分方程与相应的微分方程是相容的。这是用差分方法求解偏微分方程问题的必要条件。如果当网格步长趋于0时，差分格式的解收敛到相应微分方程定解问题的解，则称这种差分格式是收敛的。§2 椭圆型方程第一边值问题的差分解法本节以Poisson方程为基本模型讨论第一边值问题的差分方法。2.1 差分格式的建立 考虑Poisson方程的第一边值问题取分别为方向和方向的步长，如图所示，以两族平行线，将定解区域剖分成矩形网格

6、。节点的全体记为。定解区域内部的节点称为内点，记内点集为。边界与网格线的交点称为边界点，边界点全体记为。与节点沿方向或方向只差一个步长的点和称为节点的相邻节点。如果一个内点的四个相邻节点均属于，称为正则内点，正内点的全体记为，至少有一个相邻节点不属于的内点称为非正则内点，非正则内点的全体记为。问题是要求出第一边值问题在全体内点上的数值解。为简便，记，。对正则内点，由二阶中心差商公式 Poisson方程 在点处可表示为其中 为其截断误差表示式，略去，即得与方程相近似的差分方程 式中方程的个数等于正则内点的个数，而未知数则除了包含正则内点处解的近似值外，还包含一些非正则内点处的近似值，因而方程个数

7、少于未知数个数。在非正则内点处Poisson方程的差分近似不能按上式给出，需要利用边界条件得到。边界条件的处理可以有各种方案，下面介绍较简单的两种。（1）直接转移用最接近非正则内点的边界点上的值作为该点上值的近似，这就是边界条件的直接转移。例如，点为非正则内点，其最接近的边界点为点，则有 上式可以看作是用零次插值得到非正则内点处的近似值，容易求出，其截断误差为。将上式代入，方程个数即与未知数个数相等。 （2）线性插值这种方案是通过用同一条网格线上与点相邻的边界点与内点作线性插值得到非正则内点处值的近似。由点与的线性插值确定的近似值，得 其中，其截断误差为。将其与方程相近似的差分方程联立，得到方

8、程个数与未知数个数相等的方程组，求解此方程组可得Poisson方程第一边值问题的数值解。上面所给出的差分格式称为五点菱形格式，实际计算时经常取，此时五点菱形格式可化为 简记为 其中。例1 用五点菱形格式求解拉普拉斯（Laplace）方程第一边值问题其中。取。 （0,0） (1,0) (2,0) (3,0)解网格中有四个内点，均为正则内点。由五点菱形格式，得方程组 代入边界条件其解为 ，当时，对 利用点，构造的差分格式，称为五点矩形格式，简记为 其中，其截断误差为 五点菱形格式与矩形格式的截断误差均为，称它们具有二阶精度。如果用更多的点构造差分格式，其截断误差的阶数可以提高，如利用菱形格式及矩形

9、格式所涉及的所有节点构造出的九点格式就是具有四阶精度的差分格式。§3 抛物型方程的差分解法以一维热传导方程为基本模型讨论适用于抛物型方程定解问题的几种差分格式。3.1 差分格式的建立首先对平面进行网格剖分。分别取为方向与方向的步长，用两族平行直线，将平面剖分成矩形网格，节点为。为简便，记，。（一）微分方程的差分近似 在网格内点处，对分别采用向前、向后及中心差商公式一维热传导方程可分别表示为 由此得到一维热传导方程的不同差分近似 上述差分方程所用到的节点各不相同。其截断误差分别为，和。因此，它们都与一维热传导方程相容。如果将式中的用代替，则可得到又一种差分近似 差分方程用到四个节点。由

10、Taylor公式容易得出故其的截断误差为。因而不是对任意的，此差分方程都能逼近热传导方程，仅当时，才成立。综上可知，用不同的差商公式可以得到微分方程的不同的差分近似。构造差分格式的关键在于使其具有相容性、收敛性和稳定性。前面三个方程都具有相容性，而此方程则要在一定条件下才具有相容性。（二）初、边值条件的处理 为用差分方法求解定解问题初值问题 初边值问题还需对定解条件进行离散化。 对初始条件及第一类边界条件，可直接得到 其中对第二、三类边界条件 需用差分近似。下面介绍两种较简单的处理方法。（1） 在左边界处用向前差商近似偏导数，在右边界处用向后差商近似，即 则得边界条件的差分近似为 其截断误差为

11、。 （2）用中心差商近似,即 则得边界条件的差分近似为 其截断误差为。误差的阶数提高了，但出现定解区域外的节点和，这就需要将解拓展到定解区域外。可以通过用内节点上的值插值求出和，也可以假定热传导方程在边界上也成立，将差分方程扩展到边界节点上，由此消去和。 （三）几种常用的差分格式以热传导方程的初边值问题为例给出几种常用的差分格式。（1）古典显式格式令，则 可改写成 将其与初始条件及第一类边界条件 结合，我们得到求解此问题的一种差分格式 由于第0层上节点处的值已知，由此即可算出在第一层上节点处的近似值。重复使用此式，可以逐层计算出所有的，因此此差分格式称为古典显式格式。又因式中只出现相邻两个时间

12、层的节点，故此式是二层显式格式。 （2）古典隐式格式 将式 整理并与初始条件及第一类边界条件式联立，得差分格式如下 其中。虽然第0层上的值仍为已知，但不能由上式直接计算以上各层节点上的值，必须通过解下列线性方程组 才能由计算，故此差分格式称为古典隐式格式。此方程组是三对角方程组，且系数矩阵严格对角占优，故解存在唯一。 （3）Richardson格式Richardson格式是将式整理后与初始条件及第一类边界条件式联立。其计算公式为 这种差分格式中所涉及的节点出现在三层上，故为三层显式格式。Richardson格式是一种完全不稳定的差分格式，因此它在实际计算中是不能采用的。 （4）杜福特-弗兰克尔

13、（DoFort-Frankel）格式DoFort-Frankel格式也是三层显式格式，它是由式与初始条件及第一类边界条件式结合得到的。具体形式如下：用这种格式求解时，除了第0层上的值由初值条件得到，必须先用二层格式求出第1层上的值，然后再按上式逐层计算。 （5）六点隐式格式对二阶中心差商公式如果用在点与点处的二阶中心差商的平均值近似在处的值，即 同时在点处的值也用中心差商近似，即 这样又得到热传导方程的一种差分近似 其截断误差为，将上式与初始条件及第一类边界条件式联立并整理，得差分格式 此格式涉及到六个节点，它又是隐式格式，故称为六点隐式格式。与古典隐式格式类似，用六点格式由第层的值计算第层的

14、值时，需求解三对角方程组 此方程组的系数矩阵严格对角占优，故仍可用追赶法求解。 例2 用古典显式格式求初边值问题 的数值解，取。解 这里，。由格式可得到 将初值代入上式，即可算出将边界条件，及上述结果代入又可求得如此逐层计算，得全部节点上的数值解为 §4 双曲型方程的差分解法 对二阶波动方程如果令，则方程可化成一阶线性双曲型方程组 记，则方程组可表成矩阵形式 矩阵有两个不同的特征值，故存在非奇异矩阵，使得作变换，方程组可化为 方程组由二个独立的一阶双曲型方程联立而成。因此本节主要讨论一阶双曲型方程的差分解法。4.1几种简单的差分格式 考虑一阶双曲型方程的初值问题 将平面剖分成矩形网格，取方向步长为方向步长为，网格线为，。为简便，记，。 以不同的差商近似偏导数，可以得到方程的不同的差分近似 截断误差分别为，与。结合离散化的初始条件，可以得到几种简单的差分格式 其中。如果已知第层节点上的值，按上面三种格式就可