第二章 有限差分法初步

一、差商与微商第二章：有限差分法初步§1有限差分法基本概念（i）、有限差分的数学基础是用差商代替微商。有如下两种数学形式：（i）微商（导数）的定义)(xT是连续函数，则它的导数为：若xTxxTxxTdxdTxx00lim)()(lim（2.1）式（2.1）右边xT是有限的差商。x与T都不为零，而式（2.1）左边dxdT是xT当x趋于零时极限情形下的差商，称之微商。在x没有到达零之前，xT只是dxdT的近似。xT趋于dxdT的过程认为是近似向精确过渡，用xT代替dxdT就是精确向近似过渡。两者的差值dxdTxT表示差商代替微商的偏差。（ii）偏差---Taylor级数展开222!2)()()(dxTdxdxdTxxTxxTnnndxTdnx!)(（2.2）稍加整理后可写成：22!2)()(dxTdxdxdTxxTxxTxTnnndxTdnx!)(1（2.3）xT可见与只能是近似相等。dxdT偏差为：)(0x（iii）微商与差商的几何意义xx+xx-xT(x+x)T(x-x)T(x)T(x+x)-T(x-x)2xT(x+x)-T(x)xT(x)-T(x-x)xdT(x)dx图2-1差商与微商的比较图2.1表示了差商与微商之间的关系。应当指出，用不同方法得到的差商去代替微商，它们带来的偏差是不同的。向右（前）差商：（iV）差商的几种表示xxTxxTdxdT)()(（2.4）向左（后）差商：xxxTxTdxdT)()(（2.5）中心差商，取向右差商与向左差商的平均值：xxxTxTxxTxxTdxdT)()()()(21xxxTxxT2)()(（2.6）偏差分析：将Taylor级数写成：)(2)()()()(2xTxxTxxTxxT（2.7）Taylor级数还可写成：))(()(!3)(43xOxTx（2.8）))(()(!3)(43xOxTx)(!2)()()()(2xTxxTxxTxxT由式（2.7）可得)(!2)()()(xTxxTxxTxxT)(xO由式（2.8）可得)(!2)()()(xTxxTxxxTxT（2.9）)(xO（2.10）（2.9）+（2.10），得到)(!3)()(2)()(2xTxxTxxxTxxT2)(xO（2.11）比较式（2.9）、（2.10）、（2.11）可看到，用不同的差商形式去代替微商，所带来的偏差是不同的。这些偏差都是截去了Taylor级数展开式中的高阶项而引起的，常称“截断误差”。用向右差商与向左差商代替微商，其截断误差为与x同量级的小量；)(xO2)(x同量级的小量；中心差商的截断误差小于向右差商或向左差商。而用中心差商代替微商，其截断误差是与讨论：上述一阶差商一般仍是x的函数，对它们还可以求差商。这种一阶差商的差商称为二阶差商，它是二阶微商的近似，常用向右差商的向左差商来近似二阶微商，即：xxxxTxTxxTxxTdxTd)()()()(222)()()(2)(xxxTxTxxT（V）二阶差商根据式（2.7）+（2.8），可得)()()()(2)(2xTxxxTxTxxT2)(xO（2.12）由式（2.12）知，二阶差商的截断误差也为与同阶的小量。2)(x结论：由于用差商代替微商必然带来截断误差，相应地用差分方程代替微分方程也必然带来截断误差。这是有限差分法固有的。因此，在应用有限差分法进行数值解时，必须对差分的构成及其对方程造成的误差引起注意。二、从微分形式出发的差分格式图2.2给出了一个简单边界值问题。(i,j)(i+1,j)(i-1,j)(i,j-1)(i,j+1).....q’’TwThxyL1L2图2.2矩形区域离散化问题是求图2.2所示的边值问题的解，其数学表达如下，方程：02222kqyTxT（2.13）边界条件：0x20Ly)(TThxTk0y10LxqyTk1Lx20Ly0xT2Ly10LxwTT（2.14）（2.15）（2.15）（2.16）式（2.13）、（2.14）、（2.15）、（2.16）所示定解问题解法。在问题的提法已经明白之后，差分格式的构成可通过以下几步来实现：（i）区域离散法下面分别予以说明:（vi）构成差分格式（ii）建立区域内差分方程（iii）边界条件的差分形式1．区域离散化所谓离散化，就是把几何上连续的区域用一系列网格线把它划分开。一般说来，网格形式应视几何区域的不同而不同，对于矩形区域而言，用矩形的网格，如图2.2，用五条水平网线与五条垂直网线把矩形区域离散掉。网线与网线的交点称之为“节点”，节点与节点的距离称之为步长，x方向的步长表示为，y方向的步长表示为。xy节点编号：为便于计算，需对节点逐个编号。常用（i，j）表示节点位置，其中，i、j是与网线相对应的正整数。i，j的排列：可有不同的方式。习惯上，与x、y轴相一致，i由左而右逐个增长，j由下而上逐个增长。但也有，考虑到与矩阵的格式相一致，i表示行数，由上而下逐个增长，j表示列数，由左而右逐个增长。这种从上到下，从左到右的编排与一般书写习惯也是一致的，因此，在计算机上算题也常被采用。在本章中，大都采用与坐标相一致的编排方法。在区域内的节点称“内节点”，在边界上的节点称“边界节点”。图2.2所示边界是规则的,则节点或在区域内，或正好落在边界上。xy步长或可以是不变的常量，即等步长，也可以在区域内的不同处是不同的，即变步长。如果区域内各处的温度梯度变化很大，则在温度变化剧烈处，网格布得密些；在温度变化不剧烈处，网格布得疏些。至于网格布置多少，步长取多大为宜，要根据具体问题，兼顾到计算的精确度与计算的工作量等因素而定。步长:从物理方面对区域离散化可作这样的理解，即认为区域内离散的每个节点，都集中着它周围区域（尺度为步长）的热容，或者说，区域内连续分布的热容都被分别地集中到离散的节点上去了。这样，节点的温度代表着它周围区域的某种平均温度。一系列离散的节点温度值代表着连续区域内的温度分布。jiT,区域离散化物理理解:节点（i，j）处的温度表示成。2．差分方程代替微分方程在上节我们已对有限差分法的数学基础作了简要的介绍，说明了如何用差商代替微商，以及由此带来的误差。这里介绍用差商代替微商的办法来处理导热方程（2.13），得到相应的差分方程。方程（2.13）02222kqyTxT对区域内各个点都成立的，当然对任意一个内节点（i，j）也成立。或者说，在（i，j）处存在二阶偏微商与，这些二阶偏微商所对应的差商可表示成：jixT,22jiyT,222,1,122)(2xTT,TxTjijiji,ji)(2xOi=2，3，4；j=2，3，4（2.17）21,,122)(2yTT,TyTjijiji,ji)(2yO（2.18）i=2，3，4；j=2，3，4其中，与表示相应的二阶差商与二阶偏微商的差别为与的数量级。)(2xO)(2yO2)(x2)(y将式（2.17）与（2.18）代入方程（2.13），得21,,1,2,1,,1)(2)(2yTTTxTTTjijijijijiji0))()((22yxOkq（2.19）式（2.19）中去掉项，得到))()((22yxO21,,1,2,1,,1)(2)(2yTTTxTTTjijijijijiji0kq（2.20）i=2，3，4；j=2，3，4式（2.20）被称为对应于方程（2.13）的差分方程。方程（2.20）被改写成：1,2,12,22)(1)(1)(2)(2jijijiTxTxTyxkqTyTyjiji1,21,2)(1)(1（2.21）若，，则式（2.21）又被改写成：yx0q041,1,,1,1,jijijijijiTTTTT或)(411,1,,1,1,jijijijijiTTTTT（2.22）物理意义：一点（i,j）处的温度是它周围4点温度的平均值。))()((22yxO))()((22yxO由于差分方程（2.20）是从式（2.19）中去掉项得来的，称去掉的项为差分方程（2.19）的截断误差。当与趋于零时，差分方程的截断误差也趋于零，即差分方程逼近微分方程。我们称这种逼近的差分方程与相应的微分方程为“相容”。xy3．边界条件的差分形式对流换热边界条件：20,0)(LyxTThxTk（2.14）---用T对x的向前差商代替式（2.14）中的T对x的一阶偏微商，使式（2.14）变成为如下差分形式：这里介绍用差商代替微商的办法把定解问题中的各种边界条件表示成差分的形式.)(,,,1TThxTTkjijijiTkxhTTkxhjiji,1,1或（2.23）i=1；j=2，3，4热流边界条件：10,0LxyqyTk（2.15）用T对y的向前差商代替式（2.15）中T对y的一阶偏微商，使式（2.15）变成为如下差分形式:qyTTkjiji,1,kyqTTjiji1,,i=1，2，3，4，5；j=1（2.24）或绝热边界条件：2Ly0Lx0xT，变成为：0TTj1iji,,i=5；j=2，3，4（2.25）（2.16）给定温度边界条件：wjiTT,i=1，2，3，4，5；j=5（2.26）至此，我们对全部节点，包括内节点与边界节点，都用差分形式代替了原来的函数形式。对于内节点上差分形式，我们通称差分方程，因为内节点上温度都是未知的。对于边界节点的差分形式，在边界节点的温度为未知量时，它是差分方程。而对于边界节点为给定的温度时，得到的就不是差分方程了。但在实际应用中，人们往往习惯地把由内节点与边界节点建立起来的差分形式，都统称为差分方程。笼统地讲，一个节点对应一个差分方程。在边界节点的处理方面还有几点需要强调：（i）每一边界节点只应属于一种边界条件。如图2.2中，i=1，j=5的节点只属于边界条件式（2.16）。（iii）若边界节点不正好落在区域的边界上，则需对它们进行特殊处理。（ii）对应不同边界的差分方程式（2.23）、2.24）、（2.25）都是用一阶向前差商代替一阶微商得到的，也即它们的截断误差为)(xO或)(yO量级，与内节点差分方程的截断误差相比，低了一个量级。这一点也是从微分形式出发建立差分格式的弱点。4．差分格式的构成由于式（2.13）、（2.14）、（2.15）、（2.16）、（2.17）所表示的方程式与边界条件都是线性的，由此而得到的内节点与边界节点的差分方程也都是线性代数方程。由全部节点的差分方程构成一个线性代数方程组。在这个方程组里，方程式的个数等于节点的个数。针对前面讨论的例子，我们可以看到，这个有25个节点所构成的方程组，只需要5个式子，即式（2.21）、（2.23）、（2.24）、（2.25）、（2.26）就可表示，这里每个式子都表示几个节点方程。也就是说，在组成方程组时，不必把每个节点方程都写出来，而只要写出几个规格化了的方程就可以了。因此，人们把规格化了的，由内节点与边界节点全部差分方程所构成的线性代数方程组，称之为“差分格式”。一般地说，差分格式被写成如下的