第四章 复习题 1、 试简要说明对导热问题进行有限差分数值计算的基本思想与步骤。 2、 试说明用热平衡法建立节点温度离散方程的基本思想。 3、 推导导热微分方程的步骤和过程与用热平衡法建立节点温度离散方程的过程十分相似，为什么前者得到的是精确描述，而后者解出的确实近似解。 4、 第三类边界条件边界节点的离散那方程，也可用将第三类边界条件表达式中的一阶导数用差分公式表示来建立。试比较这样建立起来的离散方程与用热平衡建立起来的离散方程的异同与优劣。 5．对绝热边界条件的数值处理本章采用了哪些方法？试分析比较之． 6．什么是非稳态导热问题的显示格式？什么是显示格式计算中的稳定性问题？ 7．用高斯－塞德尔迭代法求解代数方程时是否一定可以得到收敛德解？不能得出收敛的解时是否因为初场的假设不合适而造成？

8．有人对一阶导数221,253xtttxtinininin 你能否判断这一表达式是否正确，为什么？ 一般性数值计算 4-1、采用计算机进行数值计算不仅是求解偏微分方程的有力工具，而且对一些复杂的经验公式及用无穷级数表示的分析解，也常用计算机来获得数值结果。试用数值方法对

Bi=0.1,1,10的三种情况计算下列特征方程的根:)6,2,1(nn 3,2,1,tannBi

nn



并用计算机查明，当2.02aFo时用式（3-19）表示的级数的第一项代替整个级数（计算中用前六项之和来替代）可能引起的误差。

解：Binntan，不同Bi下前六个根如下表所示： Bi μ1 μ2 μ3 μ4 μ5 μ6

0.1 0.3111 3.1731 6.2991 9.4354 12.5743 15.7143

1.0 0.8603 3.4256 6.4373 9.5293 12.6453 15.7713 10 1.4289 4.3058 7.2281 10.2003 13.2142 16.2594

Fo=0.2及0.24时计算结果的对比列于下表： Fo=0.2 x Bi=0.1 Bi=1 Bi=10 第一项的值 0.94879 0.62945 0.11866 前六和的值 0.95142 0.64339 0.12248 比值 0.99724 0.97833 0.96881 Fo=0.2 0x Bi=0.1 Bi=1 Bi=10 第一项的值 0.99662 0.96514 0.83889 前六项和的值 0.994 0.95064 0.82925 比值 1.002 1.01525 1.01163 Fo=0.24 x Bi=0.1 Bi=1 Bi=10 第一项的值 0.94513 0.61108 0.10935 前六项的值 0.94688 0.6198 0.11117 比值 0.99814 0.98694 0.98364 Fo=0.24 0x Bi=0.1 Bi=1 Bi=10 第一项的值 0.99277 0.93698 0.77311 前六项和的值 0.99101 0.92791 0.76851 比值 1.00177 1.00978 1.00598

4-2、试用数值计算证实，对方程组







5223122321321321xxxxxxxxx

用高斯-赛德尔迭代法求解，其结果是发散的，并分析其原因。 解：将上式写成下列迭代形式 







2131323213212/1252/1xxxxxxxxx

假设3,2xx初值为0，迭代结果如下： 迭代次数 0 1 2 3 4

1x 0 2.5 2.625 2.09375 2.6328125

2x 0 -0.75 0.4375 - 1.171875 1.26171825

3x 0 1.25 -0.0625 2.078125 -0.89453125

显然，方程迭代过程发散 因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值代数和。

4-3、试对附图所示的常物性，无内热源的二维稳态导热问题用高斯-赛德尔迭代法计算4321,,,tttt之值。

解：温度关系式为： 











5104/115304/130204/130404/1324413412321tttttttttttt

开始时假设取200201tt℃；150403tt℃ 得迭代值汇总于表 迭代次数 0 20 20 15 15 1 26.25 22.8125 21.5625 14.84375 2 28.59375 23.359375 22.109375 15.1171875 3 28.8671875 23.49609375 22.24607565 15.18554258 4 28.93554258 23.53027129 22.28027129 15.20263565 5 28.95263565 23.53881782 22.28881782 15.20690891 6 28.9569089 23.54095446 22.290955445 15..20797723

其中第五次与第六次相对偏差已小于410迭代终止。

4-4、试对附图所示的等截面直肋的稳态导热问题用数值方法求解节点2，3的温度。图中)./(30,25,852000KmWhCtCtf.肋高H=4cm,纵

剖面面积,42cmAL导热系数)./(20KmW。

解：对于2点可以列出：

节点2：;0)(2214321ttxhxttxtt

节点3：0)(22)(23132ttxhtthxttff。 由此得：

0)(22122321ttxhtttt，0)(2)(32332tthxtthttff，



22222312xHhtxHhtttf



2122223xhhtxhth

ttff

06.001.02002.03022xh，于是有：12.0212.0212ftttt， 

53.253.153.203.05.103.020/30103.020/302f223ffffttttttttt

＝

，代入得：

ffttttt12.053.253.112.2212

，f2123036.053.153.23636.5tttttf，

fttt8336.153.23636.412，3636.48336.153.22ffttt，

Ct8.5979.593636.484.4505.2153636.4258336.18553.22，

Ct8.3875.3853.22553.18.593。

离散方程的建立 4-5、试将直角坐标中的常物性无内热源的二维稳态导热微分方程化为显式差分格式，并指

出其稳定性条件（)yx。 解：常物性无内热源二维非稳态方程微分方程为 



2222ytxtat

 扩散项取中心差分，非稳态项取向前差分： 





211211122ytttxtttatt

ininininininini

n

所以有 inininintyxattyxat2211221

112111



稳定性条件 2/1yxFoFo 4-6、极坐标中常物性无内热源的非稳态导热微分方程为





2222211trrtrrtat

试利用本题附图中的符号，列出节点（i,j）的差分方程式。

解：将控制方程中的各阶导数用相应的差分表示式代替，可得：

1111111122222112kkikkkkkkkkjtjtjtjjtjjijijijjjtttttttttrarrrr，，，，，，，，，。

也可采用热平衡法。对于图中打阴影线的控制容积写出热平衡式得： 

111kkkkkkijijijijijijjjjttttttrrcrrrr

，，，，，，

11,22kkkkijijijijjjttttrrrrrr，，， 对等式两边同除以jrr并简化，可以得出与上式完全一样相同的结果。 4-7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷却，底面可以认为是绝热的。为用数值法确定冷却过程中柱体温度的变化，取中心角为1rad的区域来研究（如本题附图所示）。已知柱体表面发射率，自然对流表面传热系数，环境温度，金属的热扩散率，试列出图中节点（1，1），（M,1）(M,n)及（M,N）的离散方程式。在r及z方向上网格是各自均分的。

解：应用热平衡法来建立四个节点点离散方程。 节点（1，1）： 2121211211111111222282kkkkkkttttttrrzrzczr













，，，，，，

节点（m，1）：