燕山大学-传热学试题及答案

2013年《传热学》试题注：请按题目顺序在标准答题纸上作答，答在题签或草稿纸上一律无效。

本试题允许使用计算器。

一、基本概念：（8分）1、温度梯度。

2、黑度。

3、自然对流和强迫对流的概念。

4、辐射力。

二、简答题（12分）1、简述影响对流换热的因素。

2、简述求解导热问题的三类边界条件。

3、论述辐射换热与对流换热、导热有何区别。

4、简述集总参数法求解非稳态导热问题的基本思想。

三、推导题（30分）1、推导无内热源的三维非稳态导热方程τ∂∂=∂∂+∂∂+∂∂tz t y t x t a ][222222。

2、推导无内热源的三维非稳态导热方程的内节点差分方程，并说明求解的稳定性条件。

四、计算题（50分）1、蒸汽直管的外径1d =30mm ，准备外包两层厚度均为15mm 的不同材料的保温层，a 种材料的导热系数a λ=0.04 W/(m ·℃),b 种材料的导热系数b λ=0.1 W/(m ·℃)，若温差一定，从减少热损失的观点看下列两种方案：（1）a 种材料在里层，b 种材料在外层，（2）b 种材料在里层，a 种材料在外层。

试问哪一种方案好，为什么？（通过计算说明）。

（10分）2、直径为30mm ，长为200mm 的钢杆，导热系数λ=45 W/(m ·K)，将其置于恒温的流体中，流体温度t f =20℃，杆的两端保持恒定的250℃（流体与两端面不接触），流体对杆侧面的对流换热系数为 25 W/(m2·K)，问：1）计算离任一端头 100mm 处的温度。

2）计算单位时间钢杆侧面传给流体的热量。

（15分）3、两同心圆球表面，直径分别为 80mm 和 300mm ，温度分别为827℃和227℃，黑度分别为0.8和0.6。

在两圆球表面之间放入一导热热阻可以忽略不计的薄铝球壳，铝球壳的内外表面黑度均为0.28，其内、外直径可视为相同，为120mm 。

问：1）计算没有铝球壳时两球表面间的辐射换热量。

2）计算放入铝球壳后两球表面间辐射换热量。

（15分）4、 冬季室内空气温度 1f t =20℃，室外空气温温度 2f t =-25℃。

室内、外空气对墙壁的对流换热系数分别为 1α=10 W/(m2·K)和 2α= 20 W/(m2·K)，墙壁厚度为δ= 360mm ，导热系数λ=0.5W/(m ·K)，其面积F =15m2。

试计算通过墙壁的热量损失。

（10分）答案一、基本概念1. 温度梯度：当相邻等温面间的距离趋于零时，其法线方向上的温度变化率称为温度梯度，用n nt∂∂或grad t 表示，单位为K/m 。

n nt n t t n ∂∂==→ΔΔlimgrad 0Δ温度梯度是矢量，位于等温面的法线方向上，沿温度增加方向为正。

2.黑度：实际物体的辐射力与同温度下黑体的辐射力之比称为实际物体的黑度（又称发射率）。

bE E=ε 3. 自然对流换热：流体在由于流体内部各部分温度不同而引起密度差异所产生的浮升力作用下发生的流动，此时的对流换热为自然对流换热。

强迫对流换热：流体在泵、风机或其它压差作用下发生的流动，此时的对流换热为强迫对流换热。

4. 辐射力：单位时间内物体单位表面积向半球空间的一切方向发射的包含全部波长的辐射能，称为总辐射力，用符号E 表示，单位是2W/m 。

总辐射力表征物体发射辐射能力的大小，简称为辐射力。

二、简答题：1. 简述影响对流换热的因素。

答：影响对流换热的因素有： 1）. 流体运动发生的原因 2） 流体的流动状态 3）. 流体有无相变 4）. 流体的物理性质5）. 换热面的几何尺寸、形状和布置 2. 简述求解导热问题的三类边界条件。

答：(1)第一类边界条件：给定物体表面上的温度分布随时间的变化。

在特殊情况下，物体表面上各处的温度相同，=w t 常数。

(2)第二类边界条件：给定物体表面上热流密度的分布随时间的变化。

)()(2w w τλf ntq =∂∂-= 在特殊情况下，边界上的热流密度为常数，=∂∂-=w w )(ntq λ常数 (3)第三类边界条件：给定物体边界上的对流换热系数α和周围流体的温度f t 。

在特殊情况下，α和ft 为常数。

第三类边界条件可表示为：w )(Δntt ∂∂-=λα3.简述辐射换热与对流换热、导热的区别。

答：1）辐射能可以在真空中传播，而且实际上在真空中辐射能的传递最有效。

而热传导与热对流只有当存在着气体、液体、固体介质时才能进行。

因此，当两个温度不同的物体被真空隔开时，热传导与热对流均不能发生．而辐射换热照常进行。

这是辐射换热区别于热传导和热对流的根本特点。

2）辐射换热区别于热传导和热对流的另一个特点是，它不仅产生能量的转移，而且还伴随有能量形式的转换。

即发射时从热能转换为辐射能，而被吸收时又从辐射能转换为热能。

4. 简述采用集总参数法求解非稳态导热问题的基本思想。

答：当固体内部的导热热阻远小于其表面的对流换热热阻时，即11<<==λαδαλδi B 时，固体内部的温差就相当地小，以致于可以认为整个物体在同一瞬间处于同一温度下。

这时所要求解的温度仅是时间τ的一元函数，而与坐标无关，好像该物体原来连续分布的质量与热容量汇总到了一点，整个物体在同一时刻只有一个温度值，所以把这种忽略物体内部导热热阻的分析方法称为集总参数法，这是非稳态导热问题中最简单的一种情况，它是零维问题。

三、推导题1. 推导无内热源的三维非稳态导热方程 τ∂∂=∂∂+∂∂+∂∂tz t y t x t a ][222222 。

解：对微元体做能量平衡分析。

在τd 时间内，导入(导出)微元体的净热量等于微元体焓值(又称热力学能或内能)的增加或减少。

即I Q d d 1=式中，1d Q 为导热(导出)微元体的净热量；I d 为微元体焓值的增加(或减少)。

z 1d d d d Q Q Q Q y x ++=根据傅立叶定律可以确定在x 方向导入(或导出)微元体的净热量为：τλτλτd d )(d )d d d (d )d d d (d )d (d x V xtx x z y x t x xx z y q x x Q x x Q Q Q Q x x x x x ∂∂∂∂=∂∂-∂∂-=∂∂-=∂∂-=∂∂+-=同理可得y ,z 方向导入(或导出)微元体的净热量为：τλd d )(d V y ty Q y ∂∂∂∂=τλd d )(d V ztz Q z ∂∂∂∂=τλλλd d )]()()([d 1V ztz y t y x t x Q ∂∂∂∂+∂∂∂∂+∂∂∂∂=设微元体的定压比热为p c J/(kg ·K)，密度为ρ kg/m 3，则τd 时间内微元体的焓值变化为：ττρττρd d )(d )d (d V t c V t c I p p ∂∂=∂∂=则：)()()()(t c z t z y t y x t x p ρτλλλ∂∂=∂∂∂∂+∂∂∂∂+∂∂∂∂ 当λ、p c 、ρ为常数时，上式变为：τ∂∂=∂∂+∂∂+∂∂tz t y t x t a ][222222 式中，ρλp c a =称为导温系数或热扩散系数，单位为m 2/s 。

物质不同。

2.解：无内热源的三维非稳态导热方程为： τ∂∂=∂∂+∂∂+∂∂tz t y t x t a ][222222 （1）将求解区域离散成有限个网格，在时间上，从过程开始起计算，以τ∆划分时间间隔，以k 表示其序号。

因此，物体内某处某时刻的温度，即节点),,(m j i 在k τ∆时刻的温度可表示为km j i t ,,。

于是，节点在),,(m j i Q xQ x +d xQ y Q y +d yQ z Q z +d zz xy 图2-4 导热微元体在k τ∆时刻的温度对坐标的二阶导数的中心差分公式为)2()(1)(,,,,1,,12),,,(22k m j i k m j i k m j i k m j i t t t x x t -+=∂∂-+∆ (a) )2()(1)(,,,1,,1,2),,,(22k m j i k m j i k m j i k m j i t t t y y t -+=∂∂-+∆ (b) )2()(1)(,,1,,1,,2),,,(22k m j i k m j i k m j i k m j i t t t y z t -+=∂∂-+∆ (c) 对时间τ的一阶偏导数的向前差分可表示为)(1)(,,1,,),,,(k m j i k m j i k m j i t t t -=∂∂+ττ∆ (d) 若取l z y x ===∆∆∆，再将式(a)、(b)、(c)、(d)代入式（1）得：k m j i km j i k m j i k m j i k m j i k m j i t la t t t t l a t ,,2,1,,1,,,1,,121,,)61()(ττ∆∆-++++=-+-++ (2) 以准则20la F τ∆=代入上式得 k m j i k m j i k m j i k m j i k m j i k m j i t F t t t t F t ,,0,1,,1,,,1,,101,,)61()(-++++=-+-++ (3)上式就是无热源三维非稳态导热问题的内节点差分方程。

式（2）及式（3）右端第二项可能出现负值而导致温度计算值波动，出现不稳定。

为避免这种异常情况出现，在选择x ∆、y ∆、z ∆和τ∆时，必须满足下面的稳定条件：20la F τ∆=≤61(4) 四、计算题1．解：如右图：151=r mm ，302=r mm ，453=r mm ，通过单位长度圆简壁的热流量为：)ln(21)ln(21232121r r r r tlQ q l πλπλ+==∆a 材料在内时，热阻为：4.3)3045ln(1.021)1530ln(04.021)ln(21)ln(21)ln(21)ln(2123122321211=⨯+⨯=+=+=πππλπλπλπλr r r r r r r r R b ab 材料在内时，热阻为：7.2)3045ln(.04021)1530ln(.1021)ln(21)ln(21)ln(21)ln(2123122321212=⨯+⨯=+=+=πππλπλπλπλr rr r r r r r R a b可知，a 材料在内时热阻大，散热少，在第（1）方案好。

2.解：由于对称性，按一半进行计算，已知：d =30mm=0.03m, l =100mm=0.1m, λ=45, t f =20℃, t 0=250℃，α=25，则：09245.0==d U πm, 0007069.0015.022=⨯==ππr f ,524.8==fUm λα,8524.0=ml 由断面绝热的公式：)(ch )]([ch 0ml x l m -=θθ )th(0ml mf Q θλ=其中：12)0(ch )]([ch 00=+==-e e x l m , 3858.12)ch(8524.08524.0=+=-e e ml 692321.0771659.2918878.1)th(8524.08524.08524.08524.0==+-=--e e e e ml则：)(ch 100ml t t t t f f l l =--=θθ 186203858.120250)(ch 0=+-=+-=f f l t ml t t t ℃换热量为：64.8692321.0)20250(0007069.08524.0452)th(20=⨯-⨯⨯⨯⨯==ml mf Q θλW3.已知： r 1=40mm=0.04m, t 1=827℃,8.01=ε r 3=150mm=0.15m, t 3=227℃,6.01=εr 2=60mm=0.06m, 82.01=ε则：02.04211==r F πm 2，2827.04233==r F π m 2，0452.04222==r F π m 2，8301411001067.51065.5484181=⨯⨯=⨯⨯=--T E b 35435001067.51065.5484383=⨯⨯=⨯⨯=--T E b热阻分别为：5.1202.08.08.011111=⨯-=-F εε，358.22827.06.06.011333=⨯-=-F εε，89.560452.028.028.011222=⨯-=-F εε 5002.011111112113====F F F ϕϕ，12.220452.01112223===F F ϕ 没有辐射屏时换热量:1225358.2505.12354383014111333113111313,1=++-=-++--=F F F E E Q b b εεϕεεW有辐射屏是换热：396358.212.2289.5689.56505.12354383014111111333223222222112111313,1=+++++-=-++-+-++--=F F F F F F E E Q b b εεϕεεεεϕεεW4.解：由传热公式：776152015.036.01012520112121=⨯+++=++-=F t t Q f f αλδαW。