水力学教学课件 第五章 实际流体动力学基础.ppt

'第五章实际流体动力学基础§5—0实际流体动力学§5—1实际流体的运动微分方程—纳维-斯托克斯方程§5—2实际流体元流的伯努利方程§5—3实际流体总流的伯努利方程第五章实际流体动力学基础§5—5总流的动量方程 5-0实际流体动力学§5—0实际流体动力学1、实际流体动力学（fluiddynamics）的力学性质：实际流体具有粘性。在作用面上的表面力不仅有压应力即动压强，还有切应力。2、作用在一平面上M点的表面应力表面应力pn在x、y、z三个轴向都有三个分量：与平面成法向的压应力pzz，即动压强；与平面成切向的切应力τzx，和τzy。3、通过任一点在三个互相垂直的作用面上的表面应力。共有九个分量三个是压应力px、py、pz。六个是切应力τxy、τxz、τyx、τyz、τzx、τzy。压应力和切应力的第一个下标，表示作用面的法线方向，即表示应力的作用面与那一个轴相垂直；另两个下标，表示应力的作用方向，即表示应力作用方向与那一个轴相平行。 5-1实际流体的运动微分方程—纳维-斯托克斯方程5-1实际流体的运动微分方程—纳维-斯托克斯方程5-1-1以应力表示的实际流体的运动微分方程 5-1实际流体的运动微分方程—纳维-斯托克斯方程5-1实际流体的运动微分方程—纳维-斯托克斯方程5-1-1以应力表示的实际流体的运动微分方程设M点的坐标为(x,y,z)，速度、压应力、切应力、单位质量力分别为ux,uy,uz,px,py,pz,τxy,τxz,τyz,τzx,τzy,τzx,fx,fy,fz。压力为：质量力为：x轴方向受力分析：惯性力为：左右面、上下面的切力为：dzdxdy"yz"yxp"yyxzxypxxzxzypzz"xy"xzp"xxyzyxpyy"zy’zxp"zzxyz 5-1实际流体的运动微分方程—纳维-斯托克斯方程化简移项后得和上式即为以应力形式表示的实际流体的运动微分方程。同理：根据牛顿第二定律(5-1)对于不可压缩均质流体来讲，密度ρ为常数，单位质量力的分量fx、fy、fz通常是已知的，所以上式中有表面应力的九个分量和速度的三个分量，共十二个未知量。式(5-1)只有三个方程式，加上连续性微分方程也只有四个方程式，所以无法求解，需找出其他的关系式。这些其他的关系式需从分析流体质点的应力状态中获得。dzdxdy"yz"yxp"yyxzxypxxzxzypzz"xy"xzp"xxyzyxpyy"zy’zxp"zzxyz 5-1实际流体的运动微分方程—纳维-斯托克斯方程1、切应力的特性和大小5-1-2流体质点的应力状态（1）切应力的互等定律-------(5-2)在九个表面应力分量中，实际上只有六个是独立的，即为px,py,pz,τxy,τyz,τzx（2）切应力和它的应变(变形)之间的关系——广义牛顿内摩擦定律------(5-3) 5-1实际流体的运动微分方程—纳维-斯托克斯方程2、压应力的特性和大小（1）同一点上三个正交方向的压应力的平均值p是单值，它与方位无关。在实际问题中，某点压应力的各向差异并不很大，平均值p作为该点的压应力是允许的。（2）各个方向的压应力可认为等于这个平均值加上一个附加压应力：------(5-4)负号是与流体微团的拉伸和压缩相适应，即周围流体对它作用的是拉力和压力。------(5-5)则各个方向的压应力 5-1实际流体的运动微分方程—纳维-斯托克斯方程3、实际流体中任一点的应力状态讨论（1）理想流体，μ=0,px=py=pz=p。（2）实际流体的均匀直线流中ux=a(不为零的常数)，uy=O，uz=0和固体边壁处ux=uy=uz=0。（3）在实际流体中任一点的应力状态就可由一个压应力(即动压强)p和三个切应力τxy,τyz,τzx来表示。px=py=pz=p。 5-1实际流体的运动微分方程—纳维-斯托克斯方程1、纳维-斯托克斯方程5-1-3实际流体的运动微分方程—纳维-斯托克斯方程N-S方程中有四个未知数p、ux、uy、uz，因N-S方程组和连续性方程共有四个方程式，所以从理论上讲是可求解的，但实际上由于数学上的困难，N-S方程尚不能求出普遍解。上式即为不可压缩均质实际流体的运动微分方程，称纳维-斯托克斯(Navier-stokes)方程，简称N-S方程。N-S方程是不可压缩均质流体的普遍方程。----(5-6)拉普拉斯算符：考虑：不可压缩流体的连续性微分方程，切应力与主应力的关系表达式，得： 5-1实际流体的运动微分方程—纳维-斯托克斯方程式中：分别为单位质量力在坐标轴上的分量。2、在柱坐标系中，纳维-斯托克斯方程----------(5-7) 5-1实际流体的运动微分方程—纳维-斯托克斯方程3、纳维-斯托克斯方程求解条件初始条件：在起始时刻t=0时，各处的流速、压力值；对于恒定流，则不存在条件。边界条件：一般包括固体边界和自由表面等处的运动要素情况。固体边界：，不允许流体在边界上有滑移。自由表面：例5-1：圆管中的层流运动(哈根-泊肃叶流动）的流量与压降的关系（哈根-泊肃叶定律，或简称泊肃叶定律）自习 5-1实际流体的运动微分方程—纳维-斯托克斯方程例5-1：圆管中的层流运动(哈根-泊肃叶流动）的流量与压降的关系（哈根-泊肃叶定律，或简称泊肃叶定律）设实际流体在很长的水平圆管内作有压恒定(层流)运动，如图5-3所示。已知管径为dl速度ux=u(y，z)，uy=0，uz=0。试用N-S方程求解过流断面上的速度分布和通过的流量以及断面平均通度。解；由连续性方程和N-S方程可得微分方程：边界条件：(0)(1)(2)(3)固体边界无滑条件(5)(4)管轴速度梯度为零条件 5-1实际流体的运动微分方程—纳维-斯托克斯方程由(2)、(3)式知，测压管水头沿圆管断面为常数，只与x有关。因此(1)式可写成为式等号左边起y，z的函数，右边是x的函数，因此两边相等的条件是两边均为常数，即为一常数，这样，就可写成：在yoz平面内引进极坐标(r，θ)。由于管流的对称性(4)速度只是r的函数，即u(y，z)＝u(r)。或 5-1实际流体的运动微分方程—纳维-斯托克斯方程积分，得：利用边界条件（5），得C1=0积分，得利用边界条件（4），得得转化为代数方程上式表明，圆管中有压恒定层流过流断面上的速度分布是按抛物面的规律分布。 5-1实际流体的运动微分方程—纳维-斯托克斯方程流量：断面平均速度 5-1-4实际流体运动微分方程的积分5-1实际流体的运动微分方程—纳维-斯托克斯方程满足下列条件下：（1）流体是不可压缩均质的理想流体，密度为常数；（2）质量力是有势的；（3）恒定流；（4）沿流线运动积分；运动微分方程简化为：-------(5-10)第二项表示单位质量流体沿流线作微小位移时应力(主要是切应力)所作的微功。流体在作这些功的所消耗的机械能，就是流体的能量损失。（5）作用于流体上的质量力只有重力；运动微分方程简化为：-------(5-11)第二项表示对单位重量而言的应力所作的微功。以表示。 对上式沿流线由点1到点2积分，得5-1实际流体的运动微分方程—纳维-斯托克斯方程式中。----------(5-13)上式为不可压缩均质实际流体恒定流的伯努利方程式，又称伯努利方程。限于同一条流线上各点的总机械能保持不变；这和理想流体伯努利方程在有势流的应用条件不同。 5-2实际流体元流的伯努利方程5-2实际流体元流的伯努利方程5-2-1实际流体元流的伯努利方程1、用元流分析法推导出不可压缩均质实际流体恒定元流的伯努利方程：元流段的动能增量:重力所作的功为：根据动能定理压力所作的功为：得：----------(5-14)阻力所作的功为： 5-2-2实际流体元流伯努利方程的物理意义和几何意义5-2实际流体元流的伯努利方程1、实际流体元流伯努利方程的物理意义z是元流过流断面上单位重量流体从某一基准面算起所具有的位能，称单位位能。p/ρg是元流过流断面上单位重量流体从某一基准面算起所具有的位能压能，称单位压能。z+p/ρg是元流过流断面上单位重量流体从某一基准面算起所具有势能，称单位势能。物理意义：元流各过流断面上单位重量流体所具有的总机械能沿流程减小，部分机械能转化为热能或声能等而损失；同时，亦表示了各项能量之间可以相互转化的关系。u2/2g是元流过流断面上单位重量流体所具有的动能，称单位动能。h＇w是单位重量流体由过流断面1-1移动到过流断面2-2时的能量损失。是单位重量流体总机械能。 2、实际流体元流伯努利方程的几何意义5-2实际流体元流的伯努利方程z是位置水头。p/ρg是压强水头。z+p/ρg是测压管水头。几何意义：对于液体来说，元流各过流断面上总水头H(位置水头、压强水头、速度水头之和)沿流程保持减小；同时，亦表示了元流在不同过流断面上各水头水头之间可以相互转化的关系。u2/2g是速度水头。h＇w是是能量损失，称水头损失或损失水头(headloss)。是总水头（totalhead）。 水头线和水力坡度5-2实际流体元流的伯努利方程水头线：沿程水头（如总水头或测压管水头）的变化曲线。测压管水头线(piezometricheadline)总水头线(totalheadline) 水头线和水力坡度5-2实际流体元流的伯努利方程单位长度上的水头损失称能线坡度或水力坡度(hydraulicslop)。Jp可正、可负或为零。将顺流程向下的J视为正。J总是正值。单位长度上测压管水头的降低或升高，称测压管水头线坡度Jp。----------(5-15)----------(5-16)水头线：沿程水头（如总水头或测压管水头）的变化曲线。若是均匀流，则总水头线平行于测压管水头线，即J=JP。 不同固体边界下的水头线5-2实际流体元流的伯努利方程转弯转弯突扩突缩闸门V测压管水头线P-P总水头线E-EhjhfHOO出口注：出口为自由出流时，P-P线末端应落在出口断面的管轴线上。 5-3实际流体总流的伯努利方程5-3实际流体总流的伯努利方程5-3-1实际流体总流的伯努利方程总流是由无数微小流束组成的。因此，只要对式(5-14)进行积分，便可得到实际液体总流的能量方程。（1）测压管水头积分：总流的过流断面取在渐变流区域，那么，在一定的边界条件下，同一过流断面上的压强分布近似地按静压强规律分布，z+p/ρg可视为常数。----------(5-20)（2）速度水头积分：总流的过流断面取在渐变流区域，所以同一过流断面上各点速度和断面平均速度的方向可认为相同，且垂直于过流断面，但是大小不同。α——动能修正系数(kinetic-energycorrectionfactor)，是同一过流断面上各点速度不等时的实际动能与假设该过流断面上各点速度均为断面平均速度时的动能比值。α＞1.0；一般取α＝1.05-1.10之间，通常在近似计算时取α＝1.0。----------(5-23) （3）水头损失积分。令hw为单位重量流体由过流断面1-1移动到过流断面2-2能量损失的平均值5-3实际流体总流的伯努利方程整理得：上式即为实际流体总流的伯努利方程(能量方程)。方程式(5-26)各项的物理意义和几何意义类似于实际流体元流的伯努利方程式中的对应项，所不同的是指平均值。物理意义是:总流各过流断面上单位重量流体所具有的势能平均值与动能平均值之和，亦即总机械能的平均值沿流程减小，部分机械能转化为热能等而损失；同时，亦表示了各项能量之间可以相互转化的关系。几何意义是:对于液体来讲，总流各过流断面上平均总水头沿流程减小，所减小的高度即为两过流断面间的平均水头损失；同时，亦表示了各项水头之间可以相互转化的关系。。 动能及动能修正系数讨论5-3实际流体总流的伯努利方程动能（KineticEnergy）：是指物体由于机械运动而具有的能量。——单位重量流体的平均动能。——动能修正系数（层流=2.0，紊流=1.05~1.1），是实际动能与按断面平均流速计算的动能的比值，即动能修正系数是一无量纲数，取决于总流过水断面的流速分布，分布越均匀，值越小，越接近于1.0。层流流速分布紊流流速分布 5-3-2总流伯努利方程的应用条件和应用方法5-3实际流体总流的伯努利方程1、应用总流伯努利方程必须满足下列条件(1)流体运动是恒定流；(2)流体运动符合连续原理；(3)作用于流体上的质量力只有重力；(4)所取过流断面在渐变流或均匀流区域，但两断面间不必是渐变流或均匀流；(5)所取两过流断面间没有流量汇入或流量分出，亦没有能量输人或输出。否则，总流能量方程采用有分流和没有能量输人或输出形式仍适用。2、应用总流伯努利方程的步骤、方法(1)分析流动现象。弄清楚流体运动的类型，建立流体运动的流线几何图形(流谱)，判断是否能应用总流的伯努利方程；(2)选取好过流断面。所取断面须在渐变流或均匀流区域，且要垂直于流线；(3)选择好计算点和基准面。有压管流时取在管轴线上，明渠流时取在自由表面上(因压强为大气压强，可作为已知值)；基准面的选择，一般使z值为正值。(4)压强的表示方法，一般是以相对压强计，亦可用绝对压强计，但在同一方程中必须一致；所取单位要一致。(5)全面分析和考虑所取两过流断面间的能量损失。 3、流动阻力和能量损失5-3实际流体总流的伯努利方程一般将流动阻力(drag/resistance)和由于克服阻力而消耗的能量损失，按决定其分布性质的边界几何条件而分为两类。（1）沿程阻力、沿程损失式中:ξf称沿程损失系数，一般由实验测定。均匀分布在某一流段全部流阻的流动阻力称沿程阻力；克服沿程阻力而消耗的能量损失称沿程损失。单位重量流体沿程损失的平均值以表示。----------(5-29) （2）局部阻力、局部损失5-3实际流体总流的伯努利方程几何条件的急剧改变而引起对流体运动的阻力称局部阻力；克服局部阻力而消耗的能量损失称局部损失。式中:ξf称局部损失系数，一般由理论推导、实验测定。单位重量流体局部阻力损失的平均值以表示。----------(5-30)能量损失的叠加原理：任何两过流断面间的能量损失hw，在假设各损失是单独发生，且又互不干扰、影响的情况下，可视为是每个个别能量损失的简单总和。----------(5-27) 4、绘制水头线5-3实际流体总流的伯努利方程（1）一般先绘总水头线，因为在没有能量输人的情况下，它一定是沿流程下降的；然后绘测压管水头线。（4）局部损失实际上是在一定长度内发生的，但常集中地画在突变断面上。（2）已知的过流断面上的总水头端点和测压管水头端点可作为水头线的控制点(始点和终点)。（3）沿程损失则认为是均匀分布的，常画在两边界突变断面间。 预测管线真空范围；确定管线最高布置位置。5-3实际流体总流的伯努利方程∵在BC段，由图中测压管水头线知：v2=0hj出v1=0vhj进真空段BC即为真空。真空区zz+p/-p/OO虹吸管的水头线与真空区 5-3-3文丘里(Venturi)管5-3实际流体总流的伯努利方程测量恒定有压管流的文丘里管，它由渐缩段、喉道和渐扩段三部分所组成。渐缩段前的断面1-1及断面最小的喉道断面2-2处布置测压孔，并接上测压装置(压差计)。任选O-O为基准面，列1、2两断面能量方程，则得忽略水头损失，令动能修正系数＝1.0，得因为，整理得： 通过文德里流量计中的流量为5-3实际流体总流的伯努利方程当计入损失后，流量将减少，因而须对式乘以修正系数加以修正令:式中，K对给定管径是常量，称为文丘里流量计常数。于是得式中:μ称流量系数，它不是一常数，随水流情况和文丘里管的材料性质、尺寸等而变化，须对各文丘里管专门率定才能确定，一般μ≈0.95—0.98。 如遇较大压强差，可在文丘里流量计上安装水银比压计5-3实际流体总流的伯努利方程式中:γm-水银容重，γ-水容重。于是得式中:γm-水银容重，γ-水容重。 5-3-4汇流或分流的伯努利方程•有能量输入或输出的伯努利方程5-3实际流体总流的伯努利方程对汇流，每支总流建立伯努利方程为对分流，每支总流建立伯努利方程为1、汇流或分流的伯努利方程----------(5-32)----------(5-31) 2、有能量输入或输出的伯努利方程5-3实际流体总流的伯努利方程建立伯努利方程为有能量输入：Hm为单位重量的水流通过水泵后增加的能量，称管路所需的水泵扬程；Hw为全部管路中的水头损失，不包括泵内的损失。-----(5-34、37)--------(5-33、36)因为pl=p2=pa；v1,v2相对于管内流速来讲均较小，简化方程式中：z为上、下游水面高差，也称提水高度或扬水高度。有能量输出：Hm为单位重量的水流给予水轮机的能量，称水轮机的作用水头；Hw为全部管路中的水头损失，不包括水轮机内的损失。 2、有能量输入或输出的伯努利方程5-3实际流体总流的伯努利方程单位时间内原动机给予水泵的功称水泵的轴功率P。有效功率与轴功率的比值称水泵效率η,以%计：单位时间内水流从泵中实际获得的总能量为ρɡQHM，称水泵的有效功率P。-----(5-35) 例5-2设有一水位保持不变的很大的水箱，在其侧壁开小圆孔口，水从小孔口流人大气，如图5-4所示。己知水箱内自由表面到小孔口中心的水深为H，小孔口直径d0，pc值小于大气压强pa值，流体中出现真空。液体中某处的绝对压强，低于当地的汽化压强大时，形成空化现象，所以，h2值不可能无限制地增大。 （3）从虹吸管中水流流动的情况看，虹吸管顶部水平管中心点D的压强比管中点C的压强还要小，所以要限制h1值。现对A、D两点列伯努利方程：5-3实际流体总流的伯努利方程得由于vb=vD，得为了不发生汽化，h1最大的可能值是在h2=O，pD=pv的情况，由表1-1查得水温40°C时pv=7.38kPa；代人上式 5-5总流的动量方程5-5总流的动量方程5-5-1总流的动量方程1、用有限体分析法(元流分析法)，根据动量(冲量)定理来推导出不可压缩均质实际流体恒定元流的动量方程：元流段的动量增量即为1-1"段和2-2"段流体动量之差:根据动量(冲量(impulse))定律，物体在运动过程中，动量对时间的变化率，等于作用在物体上各外力的合力矢量，即：上式即为不可压缩均质实际流体恒定元流的动量方程。得：式中:Fe—作用于运动液体上的外力的合力。dp—在dt内的动量变化量。 5-5总流的动量方程积分项：对动量项积分：式中：β称动量修正系数（momentumcorrectionfactor），它可理解为是同一过流断面上各点速度不等时的实际动量与假设该过流断面上各点速度均为断面平均速度时的动量的比值。其中：对作用力项积分：式中：∑F—外力的合力。作用于控制体内流体的所有外力矢量和。该外力包括作用在该控制体内所有流体质点的质量力、作用在该控制体面上的所有表面力、以及四周边界对水流的总作用力的合力之矢量和。-----(5-42)2、总流可以看成是由流动边界内无数元流所组成，对总流过流断面面积积分，得总流的动量关系——总流的动量方程。 导出总流的动量关系:5-5总流的动量方程上式表明，单位时间内流出控制面(过流断面2-2)和流人控制面(过流断面1-1)的动量矢量差，等于作用在所取控制体内流体(总流段)上的各外力的合力矢量。上式即为不可压缩均质实际流体恒定总流的动量方程。写成分量形式：-----(5-43)-----(5-44) 动量及动量修正系数讨论5-5总流的动量方程动量（Momentum）：是物体运动的一种量度，是描述物质机械运动状态的物理量。——总流通过整个过水断面其相应的动量值为（在过水断面n方向上）。——动量修正系数（层流=4/3，紊流=1.02~1.05)指实际动量与按断面平均流速计算的动量的比值。动能修正系数是一无量纲数，取决于总流过水断面的流速分布，分布越均匀，值越小，越接近于1.0。层流流速分布紊流流速分布 5一5一2总流动量方程的应用条件和应用方法5-5总流的动量方程应用总流动量方程的步骤和方法(1)分析流动现象。首先要弄清楚流体运动的类型，建立流线几何图形(流谱)，判断是否能应用总流的动量方程。(2)选取好控制面和控制体。为计算方便，过流断面取在渐变流区域。(3)全面分析和考虑作用在控制体内流体系统(隔离体)上的外力。既要做到所有外力一个不漏，又要考虑哪些外力可以忽略不计。若所选封闭控制面中，有为流体(水)所环绕的不动固体时，固体表面亦是控制面的组成部分，固体对于环绕它的流体的反作用力，亦须包括在作用于流体的表面力之内。(4)压强的表示方法。由于在一般情况下，控制面常处于大气压强作用下，计算作用于控制面的压强时，宜采用相对压强。(5)方程式中的动量矢量差是指流出的动量减去流人的动量，两者切不颠倒。(6)正确取好外力和流速的正负号。对于己知的外力和流速的方向，凡是与选定的坐标轴方向相同者取正号；与坐标轴方向相反者取负号。对于末知待求的，则可先假定为某一方向，并按上述原则取好正负号，代入总流动量方程中的有关项。如果最后求得的结果为正值，说明假定的方向即为实际的方向；如果为负值，则说明假定的方向为实际的相反方向。 5-5总流的动量方程思考题1、实际流体区别于理想流体有何特点？理想流体的运动微分方程与实际流体的运动微分方程有何联系？1、为什么要引入动能修正系数和动量修正系数？2、拿两张薄纸，平行提在手中，当用嘴顺纸间缝隙吹气时，问薄纸是不动、靠拢、还是张开？为什么？3、恒定总流能量的限制条件有哪些？任何选取其计算断面、基准面、计算点、压强？4、什么是水头线？总水头线与测压管水头线有何区别？5、在应用恒定总流动量方程时，为什么不必考虑水头损失？为内力6、由动量方程求得的力若为负值时说明什么问题？待求未知力的大小与脱离体的大小有无关系？应用中如何选取脱离体？方向反；无关（无重力时）；计算断面与固体壁面 5-5总流的动量方程思考题1、总流能量与元流能量方程有什么不同点？（1）以断面的平均流速V代替元流中的点流速u；（2）以平均水头损失h代替元流的水头损失h´1.2；（3）各项反映的是整股水流的能量代替某一元流的能量。2、设有一水平压力管流，当不考虑水头损失的影响时，其测压管水头线沿程下降、上升或水平的条件各是怎样的？下降：压强沿程减小。例：收缩管水平：压强沿程不变。例：管径沿程不变上升：压强沿程增大。例：扩散管3、水流由等径斜管的高处A流向低处B（不考虑水头损失）。讨论压强出现以下三种情况时的情形：pA>pB；pA=pB；pApB：AB流动pA=pB：AB流动pA(z1-z2)，BA流动 本章小结本章小结一、总流分析法的三个处理方法二、能量方程：1、与元流能量方程的不同点；2、适用条件；3、总水头线、测压管水头线的绘制；4、应用。三、动量方程：1、适用条件；2、应用。'