第二十章 偏微分方程的数值解.doc

第二十章 偏微分方程的数值解 自然科学与工程技术中种种运动发展过程与平衡现象各自遵守一定的规律。这些规律的定量表述一般地呈现为关于含有未知函数及其导数的方程。我们将只含有未知多元函数及其偏导数的方程，称之为偏微分方程。 方程中出现的未知函数偏导数的最高阶数称为偏微分方程的阶。如果方程中对于未知函数和它的所有偏导数都是线性的，这样的方程称为线性偏微分方程，否则称它为非线性偏微分方程。 初始条件和边界条件称为定解条件，未附加定解条件的偏微分方程称为泛定方程。对于一个具体的问题，定解条件与泛定方程总是同时提出。定解条件与泛定方程作为一个整体，称为定解问题。 §1 偏微分方程的定解问题 各种物理性质的定常（即不随时间变化）过程，都可用椭圆型方程来描述。其最典型、最简单的形式是泊松(Poisson)方程 （1） 特别地，当时，即为拉普拉斯(Laplace)方程，又称为调和方程 （2） 带有稳定热源或内部无热源的稳定温度场的温度分布，不可压缩流体的稳定无旋流动及静电场的电势等均满足这类方程。 Poisson方程的第一边值问题为 （3） 其中为以为边界的有界区域，为分段光滑曲线，称为定解区域，分别为上的已知连续函数。 第二类和第三类边界条件可统一表示成 （4） 其中为边界的外法线方向。当时为第二类边界条件，时为第三类边界条件。 在研究热传导过程，气体扩散现象及电磁场的传播等随时间变化的非定常物理问题时，常常会遇到抛物型方程。其最简单的形式为一维热传导方程 （5） 方程（5）可以有两种不同类型的定解问题： 初值问题（也称为Cauchy问题） （6） 初边值问题 （7） 其中为已知函数，且满足连接条件 问题（7）中的边界条件称为第一类边界条件。第二类和第三类边界条件为 （8） 其中。当时，为第二类边界条件，否则称为第三类边界条件。 双曲型方程的最简单形式为一阶双曲型方程 （9） 物理中常见的一维振动与波动问题可用二阶波动方程 （10） 描述，它是双曲型方程的典型形式。方程（10）的初值问题为 （11） 边界条件一般也有三类，最简单的初边值问题为 如果偏微分方程定解问题的解存在，唯一且连续依赖于定解数据（即出现在方程和定解条件中的已知函数），则此定解问题是适定的。可以证明，上面所举各种定解问题都是适定的。 §2 偏微分方程的差分解法 差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题： （i）选取网格； （ii）对微分方程及定解条件选择差分近似，列出差分格式； （iii）求解差分格式； （iv）讨论差分格式解对于微分方程解的收敛性及误差估计。 下面我们只对偏微分方程的差分解法作一简要的介绍。 2.1 椭圆型方程第一边值问题的差分解法 以Poisson方程（1）为基本模型讨论第一边值问题的差分方法。 考虑Poisson方程的第一边值问题（3） 取分别为方向和方向的步长，以两族平行线 将定解区域剖分成矩形网格。节点的全体记为。定解区域内部的节点称为内点，记内点集为。边界与网格线的交点称为边界点，边界点全体记为。与节点沿方向或方向只差一个步长的点和称为节点的相邻节点。如果一个内点的四个相邻节点均属于，称为正则内点，正则内点的全体记为，至少有一个相邻节点不