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

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

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

3、。解：ntan n Bi ,不同Bi下前六个根如下表所示: Bi科1科2科3厂科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、试用数值计算证实,对方程组X1 2x2 2x31x1 x2x332x12x2 x3 5用高斯-赛德尔迭代法求解，其结果是发散的，并分析其原因。解：将上式写成下列迭代形式x11/2 5 2x2 x3x21/2 1 2x3 x1x33

4、x1x2假设x2，x3初值为0,迭代结果如下：迭代次数01234x10x20-1.x30-0.显然，方程迭代过程发散因为迭代公式的选择应使每一个迭代变量的系数总大于或等于式中其他变量的系数绝对值代数和。4-3、试对附图所示的常物性，无内热源的二维稳态导热问题用高斯-赛德尔迭代法计算t1，t2，t3，t4 之值。解：温度关系式为：ti1/412t34030t21/4tit42030t31/4tit43015t41/4t2t3105开始时假设取t1t220 C； t30 t40 15 C /得迭代值汇总于表迭代次数02020/1515123百度文库28345623.28.23.28.23.23.2

5、2.22.22.15.15.15.15.其中第五次与第六次相对偏差已小于10 4迭代终止。4-4、试对附图所示的等截面直肋的稳态导热问题用数值方法求解节点 23 的温度t0 850C，tf 250C，h 30W/(m2.K).肋高 H=4cm,纵剖面面积2AL 4cm，导热系数20W/(m.K) o解：对于2点可以列出:t1 t2t3t4节点2:2h x(ti t2)0；节点3:由此得:t2 t3x2h(tf3)t3)0t1t2t3t22h2(tt2t3h(tft3)x2 h-2(tf t3)0t2t1t302h xH22hxH2ht3t230 0.02220 0.01 t2 t3 30/20

6、tf2.12t2 t14.3636t 2t22.53t3h xht于是有:h x20.060.03tf30/20 0.03t2 1.53tf0.12tf2.532.53t1 1.8336tf85 1.8336 254.363659.8 1.53 25 38.752.53t2t2 1.5tf2.535.3636t2t1t2 0.12tf2 0.120.03tf t2 1.53tf2.53，代入得:2.53t1 t2 1.53tf 0.3036tft22.53tf1.8336t f4.3636215.05 45.84 59.79 59.8 C4.363638.8 C离散方程的建立4-5、试将直角坐

7、标中的常物性无内热源的二维稳态导热微分方程化为显式差分格式，并指 出其稳定性条件（x y）。/解：常物性无内热源二维非稳态方程微分方程为2t2t扩散项取中心差分，非稳态项取向前差分:t i1 t iti 2titi ti 2tititn tntn121ntn1 tn121ntn1a2/2xy所以有t；11 2a稳定性条件Fo x Fo y 1/24-6、极坐标中常物性无内热源的非稳态导热微分方程 为2tr试利用本题附图中的符号, 程式。列出节点（12t22ri,j）的差分方解：将控制方程中的各阶导数用相应的差分表示式代 替，可得：ttkt1, j t t, jakkkt t, j 1 2t t

8、, j t j 12rk1 11,rj1引，解将中谱网恃122r jrj也可米用热平衡法。tkc 对于图中打阴影线的控制容积写出热平衡式得:1 tki， j i， jtki1,j tki,jrtki 1,jtki,jrtkitki,jrjrjt" tki,jrjrj对等式两边同除以rjr并简化，可以得出与上式完全一样相同的结果。4-7、一金属短圆柱在炉内受热厚被竖直地移植到空气中冷却，底面可以认为是绝热的。为用数值法确定冷却过程中柱体温度的变化，取中心角为1rad的区域来研究（如本题附图 所示）。已知柱体表面发射率，自然对流表面传热系数，环境温度，金属的热扩散率，试列出图中节点（1,

9、 1）, （M,1）（M,n）及（M,N）的离散方程式。在 r及z方向上网格是各自 均分的。解：应用热平衡法来建立四个节点点离散方程。节点（1, 1）：k k2k kt 1,2 t 1,1 1 rt 2,1t 1,1 r Zr- - 一 一 c 一z 2 2r 228tk 一tk ,m 1,1m,11r节点(mtktkm 1, nm' nrrm m 22k1m>n1 tz-m,m 13 rm-T解m,1rm 13rmh tm tm表面传热系数采用并对如何求解这一方m, nt" ;t1 1,14-8、一个二维物体的竖直表面收液体自然对流冷却，为考虑局部表面传热系数的影响,

10、，,x 1.25/hc（t t1）来表示。试列出附图所示的稳态无内热源物体边界节点（M,n）的温度方程, 程提出你的看法。设网格均分。解：利用热平衡法:h C tM, n将h写为h一次迭代值,tf tM0.25 n tftftM , n则方程即可线性化。0.25tf ,其中tM, n为上4-9、在附图所示的有内热源的二维导热区域中，一个界面绝热，一个界面等温（包括节点4）,其余两个界面与温度为tf的流体对流换热,试列出节点1t5t15, 6,t2t1节点节点节点节点节点h均匀，内热源强度为10的离散方程式。占八、yht12:5:6:9：10:t1t2y2y2t6t2t1t5y2t9t5t6t5

11、t2t6t7t10t5t5t9t10 t9y2t10t11t10y2t6t10当x y以上诸式可简化为:t2h y-tf 2t1节点1:习就4名第图010xy210xy20yh t5t6t5t ftfy2h t9tf 0xh hiotf 0节点2:节点5:节点6:节点9:节点10：2t62t6t7t1t1t10t10t3t9t54t2tft7 4t6h2t6 t9t11 2 一维稳态导热计算4-10、h y2t5y tf0等截面直肋，高 H,厚，肋根温度为t0导热系数为 侧面）H=45cm,o将它均分成4个节点（见附图），的两种情况列出节点 22 .10mm, h 50W/（m .K）=50w

12、/,流体温度为tf,表面传热系数为 并对肋端为绝热及为对流边界条 ,3,4的离散方程式t0100 C, tf 20 C,计算节点h,肋片h同O 设2, 3,4的温度（对于肋端的两种边界条件）O解：米用热平衡法可列出节点2、3、4的离散方程为:节点2:节点3:节点其中t1t2t3 t22hx t2tft2 t3t4 t32hx t34:肋端绝热肋端对流H肋端绝热肋端对流t4tfh x t4tft4 tf 03 o将已知条件代入可得下列两方程组:t3t2t3t3t22.045t22.045t31.0225t42.045t22.045t3100.9t4 0.90.45100.9t4 0.9o00t3

13、 1.0375t4 0.8 0由此解得：肋端绝热t2 92-20C, t3 87.70C ,肋端对流t291.50Ct3 86.20Ct4t486.20C .83.80C肋端对流换热的条件使肋端温度更接近于流体温度。4-11、复合材料在航空航天及化工等工业中日益得到广泛的应用。附图所示为双层圆筒壁， 假设层间接触紧密，无接触热阻存在。已知 r1 12.5mm, r2 16mm,r3 18mm, 1 40 .2 120W /(m.K ),tf1 150 r/W/ ,。 ?22 一一.h1 1000W/(m .K)，tf2 60 C, h2 380W/(m .K)。试用数值方法确定稳态时双层圆 筒

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

15、两端可认为是绝热的。解：数值求解过程略，Q二。4-13在上题中考虑长杆与周围环境的辐射换热，其表面发射率为，环境可作为温度为t的大空间，试重新计算其导热量。解：数值求解过程略，Q二。试：(1)建立无量纲温度的控制方程；(2)在无量纲参数 t0 tf算。确定无量纲温度0.01, 0.05, 0.1,hH0.01下对上述控制方程进行数量计的分布。d2 /d20.01 2 / 5解：无量纲温度方程为：下图中，无量纲温度从肋根的1变化到肋端的。0。数值计算结果示于一维非稳态导热计算4-15、一直径为1cm,长4cm的钢制圆柱形肋片，初始温度 为25C,其后，肋基温度突然升高到 200C,同时温度为25

16、C的气流横向掠过该肋片，肋端及两侧的表面传热系数均为2100W/（m .K）。试将该肋片等分成两段（见附图）限差分法显式格式计算从开始加热时刻起相邻4温度分布（以稳定性条件所允许的时间间隔计算依据）,并用有 个时刻上的O已知52= 43W/, a 1.333 10 5m2程可按端面的对流散热与从节点 条件列出）。/ s。（提示：节点4的离散方3到节点4的导热相平衡这一解：三个节点的离散方程为: 节点2:tk1 tk2 dx/24节点3:tk tkt 4 Jx/2节点4:tk tkt 3 t 4x/2dtk2d2以上三式可化简为:x h tf tk2tk2x h tftk3dcx4k 1dxtk

17、2 J t1x4hcd3a1 -2 x4hcdtk2tkJ t2xxh tk42 tJ t4xk3 xhtf4hcdtf3a112x4hcdtk31稳定性要求3a2 x4hcd1/3a-2x4hcd o43_ 51.333 1032.258ml代入得:1/3 1.333 104 100120.0220.01 32.2581058.89877s0.099975 0.0124如取此值为计算步长，则：_ 5a 1.333 108.8987722x0.020.29664h4 100 8.89877于是以上三式化成为:cd32.258一 5 一一100.010.1103Xo2 0.2966tl 0.29

18、66t3k 0.1103tftk120.2966t0.9773tkk2 0.2966 2t40.1103tfkk30.0227tf t 4tk13时间点1234020025252520025252A2003A2004A2008.89877s之值正好使3a 4h12 一x (在上述计算中，由于因而对节点 2出现了在 0 cd时刻温度相等这一情况。如取4h0.1483cd0.05513a2 x4h 0.5cd ,于是有:kk2 0.14831 0.1483t30品0.0551tf0.14830.97732k 0.1483 2t4k 0.5t3k k k3 0.0227tft 40.0551tftk

19、12tk13为上值之半，则4.4485s时间点12/34020025/252520025252A2003A200对于相邻四个时层的计算结果如下表所示:4A2004-16、一厚为的钢板，初始温度为650 C,后置于水中淬火，其表面温度突然下降为C并保持不变。试用数值方法计算中心温度下降到450 c所需的时间。已知52.a 1.16 105m2/s。建议将平板. 8等分，取9个节点，并把数值计算的结果与按海斯勒 计算的结果作比较。解：数值求解结果示于下图中。/随着时间步长的缩小，计算结果逐渐趋向于一个恒定值，当=时，得所需时间为。J、如图所示，横轴表示时间步长从1秒，秒，秒，秒，秒，秒的变化；纵轴

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

21、算汽包内壁的温度场。试用数值方法计算：当汽包内的饱和水温度上升的速率为1C/min,3C/min时，启动后10min,20min,及30min时汽包内壁截面中的温度分布及截面中的最大温差。启动前，汽包处于100c的均匀温度。汽包可视为一无限长的圆柱体，外表面绝热，内表面与水之间/的对流换热十分强烈。汽包的内径R1 0.9m，外半径R2 1.01m，热扩散率 a 9.98 10 6m2/so/解：数值方法解得部分结果如下表所示。'汽包壁中的最大温差，K/启动后时间，min 温升速率，K/min131020/304-19、有一砖墙厚为0.3m ,=, c 1.05 10 J /(m .K)

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

23、55 t/H 1/4W/(m2.K)其中H为门高。表面发射率0.8。通过柜门的导热可看作为一维问题处理。试计算压缩机起动后2h内的冷量损失。43-解：取保温材料的 c 1 10 J / m K ,用数值计算方法得冷量损失为“5.97 10 J。34-21、一砖砌墙壁，厚度为 240mm,=,1800kg/m ,C 0.88J/kg.K、。设冬天室外温度为 24h内变化如下表所示。室内空气温度ti 15c且保持不变；外墙表面传热系数为10W/(m2.K) ,内墙为6W/(m2.K)。试用数值方法确定一天之内外墙，内墙 及墙壁中心处温度随时间的变化。取1h o设上述温度工况以 24h为周期进行变化

24、。时刻/h0:00 1 : 00 2:00 3:004: 005:00 6:00 7:008:00 9:0010:0011 :00温度0c时刻12:13:14:15:16:17:/h00000000000018:19:20:21 :22:23:000000000000温度/0C解：采用数值解法得出的结果如下表所示。时刻/h环境温度/0C外墙温外墙温度/0C墙壁中 心温度/0C内墙温度/0C时刻/h 91011121314151617环境温度/0C-7-1度/0C墙壁中心温度/0C内墙温度/0C时刻/h 181920212223环境温度/0C外墙温度/0C墙壁中/0C心温度内墙温度/0C 多维稳

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

26、) 内外壁分别维持在 10c及30 C(2) 内外壁与流体发生对流换热，且有tf1 10C, hi 20W/(m2.K) , tf2 30 C,2h24W/(m .K)。解：此题应采用计算机求解。如有墙角导热的热点模拟实验设备，则计算参数(如h, t及网格等)可以取得与实验设备的参数相一致，以把计算结果与实测值作比较。根据对称性，取1/4区域为计算区域。数值计算解出，对于给定壁温的情形，每米长通道的冷损失为，对于第三类边界条件为(取壁面导热系数0-53W/ m K )。内外表面为给定壁温时等温线分布如下图所示。第三类边界条件的结果定性上类似。4-24、为了提高现代燃气透平的进口燃气温度以提高热

27、效率,冷却通道以使叶片金属材料的温度不超过允许值,在燃气透平的叶片内部开设有为对叶片中的温度分布情况作一估算,附图a所示的截片形状简化成为附图b所示的情形。已T0 1700K,h0 1000W/(m2.K) , Ti 400Khi 250 W/(m2.K)。试计算：(1)面中最高温度及其位置；(2)单位长度通道上的热量。解：根据对称性选择1/4区域为计算区域，采用60 70网格, 取壁面15W/ m K时得单位长度的传热量为，等温线 分布如图所示。截面中最高温度发生在左上角，该处温度为 0c 。综合分析与分析、论述题4-25、工业炉的炉墙以往常用红砖和耐火科组成。 由于该两种材料的导热系数较大

28、，散热损失较严 重，为了节省能量，近年来国内广泛采用在耐火科 上贴一层硅酸纤维毡，如附图所示。今用以下的非 稳态导热简化模型来评价黏贴硅酸纤维毡的收益：设炉墙原来处于与环境平衡的状态,0s时内壁表面突然上升到 550 c并保持不变。这一非稳态导 热过程一直进行到炉墙外表面的对流，辐射热损失与通过墙壁的导热量相等为止。在炉墙升温过程中外表面的总表面传热系数由两部分组成, 即自然对流引起的部分1/3hc W/m2.K1.12 tw 0 c tf 0c及辐射部分i耳我鼾"藤、2 .hr 4 0Tm,Tm Tw Tf / 2其中：tw，Tw为外表面温度，tf,为内表面温度，1 240mm,

29、2 240mm,3 40mm。为简化计算，设三种材料的导热系数分别为11.6 W/,20.8W/,30.04 W/o 试计算每平方炉墙每平方面积上由于粘贴了硅酸纤维毡而在炉子升温过程中节省的能量。解：采用数值计算方法，详细过程从略。4-26、空气在附图所示的一长方形截面的送风管道中作充分发展的层流流动，其 z2w2x方向的动量方程简化为2w2ydz 0dP而且u v 0。上式可看成是源项为dz的一常物性导热方程。试用数值方法求解这一dP方程并计算f,Re之值。f为阻力系数，Re为特征长度为当量直径 De。计算时可任取一个 dz值，并按a/b = 及 1两种情形计算。/解：假设壁温为常数，则不同

30、a/b下换热充分发展时的fRe及Nu数的分析解为：a/bNufRe157155 C ,气体与内壁间的表面624-27、一家用烤箱处于稳定运行状态，箱内空气平均温度Li2传热系数hi40W/(m .K) o外壁面与20c的周围环境间的表面传热系数ho 10 W4m.K)。烤箱保温层厚 30mm,0.03W/,保温层两侧的护板用金属制成且很薄，分析中可不予考虑，然后，突然将烤箱调节器开大， 风扇加速， 内壁温度突然上升到185 C,设升温过程中烤箱外壁面与环-.1/4境间的表面传热系数可用 h0 ctw tf 计算，环境温度 匕仍保持为20 c, tw为烤箱外壁面温度，c之值与运行时一 样。试确定

31、烤箱内壁温度跃升后到达新的稳定状态所需时间。解：需采用数值方法求解，过程从略。小论文题目4-28、一厚为的钢管，初始温度为16C。其后，温度为572c的液态金属突然流过管内，并2经历了 10s。液态金属与内壁面间的表面传热系数h=W/(m K)。钢管可以按平壁处理，1/3其外表面的散热由对流及辐射两条路径，并分别可按hWW/m.K1.2 t 0c及hr 4。丁：计算，Tm Tf Tw/2,周围环境温度tf =20 C。试用有限差分法确定在 3液态金属开始流入后的18s时截面上的温度分布。已知钢管的41W/,7530kg/m ,c=536J/。解：在钢管壁厚方向上取 27个点，以内壁为坐标原点，

32、沿着壁厚方向为x正方向，数值计算结果如下。位置 0/cm/温度/0C位置1/cm温度/0C位置/cm温度/0C用图形表示如下1.口 -口工口1叼amo口 J92事amQ.QI9 也0 anM向”的暂寓X/cm4-29、为对两块平板的对接焊过程（见附图a）进行计算，对其物理过程作以下简化处理：钢板中的温度场仅是 x及时间 的函数；焊枪的热源作用在钢板上时钢板吸收的热流密度3r2/re2qxqme, 上为电弧有效加热半径，qm为最大热流密度；平板上下表面的散热可用q ht tf计算，侧面绝热；平板的物性为常数，熔池液态金属的物性与固体相同；固体熔化时吸收的潜热折算成当量的温升值，即如设熔化潜热为L

33、,固体比热容为C,则当固体达到熔点ts后要继续吸收相当于使温度升高（ L/C）的热量，但在这一吸热过程中该温度不变。这样，附图a所示问题就简化为附图 b所示的一维稳态导热问题。试： （1）列出该问题 的数学描写；（2）计算过程开始后内钢板中的温度场，设在开始的内有电弧的加热作用。已知 ：qm 5024 104W/m2h=W/（m2.K）3=,7800kg/m , c 670J / kg.K ,L=255kJ/kg, ts 1485 C , H=12cm,re 0.71cmobl解：取初始温度与环境温度均为200 c 。该问题的数学描写为:2tc 0<x<H,t>0t tjH ； ;q x 2h t tf />0;>0 o为了更好分辨热源附近的温度场宜采用非均分网格。 分布如下图所示。计算得出开始加热后的内钢板中的温度4-30、在壁厚为7