不可压缩流体二维边界层概述课件讲解

,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,第五章 不可压缩流体二维边界层概述,第一节 边界层的基本概念,第二节 边界层的动量积分方程,第三节 曲面边界层分离现象,卡门涡街,第四节 绕流阻力和阻力系数,在本世纪初之前，流体力学的研究分为两个分支：一是研究流体运动时不考虑黏性，运用数学工具分析流体的运动规律。另一个是不用数学理论而完全建立在实验基础上对流体运动进行研究，解决了技术发展中许多重要问题，但其结果常受实验条件限制。这两个分支的研究方法完全不同，这种理论和实验分离的现象持续了150多年，直到本世纪初普朗特提出了边界层理论为止。由于边界层理论具有广泛的理论和实用意义，因此得到了迅速发展，成为黏性流体动力学的一个重要领域。本章介绍边界层的基本概念及研究方法,第一节 边界层的基本概念,一、边界层的概念,1904年，在德国举行的第三届国际数学家学会上，德国著名的力学家普朗特第一次提出了边界层的概念。他认为对于水和空气等黏度很小的流体，在大雷诺数下绕物体流动时，黏性对流动的影响仅限于紧贴物体壁面的薄层中，而在这一薄层外黏性影响很小，完全可以忽略不计，这一薄层称为边界层。普朗特的这一理论，在流体力学的发展史上有划时代的意义。,图5-1,所示为大雷诺数下黏性流体绕流翼型的二维流动，根据普朗特边界层理论，把大雷诺数下均匀绕流物体表面的流场划分为三个区域，即边界层、外部势流和尾涡区。,图5-1 翼型上的边界层,III外部势流,II尾部流区域,I边界层,边界层外边界,边界层外边界,在边界层和尾涡区内，黏性力作用显著，黏性力和惯性力有相同的数量级，属于黏性流体的有旋流动区；在边界层和尾涡区外，流体的运动速度几乎相同，速度梯度很小，边界层外部的流动不受固体壁面的影响，即使黏度较大的流体，黏性力也很小，主要是惯性力。所以可将这个区域看作是理想流体势流区，可以利用前面介绍的势流理论和理想流体伯努里方程来研究流场的速度分布。普朗特边界层理论开辟了用理想流体理论和黏性流体理论联合研究的一条新途径。实际上边界层内、外区域并没有明显的分界面，一般将壁面流速为零与流速达到来流速度的99处之间的距离定义为边界层厚度。边界层厚度沿着流体流动方向逐渐增厚，这是由于边界层中流体质点受到摩擦阻力的作用，沿着流体流动方向速度逐渐减小，因此，只有离壁面逐渐远些，也就是边界层厚度逐渐大些才能达到来流速度。,根据实验结果可知，同管流一样，边界层内也存在着层流和紊流两种流动状态，若全部边界层内部都是层流，称为层流边界层，若在边界层起始部分内是层流，而在其余部分内是紊流，称为混合边界层，如图5-2所示，在层流变为紊流之间有一过渡区。在紊流边界层内紧靠壁面处也有一层极薄的层流底层。判别边界层的层流和紊流的准则数仍为雷诺数，但雷诺数中的特征尺寸用离前缘点的距离x表示之，特征速度取边界层外边界上的速度 ，即,(5-1),图5-2 平板上的混合边界层,层流边界层,过渡区域,紊流边界层,层流底层,对平板的边界层，层流转变为紊流的临界雷诺数为 。临界雷诺数的大小与物体壁面的粗糙度、层外流体的紊流度等因素有关。增加壁面粗糙度或层外流体的紊流度都会降低临界雷诺数的数值，使层流边界层提前转变为紊流边界层。,二、边界层的基本特征,(1)与物体的特征长度相比，边界层的厚度很小,.,(2)边界层内沿厚度方向，存在很大的速度梯度。,(3)边界层厚度沿流体流动方向是增加的，由于边界层内流体质点受到黏性力的作用，流动速度降低，所以要达到外部势流速度，边界层厚度必然逐渐增加。,(4)由于边界层很薄，可以近似认为边界层中各截面上的 压强等于同一截面上边界层外边界上的压强值。,(5)在边界层内，黏性力与惯性力同一数量级。,(6)边界层内的流态，也有层流和紊流两种流态。,第二节 边界层的动量积分方程,边界层内的流体是黏性流体的运动，理论上可以用N-S方程来研究其运动规律。但由此得到的边界层微分方程中，非线性项仍存在，因此即使对于外形很简单的绕流物体求解也是很复杂的，目前只能对平板、楔形体绕流层流边界层进行理论计算求得其解析解。但工程上遇到的很多问题，如任意翼型的绕流问题和紊流边界层，一般来说求解比较困难，为此人们常采用近似解法，其中应用的较为广泛的是边界层动量积分方程解法。,下面来推导边界层动量积分方程。假定平面边界内流动是定常的并忽略质量力，在边界层的任一处，取单位宽度、沿边界层长度为,d,的微元段作为控制体，如,图,5-3,所示。控制体的控制面由边界层的横断面,AB,与,CD,以及内边界,AD,和外边界,BC,组成。对控制体应用物理概念十分清楚的动量方程则有：通过控制面,AB,、,BC,、,CD,的动量变化率等于作用在控制面,AB,、,BC,、,CD,、,AD,上所有外力的合力。,首先计算通过边界层控制面在轴方向上的动量变化率。,单位时间流入x处控制面,AB,的动量为,从 处控制面,CD,流出的动量为,从控制面,BC,流入的动量采用下列求法，首先计算从 处控制面,AB,流入的质量流量,而从 处控制面,CD,流出的质量流量为,由不可压缩流体的连续性方程可知，通过,CD,与,AB,控制面质量流量的差值应等于由,BC,控制面流入的质量流量，于是流入,BC,控制面的质量流量与动量分别为,图5-3 推导边界层的动量积分关系式用图,整理上述单位时间内通过控制面的流体动量的通量在x方向的分量，得,下面计算作用在控制面上所有外力在x轴方向的合力。忽略质量力，故只有表面力。,作用在控制面AD上的表面力为,作用在控制面AB、CD上的表面力分别为,作用在边界层外边界控制面BC上的表面力，因摩擦应力为零，而压强可取B、C两点压强的平均值，于是有,整理上述作用在控制面上的所有表面力在x方向的代数和，并注意到略去二阶小量，得,式(5-2)又称为边界层动量积分关系式。该式是匈牙利科学家冯卡门(VonKarman),于1921年根据边界层的动量定理首先推导出来的。由于在推导过程中未加任何近似条件，从这个意义上讲，它是严格的，而且对边界层的流动性质也未加限制，因此它既可求解层流边界层，又可适用于紊流边界层。,根据动量定理，令 ，可得边界层动量积分方程为,(5-2),由于积分上限 只是 的函数，因此式(5-2)中 的可写成,又根据势流的伯努里方程,注意到式（5-3），则式（5-2）可写成,常数,则有,（5-3）,（5-4）,考察边界层的动量积分方程式（5-2）和式（5-4）可以看到，方程中含有五个未知量,：,、，其中 和 可由主流区的势流方程求得，剩下的三个未知量是 、，因此要求解边界层动量积分方程，原则上还需要补充两个方程，即,(1)满足绕流物体壁面条件和边界层外边界条件的速度分布 ；,(2)与速度分布有关的 与 的关系式。事实上，与 的关系可根据边界层内的速度分布求出,通常在求解边界层动量积分方程时，总是先选取边界层内速度分布，选取的速度分布 越接近实际，则所得结果越正确。但由于边界层运动的复杂性，而预先选定的速度分布只能满足主要的边界条件，不可能正好满足动量积分方程，这样求得的结果（、等）就都是近似的，故积分方程的解法只能是近似的解法。但这种解法有一个很大的优点，就是只要能大致选定速度分布形式，则可以得到误差并不很大的结果，而且解法较简单，因此在工程上用得较广泛。,下面列出了用,动量积分方程求得的,平板层流和紊流边界层的部分近似解。,对于层流边界层,平板上离前缘点处的边界层厚度,(5-5),在平板一个壁面上由粘滞力引起的总摩擦阻力,(5-6),摩擦阻力系数,(5-7),对于紊流边界层,平板上离前缘点处的边界层厚度,(5-8),在平板一个壁面上由粘滞力引起的总摩擦阻力,(5-9),摩擦阻力系数,(5-10),以上几式中 均匀来流速度，,m/s；,平板的宽度，,m；,平板的长度，,m；,来流的密度，,kg,/,m,3,。,第三节 曲面边界层分离现象,卡门涡街,如前所述，当不可压缩黏性流体纵向流过平板时，在边界层外边界上沿平板方向的速度是相同的，而且整个流场和边界层内的压强都保持不变。当黏性流体流经曲面物体时，边界层外边界上沿曲面方向的速度是改变的，所以曲面边界层内的压强也将同样发生变化，对边界层内的流动将产生影响。曲面边界层的计算是很复杂的，这里不准备讨论它。这一节将着重说明曲面边界层的分离现象。,一、曲面边界层的分离现象,在实际工程中，物体的边界往往是曲面（流线型或非流线型物体）。当流体绕流非流线型物体时，一般会出现下列现象：物面上的边界层在某个位置开始脱离物面，并在物面附近出现与主流方向相反的回流，流体力学中称这种现象为边界层分离现象，如,图5-4,所示。流线型物体在非正常情况下也能发生边界层分离，如图5-4(a)所示。,（a）流线形物体；（b）非流线形物体,图5-4 曲面边界层分离现象示意图,边界层,外部流动,外部流动,尾迹,外部流动,外部流动,尾迹,边界层,现以不可压缩流体绕流圆柱体为例，着重从边界层内流动的物理过程说明曲面边界层的分离现象。当黏性流体绕圆柱体流动时，在圆柱体前驻点,A,处，流速为零，该处尚未形成边界层，即边界层厚度为零。,随着流体沿圆柱体表面上下两侧绕流，边界层厚度逐渐增大。层外的流体可近似地作为理想流体，理想流体绕流圆柱体时，在圆柱体前半部速度逐渐增加，压强逐渐减小，是加速流。当流到圆柱体最高点,B,时速度最大，压强最小。到圆柱体的后半部速度逐渐减小，压强逐渐增加，形成减速流。由于边界层内各截面上的压强近似地等于同一截面上边界层外边界上的流体压强，所以，在圆柱体前半部边界层内的流动是降压加速，而在圆柱体后半部边界层内的流动是升压减速。,因此，在边界层内的流体质点除了受到摩擦阻力的作用外，还受到流动方向上压强差的作用。在圆柱体前半部边界层内的流体质点受到摩擦阻滞逐渐减速，不断消耗动能。但由于压强沿流动方向逐渐降低，使流体质点得到部分增速，也就是说流体的部分压强能转变为动能，从而抵消一部分因摩擦阻滞作用而消耗的动能，以维持流体在边界层内继续向前流动。,但当流体绕过圆柱体最高点,B,流到后半部时，压强增加，速度减小，更促使边界层内流体质点的减速，从而使动能消耗更大。当达到,S,点时，近壁处流体质点的动能已被消耗完尽，流体质点不能再继续向前运动，于是一部分流体质点在,S,点停滞下来，过,S,点以后，压强继续增加，在压强差的作用下，除了壁上的流体质点速度仍等于零外，近壁处的流体质点开始倒退。,接踵而来的流体质点在近壁处都同样被迫停滞和倒退，以致越来越多被阻滞的流体在短时间内在圆柱体表面和主流之间堆积起来，使边界层剧烈增厚，边界层内流体质点的倒流迅速扩展，而边界层外的主流继续向前流动，这样在这个区域内以,ST,线为界，如图,5-5,所示，在,ST,线内是倒流，在,ST,线外是向前的主流，两者流动方向相反，从而形成旋涡。,图5-5 曲面边界层分离现象,使流体不再贴着圆柱体表面流动，而从表面分,图5-5,曲面边界层分离现象离出来，造成边界层分离，,S,点称为分离点。形成的旋涡，不断地被主流带走，在圆柱体后面产生一个尾涡区。尾涡区内的旋涡不断地消耗有用的机械能，使该区中的压强降低，即小于圆柱体前和尾涡区外面的压强，从而在圆柱体前后产生了压强差，形成了压差阻力。压差阻力的大小与物体的形状有很大关系，所以又称为形状阻力。,二、卡门涡街,1911,年，匈牙利科学家卡门在德国专门研究了这种圆柱背后旋涡的运动规律。实验研究表明，当时黏性流体绕过圆柱体，发生边界层分离，在圆柱体后面产生一对不稳定的旋转方向相反的对称旋涡，超过,40,后，对称旋涡不断增长，至时，这对不稳定的对称旋涡，最后形成几乎稳定的非对称性的、多少有些规则的、旋转方向相反、上下交替脱落的旋涡，这种旋涡具有一定的脱落频率，称为卡门涡街，如图,5-6,所示。,图5-6 卡门涡街形成示意图,圆柱体的卡门涡街的脱落频率 与流体流动的速度 和圆柱体直径 有关，由泰勒(FTaylor)和瑞利(LRayleigh)提出下列经验公式,