流体力学粘性流体绕物体的流动

1、工程流体力学第八章 粘性流体绕物体的流动第八章 粘性流体绕物体的流动实际流动都是有粘流动，目前对粘性流动研究方法主要有：1、基于N-S方程的紊 流模拟2、流体实验流动分类 根据工程的实际情况，流动可分为： 内流和外流。 内流内流 ：如右上图。 外流：外流： 如右下图。本章的主要内容 本章主要讨论绕流问题，即外流问题。首先将介绍粘性流体的运动微分方程，然后将给出边界层的概念及其控制方程，最后针对绕流流动现象的一些具体问题进行了讨论。 空间流动三维问题，NS方程及其求解 扰流阻力及其计算 附面层的问题 第一节第一节 不可压缩粘性流体的运动微分方程不可压缩粘性流体的运动微分方程 以流体微元为分析对象

2、，流体的运动方程可写为如下的矢量形式： 这里 ： 是流体微团的加速度，微分符号： 称为物质导数或随体导数，它表示流体微团的某性质 时间的变化率。 PFDtDVVVVVtDtDiixVttDtDV(8-1)(8-2)(8-3)一、微元体的受力分析和运动微分方程的推导 如图所示，控制体的各边长分别为dx,dy,dz，微元体的体积为： （ 84） 作用在微元体上的质量力为 ，其可用 三个分量 表示为： （85）这里： （86）如果的三个分量是 ，则： （87）dxdydzdVbFkFjFiFFzyxFdxdydzFbzyxFFF,kFjFiFFbzbybxb作用在微元体上的表面力 将微元体六个面上的

3、应力分别投影到三个坐标方向上如图作用于微元体个面上的x轴方向的应力 把作用于控制体上x方向的力叠加起来，得到作用在微元体上的表面力在x方向的分量为： dxdydzzyxdzdxdyzdydxdzydxdydzxzxyxxzxyxx作用于微元体个面上的Y、Z轴方向的应力 v同理，表面力在y方向的分量为：v表面力在z方向的分量为： dxdydzxzyxyzyydxdydzyxzyzxzz作用在微元体上的表面力 v如果用 , 和 表示单位体积的表面力，则： （ 88） zPyPxPyxzPxzyPzyxPyzxzzzxyzyyyzxyxxx作用在微元体上的表面力 将上式和式（87）代入式（81）则得

4、： （89） 这就是微分形式的运动方程。 yxzFDtDvxzyFDtDvzyxFDtDvyzxzzzzxyzyyyyzxyxxxx二、本构方程v本构方程是确立应力和应变率之间关系的方程式。斯托克斯通过将牛顿内摩擦定律推广到了粘性流体的任意流动中，建立了牛顿流体的本构方程： （810） 上式也称为广义牛顿定律 zvvpyvvpxvvpzvyvxvzvyvxvzzyyxxyzyzzyzxzxxzxyyxxy232232232三、纳维斯托克斯方程（简称NS方程）v将式（810）代入式（89）可得： （ 811） 上式称纳维斯托克斯（Naver-Stokes）方程，是粘性流体运动微分方程的又一种形式

5、。2222222222222222221()()31()()31()(3yxxxxxzxyyyyyxzyyxzzzzzD vpFD txxyzxxyzD vpFD tyxyzyxyzD vpFD tzxyzzxy )zzv对于不可压流体，其连续方程为：v对于不可压缩粘性流体，粘性体膨胀应力为零，其运动方程为： （812）0zvyvxvzyx222222222222222222zvyvxvzpFDtDvzvyvxvypFDtDvzvyvxvxpFDtDvzzzzzyyyyyxxxxx三、纳维斯托克斯方程（简称NS方程） 并考虑到拉普拉斯算子： 不可压缩粘性流体的运动方程还可写为： （813） 2

6、222222zyxzzzyyyxxxvzpFDtDvvypFDtDvvxpFDtDv222三、纳维斯托克斯方程（简称NS方程）v如果质量力只有重力作用，用 代表重力加速度，不可压缩粘性流体的运动方程的矢量形式为： (8-14) 右端第一项表示单位质量的质量力；第二项代表作用于单位质量流体的压强梯度力；第三项代表黏性变形应力。g2D-pDtVgv三、纳维斯托克斯方程（简称NS方程）v对理想流动，认为流体无粘性， ，这时运动方程简化为欧拉方程： （815） 或矢量形式 （816） 0 xxyyzzD vpFD txD vpFD tyD vpFD tzpDtDFV三、纳维斯托克斯方程（简称NS方程）

7、当流体静止不动时， ，则运动方程简化为： （817） 0V000zpFypFxpFzyx三、纳维斯托克斯方程（简称NS方程）第二节第二节 蠕动流动蠕动流动v蠕动流动：雷诺数很低的流动。v特点：流动的尺度和流动的速度均很小v如：热电厂锅炉炉膛气流中绕煤粉颗粒、 油滴等的流动；滑动轴承间隙中的流 动等等。 一、蠕动流动的微分方程一、蠕动流动的微分方程v对于定常流动，忽略惯性力和质量力，在直角坐标系下，可把纳维尔斯托克斯方程（814）组简化成 ： （818）)()()(222222222222222222zvyvxvzpzvyvxvypzvyvxvxpzzzyyyxxx一、蠕动流动的微分方程一、蠕动

8、流动的微分方程如果流动是不可压缩流体，则连续性方程为： （819） 将式（818）依次求 、 、 ，然后相加，并结合连续性方程，即得： 即蠕动流动的压力场满足拉普拉斯方程。0zvyvxvzyx22xp22yp22zp02222222pzpypxp(8- 20)二、绕球的蠕动流动 v对如图所示的无穷远来流以速度 均匀平行流沿 轴绕半径为 的静止圆球流动，得速度与压 强分布为： （821） 300330030231( , )cos (1)2231( , )sin (1)443( , )cos2rrrv rrrrrv rrrrp rprVx0r二、绕球的蠕动流动 式中 为无穷远处来流的压力。 圆球以

9、很小的速度在静止流体中作等速运动时，在流场中通过x轴的平面上的流谱如图所示。 p二、绕球的蠕动流动v在圆球的前后两驻点A和B处的压强是压强的最高点和最低点，分别为：在前驻点A（ 180 ） （822） 在后驻点B（ 0）： （823）v而切应力的最大值，发生在C（ 90）为： （824） 等于A、B点处的压强与无穷远处的压强之差的绝对值。023rVppA023rVppB02/3rVc二、绕球的蠕动流动v球面上的压强和剪切应力也可根据速度分布公式算出，为： （8-25） 对上述两式积分，可分别得到作用在球面上的压强和切应力的合力。将这两个合力在流动方向的分量相加，可得到流体作用在圆球上的阻力为：

10、 （8-26） 这就是圆球的斯托克斯阻力公式。式中d=2 为圆球的直径。 00,00,03()2cos213()()sin2rr rr rrrr rvppprrvvvvrrrrr 063DFrd 0r第三节第三节 边界层的概念边界层的概念 边界层：物体壁面附近存在大的速度梯度的薄层。 我们可以用如图所示的绕平板的流动情况说明边界层的概念。边界层的定义v粘性流体绕流物体时，由于粘性的作用，在物体的表面附近，存在一速度急剧变化的薄层边界层。 例如：来流 的流体绕流平板时，在平板表面形成边界层。V边界层的定义v在平板的前部边界层呈层流状态，随着流程的增加，边界层的厚度也在增加，层流变为不稳定状态，流

11、体的质点运动变得不规则，最终发展为紊流，这一变化发生在一段很短的长度范围，称之为转捩区，转类区的开始点称为转捩点。转类区下游边界层内的流动为紊流状态。v在转捩区和紊流区的壁面附近，由于流体的质点的随机脉动受到平板壁面的限制，因此在靠近壁面的更薄的区域内，流动仍保持为层流状态，称为层流底层和粘性底层。 V边界层的特点v边界层内速度梯度很大，旋涡强度大，有旋流动惯性力和粘性具有相同的数量级，同时考虑。v边界层外部速度梯度很小，可以作为理想流体的势流处理。v边界层厚度随 的增大而增大，随 的增大而减小。v由于边界层很薄，因而可以近似认为，边界层任一截面上各点压强相等。VeR边界层的分类按流动状态，可

12、分为层流边界层和紊流边界层。判别准则雷诺准则： 平板上的临界雷诺数 = 边界层的构成： 1.层流边界层，当 较小时，边界层内全为层流，称为层流边界层。 2.混合边界层：除前部起始部分有一小片层流区，其余大部分为紊流区，称为混合边界层。eR51056103eR边界层的厚度v两个流动区域之间并没有明显的分界线。v边界层的厚度：通常，取壁面到沿壁面外法线上速度达到势流区速度的99处的距离作为边界层的厚度，以表示，这一厚度也称边界层的名义厚度。v边界层的厚度取决于惯性和粘性作用之间的关系，即取决于雷诺数的大小。雷诺数越大，边界层就越薄；反之，随着粘性作用的增长，边界层就变厚。沿着流动方向由绕流物体的前

13、缘点开始，边界层逐渐变厚。 第四节第四节 平面层流边界层的微分方程平面层流边界层的微分方程v在这一节里，将利用边界层流动的特点如流体的粘度大小、速度与温度梯度大和边界层的厚度与物体的特征长度相比为一小量等对N-S方程进行简化从而导出层流边界层微分方程。在简化过程中，假定流动为二维不可压定常流，不考虑质量力，则流动的控制方程N-S方程为： （8-27）222222221()1()0 xxxxxyyyyyxyyxvvvvpvvvxyxxyvvvvpvvvxyyxyvvxy 第四节第四节 平面层流边界层的微分方程平面层流边界层的微分方程v将上述方程组无量纲化。为此考虑如图所示的一半无穷绕流平板，假定

14、无穷远来流 的速度 ，流动绕过平板时在平板附近形成边界层，其厚度为 ，平板前缘至某点的距离为 。取 和 为特征量，可定义如下 的无量纲量： / / / / /（ ）V/xx yy /xxvvv/yyvvvpp 2v代入方程组（827），整理后得： （8-28）式中雷诺数 /2/2/2/2/2/22/2/2/2/2/1()Re111 11()11()Re1111()()011xxxxxylyyxylyxvvvvpvvxyxxyvvvvpvvxyyxyvvxy Relv l第四节第四节 平面层流边界层的微分方程平面层流边界层的微分方程 v与 相比较是很小的 ，即 或 / 1，同时注意到， 与 、

15、与 、 与 具有同一数量级，于是 、 、 和 的量级均为1，并可以得到： 1 1 1 为了估计其他各量的数量级，由连续性方程可得： 1xvvxyxvxyp/xvx2/2xvx/yvx2/2yvx2/1/yvy/xvx第四节第四节 平面层流边界层的微分方程平面层流边界层的微分方程第四节第四节 平面层流边界层的微分方程平面层流边界层的微分方程v因此 ，于是又得到： 1 通过分析方程组（828）各项的数量级，方程组（828）中第二式中各惯性项可以忽略掉 ，同时可以略 去 、 、 。于是在方程组（828）的粘性 项中只剩第一式中的一项 。 /yv/xvy2/2xvy/yvy2/2yvy/12/2xvx

16、2/2xvy2/2yvy2/2yvxv如果仅保留数量级为1的项，而将数量级比1小的各项全部略去，再恢复到有量纲的形式，便可以得到层流边界层的微分方程组为： （8-29） 沿边界层上缘由伯努利可知： 常数 上式对 求导，得： 00122yvxvypyuxpyvvxvvyxxyxx2/2bbpxbbbdpddxdx 第四节第四节 平面层流边界层的微分方程平面层流边界层的微分方程v这样，层流边界层的微分方程又可写为： （8-30） 方程组（830） 即为在物体壁面为平面的假设下得到的边界层微分方程 。 022yvxvyvdxdVVyvvxvvyxxbbxyxx第四节第四节 平面层流边界层的微分方程平

17、面层流边界层的微分方程第五节第五节 边界层的动量积分关系式边界层的动量积分关系式边界层的动量积分方程是对边界层内流动的再简化。其推导过程有两种方法：一种是沿边界层厚度方向积分边界层的方程组，一种是在边界层内直 接应用动量守恒原理。下面的推导采用第二种方法。边界层动量积分方程的推导边界层动量积分方程的推导如图所示为不可压缩流体的定常二维边界层流动 ，设物体表面型线的曲率很小。 取一个单位厚度的微小控制体，它的投影面ABDC 。用动量定理来建立该控制体内的流体在单位时间内沿x方向的动量变化和外力之间的关系。 边界层动量积分方程的推导边界层动量积分方程的推导设壁面上的摩擦应力为 根据边界层的控制方程

18、组，边界层内的压强仅近似地依赖于 而与 无关，设AB面上的压强为 ，DC上的压强为 控制面AC为边界层的外边界 其外部为理想流体的势流 ，只有与之垂直的压力 ，设AC上的压强为A，C两点压强的平均值 。作用在控制体上的表面力沿方向的合力为： wppd xx12ppdxx1sin2xwppFppdx dspdxddxxxyp边界层动量积分方程的推导v式中为边界层外边界AC与方向的夹角，由几何关系可知： ，上式经整理并略去高阶小量，得：v单位时间内沿方向经过AB流入控制体的质量和动量分别为：v经过CD面流出的质量和动量分别为：v定常流动条件下，可知从控制面AC流入控制体中的流量为：v由此引起流入的

19、动量为： sindsadxwpFd xd xx 0A Bxmd y 20ABxkdy0()xCDxvmvdx dyx2200 CDxxkv dydxv dyx0()CDABxdmmv dy dxdx0ACbxdkv dxv dydx边界层动量积分方程的推导v式中V为边界层外边界上的速度。这样，可得单位时间内该控制体内沿x方向的动量 变化为v 根据动量定理， ，则可得边界层的动量积分方程为：v （8-51） 上式也称为卡门动量积分关系式。该式是针对边界层流动在二维定常流动条件下导出的，并没有涉及边界层的流态，所以其对层流和紊流边界层都能适用。 xk200 xCDABACxbxddkkkkv dy

20、vv dy dxdxdxxxFk200wxbxdpddv dyVv dydxdxdx积分方程的求解v实际上可以把 、 和 看作已知数，而未知数只有 、 和 三个。 v再补充两个关系式：v 一、沿边界层厚度的速度分布 = (y) 二、切向应力与边界层厚度的关系式 一般在应用边界层的动量积分关系式（851）来求解边界层问题时，边界层内的速度分布是按照已有的经验来假定的。假定的 愈接近实际，则所得到的结果愈正确。所以选择边界层内的速度分布函数 是求解边界层问题的重要关键。 ( ) ( )xvv y( )v ybVdpdxxwxx第六节边界层的位移厚度和动量损失厚度 v边界层的厚度 ，表示粘性影响的范

21、围。 位移厚度 动量损失厚度v根据伯努力方程可知：v又由于：v带入（8-51）得 或 （8-52）bbdvdpvdxdx 00bbbbdvdvdpvdyv dydxdxdx 000bbxbxxdvddvv dyv v dyv dydxdxdx20000bbxxbxbwdvdvddv dyv v dyv dyv dydxdxdxdx00()()bbxxbxwdvdvv dyv vv dydxdx12边界层厚度计算式的推导因此在边界层内由于粘性影响使体积流量的减小量 ，即上式中第一项积分。 位移厚度或排挤厚度 可表示成： （8-53）同理动量损失厚度 可表示为： （8-54） 将 和 代入式（85

22、1）, 得 （8-55）0()xVv dy11001()(1)xbxbbvvv dydyvv2221()bbbwdvdvvdxdx12201()(1)xxxbxbbbvvv vv dydyvvv 2边界层厚度计算式的推导v式（8-55）是另一种形式的平面不可压缩粘性流体边界层动量积分关系式 。 、 和 都是未知数，它们决定于边界层内速度的分布规律。 v将式（855）化为无因次形式，统除以 ，得 （856） 或v式中H 。计算曲面边界层时，用上式较为方便。 12w2V022121(2)bbbdvddxvdxv022122(2)bbbdvdHdxvdxv21/第七节 平板边界层流动的近似计算 v平

23、板层流边界层的近似计算 对于式（851），如果边界层外部的压强梯度为零，方程变为： (8-57） 假定平板非常薄，对流动没有影响。边界层外层流动： 则上式可变为： (8-58） 两个补充关系式:一、冯卡门假定，二、牛顿内摩擦定律。v平板紊流边界层的近似计算 采用将边界层内的速度分布与圆管内充分发展紊流的速度分布规 律进行类比的方法。0()wxbxdv vv dydx0()wxbxdvvvdydx平板层流边界层的近似计算v选择一三次项式速度分布： (8-59) 根据下列边界条件来确定待定系数 和 . (1)在平板壁面上的速度为零，即在 处 (2)在边界层外边界上的速度等于来流速度 ,即在 处 ，

24、 (3)在边界层外边界上，摩擦切应力 为零，即在 处 ， (4)由于在平板壁面上的速度为零，即 ，由方程组（850）的第一式得 y332210yayayaavx210aaa、3a0y0 xvxvvyvxy0)(yxyv0yxvv01)(dxdpyvyx平板层流边界层的近似计算v速度分布的四个系数可确定为： v于是，层流边界层中速度的分布规律为 (8-60) 第二个补充关系式：利用牛顿内摩擦定律和式（860）得出 (8-61)式中为动力粘性系数。将速度分布方程（860）带入方程（861）并积分得：分离变量，并积分得： (8-62)00a 132va20a 332va 331( ) ( )( )

25、22xyyvyv023xwydvvdy（）23292280VdVdx4.644.64/Rexxv x平板层流边界层的近似计算式中为 运动粘性系数，为基于长度的雷诺数 。合并方程（862）和(861)得到： (8-63)如果表面摩擦系数 为： (8-64) 那么 ，为： (8-65) 根据动量损失厚度的定义式（854），并考虑式（862），可得动量损失厚度为： (8-66) 同理，位移厚度为: (8-67)上述计算结果是依赖于所假设的速度分布规律的，不同阶次的速度分布，可以得出不同的结果。表8.1 给出几种不同的情况。Rex20.323RewxvfC0212fCvfC0.646RefxC 20.

26、646Rexx11.740Rexx表8.1不同阶次的速度分布所得结果比较 548 1826 0730 0365464 1740 0646 0323584 1751 0685 0343 VvxxxRexxRe1xxRe220ReVx22yy4222yyy32123yy二、二、平板紊流边界层的近似计算v如前所述由于流动的混参以及速度和压力的波动，紊流边界层的速度分布都采用一些模型假定。普朗特建议，当边界层雷诺数 时，边界层内的速度分布可采用 次方规律，即： （8-68） 该式不能直接应用于边界层的内边界。通常认为粘性底层内的速度分布为线形分布。 雷诺数取 时的摩擦阻力系数为： 当时 普朗特和施利希

27、廷 （ H. Schlichting）采用对数速度分布，得到如下的半经验公式： 7Re10 x7117( )xyvV575 10Re10l150.074ReflC9710Re10l2.580.455(lg Re )flC层流与紊流边界层的近似计算公式汇总 q平板的层流边界层和紊流边界层的重大差别有：1)紊流边界层内沿平板壁面发向截面上的速度比层流边界层的速度增加得快 2)沿平板壁面紊流边界层的厚度比层流边界层的厚度增加得快 3)在其它条件相同的情况下，平板壁面上的切向应力 沿着壁面的减小在紊流边界层中要比层流边界层减小得慢。4)在同一 下，紊流边界层得摩擦阻力系数比层流边界层的大得多 实际情况

28、下，边界层是层流和紊流同时存在的混合边界层 xRe071yVvx5151Re37.037.0 xxxVx51Re0462. 0125. 0 xx51Re036. 01 . 0 xx512512Re0289.00289.0 xVxVV512Re036. 0lVbl4322yyyVvx21Re84. 584. 5xxVx21Re752. 13 . 0 xx21Re686. 01175. 0 xx2122Re343. 0343. 0 xVxVV212Re686. 0lVbl21Re372. 1l51Re074. 0l边界 层 内 的 流 态层 流紊 流边界层的基本特性速度分布规律边界层厚度位移厚度动量损失厚度切向应力总摩擦力摩擦阻力系数层流与紊流边界层的近似计算公式汇总 210DFfC第八节 边界层的分离与卡门涡街 一、边界层的分离一、边界层的分离v以如图所示的圆柱绕流为例 在势流流动中流体质点从D到E是加速的，为顺压强梯度;从E到F则是减速的, 为逆压强梯度流体质点由D到E过程，由于流体压能向动能的转变，不发生边界层分离E到F 段动能只存在损耗，速度减小很快，在S点处出现粘滞 ，由于压力的升高产生回流导致边界层分离，并形成尾涡。如图为边界层分离示意图。v从以上的分析中可得如下结论：粘性流体在压力降低区内流动（加速流动），决不会出现边界层的分离，只有在压力升高区内流动（减速流动），才有