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

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

2、法求解代数方程时是否一定可以得到收敛德解？不能得出收敛的解 时是否因为初场的假设不合适而造成？&有人对一阶导数x n,i3tni5tni12 x2tni你能否判断这一表达式是否正确，为什么?般性数值计算4-1、采用计算机进行数值计算不仅是求解偏微分方程的有力工具，而且对一些复杂的经验 公式及用无穷级数表示的分析解，也常用计算机来获得数值结果。试用数值方法对Bi=,1,10的三种情况计算下列特征方程的根n（n 12 ，6）：tanBin,n 1,2,3n并用计算机查明，当Fo 笃 0.22时用式（3-19）表示的级数的第一项代替整个级数（计算中用前六项之和来替代）可能引起的误差。解：ntan n

3、 Bi，不同Bi下前六个根如下表所示: Bi 1 2 3 4 5 6 /10Fo=及时计算结果的对比列于下表:Fo=XBi=Bi=1Bi=10第一项的值前六和的值比值Fo=x 0Bi= Bi=1Bi=10第一项的值前六项和的值比值Bi=Bi=1Bi=10第一项的值前六项的值比值Fo=X 0/ Bi=Bi=1Bi=10第一项的值前六项和的值/比值4-2、试用数值计算证实，对方程组Xi 2x2 2x31XiX2X32x1 2x2X35并分析其原因。用高斯-赛德尔迭代法求解，其结果是发散的, 解：将上式写成下列迭代形式X11/2 52X2X3X21/2 12X3X1X33X1X2假设X2,X3初值为

4、0,迭代结果如下:迭代次数0 1234X10X20-1.X30-0.显然，方程迭代过程发散因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值 代数和。-赛德尔迭代法计算4-3、试对附图所示的常物性，无内热源的二维稳态导热问题用高斯 讥3,上4之值。解：温度关系式为:t11/4t2t34030t21/4t1t42030t31/4t1t43015t41/4t2t3105t 0开始时假设取t1得迭代值汇总于表迭代次数0 201512020 C；20t4015 c15百度文库7345623.28.23.28.23.23.22.15.22.15.22.15.15.其中第五次与

5、第六次相对偏差已小于10 4迭代终止。4-4、试对附图所示的等截面直肋的稳态导热问题用数值方 法求解节点 2，3 的温度。图中to 850C,tf25C,h 3（W/（m2.K）.肋高日=4吓纵2剖面面积Al 4cm ,导热系数20W /（m，K）。t1 t2t3t4-2h x(t1 t2)0;节点2：xxt2 t3xh(tfX) 2h(tft3)0x2节点3：2。由此得：2h2 xh “，、t1t2t3t2 (t1t2)0t2t3(t ft3 )解：对于2点可以列出:x2 h(tf t3)t3t1t22h xH2t3t222hx230 0.02220 0.010.06于是有:t230/20t

6、f30/200.032.12tt30.03tft21.5tf4.3636t2t2t3xH2h x22.53t1t2 0.12tf2 0.12 ，t21.53tf-0.12tf 2.532.53t11.8336tf2.53 85 1.8336 254.363659.8 1.53 252.53t25.3636t22.53tf习題44附團0.03tf t21.53tf2.53，代入得:2.53t1 t21.53tf 0.3036tf1.8336t f4.3636215.05 45.844.363659.7959.8 C38.7538.8 C。离散方程的建立4-5、试将直角坐标中的常物性无内热源的二维

7、稳态导热微分方程化为显式差分格式，并指 出其稳定性条件（x y）。解：常物性无内热源二维非稳态方程微分方程为幕2t；ti2titin 1nn 12x丄i 1tn a稳定性条件Fo xti2 nyFo y1/21 2a4-6、极坐标中常物性无内热源的非稳态导热微分方程 为试利用本题附图中的符号, 程式。2tr列出节点（12t22ri,j ）的差分方解：将控制方程中的各阶导数用相应的差分表示式代 替，可得：t k i tk t1, jt t，jakkkt t，j 1 2t t，j t j 12r也可采用热平衡法。tkrjr c kk1 t t, j 1 t 1,rj2 r对于图中打阴影线的控制容积

8、写出热平衡式得: 1 . tki -i, ji, j1yj2r j1二頭1 ；（惘占tki,j1 tki, jrj对等式两边同除以rjtki1, j tkb jri,tki 1, jtki，jrrjtki，j1 tki,jrjrjr并简化，可以得出与上式完全一样相同的结果。t2t2ta 22x y扩散项取中心差分，非稳态项取向前差分:ti 1 tiinn4-7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷 却，底面可以认为是绝热的。为用数值法确定冷却过程中柱 体温度的变化，取中心角为1rad的区域来研究（如本题附图 所示）。已知柱体表面发射率，自然对流表面传热系数，环 境温度，金属的热扩散率

9、，试列出图中节点（1, 1）, （ M,1）（M,n）及（M,N ）的离散方程式。在r及z方向上网格是各自均分的。解：应用热平衡法来建立四个节点点离散方程。 节点（1, 1）：tktk11,2丄z 22 r2r8节点（m, 1）:所以有1,km,11,tkm,zztktkm,2m，1,rmr22zrm3rmnh Jtm,14oT4T4crmkm,1km,节点（m, n）：tktkt m 1, n t m, nm tm,m 1 3rmr-24-8、一个二维物体的竖直表面收液体自然对流冷却，为考虑局部表面传热系数的影响，表面传热系数采用1 25h c（t t1厂来表示。试列出附图所示的稳态无内热源

10、 物体边界节点（M,n）的温度方程，并对如何求解这一方 程提出你的看法。设网格均分。解：利用热平衡法：0m nrm 13rm40.25h C tM，n tf tM，n tf0.25将h写为h C tM，n tf tM，n tf ，其中tM，n为上 一次迭代值，则方程即可线性化。4-9、在附图所示的有内热源的二维导热区域中，一个界面绝热，一个界面等温（包括节点 4），其余两个界面与 温度为tf的流体对流换热，h均匀，内热源强度为。试列出节点1，2，5，6，9，10的离散方程式。七5 匕y；节点2：节点5：t2 t112 yh t1 tft1t2yt3 t2yx2x2t1t5yt9t5xy2y2t

11、2t6xt7t6yyxt5t9xt10 t9yy2x2t9t10yt11t10yx2x2节点6：节点9：节点10：t6 t21x x y y2t6t5t10t5xyt5 t6yxyh t5 tf 0；丄X y 丄h t9 tf 0422t6t10xyxh hi0 tf 0。当xy以上诸式可简化为:t5 t2节点1：百度文库22t6 ti t3 4t2y 0节点2:;h yh y22t6 ti t9 2 - tf 2 2 - t5 y-0节点5:2 t7 tio t5 t7 4t6 y 节点6:hy412t9y02;h y, 2 -22tio y0h yt5 tio tf 2 1节点9:h y2

12、t6 tg tii 2 tf节点io： 一维稳态导热计算h,肋片h同o 设2, 3,4-i0、一等截面直肋，高 H,厚，肋根温度为to，流体温度为tf，表面传热系数为 导热系数为。将它均分成4个节点（见附图），并对肋端为绝热及为对流边界条件（侧面）的两种情况列出节点 2，3，4的离散方程式H=45cm, i0mm,h 50W/（m2.K）, =5ow/, to i00 C, tf 20 C,计算节点 4的温度（对于肋端的两种边界条件）。4-10附冬解：采用热平衡法可列出节点2、3、4的离散方程为:节点2:节点3:节点t1 t2t3t22hx t2tft2 t3t4 t32hx t3tf4:肋端

13、绝热t4tft3 t4肋端对流Hh x t4tft4 t f93。将已知条件代入可得下列两方程组:肋端绝热t32.045t2i00.90t22.045t3t4 0.90t3i.0225t40.450肋端对流t32.045t2i00.90其中t2 2.045t3 t4 0.90百度文库15t3 1.0375t40.8由此解得：肋端绝热t292.2Ct291.5Ct3 87.70C , t4 t3 86.20C t4肋端对流-，肋端对流换热的条件使肋端温度更接近于流体温度。86.2C ；83.80C。4-11、复合材料在航空航天及化工等工业中日益得到广泛的应用。附图所示为双层圆筒壁, 设层间接触紧

14、密无接触12.5mm, r2 16mm,r318mm,2100CW/(m .K),tf2 60c, h240 W/ ,22380W/(m .K)。热阻存在。已知120W/（m.K）,tf1 150试用数值方法确定稳态时双层圆习駆All附馆h1筒壁截面上的温度分布。to解：采用计算机求解，答案从略。采用热平衡法对两层管子的各离散区域写出能量方程，进行求解；如果采用Taylor展开法列出方程，则需对两层管子单独进行， 并引入界面上温度连续及热流密度连续的条件，数值计算也需分两区进行，界面耦合。截面的温度分布定性地示于上图中。t。100 C,其表面上有自然对流散热，4-12、有一水平放置的等截面直杆

15、，根部温度t t / jni 1 / 41.75 oh Ct tf /d，其中，C 1.20W/（m . C）; d 为杆直径，m。杆高 H=10cm，直径d=1cm, = 50W/, t 25 C。不计辐射换热。试用数值方法确定长杆的散热量（需得出 与网格无关的解。杆的两端可认为是绝热的。解：数值求解过程略，Q=。的4-13在上题中考虑长杆与周围环境的辐射换热，其表面发射率为，环境可作为温度为 大空间，试重新计算其导热量。解：数值求解过程略，Q=。4-14、有如附图所示的一抛物线肋片，表面形线方程 为：肋根温度to及内热源流体温度tf为常数。定义: 试：（1 ）建立无量纲温度恒定,流体表面传

16、热系数h,t tft0 t fx/ H的控制方程；（2）在无量H2eb0.01,-0.05,纲参数t0 tfHH算。确定无量纲温度的分布。0.1,hH 0.01下对上述控制方程进行数量计d2 /d20.012 / 50。数值计算结果示于相对杆松.解：无量纲温度方程为：下图中，无量纲温度从肋根的1变化到肋端的。一维非稳态导热计算4-15、一直径为1cm长4cm的钢制圆柱形肋片，初始温度 为25C,其后，肋基温度突然升高到 200C,同时温度为25C 的气流横向掠过该肋片，肋端及两侧的表面传热系数均为2100W/（m K）。试将该肋片等分成两段（见附图），并用有 限差分法显式格式计算从开始加热时刻

17、起相邻4个时刻上的温度分布（以稳定性条件所允许的时间间隔计算依据）。已知习 JB 4-15 IT52=43W/, a 1.333 10 m /s。（提示：节点4的离散方 程可按端面的对流散热与从节点 3到节点4的导热相平衡这一 条件列出）。解：三个节点的离散方程为：节点2：k k2t 1 t 2 dx/24节点3：x/2节点4：x/2tk3tk2d2dx h tf tk2x4d2tk2tk3d2d x h tf4x4d2d2h tk4tfo44tk3tk4d2c -tk12 tk2以上三式可化简为:tk2 J t1x4hcdtf3a1 -2x4hcdtk2tkJ t2xk4xh t1稳定性要求

18、1/3a2x43a2xk3xhtf4hcd4hcdtf3a1 :2x4hcdtk351.333 103 1.333 100.0220，即32.2581/3a2x4hcd 。忖代入得:4 10010.01 32.2581058.89877s0.099975 0.0124如取此值为计算步长，则：5a 1.333 108.898772 2x0.02于是以上三式化成为:4h0.2966cd4 100 8.8987732.258510 0.010.1103o2 0.2966t1 0.2966t3k 0.1103tfkk0.296&20.2966 2t0.1103tfkk0.9773t30.0227tf

19、t 4tk12tk13时间点1234020025252520025252 2003 2004 2008.89877s3a 4h1厂在上述计算中，由于 因而对节点 2出现了在a4h厂 0.1483x之值正好使x (0 cdcd0.0551时刻温度相等这一情况。如取3a 4h2x0.5cd，于是有:kk2 0.14831 0.1483t30.5t20.1483t2k 0.1483 2t4k 0.5t3kkk0.97731 3 0.0227tft 40.0551tf0.0551tftk12tk13为上值之半，则4.4485s时间点1、 2/34020025/252520025252 2003 200

20、对于相邻四个时层的计算结果如下表所示:百度文库4 2004-16、一厚为的钢板，初始温度为650C，后置于水中淬火，其表面温度突然下降为C并保持不变。试用数值方法计算中心温度下降到450 C所需的时间。已知 a 1.1610 5m2/s。建议将平板8等分，取9个节点，并把数值计算的结果与按海斯勒计算的结果作比较。解：数值求解结果示于下图中。随着时间步长的缩小，计算结果逐渐趋向于一个恒定值，当=时，得所需时间为。如图所示，横轴表示时间步长从1秒，秒，秒，秒，秒，秒的变化；纵轴表示所需的冷却时间(用对数坐标表示)。30#4-17、一火箭燃烧器，壳体内径为400mm,厚10mm,壳体内壁上涂了一层厚

21、为2mm的包裹层。火箭发动时，推进剂燃烧生成的温度为3000 C的烟气，经燃烧器端部的喷管喷住大气。大气温度为30C。设包裹层内壁与燃气间的表面传热系数为2500 W/，外壳表面与大气间的表面传热系数为350W /(m K)，外壳材料的最高允许温度为1500C。试用数值法确定：为72使外壳免受损坏，燃烧过程应在多长时间内完成。包裹材料的=W/, a=2 10 m /s。解：采用数值方法解得 420s。4-18、锅炉汽包从冷态开始启动时，汽包壁温随时间变化。为控制热应力，需要计算汽包内壁的温度场。试用数值方法计算：当汽包内的饱和水温度上升的速率为1C /min,3 C /min时，启动后10mi

22、n,20min,及30min时汽包内壁截面中的温度分布及截面中的最大温差。启动前， 汽包处于100C的均匀温度。汽包可视为一无限长的圆柱体，外表面绝热，内表面与水之间的对流换热十分强烈。汽包的内径R .9m,外半径R21 -01m,热扩散率a 9.98 10 6m2/s。解：数值方法解得部分结果如下表所示。、汽包壁中的最大温差，K/启动后时间，min温升速率，K/mi n131020634-19、有一砖墙厚为.3mC 1.05 10 J/(m .K)室内温度为 t1 20 C百度文库2h=6W/(m K)。起初该墙处于稳定状态，且内表面温度为15C。后寒潮入侵，室外温度下 降为tf2 10 C

23、,外墙表面传热系数h2 35W/(m K)。如果认为内墙温度下降C是可 感到外界温度起变化的一个定量判据，问寒潮入侵后多少时间内墙才感知到？ 解：采用数值解法得 t=7900s。4-20、一冷柜，起初处于均匀的温度(20C)。后开启压缩机，冷冻室及冷柜门的内表面温度以均匀速度18 C /h下降。柜门尺寸为1.2m 1.2m。保温材料厚8cm,=。冰箱外表面包裹层很薄，热阻可忽略而不计。 柜门外受空气自然对流及与环境之间辐射的加热。自然对流可按下式计算：h 1.55 t/H 1/4W/(m2.K)其中H为门高。表面发射率0.8。通过柜门的导热可看作为一维问题处理。试计算压缩机起动后/2h内的冷量

24、损失。解：取保温材料的4-21、一砖砌墙壁，厚度为104 J / m3 K240mm,=4，用数值计算方法得冷量损失为5.97 10 J。31800kg/m ,c 0.88J / kg.K 。设冬天室外温度为 24h内变化如下表所示。室内空气温度 12传热系数为10W/(m .K)，内墙为及墙壁中心处温度随时间的变化。取15 C且保持不变；外墙表面6W/(m2.K)。试用数值方法确定一天之内外墙，内墙 1h。设上述温度工况以 24h为周期进行变化。时刻/h0： 001: 002 ：003： 004： 005： 006: 007:00& 009： 0010：0011 ：00温度/0C时刻12：1

25、3：14：15：16：17：18 :19：20：21:22:23:/h00000000000000000000000C温度/0c解：采用数值解法得出的结果如下表所示。时刻/h 0环境温度/0c25外墙温度/0C墙壁中心温度0C/ C内墙温度/0C时刻/h 910 11121314151617环境温度/0C-7-1外墙温度/0C墙壁中心温度内墙温度/0C时刻/h 181920212223环境温度/0C 外墙温度/0C 墙壁中心温度内墙温度/0C 多维稳态导热问题4-22、如附图所示，一矩形截面的空心电流母线的内外表面分别与温度为tfi，tf2的流体发生对流换热，表面传热系数分别为hi，h2，且各

26、自沿周界是均匀的，电流通过壁内产生均匀热源 。今欲对母线中温度分布进行数值计算，试：（1）划出计算区域（2）对该区域内的温度分布列出微分方程式及边界条件；（3）对于图中内角顶外角顶及任一内部节点列出离散方程式（x y），设母线的导热系数为常数。4-23、一个长方形截面的冷空气通道的尺寸如附图所示。假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。试用数值方法计算下列两种情况下通道壁面的温度 分布及每米长度上通过壁面的冷量损失：（1）内外壁分别维持在 10C及30C（2）内外壁与流体发生对流换热，且有tf1 10 C, X 20W/（m.K） , tf2 30 C,2h24W/（

27、m .K）。解：此题应采用计算机求解。如有墙角导热的热点模拟实验设备，则计算参数（如h, t及网格等）可以取得与实验设备的参数相一致，以把计算结果与实测值作比较。根据对称性，取1/4区域为计算区域。数值计算解出，对于给定壁温的情形，每米长通道的冷损失为，对于第三类边界条件为（取壁面导热系数0.53W/ m K ）。内外表面为给定壁温时等温线分布如下图所示。第三类边界条件的结果定性上类似。4-24、为了提高现代燃气透平的进口燃气温度以提高热效率，在燃气透平的叶片内部开设有冷却通道以使叶片金属材料的温度不超过允许值，为对叶片中的温度分布情况作一估算，把附图a所示的截片形状简化成为附图b所示的情形。

28、已知T0 17OOK,h01000W/(m2.K)， Ti 400k, h 250W/(m2.K)。试计算：(门截面中最高温度及其位置；(2)单位长度通道上的热量。解：根据对称性选择1/4区域为计算区域，采用60 70网格,综合分析与分析、论述题取壁面15W/ m K时得单位长度的传热量为，等温线 分布如图所示。截面中最高温度发生在左上角，该处温度为 0C 。4-25、工业炉的炉墙以往常用红砖和耐火砖组成。 由于该两种材料的导热系数较大，散热损失较严 重，为了节省能量，近年来国内广泛采用在耐火砖 上贴一层硅酸纤维毡，如附图所示。今用以下的非 稳态导热简化模型来评价黏贴硅酸纤维毡的收益： 设炉墙

29、原来处于与环境平衡的状态，0s时内壁表面突然上升到 550 C并保持不变。这一非稳态导 热过程一直进行到炉墙外表面的对流，辐射热损失 与通过墙壁的导热量相等为止。在炉墙升温过程中外表面的总表面传热系数由两部分组成, 即自然对流引起的部分1/3 hc W/ m2.K1.12 tw 0 c tf 0c及辐射部分TwTf /2dpdz其中：tw,Tw为外表面温度，tf,Tf为内表面温度，1 240mm, 2 240mm, 3 40mm。 为简化计算，设三种材料的导热系数分别为1 j6W/, 2 .8w/, 3 0.04 W/。试计算每平方炉墙每平方面积上由于粘贴了硅酸纤维毡而在炉子升温过程中节省的能

30、量。解：采用数值计算方法，详细过程从略。4-26、空气在附图所示的一长方形截面的送风管道中作充分 发展的层流流动，其z方向的动量方程简化为而且u v 0。上式可看成是源项为方程并计算f,Re之值。f为阻力系数, 值，并按玄巾=及1两种情形计算。Re为特征长度为当量直径 De。计算时可任取一个dpdz解：假设壁温为常数，则不同a/b下换热充分发展时的fRe及Nu数的分析解为:a/bNufRe57162ti4-27、一家用烤箱处于稳定运行状态，箱内空气平均温度2传热系数hi40W/(m .K)。外壁面与20C的周围环境间的表面传热系数ho 10 W/(m .K)。烤箱保温层厚 30mm,.3w/,

31、保温层两侧的护板用金属制成且很薄，分析中可不予考虑，然后，突然将烤箱调节器开大，风扇加速，内壁温度突然上升到185 C,设升温过程中烤箱外壁面与环1/4计算，环境温度155 C,气体与内壁间的表面境间的表面传热系数可用ho ctw tftf仍保持为20C, tw为烤箱外壁面温度，c之值与运行时 样。试确定烤箱内壁温度跃升后到达新的稳定状态所需时间。殆&解：需采用数值方法求解，过程从略。小论文题目4-28、一厚为的钢管，初始温度为16C。其后，温度为572 C的液态金属突然流过管内，并2经历了 10s。液态金属与内壁面间的表面传热系数h=W/(m K)。钢管可以按平壁处理，1/3其外表面的散热由

32、对流及辐射两条路径，并分别可按hw W/m-K1.2 t Oc及hr 4 0Tm计算，Tm Tf Tv /2，周围环境温度tf = 20C。试用有限差分法确定在3液态金属开始流入后的18s时截面上的温度分布。已知钢管的41W/,7530kg/m ,c=536J/。解：在钢管壁厚方向上取 27个点，以内壁为坐标原点，沿着壁厚方向为x正方向，数值计算结果如下。位置/cm温度/0c位置/cm0dpdz的一常物性导热方程。试用数值方法求解这一温度/C位置/cm温度/C用图形表示如下140 -H 壬口1円Amo03Q.Q19* 黑口an毎向聖的序務X/cm4-29、为对两块平板的对接焊过程（见附图a）进

33、行计算，对其物理过程作以下简化处理：钢板中的温度场仅是 x及时间 的函数；焊枪的热源作用在钢板上时钢板吸收的热流密度3r2/re2q Xqme，匚为电弧有效加热半径，qm为最大热流密度；平板上下表面的散热可用q ht tf计算，侧面绝热；平板的物性为常数，熔池液态金属的物性与固体相同；固 体熔化时吸收的潜热折算成当量的温升值，即如设熔化潜热为L，固体比热容为C,则当固体达到熔点ts后要继续吸收相当于使温度升高（ L/C ）的热量，但在这一吸热过程中该温度不 变。这样，附图a所示问题就简化为附图 b所示的一维稳态导热问题。试：（1 ）列出该问题的数学描写；（2）计算过程开始后内钢板中的温度场，设在开始的内有电弧的加热作用。已qm 5024 104W/m2h=W/(m2.K)378kg/m ,c 67J / kg.K 丄=255kJ/kg, ts 1485 C, H=12cm,re .71cm解：取初始温度与环境温度均为20 Co该问题的数学描写为：t2tq x 2h t tf /2c xc 0x0ttj0 0 x H ;t0xx500;t0xxH0。