第五章 对 流 换 热.doc

第五章 对 流 换 热本章内容要求：1 、重点内容： 对流换热及其影响因素； 牛顿冷却公式； 用分析方法求解对流换热问题的实质 边界层概念及其应用 相似原理 无相变换热的表面传热系数及换热量的计算 2 、掌握内容：对流换热及其影响因素；用分析方法求解对流换热问题的实质 3 、讲述基本的内容： 对流换热概述； 对流换热的数学描写； 对流换热的边界层微分方程组； 边界层积分方程组的求解及比拟理论； 相似原理及量纲分析； 相似原理的应用； 内部流动强制对流换热实验关联式； 外部流动强制对流换热实验关联式； 自然对流换热实验关联式在绪论中已经指出， 对流换热是发生在流体和与之接触的固体壁面之间的热量传递过程， 是发生在流体中的热量传递过程的特例。由于流体系统中流体的运动，热量将主要以热传导和热对流的方式进行，这必然使热量传递过程比单纯的导热过程要复杂得多。 本章将在对换热过程进行一般性讨论的基础上，将质量守恒、动量守恒和能量守恒的基本定律应用于流体系统，导出支配流体速度场和温度场的场方程对流换热微分方程组。 由于该方程组的复杂性，除少数简单的对流换热问题可以通过分析求解微分方程而得出相应的速度分布和温度分布之外，大多数对流换热问题的分析求解是十分困难的。因此，在对流换热的研究中常常采用实验研究的方法来解决复杂的对流换热问题。在这一章，我们将 通过方程的无量纲化和实验研究方法的介绍而得到常用的准则及准则关系式。 讨论的重点放在工程上常用的管内流动、平行流过平板以及绕流圆管的受迫对流换热，大空间和受限空间的自然对流换热，以及蒸汽凝结与液体沸腾换热。 5-1 对流换热概述本节要求：1。对流换热的概念：流体固体壁面；2对流换热中，导热核对流通式汽作用；3对流换热的影响因素：，h过程量；4对流换热系数如何确定：1 对流换热过程 对流换热是发生在流体和与之接触的固体壁面之间的热量传递过程 ，（ 直接接触是与辐射换热的区别），是宏观的热对流与微观的热传导的综合传热过程。由于涉及流体的运动使热量的传递过程变得较为复杂，分析处理较为困难。因此，在对流换热过程的研究和应用上，实验和数值分析的处理方法是常常采用的。下面我们以简单的对流换热过程为例，对对流换热过程的特征进行粗略的分析。 图 5 1 表示一个简单的对流换热过程。表示流体以来流速度 u 和来流温度 t 流过一个温度为 tw 的固体壁面。这里选取流体沿壁面流动的方向为 x 坐标、垂直壁面方向为 y 坐标。 由于固体壁面对流体分子的吸附作用 ，使得壁面上的流体是处于不流动或不滑移的状态（此论点对于极为稀薄的流体是不适用的）。又由于流体分子相互之间的穿插扩散和 ( 或 ) 相互之间的吸引造成流体之间的相互牵制。这种相互的牵制作用就是流体的黏性力，在其作用下会使流体的速度在垂直于壁面的方向上发生改变 。由于流体的分子在固体壁面上被吸附而处于不流动的状态，因而使流体速度从壁面上的零速度值逐步变化到来流的速度值。同时，通过固体壁面的热流也会在流体分子的作用下向流体扩散 ( 热传导 ) ，并不断地被流体的流动而带到下游（热对流），因而也导致紧靠壁面处的流体温度逐步从壁面温度变化到来流温度。这里，我们把流体在壁面附近的速度和温度分布也示意性地表示在图5 1中。 2 对流换热过程的分类 由于对流换热是发生在流体和固体界面上的热交换过程，流体的流动和固体壁面的几何形状以及相互接触的方式都会不同程度影响对流热交换的效果，由此也构成了许许多多复杂的对流换热过程。因此，为了研究问题的条理性和系统性，以及更便于把握对流换热过程的实质，我们按不同的方式将对流换热过程进行分类。然后再分门别类地进行分析处理。 在传热学中对流换热过程的习惯性分类方式是 : 1）按 流体运动的起因 可分为 自然对流换热 和 受迫对流换热； 2）按 流体与固体壁面的接触方式 可分为 内部流动换热 和 外部流动换热 ； 3）按 流体的运动状态 可分为 层流流动换热 和 湍流流动换热 ； 4）按流体在换热中 是否发生相变 或存在多相的情况可分为 单相流体对流换热 和 多相流体对流换热 。 湍流流动极为普遍，从自然现象看，收获季节的麦浪滚滚，旗帜在微风中轻轻飘扬，都是由空气的湍流引起的。湍流的运动服从某种统计规律，而不是杂乱无章。香烟的烟在静止的空气中上升，可以看到从层流到湍流的转化。湍流会消耗能量（同摩擦力消耗能量一样），没有湍流的世界是不可想象的。如果没有湍流，把酱油倒进汤里，花半小时酱油才能和汤混合，用汤匙一搅，依靠湍流几秒钟它们就混合在一起了。如果没有湍流的掺混，烟囱浓烟中的有害物质将长期积聚，危害人类环境。 对于实际的对流换热过程的，按照上述的分类，总是可以将其归入相应的类型之中。例如，在外力推动下流体的管内流动换热是属于受迫内部流动换热，可以为层流亦可为湍流，也可以有相变发生，使之从单相流动变为多相流动；再如，竖直的热平板在空气中冷却过程是属于外部自然对流换热（或称大空间自然对流换热），可以为层流亦可为湍流，在空气中冷却不可能有相变，应为单相流体换热；但是如果是在饱和水中则会发生沸腾换热，这就是带有相变的多相换热过程。 在本章中，我们将按照上述分类对一些典型的对流换热过程进行分析。具体步骤为，先讨论单相流体受迫对流换热，其中分层流和紊流、管内流动和掠过平板或管子的外部流动，之后讨论大空间自然对流换热，最后介绍有流体发生相变的凝结和沸腾换热。 3 换热系数和换热微分方程式 在绪论中提到对流换热的热流密度可以按照牛顿冷却公式来计算，即，式中， h 为对流换热系数（亦称表面传热系数，记为 h ），其单位是 W/(m2)。采用这样的书写形式是为了使热流的方向与流体温度的降落方向一致。如果 热流方向从固体壁面指向流体，如果 则相反。仔细分析一下这个公式，就不难看出该式只不过是定义了一个对流换热系数而已，并不能直接去解决对流换热问题。但是，利用这个定义的直接好处是，把研究复杂对流换热问题集中到研究和确定对流换热系数上，使复杂问题从形式上得到简化；同时，由于对流换热系数是表示单位时间单位换热面积在单位温度下的换热量，因而可以用来衡量各种对流换热过程换热性能的差异，这也就是对流换热系数这个定义沿用至今的道理。 对流换热系数如何确定呢 ? 分析一下流体在壁面上的特征也许会有帮助。前面已经提到，壁面上的流体分子层由于受到固体壁面的吸附是处于不滑移的状态，其流速应为零，那么通过它的热流量只能依靠导热的方式传递。由傅里叶定律传导的热流密度为，而从过程的热平衡可知，这些通过壁面流体层传导的热流量最终是以对流换热的方式传递到流体中去的，因而有 。于是得到如下关系 或 , (5 1 )式中， ， 为流体的导热系数 。 式 5 1 称为换热微分方程式 ，它给出了计算对流换热壁面上热流密度的公式，也确定了对流换热系数与流体温度场之间的关系 。它清晰地告诉我们，要求解一个对流换热问题，获得该问题的对流换热系数或交换的热流量，就必须首先获得流场的温度分布，即温度场，然后确定壁面上的温度梯度，最后计算出在参考温差下的对流换热系数。所以换热系数与流场的温度分布有关，因此，它与流速、流态、流动起因、换热面的几何因素、流体物性均有关 。 所以换热系数不是物性参数 。对流换热问题犹如导热问题一样，寻找流体系统的温度场的支配方程，并力图求解方程而获得温度场是处理对流换热问题的主要工作。由于流体系统中流体的运动影响着流场的温度分布，因而流体系统的速度分布（速度场）也是要同时确定的，这也就是说，速度场的场方程也必须找出，并加以求解。不幸的是， 对于较为复杂的对流换热问题，在建立了流场场方程之后，分析求解几乎是不可能的。此时，实验求解和数值求解是常常被采用的 。尽管如此，实验关系式的形式及准则的确定还是建立在场方程的基础上的，数值求解的代数方程组也是从场方程或守恒定律推导得出的。 下面我们将针对对流换热过程的流场从质量守恒定律、动量守恒定律和能量守恒定律出发结合傅里叶导热定律和斯托克（ Stockes ）黏性定律推导出流场的支配方程组。5-2对流换热问题的数学描写 本节要求： 1掌握对流换热问题完整的数学描写：对流换热微分方程组及定解条件；对流换热微分方程组：连续性方程+动量微分方程+能量微分方程； 2熟悉能量微分方程的推导方法及思路：对微元体应用能量守恒定律和傅里叶导热定律； 3掌握对流换热微分方程组中各项的意义。对流换热问题完整的数学描写包括对流换热微分方程组及定解条件。对流换热过程是流体中的热量传递过程，涉及流体运动造成的热量的携带和流体分子运动的热量的传导（或扩散）。因此，流体的温度场与流体的流动场（速度场）密切相关。要确立温度场和速度场就必须找出支配方程组，它们应该是，从质量守恒定律导出的连续性方程、从动量守恒定律导出的动量微分方程、和从能量守恒定律导出的能量微分方程。从一般意义上讲，推导这些方程应该尽量少的限制性条件。但是为了突出方程推导的物理实质而又不失一般性，这里选取二维不可压缩的常物性流体流场来进行微分方程组的推导工作。1. 连续性方程 图 5 2 给出了一个二维流体流场，从中选取一个微元体 ，并设定 x 方向的流体流速为 u ，而y 方向上的流体流速为 v ，流体的密度的 r 。将质量守恒定律应用于微元体，必然存在如下质量平衡关系： 单位时间流进和流出微元体的质量流量之差微元体质量随时间的变化率。 从 x 方向进入元体的质量流量为 ，流出则为 ；而 从 y 方向进入元体的质量流量为 ，流出则为 。 于是，单位时间流进和流出微元体的质量流量之差为 。 同时，微元体质量随时间的变化率应为 。从质量平衡关系可以得出连续性方程： 5 2 对于稳态流场 5 3 对于不可压常物性流场 5 4 2 动量微分方程 动量微分方程是纳维埃和斯托克斯分别于 1827 和 1845 年推导的。 流体的运动应服从动量守恒定律，对于我们所研究的二维不可压缩流场，微元体 的动量平衡关系应为： 流体流动引起的元体动量变化率 = 作用于元体上的外力之和 。 因流体运动引起的元体动量变化率 由图 5 3 可见，从 x 方向进入元体质量流量在 x 方向上的动量为 ，而从 x 方向流出元体的质量流量在 x 方向上的动量则为 ; 同时还应注意到，从 y 方向进入元体的质量流量在 x 方向上的动量为 ，而从 y 方向流出元体的质量流量在 x 方向上的动量则为 。 把离开元体的动量流量减去进入元体的动量流量，结果就是 x 方向上的动量改变量， ( 1a ) 注意，在化简过程中利用了连续性方程和忽略了高阶小量。 同理，我们可以导出 y 方向上的动量改变量， （ 1b ） 作用于微元体上的外力 作用于微元体上的力可以分为表面力和体积力。 体积力 是由于重力场、电场或磁场作用于元体上的结果。为了分析上的便利和简明，设定单位体积流体的体积力为 F ，那么相应在 x 和 y 方向上的分量分别为 F x 和 F y 。于是作用于微元体上的体积力在 x 方向为： ( 2a ) 而在 y 方向为： （ 2b ）。 表面力 是作用于微元体表面上的力。通常用作用于单位表面积上的力来表示，称之为应力。由于力作用的表面和作用力均为矢量，那么应力应该是一个二阶张量。在物理空间中面矢量和力矢量各自有三个相互独立的分量（方向），因而对应组合可构成应力张量的九个分量。于是应力张量可表示为 ，式中 为应力张量，下标 i 表示作用面的方向 , 下标 j 则表示作用力的方向。通常将作用力和作用面方向一致的应力分量称为正应力，而不一致的称为切应力。对于我们讨论的二维流场应力只剩下四个分量，记为，式中： 为 x 方向上的正应力（力与面的方向一致 ）； 为 y 方向上的正应力（力与面的方向一致 ）； 为作用于 表面上的 y 方向上的切应力；而 为作用于 y 表面上的 x 方向上的切应力。 从图 5 4 显示的微元体上的应力作用情况可以看出： 作用在 x 方向上表面力的净值为 ； (a) 而作用在 y 方向上表面力的净值为 。 (b) 由于流体黏性的作用，在应力的作用下流体的微元可以发生相应的变形。斯托克斯 提出了归纳速度变形率与应力之间的关系的黏性定律 ，即： ; ; （这是针对二维问题的形式），式中 为流体的动力黏性系数。将它们代入（ a ）、（ b ）两式，得出作用在微元体上表面力的净值表达式： 在 x 方向上为 ; ( 3a ) 在 y 方向上为 。 (3b) 现在将式（ 1a ）、（ 2a ）、（ 3a ）和 (1b) 、 (2b) 、 (3b) 分别代入动量平衡式，经整理得出分量形式的动量微分方程式： 在 x 方向上 ； 5 5a 在 y 方向上 。 4 5b 惯性力 体积力压力 粘性力 这就是二维不可压缩常物性流体的动量微分方程式，它是流场速度分布的支配方程。通过与连续性方程联立，在给定的初、边值条件下可以求出流场的速度分布和压力分布。由于动量方程产生于微元体的动量守恒，因而方程各项的物理意义是十分明确的。方程的左边表征流场的惯性力（亦为动量的当地改变量与位移改变量），方程右边的第一项表征流场的体积力、第二项表征流场静压力的变化，而最后一项表征流场的黏性力（亦为黏性扩散引起的动量变化）。方程 5 5 可以改写成， 在 x 方向上 ； 5 6a 在 y 方向上 。 5 6b 式中记号 表示流场的全导数 或称 真导数 ，表示了在流场中物理量随时间的真实改变的速率。设流场中某物理量 f ，其全微分为 ，其总的变化率也就是全导数则为 ，式中， 分别为 x,y, z 三个方向上的流体流速。从物理意义上理解， 表示物理量随奔涞牡钡乇浠剩?则表示因流体运动而造成的物理量随时间的变化率。如果是一个稳态系统则有 ；如果是一个固体系统则有 。 3 能量微分方程 流场中的温度分布无疑反映了流场能量分布的状态，受着能量守恒定律的制约。因而支配流场温度场的场方程能量微分方程可以通过对流场中微元体进行能量平衡分析而得出。对于二维不可压缩常物性流体流场而言，微元体 的 能量平衡关系式为： 式中： 为以传导方式进入元体的净的热流量； 为以对流方式进入元体的净的热流量；为元体粘性耗散功率变成的热流量； 为元体的焓随时间的变化率。 下面我们将导出微元体能量方程相应的各项。 （1）以传导方式进入元体的净的热流量 由图 5 5 可知， x 方向和 y 方向净导入元体的热流量分别为 。引入傅立叶定律有， ，代入上式可以得到总的净导入元体的热流量： （ 1 ） （2）以对流方式进入元体的净的热流量 从图 5 6 可知，从 x 方向流进元体的热流量为 ，而从 x 方向流出元体的热流量为 ，略去高次项后，从 x 方向进入元体的净热流量为： 从 y 方向流进元体的热流量为 ，而从 y 方向流出元体的热流量为 略去高次项后，从 y 方向进入元体的净热流量为： 那么，以对流方式进入元体的热流量净值为 。 （ 2 ） 应注意到，在上式的化简过程中利用了连续性方程和忽略了高阶小量。 （3）元体粘性耗散功率变成的热流量：元体粘性耗散功率变成的热流量在此不作详细的推导，有兴趣的读者可参阅相关文献。这里直接给出结果： ， （ 3 ） 也可以简记为 ， 为粘性耗散函数。 （4）元体的焓随时间的变化率 元体中流体的焓随时间的变化率 , 对于不可压缩常物性的流体，可按其定义写出： 。 （ 4 ） （5）能量微分方程 将所导出的各项代入微元体能量平衡方程，且两边消去 , 经整理得出 。 57 流体能量随时间的变化 + 对流项 = 热传导项 + 热耗散项 这就是二维不可压缩常物性流场的能量微分方程，对于稳态流场方程变为 ， 58 亦可写为 。 59 从前面的推导过程不难看出，方程 5 7 各项的物理意义是十分明确的。方程右边三项中，第一项为流体能量随时变化项，另外两项为流体热对流项；方程左边第一项为热传导（热扩散）项，第二项为黏性耗散项。当流体不流动时，流体流速为零，热对流项和黏性耗散项也为零，能量微分方程式便退化为导热微分方程式，即 。所以，固体中的热传导过程是介质中传热过程的一个特例。 4 对流换热微分方程组 上面我们基于质量守恒、动量守恒和能量守恒的原则分别导出了连续性方程、动量方程和能量方程，它们是支配对流换热过程的场方程。可以归纳如下： 5 10 。 以下就是二维、常物性、不可压缩的、无内热源条件下，牛顿流体对流换热问题的 对于给定的流场在相应的初边值条件下，联立求解连续性方程和动量方程可以获得流场的速度分布和压力分布。在速度场已知的情况下求解能量微分方程，最终可以获得流场的温度分布。此时，再引入换热微分方程 (n 为壁面的法线方向坐标 ) ，最后可以求出流体与固体壁面之间的对流换热系数，从而解决给定的对流换热问题。 十分令人遗憾的是，对于大多数对流换热问题，尤其是流体流动状态从层流转变为紊流之后的换热问题，采用直接求解微分方程的分析办法几乎是不可能的。因此，对流换热问题的求解往往是一件较为复杂的工作。通常求解对流换热问题有如下几个途径： 分析求解：主要针对一些简单问题，如二维的边界层层流流动、库特流动和管内层流流动换热等，都可以通过数学分析的办法来求解。具体的求解方法同学们可以通过阅读传热学方面的书籍而获得。 实验研究：由于对流换热的复杂性，实验研究是求解对流换热问题的主要方法，尤其是对于紊流换热问题、有相变的换热问题，或者几何结构复杂的换热问题，实验求解几乎是唯一的途径。虽然，数值分析方法得到发展，但其结果还是要通过实验来加以验证。因此，本教材将主要讨论对流换热过程的实验研究方法和所得的实验关系式的应用。 数值求解：随着计算机应用的普及和数值计算方法的发展，对流换热过程的数值分析越来越成为一种主要的求解方法，其结果的可信度越来越高。数值求解方法主要是将对流换热方程组在离散的控制体中变为代数方程组，然后编制出相应的计算机程序，通过计算机求出离散的温度分布，用于表示计算区域连续的温度分布。由于对流换热过程的数值分析较为复杂，作为一本入门的教材不可能对其进行讨论，有兴趣的同学们可以参阅流体流动和传热数值计算方面的文献。53边界层的概念及边界层微分方程组 本节重点要求： 1掌握边界层的概念：及； 2熟悉边界层的主要特点及引入边界层概念的意义，明确及（及）是同一数量级的量，且边界层内速度梯度和温度梯度很大，即都很大；3掌握二维、稳态、常物性、不可压缩、不计重力、无内热源的对流换热边界层微分方程组；4掌握特征数：、的概念及物理意义。一、边界层的概念边界层的概念是1914年普朗特提出的。他认为对流换热热阻的大小主要取决于靠近壁面附近流体的状况,因为粘滞性起作用的区域仅仅局限在靠近壁面的薄层内。在此薄层以内， u、t变化最为剧烈。 为计算的简化，因此引入了速度边界层和温度边界层的概念。1. 速度边界层 1）定义：流体流过固体壁面时，由于壁面层流体分子的不滑移特性，在流体黏性力的作用下，近壁流体流速在垂直于壁面的方向上会从壁面处的零速度逐步变化到来流速度，如图510所示。流体流速变化的剧烈程度，即该方向上的速度梯度，与流体的黏性力和速度的大小密切相关。普朗特通过观察发现，对于低黏度的流体，如水和空气等，在以较大的流速流过固体壁面时，在壁面上流体速度发生显著变化的流体层是非常薄的。因而他把在垂直于壁面的方向上流体流速发生显著变化的流体薄层定义为速度边界层，而把边界层外流体速度变化比较小的流体流场视为主流区(势流流动区域)。这样，引入边界层的概念之后，流体流过固体壁面的流场就人为地分成两个不同的区域，其一是边界层流动区，这里流体的黏性力与流体的惯性力共同作用，引起流体速度发生显著变化；其二是主流区(势流区)，这里流体黏性力的作用非常微弱，可视为无黏性的理想流体流动，也就是势流流动。2）边界层的厚度我们说边界层是壁面上方流速发生显著变化的薄层，但其边缘所在的位置人们却是模糊的。在实际分析边界层问题时通常规定，当速度变化达到主流速度的99%处的空间位置y为速度边界层的外边缘，那么从这一点到壁面的距离就是边界层的厚度。3）引入的意义：a）缩小计算区域。以温度为20的空气沿平板的流动为例，在不同来流速度下，相对于平板长度l，是一个比l小一个数量级以上的小量。在这样小的薄层内，流体的速度要从0 m/s变化到接近于主流流速，所以流体在垂直于主流方向上的速度变化是十分剧烈的。由于边界层内都很大，动量微分方程式中的粘性力和惯性力及能量微分方程中的导热和对流项都需考虑，有了边界层，可以将对流换热的研究集中在边界层内。b）边界层内的流动和换热可以利用边界层的特点进行简化，即在边界层中可以运用数量级的分析方法对纳维斯托克斯方程进行简化。例如：运用数量级分析的方法，边界层内粘性流体的稳态动量方程可简化为： 与二维稳态的纳维斯托克斯方程相比，上述运动微分方程的特点是：(i)在u方程中略去了主流方向的二阶导数项；(ii)略去了关于速度的动量方程。(iii)认为边界层中。4）临界雷诺数：流体的流动可区别为层流和湍流两种。边界层也呈现出层流和湍流。流体以的流速沿平板流动。在平扳的起始段，很薄。随着x的增加，由于壁面粘滞力的影响逐渐 向内部传递，逐渐增厚，同时黏性力对流 场的控制作用也逐步减弱，从而使边界层内 的流动变得紊乱。在某一距离xc以前一直保 持层流特征。此时流体作有秩序的分层运动， 各层之间互不干扰。这时的边界层称层流边界层。随着边界层厚度的增加，本来很平顺的层流流动状态就会变成紊乱无序的紊流流动状态。此时流体质点在沿x方向流动的同时，又作着紊乱的不规则脉动，紊流流动状态的边界层则为紊流边界层，我们把边界层从层流过渡到紊流的x值称为临界值，记为xc，其所对应的雷诺数称为临界雷诺数，即Rec。实验研究的数据表明，流体平行流过平板的临界雷诺数大约是（一般对平板：Rec500000）。当然，这一数据与来流速度的紊流程度和平板前沿的几何形状密切相关，因而临界雷诺数的数值会在一定范围内变动。 我们不难发现，要使边界层的厚度远小于流动方向上的尺度（即l），也就是所说的边界层是一个薄层，这就要求雷诺数必须足够的大。因此，对于流体流过平板，满足边界层假设的条件就是雷诺数足够大。由此也就知道，当速度很小、黏性很大时或在平板的前沿，边界层是难以满足薄层性条件。 4）边界层理论的四个基本要点 a. 当粘性流体沿固体表面流动时，流场划为主流区（势流区）和边界层区。在边界层区内，速度在垂直于壁面方向剧烈变化。而主流区速度低度几乎为零。速度边界层成立的条件是Re1。 b. 主流区的流动视为理想流体的流动，用描述理想流体的方程求解。边界层区应考虑粘性的影响，用粘性流体的边界层微分方程求解，其特点是主流方向流速的二阶导数项忽略不计。 c. 边界层的流动状态分为层流和紊流。紊流边界层分为层流底层、缓冲层和紊流核心区。 d. 边界层的厚度是与壁面尺寸相比很小的量，且远不只小于一个数量级。即：。2. 热（温度）边界层1）定义当流体流过平板而平板的温度tw与来流流体的温度t不相等时（即存在温度差），对于上述的低黏性流体，如果流体的热扩散系数也很小，在壁面上方也能形成温度发生显著变化的薄层，常称为热边界层。2）热边界层厚度 仿照速度边界层的约定规则，以过余温度为来流过余温度的99%处定义为热边界层的外边缘。（当壁面与流体之间的温差达到壁面与来流流体之间的温差的0.99倍时，即 ，此位置就是边界层的外边缘），而该点到壁面之间的距离则是热边界层的厚度,记为t 。如果整个平板都保持温度tw，那么，x=0时t(x)=0，且随着x值的增大逐步增厚。在同一位置上热边界层厚度与速度边界层厚度的相对大小与流体的普朗特数Pr有关，也就是与流体的热扩散特性和动量扩散特性的相对大小有关。理论上 tf(Re， Pr) 由此式可以看出，热边界层是否满足薄层性的条件，除了Re足够大之外还取决于普朗特数的大小，当普朗特数非常小时,热边界层相对于速度边界层就很厚，反之则很薄。二、边界层微分方程组利用上述的边界层的概念，我们可以对流体流过平板的对流换热微分方程组进行相应的简化。按照普朗特的边界层为一个薄层的假设，以及满足这一假设下的流场特征，我们运用数量级分析的方法可以对上述微分方程组在边界层中得以简化。数量级分析：比较方程中各量或各项的量级的相对大小；保留量级较大的量或项；舍去那些量级小的项，方程大大简化由边界层假设，如果设定X的数量级为1,(x 与 l 相当)，那么Y的数量级定义为（一个小量）；在设定主流方向上的速度u 的数量级为1, 温度： 的情况下，对动量微分方程进行数量级分析。对于稳态过程而言，如不考虑体积力，则动量方程为： 分析：因为，则1，；来流速度1（1）比较与：若的数量级为即1时，则的数量级为，相应地，而1（2）分析流体垂直于板面方向的速度v相当于什么数量级？即必与属于同一数量级，1可以得出v的数量级为。在边界层这样一个极薄的薄层中，流体的垂直速度从平板表面处的零，增加到边界层外缘处的，而其变化率的数量级却仍与相同，这就表明速度v的数值一定大大地低于u的数值，即v。这样v在x方向上的变化率也必远远小于u在x方向上的变化率，既然1，则。同理，既然1，。同（1）分析，也必更加大大地大于，既然，则。综上所述，动量方程的数量级可分别表示成： a 1 b 1 上两式右边括号内的两项相比，后一项要比前一项大得多，所以和均可略去。在边界层内，粘性力与其它力属于同一数量级，从上两式左边可以看出，惯性力分别相当于1和数量级，按照方程两边数量级应该一致的原则，因此粘性力也必相应地相当于1和数量级，即运动粘度2。（3）考察压力梯度的数量级：在a式中惯性力和粘性力都相当于1的数量级，假定的数量级是1，那么的数量级至多也是1，否则方程式无法成立。同理，的数量级至多也是，（注意：与不同，紧贴固体表面的流体流速u等于零，但此处的压力并不等于零。边界层断面上的压力变化并非从零变到边界层外的压力p，因而，但的数量级却至多只是。）这一结果告我们在边界层中压力不随Y的变化而变化，仅仅是X的函数。（4）比较a式和b式的数量级后，发现b式的各项要比式a的各项小一个数量级，因而b时可以舍去，于是边界层的动量微分方程就由两个变为一个，即 而边界层内的压力梯度仅沿 x 方向变化，而边界层内法向的压力梯度极小。即边界层内任一截面压力与 y 无关而等于主流压力，则 当流体纵掠平板时，边界层外主流中的流速没有变化，按照铂努力方程，主流中的压力也不变，可以得出： 同样对边界层中二维稳态能量方程进行数量级比较，有 至此，二维、稳态、无内热源的边界层换热微分方程组为：微分方程组经过在边界层中简化后，由于动量方程和能量方程分别略去了主流方向上的动量扩散项和热量扩散项 ，从而构成上游影响下游而下游不影响上游的物理特征。这就使得动量方程和能量方程变成了抛物型的非线性微分方程；且由于动量方程由两个变成为一个，而且项可在边界层的外边缘上利用伯努利方程变成的形式。于是方程组在给定的边值条件下可以进行分析求解。布劳修斯在引入相似变量下把动量方程变成常微分方程后由进行了求解，而普尔豪森也对能量方程进行了相同的处理。获得相应的速度分布和温度分布，进而求得壁面的摩擦系数和对流换热系数。普朗特从他多年从事水力学试验中所观察到的事实出发，创造性的用数量级对比法简化了原始的微分方程组，开拓了对流换热理论界的道路，成为流体力学和传热学发展史的里程碑。在他发表边界层微分方程组时，许多著名数学家都不以为然，认为他的这种舍弃微分方程中若干项以至整个方程的做法是荒谬的。（1904年发表）直到1908年，他的学生布劳修斯用边界层方程获得了外掠平板的解，并且其正确性被尔后的实验证实，普朗特的边界层理论才站住脚。所以实验研究仍是技术发展必不可少的。三、速度边界层厚度与热边界层厚的关系1。比较边界层的动量方程和能量方程，动粘性反映了流体中由于分子运动而扩散动量的能力。在忽略动量方程压力项后，当Pr=1时，动量方程与能量方程完全相同。即速度分布的解与温度分布完全相同，此时速度边界层厚度等于温度边界层厚度。当Pr1时，Pr=/a，a，粘性扩散 热量扩散，速度边界层厚度温度边界层厚度。当Pr1时，Pr=/a，a，粘性扩散 热量扩散，速度边界层厚度温度边界层厚度。这也可以从公式得出。故：湍流换热比层流换热强！层流：温度呈抛物线分布湍流：温度呈幂函数分布湍流边界层贴壁处的温度梯度明显大于层流d 与 dt 的关系：分别反映流体分子和流体微团的动量 和热量扩散的深度 2求解示例流体外掠等温平壁的对流换热：例如：对于主流场均速、均温，并给定恒定壁温的情况下的流体纵掠平板换热，即边界条件为求解上述方程组(层流边界层对流换热微分方程组)，可得局部表面传热系数的表达式 3特征方程及努塞尔特数式中：努塞尔(Nusselt)数雷诺(Reynolds)数普朗特数注意：特征尺度为当地坐标x5-4 边界层积分方程组及其求解本节要求：1掌握边界层积分方程组的求解 （1）边界层积分求解的基本思想； （2）求解2熟悉相似原理及量纲分析 （1）物理量纲相似的性质 （2）相似准则之间的相似关系 （3）判别现象相似的条件 （4）获得相似准则数的方法边界层积分方程组及比拟理论：（1）1921年，冯卡门提出了边界层动量积分方程。（2）1936年，克鲁齐林求解了边界层能量积分方程。（3）近似解，简单容易。一、 边界层动量积分方程 边界层动量积分方程是把动量定律应用于一个 控制容积 导出的。取 常物性、不可压缩流体的二维稳态强制对流为对象 作分析。在流体中划出一个如图 1 所示的控制容积，它包括 dx 一段边界层，而 z 方向为单位长度。控制容积左侧面为 ab 右侧面为 cd ，顶面为 bd ，底面为壁面的 ac 部分，即取 ac 为 dx 。由于在边界层内 y 方向上的流速很小，因此推导中 只考虑 x 方向上的动量变化 ， 不引入流速 v 。图 1 给出了速度的分布曲线。在距壁面 y 处流速为 u ，在 y 处 u=u 。 先计算单位时间内出入控制容积的动量之差。为此计算以下各项： 穿过控制面 ab 进入控制容积的动量为 而同时穿过 cd 面流出的动量为 净流出的动量为 （ a ） （ 2 ）没有流体穿过固体表面 ac 。但有流体质点穿过 bd 面。 根据质量守恒，穿过 bd 面流入控制容积的质量流量等于流出 cd 面与流入 ab 面的质量流量之差。流入 ab 面的质量流量为 ，流出 cd 面的质量流量 。于是 穿过 bd 面流入控制容积的质量流量为 相应带入控制体的动量（略去u 沿x变化引入的高阶导数项）为（ b ） 根据动量定律，在 x 方向上的动量变化必须等于 x 方向上作用在控制体表面上外力的代数和。 作用在控制体表面上 x 方向上的外力，有作用于 ac 面上的 切应力 w dx 以及 ab 和 cd 两面压力之差 于是动量定律可表达为 （ c ） 由于存在以下关系： （ d ） 于是式 (c) 可改写成为 重新组合可得 （ e ） 由伯努利方程知 代入 (e) 式，得 根据边界层理论，在边界层外的主流区 u -u=0 。改写上式积分上限得 （ 1 ） 这就是卡门在 1921 年导出的边界层动量积分方程。由积分方程求出的分析解称为近似解 ，以区别于微分方程的精确解。 二、边界层能量积分方程 把能量守恒定律应用于控制容积可推导出边界层能量积分方程。 x 方向上为 dx ， y 方向上大于流动边界层即热边界层厚度，而 z 方向上为单位长度的一个控制容积如图 2 所示。 在常物性、流速不致引起耗散热的条件下，考察控制容积的能量守恒。在边界层数量级分析中已经得出结论： 故推导中仅考虑 y 方向上的导热 。下面计算穿过控制容积各个面的热量。 （1）单位时间内穿过 ab 面进入控制容积的热量为 （f） 单位时间内穿过 cd 面带出控制容积的热量为 (g) （2）单位时间内穿过 bd 面进入控制容积的质量流量为 由它带入控制容积的热量为 （ h ） （3）穿过 ac 面，因贴壁流体层导热带出控制容积的热量为 (i) 在稳态条件下，根据能量守恒进入与带出控制容积的热量相等，于是可得 整理后得 因为在热边界层以外 t -t=0 ，上式积分上限可改为 t ，得 （ 2 ） 这就是边界层能量积分方程。它与边界层动量积分方程一起组成对流换热边界层积分方程组。 1936 年，克鲁齐林用积分方程组求解对流换热，热到了满意的结果。 三、边界层积分方程组求解示例 作为边界层积分方程组求解的示例，仍以 稳态常物性流体强制掠过平板层流时的换热作为讨论对象 。壁面具有定壁温的边界条件。 在常物性条件下。动量积分方程不受温度场的影响，可先单独求解，解出层流边界层厚度及摩擦系数，然后求解能量积分方程，解出热边界层厚度及换热系数。 1. 求解流动边界层厚度及摩擦系数 在本问题中， u 为常数，动量积分方程式（ 1 ）左边的第二项为 0 。再引入 ， 式（ 1 ）可改写成为 （ 3 ） 为求解上式，还需补充边界层速度分布函数 u=f(y) 。选用以下有 4 个任意常数的多项式作为速度分布的表达式 u=a+by+cy 2 +dy 3 式中， 4 个待定常数由边界条件及边界层特性的推论确定，即 y=0 时 u=0 且 y= 时 u= u 且 由此求得 4 个待定常数为 a= 0 c=0 于是速度分布表达式为 (4 ) 从式求得 du/dy 并代入切应力 w 的定义式得 （5） 将式（ 4 ）（ 5 ）代入（ 3 ），积分得 (j) 分离变量，并注意到 x=0 时 =0 ，可得 所以 （ 6 ） 其 无量纲表达式为 （ 7 ） 其中 Re x = u x / , 其特性尺度为离平板前缘的距离 x 。 将式（ 6 ）代入式（ 5 ）可得在 x 出的壁面局部切应力 （ 8 ） 在工程计算中常使用局部切应力与流体动压头之比，称为摩擦系数，亦称范宁摩擦系数，其表达式为 （ 9 ） 求解动量微分方程可以获得/x 、 cf 的精确解： ， 。可见，与微分方程的精确解相比，由积分方程得到的 /x 、 c f 只在常数值上略有差别。 2. 求解热边界层厚度及换热系数 先求解热边界层厚度。为从式（ 2 ）求解热边界层厚度，除 u=f(y) 已由式（ 4 ）确定外，还需要补充热边界层内的温度分布函数 t=f(y) 。对此，亦选用带 4 个常数的多项式： t=e+fy+gy 2 +hy 3 式中， 4 个待定常数由边界条件及热边界层特性的推论确定，即 y=0 时 t=t w 且 y= 时 t= t 且 由此求得 4 个待定常数为 e=t w g=0 若用以 t w 为基准点的过余温度 =t-t w 来表达，则温度分布表达式为 （ 10 ） 能量积分方程式（ 2 ）用过余温度表示为 （ 11 ） 进一步求解中， 令热边界层厚度与流动边界层厚度之比 t / = ，并假定 1 的流体显然是适用的。 用速度分布式（ 3 ）及温度分布式（ 10 ）分别求得式（ 11 ）左方积分部分及右方的 如下： （ k ） (l) 因为 1 ，包含 4 的高次方项相对于包含 2 的项可略去不计。把上列关系式代入（ 11 ）整理后得 可改写成 将式（ 5 ）的 及式（ j ）的 关系代入，整理后得 （ m ） 令 3 =Y ，上式可改写成 次微分方程的通解为 (n) 此处 C 必须等于 0 ，否则 x=0 时 Y 成为不定值，不符合物理现实。于是得 (12) 以上结果是在 Pr 1 的前提下推得的。气体类流体的 Pr 数略小于 1 。这类流体的 Pr 数的最小值约为 0.65 ，此时 =1.12 ，式 (k) 中 略去 4 项引入的误差不大。据此，式（ 12 ） 对常用气体亦可近似适用。液态金属类流体的 Pr 数的数量级是 10 -2 ， 式（ 12 ）不能适用。 其次求解局部换热系数 x 。将温度分布式（ 10 ）代入换热微分式，可得局部换热系数 x 为 将式（ 5 ）和（ 12 ）的 、关系代入上式得 （ 13 ） 其无量纲表达形式为 （ 14 ） 式中， ，称为局部努谢尔特数，上两式与波尔豪森所得的精确解相吻合。 最后求解平均换热系数 。为了求得长 l 一段平板的平均换热系数，可求 x 在 0 到 l 范围内的积分平均值，即 （ 15 ） 式中 B 代表（ 13 ）右方除 x 以外的其它量和常数。式（ 15 ）表明， 平均换热系数是 l 处局部换热系数的两倍。 于是 平均换热系数得表达式为 （ 16 ） 式中，无量纲数 Nu 、Re 中的特征尺度为平板长度l 。计算物性参数用的定性温度为边界层平均温度 。 四、边界层积