传热学讲义—第二章
第二章 稳态导热
本章重点：具备利用导热微分方程式建立不同边界条件下稳态导热问题的数学模型的能力
第一节 通过平壁的导热
1-1 第一类边界条件 研究的问题：
（1）几何条件：设有一单层平壁，厚度为δ，其宽度、高度远大于其厚度（宽度、高度是厚度的10倍以上）。这时可认为沿高度与宽度两个方向的温度变化率很小，温度只沿厚度方向发生变化。（属一维导热问题）
（2）物理条件：无内热源，材料的导热系数λ为常数。 （3） 边界条件：假设平壁两侧表面分别保持均匀稳定的温度
1w t 和2w t ，21w w t t >。（为第一类边界条件，同时说明过程是稳态的）
求：平壁的温度分布及通过平壁的热流密度值。 方法1 导热微分方程：
采用直角坐标系，这是一个常物性、无内热源、一维稳态导热问题（温度只在 x 方向变化）。
导热微分方程式为：022=dx
t
d （2-1）
边界条件为：10w x t t == ， 2w x t t ==δ （2-2）
对式（2-1）连续积分两次，得其通解： 21c x c t += （2-3）
这里1c 、2c 为常数，由边界条件确定 ，解得：??
???=-=
11221w
w w t c t t c δ （2-4）
最后得单层平壁内的温度分布为： x t t t t w w w δ
2
11--
= （2-5）
由于δ 、1w t 、2w t 均为定值。所以温度分布成线性关系，即温度分布曲线的斜率是常数（温度梯度），
const t t dx dt w w =-=δ
1
2 （2-6）
热流密度为：)(21w w t t dx dt q -=-=δ
λ
λ
2/m W （2-7） 若表面积为 A, 在此条件下 , 通过平壁的导热热流量则为 :
t A qA ?==Φδ
λ W （2-8）
考虑导热系数随温度变化的情况：
对于导热系数随温度线形变化，即)1(0bt +=λλ，此时导热微分方程为：0=??
?
??dx dt dx d λ 解这个方程，最后得：
??
?
???++-+??
? ??
+=+)(211212121121
122w w w w w w t t b x t t bt t bt t δ 或 x t
t t t b b t b t w w w w w δ
12211)(2112
2-??????+++??? ??+=??? ??+
说明：壁内温度不再是直线规律，而是按曲线变化。
对上式求导得：???
?????+??? ??-=)1/(222bt dx dt b dx t
d
因为 01>+bt ，02
>??
?
??dx dt
所以 0>b ? 02
2 d ? 曲线是向上凸的； 0dx t d ? 曲线是向上凹的。 通过平壁的导热热流密度为： () ?? ????++-=+-=-=2121211)1(00w w w w t t b t t dx dt bt dx dt q λδλλ 式中，()m w w t t b λλλλ=+=?? ????++22112 1 021 则 )(2 1 w w m t t q -= δ λ 从上式可以看出，如果以平壁的平均温度2 2 1w w m t t t +=来计算导热系数，则平壁的热流密 度仍可用导热系数为常数时的热流密度计算式： δ λ2 1w w m t t q -= 多层平壁（复合壁）的导热问题 多层壁（复合壁）：就是由几层不同材料叠加在一起组成的平壁。 以下讨论三层复合壁的导热问题，如图所示。 假设条件：层与层间接触良好，没有引起附加热阻（亦称为接触热阻）也就是说通过层间分界面时不会发生温度降。 已知各层材料的厚度为： 1δ、2δ、3δ ，导热系数为：1λ、2λ、3λ，且均为常数。多层壁的最外两侧表面分别维持均匀稳定的温度1w t 和4w t ，且41w w t t >。 求：该多层平壁中的温度分布和通过平壁的导热量。 设两个接触面的温度分别为2w t 和3w t 。 此问题是无内热源一维稳态导热。整个过程是由三个换热环节串联而成，每个环节的热流密度是相等的。 ∑=-= ++-= 3 1 ,3 ,2,1,4 14 1i i w w w w R t t R R R t t q λ λλλ （三层平壁单位面积的总热阻等于各层热阻之和） 1,12λqR t t w w -= )(2,1,3,143λλλR R q t qR t t w w w +-=+= 因为每层平壁的温度分布都是直线，各层中直线的斜率是不同的，所以多层平壁中的温度分布是一条折线。 对于n 层多层平壁，热流密度：∑=+-= n i i w w R t t q n 1 ,1 1λ 1-2 第三类边界条件 研究的问题： （1）几何条件：设有一厚度为δ的无限大平壁。。 （2）物理条件：无内热源，材料的导热系数λ为常数。 多层平壁的导热 （3）边界条件：给出第三类边界条件，即：在0=x 处，界面外侧流体的温度为1f t ，对流换热表面传热系数为1h ；在δ=x 处，界面外侧流体的温度为2f t ，对流换热表面传热系数为 2h 。 求：平壁的温度分布及通过平壁的热流密度值。 常物性、无内热源、一维稳态导热过程的导热微分方程式仍为：022=dx t d 边界条件：)(010 1==-=-x f x t t h dx dt λ )(22f x x t t h dx dt -=-==δδ λ 解得：)(1121212 1f f f f t t k h h t t q -=++-= λδ 1 111h q t t f w ? -= 2 11)1( 212h q t h q t t f f w ?+=+?-=λδ 求出1w t 、2w t ，就可得出平壁中的温度分布：x t t t t w w w δ 2 11-- = 补充：对于上述的常物性、无内热源、一维稳态导热问题，如果给定第二类边界条件，会出现什么情况？ 第二类边界条件：10C q x == 和 2C q x ==δ 由于是无内热源，稳态导热，所以21C C =，这意味着，上述两个条件是一致的，实际上就是一个条件。 根据这样一个条件，不能求出方程022=dx t d 的通解 21c x c t +=中的两个待定常数1c 和2c 。 问题的解为不定解。 所以，对于一维稳态导热问题，必须具有两个独立的边界条件才能确定出惟一的解。 第二类边界条件下的温度分布曲线： 根据 dx dt q λ-=，得 C dx dt dx dt x x == ==δ ，所以平壁内的温度分布曲线为已知斜率C 的一 簇平行直线。 t w 1 工程上会遇到这样一类平壁：无论沿宽度还是厚度方向，都是由不同材料组合而成——复合平壁。在复合平壁中，由于不同材料的导热系数不同，严格地说复合平壁的温度场是二维或三维的。 如：空斗墙、空斗填充墙、空心板墙、夹心板墙。 复合平壁中，由于不同材料的导热系数不同，严格地说复合平壁的温度场是二维或三维的。 简化处理：当组成复合平壁的各种不同材料的导热系数相差不大时，可近似当作一维导热问题处理。 复合平壁的导热量： ∑?= Φλ R t 式中，t ?——两侧表面棕温差； ∑λR ——总导热热阻。 ∑+++ +++ ++= 3 322 111 11 1 E D A E C A E B A R R R R R R R R R R λλλλλλλλλλ 工程中常用圆管作为换热壁面，如锅筒、传热管、热交换器及其外壳。圆筒受力均匀、强度高、制造方便。 3-1 第一类边界条件 研究的问题： （1）几何条件：单层圆筒壁面，内半径为1r ，外半径为2r ，长度为l ，长度l 远大于壁厚。（忽略轴向热流，热流只沿径向） （2）物理条件：无内热源，圆筒壁材料的导热系数λ为常数。 （3） 边界条件：圆筒壁内、外表面分别维持均匀稳定的温度1 w t 和2w t ，且21w w t t >。（为 第一类边界条件，同时说明过程是稳态的） 求：圆筒壁内的温度分布及通过圆筒壁的导热量。 根据以上条件知，这是一个常物性、无内热源、一维、稳态导热问题。由于温度场是轴对称的，所以采用圆柱坐标系。 导热微分方程为： 0)(=dr dt r dr d 圆筒壁边界条件为： 11 w r r t t == 22 w r r t t == 微分方程的通解为：21ln c r c t += 根据边界条件，得出： 121ln 21r r t t c w w -= 和 11 22ln ln 2 11r r r t t t c w w w -+= 则圆筒壁的温度分布为：112ln ln 211r r r r t t t t w w w --= 或 112 ln ln 211d d d d t t t t w w w --= 由此可见，圆筒壁中的温度分布呈对数曲线，而平壁中的温度分布呈线性分布。 圆筒壁的导热量 在无限大平壁中，热流密度是常数，但在圆筒壁中，不同半径处的热流密度并不相等。（dr dt q λ -=，但dr dt 不等于常数，它是r 的函数） 在稳态情况下，通过长度为l 的圆筒壁的导热量是恒定的，即： dr dt A A q λ-=?=Φ W （A 是圆筒壁的面积，在不同的r 处，有不同的A 值） 在圆筒壁内，取一个半径为r ，厚度为dr 的微圆筒壁来分析，此时，rl A π2=，则： dr dt l r ??-=Φπλ2， 而 r r r t t dr dt w w 1 ln 1 212?-= 解得：1 2ln 22 1r r t t l w w -? =Φπλ （可见，Φ与r 无关，通过整个圆筒壁面的热流量不随半径的变化而变化，在不同的r 处，通过的热流量是相等的。） 将Φ写成热阻形式，则：1 2ln 212 1r r l t t w w πλ-= Φ W 式中， 1 2 ln 21r r l πλ是长度为l 的圆筒壁的导热热阻，W K / 通过每米长圆筒壁的热流量为： 1 2 ln 212 1r r t t l q w w l πλ-= Φ= m W / 单位长度圆筒壁的导热热阻为： 1 2 ln 21r r R l πλ λ= W K m /? 多层圆筒壁的导热 多层圆筒壁：由几层不同材料紧密结合所构成的圆筒壁。 利用串联热阻跌价原理求解。该部分自学。 ∑∑=+=++-= -= n i i i i w w n i i l n w w l d d t t R t t q n 111 ,1 ,ln 21 1,11πλλ [例2-4] 自学。注意：求各层直径时，应是δ2+d 。对于圆管外，用几层材料进行保温时，应将导热系数少的材料设置在内侧。对平壁有这种要求吗？ 3-2 第三类边界条件 研究的问题： （1）几何条件：单层圆筒壁面，内半径为1r ，外半径为2r ，长度为l ，)(12r r l ->>。 （2）物理条件：无内热源，圆筒壁材料的导热系数λ为常数。 （3） 边界条件：已知1r r =一侧的流体的温度为1f t ，对流换热表面传热系数为1h ，2 r r =一侧流体的温度为2f t ，对流换热表面传热系数为2h ，且1f t >2f t 。 求：圆筒壁内的温度分布及通过圆筒壁的导热量。 这是一个常物性、无内热源、一维、稳态、导热问题。由于温度场是轴对称的，所以采用圆柱坐标系。 导热微分方程为： 0)(=dr dt r dr d 圆筒边界条件为： )(221 11 111r r f r r t t r h r dr dt ==-??=?-ππλ )(2222 2 222f r r r r t t r h r dr dt -??=?-==ππλ 在第一类边界条件中，已求出圆筒壁内的温 度变化率： r r r t t dr dt w w 1 ln 1 212?-= 根据傅立叶定律的表达式，任意r 处，单位长度圆筒壁的导热量为： dr dt r q l ? ?-=πλ2 这样，可将边界条件式改写为：)(2111 11w f r r l t t r h q -??==π )(2222 22f w r r l t t r h q -??==π 而圆筒壁的导热量为：1 2ln 212 1r r t t q w w l πλ-= 在稳态导热过程中，2 1 r r l l r r l q q q ====，可见在上述三个方程中，又三个未知数：l q 、1 w t 和2w t ，方程是可解的。 解得：22121121 ln 21212 1r h r r r h t t q f f l ππλπ?+ +?-= W/m 或 2212111 ln 2112 1d h d d d h t t q f f l ππλπ?+ +?-= W/m 也可表示为：)(21f f l l t t k q -?= 式中，l k ——传热系数，表示冷、热流体之间温差为1℃时，单位时间通过单位长度圆筒壁的传热量，W/(m ·K)。 单位长度圆筒壁的传热热阻为：2 212111 ln 2111d h d d d h k R l l ππλπ?+?+?== W K m /? 1 1111d h q t t l f w π?? -= ??? ? ??+??-=??+=121122ln 21 11 122 d d d h q t d h q t t l f l f w πλππ 根据 11 2 ln ln 211r r r r t t t t w w w -- = 可求出圆筒壁中的温度分布。 对多层圆筒壁，热流体通过圆筒壁传给冷流体的热流量为： 1 211111 ln 2112 1+=+?++?-= ∑n n i i i i f f l d h d d d h t t q ππλπ [例2-5] 课后自学。 [思考题] 若平壁与圆筒壁的材料相同，厚度相同，温度条件相同，且平壁的表面积等于圆筒的内表面积，试问：哪一种情况导热量大？ 3-3 临界热绝缘直径 工程上，为减少管道的散热损失，常在管道外侧覆盖热绝缘层或称隔热保温层。 问题：覆盖热绝缘层是否在任何情况下都能减少热损失？保温层是否越厚越好？为什么？ 分析圆管外覆盖有一保温层的情况。 对于冷、热流体之间的传热过程，给定第三类边界条件，则传热过程的热阻为： x x ins l d h d d d d d h R ππλπλπ22121 1 11 ln 21ln 211+++= 下面分析一下l R 随保温层外径x d 的变化情况。 对于一个管道进行分析时，l R 中的前两项热阻的值是确定的，在选定了保温材料后，ins λ也就确定了。这样，l R 的后两项热阻的数值随保温层的x d 而变化。 当↑x d ，???????↓↑增加，外表面的对流换热量，导热量减少x x ins d h d d ππλ221ln 21 但对l R ，随着x d 的增大（保温层加厚），先是逐渐减小，然后又逐渐增大，有一极小值。（相应地，l k 先增大，然后减小，l k 有极大值） 对l q ，随着x d 的增大，先是增大，然后减小，有一极大值。 临界热绝缘直径c d ：对于总热阻l R 为极小值时的保温层外径。 令 012112=??? ? ???-=x ins x x l d h d dd dR λπ 得 2 2h d d ins c x λ== （c d 只取决于ins λ和2h ，c d 不一定大于2d ） 从图中可见： （1）当c d d <2时，如果管道保温后的外径x d 在32~d d 之间，这时 管道的传热量l q 反而比没有保温层时更大，直到3d d x >时，才起到减少热损失的作用。 （2）当c d d >2时，l R 及l q 均是x d 的单调函数，用保温层肯定能减少热损失。 c d 的大小与ins λ和2h 有关，2h 主要取决于管道周围的环境，难以人为控制，但可以通过 选用不同的保温材料来改变c d 的值，使2d d c <，以达到只要使用保温材料就能保证减少热损失的目的。（工程上，一般2d 均会大于c d ，只有当2d 较小时，才需要注意临界绝缘直径的问题。工程上， 尽可能要求2d d c <） [思考题] 解释现象：某厂一条架空敷设的电缆使用时发现绝缘层超温，为降温特剥去一层绝缘材料，结果发现温度更高。 答：电缆外径小于了临界热绝缘直径时，导热热阻随半径增大的变化率小于对流换热随半径减小的变化率，使散热能力随半径增加而增加。剥去一层绝缘材料后，半径减小，散热能力下降，绝缘层温度更高。 第五节 通过肋壁的导热 第三类边界条件下通过平壁的一维稳态导热： A h A A h t t f f 21112 1++-= Φλδ 为了增加传热量，可以采取哪些措施？ （1）增加温差)(21f f t t -，但受工艺条件限制。 （2）减小热阻： 1）金属壁一般很薄（δ 很小）、热导率很大，故：导热热阻一般不大，可忽略。 2）增大h 1、h 2，(但提高h 1、h 2并非任意的)。 3）增大换热面积A 。 强化传热的基本思路：强化传热?降低热阻?降低串联热阻的最大项（主要矛盾） 对流换热热阻 ↑ ↓= hA R h 1 在表面传热系数不变的情况下，要降低对流换热热阻，就必须扩大对流换热面积。其方法之一就是对传热表面进行肋化（加装肋片）。 肋片：指依附于基础表面上的扩展表面 常见肋片的结构：针肋、直肋、环肋等 应用：工程应用十分广泛，如 汽车水箱、空调系统的蒸发器、冷凝器、锅炉的空气预热器、省煤器、散热器等。 肋片导热的作用及特点： ① 作用：增大对流换热面积及辐射散热面，以强化换热 ② 特点：在肋片伸展的方向上有表面的对流换热及辐射散热， 肋片中沿导热热流传递的方向上热流量是不断变化的。即： Φ≠const 。 肋片分析的任务： ① 确定沿肋片高度方向的温度分布； ② 确定肋片的散热量。 一、通过等截面直肋的导热（稳态情况下） 等截面直肋：从平直基面上伸出，而本身又具有不变截面的肋。 设肋片的高度为l ，宽度为L ，厚度为δ。肋片的横截面积为δ?=L A L ，肋片的横截面的周边长度为)(2δ+=L U 。 肋基的温度为0t =const ，金属肋片的导热系数为λ，周围流体的温度为f t ，肋片与流体的对流换热表面传热系数为h 。 求：肋片中的温度分布及通过该肋片的散热量。 采用直角坐标系，原点设在肋基处。 为简化分析，作以下假定： （1）肋片的宽度L 很大 ? 不考虑温度沿该方向的变化?肋片的温度分布是二维温度场，即),(y x f t =，此时传热情况是： 在x 方向，即沿肋片高度方向，热量从肋基以导热方式导入，随后热量继续沿 x 方向传递； 在y 方向，通过对流换热从肋片表面向周围介质散热。 （2）导热系数λ和表面传热系数h 均为常数。 （3）λ大，l <<δ? 肋片沿厚度方向的温度变化很小 ? 可认为任一横截面上的温度分布几乎是均匀 的，截面上各点温度与截面中心的温度一致?温度只是在沿肋片高度方向发生明显变化，温度分布是沿x 方向的一维温度场。 （4）肋片顶端视为绝热。 从能量平衡入手分析该问题： 在距肋基x 处，取一微元dx ，研究该微元的能量平衡。 在x 处导入的热量x Φ，应等于dx x +处导出的热量dx x +Φ和从表面传入流体的热量c Φ，即： c dx x x Φ=Φ-Φ+ 根据傅立叶定律和牛顿冷却公式有： ()f L L L t t dx hU dx dx t d A dx dt A dx dt A -??=?? ????---??? ?? -22λλλ 由于是微元段，认为各处温度是相同的 整理得：()f L t t hU dx t d A -=22λ ? ()f L t t A hU dx t d -=λ2 2 记 L A hU m λ= 1/m 则 )(22 2f t t m dx t d -= 从导热微分方程式来分析上述问题： 这里的肋片导热是常物性、有内热源、一维、稳态导热。（为什么按有内热源来处理？原因：肋片表面与周围流体的对流换热，可表示为：)(f s s t t h y t -=??-λ ，而根据分析作出的简化是：)(x f t = ? 0=??y t ，这样得出：f s t t =，这是不合理的。所以在这种情况下，无法用导热与对流换热间的关系来描述出对流换热。为了反映出这部分对流换热的热量，这时可以把对流换热看作是与导热同时存在的内热源，对流换热是从肋片带走热量，所以应为负的内热源。） 导热微分方程为： 022=+λ V q dx t d 由于肋片沿x 方向的对流换热量是变化的，所以内热源强度也沿x 方向变化。对微元段dx 分析其内热源强度V q ： 该段的对流换热量为：dx U t t h f c ??-=Φ)( c Φ是正值 则微元段 内热源强度为：)()(f L l f c V t t A hU dx A dx U t t h dV q --=???--=Φ-= 则描述等截面直肋的导热微分方程式为： ())(2 2 2f f L t t m t t A hU dx t d -=-=λ 边界条件： ?? ? ? ? ====0 00l x x dx dt t t 引入过余温度θ：f t t -=θ 则上述方程变为：??? ?? ????=====00022 2l x x dx d m dx d θθθθθ 解得等截面直肋温度分布为： [])()(0 )()(0ml ch x l m ch e e e e ml ml x l m x l m -=++=----θθθ 令l x =，得肋端的温度为：) (0 ml ch l x θθ = = 据能量守恒定律知，肋片散入外界的全部热流量都必须通过 x=0 处的肋基截面。据傅里叶定律得知通过肋片散入外界的热流量为： =-=Φx L dx d A θ λ 而 )(/)]([0ml ch x l m sh m dx d --=θθ )(00 ml th m dx d x θθ-== )()()]([000ml th A hU ml th mA ml th m A L L L θλθλθλ==--=Φ W 几点说明： （1）上述推导中忽略了肋端的散热。对于一般工程计算，尤其高而薄的肋片，可以获得较精确的结果。 若必须考虑肋端对流散热时，可采用一种简便的方法：即用假象高度2 'δ + =l l 代替实际 高度l ，然后仍认为端面是绝热时的计算式来计算肋片的散热量。这种想法是基于：为了照顾端面的散热而把端面面积铺展到侧面上去。 （2）上述分析近似认为温度场是一维的。对于肋片，当05.0/≤=λδh Bi 时，这样分析引起的误差不超过1%。对于短而厚的肋片，温度场是二维的，上述算式不适用。实际上，肋片表面上表面传热系数h 不是均匀一致的，这时需要用数值解法。 （3）敷设肋片不一定就能强化传热，只有满足一定的条件才能增加散热量。设计肋片时要注意这一点。（参考《传热学》俞佐平等编） [例2-6] 一铁制的矩形直肋，厚度mm 5=δ，高度mm l 50=，宽度m L 1=。已知肋片材料的导热系数)/(58K m W ?=λ，肋表面与周围介质之间的表面传热系数)/(122K m W h ?=，肋基的过余温度℃800=θ。求肋片的散热量和肋端的过余温度。 [解] 05.0001.058 005 .012<=?= = λ δ h Bi ，因此可以用解析法进行计算。 1 求肋片的散热量 m L m 1005.0=<=δΘ，L L U 2)(2≈+=∴δ，δL A L = 则 m h A hU m L /110.9005 .05812 22=??=== λδλ 假象的肋高：m l l 0525.00025.005.02 '=+=+ =δ 478.00525.010.9'=?=ml 4446.0)478.0()('==th ml th W ml th mA L 86.934446.080005.01.958)(0=????==Φθλ 2 求肋端过余温度 455.005.010.9=?=ml 105.1)455.0()(==ch ml ch ℃ml ch l 4.72105 .180 ) (0 == = θθ 二、肋片效率 （fin efficiency ） 由等截面直肋的导热分析知，肋片表面温度从肋基至肋端是逐渐降低的。所以肋片表面的平均温度m t 必然低于肋基温度0t 。肋片表面平均温度的高低，直接影响着肋片表面的对流换热量。于是，提出了一个如何评价换热壁面加肋后的散热效果问题。肋片效率就是衡量肋片散热有效程度的指标。 肋片效率定义：肋片的实际散热量Φ与假设整个肋表面都处于肋基温度0t 时的理想散热量0Φ的比值。 00)()(θθηm f f m f t t hUl t t hUl = --=ΦΦ= 对等截面直肋， ? ??=Φ=Φ000)(θθλhUl ml th mA L ? ml ml th l m ml th m hUl ml th m A L f ) ()()(2 0===ΦΦ=λη m θ的计算：)()]([1)()()]([11000000ml th ml x l m sh m ml lch dx ml ch x l m ch l dx l l l l m θθθθθ=-???? ??-?=-== ?? 可见，肋片效率是小于1的。 影响肋片效率的因素有：肋片的几何形状和尺寸、肋片材料的导热系数、肋片表面与周围介质的表面传热系数。 f η随ml 的变化情况： 当7.2=ml 时，99.0)(>ml th ，7.2>ml 时，)(ml th 的值变化不大，趋于1。 ? 这时f η可以认为与ml 成反比关系（ml f 1 = η），ml ↑，f η↓。当m 值一定时，随着肋片高度l 的增加，开始散热量增加迅速，后来逐渐减小，最后趋于已渐进值。（当ml 超过2.7后，增加l 对散热量没有作用，这反映出肋片高度增大到一定程度后，如果再继续增高，就会导致肋片效率的急剧降低，达不到进一步增大肋片散热量的效果） ml 数值大的肋片，其肋端的过于温度低 ? 肋片表面的m t 低 ? 肋片的效率低。 ml 数值较小时，肋片具有较高的效率（参看课本图2-16）。所以，在l 一定时，m 取较小值是有利的。 L A hU m λ=，m 小 ? ?? ? ??取小值一定时，和在较大的材料 选用L A U h λλ L A U 取决于肋片的几何形状和尺寸。 %80>f η的肋片经济适用。 （在某些场合下，采用变截面肋片，可提高f η，也减轻肋片重量） 实际上，一般将f η与ml 的关系绘成图。这里给出了矩形直肋和等厚度环肋的f η与ml 的关系图。图中横坐标已进行了变换。 2 3 22l f h l h l A hU ml L λλδλ=≈= δδδl f L A L L U L =??? ? ???=≈+=2)(2Θ是肋片的纵剖面面积。 对等截面直肋：21 23 )2(f h l ml c c λ=，其中，2δ+=l l c ，δc l f = 对等候度环肋（剖面为矩形）： 21 23 )2(f h l ml c c λ=，其中，2δ+=l l c ，c c l r r +=12，δc l f = 计算肋片散热量Φ的一般步骤： （1）根据参数计算f η（按公式或者查图表） （2）计算0Φ（假定肋表面所有温度与肋基温度相等） （3）0Φ?=Φf η 设计肋片：选择形状、计算；考虑质量、制造的难易程度、价格、空间位置的限制等。 [思考题] 1 肋片高度增加引起两种效果：肋效率下降及散热表面增加。因而有人认为，随着肋片高度的增加会出现一个临界高度，超过这个高度后，肋片导热热流量反而会下降。试分析这一观点的正确性。 答：这一观点是不正确的，计算公式表明，肋片散热量：)(0ml th A hU L θλ=Φ，与ml 的的双曲正切成正比，而双曲正切是单调增加函数，所以散热量不会随高度增加而下降。 2 对于一等截面直肋，设肋根温度为0t ，周围介质温度为f t ，且f t t >0。试定性画出沿肋高方向的温度分布，并扼要分析在设置肋片时，肋片高度是否越长越好。 答：肋片温度随高度成指数关系下降。肋片高度增加时散热表面增加，单肋效率下降，虽然换热量也在增加，但因换热器的体积、重量和成本增加，肋片并非越高越好。肋片散热量与ml 的双曲正切成正比，而双曲正切是以1为极值的单调增加函数，ml 为1.5时，其值已超过0.9。 第五节 通过接触面的导热 两个固体表面直接接触时，即使宏观上看来是非常平整的表面，他们的表面也仍是粗糙的，是点接触，而非面接触（接触只发生在一些离散的接触面上）。这样就给导热带来额外的热阻——接触热阻。 接触热阻(thermal contact resistence)：由于接触面间的不密实而产生的附加热阻。 接触热阻导致在接触面上出现温差，对传热不利。 接触热阻：Φ ?= c c t R K/W ① 在Φ不变时，↑c R ，则↑?c t ② 在c t ?不变时，↑c R ，则↓Φ 1 产生并影响接触热阻的主要因素： （1） 接触表面的粗糙度； （2） 表面接触时施加压力的大小； （3） 两接触面之间形成的空隙中气体的热物性。 2 减小接触热阻的措施： （1）减小接触表面的粗糙度； （2）增加接触压力； （3）在两接触表面之间加一层具有高导热系数和高延展性的材料； （4）在接触面之间涂以具有良好导热性的油脂。（减小气体存在空间）。 当接触热阻与接触固体的导热热阻相当时，应考虑接触热阻对导热过程的影响。 思考题： 1、一常物性、无内热源的单层长圆筒壁，内、外半径分别为1r 和2r ，其内、外表面分别维持均匀恒定的温度1w t 和2w t 。试分别就1w t >2 w t t t x A B t w t w 精确解： p=[1,5,10]; x=0:1/19:1; for i=1:1:3 for j=1:1:20 y(i,j)=(exp(p(1,i)*19*x(1,j))-1)/(exp(p(1,i)*19)-1); end plot(x,y(i,:)); hold on ; end 由题对中心差分、一阶迎风、混合格式进行模块编程： 他们之间可以通用，只需更改ae 关于p 的函数即可： 程序如下: (1)中心差分 p=[1,5,10]; for i=1:1:3 ae=1-0.5*p(1,i); x/L (Φ-ΦL )/(Φ0-ΦL ) 精确解图像 aw=p(1,i)+ae; ap=ae+aw; for i=1:1:18 for j=1:1:20 a(i,j)=0; end end for i=1:1:18 j=i; a(i,j)=aw; a(i,j+1)=-ap; a(i,j+2)=ae; end for i=1:1:17 n=i+1; for m=i:-1:1 b(1,1)=a(m,n); a(m,n)=-a(i+1,n)/a(i+1,n)*b(1,1)+a(m,n); a(m,n+1)=-a(i+1,n+1)/a(i+1,n)*b(1,1)+a(m,n+1); a(m,n+2)=-a(i+1,n+2)/a(i+1,n)*b(1,1)+a(m,n+2); end end F(1)=0; F(20)=1; F(19)=(-a(1,20)*F(20)-a(1,1)*F(1))/a(1,19); for i=2:1:18 F(i)=(-a(i,20)*F(20)-a(i,19)*F(19))/a(i,i); end x=0:1/19:1; y(1,:)=F; plot(x,y); hold on end 取步长δx=0.02。已知x=0,Φ=0;x=1,Φ=1.令k=ρｕ/Γ计算结果图表： 程序及数据结果： 追赶法： #include a[i]=2+0.02*k; b[i]=4; c[i]=2-0.02*k; f[i]=0; } tdma(a,b,c,f,x); for(i=0;i 4-1 解：采用区域离散方法A 时；网格划分如右图。内点采用中心差分 23278.87769.9 T T T === 22d T T=0dx - 有 i+1i 12 2+T 0i i T T T x ---=? 将2点,3点带入 32122 2+T 0T T T x --=? 即321 209T T -+= 432322+T 0T T T x --=?4321322+T 0T T T x --=? 即4 321 209 T T T -+-= 边界点4 （1）一阶截差 由x=1 1dT dx =，得 431 3 T T -= （2）二阶截差 11B M M q x x x T T S δδλλ -=++ 所以 434111. 1. 36311 T T T =++ 即 431 22293 T T -= 采用区域离散方法B 22d T T=0dx - 由控制容积法 0w e dT dT T x dT dT ????--?= ? ????? 所以代入2点4点有 322121011336 T T T T T ----= 即 239 028T T -= 544431011363 T T T T T ----= 即 34599 02828T T T -+= 对3点采用中心差分有 432 32 2+T 013T T T --=?? ??? 即 23499 01919 T T T -+= 对于点5 由x=1 1dT dx =，得 541 6 T T -= （1）精确解求左端点的热流密度 由 ()2 1 x x e T e e e -= -+ 所以有 ()22 20.64806911x x x x dT e e q e e dx e e λ -====- +=-=++ （2）由A 的一阶截差公式 21 0.247730.743113 x T T dT q dx λ =-=-= =?= （3）由B 的一阶截差公式 0 0.21640 0.649213 x dT q dx λ =-=-= = （4）由区域离散方法B 中的一阶截差公式： 210.108460.6504()B B T T dT dx x δ-?? ==?= ? ?? 通过对上述计算结果进行比较可得：区域离散B 有控制容积平衡法建立的离散方程与区域离散方程A 中具有二阶精度的格式精确度相当！ 4-3 解：将平板沿厚度方向3等分，如图 参考文献 [1]金国淼等.除尘设备[M].北京:化学工业出版社,2002:1-300 [2]Louis E. Stein, Alex. C. Hoffmann.旋风分离器-原理、设计和工程应用 [M].北京,化学工业出版社,2004:1-78 [3]国家环保局标准处.中华人民共和国国家标准环境空气质量标准[J],油气田环境保护,1996(04 ) [4]姚玉英,黄凤廉,陈常贵等.化工原理[M].天津:天津大学出版社,1999:138 [5]舒帆.影响旋风除尘器除尘效率的因素分析[J],粮食加工.2008, 33 (3):73-75 [6]韩占忠,王敬,兰小平.FLUENT流体工程仿真计算实例与应用[M].北京：北京理工大学出版社,2004:20 [7]魏志军,张平.旋风分离器气相流场的数值模拟[J].北京理工大学学报.2000, 20 (5):19-21 [8]嵇鹰,张红波,田耀鹏等.进口位置对旋风分离器特性影响的数值模拟[J].金属矿山,2008, 387 (3):127-129 [9]岑可法,倪明江,骆仲泱等.循环流化床锅炉理论设计与运行[M].北京:中国电力出版社, 2002:511-540 [10]陈明绍，吴光兴，张大中等.除尘技术的基本原理与应用[M].北京：中国建筑工业出版社，1981:333-518 [11]钱付平,章名耀.基于边界层理论旋风分离器分离效率的改进模型[J],中国电机工程学报.2007, 27 (5):71-74 [12]Hoffmann A C, Stein L E. Gas cyclones and twirl tubes：principles，design and operation [M]. Springer-Verlag，Berlin，Heidelberg，2002，169． [13]Leith D, Licth W. The collection efficiency of cyclone type particle collector. A new theoretical approach[J]. AIChE Symp Series,1972，126 (68)：196-206． [14]Obermair S，Woisetschlager J，Staudinger G．Investigation of the flow pattern in different dust outlet geometries of a gas cyclone by laser Doppler anemometry[J].Powder Technology，2003，2-3 (138):239-251 [15]Zhao Bingtao．Development of a new method for evaluating cyclone 【2-1】一食品冷藏室由内层为19 mm 厚的松木，中层为软木层，外层为51 mm 厚的混凝土所组成。内壁面温度为-17.8 ℃,混凝土外壁面温度为29.4 ℃。松木、软木和混凝土的平均热导率分别为, 3, W/(m ·K)，要求该冷藏室的热损失为15W/m 2。求所需软木的厚度及松木和软木接触面处的温度。 解：三层平壁的导热。 1）所需软木的厚度2b 由 ∑=-=3141i i i b T T q λ 得 151 .0019.00433.0762.0051.08.174.29152+++=b 解得： m b 128.02= 2）松木和软木接触面处的温度3T 由 151 .0019 .08.17153+==T q 解得：9.153-=T ℃ 解题要点：多层平壁热传导的应用。 【2-2】为减少热损失，在外径为150 mm 的饱和蒸汽管道外加有保温层。已知保温材料的热导率λ=+ 198 T(式中T 为℃),蒸汽管外壁温度为180 ℃,要求保温层外壁温度不超过50 ℃,每米管道由于热损失而造成蒸汽冷凝的量控制在1×10-4 kg/(m ·s)以下,问保温层厚度应为多少(计算时可假定蒸汽在180 ℃下冷凝)。 解：保温层平均热导率为： )./(126.02 501801098.1103.04K m W =+??+=-λ 由于本题已知的是蒸汽管道外壁面温度，即保温层内壁面温度，故为一层导热。 由 )()(21 221r r Ln T T L Q -=λπ 得： )()(21 221r r Ln T T L Q -=πλ （1） 式中：m W L Wr L Q /9.2011 103.20191013 4=???==- 将其及其它已知数据代入式（1）得： )075 .0()50180(126.029.2012r Ln -??=π 解得：m r 125.02= mm m 5005.0075.0125.0==-=∴δ壁厚 解题要点：单层圆筒壁热传导的应用。 【2-8】烤炉内在烤一块面包。已知炉壁温度为175 ℃，面包表面的黑度为,表面温度为100 ℃,表面积为 5 m 2，炉壁表面积远远大于面包表面积。求烤炉向这块面包辐射 传递的热量。 解：两物体构成封闭空间，且21S S <<，由下式计算辐射传热量： W T T S Q 0.65)448373(0645.085.01067.5) (448424111012-=-????=-=-εσ 负号表示炉壁向面包传递热量。 解题要点：辐射传热的应用，两个灰体构成的封闭空间。 【2-10】在逆流换热器中，用初温为20 ℃的水将1.25 kg/s 的液体［比热容为 kJ/(kg ·K)、密度为850 kg/m 3 ］由80 ℃冷却到30 ℃。换热器的列管直径为Φ25 mm ×2.5 mm,水走管内。水侧和液体侧的对流传热系数分别为850 W/（m 2·K ）和1 700W/（m 2·K ）,污垢热阻可忽略。若水的出口温度不能高于50 ℃，求水的流量和换热器的传热面积。 第二章 稳态导热 本章重点：具备利用导热微分方程式建立不同边界条件下稳态导热问题的数学模型的能力 第一节 通过平壁的导热 1-1 第一类边界条件 研究的问题： （1）几何条件：设有一单层平壁，厚度为δ，其宽度、高度远大于其厚度（宽度、高度是厚度的10倍以上）。这时可认为沿高度与宽度两个方向的温度变化率很小，温度只沿厚度方向发生变化。（属一维导热问题） （2）物理条件：无内热源，材料的导热系数λ为常数。 （3） 边界条件：假设平壁两侧表面分别保持均匀稳定的温度 1w t 和2w t ，21w w t t >。（为第一类边界条件，同时说明过程是稳态的） 求：平壁的温度分布及通过平壁的热流密度值。 方法1 导热微分方程： 采用直角坐标系，这是一个常物性、无内热源、一维稳态导热问题（温度只在 x 方向变化）。 导热微分方程式为：022=dx t d （2-1） 边界条件为：10w x t t == ， 2w x t t ==δ （2-2） 对式（2-1）连续积分两次，得其通解： 21c x c t += （2-3） 这里1c 、2c 为常数，由边界条件确定 ，解得：?? ???=-= 11221w w w t c t t c δ （2-4） 最后得单层平壁内的温度分布为： x t t t t w w w δ 2 11-- = （2-5） 由于δ 、1w t 、2w t 均为定值。所以温度分布成线性关系，即温度分布曲线的斜率是常数（温度梯度）， const t t dx dt w w =-=δ 1 2 （2-6） 热流密度为：)(21w w t t dx dt q -=-=δ λ λ 2/m W （2-7） 若表面积为 A, 在此条件下 , 通过平壁的导热热流量则为 : t A qA ?==Φδ λ W （2-8） 考虑导热系数随温度变化的情况： 对于导热系数随温度线形变化，即)1(0bt +=λλ，此时导热微分方程为：0=?? ? ??dx dt dx d λ 解这个方程，最后得： ?? ? ???++-+?? ? ?? +=+)(211212121121 122w w w w w w t t b x t t bt t bt t δ 或 x t t t t b b t b t w w w w w δ 12211)(2112 2-??????+++??? ??+=??? ??+ 说明：壁内温度不再是直线规律，而是按曲线变化。 对上式求导得：??? ?????+??? ??-=)1/(222bt dx dt b dx t d 因为 01>+bt ，02 >?? ? ??dx dt 所以 0>b ? 02 2 第二章 3.导热系数为常数的无内热源的平壁稳态导热过程，试问，若平壁两侧给定边界条件Tw1和Tw2，为什么这一导热过程的温度分布与平壁的材料无关？相同的平壁厚度，不同的平壁材料，仍给定第一类边界条件，热流密度是否相同？ （1）温度分布为 12 1w w w t t t t x δ -=- (设12w w t t >) 其与平壁的材料无关的根本原因在 coust λ=（即常物性假设） ，否则t 与平壁的材料有关 （2）由 dt q dx λ =- 知，q 与平壁的材料即物性有关 6.同上题，若已知边界条件为第三类，即已知Tf1，h1，Tf2，h2.试倒通过空心球壁热量的计算公式和球壁的传热热阻。 9.某教室有一层厚度为240mm 的砖层和一厚度为20mm 的灰泥构层。现安装空调设备，并在内表面加贴一层硬泡某塑料，是导入室内的热量比原来少了80%。已知砖的导热系数λ=0.7W/(m*k)，灰泥为λ=0.58W/(m*k)，硬泡某塑料的导热系数为λ=0.06W/(m*k)，试求出硬泡某塑料厚度。 已 知 ： 12240,20mm mm δδ==， 120.7/(),0.58/()W m k W m k λλ=?=? 3210.06/(),0.2W m k q q λ=?= 求：3δ 解： 设两种情况下的内外面墙壁温度12w w t t 和保持不变， 且12w w t t > 由题意知：12 112 12 w w t t q δδλλ-= + 12 23 12123 w w t t q δδδλλλ-= ++ 再由： 210.2q q =，有 12 12 3 12 1212 123 0.2 w w w w t t t t δδδδδλλλλλ--=+++ 22 131 3 1 2 tw 1 q tw 2 1 1 λ1 2 λ2 tw 1 tw 2 q 1 1λ1 2λ 2 3λ 3 传热学杨世铭第四版第二章答案 第二章 思考题 1 试写出导热傅里叶定律的一般形式，并说明其中各个符号的意义。 答：傅立叶定律的一般形式为：n x t gradt q ??-=λλ＝－，其中：gradt 为空间某点的温度梯度；n 是通过该点的等温线上的法向单位矢量，指向温度升高的方向；q 为该处的热流密度矢量。 2 已知导热物体中某点在x,y,z 三个方向上的热流密度分别为y x q q ,及z q ，如何 获得该点的 热密度矢量？ 答：k q j q i q q z y x ?+?+?=，其中k j i ,,分别为三个方向的单位矢量量。 3 试说明得出导热微分方程所依据的基本定律。 答：导热微分方程式所依据的基本定律有：傅立叶定律和能量守恒定律。 4 试分别用数学语言将传热学术语说明导热问题三种类型的边界条件。 答：① 第一类边界条件：)(01ττf t w =>时， ② 第二类边界条件：)()(02τλτf x t w =??->时 ③ 第三类边界条件：)()(f w w t t h x t -=??-λ 5 试说明串联热阻叠加原则的内容及其使用条件。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各串联环节热阻的和。使用条件是对于各个传热环节的传热面积必须相等。 7.通过圆筒壁的导热量仅与内、外半径之比有关而与半径的绝对值无关，而通过球壳的导热量计算式却与半径的绝对值有关，怎样理解？ 答：因为通过圆筒壁的导热热阻仅和圆筒壁的内外半径比值有关，而通过球壳的导热热阻却和球壳的绝对直径有关，所以绝对半径不同时，导热量不一样。 6 发生在一个短圆柱中的导热问题，在下列哪些情形下可以按一维问题来处理？ 答：当采用圆柱坐标系，沿半径方向的导热就可以按一维问题来处理。 8 扩展表面中的导热问题可以按一维问题来处理的条件是什么？有人认为，只要扩展表面细长，就可按一维问题来处理，你同意这种观点吗？ 答：只要满足等截面的直肋，就可按一维问题来处理。不同意，因为当扩展表面的截面不均时，不同截面上的热流密度不均匀，不可看作一维问题。 9 肋片高度增加引起两种效果：肋效率下降及散热表面积增加。因而有人认为，随着肋片高度的增加会出现一个临界高度，超过这个高度后，肋片导热热数流量反而会下降。试分析这一观点的正确性。 第二章 思考题 1 试写出导热傅里叶定律的一般形式，并说明其中各个符号的意义。 答：傅立叶定律的一般形式为： n x t gradt q ??-=λλ＝－，其中：gradt 为空间某点的温度梯度；n 是通过该点的等温线上的法向单位矢量，指向温度升高的方向；q 为该处的热流 密度矢量。 2 已知导热物体中某点在x,y,z 三个方向上的热流密度分别为y x q q ,及z q ，如何获得该点的 热密度矢量？ 答：k q j q i q q z y x ?+?+?=，其中k j i ,,分别为三个方向的单位矢量量。 3 试说明得出导热微分方程所依据的基本定律。 答：导热微分方程式所依据的基本定律有：傅立叶定律和能量守恒定律。 4 试分别用数学语言将传热学术语说明导热问题三种类型的边界条件。 答：① 第一类边界条件：)(01ττf t w =>时， ② 第二类边界条件： )()( 02τλτf x t w =??->时 ③ 第三类边界条件： )()( f w w t t h x t -=??-λ 5 试说明串联热阻叠加原则的内容及其使用条件。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各串联环节热阻的和。使用条件是对于各个传热环节的传热面积必须相等。 7.通过圆筒壁的导热量仅与内、外半径之比有关而与半径的绝对值无关，而通过球壳的导热量计算式却与半径的绝对值有关，怎样理解？ 答：因为通过圆筒壁的导热热阻仅和圆筒壁的内外半径比值有关，而通过球壳的导热热阻却和球壳的绝对直径有关，所以绝对半径不同时，导热量不一样。 6 发生在一个短圆柱中的导热问题，在下列哪些情形下可以按一维问题来处理？ 答：当采用圆柱坐标系，沿半径方向的导热就可以按一维问题来处理。 8 扩展表面中的导热问题可以按一维问题来处理的条件是什么？有人认为，只要扩展表面细长，就可按一维问题来处理，你同意这种观点吗？ 答：只要满足等截面的直肋，就可按一维问题来处理。不同意，因为当扩展表面的截面不均时，不同截面上的热流密度不均匀，不可看作一维问题。 9 肋片高度增加引起两种效果：肋效率下降及散热表面积增加。因而有人认为，随着肋片高度的增加会出现一个临界高度，超过这个高度后，肋片导热热数流量反而会下降。试分析这一观点的正确性。 答：错误，因为当肋片高度达到一定值时，通过该处截面的热流密度为零。通过肋片的热流已达到最大值，不会因为高度的增加而发生变化。 10 在式（2-57）所给出的分析解中，不出现导热物体的导热系数，请你提供理论依据。 答：由于式（2-57）所描述的问题为稳态导热，且物体的导热系数沿x 方向和y 方向的数值相等并为常数。 11 有人对二维矩形物体中的稳态无内热源常物性的导热问题进行了数值计算。矩形的一个边绝热，其余三个边均与温度为f t 的流体发生对流换热。你能预测他所得的温度场的解吗？ 答：能，因为在一边绝热其余三边为相同边界条件时，矩形物体内部的温度分布应为关于绝热边的中心线对称分布。 习题 主讲陶文铨 西安交通大学能源与动力工程学院热流中心CFD-NHT-EHT CENTER 2009年9月7日,西安 数值传热学 第一章绪论 课程简介 1. 教材－《数值传热学》第二版，2001 2. 学时－45学时理论教学；10学时程序教学 3. 考核－平时作业/计算机大作业： 考试－40/60；考查－60/40 4. 方法－开放，参与，应用 5. 助手－郭东之，周文静，李兆辉 有关的主要国外期刊 1.Numerical Heat Transfer, Part A-Applications; Part B- Fundamentals 2.International Journal of Numerical Methods in Fluids. https://www.wendangku.net/doc/ef7704082.html,puter & Fluids 4.Journal of Computational Physics 5.International Journal of Numerical Methods in Engineering 6.International Journal of Numerical Methods in Heat and Fluid Flow https://www.wendangku.net/doc/ef7704082.html,puter Methods of Applied Mechanics and Engineering 8.Engineering Computations 9.Progress in Computational Fluid Dynamics 10. Computer Modeling in Engineering & Sciences (CMES) 11.ASME Journal of Heat Transfer 12.International Journal of Heat and Mass Transfer 13.ASME Journal of Fluids Engineering 14.International Journal of Heat and Fluid Flow 15.AIAA Journal 第一章: 1-1 对于附图所示的两种水平夹层，试分析冷、热表面 间热量交换的方式有何不同？如果要通过实验来测定夹层中流体的导热系数，应采用哪一种布置？ 解：（a ）中热量交换的方式主要有热传导和热辐射。 （b ）热量交换的方式主要有热传导，自然对流和热辐射。 所以如果要通过实验来测定夹层中流体的导热系数，应采用（ a ）布置。 1-2 一炉子的炉墙厚13cm ，总面积为20m 2 ，平均导热系数为 1.04w/m 〃k ，内外壁温分别是520 ℃及50 ℃。试计算通过炉墙的热损失。如果所燃用的煤的发热量是 2.09 ×10 4 kJ/kg ，问每天因热损失要用掉多少千克煤？ 解：根据傅利叶公式 每天用煤 1-3 在一次测定空气横向流过单根圆管的对流换热实验中，得到下列数据：管壁平均温度t w = 69 ℃，空气温度t f = 20 ℃，管子外径d= 14mm ，加热段长80mm ，输入加热段的功率8.5w ，如果全部热量通过对流换热传给空气，试问此时的对流换热表面传热系数多大？ 解：根据牛顿冷却公式 1-4宇宙空间可近似的看作0K 的真空空间。一航天器在太空中飞行，其外表面平均温度为250K ，表面发射率为0.7 ，试计算航天器单位表面上的换热量？ 解：航天器单位表面上的换热量 1-5附图所示的空腔由两个平行黑体表面组成，孔腔内抽成真空，且空腔的厚度远小于其高度与宽度。其余已知条件如图。表面 2 是厚δ= 0.1m 的平板的一侧面，其另一侧表面 3 被高温流体加热，平板的平均导热系数λ=17.5w/m ? K ，试问在稳态工况下表面3 的t w3 温度为多少？ 解： 表面1 到表面2 的辐射换热量= 表面2 到表面3 的导热量 第二章: 西安交通大学西安交通大学《《数值传热学数值传热学》》课程大作业 20140114 一. 题目 （1） 百叶窗翅片的二维模型如图1 所示。在流动与换热已经进入周期性充分发展的阶段，可以取 出一个翅片单元进行传热与流动阻力的分析计算。在稳态，层流，常物性，翅片温度恒定的条件下，对于表1给定的几何尺寸，进行Re ＝10－500 范围内的数值模拟，揭示每个计算单元的平均Nu 数与阻力系数f 与Re 的关系； Nu ，f 以及Re 定为：1 12()Re ;;0.5p m m m dp dx L u L h L f Nu u νρλ==?= 其中m u 为来流平均速度；m h 为每块条片的平均换热系数。 表1 几何参数 L1/mm Tp/mm Lp/mm Delta/mm /θ 30 18.6 30 1.5 25 图1 百叶窗翅片二维模型 图2 阶梯型逼近 二. 建议建议与要求与要求 1. 为便于处理流固耦合问题，计算可对图1中打阴影线的区域进行； 2. 可采用图2 所示的阶梯型网格处理倾斜的翅片； 3. 按照《西安交通大学学报》的论文格式撰写本报告； 4. 2014年4月30号前交课程论文到东三楼204房间。 三. 参考文献 [1] 陶文铨编著，数值传热学（第二版），2001， 西安交通大学出版社，节11.2 [2] Wang L B, Tao, W Q. Numerical analysis on heat transfer and fluid flow for arrays of non-uniform plate length aligned at angles to the flow direction. Int J Numerical Methods for Heat and Fluid Flow , 1997, 7(5,6):496 [3] Gong L. Li Z Y, He Y L, Tao W Q. Discussion on numerical treatment of periodic boundary condition for temperature. Numerical Heat Transfer, Part B , 2007, 52(5):429-448 热物理过程的数值模拟Numerical Simulation of Thermophysics Process 讲稿 主讲：李隆键 第一章概论 1．1流动与传热过程的予测方法及特点 流动、传热、燃烧问题是热工类各专业和机械类动力机械专业所研究和解决的主要问题之一，燃烧问题实际上是有化学反应的流动与传热问题，推而广之，在所有热物理过程中，几乎都涉及到流动、传热问题。 预测的重要性： ①在规定设计参数的相应的结构下，热物理过程是否满足要求，达到预定的指 标？要预测； ②优化设计，不同方案的比较，要预测； ③减少设计、生产、再设计和再生产的费用； ④减少设计更改； ⑤减少试验和测量次数。 问题的核心：速度场、温度场（传热量）、浓度场等。 一、热物理问题的予测方法：理论分析法、实验测定、数值模拟 1、理论分析 以数学分析为基础，求解描述热物理过程的定解问题，获得函数形式的解，表示求解区域内物理量连续分布的场（速度场、温度场、浓度场……）。 控制方程+单值条件（数学模型）→理论解（分析解，解析解） 根据解的准确程度，又可再分为： （1）精确分析解（严格解） 特点：函数形式的解；它在求解区域精确地满足定解问题。 具体解法：直接积分法、分离变量法、积分变换法、热源法、映射法。 （2）近似分析解法 特点：函数形式的解，在求解区域上近似地满足定解问题（但在总量上满足相应的守恒原理，动量守恒、动量守恒、能量守恒、质量守恒）。 具体解法：积分法（从积分方程出发） 变分近似解法 摄动法（从微分方程出发） 2、实验测定 （1）纯实验法 （2）相似理论实验法：同类相似，减少变量数目→减少工作量，得到规律性结 数值传热学（Numerical Heat Transfer,NHT）又称计算传热学(Computational Heat Transfer,CHT)，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等），用一系列有限个离散点(称为节点)上的值的集合来代替，通过一定的原则建立起这些离散点变量值之间关系的代数方程（称为离散方程，discretization equation）,求解所建立起来的代数方程已获得求解变量的近似值。 一、数值传热学的研究作用与地位 数值传热学在最近20年中得到飞速的发展，除了计算机硬件工业的发展给它提供了坚实的物质基础外，还主要因为无论分析的方法或实验的方法都有较大的限制，例如由于问题的复杂性，既无法做分析解，也因费用的昂贵而无力进行实验测定，而数值计算的方法正具有成本较低和能模拟复杂或较理想的过程等优点。经过一定考核的数值计算软件可以拓宽实验研究的范围，减少成本昂贵的实验工作量。在给定的参数下用计算机对现象进行一次数值模拟相当于进行一次数值实验，历史上也曾有过首先由数值模拟发现新现象而后由实验予以证实的例子。在这里要指出对数值模拟结果准确度应持正确认识。计算机本身不能创造信息，发现规律，它只是把人们送入的信息按照计算者所选定的规律进行处理，加工而已。但一旦建立了实际问题合理的数学模型，数值模拟又能发挥很大的作用。由于它本身的一些固定优点，它以发展成为工业界进行CAD/CAM及过程控制的重要手段,在多种工程领域中得到广泛应用。例如：叶轮机器粘性三元流体的计算，电站锅炉炉堂内流场与温度场的模拟；大型初见凝固过程中温度场的预测； 热流问题的数值计算 物理问题数值解的基本思想 把原来在空间与时间坐标中连续的物理量的场 (如速度场,温度场,浓度场等),用一系列有限 个离散点(称为节点,node)上的值的集合来代替; 通过一定的原则建立起这些离散点上变量值之间关 系的代数方程 (称为离散方程,discretization 大规模科学计算的重要性 传热与流动问题数值计算是应用计算机求解热量传 递过程中的速度场,温度场等的分支学科,是大规模 科学计算的重要组成部分,其重要性不言而喻. 2005年美国总统顾问委员会向美国总统提出要大 力发展计算科学以确保美国在世界上的竞争能力. 波音公司实现了对航空发动机的网格数达10亿量 级的直接数值模拟,以研究所设计发动机的性能. 现代科学研究的三大基本方法及其关系 课程简介 1、二维非稳态导热微分方程：S Y T X T t T p +??+??=)2222c （λδδρ。对于时间步进（x 方向，y 方向）及空间而言，该方程为何种类型的方程？ 解： 将二维非稳态导热微分方程化为：0c 2222=+-??+??S t T Y T X T p δδρλλ （1）x 方向：0,0,a ===c b λ。则：04b 2=-=?ac ，所以该二维非稳态导热方程为抛物型方程。 （2）y 方向：0,0,a ===c b λ。则：04b 2=-=?ac ，所以该二维非稳态导热方程也为抛物型方程。 （3）对于空间而言，二维非稳态导热方程可知：,0,a b c λλ===则：2240b ac λ?=-=-<，所以该二维非稳态导热方程为椭圆型方程。 2、 （补充不可压、常物性的条件。写出守恒型和非守恒型控制方程，并推导二者关系。） 解：由题可知，该流体为不可压缩、常物性流体，而且是有内热源的二维问题。 守恒型控制方程： 质量守恒方程：0=??+??y v x u ； 由于流体自身条件，使得0==v u S S ，得 动量守恒方程：()()??? ? ????+??+??-=??+??22221y u x u v x p y vu x uu ρ ()()??? ? ????+??+??-=??+??22221y v x v v y p y vv x uv ρ ； 能量守恒方程：()()T S y T x T a y vT x uT +??? ? ????+??=??+??2222 . 非守恒型控制方程： 质量守恒方程：无非守恒型 动量守恒方程：??? ? ????+??+??-=??+??22221y u x u v x p y u v x u u ρ ??? ? ????+??+??-=??+??22221y v x v v y p y v v x v u ρ ； 考博详解与指导 考试科目代码考试科目名称参考书目： 1001英语《新世纪研究生英语教材--阅读B,C》戴炜栋，柴小平编，上海外语教育出版社 1002俄语①《基础俄语》(1-3册)北京外语学院编，外语教学与研究出版社②《大学俄语基础教程》(1-3册)张智罗，高等教育出版社 1003日语《新编日语》(1-3册)周平、陈小芬，上海外语教育出版社 1004德语①《大学德语》戴鸣钟，高等教育出版社②《新编大学德语》朱建华编，外语教学与研究出版社，2002年9月第一版 1005法语《法语》(1-3册)马晓宏，外语教育出版社 2001工程流体力学①《工程流体力学》，归柯庭汪军王秋颖，科学出版社，2004年②《工程流体力学》(第二版)，孔珑，中国电力出版社，2007年 2002传热学《传热学》杨世铭，高等教育出版社，2006年 2003计算方法《数值分析》李庆杨等编著，清华大学出版社，2008年 2004高等光学《近代光学》袁一方译，高等教育出版社，1987年 2005物理光学《物理光学》梁铨庭，机械工业出版社 2006传感器技术及应用①《传感器》强锡富主编，机械工业出版社，2004年7月第三版②《非电量电测技术》严钟豪等主编，机械工业出版社，2003年1月第二版 2007激光原理《激光原理及应用》(第1-4章，6章)清华大学出版社 2008普通物理(光学)《普通物理学》(光学部分)程守洙，人民教育出版社 2009仪器电路原理与应用①《仪器电路设计与应用》，郝晓剑等编著，电子工业出版社，2007年6月②《基于运算放大器和模拟集成电路的电路设计》，赛尔吉欧。佛朗哥著西安交通大学出版社，2004年8月第1版 传热学思考题参考答案 第一章： 1、用铝制水壶烧开水时，尽管炉火很旺，但水壶仍安然无恙。而一旦壶内的水烧干后水壶很快就被烧坏。试从传热学的观点分析这一现象。 答：当壶内有水时，可以对壶底进行很好的冷却（水对壶底的对流换热系数大），壶底的热量被很快传走而不至于温度升得很高；当没有水时，和壶底发生对流换热的是气体，因为气体发生对流换热的表面换热系数小，壶底的热量不能很快被传走，故此壶底升温很快，容易被烧坏。 2、什么是串联热阻叠加原则，它在什么前提下成立？以固体中的导热为例，试讨论有哪些情况可能使热量传递方向上不同截面的热流量不相等。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各 串联环节热阻的和。例如：三块无限大平板叠加构成的平壁。例如通过圆筒壁，对于各个传热环节的传 热面积不相等，可能造成热量传递方向上不同截面的热流量不相等。 第二章： 1、扩展表面中的导热问题可以按一维问题处理的条件是什么?有人认为,只要扩展表面细长,就可按一维问题处理,你同意这种观点吗? 答：条件：（1）材料的导热系数，表面传热系数以及沿肋高方向的横截面积均各自为常数（2）肋片温度在垂直纸面方向（即长度方向）不发生变化，因此可取一个截面（即单位长度）来分析（3）表面上的换热热阻远远大于肋片中的导热热阻，因而在任一截面上肋片温度可认为是均匀的（4）肋片顶端可视为绝热。并不是扩展表面细长就可以按一维问题处理，必须满足上述四个假设才可视为一维问题。 2、肋片高度增加引起两种效果:肋效率下降及散热表面积增加。因而有人认为随着肋片高度的增加会出现一个临界高度，超过这个高度后，肋片导热热流量会下降，试分析该观点的正确性。 答：的确肋片高度增加会导致肋效率下降及散热表面积增加，但是总的导热量是增加的，只是增加的部分的效率有所减低，所以我们要选择经济的肋片高度。 第三章： 1、由导热微分方程可知，非稳态导热只与热扩散率有关,而与导热系数无关。你认为对吗？答：错，方程的边界条件有可能与λ有关，只有当方程为拉普拉斯方程和边界条件为第一边界条件时才与λ无关。 2、对二维非稳态导热问题，能否将表面的对流换热量转换成控制方程中的内热源产生的热量？ 答：不能，二维问题存在边界微元和内边界微元，内边界微元不一定与边界换热，所以不存在源项。 第四章： 1、在第一类边界条件下，稳态无内热源导热物体的温度分布与物体的导热系数是否有关？为什么？ 答：无关，因为方程为拉普拉斯方程，边界为第一边界条件均与λ无关。 2、非稳态导热采用显式格式计算时会出现不稳定性，试述不稳定性的物理含义。如何防止这种不稳定性？ 答：物理意义：显示格式计算温度时对时间步长和空间步长有一定的限制，否则会出现不合 数值传热学第二章作业 2—1： POWER=input('POWER=?'); L1=input('L1=?'); M1=input('M1=?'); XL=input('XL=?'); YL=input('YL=?'); for i=2:L1 XF(i)=XL*((i-2)/(L1-2))^POWER; end for j=2:M1 YF(j)=YL*((j-2)/(M1-2))^POWER; end X(1)=0; for i=2:L1-1 X(i)=(XF(i)+XF(i+1))/2; end X(L1)=XF(L1); Y(1)=0; for j=2:M1-1 Y(j)=(YF(j)+YF(j+1))/2; end Y(M1)=YF(M1); for j=2:M1-1 plot(X(1),Y(j),'b.'); plot(X(1),Y(2), 'b.'); plot(X(L1),Y(j),'b.'); hold on end for i=2:L1-1 for j=1:M1 plot(X(i),Y(j),'b.'); hold on end end for i=2:L1 m=[XF(i),XF(i)]; n=[0,M1]; plot(m,n,'b-.'); hold on end for j=2:M1 m=[YF(j),YF(j)]; n=[0,L1]; plot(n,m,'b-.'); hold on end xlabel('x'); ylabel('y'); title('POWER= ') 运行结果如下： 2—3： 解：由2 2 2 1()u 2u u u x x y η ???== =???得： 原方程的守恒形式为： 2 2 2()2u u x y η ??=?? 对方程两端在t ?时间间隔内对其控制容积积分，把可积的部分积出后得： 22()t t s n e w t u u dtdy +?-? ? = 2t t e w t n s u u dtdx y y η+?????????-?? ? ????????? ? ? 选定2 u 随y 而变化的型线，这里取为阶梯式，即在控制容积内沿y 方向不变，则 2222 ()=y ()t t t t s n e w e w t t u u dtdy u u dt +?+?-?-? ? ? 选定2 u 随t 而变化的规律，这里采用阶梯式显式，则 2 2 ()t t e w t y u u dt +??-? = ()()22t t e w u u t y ??-?????? 选定 u y ??随x 而变化的型线，这里取为阶梯式，即在控制容积内沿x 方向不变，则 22t t t t e w t t n s n s u u u u dtdx x dt y y y y ηη+?+?????????????????-=?-???? ? ? ? ?????????? ???????? ? ? 选定 u y ??随t 而变化的规律，这里采用阶梯显式，则 2t t t n s u u x dt y y η+??????????-?? ? ???? ?????? = 2t t n s u u t x y y η?? ??????-???? ? ??????????? 进一步选取u 随x,y 分段线性变化，则 22 2 2 E P e u u u += ， 22 2 w 2 W P u u u +=陶文铨 数值传热学 第二版 第五章 5-2
计算传热学中国石油大学(华东)第四章大作业
数值传热学陶文铨第四章作业
旋风分离器参考文献
第二章 传热习题答案
传热学讲义设计—第二章
传热学第二章答案
最新传热学杨世铭第四版第二章答案
传热学第二章答案
数值传热学chapter_1
传热学_杨茉_部分习题与解答
西安交通大学西安交通大学《 《《《数值传热学数值传热学
热物理过程的数值模拟-计算传热学1
关于数值传热学的调研报告..
数值传热学(课件)-1
Numerical Simulations of Thermal & Fluid Problems
第一章 绪论
主讲 陶文铨
西安交通大学能源与动力工程学院 热流中心 CFD-NHT-EHT CENTER 2007年10月16日, 西安
1/88
equation);求解所建立起来的代数方程以获得所求
解变量的近似解.
2/88
3/88
理论分析
Analytical
实验研究
Experimental
数值模拟
Numerical
4/88
1. 学时- 30学时理论教学;6学时计算机作业 2. 考核- 平时作业/计算机大作业/考试: 20/30/50 3. 方法- 理解,参与,应用 努力将与数学处理相对应的物理背景联系起来理解. 4. 助手- 于乐 5. 参考教材-《计算流体力学与传热学》,中国建筑 工业出版社,1991
5/88数值传热学作业-第一章
上海理工大学博士研究生入学考试参考书目
传热学思考题参考答案(陶文铨第四版)
数值传热学第二章作业