边界层流动.ppt

,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,第四章 边界层流动,在低雷诺数的爬流流动中，由于粘性力远大于惯性力，因此惯性力项可以从运动方程中略去，从而得到斯托克斯方程。相反，对于高雷诺数的势流流动，由于惯性力远大于粘性力，可以将粘性力忽略，从而得到欧拉方程。,但欧拉方程只适用于离开壁面一定距离的流动主体，并不适用于固体壁面附近,。,为什么在势流流动中，在壁面附近不能忽略粘性力的影响？如何正确处理壁面附近大雷诺数的流体流动问题呢？,这个问题可以由,1904,年普朗特（,Prandtl,）提出的边界层理论来解决。边界层理论阐明了大雷诺数下，粘性力对流体流动的影响。流体在壁面附近的流动也称边界层流动。,问题的引入：,第一节 普兰德边界层理论,一、普兰德边界层理论的要点,流体在壁面附近存在很薄的一个流体层，称为边界层。在边界层中，流体的流速由壁面处为零变为边界层外缘的主体流速，因此，在边界层内垂直于流动方向上的速度梯度很大，剪应力也较大，所以不能忽略粘性力的作用。而在边界层以外的区域，流体的速度梯度则很小，几乎可视为零，此时粘性力可以忽略，可以将其视为理想流体的无旋流动。,可见，对于大雷诺数的流动问题，可以将整个流动区域划分为两个部分。即边界层和外流区。在边界层内，即使流体的粘度很小，但由于速度梯度很大，因此粘性力不能忽略；边界层之外称为外流区。在外流区，壁面对流动的阻滞作用大大削弱，速度梯度极小，故可将粘性力全部略去，从而可以采用欧拉方程去处理。,二、边界层的形成与发展,所谓边界层就是流体速度分布明显受到壁面影响的区域，亦即壁面附近速度梯度较大的薄流体层。,1,、边界层的形成：,u,0,u,0,u,0,A,x,c,平板上边界层的形成,u,0,u,0,u,0,A,x,c,2,、边界层的发展,由层流边界层开始转变为湍流边界层的距离称为临界距离,(,x,c,),。,x,c,的大小与壁面前缘的形状、壁面的粗糙度、流体的性质以及流速等因素有关。壁面愈粗糙、前缘愈钝，则,x,c,愈短。对于平板壁面上的流动，雷诺数的定义为,实验表明，对于光滑的平板壁面，边界层由层流开始转变为湍流的临界雷诺数范围为（,210,5,310,6,）。,u,0,u,0,u,0,层流边界层,湍流边界层,层流内层,A,x,c,平板,上边界层的发展,u,0,u,0,u,0,层流边界层,湍流边界层,层流内层,A,x,c,平板上,当一粘性流体以均匀流速流进水平圆管时，由于流体的粘性作用在管内壁面处形成边界层并逐渐加厚。在距管进口某一段距离，边界层在管中心汇合，此后便占据管的全部截面，边界层厚度即维持不变。,据此可将管内的流动分为两个区域：一是边界层汇合以前的区域，称之为进口段流动；另一是边界层汇合以后的流动，称为充分发展的流动。将入口至边界层汇合处的距离,L,称为进口段长度。,管内流动边界层的形成和发展，与平板边界层相似。如下图所示，,u,u,u,u,u,x,c,d,圆管进口处层流边界层的发展,u,u,u,u,u,x,c,d,圆管进口处层流边界层的发展,若来流速度较小，边界层在管中心汇合时流动为层流，则管内流动继续保持层流，即维持充分发展的层流流动；若来流速度较大，则在进口段内首先形成层流边界层，然后逐渐过渡到湍流边界层，再在管中心汇合后形成充分发展的湍流，如下图所示。,层流时,湍流时,在管内流动充分发展后，流体的流动型态将不再随流动距离,x,变化，此时以,x,定义的雷诺数已不再具有湍动程度的表征意义。因此对于充分发展的管内流动，判别流动型态的雷诺数定义为,式中，,d,为管内径；,u,m,为流体在管内的平均流速或主体流速。,当,Re 4000,时，管内流动为湍流。对湍流流动，进口段长度计算尚无可靠的公式，一般可用下式估计,由于湍流时边界层厚度增长较快，所以其进口段要比层流时短。近似计算时，通常取,L,e,=50,d,。,通常将流体速度为来流速度,99%,时的流层距壁面的法向距离定义为边界层的厚度，以,表示。用公式可表示为,三、边界层的厚度,y,为垂直于流动方向上的距离,边界层厚度随流体的性质（如密度与粘度）、来流速度以及流动距离而变化。在板的前缘处，,=0,；随着距离的增加，边界层逐渐增厚。,对于管内流动，在边界层未汇合以前，边界层厚度的定义和影响因素与平板壁面相同。但流动充分发展后，边界层厚度为管的内半径，即,通常，边界层厚度约在,10,-3,m,的量级。,四、边界层的基本特征,实验研究表明，对于大雷诺数下的流体流动，边界层具有以下两个基本特征：,（,2,）边界层内的粘性力与惯性力量级相同。这是因为边界层内速度梯度很大，即使流体的粘度很小，但作为速度梯度与粘度的乘积,粘性力仍然不可忽略。,（,1,）边界层厚度,要比流场流动的特征尺寸,L,小的多，即,L,。,第二节 普兰德边界层方程,一、普兰德边界层方程的推导,为了简单起见，在此仅考察不可压缩流体在无限大平板壁面上作稳态流动的情形。,u,0,u,0,u,0,层流边界层,湍流边界层,层流内层,A,x,c,平板,上边界层流动,u,0,u,0,u,0,层流边界层,湍流边界层,层流内层,A,x,c,平板上,x,y,假设流体自平板前缘至临界距离,x,c,内所形成的边界层为二维层流流动。以流动方向为,x,方向，以与壁面相垂直的方向为,y,方向。,上一章讲过对于二维的平面流动，连续性方程可以简化为,（,4-1,）,上面的两个运动方程即为普兰德边界层方程。,（,4-2,）,（,4-3,）,运动方程可以简化为,方程,(4-1),、,(4-2),和,(4-3),构成了一个二阶非线性偏微分方程组，方程组中有,3,个方程，,3,个未知数。从理论上来讲方程是可以求解的，但实际上，由于方程的非线性及原函数的复杂性，方程不经简化实际上无法直接求解。,二、普兰德边界层方程的简化,因为方程的边界条件中不出现压力项，所以可以采用以动压力梯度来表示的运动方程：,（,4-2a,）,（,4-3a,）,方程的边界条件：,下面根据边界层流动的特征，采用数量级分析（简称量阶分析）的方法对普兰德边界层方程进一步简化。在进行量级分析之前，首先作两点说明：,（,1,）,数量级分析需要预先选取标准量阶，其它物理量的量阶都是相对标准量阶而言的，当标准量阶改变以后，其它物理量的量阶也随之改变。,（,2,）,所谓量阶不是指该物理量的具体数值，而是该物理量在整个区域内相对于标准量阶的数量级。,在对边界层流动的分析中，选取如下两个标准量阶：,取流动距离,x,作为距离的标准量阶，以来流速度,u,0,作为速度的标准量阶，用符号,O,来表示，写成,O,(,x,)=1,，,O,(,u,0,)=1,，这也意味着这两个物理量的量阶相当。,取边界层厚度,作为另一个标准量阶，由于,很小，故以符号,O,(,)=,来表示。显然，标准量阶,与另外一个标准量阶,1,不在一个水平上，通常,1,是,的,10,3,倍。,当选择了标准量阶以后，可以将其它物理量的量阶与标准量阶相比较,（,4,）,y,。由于在边界层范围内，,y,由壁面处的零值变化至边界层外缘处的,，故,y,的量阶为,y,=,O,(,),。,（,7,）。,（,1,）,u,x,。,u,x,由壁面处的零值变化至边界层外缘处的,u,0,，故其量阶与,u,0,或,x,的量阶相同，即,O,(,u,x,)=1,。,（,2,）。将 写成差分形式，即,（,3,）。,（,5,）,u,y,。由不可压缩流体的二维连续性方程 可知，由于 的量阶为,O,(1),，故 的量阶亦必为,O,(1),，所以,u,y,的量阶是,O,(,),。,（,6,）。,将以上各式代入式（,4-4,），并进行量阶比较,(1)(1)(,)(1/,)(1)(1/,2,),通过量阶比较可知，上式右侧括号内第一项 的量阶远远小于第二项 的量阶，故可将第一项从方程中消去。,由于,(4-10),左侧两个惯性力的量阶均为,O,(1),，而在,边界层内粘性力和惯性力同阶，故右侧粘性力项,故这时方程可以简化为,（,4-10,）,(1)(,),(,)(1)(,2,)(,)(1/,),根据量阶分析可知，等号右边中括号中第一项的量阶远远小于第二项的量阶，因此可以将第一项忽略，这样,y,方向上的普兰德边界层方程可简化为,(1)(,),(,)(1)(,2,)(1/,),（,4-11,）,下面在来看一下,y,方向上的普兰德边界层方程中各项的量阶,由（,4-11,）式可知，动压力梯度项的量阶不可能超过,阶，即,另外，由于,y,方向上的普兰德边界层方程（,4-11,）中每一项的量阶均为,O,(,),阶，而,x,方向上普兰德边界层方程（,4-10,）中的每一项量阶均为,O,(1),阶。由此可见，,y,方向上的普兰德边界层方程可以忽略。,因此,x,方向上的运动方程（,4-10,）中，动压力梯度项就可以从方程中略去,（,4-10,）,(1)(1),2,(1),（,4-10,a,）,这样，普兰德边界层方程经简化后只剩下了,x,方向上的运动方程,但仅（,4-10,a,）一个方程解不出两个未知数,u,x,和,u,y,。,？,这时可以考虑将（,4-10,a,）与连续性方程（,4-1,）联立求解。即，,（,4-10,a,）,x,方向上简化的普兰德边界层方程,（,4-1,）,连续性方程,这是一个二阶非线性偏微分方程组，含有两个因变量,(,u,x,和,u,y,),和两个因变量,(,x,和,y,),，求解起来比较困难，因此可以考虑利用流函数,可以将其化为一个偏微分方程。,（,4-10,a,）,三、普兰德边界层方程的求解,根据流函数,的定义，,将其带入式（,4-10,a,）中，有,由于流函数,自动满足连续性方程，因此（,4-1,）就已经隐含在式（,4-12,）中了。,（,4-12,）,这样由式（,4-1,）和式（,4-10,）构成的二阶非线性偏微分方程组就简化为一个三阶非线性偏微分方程。,利用流函数的概念虽然将由,(4-1),式和,(4-10),式构成的二阶非线性偏微分方程组简化为一个三阶非线性偏微分方程,(4-12),，但要单纯利用数学方法求该方程仍然是非常困难的。,方程的边界条件：,为此，需要通过相似变换的方法将偏微分方程进一步化简为常微分方程。,下面简要介绍一下相似变换法的求解思路。,相似变换的基本思想是：,找到一个无因次位置变量,，使之与,x,和,y,两个自变量同时关联起来。这样就将,与,x,y,之间的函数关系表示为,与,之间的函数关系,f,这样就可以把,关于两个自变量,x,，,y,的偏微分方程转变成关于一个自变量,的常微分方程。,f,x,y,其中,令,这样，,的各阶导数为,（,4,13,）,（,4,14,）,（,4,15,）,（,4,16,）,（,4,17,）,将,的各阶导数带入（,4,12,），并化简得,（,4,18,）,相应的边界条件化为,这样三阶非线性偏微分方程（,4,12,）就化为了三阶非线性常微分方程（,4,18,），该方程虽然从形式上看十分简单，但由于方程的非线性，仍无法得到封闭形式的解析解。,布拉修斯采用级数衔接法近似地求出了式（,4-18,）的解，其后又许多研究者采用数值积分的方法求出了该方程的数值。在此仅给出数值积分的结果，具体求解过程参见教材,83,页。,Blasius,方程,由此可解得不同的,值所对应的,f,、,f,和,f,值，也就得到了各处的速度分布。,f,f,f,f,f,f,0,0.2,0.4,0.6,0.8,1.0,1.2,1.4,1.6,1.8,2.0,2.2,2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0,4.2,4.4,0,0.00664,0.02656,0.05974,0.10611,0.16557,0.23795,0.32298,0.42032,0.52952,0.65003,0.78120,0.92230,1.07252,1.23.99,1.39682,1.56911,1.74696,1.92954,2.11605,2.30576,2.49806,2.69238,0,0.06641,0.13277,0.19894,0.26471,0.32979,0.39378,0.45627,0.51676,0.57477,0.62977,0.68132,0.72899,0.77246,0.81152,0.84605,0.87609,0.90177,0.92333,0.94112,0.95552,0.96696,0.97587,0.33206,0.33199,0.33147,0.33008,0.32739,0.32301,0.31659,0.30787,0.29667,0.28293,0.26675,0.24835,0.22809,0.20646,0.18401,0.16136,0.13913,0.11788,0.09809,0.08013,0.06424,0.05025,0.03897,4.6,4.8,5.0,5.2,5.4,5.6,5.8,6.0,6.2,6.4,6.6,6.8,7.0,7.2,7.4,7.6,7.8,8.0,8.2,8.4,8.6,8.8,2.88826,3.08534,3.28329,3.48189,3.68094,3.88031,4.07990,4.27964,4.47948,4.67938,4.87931,5.07928,5.27926,5.47925,5.67924,5.87924,6.07923,6.27923,6.47923,6.67923,6.87923,7.07923,0.98269,0.98779,0.99115,0.99425,0.99616,0.99748,0.99838,0.99898,0.99937,0.99961,0.99977,0.99987,0.99992,0.99996,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000,1.00000,0.02948,0.02187,0.01591,0.01134,0.00793,0.00543,0.00365,0.00240,0.00155,0.00098,0.00061,0.00037,0.00022,0.00013,0.00007,0.00004,0.00002,0.00001,0.00001,0.00000,0.00000,0.00000,为应用方便，将上式各对应值列成表格形式，如下表所示。,四、普兰德边界层方程的应用,2,、求边界层的厚度,由边界层厚度的定义可知，当时,u,x,/,u,0,=0.99,时，壁面的法向距离,y,即为边界层厚度,。参见上表，当时,u,x,/,u,0,=,f,=,0.99115,时，,=5.0,。,所以有,将上式写成无因次形式,（,4-21,）,1,、求边界层速度,将流函数的定义式带入（,4-13,）和（,4-16,）式得,（,4-19,）,（,4-20,）,由此可得,3,、求流动阻力,流体的流动阻力来自于壁面的剪应力，根据牛顿粘性定律，壁面剪应力,而,根据上表中的数据，有,根据阻力系数的定义，可得距平板前缘,x,处的局部摩擦阻力系数为,于是有,（,4-22,）,（,4-23,）,将壁面剪应力的公式（,4-22,）带入，并积分得,平均阻力系数,C,D,为,上述结果称为布拉修斯解。其在层流范围内与实验数据吻合得很好，但在平板前沿处不成立。,?,不成立的原因是量级关系,x,在该处不成立。,（,4-24,）,（,4-25,）,当流体在一宽度为,b,、长度为,L,的平板壁面上流过时，流体对板面施加的总曳力,F,d,(,主要由摩擦曳力构成,),可表示为,在平板前缘处，平板边界层的阻力系数需要利用高阶边界层理论加以修正，我国力学家郭永怀研究得出的修正公式为：,该公式的适用范围：,5,Re,L,100,。,例：,25,的空气在常压下以,6m/s,的速度流过一平板壁面。试求距平壁前缘,0.15m,处的边界层厚度，并计算在该处距平壁壁面,1mm,处的,u,x,，,u,y,及速度梯度,u,x,/,y,。空气的运动粘度为,1.5510,-5,m,2,/s,，空气的密度为,1.185kg/m,3,。,解：首先计算一下距平板前缘,0.15m,处的雷诺数，以确定流型,（,1,）计算边界层的厚度,由（,4-21,）得,（,2,）计算距平壁壁面,1mm,处的,u,x,，,u,y,及速度梯度,u,x,/,y,。,先求,于是，,由此可见，,u,y,u,x,，因此,y,方向上的流动可以忽略。,查表得，当,=1.6,时，,f,=0.42,f,=0.516,f,=0.296,作业：,p164 2;,第三节 卡门边界层动量积分方程,描述边界层流动的普朗特方程虽然比奈维,-,斯托克斯方程简单，但由于方程的非线性及原函数的复杂性，使得求解过程非常复杂，并且只适用于少数几种简单的流动情形。工程中遇到问题大多是很复杂的，直接求解普兰德边界层方程相当困难，为此人们不得不采用各种近似求解的方法。,一、卡门边界层动量积分方程的推导,为简单起见，本节以不可压缩流体沿壁面作稳态流动为例进行讨论。,冯,卡门根据动量守恒定律和边界层的基本特点，避开奈维,-,斯托克斯方程，直接对边界层进行动量微分衡算，并在此基础上建立了边界层动量积分方程。,如右图所示，密度为,、粘度为,的不可压缩流体在光滑壁面上流动，设边界层外的来流速度为,u,0,，距平板前缘位置,x,处的边界层厚度为,。,在板的宽度方向取单位厚度（,z,1,）。在距壁面前缘,x,处，取微元控制体,ABCDA,。,将动量守恒定律应用于此微元控制体，有：,如果只考虑,x,方向上的受力情况，则有,（,4-24,）,（,4-25,）,下面逐一考察微元体的,4,个面的动量变化情况。,（,1,）,AB,截面：在沿壁面的法向距离,y,处，取微分高度,d,y,，则通过微元截面,d,y,1,流入的质量流率为 ，而通过该微元截面流入的动量流率为 。因此，流入整个截面质量流率和动量流率分别为：,（,4-26,）,（,4-27,）,（,2,）,CD,截面：从,CD,截面（,x,+d,x,处）流出的质量流率和动量流率可由在,AB,截面（,x,处）上的质量流率和动量流率的一阶泰勒展开得到。,（,4-28,）,（,4-29,）,（,3,）,AD,截面：因为是固体壁面，所以,AD,截面上不存在流体质量和动量的流入与流出。,（,4,）,BC,截面：根据质量守恒定律，在稳态下由此截面流入的质量流率应为,CD,截面流出的与,AB,截面流入的质量流率之差，即,由于该截面取在边界层外缘处，故此处的流体均以速度,u,0,流入控制体内，于是从该截面流入的动量流率为,（,4-30,）,（,4-31,）,微元控制体内的净动量变化速率为：“流入”“流出”。即，,（,4-32,）,F,x,=?,作用于控制体上的力有壁面摩擦力和流体静压力，流体微元在各个面上的受力分析如下：,（,2,）,AB,面：由于,AB,面与流动方向垂直，因此只受到左侧流体施加的正压力，而不受摩擦剪应力作用，压力的大小为,（,3,）,CD,面：同样，由于,CD,面垂直于流动方向，因此也只受到右侧流体施加的正压力，力的大小为,“,”,代表力的方向与流动方向相反,（,1,）,AD,面：由于是壁面，,AD,面仅受摩擦剪应力的作用，即,（,4,）,BC,面：因该截面与边界层以外的流体没有速度梯度，剪应力为零，因此也仅受到周围流体的正压力作用，作用力的大小为：,“”,代表力的方向与流动方向相反,所以流体微元受到的合外力为,（,4-33,）,将 代入（,4-34,）式得,（,4-34,）,通过比较各项的量阶，上式可以简化为,（,4-35,）,上式即为卡门边界层动量积分方程。在该方程的推导过程中并没有规定边界层内流体流动的型态，故无论对于层流边界层还是湍流边界层均适用。但求解时要分别代入层流,分布,或湍流速度,分布,方程。此外，该方程也可用于曲面物体边界层。,由于在边界层内,(1)-(1)(1)(,),3,2,(1/,),(1),二、卡门边界层动量积分方程的求解,从式（,4-35,）中可以看出，只要将速度分布，亦即,u,x,y,之间的函数关系式带入方程中，然后积分就可以得到卡门边界层动量积分方程的解析解。,从理论上讲，边界层的速度分布可以通过求解边界层的运动方程和连续性方程得到，但这样问题又回到了出发点，即普兰德边界层方程的求解问题。,为了避开这一难题，可以预先假定一个速度分布，将其带入（,4-35,）中进行求解，然后在将其结果与实验相比较。如果二者吻合，说明所假定的速度分布是正确的，这种方法求得的近似解称为实验逼近解。,边界层内的速度分布可以用,n,次多项式来逼近。多项式的级数越多，就越接近原函数,。,实际上，对于层流边界层，用四阶多项式来逼近边界层的速度分布就已经很接近了。于是，可假定边界层内的速度分布为,（,4-36,）,壁面上流体不滑脱,将上述边界条件带入（,4-36,）可以求得多项式各项的系数,于是，边界层内的速度分布可表示为,（,4-37,）,边界条件：,将速度分布方程式（,4-37,）带入边界层动量积分方程式（,4-35,）,（,4-38,）,（,4-39,）,于是有：,（,4-40,）,右侧求微分得,:,（,4-35,）,左侧求积分得,上式是一个一阶常微分方程，对上式进行积分，并将边界条件,x,=0,，,=0,带入得，,（,4-41,）,写成无因次形式为,（,4-42,）,分离变量得,（,4-21,）,与求解普兰德边界层方程得到的精确解,相比，可见二者相当接近,二、卡门边界层方程的应用,求摩擦阻力和阻力系数,摩擦阻力来自于壁面剪应力,wx,将,（,4-39),式和,（,4-41,）,式带入牛顿粘性定律得,距平板前缘,x,位置处的局部摩擦阻力系数,C,Dx,为,流体流过长度为,L,、宽度为,b,的平板壁面所受的总阻力为,所以，平均阻力系数为,其中，,注意：上述公式仅适用于流体在边界层作层流流动,?,例,:,常压下温度为,20,的空气以的,5m/s,的流速流过一块宽,1m,的平板壁面。试计算距平板前缘,0.5m,处的边界层厚度及进入边界层的质量流率，并计算该段平板壁面的阻力系数和承受的摩擦曳力。设临界雷诺数,Re,xc,=510,5,。已知空气在,20,和,1atm,下的密度和粘度分别为,1.205 kg/m,3,和,1.8110,-5,Pa,s,。,解：（,1,）判断,x,=0.5m,处边界层的流型,所以边界层流动为层流,（,2,）求边界层厚度,（,3,）计算进入边界层的质量流率,在任意位置,x,处，进入边界层内的质量流率可根据下式求出,b,为平板的宽度；,u,x,为距平板垂直距离,y,处空气的流速，将层流边界层的速度分布式代入上式并积分得,代入已知数据，得,（,4,）求阻力系数及曳力,第四节,边界层分离,边界层分离定义：当粘性流体流过曲面物体，如果物体表面曲率比较大（比如圆柱体），常常会出现,边界层与固体壁面相脱离的现象，这种现象称为边界层分离,。此时，壁面附近的流体将发生倒流并产生旋涡，导致流体出现较大的能量损失，它是粘性流体流动时能量损失的重要原因之一。,现以粘性流体绕过一长圆柱体的流动为例讨论边界层是如何形成、发展与分离的。,应予指出，由于产生边界层分离，在分离点以后的流动状况将大大改观，随之而来的是压力分布的改变，它又会使边界层分离的条件发生变化。换而言之，最终的分离点的位置将取决于最终的压力分布和速度分布，而不是取决于初始的流动条件。,由于边界层分离时，在尾流区会产生大量的旋涡，因此会消耗大量的能量。,边界层分离还是产生形体曳力的主要原因。分离点越靠前，形体曳力越大。在多数情况下，像由圆柱体这类形状凸起的物体所产生的总曳力主要是形体曳力。只有在,Re,较低时，因物体表面剪应力引起的摩擦曳力才显得重要。,工程上为了减小边界层分离带来的能量的损失，常常将物体做成流线形，如飞机的机翼、轮船的船体等。流线型的物体可以使边界层的分离点大大延后。,边界层发生分离的判据：,化工流体输送过程中，流体经管件、阀门、管路突然扩大与突然缩小以管路的进出口等局部地方，由于流向的改变和流道截面的突然变化等原因，都会出现边界层的分离现象。目前，对于因边界层分离产生的形体曳力的计算，完全依靠理论求解是很困难的，主要依靠经验方法。,作业：,p94:4,5;,