流体力学课件第9章：边界层流动

第九章边界层流动9.1

问题的提出工程中绝大多数流动都处于高Re区域，即惯性力>>粘性力，粘性力可忽略理想流体。

但固壁附近要满足无滑移条件，速度低，通常Re较小，粘性力不能忽略。普朗特（1904）提出边界层概念，边界层内为粘性流动，速度梯度较大，边界层外为理想流体势流运动。1、边界层概念的提出理想流体力学在早期较成功地解决了与粘性关系不大的一系列流动问题，诸如绕流物体的升力、波动等问题，但对绕流物体阻力、涡的扩散等问题，理想流体力学的解与实际相差甚远，且甚至得出完全相反的结论，圆柱绕流无阻力的达朗贝尔佯谬就是一个典型的例子。如何考虑流体的粘性，怎样解决扰流物体的阻力问题，在当时是一个阻碍流体力学发展的难题，1904年，普朗特（Prandtl）通过大量实验发现，虽然整体流动的Re数很大，但在靠近物面的薄层流体内，流场的特征与理想流动相差甚远，沿着法向存在很大的速度梯度，粘性力无法忽略。Prandtl把这一物面近区粘性力起重要作用的薄层称为边界层（Boundarylayer）。9.2

边界层流动的基本特征Prandtl边界层概念的提出，为人们如何计入粘性的作用开辟了划时代的途径，因此称其为粘性流体力学之父。对整个流场提出的基本分区是:（1）整个流动区域可分成理想流体的流动区域（势流区）和粘性流体的流动区域（粘流区）。（2）在远离物体的理想流体流动区域，可忽略粘性的影响，按势流理论处理。（3）粘性流动区域仅限于物面近区的薄层内，称为边界层。边界层内，粘性力的作用不能忽略，与惯性力同量级，流体质点作有旋运动。例：空气绕某一翼型的流动，整个流场可分为边界层、边界层脱离翼型物面以后形成的尾流、以及边界层和尾流以外的势流。图9-1翼型绕流实际流体绕任何形状物体的大雷诺数流动都会在物面附近形成边界层。

边界层流动可以是层流或湍流。实际中更一般地是混合边界层，即边界层前缘为层流，经过一过渡区(称为转捩区)后转变为湍流；在湍流区，紧挨物面附近还有一层流底层。图9-2所示为一均匀来流绕过平板一侧所形成的边界层流动。图9-2平板边界层流动2、边界层的流态

在湍流区，若平板表面粗糙度D大于层流底层的厚度dl，则称之为粗糙(表面)平板；否则称为光滑(表面)平板。当层流区的范围很小时，可近似地把整个边界层看成为湍流边界层。为了便于判断边界层的流态，通常假定由层流到湍流的转捩是在某一截面突变完成的，并称此截面为临界截面，它离边界层前缘的距离称为临界长度x*，临界截面边界层的厚度称为临界厚度d*。

以平板边界层为例边界层流态用临界雷诺数Re*来判断。

Re*有两种形式：

Rex*=U∞x*/v

Red*=U∞d*/v对于平板绕流，

Rex*=5105~3106，

Red*2800。

速度为U∞的来流进入前缘后，由于物体的粘附作用，低层流体的速度变为0。随着x的增加受阻滞的流体在y方向上逐渐扩大，以致行成一个有明显速度变化的区域，通常称为速度边界层。3、边界层的定义与特征

（1）边界层定义严格而言，边界层区与主流区之间无明显界线，通常以速度达到主流区速度的0.99U∞作为边界层的外缘。由边界层外缘到物面的垂直距离称为边界层名义厚度，用δ表示。注意：边界层的外边界不是流线；流线可以穿越边界层外边界进入边界层内。（2）边界层的有涡性

粘性流体运动总伴随涡旋的产生、扩散、衰减。边界层就是涡层，当流体绕过物面时，无滑移边界条件相当于使物面成为具有一定强度的连续分布的涡源。以二维流动为例说明。此时，物面上的涡源强度为：

对于不可压缩流体，二维流动的涡量输运方程为：

上式表明，由于粘性的影响，物面上的涡量一方面沿垂直流线方向扩散，另一方面，涡量沿主流方向迁移，并随之而逐渐衰减。涡量的扩散速度与粘性有关，涡量的迁移速度取决于流动速度。（3）边界层厚度的量级估计根据边界层内粘性力与惯性力同量级的条件，可估算边界层的厚度。以平板绕流为例说明。设来流的速度为U，在x方向的长度为L，边界层厚度为δ。

惯性力：粘性力：

由惯性力与粘性力同量级得到由此可见，在高Re数下，边界层的厚度远小于被绕流物体的特征长度。（4）边界层各种厚度定义a.边界层排挤厚度在边界层内，理想流体的质量流量为：

其中，U为边界层外缘速度。由于粘性的存在，实际流体通过的质量流量为：上述两项之差表示粘性存在而损失的流量，这部分流量被排挤到主流场中，相当于主流区增加了一层流体。主流区所增加的厚度称为排挤厚度：图7-3边界层位移厚度①

代表整个边界层内亏损质量流量与无粘性流动时单位厚度的质量流量之比，越大表明边界层引起的质量流量亏损越大。②

为了保持有粘性与无粘性流动的质量相等，在用无粘性理论设计管道时应将管壁向外扩大③

根据边界层速度分布渐近的概念可以导出：b.边界层动量损失厚度在边界层内，在质量流量不变的条件下，理想流体通过的动量为：由于粘性的存在，实际流体通过的动量为：上述两项之差表示粘性存在而损失的动量，这部分动量损失用外流流速U（理想流体）折算的动量损失厚度为：c.边界层能量损失厚度在边界层内，在质量流量不变的条件下，以外流速度（理想流体）通过的动能为：由于粘性的存在，实际流体通过的动能为：上述两项之差表示粘性存在而损失的动能，这部分动能损失用主流流速U（理想流体）折算的动能损失厚度为:上述各种厚度的计算公式，对于不可压缩流体而言，变为:

一般有：9.3边界层微分方程假设：不可压平面流动，不计质量力，边界层内，选取特征长度L，速度尺度ue，时间尺度t=L/ue，边界层近似假定：（1）根据边界层定义，纵向偏导数远远小于横向偏导数。（2）法向速度远远小于纵向速度。（3）边界层内的压强与外流速度的平方成正比。将这些量级关系式代入到N-S方程中，得到

无量纲化与数量级分析

比较各项后的量级后得：可见，P与y无关，即沿y均匀，可取：对于边界层问题，是已知函数,可由外部势流解得到：对定常流动，由伯努利方程得到：还原方程后得边界层方程组（普郎特边界层运动方程，1904）：定解条件：

仍然是非线性方程，只有一些特别情况可以得到解析解。9.4边界层方程的相似性解边界层方程的相似性解：任意x截面的流速u(x，y)采用无量纲坐标表示时，其速度分布均相同。定义无量纲坐标（相似变量）η是x，y的函数，h(x)函数是y坐标的特征长度。边界层方程具有相似解时，有ue(x)是外部势流速度，f(η)是某一特定函数。对二维不可压流动，存在流函数，设流函数的形式可以写为速度分量为将速度u，v的一次和二次导数都转换成关于f，h及其各阶导数的函数关系，然后代入边界层方程，可得(I)其中边界条件为当α，β为常数时，方程（I）是常微分方程，即流动具有相似性解。(I)福克勒-斯坎（Falkner-Skan）方程特征长度相似变量流函数几种存在相似解的流动：(1)平板边界层流动（β=0，m=0）(2)平面驻点流动（α=1，β=1，m=1）(3)二维楔形物体驻点附近流动（夹角πβ）(4)绕钝头柱体前驻点附近的二维流动（β=1，m=1）与平面驻点流动相同。(5)陡坡流动（β=2/3，m=1/2）

1908年，Prandtl学生Blasius利用边界层速度分布的相似性求解了平板层流边界层方程。对于零压梯度（）、定常、不可压缩流体平板层流绕流，边界层方程为:相应的边界条件为:

9.5平板边界层的Blasius精确解由于ue=U∞，所以取积分，得流函数速度分量由此得到:代入边界层方程，得到:化简后变为:边界条件为:Blasius用无穷级数进行了求解。假设：其中，为待定系数。

由边界条件，可得:（1）边界层厚度（）（2）边界层排移厚度（3）边界层动量损失厚度（4）壁面切应力:（5）壁面摩擦阻力系数（6）平均壁面摩擦总阻力系数

郭永怀（1953年）对平板前缘点的修正，得到适用范围：

平板边界层方程为三阶非线性常微分方程，无解析解，布拉修斯最早给出了级数解，后来霍华斯等人用龙格—库塔法得到了数值积分解。平板层流边界层精确解数值表图7-4平板边界层速度剖面

边界层问题评述

边界层问题求解十分麻烦，上面只是最简单的一种情形。此问题因存在相似性解（指无量纲的速度分布相同），故可采用相似变量置换，将方程组化为常微分方程。类似问题还有轴对称层流边界层。更多问题不存在相似性解（如）此时可采用:1）级数近似解法；

2）动量积分近似解法；

3）纯数值解法。湍流边界层更为复杂，只能采用:1）半经验理论；

2）数值解。

流体绕物体的流动可分为成势流区和边界层区域，而势流区可以使用位势理论求解，而求解边界层则较困难。描述边界层粘性流动的是纳维-斯托克斯方程（N-S方程）连续性方程：

9.6

边界层积分方程解法

在工程中许多近似解法有很大的实用价值，其中动量积分方程是最简单而又最实用的一种方法(冯.卡门，1921)。CD流出的动量：AB流入的动量：边界层动量积分方程AC流入的动量：控制ABCD内流体的动量变化率（物质导数）为:经AC流入控制体的质量=CD-AB（质量流量）x方向控制体的受力情况：

AB面上受力：CD面上受力：

壁面BD上的剪切应力：x方向上的总的受力：（略去高阶小量）应用动量定理，得动量积分关系式：上式中压力梯度一项可以按主流为势流，利用欧拉方程来计算，不考虑质量力时（x方向）：即AC面上的分量：（粘性切应力可忽略）又因为:得动量积分关系式：变形后得另一种形式的动量积分关系式：利用排挤厚度和动量损失厚度的定义，得动量积分方程的一般形式：上式适用于稳定流动，不可压缩流体层流或湍流时的动量方程，根据它可以计算附面边界层的厚度和流动阻力。方程式中有,为此还需要两个补充条件：边界层内的速度分布：与速度分布有关的与δ的关系:动量积分关系不仅实用于层流，对湍流也实用。例题：密度为ρ的均质不可压缩流体以速度V∞流向板面与速度方向平行的平板，设流动定常，壁面附近速度为线性分布，不计质量力，板宽为D，求流体作用在板上的切向力。

解：前缘流进控制体的动量：后缘流出的动量：由y=δ流出控制体的动量流量是：上缘流出质量流量为：∴流体作用平板的总切应力为：平板边界层的积分近似解法对于平板问题，主流为均匀流，即：故动量积分方程式简化为：要求解上式可以先假设u和τ0的表达式，再求出δ(x)的表达式，u和τ0的近似表达式越接近于实际，δ(x)越精确。假设平板边界层上的速度分布可以近似地用高次抛物线表示：壁面切应力为：多项式系数应满足一定的条件：

①在壁面上：

②在外边界：

③在外边界：

④边界层的微分方程为：时

得：

解得的4个系数是：于是层流边界层的速度表达式：根据动量厚度定义：代入：

得：

利用x=0，δ=0的边界条件，积分得排移厚度：

动量损失厚度：

壁面的切应力：长度L，单位宽度平板的摩擦力：

阻力系数：平板边界层近似解与精确解比较作业：9-159-189.7

曲面边界层

平板边界层是无压强梯度的边界层，而曲面边界层内外缘的速度和压力均有变化，这种压强梯度的存在会影响边界层内的流动，最重要的就是在一定条件上会造成边界层的分离。

边界层的分离现象

边界层中的流体质点受惯性力、粘性力和压力的作用。粘性力的作用始终是阻滞流体质点运动，使流体质点减速，失去动能；压力的作用取决于绕流物体的形状和流道形状，顺压梯度有助于流体加速前进，而逆压梯度阻碍流体运动。以圆柱绕流为例说明边界层的分离现象。

对于理想流体，流体微团绕过圆柱时，在前驻点O点压强最大，在M点压强最小，压降有利于流动，使流动逐渐加速，OM区称为顺流区，压能转化为动能；从M到F压强增加，速度逐渐减小，不利于流动，MF称逆流区。

对于粘性流体，在上述能量的转化过程中，由于粘性的作用，边界层内的流体质点将要克服粘性力作功而消耗机械能。因此微团在逆压区，不可能到达F点，而是在MF段中的某点S处微团速度降为零，以后来的质点将改道进入主流中，使来流边界层与壁面分离。在分离点下游的区域，受逆压梯度的作用而发生倒流。分离点定义为紧邻壁面顺流区与倒流区的分界点。在分离点附近和分离区，由于边界层厚度大大增加，边界层假设不再成立。边界层分离的必要条件是：逆压梯度和物面粘性的阻滞作用结果。仅有粘性的阻滞作用而无逆压梯度，不会发生边界层的分离，因为无反推力使边界层流体进入到外流区。这说明，顺压梯度的流动不可能发生边界层分离。只有逆压梯度而无粘性的阻滞作用，同样也不会发生分离现象，因为无阻滞作用，运动流体不可能消耗动能而滞止下来。气流绕翼型的流动与边界层分离现象

注意：逆压梯度和壁面粘性阻滞作用是边界层分离的必要条件，但不是充分的，也就是说只有在一定的逆压梯度下，才有可能发生分离。机翼绕流（LES）高攻角机翼边界层分离现象卡门涡街圆柱绕流问题：随着雷诺数的增大边界层首先出现分离，分离点并不断的前移，当雷诺数大到一定程度时，会形成两列几乎稳定的、非对称性的、交替脱落的、旋转方向相反