第4章-流动阻力与水头损失

1、流体力学湖南城市学院土木院第4章流动阻力与水头损失u4-1. 流动阻力和水头损失的分类u4-2. 均匀流沿程水头损失与切应力的关系u4-3.粘性流体运动的两种流态u4-4. 圆管层流u4-5. 湍流运动的特点u4- . 边界层理论简介u4-6. 湍（紊）流的沿程水头损失u4-7.局部水头损失4-1. 水头损失的两种形式l沿程水头损失l局部水头损失l水头损失的计算l均匀流的沿程水头损失与壁面切应力的关系4-1. 水头损失的两种形式水头损失的概念：水头损失的概念：单位重量的液体自一断面流至另一断面所损失的机械能。沿程水头损失局部水头损失水头损失的分类：水头损失的分类：4-1. 水头损失的两种形式一

2、、沿程水头(阻力) 损失 hf1、 定义： 水头损失沿程均有并随沿程长度增加。主要由于液体与管壁以及液体本身的内部摩擦，使得液体能量沿程降低。2、特点：（1) 沿程粘性阻力均匀地分布在整个均匀流流段上；（2) 沿程损失与管段的长度成正比。 4-1. 水头损失的两种形式二、局部水头（阻力）损失 hj由于流动边界形状突然变化，引起流线弯曲以及边界层分离而引起的水头损失。 4-1. 水头损失的两种形式EXIT常见的发生局部水头损失区域 只要局部地区边界的形状或大小改变，液流内部结构就要急剧调整，流速分布进行改组，流线发生弯曲并产生旋涡，在这些局部地区就有局部水头损失。 4-1. 水头损失的两种形式

3、液流产生水头损失的两个条件液流产生水头损失的两个条件 (1) 液体具有粘滞性。 (2) 由于固体边界的影响，液流内部质点之间 产生相对运动。 液体具有粘滞性是主要的，起决定性作用。4-1. 水头损失的两种形式三、水头损失的计算公式4-1. 水头损失的两种形式1、沿程水头损失gvdlhf22gvRlhf242圆形管道：非圆管道：沿程损失因数R水力半径AR A-过流断面的面积-湿周 2、局部水头损失 hj= v2/(2g)局部损失因数4-1. 水头损失的两种形式3、总能量损失 hw jfwhhh4-1. 水头损失的两种形式圆管恒定均匀流为例圆管恒定均匀流为例4-2. 均匀流沿程水头损失与切应力的关

4、系流体段的受力为： （1）动压力： P1= p1A，P2=p2A ； （2）重力： G= g A l ； （3）摩擦力： F = 0 l F= 0 : P1P2 + G cos F=0 p1Ap2A + gA l cos 0 l = 0或流体段的平衡方程均匀流基本方程均匀流基本方程为为除以gA，并考虑lcos=z1-z2，并考虑能量方程 p1/gp2 /g+ (z1-z2)= 0 l /(gA)= 0 l /(gR)=hf4-2. 均匀流沿程水头损失与切应力的关系各流层间均有内摩擦切应力，同理可求得由均匀流基本方程以及沿程水头损失的达西公式可得壁面切应力与沿程损失因素的关系20v8定义摩阻流速

5、v*8v0*v达西公式（4-2）（4-9） 层流 湍(紊)流4-3. 粘性流体运动的两种流态4-3. 粘性流体运动的两种流态一、雷诺试验 层流：各层质点互不掺混紊流：流体质点的轨迹曲折、混乱，各流层的流体质点相互掺杂掺杂过渡流：层流与紊流之间的流动4-3. 粘性流体运动的两种流态二、沿程水头损失与流速的关系vmkhflglglgmfkvh 层流湍流0 . 10 . 11vvkhf0 . 275. 10 . 275. 12vvkhfab段层流ef段紊流be段临界状态 4512563156020 . 11m0 . 275. 12m0 . 23m结论：流态不同，沿程损失规律不同临界速度：流态转变时的

6、速度。 下临界速度：由湍（紊）流转变为层流时的速度vc 上临界速度：由层流转变为湍流时的速度vc实验证明， vc远小于vc通过正反两种实验情况，雷诺得出如下结果： 当v vc时，流体作紊流运动； 当v vc时，流体作层流运动； 当vc v vc时，流态不稳，可能是层流也可能是紊流。4-3. 粘性流体运动的两种流态三、雷诺数与其临界值三、雷诺数与其临界值 雷诺从一系列试验中发现：1、 不同种类液体在相同直径的管中进行实验，所测得的临界速度是各不相同的；2、 同种液体在不同直径的管中实验，所得的临界速度也不同。 故判定临界速度临界速度是液体的物理性质（ ， ）和管径（ d ）的函数。4-3. 粘性

7、流体运动的两种流态液体形态的判别液体形态的判别 雷诺数： 临界雷诺数：液流型态开始转变时的雷诺数。 对圆管： 对明渠及天然河道 ddRe2300Rec500ReRc4-3. 粘性流体运动的两种流态对于非圆管，如矩形、三角形、环形管等，管道的特征尺寸是管道的当量直径（或称水力直径），即： d=4A/= 4RR为水力半径，它是过流断面面积与湿周之比，即 R= A/EXIT例例 4-1:已知水流已知水流: t=5 , qV=15 L/s，d =0.1 m 求求: Re = ？解：由解：由t=5，查表得：，查表得： = 1.519 10-6 m2/s Re = vd / = 4 qV d / (d2

8、) = 4 qV / ( d ) = 1.2571054-3. 粘性流体运动的两种流态EXIT2）降低）降低 ： 由由Rec = vcd / 得：得： = d v / Rec = 0.08 0.02 / 2320 = 6.90 10-7 m2/s 由由 = 6.90 10-7 m2/s，查表得：，查表得： t=。所以，所以，将流速将流速 v 提高到提高到 0.132 m/s 或将温度提高到或将温度提高到，均可使流态变为湍流；，均可使流态变为湍流；例例 4-2: 已知水流已知水流: t=15 , v=0.08 m/s，d =0.02 m ， = 1.14 10-6 m2/s 求求: 临界临界Re

9、c = ？怎样使流态变为湍流？怎样使流态变为湍流解：解： Re = vd / = 1404 2320 所以巷道中的空气是湍流。所以巷道中的空气是湍流。4-3. 粘性流体运动的两种流态4-4. 圆管层流1.流动特性l 流体呈层状流动，各层质点互不掺混2.切应力层流中的切应力为粘性内摩擦切应力dydudrdu其中 y=r0- r3.断面流速分布drdu牛顿内摩擦定律又rdrgJdu2积分：（a a）)(4202rrgJu旋转抛物面2020max44rgJrgJurrurdrgJdu020并考虑边界条件 r=r0 ,u=04-4. 圆管层流max2040220218824)(vugJrArgJArd

10、rrrgJAudAAQ测量圆管层流平均速度的方法（b b）平均速度速度分布又可写为：202202max1v21rrrruu4-4.圆管层流4.沿程损失系数drdrdurrxv8v4)(000又由（4-9）知20v8Re64/v64v64dd表明：圆管层流的沿程损失因数与雷诺数的一次方成反比注意：v，但hfv 例例5-4 在长度 l=1000m,直径d=300 mm的管路中输送重度为=9.31kN/m 的重油，其重量流量G=2371.6 kN/h，求油温分别为10C（运动粘度为 =25cm2 /s）和40C (运动粘度为=15 cm2 /s)时的水头损失。 0708. 0360031. 96 .

11、2371GQ13 .040708.02AQvm3 /sm /s10C时的雷诺数20001202530100Re1vd40C时的雷诺数20002001530100Re2vd平均速度解解 体积流量4-4.圆管层流该流动属层流，故可以应用达西公式计算沿程水头损失。703.908 . 923 . 012011000642Re64222121gvdlgvdlhfm油柱高同理，可计算40C时的沿程水头损失m油柱高421.548 . 923 . 0200110006422fh4-4.圆管层流 湍流的特征 湍流的处理方法 湍流的切向应力 混合长度理论4-5. 湍流运动的特点4-5. 湍流运动的特点 湍流的成因

12、 湍流发生的机理是十分复杂的，下面给出一种粗浅的描述。高速流层低速流层 任意流层之上下侧的切应力构成顺时针方向的力矩，有促使旋涡旋涡产生的倾向。xy旋涡受升力而升降，产生横向运动，引起流体层之间的混掺一、湍流的特一、湍流的特征征 紊流的基本特征是许许多多大小不等的紊流的基本特征是许许多多大小不等的涡体涡体相互混掺前进，它们的位置、形态、相互混掺前进，它们的位置、形态、流速都在时刻不断地变化。流速都在时刻不断地变化。紊流实质上紊流实质上是是非恒定流动非恒定流动。二、紊流处理方法二、紊流处理方法时均值时均值 对随机的脉动，有两种处理方法：一为对随机的脉动，有两种处理方法：一为空间平均法空间平均法；

13、二为；二为时间平均法时间平均法。4-5. 湍流运动的特点试验研究结果表明：瞬时流速虽有变化，但在足够长的时间过程中，它的时间平均值是不变的。这个平均值称为时均速度，记作 x4-5. 湍流运动的特点TxxdtuTu01紊流可根据时均流动参数是否随时间变化，分为恒定流和非恒定流 即恒定流时时间平均流速不随时间变化。瞬时速度瞬时速度ux可表示为时均速度时均速度x和脉动速度脉动速度ux 的代数和代数和，即： ux = x+ ux 即即脉动速度脉动速度的的时间平均值时间平均值 x =0。同理同理 y =z =0。 010TxxdtuTu三、紊流的切向应力4-5. 湍流运动的特点平面恒定均匀紊流按时均化方

14、法分解成时均流动和脉动流动的叠加。相应的湍流运动的切应力 看作是由两部分所组成： 第1部分为时均层流引起时均层流引起相邻两流层间相对运动产生的粘性切应力1 dyudx1第2部分为由脉动，上下层质点相互掺混产生的附加切应力（雷诺应力）yxuu2靠近壁面且Re数较小时， 1占主导地位离开壁面且Re数很大时， 2 1四、混合长度理论4-5. 湍流运动的特点普朗特认为：存在一个称之为混合长度的距离 l，质点在此距离内不与其他质点相碰，保持原有的物理属性，直至经过此行程l，才与周围质点掺混，湍流脉动速度与混合长l与（时均）速度梯度的乘积成正比，即dyudluudyudluxxyxx21,将l 1和l 2

15、合并为l2 ，湍流附加切应力为222)(dyudlxyxuu2代入（4-31）374-5. 湍流运动的特点混合长l，不受黏性影响，只与质点到壁面的距离有关，即yl（4-32）充分发展的紊流中，黏性切应力远小于雷诺应力，并认为壁面附近切应力一定，为壁面切应力 0，由（4-31）2222220)()(dyduydyudlxydyvduvdyduy*0上式变形，并根据摩阻流速的定义（P85），可得积分上式，其中 0一定，摩阻流速v*为常数cyvuln1*（4-33）普朗特-卡门对数分布律 边界层概念 边界层的两种流态 边界层分离现象 物体的绕流阻力4- . 边界层理论4- . 边界层理论 普朗特在1

16、904年的国际数学家大会上首次提出了边界层的概念。 ux= 0.99U的点规定为边界的点规定为边界层的外边界。层的外边界。边界层的厚度边界层的厚度ux=f（y） 在紧贴着物体的表面，流动受到粘性的显著影响，流速沿壁面法向的变化非常急剧，粘性切应力粘性切应力不能略去不计的极薄的一层流体，称为边界层。一、边界层概念4- . 边界层理论边界层厚度（1）物体边界附近薄层由于粘性力作用，有很大的速度梯度du/dy边界层（附面层）；（2）边界层以外的流动，粘性力作用不计理想流体 边界层中的流动也存在两种流态，从前缘起自层流开始，随x 增加，边界层越来越厚，壁面对扰动的稳定作用逐渐减弱，直至发生流态的转捩，

17、从层流转变为紊流。二、边界层的两种流态4- . 边界层理论xc的值与流体的运动粘度成正比，与边界层外边界上的流体速度成反比4- . 边界层理论UxccReRec 转捩雷诺数，通常为31053106。三、三、粘性底层（层流底层）尽管边界层内已是紊流流态，由于壁面的强烈作用，在紧贴壁面的一层中，流动仍保持层流的性质，称为粘性底层或层流底层。粘性底层厚度厚度约为xUx0Re实验表明，平板边界层内雷诺数的表达式Re9 .32vvv6 .11v6 .11*0ddd（4-26） 4- . 边界层理论粘性粘性底层的速度分布底层的速度分布 粘性底层属于层流层流流动，湍流附加切应力为零，流体受到的切应力只有粘性

18、切应力，可用牛顿内摩擦定律表示。又因为粘性底层很薄，可认为薄层内切应力为常数0dydux积分，并由边条y=0，ux=0，得到yux0或写成)v(vv0*摩阻速度yux（4-28）（4-29）黏性底层中，速度按线性分布，壁面速度为0四、边界层的分离现象4- . 边界层理论壁面最高点边界层分离点边界层边界层内边界分离区（回流区或尾流区）逆压区顺压区uA uC pA pB pC B为壁面最高点，该处过流断面最小，流速最大，压强最低。1、流体从A向B流动时，压强降低，压能减小，转化为动能，流体速度沿程增大，为顺压梯度区；2、流体从B向C流动时，流体速度沿程减小,压强增大，为逆压梯度区；3、在C点处出现

19、粘滞 ，流体速度的法向梯度为零，壁面切应力降为零。在C的下游流体发生回流。回流导致边界层分离，并形成尾涡。 C点称为分离点。4- . 边界层理论4- . 边界层理论 物体所受阻力表示为：AUCFDD220CD：阻力系数：流体密度U0：来流速度A：迎流投影面积4- . 边界层理论四. 物体的绕流阻力 4.6 湍流的沿程水头损失 一一. . 尼古拉兹实验尼古拉兹实验前面已经给出了沿程水头损失hf的计算公式gvdlhf22（4-2）圆管层流时，仅与Re有关： =64/Re由于紊流的复杂性，未能严格从理论上推导出 的理论公式尼古拉兹通过实验测定，认为 有两个影响因素)/(Re,df绝对粗糙度：管壁粗糙

20、突起的高度。 d 管径/d相对粗糙度eRlg)100lg(过渡粗过渡粗糙区糙区1.530.6r0 /601262525073.06.00.25.04.00.40.60.81.0粗糙区粗糙区过渡区过渡区光滑管光滑管层流区层流区圆管流动沿程水头损失系圆管流动沿程水头损失系数的数的尼古拉兹试验曲线尼古拉兹试验曲线 4.6 沿程损失因数的变化规律 Re2300，层流区层流区， =f (Re)=64/Re Re=23004000，层流向紊流过渡区层流向紊流过渡区， = f (Re)，该，该区范围很窄，实用意义不大。区范围很窄，实用意义不大。 4000Re 80d/，紊流光滑区紊流光滑区， =f (Re)

21、. 沿程损失沿程损失 系数仅与雷诺数有关。系数仅与雷诺数有关。 80d/Re 1140d/紊流粗糙区紊流粗糙区， =f ( /d). 沿程损失系数沿程损失系数仅与粗糙度有关。沿程水头损失将与平均流速的平方仅与粗糙度有关。沿程水头损失将与平均流速的平方成正比，通常也叫做成正比，通常也叫做阻力平方区阻力平方区。沿程损失系数的五个分区沿程损失系数的五个分区1区区2区区3区区4区区5区区 4.6 沿程损失因数的变化规律518vv*由（4-9）和（4-10）可知（4-10） 4.6 沿程损失因数的变化规律紊流光滑紊流光滑管：管：紊流粗糙紊流粗糙管：管： (Re), 4 . 0/80Re0 . 5v(5

22、. 5vln5 . 2v0*fdyux或或)/(, 6/1140Re70v(5 . 8ln5 . 2v0*dfdyux或或根据尼古拉兹实验)/(Re,1140Re80dfdd过渡粗糙区：过渡粗糙区： cyvuln1*为获得平均流速，可积分紊流流速分布，考虑到不同边界条件52将上两式写为一个通式Byux0*rln5 . 2v紊流光滑时0*vln5 . 25 . 5rB0ln5 . 25 . 8rB紊流粗糙时 4.6 沿程损失因数的变化规律dAByAAAdAAAAx0*rln5 . 21vvvu175. 1vln5 . 275. 3vv0*rB通过积分可得紊流光滑时0*vln5 . 25 . 5r

23、B0ln5 . 25 . 8rB紊流粗糙时8vv*（4-10）比较公式（4-10）EXIT 4.6 沿程损失因数的变化规律8 . 0)lg(21eR4000Re80d/ /紊流光滑紊流光滑区区紊流粗糙紊流粗糙区区74. 1lg210rRe 1140d/ /柯列勃洛克将两式改写为d7 .3lg21Re51.2lg21并将两式合并在一起)Re51. 27 . 3lg(21d此式是水力过渡粗糙湍流沿程损失因数的经验公式，称为柯列勃洛克公式层流区 =64/Re水力粗糙区 Re 1140d/-1/2 = （2 lg (d/(2) + 1.74）水力光滑区4000 Re 80d/ -1/2 = 2 lg

24、(Re1/2) 0.8 （普朗特公式）4000 Re 2300层流层流Re 2300)(4220rrJgux208rJgveR64光滑管区光滑管区5eR705eR70eR过渡粗糙管区过渡粗糙管区粗糙管区粗糙管区)5 . 8ln5 . 2(uyux)5 . 5uln5 . 2(yuux)75. 1ln5 . 2(0uruv8 . 0)lg(21eR2074. 1lg21r4/ 13164. 0eR)10(5eR流速分布流速分布圆管流动圆管流动主要公式主要公式),(0reR)75. 4ln5 . 2(u0rvEXIT三、三、莫迪图莫迪图层流层流层流区层流区过渡区过渡区粗糙区粗糙区过渡粗过渡粗糙区糙

25、区光滑管光滑管vdRe/d知: l =1000 m , d= 0.2 m， = 3.55 10-5 m2/s，Q = 3.810-2 m3/s, = 0.039 mm 求: hf =？解：v = 4qv /( d2)= 1.21 m/s Re = vd / = 6817 80(d/) = 80 0.2/0.000039 = 410256 因： 4000 Re 80(d/) ，所以为水力光滑管，且 又因为满足：4000 Re ，的确为水力光滑管，计算结果正确。例4-6 已知： l =30 m , d= 0.75 m， = 1.5710-5 m2/s，= 1.2 mm， Q= 30000 m3/h

26、 = 25/3 m3/ s 求： hf=？ 4.6 沿程损失因数的变化规律 解：v = 4Q/( d2)= 18.9 m/s Re = vd / = 9.03 105 1140(d/) = 1140 0.75/0.0012 =7.125 105 因： Re 1140(d/) ，所以为水力粗糙区。由/d=0.0012/0.75=0.0016及Re = 9.03 105查穆迪图得： =0. 022由公式： -1/2 = 2 lg (d/(2) + 1.74求得： = 0.022 hf= (l/d)v2/(2g)= 0.022(30/0.75)18.92/(2g)= 16.04m8:23:08例4-

27、7: 已知: l =300 m , d= 0.2 m， = 900 kg/m3，流量G= 245 N/s, 1 = 1.092 10-4 m2/s (冬季 )，2 = 2.10 10-5 m2/s (夏季 ) 求: hf1=？, hf2=？ 解：v = 4qv/( d2)= 4G/( g d2)= 0.884 m/s Re1 = vd / 1 = 1619 Re2 = vd / 2 = 8419 Re1 = 1619 2320，是湍流运动，但4000Re105 ，是水力光滑区，由布拉休斯公式得： =0.3164/Re1/4 =0.3164/ 84191/4 = 0.033 故： hf2 = 2

28、(l/d)v2/(2g) = 0.036(300/0.2)0.8842/(2g)= 1.97 m8:23:08 4.6 沿程损失因数的变化规律四、明渠流动的谢才公式RJC谢才公式谢才系数明渠均匀流的水力坡度即为水面线坡度C水力半径 RJ曼宁公式沿程水头损失说明 相当于我们定义的 82gCn是边界粗糙系数（糙率），其量纲不甚明确，查表后将相应数值代入公式即可。 8:23:094.7 局部水头损失 局部水头损失形成机理 局部水头损失系数 突扩圆管的局部水头损失系数 突缩和其他局部水头损失系数8:23:09 局部水头损失形成机理突然扩大突然缩小闸阀三通汇流管道弯头管道进口分离区分离区分离区分离区分离

29、区分离区分离区 有压管道恒定流遇到管道边界的局部突变 流动分离形成剪切层 剪 切 层 流 动 不 稳定，引起流动结构的重新调整，并产生旋涡 平均流动能量转化成脉动能量，造成不可逆的能量耗散。4.7 局部水头损失8:23:09v1A1A2v21122根据能量方程 认为因边界突变造成的能量损失全部产生在1-1，2-2两断面之间，不再考虑沿程损失。局部水头损失 与沿程因摩擦造成的分布损失不同，这部分损失可以看成是集中损失在管道边界的突变处，每单位重量流体承担的这部分能量损失称为局部水头损失局部水头损失。8:23:09v1A1A2v21122 上游断面1-1取在由于边界的突变， 水流结构开始发生变化的

30、渐变流段中，下游2-2断面则取在水流结构调整刚好结束，重新形成渐变流段的地方。总之，两断面应尽可能接近，又要保证局部水头损失全部产生在两断面之间。经过测量两断面的测管水头差和流经管道的流量，进而推算两断面的速度水头差，就可得到局部水头损失。 4.7 局部水头损失EXITv1A1A2v21122局部水头损失细管平均流速v1取1-1，2-2两断面如图两断面面积都为A2粗管平均流速v2 一.突扩圆管的局部水头损失系数理论结果 F=p1A1p2A2+ p1(A2A1) + glA2 cos Q = A1v1= A2v24.7 局部水头损失EXIT动量方程v1A1A2v21122代入局部水头损失的表达式020121,均取为1.0EXIT4.7 局部水头损失EXIT 突扩圆管局部水头损失之所以能够导出上述解析表达式是因为： 我们假设1-1断面上的测管水头为常数； 旋涡区全部在粗管中，所以1-1，2-2两断面都取在粗管中，它们的面积相等。EXITv1A1A2v211224.7 局部水头损失 当流体在淹没出流条件下