传热学第四版课后题答案第四章

第四章复习题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nBinn并用计算机查明，当2.02aFo时用式（3-19）表示的级数的第一项代替整个级数（计算中用前六项之和来替代）可能引起的误差。解：Binntan，不同Bi下前六个根如下表所示：Biμ1μ2μ3μ4μ5μ60.10.31113.17316.29919.435412.574315.71431.00.86033.42566.43739.529312.645315.7713101.42894.30587.228110.200313.214216.2594Fo=0.2及0.24时计算结果的对比列于下表：Fo=0.2xBi=0.1Bi=1Bi=10第一项的值0.948790.629450.11866前六和的值0.951420.643390.12248比值0.997240.978330.96881Fo=0.20xBi=0.1Bi=1Bi=10第一项的值0.996620.965140.83889前六项和的值0.9940.950640.82925比值1.0021.015251.01163Fo=0.24xBi=0.1Bi=1Bi=10第一项的值0.945130.611080.10935前六项的值0.946880.61980.11117比值0.998140.986940.98364Fo=0.240xBi=0.1Bi=1Bi=10第一项的值0.992770.936980.77311前六项和的值0.991010.927910.76851比值1.001771.009781.005984-2、试用数值计算证实，对方程组5223122321321321xxxxxxxxx用高斯-赛德尔迭代法求解，其结果是发散的，并分析其原因。解：将上式写成下列迭代形式2131323213212/1252/1xxxxxxxxx假设3,2xx初值为0，迭代结果如下：迭代次数012341x02.52.6252.093752.63281252x0-0.750.4375-1.1718751.261718253x01.25-0.06252.078125-0.89453125显然，方程迭代过程发散因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值代数和。4-3、试对附图所示的常物性，无内热源的二维稳态导热问题用高斯-赛德尔迭代法计算4321,,,tttt之值。解：温度关系式为：5104/115304/130204/130404/1324413412321tttttttttttt开始时假设取200201tt℃；150403tt℃得迭代值汇总于表迭代次数020201515126.2522.812521.562514.84375228.5937523.35937522.10937515.1171875328.867187523.4960937522.2460756515.18554258428.9355425823.5302712922.2802712915.20263565528.9526356523.5388178222.2888178215.20690891628.956908923.5409544622.29095544515..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2122223xhhtxhthttff06.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扩散项取中心差分，非稳态项取向前差分：211211122ytttxtttattinininininininin所以有inininintyxattyxat2211221112111稳定性条件2/1yxFoFo4-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）：11,11112111;22222kkkkkkkkmmmmmmmmmmmmttttttttzzzzzrrrrcrrrrz，，，，，，，节点（m，n）：1,4411111003332242242422kkkkkkmnmnmnmnmnmnmmmmmmmmmmnmnttttttrrrrrrrzrzrrzrrhttTTcrzr，，，，，，，。4-8、一个二维物体的竖直表面收液体自然对流冷却，为考虑局部表面传热系数的影响，表面传热系数采用25.11)(ttch来表示。试列出附图所示的稳态无内热源物体边界节点（M,n）的温度方程，并对如何求解这一方程提出你的看法。设网格均分。解：利用热平衡法：0.25MnfMnfhctttt，，，将h写为0.25MnfMnfhctttt，，，其中Mnt，为上一次迭代值，则方程即可线性化。4-9、在附图所示的有内热源的二维导热区域中，一个界面绝热，一个界面等温（包括节点4），其余两个界面与温度为ft的流体对流换热，h均匀，内热源强度为。试列出节点1，2，5，6，9，10的离散方程式。解：节点1：512111102242fttttxyxyyhttyx；节点2：32621210222ttttttyyxxyxxy；节点5：159565510222fttttttyxyxyyhttyyx；节点6：2676105560ttttttttxyxyxyyxyx；节点9：5910991022422fttttxyxyxyhttyx；节点10：91011106101010222fttttttyyxxyxhhtxxy。当xy以上诸式可简化为：节点1：252112202fhyhytttty；节点2：26132240tttty；节点5：2619522220fhyhyttttty节点6：271057640ttttty；节点9：25109122102fhyhytttty