第四章不可压缩粘性流体一元流动优秀PPT.ppt

第四章不可压缩粘性流体一元流动第一页，本课件共有61页第四章第四章 不可压缩粘性流体一元流动不可压缩粘性流体一元流动实际流体具有粘粘性性，粘性的存在会使流体内部及流动壁面处产生粘粘性性阻阻力力，流体运动为了克服这部分阻力必然要使流动流体的部分机械能不可逆转地转化为热能，造成能量损失。这一能量损失就是实际流体伯努利方程中的hw这种能量损失也可以压力损失的形式（g hw）来表示。管内流动的能量损失（水头损失）hw一般可分为两部分：沿程阻力损失hf：局部阻力损失hj：第二页，本课件共有61页沿程阻力损失沿程阻力损失：当流体沿管径不变的管道流动时，此时过流断面上的速度分布沿流程不变，粘性切应力也沿流程不变，管道沿程越长，能量损失越大。这种能量损失称为沿程阻力损失hf。沿程阻力损失可用下式计算：上式也称为达西（达西（DarcyDarcy）公式）公式。称为沿程阻力损失系数，其与流体的流动状态有关局部阻力损失：局部阻力损失：当流动边界发生急剧变化（如流体流经弯管、管截面突然扩大或缩小、各种阀门、流量计等等）情况下，流线发生变形，产生旋涡，这也会造成部分机械能的损耗。这种由管路局部原因造成的能量损失（水头损失）称为局部阻力损失hj。局部损失可用下式计算：称为局部阻力损失系数。的大小一般由实验测的，不同不同形式的局部损失，其大小是不同的。具体计算时可查有关手册。第三页，本课件共有61页对于一个管路系统，一般都有若干管段及各种管件组成。故管路系统总的能能量损失量损失为：研究阻力损失的目的：研究阻力损失的目的：一、正确计算管路或液压系统中的流动阻力二、找出减小流动阻力的途径。以下将介绍不同流动情况下管路内压力损失的计算，为管路系统水力分析、设计、节能打下理论基础。第四页，本课件共有61页4.1流体的流动状态流体的流动状态流体流动的能量损失与流体的流流动动状状态态密切相关。因此在讨论流体流动能量损失之前必须对流体的流动状态进行分析。4.1.1 雷诺实验雷诺实验1883年英国物理学家雷诺通过实验证明流体流动存在两种不同的流动状态，确定了流态的判别及其与能量损失的关系。以下介绍一下雷诺的实验过程及结论。实验配置：1：盛有无色水的大容器，液位高保持不变。2：透明等直径管。3：盛有有色水的小容器。4：阀门。5、6：两个测压管。第五页，本课件共有61页实验过程：实验过程：（1）打开阀门4，管内流速由小变大，当水以较低速度流过透明管。管中有色水近似呈直线状流动，说明管中水流的质点以一种有规则的、互不混杂的形式作分层流动。这种流动状态称为层流层流。（2）继续增大管内流速，在一定范围内，管内仍保持为层流。但到流速增大到一定值时，有色水开始波动。说明此时流体质点已出现横向运动。第六页，本课件共有61页（3）近一步增大流速，有色水破裂，完全混杂在周围的无色水中。说明此时流动质点呈现一种相互混杂的、无规则的流动。这种流动状态称为紊紊流流。继续增大流速，管内流动将一直处于紊流。（4）由紊流状态，逐渐减小流速V，流动状态将由紊流又转换为层流。但实验发现：流动由紊流转换为层流时的临界速度与层流转换为紊流时的临界速度并不相同。前者称为下下临临界界速速度度（紊紊流流-层层流流）；后者称为上上临临界界速速度度（层层流流-紊紊流流）。第七页，本课件共有61页 流动状态的判别：流动状态的判别：前面实验表明：流体流动存在两种不同的流动状态：层层流流和和紊紊流流。当流速达到一定值时（临界速度），两种流态将会互相转换。雷诺用不同的管径、不同的流体进行了实验，结果表明：流动状态实际与一个无量纲数，即雷诺数Re有关：即雷诺数Re才是判别流态的准则。临界雷诺数也有上、下临界之分：上临界雷诺数：上临界雷诺数：对应于上临界速度下临界雷诺数：下临界雷诺数：对应于下临界速度。上临界雷诺数约为12000-40000，下临界雷诺数较稳定，约为2320。工程上一般认为临界雷诺数Recr=2320。即：层流层流紊流紊流第八页，本课件共有61页沿程损失沿程损失hf与流态的关系：与流态的关系：层流和紊流情况下，管内沿程损失hf的变化规律是不同的。层流时hf与流速V的一次方成正比；而紊流时hf与流速V的1.75-2次方成正比。如下图所示：Vcr，Vcr分别为上、下临界速度。第九页，本课件共有61页4.2 圆管中的层流圆管中的层流4.2.1截面速度分布截面速度分布粘性流体在圆管内作层流运动时，管截面上的速度分布是不均匀的。在管壁上，由于流体的粘附作用，流体速度为零，越接近管轴线，速度越大。在轴线上，速度达到最大值。以下推导一下层流速度分布：在一个足够长的等截面圆管中，取一个以管轴线为中心的流体圆柱，半径为r，长为L。由于管道足够长，任意两个截面的速度分布相同。如果流动为定常，则加速度为0，作用在圆柱体上流体的外力（表面力和质量力）平衡，即：第十页，本课件共有61页或可写为：由上式可看出，切应力与r成正比。在管壁上，r=d/2，切应力最大：由伯努利方程可知，管流的沿程水头损失：又由达西公式：比较以上三式可得：第十一页，本课件共有61页上式V为截面平均速度。引入摩擦速度u*：式（4-10）可写为：注注意意：上面各式均没有涉及到流态。故式（4-12）、式（4-10）对层流和紊流均可适用。如果流体作层流运动，式（4-6）中的粘性切应力可由牛顿内摩擦定律表示为：上式中的负号是由于，沿径向r方向速度是减小的。将式（4-13）代入式（4-6），并积分可得：第十二页，本课件共有61页当r=r0时，速度u=0，由此可求得上式积分常数C=-r02。速度分布可确定如下：上式表明：圆管内层流时，速度分布为旋转抛物面。最大速度在管轴线上，并且为：由截面速度分布，还可计算出通过管截面上的体积流量Q。第十三页，本课件共有61页管截面平均速度V：故截面速度分布也可用平均速度来表示：通过以上速度分布，还可计算层流时的动能修正系数及动量修正系数。第十四页，本课件共有61页4.2.2 沿程阻力损失系数沿程阻力损失系数达西公式：由式（4-18）得：将式（4-21）代入达西公式（4-20）得沿程损失阻力系数：沿程损失阻力系数：上式即为圆管内层流时的沿程阻力系数的计算公式。由上式计算出，再代入达西公式即可计算出圆管内层流时的沿程水头损失hf（或沿程压力损失pf=g hf）。第十五页，本课件共有61页对于非圆截面管道，达西公式仍然适用，但此时应用水力半径R来计算。水水力力半径半径R的定义如下：式中A：管道截面积；S：湿周（流体湿润的周界长度）对于圆管R=d/4，故非圆管道的达西公式为：相应的雷诺数：第十六页，本课件共有61页4.3 圆管中的紊流圆管中的紊流4.3.1 紊流特征紊流特征（1）流体作紊流流动时不同流层的流体质点相互混杂，作无定向、无规则的运动。表现在运动参数（如速度、压强、密度等）上具有随随机机性性变变化化。即紊流的主要特征是在时间与空间上具有脉动性脉动性。（2）流体的物理参数虽然是随机脉动的，但又服从统计规律。数学处理上可用统计统计平均平均方法。即：上式中u：某点的瞬时速度。：平均速度。：脉动速度。第十七页，本课件共有61页类似上式，对于其他物理量同样有：对于上面物理量的平均处理方法，主要有时间平均法、空间平均法、概率平均法等。其中常用的为时间平均法时间平均法。时间平均法（时均法）定义时间平均法（时均法）定义：上 式T应足够长，以获得比较稳定的平均值。第十八页，本课件共有61页（3）紊流具有脉动速度（u，v），脉动速度会引起紊流的动动量量交交换换。动量交换的结果又会产生紊流附加切应力紊流附加切应力。以下简单说明紊流附加切应力产生的原因及表达式：取控制体，设控制体底面积为A，由于脉动速度，下面流体会进入控制体，并代入一定动量。根据动量方程，作用在控制体上外力为动量流出动量为0，故：第十九页，本课件共有61页上式两边同除A，并取时均值，紊流附加切应力为：因此，紊流时流体的切应力除了有粘粘性性切切应应力力，还有因脉动速度（u，v）引起动量交换而产生的紊流附加切应力（雷诺应力紊流附加切应力（雷诺应力）。即：对于层流，没有脉动速度，只有粘性切应力。在紊流中，脉动速度引起的动量交换非常强烈，紊流切应力要远大于粘性切应力。第二十页，本课件共有61页4.3.2 普朗特混合长理论及紊流附加切应力普朗特混合长理论及紊流附加切应力式（4-29）虽给出了紊流附加切应力的表达式。但脉动速度为多大，它与平均速度又有何关系，这些至今还不能以严格的数学推理给以确定。为此许多学者提出许多假设。在工程上目前以普普朗朗特特提出的半经验理论-混混合合长长度度理理论论应用最为广泛。故我们就主要介绍一下普朗特的混合长理论。普朗特混合长理论：普朗特混合长理论：为了确定脉动速度，普朗特认为流体质点在y方向的脉动，即由一层跳入另一层，要要经经过过一一段段不不与与其其他他流流体体质质点点相相碰碰撞撞的的距距离离l。然后以自己原来的动量和新位置周围的流体质点相混，完成动量交换。流体质点从一层跳入另一层，并且不于其他质点相碰的这段距离l称为混混合合长长度度。它是流体质点横向混杂运动中自由行程的平均值。第二十一页，本课件共有61页在上述基本假设的基础上，普朗特还作出以下假设：（1）流体质点的纵向脉动速度u近似等于两层流体的时均速度之差。即（2）横向脉动速度与纵向脉动速度成比例。即：根据以上假设，普朗特获得了脉动速度与时均速度的关系，进而紊流附加切应力可表达为：为了书写简便，上式时均速度 ，后面用u来代替。第二十二页，本课件共有61页4.3.3 紊流速度分布紊流速度分布（1）紊流结构及粘性底层）紊流结构及粘性底层图示为不同雷诺数下管截面速度分布。层流时速度分布不均匀，紊流时分布均匀。雷诺数越大，速度分布越趋向于均匀。第二十三页，本课件共有61页对于紊流，由于管壁限制，紧贴管壁的流体速度为0，在靠近管壁很薄的一个流层内，流体速度比较小，流体的紊动程度几乎为0，流动处于层流，其切应力主要表现为粘粘性性切切应应力力。这层流体称为粘粘性性底底层层（近近壁壁层层流流区区）。在远离壁面的很大区域内，流速较大，并且速度分布比较均匀，称为紊紊流流核核心心区区。粘性底层与紊流核心区之间还有一个过渡区域，一般将过渡取也划入紊流核心区。如下图所示。第二十四页，本课件共有61页粘性底层的厚度非常薄（大约只有几分之一或几十分之一毫米），其厚度可用下列经验公式计算：或：两个概念：水力光滑管水力光滑管与与水力粗糙管水力粗糙管实际管壁是凹凸不平的，管壁粗糙物凸出的平均高度称为壁壁面面绝绝对对粗粗糙度糙度。/d则称为相对粗糙度。当 ，称为水力光滑管；时，称为水力粗糙管。第二十五页，本课件共有61页（2）速度分布）速度分布以上表明，管内紊流结构可分为粘性底层与紊流核心区，两个区的速度分布是不同的，因此应分别讨论。对对于于粘粘性性底底层层：紊流切应力很小，只有粘性切应力对流动起作用。实验表明：在很薄的粘性底层，切应力变化不大，可近似认为常数0（壁面上的剪切应力）。根据牛顿内摩擦定律，有：引入摩擦速度u*（0=u*2），代入上式积分，得：或：上式表明，粘性底层的速度呈线性分布。第二十六页，本课件共有61页对对于于紊紊流流核核心心区区：紊流切应力远大于粘性切应力，故粘性切应力可略去。实验还表明：紊流切应力的大小与壁面处的粘性切应力近似相等。故：普朗特还认为，混合长l与距壁面的距离y成正比，即：式（4-36）可变为：上式积分，得：第二十七页，本课件共有61页上式中k称为卡卡门门（Von Karman）常常数数，k=0.4。式中的积分常数C可这样求得：当y=时，用式（4-38）与前面粘性底层速度分布（4-35）求得的速度相等，从而可求得积分常数。故式（4-38）变为：根据大量实验得出B=5.5，而k=0.4，故水力光滑管速度分部为：可见，紊流区速度分布比粘性底层速度分布要平缓的多。利用上述速度分布，还可计算出管截面平均速度V。由于粘性底层很薄，计算流量时可只计算紊流区流量。将式（4-40）速度代入式（4-41），并注意：第二十八页，本课件共有61页得：该对数形式速度分布（4-42）比较复杂，应用起来不太方便。普朗特根据实验，给出了一种指指数数形形式式的的速度分布速度分布：指数n的大小与雷诺数有关，如下表Re41032.31041.11051.11062106n1/61/6.61/71/8.81/10第二十九页，本课件共有61页4.3.4 圆管内紊流沿程损失系数圆管内紊流沿程损失系数（1）水力光滑管）水力光滑管：利用上式（4-42）及式（4-12），并进行对数换底，得：根据实验，将上式修正为：上式即为紊流水力光滑管的紊流水力光滑管的-Re的关系式的关系式，也称卡门卡门-普朗特方程普朗特方程。第三十页，本课件共有61页（2）水力粗糙管）水力粗糙管水力粗糙管的速度分布与管壁粗糙度有关，其平均速度分布：同样方法可得：上式表明：水力粗糙管的沿程损失系数上式表明：水力粗糙管的沿程损失系数与雷诺数与雷诺数Re无关。无关。（3）过渡区（水力光滑与水力粗糙）过渡区（水力光滑与水力粗糙）对于水力光滑与水力粗糙的过渡区域，沿程损失系数可用下列柯列勃洛克（C.F.Colebrook）公式计算。由上式可见，在该区域内，沿程损失系数与雷诺数及管壁粗糙度均有关系。由上式可见，在该区域内，沿程损失系数与雷诺数及管壁粗糙度均有关系。第三十一页，本课件共有61页上面各公式中均涉及到管壁粗糙度，它与管道种类有关，具体计算时可参照下表来选取：常用工业管道当量粗糙度常用工业管道当量粗糙度管材种类（mm）新聚氯乙烯管、玻璃管，铜管、铅管0.0015-0.01新无缝钢管0.014旧无缝钢管0.2新焊接钢管0.06旧焊接钢管1.0新铸铁管0.3旧铸铁管1.0-1.2新镀锌钢管0.15旧镀锌钢管0.5水泥管0.5第三十二页，本课件共有61页工程设计中，求取沿程水头损失（或压力损失）的方法：工程设计中，求取沿程水头损失（或压力损失）的方法：达西公式：其中沿程损失系数其中沿程损失系数的确定方法可归纳以下两种：的确定方法可归纳以下两种：1.通过前面的公式来计算。即（1）层流区：（2）水力光滑区：在该区域，假如雷诺数不太大（Re=4000-105），还可用简单的勃拉修斯公式计算：第三十三页，本课件共有61页（3）过渡区：（4）水力粗糙区：柯列勃洛克（C.F.Colebrook）公式：第三十四页，本课件共有61页计算可用牛顿迭代法牛顿迭代法或其他迭代法。牛顿迭代法用于求解超越方程f(x)=0：第三十五页，本课件共有61页2.查表法：查表法：在有关工程应用手册给出了相关的，Re之间的关系图表（穆迪图）可根据实际的/d，Re直接在图表中查出的大小。第三十六页，本课件共有61页4.4 局部水头损失局部水头损失流体流经固体壁面急剧变化的部位（截面突然扩大、突然缩小，弯管，阀门等）时，流体微团相互碰撞和产生旋涡，引起该局部区域较大的机械能损失，称为局部损失局部损失。局部损失产生的原因是：（1）边界面急剧变化使主流与边界面间形成旋涡区，流体质点相互撞击；（2）产生与主流方向相正交的流动，即二次流。以下为几种局部损失：第三十七页，本课件共有61页弯管：流经弯管时的局部损失：流体流经弯管时，不但会产生旋涡区，而且会产生二此流。流体进入弯管时，流线发生弯曲，流体质点由于受到离心力而压向外壁面，使外壁面压力升高，流速减小。大约在B点以后压力才逐渐降低直至弯管终止为止。与此对应，在管内侧将出现压力降低随后升高的现象。因此在相应的AB，BC区域都会因升升压压减减速速脱离壁面而形成旋涡区，造成旋旋涡涡损损失失。流经弯管的另一个局部损失的原因是二次流。由于粘性的作用使管壁附近的液流速度低于内部，速度差将造成中心与外侧离心力的不同，其结果是在离心力差值的作用下使流体质点压向外侧，形成图示的自中心向外侧的二此流。径向方向的二次流与主流结合会形成双螺旋流动。这一流动情况在紊流情况下将大于100倍管径。弯管曲率半径小，内外压差将增大；管子直径大，二次流的范围会增加，其结构都将使局部能量损失增大。第三十八页，本课件共有61页造成局部损失的形式是各种各样的，但起产生阻力、形成能量损失的本质基本相同。故对各种局部水头损失，均可用下面形式计算：上式中：局部损失系数。对于不同的局部损失，其大小也是不同的。目前对于，除了少数几种情况下可理论推导外，其他形式局部损失系数均需实验获得。（1）截面扩大时的）截面扩大时的由截面1-1、2-2建立伯努利方程，动能修正系数取为1：可得：第三十九页，本课件共有61页为求上式中的p1-p2，还必须利用动量方程。取1-1、2-2截面间区域为控制体。建立动量方程：即：将上式代入前面局部损失，得：上面局部损失水头可用上游速度计算，也可用下游速度计算。即：第四十页，本课件共有61页（2）截面突然缩小的）截面突然缩小的（3）渐扩圆管的）渐扩圆管的8o10o12o15o20o25ok0.140.160.220.30.420.62第四十一页，本课件共有61页（4）圆管弯头的）圆管弯头的第四十二页，本课件共有61页（5）圆管进口的）圆管进口的1：容器上的圆管进口：2：容器上的圆管平滑进口：3：容器上的内伸进口：/Db/D00.0020.010.050.500.50.570.630.81.00.0080.50.530.580.740.880.0160.50.510.530.580.770.0240.50.50.510.530.680.030.50.50.510.520.610.050.50.50.510.50.53第四十三页，本课件共有61页（6）闸阀的）闸阀的h/d全开6/84/83/82/81/80.110.262.065.5217.097.8第四十四页，本课件共有61页(8)渐缩圆管渐缩圆管10o20o40o80o100o140ok0.40.250.20.30.40.6第四十五页，本课件共有61页（9）折管（急剧弯管）：折管（急剧弯管）：（10）孔板）孔板第四十六页，本课件共有61页（11）流入静止大容器的管路出口）流入静止大容器的管路出口：这种情况可看作为截面突然扩大的极限情况。即：（12）截止阀全开：）截止阀全开：（13）蝶阀全开：）蝶阀全开：其他阀可参阅有关手册。第四十七页，本课件共有61页减小局部损失的措施减小局部损失的措施：对于实际的装置，只能根据其具体的结构条件来确定减阻措施。（1）对于管截面需要变化的位置，如果条件许可，可以用截面逐渐扩大或截面逐渐缩小来代替突然扩大或突然缩小。（2）容器上的圆管进口尽量改为平滑进口。（3）在可能的情况下，尽量增大阀门的开度。（4）对容器上的内伸进口，尽量缩短内伸长度。第四十八页，本课件共有61页（5）对弯管（如蒸汽或风送管道），弯曲部分将产生较大的局部阻力。可安装适当形式的导流板。这样可避免弯曲部分的内、外侧出现大范围的旋涡区，也可减小二次流的产生和影响范围。实验证明：实验证明：不装导流板时直角弯头的损失系数安装月牙形导流板时的损失系数：总总的的原原则则：尽量减小旋涡区，减小二次流的产生及波及范围，从而减小流体质点的撞击损失和减小速度重新分布时的动量交换。第四十九页，本课件共有61页4.5 工程应用工程应用（1）孔口出流）孔口出流恒定水位的容器，容器侧面的小孔面积为A。求出流流量。流体流出后收缩，在喉部的截面积为Ac称为面积收缩系数。由液面0-0，及截面c-c建立伯努利方程。第五十页，本课件共有61页流量为：上式中称为流量系数，显然：实验测得局部损失系数：面积收缩系数：因此流量系数：第五十一页，本课件共有61页（2）管嘴出流）管嘴出流将上例孔口改为一个短管。按照同样的方法可求得：但此时管嘴的流量系数为：第五十二页，本课件共有61页（3）串联管路的水力计算）串联管路的水力计算串联管路是各种不同直径的管路串联而成。如下图所示：由液面0-0及出口截面3-3建立伯努利方程，有：第五十三页，本课件共有61页 故：对于无外泄漏的串联管路，各段流量应相等，即：当有外泄漏时，各段流量应为：第五十四页，本课件共有61页对于长管（管路系统中局部损失相对沿程损失来说，可以忽略）来说，管路系统中的总水头损失等于各段沿程损失之和。即：假如引入流量模数流量模数K的概念（K具有流量的单位）：并且假定各管段无外泄漏，则总水头损失还可表达为：显然，此时流量Q：第五十五页，本课件共有61页（4）并联管路水力计算）并联管路水力计算并联管路：并联管路：由若干个简单长管并联而成。并联管路有下列特点：A：并联管路各支路上的阻力损失相同（因为在共同的分流处及汇流处，对于各支路来说的有共同的压强）。即：因此，对整个管路的损失，有：B：总流量等于各支路流量之和。即：第五十六页，本课件共有61页（5）环状管网的水力计算）环状管网的水力计算一般用逐次迭代法，计算量很大，一般要计算机编程计算。计算时一般可利用下列两个原则：A：节点连续性条件：即每个节点上流入 的流量等于流出的流量。假如流入为 正，流出为负，则总和应为0。即：B：沿环状管路，每条支路的流动方向不一定相同，我 们沿环路可作一逆时针曲线，当管路流向与逆时针曲线方向相同时，水头损失为正，相反时为负。沿整个环路，总的水头损失之和应等于0。即：迭代计算时，首先假定每一个管道的流量及流向，然后计算水头损失。迭代计算时，首先假定每一个管道的流量及流向，然后计算水头损失。看看是是否否满满足足 ，如如不不满满足足，则则调调节节流流量量，继继续续计计算算，直直到到满满足足为为止止。一一般般要要编编程程计算。计算。第五十七页，本课件共有61页4.6 本章基本理论应用举例本章基本理论应用举例例例一一：如图，用虹吸管从水箱中取水，虹吸管管径d=5cm，管总长10.5m，管道沿程损失系数=0.03，弯管局部损失系数为0.2，LAB=3.5M，试求管中流量和管内最大真空度。解：以液面及管出口处建立伯努利方程：可求得：第五十八页，本课件共有61页虹吸管最高点B处真空度为最大，故由液面及管道B点处建立伯努利方程：求得：真空度为：第五十九页，本课件共有61页例二例二：如图所示，串联管路，已知：计及截面突变时的局部损失，试求流量Q。解：由液面及管出口截面处建立 伯努利方程：第六十页，本课件共有61页由查表可得：代入已知参数可得：上面方程有三个未知参数，再用连续性方程：从而可求得：第六十一页，本课件共有61页