第九章对流传热

1、1第九章 对流传热Chapter 9 Heat Convection29.1 对流传热的基本概念及传热方程Basic Concept for Convection一、特点(Features)n对流传热是研究有流体（气、液及其混合物）存在的传热物质体系中，通过流体的流动（质量团的宏观迁移）产生的热量传热现象及热量传输速率和温度分布的定量分析。n在工程中常见到的也是具有重要的工程意义的对流传热情况就是某种流动流体与固体壁之间的界面对换热。3n我们在绪言中已讲到流体与流体器壁之间的对流换热与许多因素有关，除了与固体和流体的自身的导热能力等因素之外，更重要的还与固体器壁前沿流体的流动强度（速度大小）和

2、流动方式（流动方向、涡流密度，尤其是层流还是紊流）密切相关。而决定或者说影响固体器壁面前沿流体流动强度和流动状态在很大程度上取决于引起流体流动的原因或驱动力。n因此，工程上根据引起流体流动的原因或流动驱动力不同，通常将流动流体与固体器壁之间的对流换热分为强制对流换热和自然对流换热 。4n强制对流换热(Heat Transfer with Forced Convection)： 流动靠外力，如水泵、风机等来驱动 特点：流动速度大、常常是紊流对流换热能力强、hn自然对流换热(Heat Transfer with Natural Convection) ： 由于温度浓度分布不均匀(T, CL )不均

3、匀，在重力场下，会上浮或下沉流动 特点：流动强度直接与 T/n, c/n有关，一般强度不大，通常为层流，对流换热相对弱些，h5 无论是强制对流换热，还是自然对流换热，流动流体与固体壁之 间的对流换热强度，均由牛顿换热来计算，即： q = h (T - Ts) (9-2) 式中： q对流换热热流密度（通量：W/） h对流换热系数（W/） T，Ts流体和固体表面温度（） 6 换热系数h 的影响因素 Factors Influencing Heat Transfer Coefficient1、流体的流速:n流速层流底尽c变薄，热阻变小，导热增强，hn流体内部相对运动加剧，hn在强制对流时，加强流速，

4、耗更多能量。应选取适当流速。72、流体的物性量n，热阻 / 小，h h水= 20h空气 nc载热能力强，热交换强，hn滞止作用大，c厚 ，减弱对流，hn需要综合考虑：如水，粘性大，但c、也大。因此比空气（粘性小）的换热系数大的多。83、壁面几何尺寸、形状、位置：n垂直放置h，水平放置h（顺流动方向放置）9二、热边界层概念 Thermal Boundary Layer如图所示： 当紊流整体速度为V的流体流过一平板时，除产生一速度边界层(x)外，类似的：当紊流流体温度T 与平板表面温度Ts不相等时，由于传热还会产生温度边界层 T(x) 。n解释：紊流温度 T(x) Ts ，接平板后，接触平板处，传

5、热T下降，随流动 x 增加，温度影响层增大，(x)增大。1011n温度边界层的厚度 T 是如下定义的: 当流体（其温度分布不均匀）温度为整体温度 T 的99% 时，即： T (x, y, z) = 0.99T时所对应的距离平板的高度( T =y)，定义为温度边界层厚度。n由于随着流体沿平板的流动距离 x 增加，流体与平板接触时间加长，增加了流体与平板之间的传热量。所以受平板较低温度的影响范围增加。即T(x)随 x 增加，从而T 增加。 n可见为确定T(x)和 h。均需要知道沿平板流动流体中的温度分布T (x, y)。（稳态的）121314 将热边界层三个区间简化成一层流区和一紊流区，在 层流底

6、层只有导热。故： （ 有效热边界层厚度）()sTTTAQ()sQh TTAThT15 三、对流换热系数h的数学计算式Calculation of Heat Transfer Coefficient (92)是计算某一对流换热条件下，传热热通量 q 的基本计算式。可见(92)式的数学形式是十分简单的，但通过前面对固体壁与流动流体对流换热机理方面分析知道，这一过程是十分复杂的 ： n 固相中的导热n 流体中的导热n 流体的流动强度和流动状态n 固体表面几何状况16n因此，尽管(9-2)式的数学形式简单，但把复杂的影响因素都归结到对流换热系数 h 中，即： h = f（ ，流动状态, s, ,cp,

7、界面条件等）n通常影响 h 的因素是十分复杂的。但另一方面，要想应用(9-2)式准确计算 q ，又需要知道某具体对流传热条件下的 h 的精确值。 v17n解决这一问题的办法就是针对某一具体对流换热条件下，分析计算或测出靠近固相壁附近的温度场 T(y)，基于下式计算出对流换热系数 h： q = h (Tf - Ts) 式中：流体导热系数， Ty固体壁附近的的流体温度分布 。 上式是根据粘性流动流体在固相表面上的状态导出的 。0|yTqy 0|()ysThTTy 18n由于流体的粘性0，在固相表面总存在一薄层流体相对固体静止，这样流体与固相壁之间的传热热流量q，毕竟要通过这层相对静止的流体，而且是

8、通过导热方式进行。由傅立叶定律0|yTqy 导19n而 q导 应该等于(9-2)式中的 q，即：n从而得到(9-3)公式：0|()ysTqh TTy 导0|()ysThTTy 209.2 、热量传输方程 傅立叶-克希荷夫导热微分方程 本节将用微元体法导出含有对流条件下的流体中的热量传出方程。n做下列假设： 1) 没有内热源（如化学反应热效应）产生 2) 流体流速不高，由粘性引起的耗散热可忽 略 3) 流体的物性参数、Cp、等不随温度和压 力而变21( VxH )x( VxH )x+ x( VyH )y( VyH )y+ y( VzH )z( VzH )z+ z 图图9-6 进出微元体的热能进出

9、微元体的热能22n根据能量守恒 热量收入热量支出热焓变化1). x 方向的由传导引起的热量输入： y z t qx2). x 方向的热量输出： y z t q x+ x3). y 方向的热量输入： z x t qy4). y 方向的热量输出： z x t q y+ y 5). z 方向的热量输入： y x t qz6). z 方向的热量输出： y x t q z+ z237). 由对流引起的 x 方向的热量输入： y z t ( VxH )x H=CpT单位重量物质在T温度下的热焓8). 由对流引起的 x 方向的热量输出： y z t ( VxH )x+ x9). 由对流引起的 y 方向的热量

10、输入： z x t ( VyH )y10). 由对流引起的 y 方向的热量输出： z x t ( VyH )y+ y11). 由对流引起的 z 方向的热量输入： y x t ( VzH )z12). 由对流引起的 z 方向的热量输出： y x t ( VzH )z+ z2413). 热焓的变化： x y z ( H )t+ t ( H )t y z t qxq x+ x + z x t qyq y+ y + x y t qzq z+ z + y z t ( VxH )x( VxH )x+ x + z x t ( VyH )y( VyH )y+ y + x y t ( VzH )z( VzH )

11、z+ z = x y z ( H )t+ t ( H )t25同除x y z ttHHzHVHVyHVHVxHVHVzqqyqqxqqtttzzzzzyyyyyxxxxxzzzyyyxxx)()()()()()(26222222()()()() C- xyyxxzzpxyxzxyzqV HqV HqV HHxyzxyztTdHdT qVVVTTTHHHHHHHVVVxyzxyzxyzt将代入27222222 0 ()yxzpxyzppxyzVVVxyzTTTTTTTC VVVCxyzxyztTTTTCVVVxyzt代入连续方程28 此式称为傅立叶克希荷夫导热微方程 适用：一切传导、对流的稳定和

12、不稳定传热 222222xyzpTTTTTTTVVVtxyzCxyz29n对固体中的导热中的傅立叶第二定律是v=0时的特例： 固体稳态时的导热是 时的特例： 222222zTyTxTtT0tT0 02222222TzTyTxT即：30n Vx, Vy, Vz, T, p 五个方程，五个未知数 联立求解，可用于各种稳定于非稳定导热222222222222222222xyzxxxxxxxyzxyyyyyyxyzyxTTTTTTTVVVtxyzxyzVVVVVVVpVVVgtxyzxyzxVVVVVVVpVVVgtxyzxyzyVVVt2222220zzzzzzyzzyxzVVVVVpVVgxyzx

13、yzzVVVxyz31n（9-5）方程可写成：n书上191页给出柱面坐标系的对流传热方程的形式。即（9-7）式： 2TVTTt2222211()rzVTTTTTTTVVrtrrzr rrrz329.3 沿平板强制层流的对流换热Heat Transfer with Laminar Forced Convection over a Flat Platen工程中常遇到运动流体与固体壁之间的对流换热情形，因而，定量分析和计算流动流体与固体壁间的对流换热具有工程意义。 如前所述，牛顿换热计算式: q = h (Ts - T ) 为定量计算对流流体与固体壁之间的换热热流密度(通量)的基本表达式。33 又由

14、(9-3)式知n为准确确定对流换热系数h，不仅需要知道T 和Ts，而且还需要知道流动流体内部（至少是固体壁附近的）的温度分布 T (t,x,y,z)。n(9-5)式给出了定量描述流体内部对流传热及温度场的基本微分方程，给出具体对流换热的初始条件和边界条件即给出具体问题的定解方程组（包括流动方程和连续性方程）。n本节首先研究沿一水平大平板强制层流对流换热的情况。0ST -|y(T -T)yh34一、平板边界层流流动对流传热的定解描述方程Governing Equations for Heat Transfer with Laminar Forced Convection over a Flat

15、Plate 本问题为流动流体中的二维稳态对流传热，其传热方程可由(9-5)式直接简化而得到，为： (9-9) 式中： 热扩散系数； 项可忽略。22xyTTTVVxxypc/22xT35说明： 上式(9-9)中，考虑了 Vx和Vy对流传热的影响； 方程右端为扩散项，对于沿平板流动流体与板间的对流传热，以垂直平板方向传热为主，在靠近平板处（V0），以导热为主； 在 x方向上，由于对流，温度分布近于均匀T/x=0 ， 所以在(9-9)式中忽略了2T/x2 项。36 对流传热方程(9-9)式中，含流动速度分量Vx和Vy，所以还需要流动的方程组（N-S方程，和连续方程） 强制流动动量方程（作用力忽略不计

16、）： (1) 连续方程： (2) 方程(9-9)，(1)和(2)构成沿平板强制对流传热的微分方程。还需要边界条件，才能构成相应的问题的定解方程组。 22xxxxyVVVVVxxy0yVxVyx37 边界条件：n传热边界条件 ： (1) (2)n流动边界条件： ,yTT 0,SyTT(1)0,0,0(2),0 xyxyyvvyvvv 38n分析(9-9)式与(1)式可见，这两个方程的形式是相似的（不同之处在于两个方程中的扩散传输项系数 和 ）。 在前面已经给出沿平板流动速度Vx的解的形式及(x)的表达式，显然，流动流体中的温度分布 T (x，y)及T(x) 的解的形式应该与Vx和(x) 的形式是

17、相似的。39n其相似程度由如下无量纲数群（也称相似准则数）决定。该无量纲数群，称为普郎特(Prandtl)数，为： 显然Prandtl数为流体的物性参数，它的物理意义是： Pr=流体的动量扩散能力流体的热量扩散能力 当 Pr=/=1时，方程(9-9)与(1)式等效。T(x)应与v(x)重合。prCP40二、精确解的结果及努塞尔数Nu Exact Solution and Nusselt Numbern利用上述的传热边界条件，求解对流传热方程（9-9）的方法和步骤是类似的，（在方程中已包括了一些近似，但仍称为精确解）其精确解的结果，表示为无量纲温度 = (T - Ts)/(T - Ts)，与曲线

18、坐标 之间的不同Pr数的曲线簇。 见图(9-8)各种Pr下的平板上层流边界层内的无量纲温度分布。 /()y Vx4142n已知温度分布，就可求得对流传热系数h，对于Pr 0.6的情形，根据解的结果，图9-8求得局部换热系数为： (9-10) 引入无量纲数群 Nux表示上式(9-10)，则有： (9-11)130/0.332 PxsryVhTyTx 10.530.332uxrexNP R43n式中： 称为 x 处的局部努塞系数，是描述流体与平板之 间对流传热的重要的无量纲数群，其物理含义为： 流体的边界导热热阻/流体的对流换热热阻 n对于整块平板上的平均对流传热系数，可对 x=0,L区域对hx求积分平均值即可： xuxhxN/1/xuxxhxxNh13010.664LLxrVhh dxPLL44n最后得到的整块平板的平均数 Nu 为： (9-13)n由Nu的定义式知，求 h与求 Nu是等价的（Nu 已知，h 便可知）n由(9-11)式及(9-13)式可知，对于强制对流换热来说，其 Nu数为 Re 和Pr 的系数，即： Nu = f ( Re , Pr)10.530.6