第二章-有限差分法初步.pptVIP

结论： 由于用差商代替微商必然带来截断误差，相应地用差分方程代替微分方程也必然带来截断误差。这是有限差分法固有的。因此，在应用有限差分法进行数值解时，必须对差分的构成及其对方程造成的误差引起注意。 二、从微分形式出发的差分格式 图2.2给出了一个简单边界值问题。 问题是求图2.2所示的边值问题的解，其数学表达如下，方程： 边界条件： 在问题的提法已经明白之后，差分格式的构成可通过以下几步来实现： （i） 区域离散法 1．区域离散化 所谓离散化，就是把几何上连续的区域用一系列网格线把它划分开。一般说来，网格形式应视几何区域的不同而不同，对于矩形区域而言，用矩形的网格，如图2.2，用五条水平网线与五条垂直网线把矩形区域离散掉。网线与网线的交点称之为“节点”，节点与节点的距离称之为步长，x方向的步长表示为 ，y方向的步长表示为 。 节点编号：为便于计算，需对节点逐个编号。常用（i，j） 表示节点位置，其中，i、j是与网线相对应的正整数。 在区域内的节点称“内节点”，在边界上的节点称“边界节点”。图2.2所示边界是规则的,则节点或在区域内，或正好落在边界上。 从物理方面对区域离散化可作这样的理解，即认为区域内离散的每个节点，都集中着它周围区域（尺度为步长）的热容，或者说，区域内连续分布的热容都被分别地集中到离散的节点上去了。这样，节点的温度代表着它周围区域的某种平均温度。一系列离散的节点温度值代表着连续区域内的温度分布。 2．差分方程代替微分方程 在上节我们已对有限差分法的数学基础作了简要 的介绍，说明了如何用差商代替微商，以及由此带来的误差。这里介绍用差商代替微商的办法来处理导热方程（2.13），得到相应的差分方程。 或者说，在（i，j）处存在二阶偏微商 与 ，这些二阶偏微商所对应的差商可表示成： 将式（2.17）与（2.18）代入方程（2.13），得 式（2.19）中去掉 项，得到 3．边界条件的差分形式 对流换热边界条件： 用T 对 y 的向前差商代替式（2.15）中T 对 y的一阶偏微商，使式（2.15）变成为如下差分形式: 至此，我们对全部节点，包括内节点与边界节点，都用差分形式代替了原来的函数形式。对于内节点上差分形式，我们通称差分方程，因为内节点上温度都是未知的。对于边界节点的差分形式，在边界节点的温度为未知量时，它是差分方程。而对于边界节点为给定的温度时，得到的就不是差分方程了。但在实际应用中，人们往往习惯地把由内节点与边界节点建立起来的差分形式，都统称为差分方程。笼统地讲，一个节点对应一个差分方程。 4．差分格式的构成 也就是说，在组成方程组时，不必把每个节点方程都写出来，而只要写出几个规格化了的方程就可以了。因此，人们把规格化了的，由内节点与边界节点全部差分方程所构成的线性代数方程组，称之为“差分格式”。 一般地说，差分格式被写成如下的形式： 方程组（2.27）可被进一步写成矩阵的形式 : 如果在方程组中去掉其中已知温度节点的那些方程，由此构成的线性代数方程组中，方程的个数等于温度未知的节点个数，也即方程组的未知数。这样的方程组也是差分格式。也可写成式（2.28）。 综上所述，用有限差分法对式（2.13）、（2.14）、（2.15）、（2.16）、（2.17）所组成的边值问题的数值处理，最终归结成求解线性代数方程组（2.28），方程组的解即各节点的温度。如果整个区域的节点足够多，那么，离散节点的温度分布就近似代替了区域内的连续温度分布。 为便于讨论各种差分格式的优缺点，最好把方程组（2.28）中系数矩阵具体地写出来。但当我们着手书写由式（2.21）、（2.23）、（2.24）、（2.25）、（2.26）组成的代数方程组时，发现它所占的版面太大，造成印刷的困难。 所以，为便于书写，采用图2.3所示的网格 得到由全部节点组成的线性代数方程组为： 表示成矩阵形式： 式（2.29）或（2.30）构成差分格式。若在方程组中去掉已知温度节点所对应的方程，即第1至第3个方程，则式（4.2.24）被改写成： 将式（2.31）与式（4.2.28）进行对照，即可得到矩阵 、 、 的各个元素。 这里特别提醒读者注意，在式（2.31）与（2.28）中温度T下角码的对应关系，在讨论二维稳定导热问题时，人们习惯用二个角码（i，j）来表示节点位置及对应的曾度 。 （二）解线性代数方程组的直接法 这里介绍计算机上常用的解线性代数方程组的直接法。 大家知道，线性