流体力学-物体绕流流动分析课件.ppt

1、2022-6-8流体绕流流动1内蒙古工业大学工程流体力学电子课件内蒙古工业大学工程流体力学电子课件10物体绕流流动物体绕流流动10.1 边界层理论及基本概念10.1.1 边界层理论本章讨论大雷诺数情形下的流动问题,着重介绍普朗特的边界层理论。自1904年普朗特创立边界层理论以来,由于它的应用范围极为广泛,发展非常迅速,早已成为粘性流体力学的主要发展方向之一。边界层学说还与传热过程和传质过程有密切关系。边界层理论的主要任务是研究物体在流体中运动时所受到的摩擦阻力和物体与流体间的热交换。2022-6-8流体绕流流动2例例1 1：空气运动粘度：空气运动粘度51 4 102. ms6522 1 52

2、4 101 4 10Vh.Re.大大Re数流动是常见现象数流动是常见现象.设设汽车汽车1 58022h. m , V km h m s例例2 2：水运动粘度：水运动粘度61 102 ms 76108 . 2101108 . 2ReVl设船设船10102 8lm , V km h. m s2022-6-8流体绕流流动3用用N-S方程可以得到小雷诺数流动条件下的近似解，方程可以得到小雷诺数流动条件下的近似解，工程上涉及到大雷诺数流动，要寻求新的近似方法。工程上涉及到大雷诺数流动，要寻求新的近似方法。若采用欧拉方程,同时在固体壁面上采用滑移条件(而不是无滑移粘附条件),这就是理想流体的模型。在理想流

3、体模型的范围内,算出的物体表面的压力分布,在流动不分离或在接近尾缘处有小分离区的情况下与实测结果比较符合。但无法解决阻力问题。在实际流体绕流固体时，固体边界上的流速为在实际流体绕流固体时，固体边界上的流速为0，在，在固体边界的外法线方向上的流体速度从固体边界的外法线方向上的流体速度从0迅速增大，在迅速增大，在边边界附近的流区存在相当大的速度梯度，在这个流区内粘界附近的流区存在相当大的速度梯度，在这个流区内粘性作用不能忽略，性作用不能忽略，边界附近的流区称为边界附近的流区称为边界层（或附面边界层（或附面层），层），边界层外流区，粘性作用可以忽略，当作理想流边界层外流区，粘性作用可以忽略，当作理想

4、流体来处理体来处理。 2022-6-8流体绕流流动4如图，平板前方均匀来流的速度如图，平板前方均匀来流的速度v，从平板前缘开始从平板前缘开始形成边界层，其厚度沿流增加。在边界层外缘附近流形成边界层，其厚度沿流增加。在边界层外缘附近流速渐近于当地外流速度。认为边界层厚度是沿表面法速渐近于当地外流速度。认为边界层厚度是沿表面法线方向从到的一段距离。线方向从到的一段距离。边界层定义：边界层定义：绕流物体表面上一层厚度很小且其中的绕流物体表面上一层厚度很小且其中的流动具有很大法向速度梯度的流动区域。流动具有很大法向速度梯度的流动区域。2022-6-8流体绕流流动5整个流场可以明显地分成性质很不相同的两

5、个区域：(1)(1)紧贴物面非常薄的一层区域称为边界层。在该区域内紧贴物面非常薄的一层区域称为边界层。在该区域内, ,速度速度分量分量u ux x沿物面的法向变化非常迅速沿物面的法向变化非常迅速, ,它比沿切向的变化高一个数它比沿切向的变化高一个数量级。即量级。即 甚大。虽然在大甚大。虽然在大 数情况数情况, , 很小很小, ,但因但因 很大很大, ,故粘性应力故粘性应力 仍然可以达到很高的数值。仍然可以达到很高的数值。2)2)边界层外的整个流动区域称为外部流动区域。在该区域内边界层外的整个流动区域称为外部流动区域。在该区域内 , , 很小很小, ,因此粘性应力因此粘性应力 在大在大 ReRe

6、数情况下的确比惯性数情况下的确比惯性力小得多力小得多, ,可以将粘性力全部略去可以将粘性力全部略去, ,因而把流体近似地看成是理因而把流体近似地看成是理想的。对于均匀来流绕过物体的流动而言想的。对于均匀来流绕过物体的流动而言, ,在整个外部流动区域在整个外部流动区域中不仅可把流体视为理想的中不仅可把流体视为理想的, ,而且可视为运动是无旋的。而且可视为运动是无旋的。 Reyuxyuxyux注意：1. 对于平板绕流，边界层外缘，对于弯曲固壁，边界层外缘。2. 边界层的外边界线与流线不重合，外流区域中的流体质点可以连续地穿过边界层的外缘 进入边界层内。2022-6-8流体绕流流动6一、边界层特点一

7、、边界层特点普朗特理论：边界层内惯性力与粘性力量级相等。普朗特理论：边界层内惯性力与粘性力量级相等。1. 边界层很薄边界层很薄22 yuxuu22 UlUUll 22Re1 l当当610 , 0 001Rel.2. 边界层厚度增长边界层厚度增长Uxxx )(22x( x) U10.1.2 边界层基本概念2022-6-8流体绕流流动7二、边界层厚度3. 边界层内流态边界层内流态 实验测量表明边界层内层流实验测量表明边界层内层流态向湍流态转态向湍流态转捩捩的雷诺数为的雷诺数为53 2 10 xcrRe.1名义厚度名义厚度定义为速度达到外流速度定义为速度达到外流速度99%的厚度。的厚度。对对平板层流

8、边界层平板层流边界层Ux0 . 52022-6-8流体绕流流动82.排挤厚度（位移厚度）排挤厚度（位移厚度）将由于不滑移条件造成的质量亏将由于不滑移条件造成的质量亏损折算成无粘性流体的流量相应损折算成无粘性流体的流量相应的厚度的厚度d 。又称为。又称为质量流量亏损厚度质量流量亏损厚度uouo yuddyuuod013. 动量损失厚度动量损失厚度m将由于不滑移条件造成的动量流将由于不滑移条件造成的动量流量亏损折算成无粘性流体的动量量亏损折算成无粘性流体的动量流量相应的厚度流量相应的厚度m 。dyuuuuom001 动量损失厚度动量损失厚度排挤厚度排挤厚度2022-6-8流体绕流流动910.2 平

9、板边界层流动10.2.1 普郎特边界层方程51011022222222yuxuvypyuuxuuyuxuvxpyuuxuuyuxuyyyyyxxxxyxxyx2022-6-8流体绕流流动10今以来流速度 为特征速度,平板长度 为特征长度。引进如下无量纲量 0uLLyy,Lxxupp,uuu,uuu*20*0y*y0 x*x无量纲的基本方程组： 2*22*2*2*22*2*110yuxuReypyuuxuuyuxuRexpyuuxuuyuxuyyLyyyxxxLxyxxyx2022-6-8流体绕流流动11在大 数情况下边界层内的流动有如下两个主要性质：LRe(1)边界层的厚度较物体的特征长度 L

10、小得多,即 是一个小量。(2)在边界层内粘性力和惯性力同数量级。L*对“量级”的两点说明：(1)估计量级必须有个标准。在边界层问题中,我们取*作为估计量级的标准。并采用符号 ,例如（*）表示和* 同量级。(2)所谓量级不是指该物理量或几何量的具体数值,而是指该量在整个区域内相对于标准小参数(在边界层问题中,即相对于*)而言的平均水平。所以允许(1)的量(即与1同量级的量)甚至比1大一个量级的量在个别点上或局部小区域内取较小的数值,甚至等于零,重要的是它的平均水平是与1同量级,或比1大一个量级就行了。 2022-6-8流体绕流流动121. 及其各阶导数 的量级*xu2*2*2*2*,xuxuyu

11、yuxxxx 1*Oux*1Oyux2*2*21Oyux 1*Oxux 12*2Oxux2. 及其各阶导数 的量级 *yu2*2*2*2*,xuxuyuyuxyyy由连续方程得 1*Oxuyuxy积分得*0*yyyOdyyuu2022-6-8流体绕流流动13*Oxuy*2*2Oxuy*2*21Oyuy3. 及 的量级 *xp*yp压力梯度是起调节作用的被动的力。它们的量级由方程中其它类型力中的最大量级所决定 1*Oxp *Oyp2022-6-8流体绕流流动14*2*2*22*2*2*2*2*22*2*1)()(1)()1()()()1(1)1()(1)1(1)()1()1()1()1(0OOO

12、yuxuReOOOOOypyuuxuuOOOyuxuReOOOOOxpyuuxuuOOyuxuyyLyyyxxxLxyxxyx2022-6-8流体绕流流动15无量纲形式的普朗特层流边界层方程 2*2*1dd0yuRexpyuuxuuyuxuxLxyxxyx2022-6-8流体绕流流动16有量纲形式的普朗特层流边界层方程为： 22dd/10yuvxpyuuxuuyuxuxxyxxyx边界条件 在物面上 , 0y0yxuu 在边界层外边界 上, )(xy099. 0uux2022-6-8流体绕流流动17内蒙古工业大学工程流体力学电子课件内蒙古工业大学工程流体力学电子课件10.2.2 平板层流边界层

13、的精确解因为平板没有厚度,当理想流体沿平板方向流过平板时,平板对流动没有扰动,因此外流的速度场是均匀的且等于常数 。根据伯努利方程,压力也均匀 0u常数, p0ddxp2022-6-8流体绕流流动180yuxuyx22yuvyuuxuuxxyxx普朗特边界层方程式 （10-10）边界条件是： x0， 0:0yxuuy0:uuyx0yp2022-6-8流体绕流流动19说明说明：第三式表明边界层内第三式表明边界层内y方向压强梯度为零，说明外部压强可穿方向压强梯度为零，说明外部压强可穿透边界层直接作用在平板上。外部压强由势流决定透边界层直接作用在平板上。外部压强由势流决定ddddpUUxx 第二式右

14、边得到简化（第二式右边得到简化（x方向二阶偏导数消失），有利于数值方向二阶偏导数消失），有利于数值计算。利用该方程就可计算壁切应力和流动阻力，具有里程碑计算。利用该方程就可计算壁切应力和流动阻力，具有里程碑式意义。式意义。2022-6-8流体绕流流动20无量纲化方程和边界条件为： 2*2*0yuyuuxuuyuxuxxyxxyx1:0:0, 0*xyxuyuuyx2022-6-8流体绕流流动21返回有量纲形式为 Lvuy,Lxuuu0*x0 xLvuyLxuvuLuyy0*0,由于讨论的是半无限长平板上边界层内的流动,不存在特征长度 ,因此在最后的解中不应该出现特征长度 ,即 和 应该与 无关

15、。这就要求自变量以下列组合LLxuyuLvxuyLxLvuyxy00*)(*fux)(1*Fxuy2022-6-8流体绕流流动22*2xx*1xy)(2*fxxuxuxx)(1*Fxyuyuyy)(1*fxyux)(1*2*2fxyux )(2)(fFdfF)(21)(dff)()(21)()(f令)()(21)(F则2022-6-8流体绕流流动2317100)()(21)( 1)(:0)(, 0)(:0这种原有两个自变量的偏微分方程组,若其解只依赖于一个组合变量,能使偏微分方程变为常微分方程,则称此方程式具有相似性解。（） 所满足的方程式是一个非线性的三阶常微分方程,形式虽然简单,但却无法找

16、出封闭形式的解析解来。 布拉修斯于1908年用级数衔接法求出此问题的近似解,而后托普费尔(Topfer)、哥德斯坦(Goldstein)、霍华斯(Howarth)、哈脱利(Hartree)等人分别用数值方法求出了此问题的解。 2022-6-8流体绕流流动24由数值解绘制的无量纲速度廓线由数值解绘制的无量纲速度廓线与尼古拉兹实验测量结果吻合。与尼古拉兹实验测量结果吻合。布拉修斯解边界层厚度21Re96. 4xx排挤厚度21Re73. 1xdx动量损失厚度21Re664. 0 xmx2022-6-8流体绕流流动25布拉修斯的精确解有两个重要的前提条件：边界层内的流动是层流；沿 平板的压力梯度为零。

17、工程中常遇到的问题要复杂些，如边界层是湍流，沿 壁面的压力梯度不为零等等，这就促使人们去寻求近似解法。在近似解法 中应用最广泛的是冯卡门的边界层近似积分法。 10.2.310.2.3 边界层动量方程边界层动量方程2022-6-8流体绕流流动26 流体绕流中作用在物体上的力可以分为垂直于来流体绕流中作用在物体上的力可以分为垂直于来流方向的升力和平行于来流方向的阻力，绕流阻力可流方向的升力和平行于来流方向的阻力，绕流阻力可以分成摩擦阻力与形状阻力，都与边界层有关。绕流以分成摩擦阻力与形状阻力，都与边界层有关。绕流阻力作用表现在于边界层内流速的降低，引起动量的阻力作用表现在于边界层内流速的降低，引起

18、动量的变化。通过建立边界层的动量方程来研究摩擦阻力。变化。通过建立边界层的动量方程来研究摩擦阻力。 沿物体的曲面取沿物体的曲面取x轴，沿物体表面法线取轴，沿物体表面法线取y轴，在物轴，在物体表面取边界层微元段体表面取边界层微元段ABCD，把它放大，把它放大，x轴便成为轴便成为直线，线段直线，线段BD长为长为dx，AC为边界层外边界，为边界层外边界，AB、CD垂直于物体表面。垂直于物体表面。0u0u2022-6-8流体绕流流动27假设：假设： 不计质量力不计质量力 流动为定常流动流动为定常流动 dx无限小，无限小，BD、AC可看成直线可看成直线由动量方程由动量方程 由控制面由控制面AB沿沿x方向

19、流入动量方向流入动量 （1）由控制面由控制面CD沿沿x方向流出动量方向流出动量 由控制面由控制面AC沿沿x方向流入动量方向流入动量 xACABCDFMMM02dyvMxAB)2()(0202dxdyvxdyvdxxMMMxxABABCD)3()(0dxdyvxuMxoAC 4sin)21()(0dxdsdxxppddxxpppFx2022-6-8流体绕流流动28因为因为 ，所以，所以 边界层内边界就是物体表面，其流速为边界层内边界就是物体表面，其流速为0，其压强等于边，其压强等于边界层外边界的压强，即沿物体表面的法线界层外边界的压强，即沿物体表面的法线y方向压强不变，方向压强不变，p与与y 无

20、关，可用全微分代替偏微分，上式可写作无关，可用全微分代替偏微分，上式可写作 （5）将（将（2）、（）、（3）、（）、（4）、（）、（5）代入（）代入（1）得到）得到 （6）方程（方程（6）就是边界层积分方程，由冯）就是边界层积分方程，由冯卡门首先推导出卡门首先推导出来来的，称作卡门动量积分方程。的，称作卡门动量积分方程。ddssindxdxdxpdxxpFx021dxdxdxdpFx00002dxdpdyudxddyudxduxxo2022-6-8流体绕流流动29讨论：（1）如果dp/dx=0，um为常数（2）动量损失厚度020dyuuudxdxxo20omudxd0002dxdpdyudxd

21、dyudxduxxo2022-6-8流体绕流流动3010.2.4 平板层流边界层的近似解绕平板的稳态不可压缩流动：x方向的压力梯度dpdx取边界层外边界上外流(即势流)的只值；对于平扳来说，边界层外边界上的速度等于来流速度uo ；假定一个合理的速度分布通常认为沿y方向的速度分布曲线在任 意x处类似，即认为u是y的函数；根据牛顿剪切定律，用 代替。 0/yyu2022-6-8流体绕流流动311. 速度剖面族的选取我们选取如下三次多项式3322100)()()()(xaxaxaxauux)(xy其中依据条件 , 可得： 00yxu0)(0 xa由条件 ,可得0yxuu1321aaa由条件 ,可得0

22、yxyu032321aaa设定平板上为层流边界层，首先设定平板上为层流边界层，首先补充边界层流速分布关系式，补充边界层流速分布关系式，2022-6-8流体绕流流动320022yxyu由条件 ,可得 02a21,2331aa选定的速度剖面为 302123uux2. 单参数（x）的确定0100000d1d1uuuuyuuuuxxxx因为d2123121233103)(28039x2022-6-8流体绕流流动33xyuuuuxxxdd28039d1dd000所以00200232323uuyuyx又020020023uvuyuuyx故023dd28039uvxxuvd13140d02022-6-8流体

23、绕流流动34102Cuvx131402积分得0 x001C21064. 464. 4xxReuvx（ ）21048. 548. 5xxReuvx3 壁面摩擦阻力xReu200323. 0 xuRex02022-6-8流体绕流流动35xfReuC1646. 021200局部阻力系数）（xfReuC1730. 021200 xfReuC1664. 021200布拉修斯的精确解则对长为 宽为 且两边浸没在流体中的平板所受到的总摩擦阻力为：,Lb10200292. 12LubdxbFd2022-6-8流体绕流流动36LddRebLuFC292. 122120总阻力系数为LddRebLuFC328. 1

24、22120（ ）10.2.5 平板湍流边界层的近似解平板湍流边界层的近似解为了方便,在工程上往往采用幂次公式作为近似速度剖面族。例如采用经验公式(以下表示时均值的上标“一”略去,并用 代替) ：uxu一、湍流边界层内速度剖面的选取一、湍流边界层内速度剖面的选取 nnvyuCuu1*2022-6-8流体绕流流动37vLuL/Re0当 在510 5到107范围内, 随着 ReL增加, 1/n下降 ,711,74. 8nCn1011, 5 .11nCn71*74. 8vyuuu如取0yuu在湍流外边界处710*74. 8vuuu（*）2022-6-8流体绕流流动38710yuu式中（x）是依赖于x的

25、未知函数 =f（x）。为了完全确定速度剖面,还需要确定（x） 。 为了确定 （x） ,需要应用式 二、单参数二、单参数 （x） 的确定的确定 0200001udyuuuudxdxx由（*）式得又因410202*uvu0225. 0u2*uw2022-6-8流体绕流流动39所以410202*2000225. 0uvuuu将以上代入卡门积分动量方程式，得4100225. 0727uvdxdCxuv41045289. 0通过采用近似方法,认为湍流边界层从平板前缘就已形成,即认为：00 x, 0C51Re37. 0)(xxx2022-6-8流体绕流流动40三、摩擦阻力三、摩擦阻力 514120037.

26、 0045. 021xfReuC510577. 0 xRe局部阻力系数长为L,宽为b的平板受到的总阻力系数为： 510Re072. 01LLfddxCLC同实验结果比较表明,若把上式中的0.072改为0.074,即 ）（5110074. 051LdReC7510105LRe适用范围2022-6-8流体绕流流动41当 时,式（10-51）就不太准确了。需要采用对数速度剖面。在 的范围内,史里希丁用对数速度分布与积分关系式联合求解得出的摩擦阻力系数公式为：710LRe971010LRe它表明阻力和来流速度的1.8次方成正比,而层流时阻力和来流速度 的1.5次方成正比,故湍流边界层的摩阻比层流边界层

27、的摩阻大。0u58. 2Relg455. 0LdC 边界层内流动形态转变的典型情况表示在图中。： 2022-6-8流体绕流流动42四、平板混合边界层的近似计算四、平板混合边界层的近似计算前面假定整个平板上是层流或紊流边界层，实际上，前面假定整个平板上是层流或紊流边界层，实际上，当当Re增大到一定数值时，平板长度达到一定长度，即增大到一定数值时，平板长度达到一定长度，即 Lxer时，平板前部是层流边界层，后部是紊流边界层，时，平板前部是层流边界层，后部是紊流边界层，中间有一过渡段，这种边界层称为混合边界层。中间有一过渡段，这种边界层称为混合边界层。计算时引入假设：计算时引入假设：（1）层流边界层

28、转变为紊流边界层是在处突然发生，无）层流边界层转变为紊流边界层是在处突然发生，无 过渡段；过渡段；（2） 混合边界层的紊流边界层可以看作是从平板的首端混合边界层的紊流边界层可以看作是从平板的首端 开始的紊流边界层的一部分开始的紊流边界层的一部分。普朗特建议：在边界层转捩位置 以前采用层流的摩擦阻力系数 ,在其后采用的摩擦阻力系数 ,于是混合边界层的总的阻力系数 为erxflCftCdC2022-6-8流体绕流流动43ererxLxftfddxCdxCLC011erxLftftfdxCdxCCL001)(1LxfterfererftdxCdxCxLxdxCLer000111LdtReACerer

29、erxxxReA5 . 0328. 12 . 074. 0ReRe5 . 08 . 0328. 1074. 0ererxxReRe其中2022-6-8流体绕流流动44近似计算方法2022-6-8流体绕流流动4510.3 边界层分离与压差阻力10.3.1 边界层分离现象 流体绕过非线型钝头物体时，较早脱离物体表面，流体绕过非线型钝头物体时，较早脱离物体表面，在物体后部形成较宽阔的尾流区，在物体后部形成较宽阔的尾流区，在边界层内，流体在边界层内，流体质点在某些情况下向边界层外流动的现象称为边界层质点在某些情况下向边界层外流动的现象称为边界层从固体分离。从固体分离。 以圆柱绕流为例，虚线为边界层外边

30、界。以圆柱绕流为例，虚线为边界层外边界。2022-6-8流体绕流流动46注意：注意：C C点的位置点的位置，这是由于在加速减压和减速增压的过，这是由于在加速减压和减速增压的过程中，还存在克服流动阻力所消耗的能量损失程中，还存在克服流动阻力所消耗的能量损失 由伯努利方程知，愈靠近圆柱，流速越小，压强由伯努利方程知，愈靠近圆柱，流速越小，压强越大，在贴近圆柱面越大，在贴近圆柱面A处流速为处流速为0，压强最大，压强最大，A点称点称为驻点。由于液体不可压缩，继续流来的液体质点在为驻点。由于液体不可压缩，继续流来的液体质点在驻点的压强的作用下，将压能转化为动能，从而改变驻点的压强的作用下，将压能转化为动

31、能，从而改变流向，沿圆柱面两侧继续向前流动。由于圆柱面的阻流向，沿圆柱面两侧继续向前流动。由于圆柱面的阻滞作用，在表面产生边界层，从滞作用，在表面产生边界层，从A点经点经1/4圆周到圆周到B点之点之前，柱面向外凸出，流线趋于密集，边界层内流体处前，柱面向外凸出，流线趋于密集，边界层内流体处在加速减压情况，在加速减压情况， ，这时由于压能减小部分还能够，这时由于压能减小部分还能够补偿动能增加和由于克服流动阻力而消耗的能量损失补偿动能增加和由于克服流动阻力而消耗的能量损失，因此此时因此此时B点处边界层内流体质点速度不为点处边界层内流体质点速度不为0。0 xp2022-6-8流体绕流流动47 过过B

32、 B点之后，流线逐渐疏散，边界层内流体处于减速点之后，流线逐渐疏散，边界层内流体处于减速增压的情况，动能转化成压能，同时也用以克服流动阻增压的情况，动能转化成压能，同时也用以克服流动阻力而消耗的能量。在力而消耗的能量。在C C点处边界层内流体质点速度下降为点处边界层内流体质点速度下降为0 0。流体质点在流体质点在C C点停滞下来，形成新的停滞点，继续流点停滞下来，形成新的停滞点，继续流来的流体质点将脱离原来的流线，沿另一流线来的流体质点将脱离原来的流线，沿另一流线CECE流去，流去，从而使边界层脱离了圆柱面，这样就形成了边界层的分从而使边界层脱离了圆柱面，这样就形成了边界层的分离现象，离现象，

33、C C点为点为分离点分离点。分离点的位置与绕流物的形状、。分离点的位置与绕流物的形状、粗糙程度、流动的粗糙程度、流动的ReRe数和来流与物体的相对方向有关。数和来流与物体的相对方向有关。 边界层分离后，边界层和圆柱面之间，由于分离点下边界层分离后，边界层和圆柱面之间，由于分离点下游压强大，从而使流体发生反向回流，游压强大，从而使流体发生反向回流，形成旋涡区形成旋涡区。2022-6-8流体绕流流动48边界层分离：边界层脱离壁面边界层分离：边界层脱离壁面 边界层分离边界层分离2.2.分离的原因分离的原因 粘性粘性圆柱后部：猫眼圆柱后部：猫眼1.分离现象分离现象在顺压梯度区（在顺压梯度区（BC)：流

34、体加速流体加速 在逆压梯度区（在逆压梯度区（CE)：CS段减速段减速 S点停止点停止SE段倒流。段倒流。3.3.分离的条件分离的条件 逆压梯度逆压梯度4.4.分离的实际发生分离的实际发生 微团滞止和倒流微团滞止和倒流2022-6-8流体绕流流动49分离实例分离实例从静止开始边界层发展情况从静止开始边界层发展情况扩张管扩张管（上壁有抽吸）（上壁有抽吸） 边界层分离边界层分离2022-6-8流体绕流流动5010.3.2 压差阻力DpDDLDFFFFFFf压差阻力摩擦阻力总阻力升力总阻力总作用力CD=CDf+CDp形状阻力（压差阻力）：粘性流形状阻力（压差阻力）：粘性流体绕流时，在物体表面上所作用体

35、绕流时，在物体表面上所作用的压力的合力在流动方向上的投的压力的合力在流动方向上的投影。影。对非流线型物体，是由于边界层对非流线型物体，是由于边界层的分离，在物体尾部形成旋涡，的分离，在物体尾部形成旋涡，旋涡区的压强较物体前部低，在旋涡区的压强较物体前部低，在流动方向上产生了压强差，形成流动方向上产生了压强差，形成了作用于物体上的阻力，称为压了作用于物体上的阻力，称为压差阻力。压差阻力主要取决于物差阻力。压差阻力主要取决于物体的形状。体的形状。2022-6-8流体绕流流动51一一 摩擦阻力摩擦阻力 是由于流体的粘性引起的，当流体绕流物体时，在表是由于流体的粘性引起的，当流体绕流物体时，在表面上形

36、成了边界层，边界层内速度梯度大，粘性的牵面上形成了边界层，边界层内速度梯度大，粘性的牵制作用使物体受到阻力。阻力发生在运动物体表面上。制作用使物体受到阻力。阻力发生在运动物体表面上。 摩擦阻力特点摩擦阻力特点1) 阻力系数强烈地依赖于雷诺数；阻力系数强烈地依赖于雷诺数；2) 对相同雷诺数，层流态的阻力明显低于湍流态；对相同雷诺数，层流态的阻力明显低于湍流态；3) 对湍流边界层，光滑壁面的阻力最小，粗糙度增加使阻力系对湍流边界层，光滑壁面的阻力最小，粗糙度增加使阻力系数增大；数增大；4) 摩擦阻力与壁面面积成正比。摩擦阻力与壁面面积成正比。2022-6-8流体绕流流动52二二 压差阻力压差阻力

37、与边界层的分离现象密切相关。当流体流过一个圆与边界层的分离现象密切相关。当流体流过一个圆头尖尾的回转体时，在物体前端形成减速区，在前端头尖尾的回转体时，在物体前端形成减速区，在前端顶点顶点A形成驻点，流体压强随流速变化而变化，在驻形成驻点，流体压强随流速变化而变化，在驻点处最大，离开驻点，压强逐渐减小，从点处最大，离开驻点，压强逐渐减小，从B点处开始点处开始变成负值，过最大速度点变成负值，过最大速度点C后，流速减小，压强上升，后，流速减小，压强上升，压强又变成正值。压强又变成正值。2022-6-8流体绕流流动53 压强分布如实线所示，虚线理想压强分布。压强分布如实线所示，虚线理想压强分布。 从

38、图中可以看出，前端的正压强产生一个向后的水平从图中可以看出，前端的正压强产生一个向后的水平合力，后端的正压强产生一个向前的水平合力，中段压强合力，后端的正压强产生一个向前的水平合力，中段压强为负值，产生吸力，其前半部合成一向前的水平力，后半为负值，产生吸力，其前半部合成一向前的水平力，后半部合成一向后的水平力，这两者数值相差不大，几乎相互部合成一向后的水平力，这两者数值相差不大，几乎相互抵消。因此，物体所受的水平合力取决于前端正压强造成抵消。因此，物体所受的水平合力取决于前端正压强造成的向后的较大的力与后端正压强造成的向前的较小的力，的向后的较大的力与后端正压强造成的向前的较小的力，相互抵消后

39、，还剩下向后的反物体前进的力，即压差阻力。相互抵消后，还剩下向后的反物体前进的力，即压差阻力。物体形状物体形状后部逆压梯度后部逆压梯度压强分布压强分布压强合力压强合力用实验方法确定形状阻力用实验方法确定形状阻力阻力曲线阻力曲线2022-6-8流体绕流流动54内蒙古工业大学工程流体力学电子课件内蒙古工业大学工程流体力学电子课件2022-6-8流体绕流流动55内蒙古工业大学工程流体力学电子课件内蒙古工业大学工程流体力学电子课件2022-6-8流体绕流流动56圆柱体绕流问题的重点是考察圆柱体表面及其附近的流动。一、低雷诺数下的绕流当 时，整个流场呈稳定的层流状态，且上下流场是对称的。低雷诺数下，圆柱

40、体对流场的影响区域较大，在距离圆柱体数倍柱体直径的地方，流体的速度仍然与来流速度不同，此时圆柱体受到的阻力仅为摩擦阻力。1Re2022-6-8流体绕流流动57内蒙古工业大学工程流体力学电子课件内蒙古工业大学工程流体力学电子课件二、中等雷诺数下的绕流中等量诺数下的绕流 随着雷沿数的增大，上下游对称消失，迎流面的 流动与理想流体相似，但背流面出现边界层分离、产生尾迹流。在35 Re30-40的范围内，尾迹区有较弱的对称旋涡。如图l011(a)所示。 2022-6-8流体绕流流动581. 圆柱表面压强系数分布圆柱表面压强系数分布2. 阻力系数随阻力系数随Re数的变化数的变化DCf Re10.4 绕圆

41、柱体的流动分析10.4.1 绕圆柱体的流动2022-6-8流体绕流流动595）565 103 10Re 6）63 10Re 1）1Re （图图(a)）2）1500Re（图图(b)(c)）3）55002 10Re （图（图(d)）4）552 105 10Re（图图(e)）2022-6-8流体绕流流动60 大雷诺数下的绕流 当150300Rel.9105时，在圆柱体的迎流 而上形成的边界层为层流，边界层与物面的分离点在迎流面，分离点与来流 的夹角为85左右，这种情况称为亚临界状态亚临界状态；此时物体的阻力以压差阻力 为主。在雷诺数大于1.9105的条件下，边界层小的流动逐渐向湍流过渡。 当 Re6

42、.7105 时，边界层在分离前已由层流转变为湍流，且分离点向后 移动至背流而，分离点与来流的夹角为135左右，这种状态称为超临界状 态。这时摩擦阻力有所增大，而压差阻力有所减小。但摩擦阻力与压差阻力 相比可忽略2022-6-8流体绕流流动612022-6-8流体绕流流动623. 卡门涡街卡门涡街1）定义：在圆柱绕流中，）定义：在圆柱绕流中， 涡旋从圆柱上交替脱涡旋从圆柱上交替脱 落，在下游形成有一落，在下游形成有一 定规则，交叉排列的定规则，交叉排列的 涡列。涡列。2）Re范围：范围：60-50003）Sr（斯特劳哈尔）数：斯特劳哈尔）数：19 70 198 1fd.Sr.URe2022-6-8流体绕流流动6310.4.2 不同形状物体的阻力系数不同形状物体的阻力系数1. 二维钝体二维钝体（1）光滑圆球阻力曲线）光滑圆球阻力曲线Re1时时24DCRe,3DFdU（2）粗糙圆球阻力曲线）粗糙圆球阻力曲线4. 钝体绕体阻力特点：钝体绕体阻力特点： (1) 头部形状头部形状5. 流线型体流线型体2. 三维钝体三维钝体3. 圆球：圆球：(2) 后部形状后部形状(3) 物体长度物体长度(4) 表面粗糙度表面粗糙度2022-6-8流体绕流流动64不同形状物体的阻力系数不同形状物体的阻力系数 光滑圆球阻力曲线光滑圆球阻力曲线粗糙圆球阻力曲线粗糙圆球阻力曲线