研一 流体力学课件 传热与流体流动的数值计算.pptVIP

* * §10 小结 本章完成了我们传热与流体流动数值计算方法的最后一步。对于完整的非稳态对流扩散问题已有了全貌的认识。下一步将探讨一下求解方法和一些专题及应用。 ? 作业：3、 4、5、6 * * 第七章 解法与最后的完善 引言：回顾与回忆 重点内容： ① 线性代数方程的解． ② 代数方程的解． ③ 超松弛与欠松弛． ④ 最后的修饰． * * 仍对一维问题的所有未知温度点： c1=bN=0，边界条件 ? 这种方程可由标准的Gauss消元法（列消元、全主元消元法等）。这种算法有时称为Thomas算法或三对角矩阵法、TDMA方法。矩阵形式如下： i∈[1, N] (7.1) §1.线性代数方程的解 P W E 1 i-1 i i+1 N * * A T=d c1=0，T0不起作用，归入d1 bN=0，TN+1不起作用，归入dN。? * * 采用边界条件下，表明T1是T2的已知函数。 如i=2时，方程表示T1、T2与T3间的关系式，但T1又可表示为T2的函数关系，所以T1、T2、T3之间可最终表示成T2、T3的关系，依次类推直至TN表示成TN+1的函数关系，但TN+1不存在，由此得到TN的值，之后开始“回代”过程，即 由TN TN-1 TN-2......T2 最后得到 T1。 * * 计算过程如下： 前 向 代 入 T1=f(T2) T2=f(T3) ... ... Ti=f(Ti+1) ... ... TN=c 回 代 过 程 具体做法： 1)前向代入时，得到： (7.3) * * 2)将Ti-1代入(7.1)得到： (7.4) 这样，即可得到前向过程递推的Pi、Qi与 Pi-1、Qi-1的关系? 3)再来看在递推过程开始时 对i=1, 方程(7.1) (7.1) (7.1) (7.1) (7.6) * * c1=0 (7.7) [或直接在方程(7.6)中代入c1=0 可得到方程(7.7)的表达式] 在Pi、Qi序列的另一端，bN=0 PN=0 (7.8) * * 4)此后，开始回代过程 ? 回代过程概要： (1)由（7.7）得到P1、Q1 (2)由（7.6）递推得到i=2, 3, 4, ......, N 得到 Pi, Qi, i=2, 3, 4, ......, N (3)令TN=QN (4)对i=N-1, N-2, N-3, ......, 3, 2, 1，应用（7.4）得到TN-1, TN-2, TN-3, ......, T3, T2, T1 说明：TDMA方法是一个非常有效的求解方法。所需计算机存储量及计算机时仅正比于N而不是N2或N3。 * * 对离散化方程可采用任何合适的方法求解 TDMA 方法不能直接求解2－D、3－D问题。由于TDMA直接法是非常复杂且需要相当大量的计算机容量和时间。 对非线性问题，系数矩阵随时在变，采用直接解法不经济。 迭代方法：如前所述，由某一估计值开始求解代数方程，获得改进的新值。重复这一过程最后得到一个充分接近代数方程精确解的解。对非线性问题，迭代方法特别有吸引力。 §2. 代数方程的解 * * 但要注意：计算系数所付出的代价与求解方程所需的时间之间应当存在某一平衡。计算出一组系数后，应充分迭代从当前系数上得到足够的好处，但不要过于多地利用这些中间值的系数。 1.Gauss-Seidel迭代，逐点计算法：最简单的迭代法。 过程：按一定顺序逐个访问每一个结点，计算该处变量值。计算机内只存储一组T值。开始，这一组T值代表最初的估计值或上一次迭代得到的值。 (7.21) * * 实例： T1=0.4T2+0.2 T2=T1+1 迭代过程： 迭代次数 0 1 2 3 4 5 6 ∞ T1 0 0.2 0.68 0.872 0.949 0.980 0.992 1.0 T2 0 1.2 1.68 1.872 1.949 1.980 1.992 2.0 * * 式(6.32)即为速度修正公式，它也可写作： 这就表明了对压力的修正，ue*如何进行修正以产生ue，或ue与p的相应关系。 (6.34) 类似地可得到： 至此，就逐渐得到了p的离散化方程，下面来具体解决这个问题。 连续性方程 p 方程 u，v，w满足连续性方程 (6-35) (6-36) * * §6