第四章工程传热学ppt课件

1、2022-1-3014-1 对流换热概述对流换热概述4-2 层流流动换热的微分方程组层流流动换热的微分方程组4-3 对流换热过程的相似理论对流换热过程的相似理论4-4 边界层理论边界层理论4-5 紊流流动换热紊流流动换热2022-1-3021 对流换热过程对流换热过程对流换热定义：流体和与之接触的固体壁对流换热定义：流体和与之接触的固体壁面之间的热量传递过程，是宏观的热对流与面之间的热量传递过程，是宏观的热对流与微观的热传导的综合传热过程。微观的热传导的综合传热过程。对流换热与热对流不同，既有热对流，也有对流换热与热对流不同，既有热对流，也有导热；不是基本传热方式导热；不是基本传热方式对流换热

2、实例：对流换热实例：1) 暖气管道暖气管道; 2) 电子器件冷电子器件冷却却2022-1-303对流换热的特点：对流换热的特点：(1) 导热与热对流同时存在的复杂热传递过程导热与热对流同时存在的复杂热传递过程(2) 必须有直接接触流体与壁面和宏观运动；必须有直接接触流体与壁面和宏观运动；也必须有温差也必须有温差特征：以简单的对流换热过程为例，对对流特征：以简单的对流换热过程为例，对对流换热过程的特征进行粗略的分析。换热过程的特征进行粗略的分析。2022-1-304图表示一个简单的对流换热过程。流体以来图表示一个简单的对流换热过程。流体以来流速度流速度u和来流温度和来流温度t流过一个温度为流过一

3、个温度为tw的固体壁面。选取流体沿壁面流动的方向为的固体壁面。选取流体沿壁面流动的方向为x坐标、垂直壁面方向为坐标、垂直壁面方向为y坐标。坐标。y t u tw qw x ttw2022-1-305壁面对流体分子的吸附作用，使得壁面上壁面对流体分子的吸附作用，使得壁面上的流体是处于不滑移的状态此论点对于的流体是处于不滑移的状态此论点对于极为稀薄的流体是不适用的）。极为稀薄的流体是不适用的）。又由于粘性力的作用，使流体速度在垂直于又由于粘性力的作用，使流体速度在垂直于壁面的方向上发生改变。流体速度从壁面上壁面的方向上发生改变。流体速度从壁面上的零速度值逐步变化到来流的速度值。的零速度值逐步变化到

4、来流的速度值。同时，通过固体壁面的热流也会在流体分子同时，通过固体壁面的热流也会在流体分子的作用下向流体扩散的作用下向流体扩散(热传导热传导)，并不断地被，并不断地被流体的流动而带到下游热对流），也导致流体的流动而带到下游热对流），也导致紧靠壁面处的流体温度逐步从壁面温度变化紧靠壁面处的流体温度逐步从壁面温度变化到来流温度。到来流温度。 2022-1-3062 对流换热的分类对流换热的分类自然强制hh对流换热：导热对流换热：导热 + 热对流；壁面热对流；壁面+流动流动 流动起因流动起因自然对流：流体因各部分温度不同而引起的自然对流：流体因各部分温度不同而引起的密度差异所产生的流动密度差异所产生

5、的流动Free convection）强制对流：由外力如：泵、风机、水压头强制对流：由外力如：泵、风机、水压头作用所产生的流动作用所产生的流动Forced convection）2022-1-307 流动状态流动状态层流：整个流场呈一簇互相平行的流线层流：整个流场呈一簇互相平行的流线Laminar flow）湍流：流体质点做复杂无规则的运动湍流：流体质点做复杂无规则的运动（Turbulent flow）紊流流动极为普遍紊流流动极为普遍自然现象：收获季节的麦浪滚滚，旗帜在微自然现象：收获季节的麦浪滚滚，旗帜在微风中轻轻飘扬，以及袅袅炊烟都是由空气的风中轻轻飘扬，以及袅袅炊烟都是由空气的紊流引起的

6、。紊流引起的。 层流湍流hh2022-1-3082022-1-309 流体有无相变流体有无相变单相换热单相换热相变换热：凝结、沸腾、升华、凝固、融化相变换热：凝结、沸腾、升华、凝固、融化 流体运动是否与时间相关流体运动是否与时间相关非稳态对流换热：与时间有关非稳态对流换热：与时间有关稳态对流换热：与时间无关稳态对流换热：与时间无关单相相变hh2022-1-3010 流体与固体壁面的接触方式流体与固体壁面的接触方式 内部流动对流换热：管内或槽内内部流动对流换热：管内或槽内外部流动对流换热：外掠平板、圆管、管束外部流动对流换热：外掠平板、圆管、管束管管内内管管外外热热面面朝朝上上 、朝朝下下202

7、2-1-3011管内沸腾管内沸腾对对流流换换热热有有相相变变无无相相变变强制对流强制对流内部流动内部流动圆管内强制对流换热圆管内强制对流换热其它形状管道的对流换热其它形状管道的对流换热外部流动外部流动外掠单根圆管的对流换热外掠单根圆管的对流换热外掠圆管管束的对流换热外掠圆管管束的对流换热外掠平板的对流换热外掠平板的对流换热外掠其它截面柱体的换热外掠其它截面柱体的换热射流冲击换热射流冲击换热自然对流自然对流大空间自然对流大空间自然对流有限空间自然对流有限空间自然对流混合对流混合对流沸腾换热沸腾换热凝结换热凝结换热大空间沸腾大空间沸腾管内凝结管内凝结管外凝结管外凝结2022-1-30123 对流换

8、热系数与对流换热微分方程对流换热系数与对流换热微分方程当流体与壁当流体与壁面温度相差面温度相差1时、时、每单位壁面面积每单位壁面面积上、单位时间内上、单位时间内所传递的热量所传递的热量.对流换热系数对流换热系数( (表面传热系数表面传热系数) )C)(mW)( 2 ttAhw确定确定h及增强换热的措施是对流换热的核心问题及增强换热的措施是对流换热的核心问题2022-1-3013 对流换热过程微分方程式对流换热过程微分方程式壁面上的流体分子层由于受到固体壁面的吸附壁面上的流体分子层由于受到固体壁面的吸附是处于不滑移的状态，其流速应为零，那么通是处于不滑移的状态，其流速应为零，那么通过它的热流量只

9、能依靠导热的方式传递。过它的热流量只能依靠导热的方式传递。 y t u tw qw x由傅里叶定律由傅里叶定律 0ywytq通过壁面流体层传导的热流量最终是以对流通过壁面流体层传导的热流量最终是以对流换热的方式传递到流体中换热的方式传递到流体中 2022-1-3014cwqq0ywcyttthq0yytth或或对流换热过程微对流换热过程微分方程式分方程式h 取决于流体热导率、温度差和贴壁流体的取决于流体热导率、温度差和贴壁流体的温度梯度温度梯度温度梯度或温度场与流速、流态、流动起因、温度梯度或温度场与流速、流态、流动起因、换热面的几何因素、流体物性均有关。换热面的几何因素、流体物性均有关。速度

10、场和温度场由对流换热微分方程组确定：速度场和温度场由对流换热微分方程组确定：连续性方程、动量方程、能量方程连续性方程、动量方程、能量方程2022-1-3015为便于分析，只限于分析二维对流换热为便于分析，只限于分析二维对流换热假设：假设：a) 流体为不可压缩的牛顿型流体，流体为不可压缩的牛顿型流体，（即：服从牛顿粘性定律的流体；而油漆、（即：服从牛顿粘性定律的流体；而油漆、泥浆等不遵守该定律，称非牛顿型流体）泥浆等不遵守该定律，称非牛顿型流体）b) 所有物性参数（所有物性参数（、cp、）为）为常量常量yu2022-1-30164个未知量：速度个未知量：速度 u、v；温度；温度 t；压力；压力

11、p需要需要4个方程个方程: 连续性方程连续性方程1）; 动量方程动量方程2）;能量方程能量方程1）1 连续性方程连续性方程流体的连续流动遵循质量守恒规律。流体的连续流动遵循质量守恒规律。从流场中从流场中 (x, y) 处取出边长为处取出边长为 dx、dy 的微元的微元体，并设定体，并设定x方向的流体流速为方向的流体流速为u，而，而y方向上方向上的流体流速为的流体流速为v 。 M 为质量流量为质量流量 kg/s2022-1-3017balance mass changedmassoutmassinmassudyyvv vdxdy单位时间内流入微元体的净质量单位时间内流入微元体的净质量 = 微元体

12、内流微元体内流体质量的变化。体质量的变化。 单位时间内、沿单位时间内、沿x轴方轴方向流入微元体的净质量：向流入微元体的净质量：dxxxMM dxdyxu )( dydxxuuudy dxxuu 2022-1-3018单位时间内、沿单位时间内、沿y轴方向流入微元体的净质量：轴方向流入微元体的净质量：dxdyyv )( dyyyMM dxdyyvvvdx 单位时间内微元体内流体质量的变化单位时间内微元体内流体质量的变化:dxdydxdy )(2022-1-3019dxdyxu )( dxdyyv )( dxdy 单位时间：流入微元体的净质量单位时间：流入微元体的净质量 = = 微元体内微元体内流体

13、质量的变化流体质量的变化0yvxu连续性方程：连续性方程：对于二维、稳定、常物性流场对于二维、稳定、常物性流场 ：xu0yv 2022-1-30202 动量微分方程动量微分方程作用力作用力 = 动量变化率动量变化率 F=d(mc)/d动量微分方程式描述流体速度场动量微分方程式描述流体速度场动量守恒动量守恒动量微分方程是纳维埃和斯托克斯分别于动量微分方程是纳维埃和斯托克斯分别于1827和和1845年推导的。年推导的。 Navier-Stokes方程方程N-S方程）方程） 牛顿第二运动定律：作用在微元体上各外牛顿第二运动定律：作用在微元体上各外力的总和等于控制体中流体动量的变化率力的总和等于控制体

14、中流体动量的变化率控制体中流体动量的变化率控制体中流体动量的变化率2022-1-3021从从x x方向进入元体质量流量方向进入元体质量流量在在x x方向上的动量方向上的动量 ：uudy1从从x x方向流出元体的质量流方向流出元体的质量流量在量在x x方向上的动量方向上的动量 dxxuudydxxuu1从从y方向进入元体的质量流量在方向进入元体的质量流量在x方向上的动量为方向上的动量为 ：uvdx1从从y方向流出元体的质量流量在方向流出元体的质量流量在x方向上的动量：方向上的动量： dyyuudxdyyvv1vudxdydxxuu dyyvv 2022-1-3022x方向上的动量改变量方向上的动

15、量改变量 ：1 dxdyyuvxuu化简过程中利用了连续性方程和忽略了高阶化简过程中利用了连续性方程和忽略了高阶小量。小量。 同理，导出同理，导出y方向上的动量改变量方向上的动量改变量 ：1 dxdyyvvxvu作用于微元体上的外力作用于微元体上的外力 作用力：体积力、表面力作用力：体积力、表面力2022-1-3023体积力：重力、离心力、电磁力体积力：重力、离心力、电磁力设定单位体积流体的体积力为设定单位体积流体的体积力为F，相应在，相应在x和和y方向上的分量分别为方向上的分量分别为Fx和和Fy。 在在x方向上作用于微元体的体积力：方向上作用于微元体的体积力：在在y方向上作用于微元体的体积力

16、：方向上作用于微元体的体积力：1 dxdyFx1 dxdyFy表面力表面力:作用于微元体表面上的力。作用于微元体表面上的力。通常用作用于单位表面积上的力来表示，称通常用作用于单位表面积上的力来表示，称之为应力。包括粘性引起的切向应力和法向之为应力。包括粘性引起的切向应力和法向应力、压力等。应力、压力等。法向应力法向应力 中包括了压力中包括了压力 p 和法向粘性应和法向粘性应力力 。2022-1-3024在物理空间中面矢量和力矢量各自有三个相互在物理空间中面矢量和力矢量各自有三个相互独立的分量方向），因而对应组合可构成应独立的分量方向），因而对应组合可构成应力张量的九个分量。于是应力张量可表示为

17、力张量的九个分量。于是应力张量可表示为 333231232221131211ij式中式中 为应力张量，下标为应力张量，下标i表表示作用面的方向，下标示作用面的方向，下标j则表示作用力的方向则表示作用力的方向 3 , 2, 1; 3 , 2, 1,jiij通常将作用力和作用面方向一致的应力分量通常将作用力和作用面方向一致的应力分量称为正应力，而不一致的称为切应力。称为正应力，而不一致的称为切应力。 2022-1-3025对于我们讨论的二对于我们讨论的二维流场应力只剩下维流场应力只剩下四个分量，记为四个分量，记为 yyxxyxx为为x方向上的正应力力与面方向一致）；方向上的正应力力与面方向一致）；

18、 y为为y方向上的正应力力与面方向一致）；方向上的正应力力与面方向一致）； xy为作用于为作用于x表面上的表面上的y方向上的切应力；方向上的切应力； yx为作用于为作用于y表面上的表面上的x方向上的切应力。方向上的切应力。 2022-1-3026作用在作用在x方向上表面力的净值为方向上表面力的净值为 :11dxdyxdxdyyxyx作用在作用在y方向上表面力的净值为方向上表面力的净值为 11dxdyydxdyxyxyxvyuyxxy斯托克斯提出了归纳速斯托克斯提出了归纳速度变形率与应力之间的度变形率与应力之间的关系的黏性定律关系的黏性定律 xupx2yvpy22022-1-3027得出作用在微

19、元体上表面力的净值表达式：得出作用在微元体上表面力的净值表达式： x方向上方向上 12222dxdyyuxuxpy方向上方向上 12222dxdyyvxvyp动量微分方程式动量微分方程式在在x方向上方向上 2222yuxuxpFyuvxuuuxy方向上方向上 2222yvxvypFyvvxvuvy2022-1-3028对于稳态流动：对于稳态流动：0 0vu；只有重力场时：只有重力场时：yyxxgFgF ;3 能量微分方程能量微分方程能量微分方程式描述流体温度场能量微分方程式描述流体温度场能量守恒能量守恒 导入与导出的净热量导入与导出的净热量 + + 热对流传递的净热对流传递的净热量热量 + +

20、内热源发热量内热源发热量 = = 总能量的增量总能量的增量 + + 对外作膨胀功对外作膨胀功 2022-1-3029Q = E + W内热源内热源对流对流导热导热QQQQ （动能）（动能）热力学能热力学能K UUE W 体积力体积力( (重力重力) )作的功作的功表面力作的功表面力作的功（1压力作的功：压力作的功： a) 变形功；变形功；b) 推动功推动功（2表面应力作的功：表面应力作的功：a) 动能；动能；b) 2022-1-3030Q = E + W内热源内热源对流对流导热导热QQQQ （动能）（动能）热力学能热力学能K UUE W 体积力体积力( (重力重力) )作的功作的功表面力作的功

21、表面力作的功一般可忽略一般可忽略（1压力作的功：压力作的功：a) 变形功；变形功；b) 推动功推动功 （2表面应力法向表面应力法向+切向作的功：切向作的功：a) 动能；动能；b) 耗散热耗散热 假设：（假设：（1流体的热物性均为常量流体的热物性均为常量变形功变形功=0Q内热源内热源=0（2流体不可压缩流体不可压缩 （3一般工程问题流速低一般工程问题流速低 （4无化学反应等内热源无化学反应等内热源2022-1-3031Q Q导热导热 + Q + Q对流对流 = = U U热力学能热力学能 + + 推动功推动功 = = H H耗散热（耗散热（ ）：由表面粘性应力产生的摩）：由表面粘性应力产生的摩擦

22、力而转变成的热量。擦力而转变成的热量。对于二维不可压缩常物性流体流场而言，微元对于二维不可压缩常物性流体流场而言，微元体的能量平衡关系式为：体的能量平衡关系式为： HQQQ 321Q1Q1为以传导方式进入元体的净的热流量；为以传导方式进入元体的净的热流量；Q2Q2为以对流方式进入元体的净的热流量；为以对流方式进入元体的净的热流量； Q3Q3为元体粘性耗散功率变成的热流量；为元体粘性耗散功率变成的热流量；HH为元体的焓随时间的变化率。为元体的焓随时间的变化率。 2022-1-3032以传导方式进入元体的净热流量以传导方式进入元体的净热流量 单位时间沿单位时间沿x轴方向导入与导出微元体净热量：轴方

23、向导入与导出微元体净热量：dyQQdxxx)(单位时间沿单位时间沿y轴方向导入与导出微元体净热量：轴方向导入与导出微元体净热量：dxdyytdxQQdyyy22)(dydxytdxdyxtQ2222 导热导热dydxxQQQxxxdxdyxt222022-1-3033以对流方式进入元体的净热流量以对流方式进入元体的净热流量 单位时间沿单位时间沿 x 方向热对流传递到微元体净热量方向热对流传递到微元体净热量dxxxQQ单位时间沿单位时间沿y 方向热方向热对流传递到微元体对流传递到微元体的净热量：的净热量：dydxyvtcdyyQpy)(dxxQQQxxxdxxQxdxdyxutcp)(2022-

24、1-3034dydxyvtcdxdyxutcQpp)( )(对流元体粘性耗散功率变成的热流量元体粘性耗散功率变成的热流量 13dxdyQ122223dxdyxvyuyvxuQ单位时间内、微元体内焓的增量：单位时间内、微元体内焓的增量：dxdytctdxdyctmcppp2022-1-3035能量微分方程能量微分方程 2222ytxtytvxtutcp当流体不流动时，流体流速为零，热对流项和当流体不流动时，流体流速为零，热对流项和黏性耗散项也为零，能量微分方程式便退化为黏性耗散项也为零，能量微分方程式便退化为导热微分方程式，导热微分方程式， 2222ytxttcp所以，固体中的热传导过程是介质中

25、传热过程所以，固体中的热传导过程是介质中传热过程的一个特例。的一个特例。 2022-1-30364层流流动对流换热微分方程组层流流动对流换热微分方程组（常物性、无内热源、二维、不可压缩牛顿（常物性、无内热源、二维、不可压缩牛顿流体）流体）2222ytxtytvxtutcp)()()22222222yvxvypFyvvxvuvyuxuxpFyuvxuuuyx（xu0yv 4个方程，个方程，4个未知量个未知量 ， 可求速度场和温度场可求速度场和温度场2022-1-3037再引入换热微分方程再引入换热微分方程 (n为壁面为壁面的法线方向坐标的法线方向坐标)，最后可以求出流体与固体，最后可以求出流体与

26、固体壁面之间的对流换热系数，从而解决给定的壁面之间的对流换热系数，从而解决给定的对流换热问题。对流换热问题。 0nntth5 求解对流换热问题的途径求解对流换热问题的途径 分析求解。分析求解。实验研究。实验研究。数值求解。数值求解。6 对流换热单值性条件对流换热单值性条件2022-1-3038单值性条件：能单值反映对流换热过程特点的单值性条件：能单值反映对流换热过程特点的条件条件完整数学描述：对流换热微分方程组完整数学描述：对流换热微分方程组 + 单值性单值性条件条件单值性条件包括：几何、物理、时间、边境单值性条件包括：几何、物理、时间、边境 几何条件：说明对流换热过程中的几何形几何条件：说明

27、对流换热过程中的几何形状和大小，平板、圆管；竖直圆管、水平圆管；状和大小，平板、圆管；竖直圆管、水平圆管；长度、直径等长度、直径等物理条件：说明对流换热过程物理特征，如：物理条件：说明对流换热过程物理特征，如：物性参数物性参数 、 、c c 和和 的数值，是否随的数值，是否随温度温度 和压力变化；有无内热源、大小和分布和压力变化；有无内热源、大小和分布2022-1-3039时间条件：说明在时间上对流换热过程的特时间条件：说明在时间上对流换热过程的特点，稳态对流换热过程不需要时间条件点，稳态对流换热过程不需要时间条件 与与时间无关时间无关边界条件：说明对流换热过程的边界特点，边界条件：说明对流换

28、热过程的边界特点，边界条件可分为二类：第一类、第二类边界条边界条件可分为二类：第一类、第二类边界条件件（1第一类边界条件：已知任一瞬间对流换第一类边界条件：已知任一瞬间对流换热过程边界上的温度值热过程边界上的温度值（2第二类边界条件：已知任一瞬间对流换第二类边界条件：已知任一瞬间对流换热过程边界上的热流密度值热过程边界上的热流密度值2022-1-3040由于对流换热是复杂的热量交换过程，所涉由于对流换热是复杂的热量交换过程，所涉及的变量参数比较多，常常给分析求解和实及的变量参数比较多，常常给分析求解和实验研究带来困难。验研究带来困难。人们常采用相似原则对换热过程的参数进行人们常采用相似原则对换

29、热过程的参数进行归类处理，将物性量，几何量和过程量按物归类处理，将物性量，几何量和过程量按物理过程的特征组合成无量纲的数，这些数常理过程的特征组合成无量纲的数，这些数常称为准则称为准则2022-1-3041首先选取对流换热过程中有关变量的特征值，首先选取对流换热过程中有关变量的特征值，将所有变量无量纲化，进而导出无量纲形式将所有变量无量纲化，进而导出无量纲形式的对流换热微分方程组。的对流换热微分方程组。出现在无量纲方程组中的系数项就是我们所出现在无量纲方程组中的系数项就是我们所需要无量纲数或称：无因次数），也就是需要无量纲数或称：无因次数），也就是无量纲准则，它们是变量特征值和物性量的无量纲准

30、则，它们是变量特征值和物性量的某种组合。某种组合。流场中的任一无量纲变量均可表示为其余无流场中的任一无量纲变量均可表示为其余无量纲变量和无量纲准则的函数形式。量纲变量和无量纲准则的函数形式。 2022-1-3042y u tPin Pout0 L x以流体流过平板以流体流过平板的对流换热问题的对流换热问题为例来进行换热为例来进行换热过程的相似分析。过程的相似分析。 流体平行流过平板的对流换热过程如图所示，流体平行流过平板的对流换热过程如图所示，来流速度为来流速度为u，来流温度，来流温度t，平板长度，平板长度L, 平平板温度板温度tW ，流体流过平板的压力降为，流体流过平板的压力降为 p。如果为

31、二维、稳态、流体物性为常数，且忽如果为二维、稳态、流体物性为常数，且忽略黏性耗散项和体积力项，按图中所示的坐略黏性耗散项和体积力项，按图中所示的坐标流场的支配方程为标流场的支配方程为 2022-1-3043y u tPin Pout0 L x2222ytxtytvxtutcp)()()22222222yvxvypFyvvxvuvyuxuxpFyuvxuuuyx（xu0yv 0nntth2022-1-3044y u tPin Pout0 L x今选取板长今选取板长L，来，来流流速流流速u，温度，温度差差t=tw-t 和压和压力降力降 p=pin-pout为变量的特征值为变量的特征值;/;/uvv

32、uuu用这些无量纲变量去取代方程组中的相应用这些无量纲变量去取代方程组中的相应变量，可得出无量纲变量组成的方程组。变量，可得出无量纲变量组成的方程组。 ;/;/LyyLxx;/ppp)/()(ttttw2022-1-3045y u tPin Pout0 L x；0yvxuLxuuLu0)(yvxuLu；2222yuxuxpyuvxuuLyuuLuuvuLxuuLuuuu；2222yuxuLuxpLpyuvxuuLu；0LyuvLuLxppLp；）（）（2222LuuLuLxuuLu2022-1-3046；2222yvxvLuypLpyvvxvuLu；2222yxLtyvxuLtucp0yyNu

33、惯性力惯性力粘性力粘性力热对流热量热对流热量热传导热量热传导热量2022-1-3047对方程整理，可以得到无量纲化的方程组。对方程整理，可以得到无量纲化的方程组。 0yvxu2222Re1yuxuxpEuyuvxuu22222yuxuLuxpupyuvxuu2222Re1yvxvypEuyvvxvu2222PrRe1yxyvxu2222yxLuayvxu0yyNu2022-1-30482 无量纲准则的表达式和物理意义无量纲准则的表达式和物理意义)/(2upEu定义为欧拉数定义为欧拉数EulerEuler），它反映了流场压力），它反映了流场压力降与其动压头之间的相对关系，体现了在流降与其动压头之

34、间的相对关系，体现了在流动过程中动量损失率的相对大小。动过程中动量损失率的相对大小。 2222Re1yuxuxpEuyuvxuu22222yuxuLuxpupyuvxuu2022-1-3049LuLuRe称为雷诺数，表征了给定流场的惯性力与称为雷诺数，表征了给定流场的惯性力与其黏性力的对比关系，也就是反映了这两其黏性力的对比关系，也就是反映了这两种力的相对大小。种力的相对大小。利用雷诺数可以判别一个给定流场的稳定性，利用雷诺数可以判别一个给定流场的稳定性，随着惯性力的增大和黏性力的相对减小，雷随着惯性力的增大和黏性力的相对减小，雷诺数就会增大，而大到一定程度流场就会失诺数就会增大，而大到一定程

35、度流场就会失去稳定，而使流动从层流变为紊流。去稳定，而使流动从层流变为紊流。2022-1-30502222PrRe1yxyvxuRe PrpPec u Lu L a称为贝克莱称为贝克莱Peclet准绳，记为准绳，记为Pe，它反映，它反映了给定流场的热对流能力与其热传导能力的对了给定流场的热对流能力与其热传导能力的对比关系。它在能量微分方程中的作用相当于雷比关系。它在能量微分方程中的作用相当于雷诺数在动量微分方程中的作用。诺数在动量微分方程中的作用。 其中其中 :aPr称为普朗特称为普朗特Prandtl数，它反映了流体的数，它反映了流体的动量扩散能力与其热扩散能力的对比关系。动量扩散能力与其热扩

36、散能力的对比关系。 2022-1-3051hLNu 努塞尔努塞尔Nusselt准绳，它反映了给定流场的准绳，它反映了给定流场的换热能力与其导热能力的对比关系。这是一个换热能力与其导热能力的对比关系。这是一个在对流换热计算中必须要加以确定的准则。在对流换热计算中必须要加以确定的准则。 uchNuStpPrRe斯坦顿斯坦顿Stanton数，修正的努塞尔数，流数，修正的努塞尔数，流体实际的换热热流密度与可传递之最大热流密体实际的换热热流密度与可传递之最大热流密度之比。度之比。2022-1-3052努谢尔特准则与非稳态导热分努谢尔特准则与非稳态导热分析中的毕欧数形式上是相似的。析中的毕欧数形式上是相似

37、的。但是，但是，Nu中的中的Lf为流场的特为流场的特征尺寸，征尺寸，f为流体的导热系数；为流体的导热系数；流体侧流体侧固体侧固体侧LsLfffssNuBi而而Bi中的中的Ls为固体系统的特征尺寸，为固体系统的特征尺寸，s为固体为固体的导热系数。的导热系数。它们虽然都表示边界上的无量纲温度梯度，它们虽然都表示边界上的无量纲温度梯度，但一个在流体侧一个在固体侧。但一个在流体侧一个在固体侧。2022-1-3053在运用相似理论时，应该注意：只有属于同在运用相似理论时，应该注意：只有属于同一类型的物理现象才有相似的可能性，也才一类型的物理现象才有相似的可能性，也才能谈相似问题。能谈相似问题。所谓同类现

38、象，就是指用相同形式和内容的所谓同类现象，就是指用相同形式和内容的微分方程控制方程微分方程控制方程+单值性条件方程所描单值性条件方程所描述的现象。述的现象。电场与温度场：电场与温度场： 微分方程相同；内容不同微分方程相同；内容不同强制对流换热与自然对流换热：微分方程的强制对流换热与自然对流换热：微分方程的形式和内容都有差异形式和内容都有差异外掠平板和外掠圆管：控制方程相同；单值外掠平板和外掠圆管：控制方程相同；单值性条件不同性条件不同2022-1-3054判断两个现象是否相似的条件：凡同类现象、判断两个现象是否相似的条件：凡同类现象、单值性条件相似、同名已定特征数相等，那单值性条件相似、同名已

39、定特征数相等，那么现象必定相似。据此，如果两个现象彼此么现象必定相似。据此，如果两个现象彼此相似，它们的同名准则数必然相等。相似，它们的同名准则数必然相等。2022-1-30550yvxu2222Re1yuxuxpEuyuvxuu2222Re1yvxvypEuyvvxvu2222PrRe1yxyvxu0yyNu, , ,Re, , , ,Re,yxpEufvyxpEufuvuRe, ,epfEuyxEufp, , , Pr,Re,yxvuf, Pr,Re,yxfxyNuxfyPr,Re,0PrRe,fNux2022-1-3056从上式不难看出，在计算几何形状相似的流从上式不难看出，在计算几何形

40、状相似的流动换热问题时，如果只是求取其平均的换热动换热问题时，如果只是求取其平均的换热性能，就可以归结为确定几个准则之间的某性能，就可以归结为确定几个准则之间的某种函数关系，最后得出平均的表面传热系数种函数关系，最后得出平均的表面传热系数和总体的换热热流量。和总体的换热热流量。由于无量纲准则是由过程量、几何量和物性由于无量纲准则是由过程量、几何量和物性量组成的，从而使实验研究的变量数目显著量组成的，从而使实验研究的变量数目显著减少，这对减少实验工作量和实验数据处理减少，这对减少实验工作量和实验数据处理时间是至关重要的。时间是至关重要的。2022-1-3057对流动换热微分方程组进行无量纲化时，

41、选对流动换热微分方程组进行无量纲化时，选定了对应变量的特征值，然后进行无量纲化定了对应变量的特征值，然后进行无量纲化的工作，这些特征参数是流场的代表性的数的工作，这些特征参数是流场的代表性的数值，分别表征了流场的几何特征、流动特征值，分别表征了流场的几何特征、流动特征和换热特征。和换热特征。 特征尺寸，它反映了流场的几何特征，对于特征尺寸，它反映了流场的几何特征，对于不同的流场特征尺寸的选择是不同的。如，不同的流场特征尺寸的选择是不同的。如，对流体平行流过平板选择沿流动方向上的长对流体平行流过平板选择沿流动方向上的长度尺寸；管内流体流动选择垂直于流动方向度尺寸；管内流体流动选择垂直于流动方向的

42、管内直径；对于流体绕流圆柱体流动选择的管内直径；对于流体绕流圆柱体流动选择流动方向上的圆柱体外直径。流动方向上的圆柱体外直径。 2022-1-3058特征流速，它反映了流体流场的流动特征。特征流速，它反映了流体流场的流动特征。不同的流场其流动特征不同，所选择的特征不同的流场其流动特征不同，所选择的特征流速是不同的。流速是不同的。如，流体流过平板，来流速度被选择为特征如，流体流过平板，来流速度被选择为特征尺寸；尺寸；流体管内流动，管子截面上的平均流速可作流体管内流动，管子截面上的平均流速可作为特征流速；为特征流速；流体绕流圆柱体流动，来流速度可选择为特流体绕流圆柱体流动，来流速度可选择为特征流速

43、。征流速。2022-1-3059定性温度，无量纲准则中的物性量是温度的函定性温度，无量纲准则中的物性量是温度的函数，确定物性量数值的温度称为定性温度。对数，确定物性量数值的温度称为定性温度。对于不同的流场定性温度的选择是不同的。于不同的流场定性温度的选择是不同的。外部流动常选择来流流体温度和固体壁面温度外部流动常选择来流流体温度和固体壁面温度的算术平均值，称为膜温度；的算术平均值，称为膜温度；内部流动常选择管内流体进出口温度的平均值内部流动常选择管内流体进出口温度的平均值算术平均值或对数平均值），当然也有例外。算术平均值或对数平均值），当然也有例外。2)(twfmtt2022-1-3060由于

44、对流换热问题的复杂性，实验研究是解决由于对流换热问题的复杂性，实验研究是解决换热问题的主要方法。换热问题的主要方法。在工程上大量使用的对流换热准则关系式都是在工程上大量使用的对流换热准则关系式都是通过实验获得的。通过实验获得的。 我们从无量纲微分方程组推出了一般化的准则我们从无量纲微分方程组推出了一般化的准则关系式关系式 。但这是一个原则性的式。但这是一个原则性的式子，要得到某种类型的对流换热问题在给定范子，要得到某种类型的对流换热问题在给定范围内的具体的准则关系式，在多数情况下还必围内的具体的准则关系式，在多数情况下还必须通过实验的办法来确定。须通过实验的办法来确定。 PrRe,fNux20

45、22-1-3061twqLBt u图中给出了平板在风图中给出了平板在风洞中进行换热实验的洞中进行换热实验的示意图。示意图。 为了得出该换热问题的准则关系式，必须测量为了得出该换热问题的准则关系式，必须测量的物理量有：流体来流速度的物理量有：流体来流速度u,来流温度来流温度t,平平板表面温度板表面温度tw,平板的长度平板的长度L和宽度和宽度B，以及平，以及平板的加热量板的加热量Q通过测量电加热器的电流通过测量电加热器的电流I和电和电压压V而得出）。而得出）。 可由可由 得到得到 wQI Vh ttLBwhIVttLB必须在不同的工况下获得不同的换热系数值必须在不同的工况下获得不同的换热系数值 。

46、2022-1-3062123123NNhuhuhuhu123111222333ReReReReNNNNNuh Lu LNuh LuLNuh LuLNuh LuL如果认为准则关系式有如果认为准则关系式有 这样的形式。这是一种先验的处理办法，但这样的形式。这是一种先验的处理办法，但是，这给拟合准则关系式带来较大的方便。是，这给拟合准则关系式带来较大的方便。最小二乘法是常用的线性拟合方法最小二乘法是常用的线性拟合方法 。采用几何作图的方法亦可以求解采用几何作图的方法亦可以求解 。 mnmccNuRePrRe12022-1-3063n=tglogNuLogc1 logRe对于几何结构比较复杂对于几何结

47、构比较复杂的对流换热过程，特征的对流换热过程，特征尺寸无法从已知的几何尺寸无法从已知的几何尺度中选取，通常的做尺度中选取，通常的做法是采用当量尺寸。如法是采用当量尺寸。如异型管槽内的流动换热，异型管槽内的流动换热，其当量直径定义为其当量直径定义为Pfde4Pf式中式中f为流体流通面积；为流体流通面积；P为流体的润湿周边。为流体的润湿周边。 2022-1-30642022-1-30652022-1-3066边界层的概念是边界层的概念是1904年德国科学家普朗特提年德国科学家普朗特提出的。出的。1 边界层定义边界层定义速度边界层速度边界层(a) 定义定义流体流过固体壁面时，由于壁面层流体分子的流体

48、流过固体壁面时，由于壁面层流体分子的不滑移特性，在流体黏性力的作用下，近壁流不滑移特性，在流体黏性力的作用下，近壁流体流速在垂直于壁面的方向上会从壁面处的零体流速在垂直于壁面的方向上会从壁面处的零速度逐步变化到来流速度。速度逐步变化到来流速度。 2022-1-3067垂直于壁面的方向上流体流速发生显著变化垂直于壁面的方向上流体流速发生显著变化的流体薄层定义为速度边界层。的流体薄层定义为速度边界层。 普朗特通过观察发现，对于低黏度的流体，普朗特通过观察发现，对于低黏度的流体，如水和空气等，在以较大的流速流过固体壁如水和空气等，在以较大的流速流过固体壁面时，在壁面上流体速度发生显著变化的流面时，在

49、壁面上流体速度发生显著变化的流体层是非常薄的。体层是非常薄的。 2022-1-3068流体流过固体壁面的流体流过固体壁面的流场就人为地分成两流场就人为地分成两个不同的区域。个不同的区域。tw t u t 0 x其一是边界层流动区，这里流体的黏性力与其一是边界层流动区，这里流体的黏性力与流体的惯性力共同作用，引起流体速度发生流体的惯性力共同作用，引起流体速度发生显著变化；显著变化；其二是势流区，这里流体黏性力的作用非常其二是势流区，这里流体黏性力的作用非常微弱，可视为无黏性的理想流体流动，也就微弱，可视为无黏性的理想流体流动，也就是势流流动。是势流流动。2022-1-3069(b)边界层的厚度边

50、界层的厚度当速度变化达到当速度变化达到 时的空间位置为时的空间位置为速度边界层的外边缘，那么从这一点到壁面速度边界层的外边缘，那么从这一点到壁面的距离就是边界层的厚度的距离就是边界层的厚度 99. 0uu x 小：空气外掠平板，小：空气外掠平板，u u =10m/s=10m/s：mm5 . 2 ;mm8 . 1200100mmxmmx热温度边界层热温度边界层(a) 定义定义当流体流过平板而平板的温度当流体流过平板而平板的温度tw与来流流体与来流流体的温度的温度t不相等时，在壁面上方也能形成温不相等时，在壁面上方也能形成温度发生显著变化的薄层，常称为热边界层。度发生显著变化的薄层，常称为热边界层

51、。 2022-1-3070(b)热边界层厚度热边界层厚度 当壁面与流体之间的温差达到壁面与来流流体当壁面与流体之间的温差达到壁面与来流流体之间的温差的之间的温差的0.99倍时，即倍时，即 ，此位置就是边界层的外边缘，而该点到壁面之此位置就是边界层的外边缘，而该点到壁面之间的距离则是热边界层的厚度间的距离则是热边界层的厚度,记为记为 99. 0)/()(ttttww xt层流：温度呈抛物层流：温度呈抛物线分布线分布湍流：温度呈幂函湍流：温度呈幂函数分布数分布湍流边界层贴壁处湍流边界层贴壁处温度梯度明显大温度梯度明显大湍流换热比层流换热强！湍流换热比层流换热强！2022-1-3071引入边界层概念

52、可使换热微分方程组得以简化引入边界层概念可使换热微分方程组得以简化数量级分析数量级分析order of magnitude ：比较方程中：比较方程中各量或各项的量级的相对大小；保留量级较各量或各项的量级的相对大小；保留量级较大的量或项；舍去那些量级小的项，方程大大的量或项；舍去那些量级小的项，方程大大简化大简化无量纲形式的微分方程组对于流体平行流过平无量纲形式的微分方程组对于流体平行流过平板形成的边界层流动换热问题也是同样适用的。板形成的边界层流动换热问题也是同样适用的。 2022-1-30720yvxu2222Re1yuxuxpEuyuvxuu2222Re1yvxvypEuyvvxvu222

53、2PrRe1yxyvxu0yyNu2022-1-30735个基本量的数量级：个基本量的数量级：主流速度：主流速度：);1 (0u温度：温度：);1 (0t壁面特征长度：壁面特征长度：);1 (0L边界层厚度：边界层厚度：)(0 );(0tx与与L相当，即：相当，即：);1 (0 lx)(0 0yy0(1)、0()表示数量级为表示数量级为1和和 ，1 。 “” 相当于相当于2022-1-3074u沿边界层厚度由沿边界层厚度由0到到u:) 1 (0uu主流方向上的无量纲速度主流方向上的无量纲速度 的数量级为的数量级为1 uuu0yvxu2222Re1yuxuxpEuyuvxuu2222Re1yvx

54、vypEuyvvxvu由连续性方程由连续性方程 :0yvxu可以得出可以得出v的数量级为的数量级为 2022-1-3075 0yvxu 11 2222Re1yuxuxpEuyuvxuu 1 11 1 )(22xuxxu)(22yuyyu221 11 211 1 x方向上的动量方程变为：方向上的动量方程变为： 22Re1yuxpEuyuvxuu2022-1-30762222Re1yvxvypEuyvvxvu2222PrRe1yxyvxu) 1 ( 1 1222）（221 11 1 11 1 2t2022-1-3077这就使得动量方程和能量方程变成了抛物型这就使得动量方程和能量方程变成了抛物型的非

55、线性微分方程；的非线性微分方程；0yvxu22Re1yuxpEuyuvxuu22PrRe1yyvxu0yyNu微分方程组经过在边界层中简化后，由于动微分方程组经过在边界层中简化后，由于动量方程和能量方程分别略去了主流方向上的量方程和能量方程分别略去了主流方向上的动量扩散项和热量扩散项，从而构成上游影动量扩散项和热量扩散项，从而构成上游影响下游而下游不影响上游的物理特征。响下游而下游不影响上游的物理特征。 2022-1-3078由于动量方程由两个变成为一个，而且由于动量方程由两个变成为一个，而且项可在边界层的外边缘上利用伯努利方程求项可在边界层的外边缘上利用伯努利方程求解，于是方程组在给定的边值

56、条件下可以进解，于是方程组在给定的边值条件下可以进行分析求解，所得结果为边界层的精确解。行分析求解，所得结果为边界层的精确解。 对于外掠平板的层流流动，主流场速度是均对于外掠平板的层流流动，主流场速度是均速速u u ，温度是均温，温度是均温t t ；并假定平板为；并假定平板为恒温恒温twtw。 xp wttvuy, 0:0ttuuy,:3/12/1332. 0axuxhx3/12/1PrRe332. 0 xxNu 留意：层流留意：层流2022-1-3079 比较边界层无量纲的动量方程和能量方程：比较边界层无量纲的动量方程和能量方程：22Re1yuxpEuyuvxuu22PrRe1yyvxu在忽

57、略动量方程压力项后，温度边界层的厚在忽略动量方程压力项后，温度边界层的厚度与速度边界层的厚度的相对大小则取决于度与速度边界层的厚度的相对大小则取决于普朗特数的大小。普朗特数的大小。当当Pr=1时，动量方程与能量方程完全相同。时，动量方程与能量方程完全相同。即速度分布的解与温度分布完全相同，此时即速度分布的解与温度分布完全相同，此时速度边界层厚度等于温度边界层厚度。速度边界层厚度等于温度边界层厚度。 2022-1-3080当当Pr1时，时，Pr= ，a，粘性扩散，粘性扩散 热量扩热量扩散，速度边界层厚度散，速度边界层厚度温度边界层厚度。温度边界层厚度。当当Pr1时，时，Pr=/a，a，粘性扩散，

58、粘性扩散 热量扩热量扩散，速度边界层厚度散，速度边界层厚度温度边界层厚度。温度边界层厚度。也可从公式得出也可从公式得出 31Pr026. 11xxtTuTx0ttux0tt(a)Pr1a2022-1-3081边界层积分方程组边界层积分方程组1921年，冯年，冯卡门提出了边界层动量积分方程。卡门提出了边界层动量积分方程。1936年，克鲁齐林求解了边界层能量积分方程。年，克鲁齐林求解了边界层能量积分方程。所得的结果称为边界层问题的近似解所得的结果称为边界层问题的近似解 。边界层积分方程一般可由两种方法获得：其一边界层积分方程一般可由两种方法获得：其一是将动量守恒定律和能量守恒定律应用于控制是将动量

59、守恒定律和能量守恒定律应用于控制体；其二是对边界层微分方程直接进行积分。体；其二是对边界层微分方程直接进行积分。2022-1-3082对一固定对一固定x，将能量方程从，将能量方程从y = 0到到y = 积分积分得：得： 采用对微分方程积分得到积分方程采用对微分方程积分得到积分方程能量方程为：能量方程为：22ytaytvxtu00220dyytadyytvdyxtu2022-1-308300220dyytadyytvdyxtu0, 0,22ytytxtttyt故时由分部积分：由分部积分：tttttdyyvttvdyyvtvtdyytv0000（b）tttdyytadyytvdyxtu002202

60、022-1-3084将将v转化为转化为u，利用，利用0yvxuttdyxutdyyvt00tttdyxudyyvv00tttdyxutdyxutdyytv0002022-1-3085式中的扩散项为：式中的扩散项为：代入代入(b)式得：式得： 上式左边可进一步简化为：上式左边可进一步简化为： 00022)(yytaytadyytatttttyytadyxutdyxutdyxtu0000)(ttttdyututdxddyututxdyxutdyxutxtu0000)(最后能量积分方程为：最后能量积分方程为：00yytadyttudxdt2022-1-3086边界层积分方程组求解示例边界层积分方程组