《数值传热学作业答案第7章题解》.doc

第7章题解 7-1 试用计算证实，对下列方程组 采用Jacobi点迭代时是收敛的，而采用G-S点迭代时则是发散的。 解：先采用Jacobi点迭代方法 将方程组写成Jacobi点迭代形式 其矩阵形式为 用雅可比迭代法计算方程组的解。 表7-1-1 Jacobi迭代法的各次结果 迭代次数 0 0 0 0 1 1 3 5 2 5 -3 -3 3 1 1 1 4 1 1 1 再采用G-S点迭代方法 将方程组写成G-S点迭代形式 其矩阵形式为 用G-S迭代法计算方程组的解。 表7-1-2 G-S迭代法的各次结果 迭代次数 0 0 0 0 1 1 2 -1 2 -5 9 -3 3 -23 29 -7 4 -71 81 -15 5 -191 209 -31 6 -479 513 -63 7 -1151 1217 -127 8 -2687 2817 -255 从表7-1-2中明显可以看出采用G-S方法计算该方程组时发散。 7-4 一正方形导热物体的底边是绝热的，其余三边上的温度如图7-14所示。试确定该正方形内节点1、2、3、4的温度，导热物体无内热源，物性为常数。 图7-14 习题7-4插图 解：导热物体无内热源，物性为常数的控制方程可以写成 7-4-1 考虑到空间步长相同时，离散方程可以写成： 7-4-2 式中 从而有： 7-4-3 对节点1列方程，并整理： 即 7-4-4 同理，对节点2列方程，并整理： 7-4-5 考虑到下底面是绝热面，所以有，这里可以认为下底面上的节点温度对应节点，的温度。 所以对节点3，4列方程，有： 即 7-4-6 7-4-7 对上面各节点的离散方程进行联立求解。 下面列出计算结果，见表7-4。 表7-4 各节点的计算结果 28.737 24.263 20.684 18.316 7-6 对如图7-15所示的正方形内稳态无内热源的常物性导热问题，试分别用G-S点迭代及线迭代求解节点1、2、3、4的温度，比较其收敛速度。将求解结果与例1进行对比，并解释所观察到的事实。 图7-15 习题7-6插图 解：导热物体无内热源，物性为常数的控制方程可以写成 7-6-1 考虑到空间步长相同时，离散方程可以写成： 7-6-2 式中 从而有： 7-6-3 列出节点1，2，3，4的温度方程，有： 7-6-4 7-6-5 7-6-6 7-6-7 对上面提出的方程组，采用采用G-S点迭代方法以及线性迭代法进行求解。 表7-6-1 G-S迭代法的各次结果 迭代次数 0 35 35 35 35 1 30 38.75 33.75 45.625 2 30.625 41.563 36.563 47.031 3 32.031 42.266 37.266 47.383 4 32.383 42.441 37.441 47.471 5 32.471 42.485 37.485 47.493 6 32.493 42.496 37.496 47.498 7 32.498 4