第四版传热学习题解答

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

2、代数方程时是否一定可以得到收敛德解？不能得出收敛的解- 3tn' 5tn d -'tn 2n,i时是否因为初场的假设不合适而造成？;:t8 .有人对一阶导数 ：X你能否判断这一表达式是否正确，为什么?一般性数值计算4-1、采用计算机进行数值计算不仅是求解偏微分方程的有力工具，而且对一些复杂的经验 公式及用无穷级数表示的分析解，也常用计算机来获得数值结果。试用数值方法对Bi=0.1,1,10的三种情况计算下列特征方程的根(n =1,2，6)：x . Bitan 七=-,n =1,2,3na.Fo = = _0.2并用计算机查明，当 6 时用式(3-19)表示的级数的第一项代替整个

3、级数(计算中用前六项之和来替代)可能引起的误差。解：*n tan *n = Bi，不同Bi下前六个根如下表所示：Bi(1 12H 3H 4H 5H 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.2x = 5Bi=0.1Bi=1Bi=10第一项的值0.948790.629450.11866前六和的值0.951420.643390.12248

4、比值0.997240.978330.96881Fo=0.2X = 0Bi=0.1Bi=1Bi=10第一项的值0.996620.965140.83889前六项和的值0.9940.950640.82925比值1.0021.015251.01163Fo=0.24X = §Bi=0.1Bi=1Bi=10第一项的值0.945130.611080.10935前六项的值0.946880.61980.11117比值0.998140.986940.98364Fo=0.24 X = 0Bi=0.1Bi=1Bi=10第一项的值0.992770.936980.77311前六项和的值0.991010.9279

5、10.76851比值1.001771.009781.005984-2、试用数值计算证实，对方程组又 + 2x2 2X3 = 1'X1 + x2 + x3 = 3.2x1 +2x2 +x3 =5用高斯-赛德尔迭於法求解，其结果是发散的，并分析其原因。解：将上式写成下列迭代形式Xi =1/2(5 2X2 X3 )灯2 =1/2(1 +2x3 -X1 )，X3 = 3 Xi X2假设X2，X3初值为0,迭代结果如下：迭代次数01234X102.52.6252.093752.6328125X20-0.750.4375 -1.1718751.26171825X301.25-0.06252.078

6、125-0.89453125显然，方程迭代过程发散因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值 代数和。4-3、试对附图所示的常物性，无内热源的二维稳态导热问题用高斯-赛德尔迭代法计算打,灼3,七4之值。解：温度关系式为：b =1/她 +t3 +40 + 30 jt2 =1/4(54 +20+30),t3 =1/4t1 t4 30 15t4 =1/4t2 t3 10 5开始时假设取 t11220 c ； t )t )15 c得迭代值汇总于表迭代次数020201515126.2522.812521.562514.84375228.5937523.35937522

7、.10937515.1171875328.867187523.4960937522.24607565159355425823.5302712922.2802712915.20263565528.9526356523.5388178222.2888178215.20690891628.956908923.5409544622.29095544515.20797723其中第五次与第六次相对偏差已小于10圣迭代终止。4-4、试对附图所示的等截面直肋的稳态导热问题用数值方 &山法求解节点 2 ,3 的温度。图中 ，,0 0 2/ J L *2> 3 | A.tt

8、t0=85C，tf =25C，h=3(W/(m .K).肋高心冲纵 顼 一 剖面面积Al = 4cm2,导热系数L=20W/(m.K)。*解：对于2点可以列出：习JSi恤，、.，、. h_ 2hx(t f) = 0；节点2-X -Xto tax23 h(tf -t1) 2h (tf -t3) =0节点 3:A/22。由此得：ti,2h ：x2 ,，、八-t2 '七3 -t? '(ti "2) =0aS6hAx h2 人 、ct 2 - 七3 (t f - t3 )(t f - t3) = 0/uO/uOt2t3,，2hAxH 2 1 /L 2hAxH 2 22h ：x

9、 _ 30 0.02:. 一 20 0.01= 0.06于是有:t2t1t20.12tf2 0.12t3t2 - 30/20 tf 0.03tf1 30/20 0.03t2 1.5tf 0.03tf t2 1.53tf=2.532.532.12t2 =t上g2.530.12tf5.3636t2 =2.53t t2 1.53tf0.3036tf4.3636t2 =2.53t 1.8336tf t22.53tf 1.8336tf4.3636t22.53 85 1.8336 25215.05 45.84=59.79 三 59.8 C4.36364.3636t359.8 1.53 25 =38.75

10、三 38.8 C2.53离散方程的建立4-5、试将直角坐标中的常物性无内热源的二维稳态导热微分方程化为显式差分格式，并指 出其稳定性条件(皈# Ay)。解：常物性无内热源二维非稳态方程微分方程为扩散项取中心差分，非稳态项取向前差分:t； 1弋 =aAt所以有1' 2t' +1' ti 2t' +1' 蚓书2tn Tln= 土 In 千一 2tn Tin3x2稳定性条件tn.i 1.tn =a7LnJ.产 1/2习射陌 根*等中斡阿怡4-6、极坐标中常物性无内热源的非稳态导热微分方程为，约 1 a 1 c2t =a z- + -与y槌r2r济r2 &quo

11、t;2,试利用本题附图中的符号，列出节点（ i,j）的差分方程式。.:tCT解：将控制方程中的各阶导数用相应的差分表示式代替，可得：k i kkk kkk【1, jTt,j1 t t, j+2t t, jt j 11 【t,j¥l1,j A=a 2-rrj2、：r1-2"r jti，jk kk-2t i, j r ”, j也可采用热平衡法。对于图中打阴影线的控制容积写出热平衡式得：k 1kkkkr 土r iMc，j 一，j = 心门 一、j .汀 i 1, j江".,t：n-t'jj 2r对等式两边同除以r"；：；：r并简化，可以得出与上式完全一

12、样相同的结果。tL.-1"r-tk . rrj：j-三j 24-7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷 却，底面可以认为是绝热的。为用数值法确定冷却过程中柱 体温度的变化，取中心角为1rad的区域来研究（如本题附图 所示）。已知柱体表面发射率，自然对流表面传热系数，环 境温度，金属的热扩散率，试列出图中节点（1, 1）, （M,1）（M,n）及（M,N ）的离散方程式。在 r及z方向上网格是各自 均分的。解：应用热平衡法来建立四个节点点离散方程。节点（1, 1）:.k , k2, k , k-t 1,2 -1 1,1 1 r I ) -+ - t 2,1 -1 1,1”A

13、z2 12 J"Ar仔r仔r2 Az'tTj皆f%附图I I = ic l22j 、82 J Attk ” Tk ,>> tk tk-z z tk.L mJJ r m,1_±_± . ：. ' m 11_1 m,1 .二 . :. ' m,2节点（m, n）:rrm22r 繇 224-8、一个二维物体的竖直表面收液体自然对流冷却，为考虑局部表面传热系数的影响，表1.25面传热系数采用h = c(t -&)来表示。试列出附图所示的稳态无内热源物体边界节点(M,n)的温度方程，并对如何求解这一方程提出你的看法。设网格均分。解

14、：利用热平衡法：0.25h =C(tM , n -tf Jm , n Tf )0.25将h写为h =C(tM，n _tf 2m，n _tf )，其中tM , n为上一次迭代值，则方程即可线性化。习日4丹尉U4-9、在附图所示的有内热源的二维导热区域中， 一个界 面绝热，一个界面等温(包括节点 4),其余两个界面与 温度为tf的流体对流换热，h均匀，内热源强度为 $。 试列出节点1, 2, 5, 6, 9, 10的离散方程式。解：节 点1:tV f r"y -捉&土 =0土q包一七或包 , 土也*l.*y =0节点2:节点5:节点6:Ax(2 JAx(2 JAy2；'土

15、 心三，W ： . x ：y - ：yh t5-tf =0 y 2y 2x 2, ；冰,Ly-°：匚乂-:工y 广 xy= 0yAxAyAx、t5 t9 fAx)t0 T9 fAy " 1 相Ax Ay ,、k. + k. l + AxAyW .+ h(t9 tf )=0节点 9:也 一2 ) 芯 12 J 4V22 J-E*板项-由 hw-tf =0节点 10：* I 2 .jAx 12 jAy2当板=Ay以上诸式可简化为:t5 t22 tf _2 2 2 t1 1：y2 =0节点1：节点2：节点5：节点6:222t6 t1 t3 -4t2y2 =0vJ ；2t6 t1

16、t9 2 亶 tf -2 2 4 t5：y2 - =0V. X ) k )-J2t7 t10 t5 t7 -4t6 Ny =0节点9： 5”+普卜4+却9E&）=0；2t6t9 til2 i 空 tf 2 i 2 H tioy2 =0节点 10：H I A ） E。一维稳态导热计算4-10、一等截面直肋，高 H,厚8 ,肋根温度为to，流体温度为tf，表面传热系数为 h,肋片 导热系数为九。将它均分成4个节点（见附图），并对肋端为绝热及为对流边界条件（ h同 侧面）的两种情况列出节点 2 ,3 ,4 的离散方程式。设H=45cm, & =10mm，h =5网/前.心，*=50w

17、/（m.K）, to =100 c, tf =20 c,计算节点 2, 3, 4的温度（对于肋端的两种边界条件）。解：采用热平衡法可列出节点2、3、4的离散方程为:节点节点't 土 /，什3 -12 '2h" t2 tf2：xx't2 -匕顶 / 14 -t3 '2h ：x t3tf3：x节点x t3 - t43 -h ：x t4 -tf ）：=04:肋端绝热Ax,十 I t3 - t4 ）、打5、力：-h x n -tf -卜 a -tf）=°肋骊对流Ax。其中3。将已知条件代入可得下列两方程组:肋端绝热t2.045t2 10°.

18、9=0t2 -2.045t3 t4 0.9 =0t3 -1.0225t4 0.45 =0肋端对流t2.045t2 100.0t2 -2.045t3 - t4 0.9 =0t3 -1.0375t4 0.8 =0由此解得：肋端绝热 t2 =92.20C ,乜=87.70。，t4 =86.20C ；肋端对流 &=91.50。,，86.& , %=83.此。肋端对流换热的条件使肋端温度更接近于流体温度。4-11、复合材料在航空航天及化工等工业中日益得到广泛的应用。附图所示为双层圆筒壁,假设层间接触紧密，无接触热阻存在。已知 r1 =12.5mm,r2 = 16mm, e =18mm,兀

19、1 = 40 w/(m K)舄2 =120W /(m.K), tf 1 = 150 h =100CW/(m2.K),tf2 =60 c, h? =380W/(m2.K)。试用数值方法确定稳态时双层圆 筒壁截面上的温度分布。解：采用计算机求解，答案从略。采用热平衡法对两层管子的各离散区域写出能量方程，进行求解；如果采用Taylor展开法列出方程，则需对两层管子单独进行， 并引入界面上温度连续及热流密度连续的条件，数值计算也需分两区进行，界面耦合。截面的温度分布定性地示于上图中。4-12、有一水平放置的等截面直杆，根部温度t。=100 C,其表面上有自然对流散热，h=c(ttf yd】/4,其中，

20、c=1.20W/(m175.°C)；d 为杆直径，m。杆高 H=10cm ,直径 d=1cm,赤=50W/(m.K) , t*-25C。不计辐射换热。试用数值方法确定长杆的散热量(需 得出与网格无关的解。杆的两端可认为是绝热的。解：数值求解过程略，Q=2.234W。4-13在上题中考虑长杆与周围环境的辐射换热，其表面发射率为0.8,环境可作为温度为too的大空间，试重新计算其导热量。解：数值求解过程略， Q=3.320W。4-14、有如附图所示的一抛物线肋片，表面形线方程为：y (x )=。+ (b - e (1 - x/ H f 电2匕b肋根温度0及内热源$恒定，流体表面传热系数

21、h, t -tf 。=流体温度tf为常数。定义：to试：(1)建立无量纲温度 O的控制方程；>H 2eb=0.01,0.05,纲参数t0 -tfHH算。确定无量纲温度 &的分布。= x/H(2)在无量hH= 0.1,0.01下对上述控制方程进行数量讨d20/d2E+0.01 2。/(5 + 5(1 匚)=0 山解：无量纲温度方程为：')。数值计算结果示于下图中，无量纲温度从肋根的1变化到肋端的0.852。习JB4以附甜.k i.k.2 .k . kt i t 2 顷d ) t 3 T 2/ + 九x/2 、 4 jAx节点3:节点tk _tk I 4 I 3/.x/24:

22、二 d2tk-二 : d2x/2以上三式可化简为:k 1t 2h(tk4-tfJd24二 d x h tf -tk3 =七 : d2,k虫t 2k2t 3=2 W t" 土 t3 x2 1x2 3=TWx )t +2Wt + t2 2W44h *Il fcd f4h ":cd3aA t 4hA 7x2Pcd j3a 4htk2、tk3cd一维非稳态导热计算4-15、一直径为1cm,长4cm的钢制圆柱形肋片，初始温度 为25 C,其后，肋基温度突然升高到 200 C,同时温度为25 C 的气流横向掠过该肋片，肋端及两侧的表面传热系数均为 100W/（m2.K）。试将该肋片等分

23、成两段（见附图） ，并用有 限差分法显式格式计算从开始加热时刻起相邻4个时刻上的温度分布（以稳定性条件所允许的时间间隔计算依据）。已知舄=43W/（m.K） , a =1.333=<10Mm2/s。（提示：节点 4 的离 散方程可按端面的对流散热与从节点3到节点4的导热相平衡这一条件列出）。解：三个节点的离散方程为:节点2:2 xh tk4 =2 tk3xhtf,3a." 4h."1 - 0: _1/稳定性要求&Pcd,即3a 4hI Ax2Pcd八435=32.258 1054 1001.333x10,代入得:=8.89877s0.099975 0.0124

24、3x1.333x102+ 50.0220.01 x 32.258k 10"如取此值为计算步长，贝U：a."1.333 10 直 8.898774h"0.2966-,Pcd2 0.2966t1 0.2966t3k 0.1103tf =tk 120.0224 100 8.898775 = 0.110332.258 100.01于是以上三式化成为:kkk 10.296&20.2966 2t40.1103tf=t 30.9773t3k 0.0227tf =tk4：=8.89877s时间点12340200252525 T200128.8125252A T200128

25、.8155.8055.093A T200137.9573.6472.544 T200143.0486.7085.30/ 3az 4h-1 - -: = 0在上述计算中，由于 Ae之值正好使皈Pcd,因而对节点 2出现了在 At及2 At时刻温度相等这一情况。如取At为上值之半，则a. ： .4h3a " 4h ：厦=0.1483 ,冠=0.0551 ,"寸一岳=0.5,于是有：kkk 12 0.14831 0.148&3 0星 0.0551tf =t 2 10.1483t20.1483 2t40.%0.0551tf =t 30.9773 k3 0.0227tf =也

26、对于相邻四个时层的计算结果如下表所示：3 =4.4485s)时间点12340200252525 720076.9125252A 1200102.8632.7032.533 T200116.9842.6342.234 7200125.5152.5751.944-16、一厚为2.54cm的钢板，初始温度为 650 C ,后置于水中淬火，其表面温度突然下 降为93.5 C并保持不变。试用数值方法计算中心温度下降到450 C所需的时间。已知a =1.16*10'm2/s。建议将平板8等分，取9个节点，并把数值计算的结果与按海斯勒 计算的结果作比较。解：数值求解结果示于下图中。随着时间步长的缩小

27、，计算结果逐渐趋向于一个恒定值，当A=0.00001s时，得所需时间为 3.92s。如图所示，横轴表示时间步长从1秒，0.1秒，0.01秒，0.001秒，0.0001秒，0.00001秒的变化；纵轴表示所需的冷却时间(用对数坐标表示)。4-17、一火箭燃烧器，壳体内径为400mm，厚10mm,壳体内壁上涂了一层厚为2mm的包裹层。火箭发动时，推进剂燃烧生成的温度为3000 C的烟气，经燃烧器端部的喷管喷住大气。大气温度为30C。设包裹层内壁与燃气间的表面传热系数为2500 W/(m.K),外壳表面与大气间的表面传热系数为 350W/(m2.K)，外壳材料的最高允许温度为1500 C。试用数值法

28、确定：为使外壳免受损坏，燃烧过程应在多长时间内完成。包裹材料的舄=0.3 W/(m.K),a=2x10'm2/s。解：采用数值方法解得 t =420s。4-18、锅炉汽包从冷态开始启动时，汽包壁温随时间变化。为控制热应力，需要计算汽包内 壁的温度场。试用数值方法计算：当汽包内的饱和水温度上升的速率为1C /min,3 C /min时，启动后10min,20min,及30min时汽包内壁截面中的温度分布及截面中的最大温差。启动前， 汽包处于100C的均匀温度。汽包可视为一无限长的圆柱体，外表面绝热，内表面与水之间的对流换热十分强烈。汽包的内径R = 0.9m,外半径R2 =dm,热扩散率

29、 .一 .-62 a = 9.98 乂10 m /s。解：数值方法解得部分结果如下表所示。汽包壁中的最大温差，K启动后时间，min温升速率，K/min13107.13621.41209.46328.393010.1930.574-19、有一砖墙厚为 &=0.3m，丸=0.85W/(m.K) , % = 1.05 乂 106 J/(m3.K)室内温度为 G =20c, h=6W/(m2.K)。起初该墙处于稳定状态，且内表面温度为15C。后寒潮入侵，室外温度下降为tf2 =T°c，外墙表面传热系数h2 =35W/(m2.K)。如果认为内墙温度 下降0.1 C是可感到外界温度起变化

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

31、 K )，用数值计算方法得冷量损失为5.97x104J。34-21、一砖砌墙壁，厚度为 240mm，丸=0.81W/(m.K), P = 1800kg / m , c = 0.88J/(kgK )。设冬天室外温度为 24h内变化如下表所示。室内空气温度 =15c且保持不变；外墙表面传热系数为10W/(m2.K)，内墙为6W/(m2.K)。试用数值方法确定一天之内外墙，内墙 及墙壁中心处温度随时间的变化。取At =1h。设上述温度工况以 24h为周期进行变化。时刻10：11 ：/h0： 001 : 002 : 003: 004: 005: 006: 007: 008: 009: 000000温度

32、/°C-5.9-6.2-6.6-6.7-6.8-6.9-7.2-7.7-7.6-7.0-4.9-2.3时刻12:13:14:15:16:17:18 :19:20:21 :22:23:/h000000000000000000000000温度/0C-1.02.41.81.81.60.5-1.6-2.8-3.5-4.3-4.8-5.3解：采用数值解法得出的结果如下表所示。时刻/h012345678环境温度/0C-5.9-6.2-6.6-6.7-6.8-6.9-7.2-7.7-7.6外墙温度/0C-1.70-2.19-2.44-2.76-2.85-2.93-3.01-3.26-3.67/0C

33、内墙温，0八 度/ C3.658.993.328.823.158.732.928.612.878.582.818.552.758.522.598.43时刻/h910111213141516环境温度/0c-7-4.9-2.3-12.41.81.81.6外墙温度/0C墙壁中心温度-3.58-3.07-1.340.781.874.634.154.14/0C2.362.703.875.326.057.957.627.62内墙温度/0C8.318.499.119.8710.2611.2611.1011.10墙壁中心温度2.318.28170.53.977.5111.10时刻/h 181920212223

34、环境温广 °C -1.6-2.8-3.5-4.3-4.8-5.3外墙温度/0C3.061.340.36-0.22-0.87-1.29墙壁中心温度6.095.735.054.664.213.93/0C内墙温度/0C10.7110.109.739.539.309.14多维稳态导热问题4-22、如附图所示，一矩形截面的空心电流母线的内外表面分别与温度为tf1,tf2的流体发生对流换热，表面传热系数分别为h1，h2,且各自沿周界是均匀的，电流通过壁内产生均匀热源。今欲对母线中温度分布进行数值计算，试：(1) 划出计算区域(2) 对该区域内的温度分布列出微分方程式及边界条件；(3) 对于图中内

35、角顶外角顶及任一内部节点列出离散方程式(Ax#Ay),设母线的导热系数兀为常数。4-23、一个长方形截面的冷空气通道的尺寸如附图所示。假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。试用数值方法计算下列两种情况下通道壁面的温度 分布及每米长度上通过壁面的冷量损失：(1)内外壁分别维持在 10C及30 C(2)内外壁与流体发生对流换热，且有tf1 =10 c, h1 =20W/(m2.K) , tf2=30c, h2 =4W/(m2.K)。(B)解：此题应采用计算机求解。如有墙角导热的热点模拟实验设备，则计算参数(如h, At及网格等)可以取得与实验设备的参数相一致，以把计算结

36、果与实测值作比较。根据对称性，取1/4区域为计算区域。数值计算解出，对于给定壁温的情形，每米长通道的冷损失为39.84W ,对于第三类边界条件为30.97W (取壁面导热系数 九=0.53W/(m K )。内外表面为给定壁温时等温线分布如-F图所示。第三类边界条件的 结果定性上类似。4-24、为了提高现代燃气透平的进口燃气温度以提高热效率，在燃气透平的叶片内部开设有冷却通道以使叶片金属材料的温度不超过允许值，为对叶片中的温度分布情况作一估算，附图 a所示的截片形状简化成为附图 b所示的情形。已To =17OOK,ho =1000W/(m2.K) , Ti = 400k,瓦=250 W/(m2.

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

38、s时内壁表面突然上升到 550 C并保持不变。这一非稳态导 热过程一直进行到炉墙外表面的对流，辐射热损失 与通过墙壁的导热量相等为止。在炉墙升温过程中外表面的总表面传热系数由两部分组成， 即自然对流引起的部分kc)W/(m2.K )=1.12(妃°c -I Lc J/'及辐射部分hr =4Fn2,Tm = Tw Tf /2其中：tw,Tw为外表面温度，tf'Tf为内表面温度，鬲= 240mm,&2 = 240mm03 =40mm。 为简化计算，设三种材料的导热系数分别为*1 =1.6 W/(m.K),七2 =。.8 W/(m.K), ，一3 =0.04W/(m

39、.K)。试计算每平方炉墙每平方面积上由于粘贴了硅酸纤维毡而在炉子升温 过程中节省的能量。解：采用数值计算方法，详细过程从略。4-26、空气在附图所示的一长方形截面的送风管道中作充分 发展的层流流动，其 z方向的动量方程简化为/ _2_2n c w + G wdp而且u =v = 0。上式可看成是源项为 dz的一常物性导热方程。试用数值方法求解这一dp 方程并计算f,Re之值。f为阻力系数，Re为特征长度为当量直径 De。计算时可任取一个 dz 值，并按a/b= 0.5及1两种情形计算。解：假设壁温为常数，贝U不同a/b下换热充分发展时的fRe及Nu数的分析解为：a/bNufRe12.98570

40、.53.39624-27、一家用烤箱处于稳定运行状态，箱内空气平均温度G =155c, 气体与内壁间的表面传热系数加=40W/(m2.K)。外壁面与20C的周围环境间的表面传热系数ho=1OW/(m.K)。烤箱保温层厚30mm,舄=0.03w/(m.K),保温层两侧的护板用金属制成且很 薄，分析中可不予考虑，然后，突然将烤箱调节器开大，风 扇加速，内壁温度突然上升到185C,设升温过程中烤箱外壁面与环境间的表面传热系数可用加=dtw -tf 1/4计算，环境温度tf仍保持为20C,九为烤箱外壁面温度，c之值与 运行时一样。试确定烤箱内壁温度跃升后到达新的稳定状态 所需时间。解：需采用数值方法求

41、解，过程从略。小论文题目4-28、一厚为2.54cm的钢管，初始温度为16 C。其后，温度为 572 C的液态金属突然流过管内，并经历了 10s。液态金属与内壁面间的表面传热系数 h=2.84W/(m2.K)。钢管可以按 平壁处理，其外表面的散热由对流及辐射两条路径，并分别可按)=1.2乂灿% )及阵=4气3计算，Tm=(Tf +Tw2，周围环境温度tf =20C。试用有限差分法确定在液态金属开始流入后的18s时截面上的温度分布。已知钢管的次=41W/(m.K) , P = 7530kg / m' , c=536J/(kg.K)。解：在钢管壁厚方向上取 算结果如下。27个点，以内壁为坐

42、标原点，沿着壁厚方向为x正方向，数值计位置/cm00.10.20.30.40.50.60.70.80.9温度/0C216.0215.6214.6213.0210.7207.9204.6200.8196.6192.1位置/cm11.11.21.31.41.51.61.71.81.9温度/0C187.3182.3177.2172.0166.9161.8157.0152.5148.2144.4位置22.12.22.32.42.52.54/cm温度0 -141.0138.1135.7133.9132.6132.0131.9/°C用图形表示如下14 Q -口口1 amo尊eg q.qis 香致

43、口 an坷却H的魔野X/cm4-29、为对两块平板的对接焊过程（见附图a）进行计算，对其物理过程作以下简化处理：钢板中的温度场仅是 x及时间T的函数；焊枪的热源作用在钢板上时钢板吸收的热流密度 q（x） = qme°r，。）, re为电弧有效加热半径，qm为最大热流密度；平板上下表面的散热可 用q = h（ttf 十算，侧面绝热；平板的物性为常数，熔池液态金属的物性与固体相同；固 体熔化时吸收的潜热折算成当量的温升值，即如设熔化潜热为L,固体比热容为c,则当固体达到熔点匕后要继续吸收相当于使温度升高（ L/c）的热量，但在这一吸热过程中该温度不 变。这样，附图a所示问题就简化为附图b所示的一维稳态导热问题。试：（1）列出该问题的数学描写；（2）计算过程开始后 3.4s内钢板中的温度场，设在开始的0.1s内有电弧的加热作用。已知：qm=502U104W/m2 ,h=12.6 W/（m2.K）,舄=41.9W/（m.K）,=7800kg/m3,c =670J/（kgK）L=255kJ/kg, ts = 1485 C , H=12cm, re = 0.71cm。解