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

WRITE (6, 122) C1(1, J), I=1, N1), J=1, MX) ROFMAT (3X, ‘I=’, I2, 3X, J=’, 12, 3X, 122 $‘C1 (I,J)=’, F10.5/) stop end The END ! 思考题: 上述几何区域上两维扩散方程的有限差分法求解. 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * * * 若 ， ，则式（2.22）又被 改写成： 或 （2.23） 物理意义： 一点（i, j）处的温度是它周围4点温度的平均值。 由于差分方程（2.21）是从式（2.20）中去掉项 得来的，称去掉的项 为差分方程（2.21）的截断误差。当 与 趋于零时，差分方程的截断误差也趋于零，即差分方程逼近微分方程。我们称这种逼近的差分方程与相应的微分方程为“相容”。 相容性问题 3．边界条件的差分形式 对流换热边界条件： （2.14） ---用T 对 x 的向前差商代替式（2.14）中的 T 对 x 的一阶偏微商，使式（2.14）变 成为如下差分形式 ： 这里介绍用差商代替微商的办法把定解问题中 的各种边界条件表示成差分的形式. 或 （2.24） （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.25） 或 绝热边界条件： 变成为： （i=5；j=2，3，4） （2.26） （2.16） 给定温度边界条件： （i=1，2，3，4，5；j=5 （2.27） 至此，我们对全部节点，包括内节点与边界节点，都用差分形式代替了原来的函数形式。对于内节点上差分形式，我们通称差分方程，因为内节点上温度都是未知的。对于边界节点的差分形式，在边界节点的温度为未知量时，它是差分方程。而对于边界节点为给定的温度时，得到的就不是差分方程了。但在实际应用中，人们往往习惯地把由内节点与边界节点建立起来的差分形式，都统称为差分方程。笼统地讲，一个节点对应一个差分方程。 在边界节点的处理方面还有几点需要强调： （i） 每一边界节点只应属于一种边界条件。如图2.2 中，i=1，j=5的节点只属于边界条件式（2.16）。 （iii）若边界节点不正好落在区域的边界 上，则需对它们进行特殊处理。 （ii）对应不同边界的差分方程式（2.24）、2.25）、 （2.26）都是用一阶向前差商代替一阶微商得到的， 也即它们的截断误差为 或 量级，与内节点差分方程的截断误差相比， 低了一个量级。这一点也是从微分形式出发 建立差分格式的弱点。 4．差分格式的构成 由于式（2.13）~（2.17）所表示的方程式与边界条件都是线性的，由此而得到的内节点与边界节点的差分方程也都是线性代数方程。由全部节点的差分方程构成一个线性代数方程组。在这个方程组里，方程式的个数等于节点的个数。针对前面讨论的例子，我们可以看到，这个有25个节点所构成的方程组，只需要5个式子，即式（2.21）、（2.24）~（2.27）就可表示，这里每个式子都表示几个节点方程。 也就是说，在组成方程组时，不必把每个节点方程都写出来，而只要写出几个规格化了的方程就可以了。因此，人们把规格化了的，由内节点与边界节点全部差分方程所构成的线性代数方程组，称之为“差分格式”。 一般地说，差分格式被写成如下的形式： （2.28） 其中，n是节点数，也即方程个数，每个方程对应一个节点。 和bi（i=1，2，…，n；j=1，2，…，n）都是常数 ， 方程组（2.28）可被进一步写成矩阵的形式 : （2.29） 其中 如果在方程组中去掉其中已知温度节点的那些方程，由此构成的线性代数方程组中，方程的个数等于温度未知的节点个数，也即方程组的未知数。这样的方程组也是差分格式。也可写成式（2.29）。 综上所述，用有限差分法对