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

若令 ，则系数 和 的计 算公式应为： 类似地，分别从上述等价方程组的第三个方程减去第二个方程乘以 ，第四个方程减去第二个方程乘以 ，如此等等，即可进一步消去后面n-2个方程中的未知量T2，而将方程(2.39)变为如下等价形式： 3．边界条件的差分形式 对流换热边界条件： （2.14） ---用T 对 x 的向前差商代替式（2.14）中的 T 对 x 的一阶偏微商，使式（2.14）变 成为如下差分形式 ： 这里介绍用差商代替微商的办法把定解问题中 的各种边界条件表示成差分的形式. 或 （2.23） i=1；j=2，3，4 热流边界条件： （2.15） 用T 对 y 的向前差商代替式（2.15）中T 对 y的一阶偏微商，使式（2.15）变成为如下差分形式: i=1，2，3，4，5；j=1 （2.24） 或 绝热边界条件： 变成为： i=5；j=2，3，4 （2.25） （2.16） 给定温度边界条件： 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）都是用一阶向前差商代替一阶微商得到的， 也即它们的截断误差为 或 量级，与内节点差分方程的截断误差相比， 低了一个量级。这一点也是从微分形式出发 建立差分格式的弱点。 4．差分格式的构成 由于式（2.13）、（2.14）、（2.15）、（2.16）、（2.17）所表示的方程式与边界条件都是线性的，由此而得到的内节点与边界节点的差分方程也都是线性代数方程。由全部节点的差分方程构成一个线性代数方程组。在这个方程组里，方程式的个数等于节点的个数。针对前面讨论的例子，我们可以看到，这个有25个节点所构成的方程组，只需要5个式子，即式（2.21）、（2.23）、（2.24）、（2.25）、（2.26）就可表示，这里每个式子都表示几个节点方程。 也就是说，在组成方程组时，不必把每个节点方程都写出来，而只要写出几个规格化了的方程就可以了。因此，人们把规格化了的，由内节点与边界节点全部差分方程所构成的线性代数方程组，称之为“差分格式”。 一般地说，差分格式被写成如下的形式： （2.27） 其中，n是节点数，也即方程个数，每个方程对应一个节点。 和bi（i=1，2，…，n；j=1，2，…，n）都是常数 ， 方程组（2.27）可被进一步写成矩阵的形式 : （2.28） 其中 如果在方程组中去掉其中已知温度节点的那些方程，由此构成的线性代数方程组中，方程的个数等于温度未知的节点个数，也即方程组的未知数。这样的方程组也是差分格式。也可写成式（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） 式（2.29）或（2.30）构成差分格式。若在方程组中去掉已知温度节点所对应的方程，即第1至第3个方程，则式（4.2.24）被改写成： （2.31） 将式（2.31）与式（4.2.28）进行对照，即可得到矩阵 、 、 的各个元素。