第5章导热问题数值解.pptVIP

梁秀俊 高等传热学 2．边界节点离散方程的建立 采用热平衡的方法建立。 梁秀俊 高等传热学 热平衡法的优点 （1）物理意义明确； （2）适用于内节点和边界节点。 对于其它外节点，如外部角点和内部角点等可按同样的方法推导，区别是其所代表的控制体体积不同。 梁秀俊 高等传热学 写出所有内节点和边界节点的温度差分方程n个未知节点温度，n个代数方程式： 代数方程组的求解方法：直接解法、迭代解法 五、节点方程组的求解 梁秀俊 高等传热学 直接解法：通过有限次运算获得代数方程精确解; 如：矩阵求逆、高斯消元法等。 迭代解法：先对要计算的场作出假设、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。 缺点：所需内存较大、方程数目多时不便、不适用于非线性问题。 迭代解法有多种：高斯-赛德尔迭代是最常用的。高斯-赛德尔迭代的特点：每次迭代时总是使用节点的最新值。 梁秀俊 高等传热学 在计算后面的节点温度时应按下式（采用最新值） 例如：根据第 k 次迭代的数值 可以求得节点温度： 梁秀俊 高等传热学 判断迭代是否收敛的准则： 6 3 10 ~ 10 - - 取 值一般 相对偏差 允许的偏差； — e e k及k+1表示迭代次数； —第k次迭代得到的最大值 梁秀俊 高等传热学 对于一个代数方程组，若选用的迭代方式不合适，有可能导致发散，即称迭代过程发散； 对于常物性导热问题，组成的差分方程组，迭代公式的选择应使一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值的代数和，此时，结果一定收敛。 这一条件数学上称主对角线占优（对角占优）； 梁秀俊 高等传热学 例题5-1 有一根截面尺寸为1m×1m的很长的方柱，导热系数为1 W / (m·K)，在稳态条件下，柱体的三个面表维持在500℃，另一个表面暴露在温度是300℃的气体中，气体与柱体间的表面传热系数为10 W / (m2·K)。利用 的网格，确定柱体中的二维温度分布。 梁秀俊 高等传热学 梁秀俊 高等传热学 梁秀俊 高等传热学 计算次数 t1 t2 t3 t4 t5 t6 t7 t8 0 480.0 470.0 440.0 430.0 400.0 390.0 370.0 350.0 1 477.5 471.3 451.9 441.3 428.0 411.8 355.2 337.3 2 480.8 475.7 462.5 453.1 432.6 413.9 355.8 337.7 3 484.6 480.6 467.6 457.4 434.3 415.9 355.2 338.3 4 487.0 482.9 469.7 459.6 435.5 417.2 355.6 338.6 5 488.1 484.0 470.8 460.7 435.1 417.9 355.7 338.8 6 488.7 484.5 471.4 461.3 435.5 418.3 355.9 338.9 7 489.0 484.8 471.7 461.6 435.7 418.5 355.9 339.0 8 489.1 485.0 471.9 461.8 435.8 418.6 355.9 339.0 梁秀俊 高等传热学 练习：试用数值计算求解方程组 用高斯-赛德尔迭代法求解，其结果是发散的还是收敛的？为什么？ 梁秀俊 高等传热学 　 x1 x2 x3 0 0 0 0 1 31 -42 -44 2 373 -550 -650 3 4731 -6842 -8644 4 60173 -85750 -112050 5 764931 -1081642 -1437244 6 9719973 -1.4E+07 -1.8E+07 7 1.2E+08 -1.7E+08 -2.3E+08 8 1.6E+09 -2.2E+09 -3E+09 9 2E+10 -2.8E+10 -3.8E+10 10 2.5E+11 -3.6E+11 -4.8E+11 11 3.2E+12 -4.5E+12 -6.1E+12 12 4.1E+13 -5.7E+13 -7.7E+13 　 x1 x2 x3 0 0 0 0 1 4.25 5.35 1.5375 2 1.8594 4.9056 2.8439 3 1.6994 4.1538 3.1118 4 1.9144 3.95 3.0553 5 2.0049 3.9658 3.0061 6 2.0113 3.9941 2.9958 7 2.0033 4.0018 2.9979 8 1.9998 4.0013 2.9998 9 1.9996 4.0002 3.0002 10