第四版传热学第四章习题解答讲诉

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

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

3、=bi ,不同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.2x =、bi=0.1bi=1bi=10第一项的值0.948790.629450.11866前六和的值0.951420.643390.12248比值0.997240.978330.96881fo=0.2x = 0bi=0.1bi

4、=1bi=10第一项的值0.996620.965140.83889前六项和的值0.9940.950640.82925比值1.0021.015251.01163fo=0.24 x =、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.927910.76851比值1.001771.009781.005984-2、试用数值计算证实，对方程

5、组x1 +2x2 -2x3 = 1 (x1 +x2 + x3 = 332x1 +2x2 +x3 =5用高斯-赛德尔迭法求解，其结果是发散的，并分析其原因。解：将上式写成下列迭代形式x1= 1/2(5 -2x2- x3 )%= 1/2(1 +2x3- x1卜|x3=3 x 一 x2-假设x2,x3初值为0,迭代结果如下：迭代次数01234x102.52.6252.093752.6328125乂20-0.750.4375 -1.1718751.26171825x301.25-0.06252.078125-0.89453125显然，方程迭代过程发散因为迭代公式的选择应使每一个迭代变量的系数总大于或等

6、于式中其他变量的系数绝对值 代数和。-赛德尔迭代法计算4-3、试对附图所示的常物性，无内热源的二维稳态导热问题用高斯 tl, t2,t3,t4 之值。解：温度关系式为:126.2522.812521.5625228.5937523.35937522.109375328.867187523.4960937522.24607565tlt2t3t4开始时假设取 1112 cz ； t3 t41 cz得迭代值汇总于表迭代次数0202015= 1/4(t2 +t3 +40 +30 j =1/4(t1 +t4 +20 +30)、 = 1/4(ti +t4 +30 +15 ) =1/4(t2 +t3 +10

7、 +5)1514.8437515.1171875159355425823.5302712922.28027129528.9526356523.5388178222.28881782628.956908923.5409544622.290955445其中第五次与第六次相对偏差已小于10工迭代终止。4-4、试对附图所示的等截面直肋的稳态导热问题用数值方 法求解节点 2,3 的温度。图中t =85匕tf =25th =30w/(m2.k).肋高 +4cm,纵2剖面面积al =4cm，导热系数九=20w/(m.k)。解：对于2点可以列出:、21、-2hx(ti -12)= 0

8、;节点2:xxto 一 八. ：x h(tf -t1)2h (tf -t3) = 0x2节点3:2由此得:15.2026356515.2069089115.20797723(tf -t3) =0t1,2h. ：x2 ,.、-,-% t3 - t2(t1 - t2) = 0 t2 - t3(tf y)x2 hx 2t2=t1 t32h、xh2c 2h：xh22 t3t2-tf22h x 30 0.0220 0.01= 0.06t3t2 30 / 20 t f0.03tf1 30/20 0.032.12t2 =t1t2 1.53tf2.53l 湎)1+0+比22a6 )11t2 0.12tft2

9、:于是有：2十0.12,0.12tf4.3636t2 =2.53tl 1.8336tft2t21.5tf 0.03tft2 1.53tf2.532.53，代入得:5.3636t2 =2.53t1t2 1.53tf0.3036tf2.53tf 1.8336t f4.3636t22.53 85 1.8336 254.3636215.05 45.84 =59.79 = 59.8 c4.3636t359.8 1.53 252.53=38.75 三 38.8 c离散方程的建立4-5、试将直角坐标中的常物性无内热源的二维稳态导热微分方程化为显式差分格式，并指出其稳定性条件( 2手)。解：常物性无内热源二维

10、非稳态方程微分方程为2 22日 ,+硝=a 十f a 4 24 2心七 、yx cy扩散项取中心差分，非稳态项取向前差分:tni 1 -tni a所以有j 1tn =a稳定性条件1/2tn 1 - tnl-2tntnd tnx2j习用46对四 横法行中的两府4-6、极坐标中常物性无内热源的非稳态导热微分方程为笈2 e2t1a1e2t 丁 二 a.2+ : +2t2ct(frrcrr )试利用本题附图中的符号，列出节点( i,j)的差分方程式。解：将控制方程中的各阶导数用相应的差分表示式代替，可得：k 1kkkkkkt1, j -t t, j t t, j 平 一4 t, j tt j1 tt,

11、 j 4t1, j二 a 2-匚三l匚rrj2r也可采用热平衡法。对于图中打阴影线的控制容积写出热平衡式得：12r j,jk-2tki,jk. r i，jtrj - - r ii :c -k 1ki, j i i, jattki,j -匕，jjtwjlr -rj：j i，jlt i，j i； +&1中 + j i jt i，j i -甑温rv j 2 jvj 2 j对等式两边同除以j *-r并简化，可以得出与上式完全一样相同的结果。4-7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷却，底面可以认为是绝热的。为用数值法确定冷却过程中柱体温度的变化，取中心角为1rad的区域来研究(如本题附图

12、所示)。已知柱体表面发射率，自然对流表面传热系数，环境温度，金属的热扩散率，试列出图中节点(1, 1), (m,1)(m,n)及(m,n)的离散方程式。在 r及z方向上网格是各自 均分的。解：应用热平衡法来建立四个节点点离散方程。节点(1, 1)：a工k k x k 二 / 2 2 jxkx k/ a儿/ i 21，2 t 1，1包 | +/ t 2，1 t 1，1 i 包，包=pc 生z2l2)j ar 122)18tk-tk ,m _l 1m,1a-7-rrm 22m,1节点(m, n)：tk.,tkm 1, n m, n1.r-tkn 1 m nt中孩t贲htnv 5习裔4电附田4-8、

13、一个二维物体的竖直表面收液体自然对流冷却，为考虑局部表面传热系数的影响，表 1.25面传热系数采用 h = c（t -d 来表示。试列出附图所示的稳态无内热源物体边界节点 （m,n）的温度方程，并对如何求解这一方程提出你的看法。设网格均分。解：利用热平衡法:0.25 h c (tm , n tf jftm, n tf )0.25将h写为h =c（tm，n 一力xtm，口一匕），其中i, n为上一次迭代值，则方程即可线性化。4-9、在附图所示的有内热源的二维导热区域中，一个界面绝热，一个界面等温（包括节点4）,其余两个界面与温度为tf的流体对流换热, 试列出节点1, 2, 5, 6,h均匀，内热

14、源强度为10的离散方程式。解 ：节t5 -11x t2 - t1a i 十 /u占八、y j1 1 x y - yh t1-tf =0:x 242节点2:t6 -t21 x x y =0y2；节点5:节点6:t1 _t5yt2 -t6 yt5 -t9lx1. y x y ，.yh t5-tf = 02.汉 ：y - t-5 x 左士 ；：y：x ：y =0xayax；节点9:”一节点10:x2x 2y1.一 x - =x y - xh hi0-tf =02o节点2:2t6 +3+t34t2 + ay2:节点2t6 t1t9 2 i 凹 tf -2 2节点5:.一=6:17 t10t5t76 w

15、 1当ax = ay以上诸式可简化为：节点1:15 t2四 tf -2 2t5 . t10 . 2节点9:2 1.h . ：y ,1.2t9y i 一2hay 2t6 +t9 +t1 +2 - 节点10:l九j一维稳态导热计算tf 212 +50wilk。4-10、一等截面直肋，高 h,厚6 ,肋根温度为t0 导热系数为 九。将它均分成4个节点（见附图）， 侧面）的两种情况列出节点 2,流体温度为tf ,表面传热系数为 h,肋片 并对肋端为绝热及为对流边界条件（ h同 ,3,4 的离散方程式。设o可期4-10的至节点2:：_x-2h x t2 -tf=0节点3:-2h ：x t3 -tf=0节

16、点其中4:肋端绝热 xt3 -14 肋端对流 xhd =- h-x 114 - tf - h： 114-1 f - 0o肋端绝热3 o将已知条件代入可得下列两方程组: t3 -2.045t2 100.9 = 0h=45cm,5=10mm, h=50w/(m2.k)，九=50w/(m.k), to =100c, tf =20c,计算节点2, 3, 4的温度（对于肋端的两种边界条件）解：采用热平衡法可列出节点 2、3、4的离散方程为:t2 -2.045t3 t4 0.9 = 0t3 -1.0225t4 0.45 =0肋端对流t3 -2.045t2 100.9 = 0t2 -2.045t3 t4 0

17、.9 = 0t3 -1.0375t4 0.8=0由此解得：肋端绝热t2 =92.2 ct3 = 87.7c, t4=86.2c；肋端对流.眈，86.20 t4 =83.80c。肋端对流换热的条件使肋端温度更接近于流体温度。4-11、复合材料在航空航天及化工等工业中日益得到广泛的应用。附图所示为双层圆筒壁,假设层间接触紧密，无接触热阻存在。已知r1 =12.5mm,r2 =16mm,r3 =18mmj” = 40w/(mk)九2 = 12cw/(m.k),tf1 = 150 c,2h1 =1000w/(m .k),tf2 =60c, h2 =380w/(m .k)o试用数值方法确定稳态时双层圆筒

18、壁截面上的温度分布。解：采用计算机求解，答案从略。采用热平衡法对两层管子的各离散区域写出能量方程，进行求解；如果采用taylor展开法列出方程，则需对两层管子单独进行， 并引入界面上温度连续及热流密度连续的条件，数值计算也需分两区进行，界面耦合。截面的温度分布定性地示于上图中。4-12、有一水平放置的等截面直杆，根部温度to =100c,其表面上有自然对流散热，h=c&_tf /d 1 ,其中，cml2wkjc)为杆直径，m。杆高 h=10cm ,直径 d=1cm,九=50w/(m.k) , too - 25 c。不计辐射换热。试用数值方法确定长杆的散热量(需 得出与网格无关的解。杆的两端可认

19、为是绝热的。解：数值求解过程略，q=2.234w。4-13在上题中考虑长杆与周围环境的辐射换热，其表面发射率为0.8,环境可作为温度为tg的大空间，试重新计算其导热量。解：数值求解过程略，q=3.320w。口为:4-14、有如附图所示的一抛物线肋片，表面形线方程y (x )= + b - e (1 - x/ h 2 ip_2卜肋根温度t0及内热源 中恒定，流体表面传热系数 h,t -tf流体温度tf为常数。定义:t0 -tf=x/ h试：(1)建立无量纲温度 。的控制方程；:h2ev b =0.01,0.05,一纲参数t0 -tfhh(2)在无量hh= 0.1,0.01九下对上述控制方程进行数

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

21、） , a =1.33310 m /s。（提不：节点 4 的离散方程可按端面的对流散热与从节点3到节点4的导热相平衡这一条件列出）。解：三个节点的离散方程为:节点2:tk1 -tk2x/2 、4节点3:tk -tk t 3 t 2x，42、 nd- ld x h tf -tk2 二：cxrk 1 k2 - t 2at节点tk4 -tk3x/24:tk3-tk4x/2nd2tk -tk t 2 t 3xnd2二 d x h tf -tk3 = :?c二 d2gd2 fnd2 4八4h(tk4-tf以上三式可化简为:(x )t +2l，aa +2324ht一 %d jtf4h:cdf3a x24h

22、：cdtk24h k t 3cd2 xh tk4 = 2 tk3xhtf, 3a.：.4h.：.1 2- 0： _1/稳定性要求 axpcd,即3a4h2心x pcd八pc = 一 a4355=32.258 1051.333x10-，代入得:3x1.333x104x10025 0.0220.01父 32.258父 1050.099975 0.0124= 8.89877s如取此值为计算步长，则：a.：.1.333 10 为 8.898770.0224h;=0.2966-:cd4 100 8.89877一 _ 一 5 一一32.258 100.01=0.1103于是以上三式化成为:2 0.2966

23、t1 0.2966t3k 0.1103tf =tk 12kkk 10.2966t20.2966 2t40.1103tf=t 30.9773t3k 0.0227tf =tk4(al =8.89877s时间点12340200252525 t200128.8125252a t200128.8155.8055.093a t200137.9573.6472.544a弋200143.0486.7085.30, 3a； 4h ：-1 0在上述计算中，由于之值正好使取pcd7为上值之半，则因而对节点 2出现了在 丁及24工时刻温度相等这一情况。如取a -.r =0.1483x,4h 二 =0.05511:cd

24、3a； 4h ： .八2- - = 0.5xpcd,于是有:kkk 12 0.14831 0.1483t30显 0.0551tf =t 2kkkk 10.148320.1483 2t40.5t30.0551tf=t 30.9773 k3 0.0227tf =tk4对于相邻四个时层的计算结果如下表所示：3 =4.4485s )时间点1234020025252520076.9125252a t200102.8632.7032.533a t200116.9842.6342.234a t200125.5152.5751.944-16、一厚为2.54cm的钢板，初始温度为 650 c ,后置于水中淬火，

25、其表面温度突然下 降为93.5 c并保持不变。试用数值方法计算中心温度下降到450c所需的时间。已知a =1.16m10*m2/s。建议将平板8等分，取9个节点，并把数值计算的结果与按海斯勒 计算的结果作比较。解：数值求解结果示于下图中。随着时间步长的缩小，计算结果逐渐趋向于一个恒定值，当7二0.00001s时，得所需时间为 3.92s。如图所示，横轴表示时间步长从1秒，0.1秒，0.01秒，0.001秒，0.0001秒，0.00001秒的变化；纵轴表示所需的冷却时间(用对数坐标表示)。4-17、一火箭燃烧器，壳体内径为400mm，厚10mm,壳体内壁上涂了一层厚为2mm的包裹层。火箭发动时，

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

27、20min,及30min时汽包内壁截面中的温度分布及截面中的最大温差。启动前， 汽包处于100c的均匀温度。汽包可视为一无限长的圆柱体，外表面绝热，内表面与水之间 的对流换热十分强烈。汽包的内径r1 = 0.9m,外半径r2 = 1.01m,热扩散率 _62a =9.98父10 m /s。解：数值方法解得部分结果如下表所示。汽包壁中的最大温差，k启动后时间，min温升速率，k/min13107.13621.41209.46328.393010.1930.574-19、有一砖墙厚为 6=0.3m,九=0.85w/(m.k),牝=1.05父 106 j/(mk)室内温度为 2t1 =20c, h=

28、6w /(m .k)。起初该墙处于稳定状态，且内表面温度为15co后寒潮入侵，室外温度下降为tf2 =tc,外墙表面传热系数 h2 =35w/(m2.k)。如果认为内墙温度下降0.1 c是可感到外界温度起变化的一个定量判据，问寒潮入侵后多少时间内墙才感知到？解：采用数值解法得 t=7900so4-20、一冷柜，起初处于均匀的温度(20c)。后开启压缩机，冷冻室及冷柜门的内表面温度以均匀速度18 c/h下降。柜门尺寸为1.2m m 1.2m。保温材料厚8cm,九=0.02w/(m.k)。 冰箱外表面包裹层很薄，热阻可忽略而不计。柜门外受空气自然对流及与环境之间辐射的加 热。自然对流可按下式计算：

29、h =1.55 ：t/h 1/4 w/(m2.k)其中h为门高。表面发射率 =0.8。通过柜门的导热可看作为一维问题处理。 机起动后2h内的冷量损失。4 ,3 .解：取保温材料的pc=1x10 j/(m k ),用数值计算方法得冷量损失为试计算压缩5.97x104j 。4-21、3砖砌墙壁,厚度为 240mm, k = 0.81w/(m.k), p = 1800kg/m，c =0.88j /(kg.k )。设冬天室外温度为 24h内变化如下表所示。室内空气温度 ti =15c且保持不变；外墙表面2 .传热系数为10w/(m .k),内墙为及墙壁中心处温度随时间的变化。取2 .6w/(m .k)

30、 o试用数值方法确定一天之内外墙，内墙 t =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:0023:00温度/0c-1.02.41.81.81.60.5-1.6-2.8-3.5-4.3-4.8-5.3解：采用数值解法得出的结果如下表所示。时刻/h环

31、境温度/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墙壁中 心温度3.653.323.152.922.872.812.752.592.31/0c内墙温度/0c8.998.828.738.618.588.558.528.438.28时刻 /h91011121314151617环境温度/0c-7-4.9-2.3-12.41.81.81.60.5外墙温度/0c-3.58-3.07-1.340.781.874.634.154.143.97墙壁中心温度2.362.703

32、.875.326.057.957.627.627.51/0c内墙温度/0c8.318.499.119.8710.2611.2611.1011.1011.10时刻 /h181920212223环境温0-1.6-2.8-3.5-4.3-4.8-5.3度/0c外墙温度/0c3.061.340.36-0.22-0.87-1.29墙壁中心温度6.095.73/0c5.054.664.213.93内墙温度/0c10.7110.109.739.539.309.14多维稳态导热问题4-22、如附图所示，一矩形截面的空心电流母线的内外表面分别与温度为tf1,tf2的流体发生对流换热，表面传热系数分别为h1,h2

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

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

35、况作一估算，把附图a所示的截片形状简化成为附图b所示的情形。已知to =1700k,h。=1000w/(m2.k), = 400k,hi = 250 w/(m2.k)。试计算：(1)截 面中最高温度及其位置；(2)单位长度通道上的热量。解：根据对称性选择1/4区域为计算区域，采用60 m 70取壁面九=15w/(m k )时得单位长度的传热量为98： 等温线分布如图所示。截面中最高温度发生在左上角，该处 温度为1419.90c。综合分析与分析、论述题4-25、工业炉的炉墙以往常用红砖和耐火科组成。由于该两种材料的导热系数较大，散热损失较严 重，为了节省能量，近年来国内广泛采用在耐火科 上贴一层

36、硅酸纤维毡，如附图所示。今用以下的非 稳态导热简化模型来评价黏贴硅酸纤维毡的收益：设炉墙原来处于与环境平衡的状态，0 = 0s时内壁表面突然上升到 550 c并保持不变。这一非稳态导 热过程一直进行到炉墙外表面的对流，辐射热损失 与通过墙壁的导热量相等为止。在炉墙升温过程中外表面的总表面传热系数由两部分组成， 即自然对流引起的部分hc，w/m2.k =1.12 %- tf , 0c及辐射部分hr =4；二 0tm2,tm = tw tf /2其中：tw,tw为外表面温度，tf,tf为内表面温度，2=240mml 3 =240mml2=40mm。为简化计算，设三种材料的导热系数分别为储=1.6

37、w/(m.k) , % = 0.8 w/(m.k),1=0.04w/(m.k) o试计算每平方炉墙每平方面积上由于粘贴了硅酸纤维毡而在炉子升温 过程中节省的能量。解：采用数值计算方法，详细过程从略。4-26、空气在附图所示的一长方形截面的送风管道中作充分发展的层流流动，其 z方向的动量方程简化为2 a2二2 、.虫=0n 2 口 2ucy j dzdp而且u = v = 0。上式可看成是源项为dz的一常物性导热方程。试用数值方法求解这一dp方程并计算f,re之值。f为阻力系数，re为特征长度为当量直径 de。计算时可任取一个 dz 值，并按a/b=0.5及1两种情形计算。解：假设壁温为常数，则

38、不同a/b下换热充分发展时的fre及nu数的分析解为：a/bnufre12.98570.53.39624-27、一家用烤箱处于稳定运行状态，箱内空气平均温度ti =155c,气体与内壁间的表面2 ,八传热系数hi =40w/(m .k) o外壁面与20c的周围环境间的表面传热系数ho =10 w/(m .k)。烤箱保温层厚 30mm,九= 0.03w/(m.k),保温层两侧的护板用金属制成且很 薄，分析中可不予考虑，然后，突然将烤箱调节器开大，风 扇加速，内壁温度突然上升到185c,设升温过程中烤箱外壁面与环境间的表面传热系数可用h0 =ctw tf i”计算，环境温度tf仍保持为20 c,

39、tw为烤箱外壁面温度，c之值与 运行时一样。试确定烤箱内壁温度跃升后到达新的稳定状态 所需时间。解：需采用数值方法求解，过程从略。小论文题目4-28、一厚为2.54cm的钢管，初始温度为16 c。其后，温度为 572c的液态金属突然流过2管内，并经历了 10s。液态金属与内壁面间的表面传热系数h=2.84w/(m .k) o钢管可以按平壁处理，其外表面的散热由对流及辐射两条路径，并分别可按 (hwl(m.k 尸.2父&% a及 =4m0tm计算，tm =tf十丁卬v2 ,周围环境温度tf = 20co试用有限差分法确定在液态金属开始流入后的18s时截面上的温度分布。已知钢管的九=41w/(m.

40、k) , p = 7530kg/m3 , c=536j/(kg.k)。解：在钢管壁厚方向上取 27个点，以内壁为坐标原点，沿着壁厚方向为 x正方向，数值计 算结果如下。位置/cm0.10.20.30.40.50.60.70.80.9温度0c 216.0215.6214.6213.0210.7207.9204.6200.8196.6192.1位置1 /cm1.11.21.31.41.51.61.71.81.9温度/0c187.3182.3177.2172.0166.9161.8157.0152.5148.2144.4位置22.12.2/cm2.32.42.52.54温度/0c141.0138.1

41、135.7133.9132.6132.0131.9用图形表示如下4-29、为对两块平板的对接焊过程（见附图a）进行计算，对其物理过程作以下简化处理：钢板中的温度场仅是 x及时间t的函数；焊枪的热源作用在钢板上时钢板吸收的热流密度 q（x） = qmegr ”e i 1为电弧有效加热半径，qm为最大热流密度；平板上下表面的散热可 用口=卜6-3 1十算，侧面绝热；平板的物性为常数，熔池液态金属的物性与固体相同；固 体熔化时吸收的潜热折算成当量的温升值，即如设熔化潜热为l,固体比热容为c,则当固体 达到熔点ts后要继续吸收相当于使温度升高（ l/c）的热量，但在这一吸热过程中该温度不变。这样，附图a所示问题就简化为附图 b所示的一维稳态导热问题。试： （1）列出该问题 的数学描写；（2）计算过程开始后 3.4s内钢板中的温度场，设在开始的 0.1s内有电弧的加422热作用。已知：qm =5。24 灯0 w/m ,h=12.6 w/（m .k） ,九=41.9w/（m.k）,7800kg/m3,c=670j/（kg