第六章粘性流体的一维定常流动

1、2022-3-7工程流体力学第六章 粘性流体的一维定常流动第一节 黏性流体总流的伯努利方程第二节 黏性流体的两种流动型态 第三节 流动损失分类第四节 圆管中流体的层流流动第五节 圆管中流体的紊流流动第六节 沿程阻力系数的实验研究第七节 非圆形截面管道沿程损失的计算第八节 局部损失的计算第九节 管 道 水 力 计 算第十节 水击现象2022-3-7工程流体力学 在第三章中，通过对理想流体运动的基本规律的讨论，得到了流场中任一空间点上、任一时刻流体微团的压强和速度等流动参数之间的关系式，但在推导流体微团沿流线运动的伯努利方程中，仅局限于微元流束的范围内。而在工程实际问题中要研究实际流体在整个流场中

2、的运动，其中大量的是在管道和渠道中的流动问题。所以除了必须把所讨论的范围从微元流束扩展到整个流场（如管道）外，还需考虑黏性对流体运动的影响，实际流体都具有黏性，在流动过程中要产生摩擦阻力，为了克服流动阻力以维持流动，流体中将有一部分机械能不可逆地损失掉。由此可见，讨论黏性流体流动的重点就是讨论由于黏性在流动中所造成的阻力问题，即讨论阻力的性质、产生阻力的原因和计算阻力的方法。 2022-3-7工程流体力学第一节 黏性流体总流的伯努利方程一、黏性流体微元流束的伯努利方程一、黏性流体微元流束的伯努利方程在第三章中已经得到了理想不可压缩流体作定常流动时，质量力仅为重力情况下的微 元流束的伯努利方程，

3、该式说明流体微团沿流线运动时总机械能不变。但是对于黏性流体， 在流动时为了克服由于黏性的存在所产生的阻力将损失掉部分机械能，因而流体微团在流 动过程中，其总机械能沿流动方向不断地减少。如果黏性流体从截面1流向截面2，则截面2处的总机械能必定小于截面1处的总机械能。若以 表示单 位重量流体自截面1到2的流动中所损失的机械能（又称为水头损失），则黏性流体微元流束的伯努利方程为 (6-1)式（6-1）的几何解释如图6-1所示，实际总水头线沿微元流束下降，而静水头线则随流束的形状上升或下降。 WhwhgVgpzgVgpz22222221112022-3-7工程流体力学图6-1 伯努利方程的几何解释20

4、22-3-7工程流体力学二、黏性流体总流的伯努利方程二、黏性流体总流的伯努利方程 流体的实际流动都是由无数微元流束所组成的有效截面为有限值的总流流动，例如流体在管道中和渠道中的流动等。 微元流束的有效截面是微量，因而在同一截面上流体质点的位置高度 、压强 和流速 都可认为是相同的。而总流的同一有效截面上，流体质点的位置高度 、压强和流速 是不同的。总流是由无数微元流束所组成的。因此，由黏性流体微元流束的伯努利方程来推导总流的伯努利方程，对总流有效截面进行积分时，将遇到一定的困难，这就需要对实际流动作某些必要的限制。为了便于积分，首先考虑在什么条件下总流有效截面上各点的 常数？这只有在有效截面附

5、近处有缓变流动时才能符合这个要求。zzppVVgpz2022-3-7工程流体力学 由于流线几乎是平行直线，则各有效截面上相应点的流速几乎不变，成为均匀流，由于速度的变化很小即可将惯性力忽略不计，又由于流线的曲率半径很大，故向心力加速度很小，以致可将离心力忽略。于是缓变流中的流体微团只受重力和压强的作用，故缓变流的有效截面上各点的压强分布与静压强分布规律一样，即在同一有效截面上各点的 常数。当然在不同的有效截面上有不同的常数值。 掌握了缓变流动的特性之后，就可以将黏性流体微元流束的伯努利方程应用于总流，从而推导出适用于两个缓变流有效截面的黏性流体总流的伯努利方程。gpz2022-3-7工程流体力

6、学 以总流中每一微元流束的任意两个截面可以写出 则通过该微元流束的总能量在截面1与截面2之间的关系式为积分上式，则得总流在有效截面1和有效截面2之间的总能量关系式 （6-2）whgVgpzgVgpz2222222111VwVVqghqggVgpzqggVgpzdd2d222222111VVVqVwqVqVqghqggVgpzqggVgpzdd2d2222221112022-3-7工程流体力学若有效截面1和有效截面2处的流动都是缓变流动，则 和 ， 和 是两个不同的常数，于是式（6-2）可写成 （6-3）对于不可压缩流体，以 通除式（6-3）各项得 （6-4）用有效截面上的平均流速 代替真实流速

7、 ，则可将式（6-4）中总流的平均单位重量 流体的动能项改写为 （6-5）式中 总流的动能修正系数 （6-6）111Cgpz222Cgpz1C2CVVVVVqVwqVqVqVqVqghqggVqggpzqggVqggpzdd2dd2d22222111VqVVgqqgdVVVqVwVqVVqVVqhqqgVqgpzqgVqgpzd1d21d2122222111VVVqAAVVgVAgVVVAAVgVVAqgVq2d21d21d21223222AAVVAd132022-3-7工程流体力学以 表示总流有效截面1和有效截面2之间的平均单位重量流体的能量损失，即 （6-7）将式（6-5）和式（6-7）代

8、人式（6-4）中得： （6-8） 这就是黏性流体总流的伯努利方程。适用范围是：重力作用下不可压缩黏性流体定常流动的任意两个缓变流的有效截面，至于两个有效截面之间是否是缓变流则无关系。由式(6-8)可以看出，如同黏性流体沿微元流束的流动情况一样，为了克服流动阻力，总流的总机械能即实际总水头线也是沿流线方向逐渐减少的，如图6-2所示。WhVqVVqhqhd1WWwhgVgpzgVgpz2222222211112022-3-7工程流体力学图6-2 总流总水头线2022-3-7工程流体力学 动能修正系数 是由于截面上速度分布不均匀而引起的，它可按式（6-6）根据有效截面上的速度分布规律而求得。 是个大

9、于1的数，有效截面上的流速越均匀， 值越趋近于1。在实际工业管道中，通常都近似地取 。以后如不加特别说明，都假定 ，并以 代表平均流速。而对于圆管层流流动 。0 . 112V2022-3-7工程流体力学【例例6-1】 有一文丘里管如图6-3所示，若水银差压计的指示为360mmHg，并设从截面A流到截面B的水头损失为0.2mH2O， =300mm， =150mm，试求此时通过文丘里管的流量是多少?图6-3 文丘里管AdBd2022-3-7工程流体力学【解解】 以截面A为基准面列出截面A和B的伯努利方程由此得 （a）由连续性方程所以 （b）w2BB2AA276. 020hgVgpgVgp2 . 0

10、76. 0222A2BBAgVgVgpgpBBAAAVAV2ABBABBAddVAAVV2022-3-7工程流体力学水银差压计11为等压面，则有由上式可得 （c）将式（b）和式（c）代入（a）中解得 （m/s） （m3/s）ggzpgzpHgBA36. 076. 036. 0）（）（）（OmmH3 . 5980613340036. 040. 0g36. 036. 076. 02HgBAggpgp96. 0123 . 542ABBddgV53. 93001501)96. 03 . 5(806. 921)96. 03 . 5(244ABBddgV168. 015. 0453. 9422BBdVqV

11、2022-3-7工程流体力学【例例6-2】 有一离心水泵装置如图6-4所示。已知该泵的输水量 m3/h，吸水管内径 150mm，吸水管路的总水头损失 mH2O，水泵入口22处，真空表读数为450mmHg，若吸水池的面积足够大，试求此时泵的吸水高度 为多少? 60Vqd5 . 0whgh图6-4 离心泵装置示意图2022-3-7工程流体力学【解解】 选取吸水池液面l1和泵进口截面22这两个缓变流截面列伯努利方程，并以11为基准面，则得因为吸水池面积足够大，故 。且 （m/s） 为泵吸水口截面22处的绝对压强，其值为将和值代入上式可得 （mH2O）w22221a220hgVgphgVgpg01V9

12、4. 015. 014. 336006044222dqVV45.01330002app2pw22245. 0133000hgVghg5 . 0806. 9294. 0980645. 0133000256.52022-3-7工程流体力学第二节第二节 黏性流体的两种流动型态黏性流体的两种流动型态从上节式（6-8）的黏性流体总流的伯努利方程可以看出，要想应用此关系式计算有关工程实际问题，必须计算能量损失 项，由于流体流动的能量损失与流动状态有很大关系，因此，我们首先讨论黏性流体流型。wh黏性流体的流动存在着两种不同的流型，即层流和紊流，这两种流动型态由英国物理学家雷诺（Reynolds）在1883年

13、通过他的实验（即著名的雷诺实验）大量观察了各种不同直径玻璃管中的水流，总结说明了这两种流动状态。2022-3-7工程流体力学 一、雷诺实验一、雷诺实验 雷诺实验装置如图6-5所示。实验的步骤如下： (1) 首先将水箱A注满水，并利用溢水管H保持水箱中的水位恒定，然后微微打开玻璃管末端的调节阀C，水流以很小速度沿玻璃管流出。再打开颜色水瓶D上的小阀K，使颜色水沿细管E流入玻璃管B中。当玻璃管中水流速度保持很小时，看到管中颜色水呈明显的直线形状，不与周围的水流相混。这说明在低速流动中，水流质点完全沿着管轴方向直线运动，这种流动状态称为层流，如图6-6(a)所示。图6-5 雷诺实验图6-6 层流、紊

14、流及过渡状态2022-3-7工程流体力学(2) 调节阀C逐渐开大，水流速度增大到某一数值时颜色水的直线流将开始振荡，发生弯曲，如图6-6(b)所示。(3) 再开大调节阀C，当水流速度增大到一定程度时，弯曲颜色水流破裂成一种非常紊乱的状态，颜色水从细管E流出，经很短一段距离后便与周围的水流相混，扩散至整个玻璃管内，如图6-6(c)所示。这说明水流质点在沿着管轴方向流动过程中，同时还互相掺混，作复杂的无规则的运动，这种流动状态称为紊流（或湍流）。如果将调节阀C逐渐关小，水流速度逐渐减小，则开始时玻璃管内仍为紊流，当水流速度减小到另一数值时，流体又会变成层流，颜色水又呈一明显的直线。但是，由紊流转变

15、为层流时的流速要比由层流转变为紊流时的流速小一些。我们把流动状态转化时的流速称为临界流速，由层流转变为紊流时的流速称为上临界流速，以cVcVccVV以表示。则表示。由紊流转变为层流时的流速称为下临界速，2022-3-7工程流体力学雷诺实验表明：当流速大于上临界流速时为紊流；当流速小于下临界流速时为层流；当流速介于上、下临界流速之间时，可能是层流也可能是紊流，这与实验的起始状态、有无扰动等因素有关，不过实践证明，是紊流的可能性更多些。在相同的玻璃管径下用不同的液体进行实验，所测得的临界流速也不同，黏性大的液体临界流速也大；若用相同的液体在不同玻璃管径下进行试验，所测得的临界流速也不同，管径大的临

16、界流速反而小。2022-3-7工程流体力学二、雷诺数二、雷诺数综上可知，流体的流动状态是层流还是紊流，与流速、管径和流体的黏性等物理性质有关。雷诺根据大量的实验数据证明，流体的临界流速cVddVc他引出一个比例系数cRedRedReVccc或dVRecc （6-9） 这个比例系数cRe与流体的动力黏度成正比，与管内径和流体的密度成反比，即，上式可写成等式 称为临界雷诺数，是一个无量纲数。2022-3-7工程流体力学 经过雷诺实验和他以后的许多学者如席勒（Ludwig Schiller）的精密实验结果指明，对于非常光滑、均匀一致的直圆管，下临界雷诺数 等于2320。但对于一般程度的粗糙壁管 值稍

17、低，约为2000，所以在工业管道中通常取下临界雷诺数 。上临界雷诺数 不易测得其精确数值，一般取为13800。于是得cRecRe2000cReceR 2000dVRecc13800dVeRcc无数实验证明，不管流速多少、管内径多大、也不管流体的运动黏度如何，只要雷诺数相等，它们的流动状态就相似。所以雷诺数是判别流体流动状态的准则数，即：2022-3-7工程流体力学当流体流动的雷诺数 时，流动状态为层流；当时 ，则为紊流；当 时，流动状态可能是层流，也可能是紊流，处于极不稳定的状态，任意微小扰动都能破坏稳定，变为紊流。 显然，上临界雷诺数在工程上一般没有实用意义，故通常都采用下临界雷诺数 作为判

18、别流动状态是层流或紊流的准则数。即：cceRReRecReRe ceRRecReVdRe VdRe 20002000是层流是紊流2022-3-7工程流体力学工程中实际流体（如水、空气、蒸汽等）的流动，几乎都是紊流，只有黏性较大的液体（如石油、润滑油、重油等）在低速流动中，才会出现层流。流体在任意形状截面的管道中流动时，雷诺数的形式是eVdRe (6-10) 式中ed雷诺数之所以能作判别层流和紊流的标准，可根据雷诺数的物理意义来解释。黏性流体流动时受到惯性力和黏性力的作用，这两个力用量纲可分别表示为22lVdtdVmVlAdydV黏性力惯性力VllVVl22Re为当量直径。惯性力黏性力2022-

19、3-7工程流体力学 由此可知雷诺数是惯性力与黏性力的比值。雷诺数的大小表示了流体在流动过程中惯性力和黏性力哪个起主导作用。雷诺数小，表示黏性力起主导作用，流体质点受黏性的约束，处于层流状态；雷诺数大表示惯性力起主导作用，黏性不足以约束流体质点的紊乱运动，流动便处于紊流状态。2022-3-7工程流体力学三、能量损失与平均流速的关系三、能量损失与平均流速的关系 如果将两根测压管接在雷诺实验装置中玻璃管B的前后两端，如图6-7所示，可测出有效截面1-1和2-2间的能量损失，并找出管中平均流速与能量损失之间的关系。列截面1-1和2-2的伯努利方程f222222111122hgVgpzgVgpz由于玻璃

20、管是等截面管，所以 ，21VV 2121zz gpphf21可见，测压管中的水柱高差即为有效截面1-1和2-2间的压头损失。并令，另外玻璃管是水平放置的，即，于是上式可写成2022-3-7工程流体力学将测得的平均流速和相应的压头损失，在对数坐标上表示出，如图4-8所 示。先做层流到紊流的试验，当流速逐渐增加时， 与 成正比增大，如图中的OAB直线。当流速增加到一定程度时层流变为紊流， 突然从B点上升到C点。以后再增大流速时， 要比 增加得快，如图中的CD线，其斜率比OAB线的斜率大，此后若将流速逐渐减小，则 与 的关系曲线沿DCAO线下降。A点和B点各为相应的下临界流速 和上临界流速 ，ABC

21、为过渡区。fhVfhfhVfhVcVcV2022-3-7工程流体力学图6-7 水平等直管道中水头损失 图6-8 层流和紊流的与的关系曲线 2022-3-7工程流体力学由实验所得的图6-8可知，当 时，即层流时， 与 的一次方成正比；当 时，即紊流时， 与 成正比。 值与管壁粗糙度有关：对于管壁非常光滑的管道 ；对于管壁粗糙的管道 .所以紊流中的压头损失比层流中的要大。cVV fhVcVVfhmVm75. 1m2m从上述讨论可以得出，流型不同，其能量损失与速度之间的关系差别很大，因此，在计算管道内的能量损失时，必须首先判别其流态（层流，紊流），然后根据所确定的流态选择不同的计算方法。2022-3

22、-7工程流体力学【例例6-36-3】 管道直径 100mm，输送水的流量 m3/s，水的运动黏度 m2/s，求水在管中的流动状态？若输送 m2/s的石油，保持前一种情况下的流速不变，流动又是什么状态？d01. 0Vq610141014. 1【解解】 （1）雷诺数 VdRe27. 11 . 014. 301. 04422dqVV20001027. 11011 . 027. 1Re56（m/s） 故水在管道中是紊流状态。2022-3-7工程流体力学 （2） 200011141014. 11 . 027. 1Re4Vd故油在管中是层流状态。2022-3-7工程流体力学第三节 流动损失分类实际流体在管

23、内流动时，由于黏性的存在，总要产生能量损失。产生能量损失的原因和影响因素很复杂，通常可包括黏性阻力造成的黏性损失fhjh 一、沿程阻力与沿程损失一、沿程阻力与沿程损失黏性流体在管道中流动时，流体与管壁面以及流体之间存在摩擦力，所以沿着流动路程，流体流动时总是受到摩擦力的阻滞，这种沿流程的摩擦阻力，称为沿程阻力。流体流动克服沿程阻力而损失的能量，就称为沿程损失。沿程损失是发生在缓变流整个流程中的能量损失，它的大小与流过的管道长度成正比。造成沿程损失的原因是流体的黏性，因而这种损失的大小与流体的流动状态（层流或紊流）有密切关系。两部分。和局部阻力造成的局部损失2022-3-7工程流体力学单位重量流

24、体的沿程损失称为沿程水头损失，以 表示，单位体积流体的沿程损失，又称为沿程压强损失，以 表示 。fhfpffghp在管道流动中的沿程损失可用下式求得gVdlh22f22fVdlp (6-11) (6-11a)式中 ldV沿程阻力系数，它与雷诺数和管壁粗糙度有关，是一个无量纲的系数，将在本章第六节进行讨论；式（6-11）称为达西-威斯巴赫（Darcy-Weisbach）公式。管道长度，m； 管道内径，m；管道中有效截面上的平均流速，m/s。2022-3-7工程流体力学二、局部阻力与局部损失二、局部阻力与局部损失 在管道系统中通常装有阀门、弯管、变截面管等局部装置。流体流经这些局部装置时流速将重新

25、分布，流体质点与质点及与局部装置之间发生碰撞、产生漩涡，使流体的流动受到阻碍，由于这种阻碍是发生在局部的急变流动区段，所以称为局部阻力。流体为克服局部阻力所损失的能量，称为局部损失。单位重量流体的局部损失称为局部水头损失，以 表示，单位体积流体的局部损失，又称为局部压强损失，以 表示 。jhjpjjghp在管道流动中局部损失可用下式求得gVhj222 2Vpf (6-12)(6-12a) 式中 局部阻力系数。 局部阻力系数 是一个无量纲的系数，根据不同的局部装置由实验确定。在本章第八节进行讨论。2022-3-7工程流体力学三、总阻力与总能量损失三、总阻力与总能量损失在工程实际中，绝大多数管道系

26、统是由许多等直管段和一些管道附件连接在一起所组成的，所以在一个管道系统中，既有沿程损失又有局部损失。我们把沿程阻力和局部阻力二者之和称为总阻力，沿程损失和局部损失二者之和称为总能量损失。总能量损失应等于各段沿程损失和局部损失的总和，即jfwhhh jfwwppghp (6-13) (6-13a) 上述公式称为能量损失的叠加原理。 2022-3-7工程流体力学第四节第四节 圆管中流体的层流流动圆管中流体的层流流动黏性流体在圆形管道中作层流流动时，由于黏性的作用，在管壁上流体质点的流速等于零，随着流层离开管壁接近管轴时，流速逐渐增加，至圆管的中心流速达到最大值。本节讨论流体在等直径圆管中作定常层流

27、流动时，在其有效截面上切应力和流速的分布规律。一、数学模型一、数学模型图6-9 等直径圆管中的定常层流流动流体在等直径圆管中作定常层流流动时，取半径为 ，长度为 的流段1-2为分析对象，如图6-9所示。作用在流段12上的力有：截面1-1和2-2上的总压力 和 ，在这里是假设截面1-1和2-2上的压强分布是均匀的；流段1-2的重力 ；作用在流段侧面上的总摩擦力 ，方向与流动方向相反。rlApP11ApP22gAlGrlT22022-3-7工程流体力学图6-9 等直径圆管中的定常层流流动由于流体在等直径圆管中作定常流动时加速度为零，故不产生惯性力。根据平衡条件，写出作用在所取流段上各力在流动轴线上

28、的平衡方程：0sin221gAlrlApAp式中： 21sinzzl2rA2022-3-7工程流体力学以 除以上式各项，整理得gAlGlgrgpzgpz22211(6-14) 对截面1-1和2-2列出伯努利方程得f222222111122hgVgpzgVgpz在等直径圆管中 ， ，故2121VV ，gpzgpzhf2211(6-15) 2022-3-7工程流体力学将式（6-15）代入式（6-14）中得lgrhf2(6-16) 在层流中切应力 可用牛顿内摩擦定律来表示，即rudd(6-17) 由于流速 随半径 的增加而减小，即 是负值，为了使 为正值，式（6-17）等号在右端取负号。urrudd

29、2022-3-7工程流体力学二、速度分布二、速度分布为了求出速度分布，现将式（6-17）代入式（6-16）中整理得rrlprrlhguffd 2d2d积分上式得Crlpuf2 4根据边界条件确定积分常数 ，在管壁上 ， ，则C0rr 0u20 4rlpCf代入上式得)( 4220rrlpuf (6-18)2022-3-7工程流体力学式（6-18）表明在有效截面上各点的流速 与点所在的半径 成二次抛物线关系，如图6-10所示。在 的管轴上，流速达到最大值：ur0r20max 4rlpuf (6-19) 图6-10 圆管中层流的速度分布 三、流量及平均流速三、流量及平均流速 现求圆管中层流的流量：

30、取半径 处厚度为d 的一个微小环形面积，每秒通过这环形面积的流量为rrrdrudqV22022-3-7工程流体力学由通过圆管有效截面上的流量为ArrfVVrrrrlprruqq0000220d2)( 4d2d400220 8d)( 20rlprrrrlpfrf (6-20) 这就是层流管流的哈根-普索勒（Hagen-Poiseuille）流量定律。该定律说明：圆管中流体作层流流动时，流量与单位长度的压强降和管半径的四次方成正比。2022-3-7工程流体力学圆管有效截面上的平均流速202040 8 8rlprlrpAqVffV (6-21) 比较式（6-19）和式（6-21）可得max21uV

31、(6-22) 即圆管中层流流动时，平均流速为最大流速的一半。工程中应用这一特性，可直接从管轴心测得最大流速从而得到管中的流量 ，这种测量层流的流量的方法是非常简便的。AuqVmax212022-3-7工程流体力学四、切应力分布四、切应力分布由牛顿内摩擦定律可得到切应力在有效截面上的分布规律。lrprrlpdrddrduff2)( 4220(6-23) 在管壁处 ， ，故式（6-23）成为0rr 0lrpf200 (6-24) 由式（6-23）和式（6-24）得00rr(6-25) 式（6-25）表明，在圆管的有效截面上，切应力 与管半径 的一次方成比例，为直线关系，在管轴心处 时 ，如图6-1

32、1所示。r0r0图6-11 圆管有效截面上的切应力2022-3-7工程流体力学五、沿程损失五、沿程损失 fh流体在等直径圆管中作层流流动时，流体与管壁及流体层与层之间的摩擦，将引起能量损失，这种损失为沿程损失。由式（6-21）可得沿程损失20ff8ggrlVph由此可见，层流时沿程损失与平均流速的一次方成正比。 由于 ，代入上式得 gVdlRegVdlVdgrlVh264223282220f2022-3-7工程流体力学令Re64 为沿程阻力系数，在层流中仅与雷诺数有关。于是得 gVdlh22f该式与式（6-11）的形式相同。 六、动能修正系数六、动能修正系数 已知黏性流体在圆管中作层流流动时的

33、速度分布规律，便可求出黏性流体总流伯努利方程中的动能修正系数 ，将式（6-18）和式（6-21）代入到式（6-6）得：22121130202030rdrrrrdAVuAAr(6-27) （6-26）2022-3-7工程流体力学【例例6-46-4】 圆管直径 mm，管长 m，输送运动黏度 cm2/s的石油，流量 m3/h，求沿程损失。200d1000l6 . 1144Vq【解解】 判别流动状态20005 .1587106 . 12 . 027. 1Re4Vd为层流 式中 27. 12 . 014. 336001444422dqVV（m/s） 由式（6-6） 57.16806. 9227. 12

34、. 010005 .1587642642222fgVdlRegVdlh（m 油柱） 2022-3-7工程流体力学【例例6-56-5】 输送润滑油的管子直径 8mm，管长 15m，如图6-12所示。油的运动黏度 m2/s，流量 12cm3/s，求油箱的水头 （不计局部损失）。 dl61015Vqh图6-12 润滑油管路 239. 0008. 014. 3101244242dqVV（m/s） 雷诺数 20005 .1271015008. 0239. 06VdRe为层流列截面1-1和2-2的伯努利方程f222211202hgVgpgVgphaa2022-3-7工程流体力学认为油箱面积足够大，取01V

35、gVdlRegVh264222222f806. 92239. 0008. 0155 .12764806. 92239. 022275. 2(m) ，则2022-3-7工程流体力学第五节第五节 圆管中流体的紊流流动圆管中流体的紊流流动 从本章第二节中的雷诺实验可知，当ceRRe 一、紊流脉动现象与时均速度一、紊流脉动现象与时均速度 流体质点在运动过程中，不断地互相掺混，引起质点间的碰撞和摩擦，产生了无数旋涡，形成了紊流的脉动性，这些旋涡是造成速度等参数脉动的原因。紊流是一种不规则的流动状态，其流动参数随时间和空间作随机变化，因而本质上是三维非定常流动，且流动空间分布着无数大小和形状各不相同的旋涡

36、。因此，可以简单地说，紊流是随机的三维非定常有旋流动。流动参数的变化称为脉动现象。 时，管内流动便会出现杂乱无章的紊流，流体运动的参数，如速度、压强等均随时间不停地变化。在紊统流动时，其有效截面上的切应力、流速分布等与层流时有很大的不同。2022-3-7工程流体力学在流场中的某一空间点如用高精度的热线热膜风速仪来测量流体质点的速度，则可发现速度是随时间而脉动的，如图6-13所示。从图中可见紊流中某一点的瞬时速度随时间的变化极其紊乱，似乎无规律可循。但是在一段足够长时间 内，即可发现这个变化始终围绕着某一平均值，在其上下脉动，这就反映了流体质点掺混过程中脉动现象的实质，揭示了紊流的内在规律性。

37、1t图6-13 脉动速度时间1t101d1ttutu（6-28） 内，速度的平均值称为时均速度，定义为2022-3-7工程流体力学于是流场的紊流中某一瞬间，某一点瞬时速度可用下式表示。uuu(6-29) 其中， 称为脉动速度，由于 流体质点在紊流状态下作不定向的杂乱无章的流动，脉动速度 有正有负。但是在一段时间内，脉动速度的平均值为零，即 。uuu0u紊流中的压强和密度也有脉动现象，同理 和 也同样可写成 pppp（6-30）在实际工程和紊流试验中，广泛应用的普通动压管只能测量它的时均值，所以在研究和计算紊流流动问题时，所指的流动参数都是时均参数，如时均速度 ，时均压强 等。为书写方便起见，常

38、将时均值符号上的“一”省略。我们把时均参数不随时间而变化的流动，称为准定常紊流。up2022-3-7工程流体力学 二、紊流中的切向应力二、紊流中的切向应力在黏性流体层流流动时，切向应力表现为由内摩擦力引起的摩擦切向应力。在黏性流体紊流流动中，与层流一样，由于流体的黏性，各相邻流层之间时均速度不同，从而产生摩擦切向应力 。t1.1.摩擦切向应力摩擦切向应力另外，由于流体有横向脉动速度，流体质点互相掺混，发生碰撞，引起动量交换，因而产生附加切应力 ，向应力是由摩擦切向应力和附加切应力两部分组成。因此紊流中的切摩擦切向应力可由牛顿内摩擦定律式（1-10）求得yudd2022-3-7工程流体力学2 2

39、附加切向应力附加切向应力 附加切向应力可由普朗特混合长度理论推导出来。 设管内紊流时均速度 的分布如图6-14所示，在流层1上某一流体质点有轴向脉动速度 和横向脉动速度 。横向脉动速度 使流体质点从流层1运动一个微小距离 到另一流层2。普朗特假定 相当于气体分子的平均自由行程。流层1上的流体的时均速度为 ，则流层2上的时均速度为 。uullulyuudd图6-14 紊流时均速度分布在 时间内，由流层1经微小面积d 流向流层2的流体质量为 t dAtAmddd2022-3-7工程流体力学质量 的流体到流层2后与该层上的流体互相碰撞，发生动量交换。在 时间内动量变化为mdt dlyutAulyuu

40、mddddddd根据动量定理，动量变化等于作用在 流体上外力的冲量。这个外力就是作用在 上的水平方向的附加阻力 ，于是得mdAdFdlyutAtFdddddd式中 表示与X轴平行的流层之间作用在面积 上的总切力。则单位面积上的附加切应力为FdAdlyuAFtdddd（6-31） 2022-3-7工程流体力学假设脉动速度uyudd与时均速度的增量 成正比，即yukdd代入式（6-31），得到紊流的附加切应力2222dddd yulyulkt式中 22lkl 普朗特将 l称为混合长度，并认为它与 y成正比，即 kyl 式中 k比例常数，由实验确定 2022-3-7工程流体力学所以，紊流中的总切向应

41、力等于22ddddyulyut摩擦切应力t不同的，例如在接近管壁的地方黏性摩擦切应力起主要作用，等号右边的第二项可略去不计；在管道中心处，流体质点之间混杂强烈，附加切应力起主要作用，故可略去等号右边的第一项。的影响在有效截面上的各处是和附加切应力三、紊流结构、三、紊流结构、“光滑管光滑管”和和“粗糙管粗糙管” 1 1紊流结构分析紊流结构分析 由上节可知，黏性流体在管内作层流流动时，有效截面上的速度分布为抛物线分布。2022-3-7工程流体力学黏性流体在管中作紊流流动时，管壁上的流速为零，从管壁起流速将从零迅速增大，在紧贴管壁处一极薄层内，速度梯度很大，黏性摩擦切应力起主要作用，处于层流状态，称

42、为层流底层，距管壁稍远处有一黏性摩擦切应力和紊流附加切应力同样起作用的薄层，称为层流到紊流的过渡区；之后便发展成为完全紊流，称为紊流核心。如图6-15所示。层流底层的厚度在紊流水流中通常只有十分之几毫米。层流底层的厚度 可由下列两个半经验公式计算管道中 mm (6-33)Red8 .32明渠中 ReRh8 .32mm （6-34） 图6-15 紊流结构1层流底层；2过渡区；3紊流核心式中 管道直径，mm； 水力半径，mm； 沿程阻力系数 dhR2022-3-7工程流体力学从上式可以看出，层流底层的厚度取决于流速的大小，流速越高，层流底层的厚度越薄，反之越厚。层流底层虽然很薄，但是它对紊流流动的

43、能量损失以及流体与管壁之间的热交换起着重要的影响。例如层流底层的厚度越薄，换热就越强，流动阻力也越大。任何管子由于材料、加工、使用条件和年限等影响，管道内壁总是凹凸不平，其管壁粗糙凸出部分的平均高度 称为管壁的绝对粗糙度，而把 与管内径 的比值 称为管壁的相对粗糙度。常用管道绝对粗糙度见表6-1和表6-2。dd2022-3-7工程流体力学2022-3-7工程流体力学 2“光滑管”和“粗糙管” 从式（6-33）可知，层流底层的厚度 随着 的减小而增厚，当 时，则管壁的粗糙凸出的高度完全被层流底层所掩盖，如图6-16(a)所示。这时管壁粗糙度对流动不起任何影响，液体好象在完全光滑的管道中流动一样。

44、这种情况下的管道称为“水力光滑”管，简称为“光滑管”。 当 时，即管壁的粗糙凸出部分突出到紊流区中，如图6-16(b)所示。当流体流过凸出部分时，在凸出部分后面将引起旋涡，增加了能量损失，管壁粗糙度将对紊流流动发生影响。这种情况下的管道称为“水力粗糙”管，简称“粗糙管”。 在这里需要说明的是，对同一绝对粗糙度 的管道，当流速较低时，其层流底层厚度 可能大于 ，当流速较高时，其层流底层厚度 可能小于 ，因此同一根管道，在不同的流速下，可能是光滑管也可能是粗糙管。 Re2022-3-7工程流体力学图6-16 水力光滑和水力粗糙（a）“光滑管”；（b）“粗糙管”2022-3-7工程流体力学四、圆管中

45、紊流有效截面上的切应力分布和速度分布四、圆管中紊流有效截面上的切应力分布和速度分布 1切应力分布 紊流在半径 的管内流动，轴向时均速度为 ，切向应力在管长为 的管段上产生的能量损失，即压强损失 。若用管壁上的切向应力 来计算，则 (6-35) 如果在二个有效截面之间取半径为 （ ）的流管，则流管表面上切应力 可表示为 (6-36)0rulp0)(2212000pprlrlpr200r0rrlpr22022-3-7工程流体力学 因此，在有效截面上的切应力分布为 (6-37) 上式说明，紊流切向应力分布也与层流一样，与管半径 的一次方成比例，为直线关系，在 处切应力为零，如图6-17所示，从图中可

46、以看出，层流(a)的 与紊流(b)的 是不同的，两者的斜率不一样。 在紊流中切应力是指摩擦切应力和附加切应力，这两种切应力在层流底层和紊流核心所占比例不一样，在层流底层中，摩擦切应力 占主要地位，在紊流核心中附加切应力 占主要地位，根据对光滑管紊流实验，如图6-17(b)中的斜线部分为摩擦切应力，在 处附加切应力最大，当 摩擦切应力占主要，而在 范围内，摩擦切应力几乎为零，是以附加切应力为主的紊流核心区。rr00r0r00t095. 0rr095. 0rr07 . 0 rr2022-3-7工程流体力学图6-17 切应力分布（a）层流；（b）紊流2022-3-7工程流体力学 2速度分布在层流底层

47、（ ）中的切向应力为所以令 ，由于它具有速度的量纲，故称其为切应力速度，则有或 (6-38)由此可知，层流底层中的速度是按直线规律分布的。yyuyyuuyuu2*yuuu2022-3-7工程流体力学 在紊流区（ ）中假定切应力不变，令 （ 为管壁上的切向应力），则 常数或 (6-39)由式（6-39）可得积分得 (6-40)式中的积分常数 可根据层流底层与紊流区交界处（ ）的速度 相等的条件来确定，即 或 （6-41） （6-42） y0022220dd)(ddyukyyulkyukyyu*01ddkyyuuddCykuuln1*Cy*uuu2*)(uuCkuuln1*2022-3-7工程流体

48、力学由式（6-41）得或 （6-43）式中 层流底层的雷诺数， 将式（6-41）和（6-43）代入式（6-42）得 （6-44）将式（6-44）代入式（6-40）得再令 ，整理上式得 （6-45）uuReuuuuu*Reuu*ReuRe ln1*kuuC2*)(ln1ReuukReln1Reln1*kukRekReukykuuln1ln1ln1*RekReCln111*ln1Cyukuu2022-3-7工程流体力学尼古拉兹（Nikuradse）对光滑圆管中的紊流进行试验的结果得到： ， 。代入式（6-45）得 （6-46）式（6-46）即为圆管紊流速度分布的对数规律，此式只适用于光滑圆管。在圆

49、管的轴线处（ ）， ，代入式（6-46）得 （6-47）将式（6-47）与式（6-46）相减后得到 （6-48）40. 0k5 . 51C5 . 5log75. 55 . 5ln50. 2*yuuuyuuu0ry maxuu5 . 5log75. 55 . 5ln50. 2*0*max*0*maxuruuuruuyruuuyruuu0*max0*maxlog75.5ln50.22022-3-7工程流体力学式（6-48）称为普朗特公式。由于消去了常数项5.5，并经大量实验证明，此式对光滑管和粗糙管都适用。圆管紊流流速分布还可以近似地用一个简单的指数规律示之，即 （6-49）则平均流速 与最大流速

50、 之比，可由下式求得 即 （6-50）指数 n 随雷诺数 而变化，在不同指数n下的 与 的此值见表6-3。由表6-3知，当 =1.1105时， n=7。由式（6-49）则有 （6-51）这就是紊流的七分之一次方规律公式。nryuu10maxVmaxu0000maxmax0020d)(2d)(2rryyruuuyyrurV00010maxd)(2rnyyrryu) 12)(1(22maxnnnuVReVmaxuRe7/10maxryuu2022-3-7工程流体力学 表6-3中列出了平均流速 与最大流速 在不同雷诺数 下的比值。因而可用测定管轴处最大流速，用表6-3内的比值换算出平均流速，即可求出

51、流量。利用这种方法求管道有效截面上的平均流速及流量是非常简便的。 从以上分析可知，层流底层中的速度是按直线规律分布的，在紊流的核心区速度是按对数规律分布的，在核心区速度分布的特点是速度梯度较小，速度比较均匀，如图6-18所示，这是由于紊流时质点脉动掺混，动量交换强烈的结果 VmaxuRe2022-3-7工程流体力学2022-3-7工程流体力学 五、紊流流动中沿程损失的计算五、紊流流动中沿程损失的计算 式（6-11）也适用于对紊流流动沿程损失的计算，关键要确定紊流中的沿程阻力系数 。在一般情况下 ，即 值不仅取决于雷诺数 ，而且还取决于管壁相对粗糙度 ，情况比较复杂。紊流流动中的沿程阻力系数 的

52、计算公式，要在大量实验的基础上，对实验结果进行归纳分析，得出在不同条件下的经验公式，下节将详细讨论。/d)Ref，(Red图6-18 紊流速度分布2022-3-7工程流体力学第六节 沿程阻力系数的实验研究 层流流动的沿程阻力系数的计算公式已在第四节中用理论分析的方法推导出。由于紊流流动的复杂性，管壁粗糙度又各不相同，所以紊流流动的沿程阻力系数 值还不能与层流 一样完全从理论上来求得，而依靠对实验测得的数据进行整理归纳，得到经验公式。有许多学者和工程师做过 值的实验研究工作，在这类实验研究中，以德国尼古拉兹（JNikuradse）实验最有系统、范围最广，具有一定的代表性。 2022-3-7工程流

53、体力学 一、尼古拉兹实验一、尼古拉兹实验 各种管道的管壁都有一定的粗糙度，但管壁的粗糙度是一个既不易测量也无法准确确 定的数值。为了避免这个困难，尼古拉兹采用人工方法制造了各种不同粗糙度的圆管，即用漆胶将颗粒大小一样的砂粒均匀地贴在管壁上，砂粒直径表示管壁粗糙突出高度 。实验时采用砂粒直径 （即管壁的绝对粗糙度）与圆管半径 之比 表示以半径计算的管壁的相对粗糙度，用三种不同管径的圆管（25mm、50mm、l00mm）和六种不同的 值（15、30.6、60、126、252、507）在不同的流量下进行实验。对每一个实验找出沿程阻力系数且与雷诺数 和 之间的关系曲线。为了便于分析起见，将所有的实验结

54、果画在同一对数坐标纸上，以 为横坐标，以100 为纵坐标，并以 为参变数，即属于同一 的实验点用线连起来。 从6102106，包括层流在内，这个实验结果反映了圆管流动中的全部情况，如图6-19所示。现在将尼古拉兹实验曲线分成五个区域加以分析：rrrrrReReRer2022-3-7工程流体力学图6-19 尼古拉兹实验曲线2022-3-7工程流体力学 1层流区 当 2300时，所有六种不同的 的实验点都落在同一条直线上。这说明在层流流动时，沿程阻力系数 与管壁相对粗糙度 无关，而仅与雷诺数 有关，即 图6-19中的直线1恰好满足此方程，说明沿程损失 与有效截面平均流速 一次方成正比，实验进一步证

55、实了层流理论分析的正确性。 2层流到紊流的过渡区 2300 4000，当雷诺数超过2300时，流动状态开始发生变化，各种 的实验点离开1线，集中在一个很狭小的三角形区域内，这区域就是上、下临界雷诺数之间的不稳定区域，也就是层流到紊流的过渡区。rReRerRe）（ RefRe64fhVr2022-3-7工程流体力学 3紊流水力光滑管区 4000 59.6，各种不同 的实验点都落在同一倾斜直线2上，在这区域内沿程阻力系数 仍与相对粗糙度 无关，而仅与 有关，即 。这是由于层流底层的厚度还较大，足以掩盖粗糙突出高度 的影响，这区域就是紊流水力光滑管区。但是不同的 所占该直线上区段的长短也不同， 值越

56、小所占区段越短， 值越大所占区段越长。 =30.6的曲线几乎没有紊流光滑管区。这是由于在相同的雷诺数下，即在同样的层流底层厚度的情况下，较大的粗糙突出高度 先露出层流底层，变为水力粗糙管。ReRerrrr）（Refrr2022-3-7工程流体力学对于4103 105的这段直线2，勃拉休斯（H.Blasius）归纳了大量的实验数据，得出下列计算式 （6-52）在105 3106范围内，尼古拉兹结合普朗特的理论分析得到的公式为 （6-53）这就是光滑管的普朗特阻力公式，即图6-19中的曲线3。 若将式（6-52）代入式（6-11）中，可以证明沿程损失 与平均流速 成正比。ReRe25.03164.

57、0Re8.0)log(21Refh75. 1V2022-3-7工程流体力学 4紊流水力粗糙管过渡区 。随着雷诺数继续增加，各种相同 的实验点所连成的线先后进入区域4后部的5区域，所有的线都是平行于横坐标的直线，也就是说同一相对粗糙度的圆管有相同的 值，而与 无关，仅与相对粗糙度 （或 ）有关，即 ，这是因为此时层流底层的厚度已经非常薄，管壁粗糙度的作用已大大超过了层流底层内流体的黏性作用。 由式（6-11）可知，在水力粗糙区 值仅是 的函数，而同一 的圆管中 值是一个常数，沿程损失与平均流速的平方成正比，所以这个区域称为平方阻力区。 平方阻力区的 值可按尼古拉兹归纳的公式进行计算，即 （6-5

58、7）ReRer85. 0)(4160rrd）（ df rr2)/log274. 1 (r2022-3-7工程流体力学 由式（6-57）可知，在这区域中，要使两个流动的沿程阻力系数 值相等，只要使这两个流动（模型与实型）的相对粗糙度 相等即可，无需雷诺数 相等。因此紊流粗糙管平方阻力区又称为“自动模化区”，简称“自模区”。式（6-57）根据不同的 的计算结果列成表6-5。 以上介绍了尼古拉兹用人工粗糙度的管子所进行的实验。由此实验可知，流动在图6-19中不同的区域里，沿程阻力系数 值的计算公式是不同的。因此在计算沿程损失时，应先判别流动处在哪个区域，然后采用相应的公式去计算 值。rrRe2022

59、-3-7工程流体力学 二、莫迪图二、莫迪图 尼古拉兹的实验曲线是用各种不同的人工均匀砂粒粗糙度的圆管进行实验得到的，这 与工业管道内壁的自然不均匀粗糙度有很大差别。因此在进行工业管道的阻力计算时，不 能随便套用图6-19去查取 值。莫迪（F.Moody）根据光滑管、粗糙管过渡区和粗糙管平方阻力区中计算 的公式绘制了莫迪实用曲线，如图6-20所示。该图按对数坐标绘制，表示 与 、 之间的函数关系。整个图线分为五个区域，即层流区、临界区（相当于尼古拉兹曲线的过渡区）、光滑管区、过渡区（相当于尼古拉兹曲线的紊流水力粗糙管过渡区）、完全紊流粗糙管区（相当于尼古拉兹曲线的平方阻力区）。利用莫迪曲线图确定

60、沿程阻力系数 值是非常方便的。在实际计算时根据 和 ，从图6-20中查得 值，即能确定流动是在哪一区域内。 ddReRe2022-3-7工程流体力学图6-20 莫迪图2022-3-7工程流体力学【例例6-6】 输送石油的管道长 5000m，直径 250mm的旧无缝钢管，通过的质量流量 100t/h，运动黏度在冬季 =1.0910-4m2/s，夏季 =0.3610-4m2/s，若取密度 885kg/m3，试求沿程水头损失各为多少？ 解析解析【例【例6-7】 输送空气（输送空气（t=20）的旧钢管道，取管壁绝对粗）的旧钢管道，取管壁绝对粗糙度糙度 lmm，管道长，管道长 400m，管径，管径 25