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

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

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

3、用前六项之和来替代）可能引起的误差。解：ttan t = Bi不同Bi下前六个根如下表所示：Bia 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.2X =、Bi=0.1Bi=1Bi=10第一项的值0.948790.629450.11866前六和的值0.951420.643390.12248比值0.997240.9783

4、30.96881Fo=0.2X = 0Bi=0.1Bi=1Bi=10第一项的值0.996620.965140.83889前六项和的值0.9940.950640.82925比值1.0021.015251.01163Bi=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.927910.76851比值1.001771.009781.005984

5、-2、试用数值计算证实，对方程组X1 +2x2 -2x3 =1、« x1 +x2 + x3 = 3，2x1 +2x2 + x3 =5用高斯-赛德尔迭入法求解，其结果是发散的，并分析其原因。解：将上式写成下列迭代形式x1 =1/2(5-2x2 -x3 J<x2 =1/ 2(1 +2x3 -x1 ).x3 3 x1 一 x2假设x2,x3初值为0,迭代结果如下：迭代次数01234x102.52.6252.093752.6328125x20-0.750.4375-1.1718751.26171825x301.25-0.06252.078125-0.89453125显然，方程迭代过程发

6、散因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值 代数和。-赛德尔迭代法计算4-3、试对附图所示的常物性，无内热源的二维稳态导热问题用高斯 t1，t2，t3，t4 之值。解：温度关系式为:ti =1/4d +t3 +40+30),t2 =1/4(ti +t4 +20 +30)，t3 =1/4t1 t4 30 15t4 =1/4t2 t3 10 5开始时假设取 t*0)=t20=20 c； t30)=t40)=15 c得迭代值汇总于表迭代次数0202015126.2522.812521.5625228.5937523.35937522.1093751514.843

7、7515.1171875345628.867187528.9355425828.9526356528.956908923.4960937523.5302712923.5388178223.5409544622.2460756522.2802712922.2888178222.29095544515.1855425815.2026356515.2069089115.20797723习题用图其中第五次与第六次相对偏差已小于10 -迭代终止。4-4、试对附图所示的等截面直肋的稳态导热问题用数值方法求解节点 2,3 的温度。图中t0 =85°C,tf =25°C,h =3CW/(m

8、2.K) E0, f,'，.肋身 H=4cm,纵2剖面面积Al =4cm，导热系数九=20W/(m.K)。解：对于2点可以列出:节点2:2 2h 二 x(t1 -12) = 0;xxt2 -t3节点3:由此得:x2.x、h(tf -3) 2h (tf -t3)=02-t2t3 - t2t2t32h xH2，、c，、h ，、(t -t2)=0t2- t3 (t f- t3 )2 2h xH2x2 h2(tf -t3)=0t3t2-tfh x22 、.h.：x230 0.022L t2 0.12tf=0.06t2 =谪 20 Mo.01，于是有：2+0.12,t230/20 tf 0.03

9、tft2 1.5tf 0.03tft2 1.53tf2.532.53，代入得：5.3636t2 =2.533 t2 1.53tf0.3036tf_ 2.53tf 1.8336tf4.3636,13 二二1 30/20 0.03t21.53tf2.12t2 =L 0.12tf2.534.3636t 2 =2.533 1.8336tf t2t2t32.53 85 1.8336 25215.05 45.84=59.79 与 59.8 C4.36364.363659.8 1.53 252.53 =38.75 三 38.8 C离散方程的建立4-5、试将直角坐标中的常物性无内热源的二维稳态导热微分方程化为

10、显式差分格式，并指 出其稳定性条件(2手')。解：常物性无内热源二维非稳态方程微分方程为扩散项取中心差分，非稳态项取向前差分:.t；.i -2t；七t厂）中心-2tn +6；二a 七IAx所以有ti 1 nkkkkt i，j t i，j . t i 1, j t i，j .r 1 r7：k 1krj - -r)iTci-1-ji稳定性条件Fo x FO y 1/24-6、极坐标中常物性无内热源的非稳态导热微分方程 为笈j d2t + 1 笈1 ctir2 -9a2 r mr2 /2 ,试利用本题附图中的符号，列出节点（ i,j）的差分方程式。解：将控制方程中的各阶导数用相应的差分表示式

11、代 替，可得：h尸tk t, j/2tkt, j +t1 1 tk一匕j JKt糙甲Y曲舞即制巴附'=a + + At（ir2rj2Arr 2j,也可采用热平衡法。对于图中打阴影线的控制容积写出热平衡式得：对等式两边同除以rj ；：r并简化，可以得出与上式完全一样相同的结果。4-7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷却，底面可以认为是绝热的。为用数值法确定冷却过程中柱体温度的变化，取中心角为1rad的区域来研究（如本题附图 所示）。已知柱体表面发射率，自然对流表面传热系数，环境温度，金属的热扩散率，试列出图中节点（1, 1）, （M,1）（M,n）及（M,N）的离散方程式。

12、在 r及z方向上网格是各自 均分的。解：应用热平衡法来建立四个节点点离散方程。节点（1, 1）：k kI 1,2 I 1,1z1色j心,1 -匕肉Az p伍r2 AzyZ驾f 1用甥_ L /l尸尸 C I _2l2jj 12 2)18 2 J Att m 工1 -t m,1 rzz.zz.t'7 -m5 T "节点(m, n)：m 1,1kkm 十一t m,1:rm,1km,1tk tkm 1, n m, n I7r-rm/-：/ztkn 1 一 m> n飞平杯m拳U /htmlm,4-8、一个二维物体的竖直表面收液体自然对流冷却，为考虑局部表面传热系数的影响，表面传

13、热系数采用,，1 1.25h = c（t -11） 来表示。试列出附图所示的稳态无内热源 物体边界节点（M,n）的温度方程，并对如何求解这一方 程提出你的看法。设网格均分。解：利用热平衡法：W3mJ-,1.0 T40 J4m, n习名第图中财酊节点2:节点5:节点6:节点9:0.25h =C(tM, n 一tf 1tM, n -tf )0.25将h写为h =C(tM，n T K", n T )，其中tM, n为上一次迭代值，则方程即可线性化。4-9、在附图所示的有内热源的二维导热区域中， 一个界 面绝热，一个界面等温（包括节点 4）,其余两个界面与 温度为tf的流体对流换热，h均匀，

14、内热源强度为 6。 试列出节点1, 2, 5, 6, 9, 10的离散方程式。解：节 点1'号 lx '号3 4 x：y -2舛 ti =o,上.，X ：x 1以旷=0x 12 JAx 12 JAy2；r t1 15Ay )t9 15 rAx )t6 15 小 *1八八力 a u/x , c 1 +九(iy )+-Axiy?-iyh(t5-tf )=0y 12 J 句 12 J Ax 2 t26 ：x 76 y - tt5 x 56 ：y :x y =0yAxAyZx；，2.X Jy J/ 2 T ht9Tf =0Ayl2jAxl2，4122，；节点10,e件"若停卜

15、号31+23-，尸。当&X = Ay以上诸式可简化为：t5 t -i h-y tf 2 2 -h-y t L ：y2 - =0 节点 1： 52. ' f'12 y / 2c2t6 tl t3 - 4t2 y =0节点2:I'J ;2t6 ti t9 2 i 皿 tf -2 12 h-:y t5 . ：y2 - = 0节点5:'''2人t7 ti0 t5 t7 - 4t6 -N =0节点6:VU ;15 t10 2四 tf-2i 3 t9 ；：y2 =0节点9:八 I Z J 2 l?J；2t6 +t9 +tii +2 也 Itf -2

16、.'2+y 110 +Ay2 '-）=0节点i0：- l 九J n 。一维稳态导热计算4-10、一等截面直肋，高 H,厚S ,肋根温度为t0 ,流体温度为tf ,表面传热系数为 h,肋片 导热系数为九。将它均分成4个节点（见附图），并对肋端为绝热及为对流边界条件（ h同 侧面）的两种情况列出节点 2 ,3 ,4 的离散方程式。设H=45cm, 5 =10mm, h =50W /（m'.K），九=50W/（m.K）, t0 = 100 c , tf = 20 C ,计算节点 2, 3, 4的温度（对于肋端的两种边界条件）。5»4-10舒至解：采用热平衡法可列出节

17、点 2、3、4的离散方程为:节点2:-2h x t2 -tf)=0;节点3:-2h ：x t3-tf工XLX, t3 -4:;：-h：xt4-tf =0节点4：肋端绝热Ax,- t3 _ t4 I-1-hxt4-tf -h、t4-tf )=0 肋端对流4x。其中.：x=H3。将已知条件代入可得下列两方程组:肋端绝热g2.045t2 100.9=0t2 2.045t3 t4 0.9 -0 t3 -1.0225t4 0.45 =0肋端对流t3 -2.045t2 . 100.9=0t2 2.045t3 t4 0.9 -0t3 -1.0375t4 0.8 =0由此解得：肋端绝热 t2 =92.20C,

18、 t3 =87.70C, t，=86.20C；肋端对流.支，=86.2七，33.8。肋端对流换热的条件使肋端温度更接近于流体温度。4-11、复合材料在航空航天及化工等工业中日益得到广泛的应用。附图所示为双层圆筒壁， 假设层间接触紧密，无接触热阻存在。已知r1 =12.5mm,r2 = 16mm,r3 =18mm,% = 40A/m九2 = 120W/(m.K),tf 1 = 1500vv/(m.K),。,h1 =1000W/(m2.K),tf2 =60c, h2 =380W/(m2.K)。试用数值方法确定稳态时双层圆筒壁截面上的温度分布。为：肋根温度流体温度ty(x ) = - + (b-e(

19、1-x/H f $2: rt0及内热源&恒定，流体表面传热系数 h, 上O =f为常数。定义：t -tf, 二x / Ht0 -tf解：采用计算机求解，答案从略。采用热平衡法对两层管子的各离散区域写出能量方程，进行求解；如果采用Taylor展开法列出方程，则需对两层管子单独进行， 并引入界面上温度连续及热流密度连续的条件，数值计算也需分两区进行，界面耦合。截面的温度分布定性地示于上图中。4-12、有一水平放置的等截面直杆，根部温度t0 =100 C,其表面上有自然对流散热，h=c&_tf yd1/4,其中，c=1.20W/(m1.75.oC)；d 为杆直径，m。杆高 H=10c

20、m,直径 d=1cm,人= 50W/(m.K) , tm = 25c。不计辐射换热。试用数值方法确定长杆的散热量(需 得出与网格无关的解。杆的两端可认为是绝热的。解：数值求解过程略，Q=2.234W。4-13在上题中考虑长杆与周围环境的辐射换热，其表面发射率为0.8,环境可作为温度为tm的大空间，试重新计算其导热量。解：数值求解过程略，Q=3.320W 04-14、有如附图所示的一抛物线肋片，表面形线方程试：(1)建立无量纲温度 。的控制方程；(2)在无量中H 2纲参数t0 -tf= 0.01,£ =0.05上=0.1,hH=0.01下对上述控制方程进行数量计算。确定无量纲温度 O的

21、分布。解：无量纲温度方程为:_ 22 ._ _.2_d。/d0.01-20/ 5 5 1)')=0数值计算结果示于卜图中，无量纲温度从肋根的1变化到肋端的0.852。相对杵松，叫一维非稳态导热计算4-15、一直径为1cm,长4cm的钢制圆柱形肋片，初始温度为25C,其后，肋基温度突然升高到 200 C,同时温度为25 c的气流横向掠过该肋片，肋端及两侧的表面传热系数均为2100W/（m .K）。试将该肋片等分成两段（见附图），并用有限差分法显式格式计算从开始加热时刻起相邻4个时刻上的温度分布（以稳定性条件所允许的时间间隔计算依据）。已知九52 ,= 43W/（m.K） , a =1.3

22、33 :<10 m /s。（提不：节点 4 的离散方程可按端面的对流散热与从节点3到节点4的导热相平衡习4上5的引这一条件列出）解：三个节点的离散方程为:节点2:.tk1 -tk2二 d2x/2 4节点3:tk -tk,t 4 t 3x/2节点4:tk _tk,t 3 t 4Aux/2二 d2二 d2一dx h tf -tk2 - :?cd2< 4.AxP犬k 1- k2 - t 2Attk _tk t 2 t 3xh tk4 -tfo二 d x h tf -tk3 = ;?c二 d2以上三式可化简为:/ 1 八 a匚.；i'a='tk 12 =2亍t1亍，x2;

23、x24h _ a _ 1-rI,：x2a- 1t2 2 x2I Pcd J4hi tItf + 13a 1 14h. x2：'cd3a 4h-Pcdx2Pcdtk2kk2 xh t 4 =2 t 3：xhtf, 3aLi4h1 - -2- 0 .三 1/稳定性要求AxPcd,即3a + ©x24hPcd J。% = a435-5 = 32.258 101.333 m 10，代入得:3 1.333 1014 1000.0220.01m 32.258父 105 ,0.099975 0.0124=8.89877s如取此值为计算步长，则：a.：.1.333 10 上 8.89877x

24、2_20.0224h.：.=0.2966,Pcd4 100 8.89877一 _ 一 5 一一32.258 100.01=0.1103kk 1于是以上三式化成为:2 0.2966L 0.2966t30.1103tf=t 2kkk 1 10.2966t20.2966 2t40.1103tf =t 30.9773t3k 0.0227tf =t)=8.89877s时间点12340200252525 T200128.8125252A T200128.8155.8055.093A T200137.9573.6472.544A 7200143.0486.7085.303a； 4h1 -2-二 0在上述计算

25、中，由于 Ae之值正好使取Pcd因而对节点 2出现了在 丁及2 At时刻温度相等这一情况。如取 七为上值之半，则= 0.14834h；3a； 4h0.05511 -r -=0.5PcdAx Pcd,于是有:kk2 0.14831 0.1483t30品0.0551tf =t0.14832k 0.1483 2t4k 0.5t3k 0.0551tf =t kk0.9773 3 0.0227tf =t 4时间点12340200252525三20076.9125252A T200102.8632.7032.533A T200116.9842.6342.23. 4 4.4485s对于相邻四个时层的计算结果

26、如下表所示:200125.5152.5751.944-16、一厚为2.54cm的钢板，初始温度为 650 C ,后置于水中淬火，其表面温度突然下 降为93.5 C并保持不变。试用数值方法计算中心温度下降到450c所需的时间。已知a =1.16 M10m2/s。建议将平板8等分，取9个节点，并把数值计算的结果与按海斯勒 计算的结果作比较。解：数值求解结果示于下图中。随着时间步长的缩小，计算结果逐渐趋向于一个恒定值，当 T=0.00001s时，得所需时间为 3.92s。如图所示，横轴表示时间步长从1秒，0.1秒，0.01秒，0.001秒，0.0001秒，0.00001秒的变化；纵轴表示所需的冷却时

27、间(用对数坐标表示)。4-17、一火箭燃烧器，壳体内径为400mm，厚10mm,壳体内壁上涂了一层厚为2mm的包裹层。火箭发动时，推进剂燃烧生成的温度为3000c的烟气，经燃烧器端部的喷管喷住大气。大气温度为30Co设包裹层内壁与燃气间的表面传热系数为2500 W/(m.K),外壳表面与大气2间的表面传热系数为 350W/(m .K),外壳材料的最高允许温度为1500 C。试用数值法确定：为使外壳免受损坏，燃烧过程应在多长时间内完成。包裹材料的九=0.3 W/(m.K),72a=2M10 m2 /So解：采用数值方法解得 t = 420s。4-18、锅炉汽包从冷态开始启动时，汽包壁温随时间变化

28、。为控制热应力，需要计算汽包内壁的温度场。试用数值方法计算： 当汽包内的饱和水温度上升的速率为1C/min,3C/min时，启动后10min,20min,及30min时汽包内壁截面中的温度分布及截面中的最大温差。启动前，汽包处于100c的均匀温度。汽包可视为一无限长的圆柱体，外表面绝热，内表面与水之间的对流换热十分强烈。汽包的内径R1 = 0.9m,外半径R2 = 1.01m,热扩散率a =9.98M10"m2/s。解：数值方法解得部分结果如下表所示。汽包壁中的最大温差，K启动后时间，min温升速率，K/min13107.13621.41209.46328.393010.1930.5

29、763.-、4-19、有一砖墙厚为 6= 0.3m,九=0.85W/(m.K) , pc=1.05x10 J /(m .K)室内温度为 2t1 =20C, h=6W/(m .K)。起初该墙处于稳定状态，且内表面温度为15C。后寒潮入侵，室外温度下降为tf2 =10C,外墙表面传热系数 h2 =35W/(m2.K)。如果认为内墙温度 下降0.1 C是可感到外界温度起变化的一个定量判据，问寒潮入侵后多少时间内墙才感知 到？解：采用数值解法得 t=7900so4-20、一冷柜，起初处于均匀的温度(20C)。后开启压缩机，冷冻室及冷柜门的内表面温度以均匀速度18 C/h下降。柜门尺寸为1.2m M 1

30、.2m。保温材料厚8cm,九=0.02W/(m.K)。 冰箱外表面包裹层很薄，热阻可忽略而不计。柜门外受空气自然对流及与环境之间辐射的加 热。自然对流可按下式计算：1/42h =1.55 ：t/HW/(m2.K)其中H为门高。表面发射率 君=0.8。通过柜门的导热可看作为一维问题处理。试计算压缩 机起动后2h内的冷量损失。解：取保温材料的pc = 1x104J/(m3 K ),用数值计算方法得冷量损失为5.97M104J 。34-21、，砖砌墙壁,厚度为 240mm,儿=0.81W/(m.K), * =1800kg/m , c = 0.88J /(kg.K )o2 ,传热系数为10W/(m .

31、K) ,内墙为及墙壁中心处温度随时间的变化。取设冬天室外温度为 24h内变化如下表所示。室内空气温度ti = 15 C且保持不变；外墙表面2 ,6W/(m .K) o试用数值方法确定一天之内外墙，内墙Ae =1h。设上述温度工况以 24h为周期进行变化。时刻/h0: 001 : 002:003: 004: 005: 006: 007:008: 009:0010:0011 :00温度/0C-5.9-6.2-6.6-6.7-6.8-6.9-7.2-7.7-7.6-7.0-4.9-2.3时刻/h12:0013:0014:0015:0016:0017:0018:0019:0020:0021 :0022

32、:0023:00温度/0C-1.02.41.81.81.60.5-1.6-2.8-3.5-4.3-4.8-5.3解：采用数值解法得出的结果如下表所示。时刻/h 0环境温度/0C-5.9-6.2-6.6-6.7外墙温度/0C墙壁中 心温度-1.70-2.19-2.44-2.76/0C内墙温3.653.323.152.92度/0C8.998.828.738.61时刻/h环境温9101112度/0C外墙温-7-4.9-2.3-1度/0C墙壁中 心温度-3.58-3.07-1.340.78/0C内墙温2.362.703.875.320度/ C8.318.499.119.87-6.8-6.9-7.2-7

33、.7-7.6-2.85-2.93-3.01-3.26-3.672.872.812.752.592.318.588.558.528.438.2813141516172.41.81.81.60.51.874.634.154.143.976.057.957.627.627.5110.2611.2611.1011.1011.10时刻/h 181920212223环境温度/0C-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内墙温0 、度/ C10.7110.109.73

34、9.539.309.14多维稳态导热问题4-22、如附图所示，一矩形截面的空心电流母线的内外表面分别与温度为tf1,tf2的流体发生对流换热，表面传热系数分别为由2,且各自沿周界是均匀的，电流通过壁内产生均匀热源不。今欲对母线中温度分布进行数值计算，试：(1)划出计算区域(2)对该区域内的温度分布列出微分方程式及边界条件；(3)对于图中内角顶外角顶及任一内部节点列出离散方程式(Ax# &y),设母线的导热系数人为常数。4-23、一个长方形截面的冷空气通道的尺寸如附图所示。假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。试用数值方法计算下列两种情况下通道壁面的温度 分布

35、及每米长度上通过壁面的冷量损失：(1)内外壁分别维持在 10c及30 C(2)内外壁与流体发生对流换热，且有tf1=10C, hi =20W/(m2k) , tf2 =30c,2h2 =4 W/(m .K) o解：此题应采用计算机求解。 如有墙角导热的热点模拟实验设备，则计算参数(如h, At及网格等)可以取得与实验设备的参数相一致，以把计算结果与实测值作比较。根据对称性，取1/4区域为计算区域。数值计算解出，对于给定壁温的情形，每米长通 道的冷损失为39.84W ,对于第三类边界条件为30.97W (取壁面导热系数0 =0.53W/(m K内外表面为给定壁温时等温线分布如一F图所示。第三类边

36、界条件的结果定性上类似。4-24、为了提高现代燃气透平的进口燃气温度以提高热效率，在燃气透平的叶片内部开设有冷却通道以使叶片金属材料的温度不超过允许值，为对叶片中的温度分布情况作一估算，把附图 a所示的截片形状简化成为附图 b所示的情形。已知To =17OOK,ho =1000W/(m2.K) , Ti =400k,。=250 W/(m2.K) o 试计算：(1)截 面中最高温度及其位置；(2)单位长B通道上的热量。解：根据对称性选择1/4区域为计算区域，采用60 M 70网格,取壁面 九= 15W/(m K )时得单位长度的传热量为987.8W,等温线分布如图所示。截面中最高温度发生在左上角

37、，该处 温度为1419.90C 。综合分析与分析、论述题4-25、工业炉的炉墙以往常用红砖和耐火科组成。 由于该两种材料的导热系数较大，散热损失较严 重，为了节省能量，近年来国内广泛采用在耐火科 上贴一层硅酸纤维毡，如附图所示。今用以下的非 稳态导热简化模型来评价黏贴硅酸纤维毡的收益：设炉墙原来处于与环境平衡的状态，7=0s时内壁表面突然上升到 550 c并保持不变。这一非稳态导 热过程一直进行到炉墙外表面的对流，辐射热损失与通过墙壁的导热量相等为止。在炉墙升温过程中外表面的总表面传热系数由两部分组成, 即自然对流引起的部分%c)W/(m2.K )=1.12(4w 0 c - 4f )0c f

38、/3及辐射部分hr =4，Tn2,Tm = Tw Tf /2其中：tw,Tw为外表面温度，tf,Tf为内表面温度，a = 240mm力2 = 240mm/3 =40mm。 为简化计算，设三种材料的导热系数分别为% =1.6 W/(m.K),工2 = 0.8 W/(m.K), 3 =。.。4W/(m.K)。试计算每平方炉墙每平方面积上由于粘贴了硅酸纤维毡而在炉子升温过程中节省的能量。解：采用数值计算方法，详细过程从略。4-26、空气在附图所示的一长方形截面的送风管道中作充分 发展的层流流动，其 z方向的动量方程简化为dp而且u = v = 0。上式可看成是源项为 一dz的一常物性导热方程。试用数

39、值方法求解这一方程并计算f,Re之值。f为阻力系数，Re为特征长度为当量直径 De。计算时可任取一个 情形计算。dpdz值，并按a/b= 0.5及1两种解：假设壁温为常数，则不同a/b下换热充分发展时的fRe及Nu数的分析解为：a/bNufRe12.98570.53.39624-27、一家用烤箱处于稳定运行状态，箱内空气平均温度ti =155C,气体与内壁间的表面2,、传热系数hi =40W/(m .K)。外壁面与20c的周围环境间2的表面传热系数h。=10 W/(m .K)。烤箱保温层厚30mm,九= 0.03w/(m.K),保温层两侧的护板用金属制成且很 薄，分析中可不予考虑，然后，突然将

40、烤箱调节器开大，风 扇加速，内壁温度突然上升到185C,设升温过程中烤箱外壁面与环境间的表面传热系数可用h0 =ctw 一31/4计算，环境温度tf仍保持为20 C, 均为烤箱外壁面温度，c之值与 运行时一样。试确定烤箱内壁温度跃升后到达新的稳定状态 所需时间。解：需采用数值方法求解，过程从略。小论文题目4-28、一厚为2.54cm的钢管，初始温度为16 C。其后，温度为 572c的液态金属突然流过2管内，并经历了 10s。液态金属与内壁面间的表面传热系数 h=2.84W/(m .K)。钢管可以按 平壁处理，其外表面的散热由对流及辐射两条路径，并分别可按 Gw%/”厂1.2父(5% )及hr=

41、4股。丁：计算，Tm =(Tf *Tw>2 ,周围环境温度tf = 20Co试用有限差分法确定在液态金属开始流入后的18s时截面上的温度分布。 已知钢管的3九=41W/(m.K) , * = 7530kg /m , c=536J/(kg.K)。解：在钢管壁厚6向上取 27个点，以内壁为坐标原点，沿着壁厚方向为x正方向，数值计算结果如下。位置00.10.20.30.40.50.60.70.80.9/cm温度/0C216.0215.6214.6213.0210.7207.9204.6200.8196.6192.1位置/cm1.11.21.31.41.51.61.71.81.9温度/0C187

42、.3182.3177.2172.0166.9161.8157.0152.5148.2144.4位置/cm2.12.22.32.42.52.54温度/0C141.0138.1135.7133.9132.6132.0131.9用图形表示如下4-29、为对两块平板的对接焊过程（见附图a）进行计算，对其物理过程作以下简化处理：钢板中的温度场仅是 x及时间T的函数；焊枪的热源作用在钢板上时钢板吸收的热流密度a2/%2q（x）=qme1 re为电弧有效加热半径，qm为最大热流密度；平板上下表面的散热可用q =hft -tf狂算，侧面绝热；平板的物性为常数，熔池液态金属的物性与固体相同；固 体熔化时吸收的潜热折算成当量的温升值，即如设熔化潜热为L,固体比热容为C,则当固体达到熔点ts后要继续吸收相当于使温度升高（ L/C）的热量，但在这一吸热过程中该温度不 变。这样，附图a所示问题就简化为附图b所示的一维稳态导热问题。试： （1）列出该问题的数学描写；（2）计算过程开始后 3.4s内钢板中的温度场，设在开始的0.1s内有电弧的加422热作用。已知：qm =5024 10 W/m h=12.6 W/（m.K），九= 41.9W/（m.K）,=7800kg/m3,c =670J "kg.K）L=255kJ/kg, ts =1485 C,