工程流体力学第三版课件B.ppt

《工程流体力学第三版课件B.ppt》由会员分享，可在线阅读，更多相关《工程流体力学第三版课件B.ppt（358页珍藏版）》请在三一办公上搜索。

1、相似原理的提出,流体力学理论的检验依赖于流体力学试验；,流体力学的模型实验：工程实际需要的流体力学试验一般很难在实物上进行。,飞机风洞试验,汽车风洞试验,水利大坝试验,轮船水洞试验,飞机风洞试验,返回,汽车风洞试验,返回,水力大坝试验,返回,轮船水洞试验,返回,怎么做模型试验？,1.如何根据实物正确的设计和布置模型实验？,2.模型实验的结果如何推广到原型上去，并进行推广应用？,相似原理,第一节 流动的力学相似,一个流动某点的运动参数由另一个流动相应点的同名参数乘以对应点均相同的因子得到，称两流动相似。,一、流动相似的概念,具体的说两流动相似应满足几何相似、运动相似和动力相似三个条件。,B,1.

2、几何相似,几何相似是指模型与原型的全部对应线性长度成比例，且对应的特征角度相等。,长度比例系数：,面积比例系数：,体积比例系数：,2.运动相似（时间相似）,运动相似是指：模型与原型的流场中所有对应点上对应时刻的流速方向相同，且对应流速的大小的比例相等，即它们速度场相似。,时间比例系数：,加速度比例系数：,速度比例系数：,由上述比例系数可推出：,3.动力相似,动力相似是指模型与原型的流场所有对应点上作用在流体微团上的各种力彼此方向相同，且它们大小的比例相等，即它们的动力场相似。,力比例系数：,也可写成：,综上所述：,在做模型试验时，要想使两个流动相似必须在几何相似、运动相似和动力相似三个方面都得

3、到满足。,实际应用中，并不能用定义来检验流动是否相似，因为通常原型的流动是未知的。,相似准则,第二节 动力相似准则,什么是相似准则？,两个矩形要相似必满足：,相似准则数,流体的运动必须符合牛顿第二律，对模型和原型流场中的流体微团应用牛顿第二定律，并根据动力相似，各种力大小的比例相等，可得：,或,动力相似准则的推导,令：,Ne称为牛顿数，它是作用力与惯性力的比值。,Ne称为牛顿数，它是某种作用力与惯性力的比值，是无量纲数。由此可知，模型与原型的流场动力相似，它们的牛顿数必相等。,作用在流体微团上的作用力有各种性质的力，如重力、粘滞力、压力等，根据上式可导出单项力相似准则。,在粘滞力作用下相似的流

4、动，其粘滞力场相似。,Re雷诺数，惯性力与粘滞力的比值。,一、粘滞力相似准则,代入,在重力作用下相似的流动，其重力场相似。,二、重力相似准则,Fr-弗劳得数，惯性力与重力的比值。,代入,在压力作用下相似的流动，其压力场相似。,代入,Eu-欧拉数，压力与惯性力的比值。,压力相似准则,三、其它的相似准则数,弹性力相似准则 对于可压缩流体的模型试验，由压缩引起的弹性力场相似。（Ca柯西数 Ma马赫数，惯性力与弹性力的比值）。,非定常相似准则 对于非定常流动的模型试验，模型与原型的流动随时间的变化必相似。（Sr 斯特劳哈尔数，当地惯性力与迁移惯性力的比值）。,3.应用举例,1）如果模型比例尺为1:20

5、，考虑粘滞力相似，采用模型中流体与原型中相同，模型中流速为50m/s，则原型中流速为多少？,查看答案,解：由粘滞力相似准则知模型与原型中的雷诺数应相等：,所以：,雷诺数：,因为：,返回,2）设模型比例尺为1:100，符合重力相似准则，如果模型流量为100 cm3/s，则原型流量为多少 cm3/s?,A：0.01 B：1000 C：10 D：10000,答案:c,1.相似理论的提出：相似理论的作用是什么？,2.流动相似的概念：两流动相似的条件是什么？,3.相似准则：为什么应用相似准则？怎么用？,小结,第三节 流动相似条件,流动相似：在对应点上、对应瞬时，所有物理量 都成比例。,相似流动必然满足以

6、下条件：,1任何相似的流动都是属于同一类的流动，相似流场对应 点上的各种物理量，都应为相同的微分方程所描述；2相似流场对应点上的各种物理量都有唯一确定的解，即 流动满足单值条件；3由单值条件中的物理量所确定的相似准则数相等是流动 相似也必须满足的条件。,模型实验主要解决的问题：,1根据物理量所组成的相似准则数相等的原则去设计模 型，选择流动介质；2在实验过程中应测定各相似准则数中包含的一切物理量；3用数学方法找出相似准则数之间的函数关系，即准则方程 式。该方程式便可推广应用到原型及其他相似流动中去。,油池模型,第四节 近似模拟试验,完全相似和不完全相似 动力相似可以用相似准则数表示，若原型和模

7、型流动动力相似，各同名相似准数应均相等，如果满足则称为完全的动力相似。但是事实上，不是所有的相似准数之间都是相容的，满足了甲，不一定就能满足乙。所以通常考虑主要因素忽略次要因素，只能做近似的模型实验。,例如：粘滞力相似：由 得重力相似：由 得,由此可以看出，有时要想做到完全相似是不可能的，只能考虑主要因素做近似模型实验。,以相似原理为基础的模型实验方法，按照流体流动相似的条件，可设计模型和安排试验。这些条件是几何相似、运动相似和动力相似。前两个相似是第三个相似的充要条件，同时满足以上条件为流动相似，模型试验的结果方可用到原型设备中去。,在工程实际中的模型试验，好多只能满足部分相似准则，即称之为

8、局部相似。如上面的粘性不可压定常流动的问题，不考虑自由面的作用及重力的作用，只考虑粘性的影响，则定性准则只考虑雷诺数Re，因而模型尺寸和介质的选择就自由了。,有压粘性管流中，当雷诺数大到一定数值时，继续提高雷诺数，管内流体的紊乱程度及速度剖面几乎不再变化，沿程能量损失系数也不再变化，雷诺准则已失去判别相似的作用。称这种状态为自模化状态，称自模化状态的雷诺数范围为自模化区。,一、物理方程量纲一致性原则,第五节 量纲分析,1、量纲,量纲是物理量的一种本质属性，是同一物理量各种不同单位的集中抽象。,如：s单位：km，m，cm，mm 等 t单位：hour，min，second 等s-具有长度的量纲L

9、t-具有时间的量纲T V-具有速度的量纲,同时还有,如质量量纲M，力的量纲F等。,基本量纲-相互独立，不相互依赖，如M，L，T等。导出量纲-由基本量纲导出,如,一个合理的物理方程等号两端的量纲必须相同。,2、方程量纲一致性,-方程两端具有相同量纲,量纲式中各基本量纲指数均为零-无量纲量。,二、瑞利法,1.定义：根据量纲量一致性原则，确定相关量的函数关系。,假定物理量y是x1、x2等的函数。则,关键的问题是怎么根据量纲一致性原则确定各个x的指数。,2.举例：,三角堰,三、定理,由于方程量纲一致性，用其中任意一项去除方程两端，都可将有量纲函数式变成无量纲函数式。,但这些无量纲函数式可以互相推导：,

10、它们所包含的独立无量纲量数目（不能相互推导）是多少？,共有两个，它们可以表示为：,为什么？,与匀变速运动有关物理量的量纲指数：,设四个物理量组成的函数形式：,如果是无量纲量其基本量纲指数为零。,-只有那些线性无关的解组（基础解系）才能组成彼此独立的无量纲量。,当,当,解向量,解向量,独立无量纲量(m)=与问题有关的物理量(n)-量纲指数矩阵秩(r),柏金汉（Buckingham）定理：如果与某一物理问题有关的物理量有n个，这n个量的量纲指数矩阵的秩为r，则表达这一问题规律性的n个量之间的函数关系，可以用这n个量组成的n-r个独立无量纲量的函数式来表示：,1.选取影响流动的 n 个物理量并写出下

11、述函数关系如,求解过程：,2.选择 m 个独立变量，原则是要既相互独立，又包含三个基本量纲.一般选:,几何尺度,几何尺度,速度,质量,3.用 n m 个无量纲量写出准则方程,4.求,5.将 带入准则方程式求得 结果。,例.流体通过孔板流量计的流量qv与孔板前、后的压差P、管道的内径d1、管内流速v、孔板的孔径d、流体密度和动力粘度有关。试用定理导出流量qv的表达式。,（dimP=ML-1T-2，dim=ML-1T-1）。,查看答案,解：设qv=f(P,d1,v,d,),选d,v,为基本变量,上述方程的量纲方程为：,由量纲一致性原则，可求得：,a1=0 a2=1 a3=0 a4=1b1=1 b2

12、=2 b3=0 b4=1c1=2 c2=0 c3=1 c4=1,第六章,管内流动和水力计算 液体出流,第一节管内流动的能量损失,一、沿程损失,-沿流程上流体与壁面以及流体本身内部摩擦而产生的能量损失（用hf来表示）。,沿程损失，是发生在缓变流整个流程中的能量损失，是由流体的粘滞力造成的损失。,L：管长，d：管径，V：管断面平均速度，：沿程阻力系数。,影响因素,二、局部损失,-流动中，由于边界急剧变化（如管径突然变大或变小；弯管引起流速方向改变；或阀门、三通等）而产生的局部能量损失(一般用hj表示)。,局部损失：是发生在流动状态急剧变化的急变流中的能量损失。是主要由流体微团的碰撞、流体中的涡流等

13、造成的损失。,变径管,发生位置,弯头,阀门,渐缩,渐扩,突缩,突扩,计算公式：,V：断面平均速度，：局部阻力系数。,若为管路系统，能量损失应是各段沿程损失和局部损失之和，即,局部阻力系数由试验确定。,第二节黏性流体的两种运动状态,一、雷诺实验 两种流态,流体分层运动，各层间互不干扰、互不相混的流动状态。,1.层流,流体质点运动彼此混杂、互相干扰，完全无规则的流动状态。,2.紊流,3.上临界速度和下临界速度：,随着水流速度的增大，水流将由层流状态过渡到紊流状态。由层流过渡到紊流的临界状态下的流体速度称为上临界速度，用Vcr表示。,当玻璃管内的水流已经是紊流运动，此时逐渐关小阀门K，使水流速度逐渐

14、减小，当水流速度减小到一定程度时，紊乱的红色液体又将重新成为一条明晰的红色直线流，即紊流又转变为层流。但是，由紊流转变为层流的临界速度比上临界速Vcr更低，称为下临界速度，用Vcr表示。,实验表明，这两种情况下的流动状态都不稳定，并且取决于实验的起始状态有无扰动等因素。,说明,(1)当流体的流速超过上临界速度(VVcr)，管内水流一定是紊流状态；,(2)当流体的流速低于下临界速度时(VVcr)，管内水流一定是层流状态；,(3)当流体的流速介于上临界速度和下临界速度之间时(VcrVVcr)，管内水流可能是层流，也可能是紊流。如果流速是由小增大时，流动是层流，如果流速是由大变小时，则流动是紊流。,

15、二、沿程损失和平均流速的关系,层流状态,紊流状态,m=1,m=1.752,可能是层流，也可能是紊流,沿程损失和平均流速的关系图,三、流态的判别临界雷诺数,实验发现，判别流体的流动状态，仅靠临界速度很不方便，因为随着流体的粘度、密度以及流道线尺寸的不同，临界速度在变化，很难确定。雷诺根据大量的实验归纳出一个无因次综合量作为判别流体流动状态的准则，称为雷诺准则或雷诺准数，简称雷诺数，用Re表示，即,式中V为流体的特征流速，d为流体通道的特征尺寸。对于直径为d的圆截面管道，有,对应于临界速度的雷诺数称为临界雷诺数，用Recr表示，,流体的流动状态是层流还是紊流，对于流场的速度分布、产生阻力的方式和大

16、小，以及对传热传质过程和动量传递规律等都各不相同，所以在研究这些问题之前，首先需要判别流体的流动是属于哪一种状态。,实验结果表明，对于光滑的圆截面直管，不论流体的性质和管径如何变化，其下临界雷诺数一般均为Rec=21002300，而上临界雷诺数Recr可达1200013800，甚至更高些，但这时流动处在极不稳定的状态，稍有扰动层流瞬即被破坏而转变为紊流。因此，上临界雷诺数在工程上没有实用意义，通常用下临界雷诺数来判别流体的流动状态，即取圆管内流动的临界雷诺数为Rec=2300。对于圆截面管道，当Re2300时为层流，Re2300时为紊流。,说明,三、流态分析,雷诺数之所以能判别流态，是因为它反

17、映了流体运动时惯性力与粘滞力的对比关系：,当ReRecr，惯性力起主导作用，粘性力控制减弱，不足以控制和约束外界扰动，惯性力将微小扰动不断扩大，形成紊流。,当Re较小时，粘性力作用大，对质点运动起约束作用，流体质点表现为有秩序互不掺混的层流状态；,第三节管道进口段中粘性流体的流动,一、圆管内层流流动的起始段,由于流体的粘性作用，自圆管入口起，在管壁附近形成一层有速度梯度存在的流体薄层，该流体薄层内壁面上流体的速度为零，薄层外边界上的流速为u(x)。这一有速度梯度存在的流体层称为附面层或边界层。,从管进口到附面层在管中心汇合处的截面间的一段距离L*称为层流(紊流）的起始段。以下将证明，在起始段以

18、后的各管截面上的速度分布均为抛物线分布(对数曲线)。起始段以后的管段称为层流（紊流）的充分发展段。,层流:,兰哈尔,L*0.058dRe,L*经验公式,紊流:,L*（2540）d,L*(层流)L*（紊流）,实验发现，圆管层流流动起始段的长度L*是雷诺数Re的函数，可按下式确定：,希累尔,L*0.2875dRe,布西内斯克,L*0.065dRe,第四节 圆管中的层流流动,一、圆管有效截面上的切应力分布.,1.取微元体：如图.半径,长 中心线和轴重合.,受力分析,两截面压力：,重力：,切向力：,3.在流动方向上的平衡方程.,由：,方程两边同除 得：,不随r变化,粘性流体在圆管中作层流流动时，同一截

19、面上的切向应力的大小与半径成正比。,注：此式同样适用于圆管中的紊流流动.,对水平管道：,在管壁上：,没有负号,由前述：,代如上式得：,二、速度分布.,根据牛顿内摩擦定律：,所以,旋转抛物面,三、最大流速:,四、平均流速:,五、流量:,圆管中的流量：,哈根一泊肃叶公式,选取管径的问题,经济流速,对于水平圆管：,或,六、达西公式:,由前述沿程损失公式：,可见，层流流动的沿程损失与平均流速的一次方成正比。,得,例 在管径，管长 的圆管中，冷冻机润滑油作层流运动，测得流量，水头损失，试求油的运动粘滞系数？,解：管中润滑油的平均流速,沿程阻力系数为,是层流,第五节粘性流体的紊流流动,一、紊流脉动与时均法

20、,紊流流动是极不规则的流动，这种不规则性主要体现在紊流的脉动现象。所谓脉动现象，就是诸如速度、压强等空间点上的物理量随时间的变化作无规则的随机的变动。在作相同条件下的重复实验时,所得瞬时值不相同,但多次重复实验的结果的算术平均值趋于一致。,时均速度,脉动速度,瞬时速度,同理,瞬时轴向速度与时均速度图,从工程应用的角度看,关心流体主流的速度分布、压强分布以及能量损失,流体主流的速度和压强，指的正是时均速度和时均压强,普通测速管的测量值均为平均值,空间各点的时均速度不随时间改变的紊流流动也称为定常流动或准定常流动,1.紊流中的切向应力,二、紊流中的切向应力 普朗特混合长度,对于,2.惯性切应力（雷

21、诺应力）,如图，在恒定流中时均速度沿x方向，脉动速度沿x和y方向的分量分别为 和,任取一水平截面A-A，设在某一瞬时，原来位于低流速层a点处的质点，以脉动速度 向上流动，穿过A-A截面到达 点。,则：,1）单位时间内通过A-A截面单位面积的流体质量为。,2）单位时间内通过单位面积的动量为,3）由动量定律，动量的变化率等于作用力。此时，动量变化率 通过截面A-A的动量流量。,作用力 沿x方向单位面积上的切向作用力惯性切应力。,（对取时均值）,3.惯性切应力的正负,当质点由下往上脉时，。由于a处X方向的时均速度 处x方向的时均速度，故，当质点由 时，会对该处原有的质点的运动起阻滞作用，产生负的沿x

22、方向的脉动流速；相反从上到下层会产生。但；无论哪一种情况：，为保证切应力非负：,2、普朗特混合长度理论,1）掺混类似于气体分子运动，而流体微团以 的速度自由的经过一段路程L，才与该层其他微团碰撞掺混。,Y,2）流体微团的纵向脉动速度 与横向脉动速度 的大小是属于同一个数量级。,3）脉动速度 与 与流层时均速度差 成正比。,观看录像,式中：为和的比例系数,令（L-混合长度）,则,三、圆管中紊流的速度分布和沿程损失,1.圆管中的紊流区划，粘性底层，水力光滑与水力粗糙,1）区域划分,2)水力光滑与水力粗糙,管壁粗糙凸出部分的平均高度叫做管壁的绝对粗糙度()/d 称为相对粗糙度,水力光滑,水力粗糙,光

23、滑管,粗糙管,2.圆管中紊流的速度分布,1)紊流光滑管,切应力常数,或：,2)紊流粗糙管,3.圆管中的沿程损失,紊流光滑管,紊流粗糙管,第六节沿程损失的试验研究,一、沿程阻力系数影响因素,研究沿程阻力系数，首先分析影响的因素:层流=64/Re，仅与Re有关，与管壁粗糙无关。,壁面粗糙在一定条件下成为产生惯性阻力的主要外因。,二尼古拉兹实验及尼古拉兹曲线,确定阻力系数与雷诺数Re及相对粗糙度/d 之间的关 系，具体关系要由实验确定，最著名的是尼古拉茨于19321933年间做的实验。,1.实验方法：人为造出六种不同的相对粗糙度的管；对不同的管径通过改变流量来改变雷诺数；测出沿程阻力损失，由 求阻力

24、系数.,2.实验结果:,观看动画,3.阻力分区:,1）层流区：（Re2320）不论 如何变化，都集中在一条直线上。-表明仅随Re，与相对粗糙度无关。（此为层 流运动，证明了理论推导的结果）,2）过渡区（2320Re4000）实验点比较分散，层流向紊流过渡的不稳定区域。,不同相对粗糙点，起初都集中在一条直线上-紊流光滑区。（当Re，逐渐偏离，较小，Re较大时才偏离）,5）紊流粗糙管平方阻力区,2）在过渡区，层流底层变薄，粗糙开始影响到核心区内流动，加大了核心区紊流强度，因此增加了阻力和能量损失，,1）在光滑区，粗糙突起高度k比层流底层小得多，,说明,3）紊流粗糙区，层流底层更薄，粗糙突起高度几乎

25、全部暴露在紊流核心中，,尼古拉兹实验比较完整地反映了沿程损失系数的变化规律，揭示了影响变化的主要因素，对和断面流速分布的测定，推导紊流的半经验公式提供了可靠的依据。,三、莫迪图（用于计算新的工业管道）,根据普朗特的半经验理论，以及尼古拉兹实验曲线得到。,莫迪图对计算新的工业管道的沿程损失系数很方便。,柯列布茹克公式,柯氏公式是在合并两个半经验公式的基础上获得的，可以认为该公式是普朗特理论的尼古拉兹实验结合后进一步发展到工程应用阶段的产物，该公式在国内外得到了极为广泛的应用。柯氏公式的求解相对复杂，一般采用计算机数值计算方式。为了简化计算,莫迪以柯氏公式为基础绘制出反映Re、k/d与对应关系的莫

26、迪图,在图上可根据Re、k/d直接查出（如下图）：,例 在直径，相对粗糙度 的工业管道内，运动粘滞系数，的水以的速度运动。试求：管长的管道内的沿程水头损失。解：1）查莫迪图：流动处于紊流粗糙区2）用尼古拉兹粗糙区公式,两种方法较为接近,例沿程损失：已知管道和流量求沿程损失,已知:d200mm,l3000m 的旧无缝钢管,900 kg/m3,Q90T/h,在冬天为1.092 10-4 m2/s,夏天为0.355 10-4 m2/s。,求：冬天和夏天的沿程损失hf,解：,在夏天，查旧无缝钢管等效粗糙度=0.2mm,/d=0.001查穆迪图2=0.0385,例 沿程损失：已知管道和压降求流量,已知:

27、d10cm,l400m 的旧无缝钢管比重为0.9,=10-5 m2/s 的油，,求：管内流量qv,解：,Moddy图完全粗糙区的0.025,设10.025,由达西公式,查Moddy图得20.027,重新计算速度,查Moddy图得20.027,例 沿程损失：已知沿程损失和流量求管径,求：管径d 应选多大,解：,由达西公式,由/d=0.2/98.4=0.002，查Moody图得2=0.027,d 2=(3.6910 4 0.027)1/5=0.0996(m),Re2=4000/0.0996=4.01104,/d=0.2/99.6=0.002，查Moody图得3=0.027,取d=0.1m。,用迭代

28、法设1=0.025,第七节非圆管的沿程损失,输送流体的管道不一定都是圆形截面。,对于这些非圆形管道的沿程损失计算问题，达西公式和雷诺数的计算公式仍然可以应用。,但要把公式中的直径d用当量直径D来代替。,圆管：当量直径公式则矩形：，其矩形当量直径同样，非圆管道Re和k/d分别为：,矩形：,此时，,说明：,1.试验表明，在使用当量直径原理计算时，对矩形 三角行 方形 的计算结果和试验结果较接近，在对和圆形差别比较大的形状，计算结果就不可靠。,2.由于层流和紊流的流速分布不同，沿程损失不像紊流那样集中在管壁附近，这样单纯用湿周大小作为影响能量损失的主要外因条件，对层流来说就不充分了。,第八节 局部损

29、失,流体经过阀门、弯管、突扩和突缩等管件,流体经过这些局部件时，由于通流截面、流动方向的急剧变化，引起速度场的迅速改变，增大流体间的摩擦、碰憧以及形成旋涡等原因,从而产生局部损失。,损失均按 计算，关键是如何确定,一、突然扩大,1.损失机理,速度分布变化附加摩擦碰撞 漩涡,根据连续方程有:,根据动量方程有:,2.求局部损失系数,由伯努利方程,比较得,整理得,求,由,得,讨论：若（如管道流入很大的容器（水池）或气体流入大气）（速度头完全损失）。,二、其它局部阻力系数1）其它截面变化(如突然缩小,渐扩管)2）弯管3）管道阀件4）三通（如合流，分流，Y型，T型等）5)其它(如过滤网,波纹管),计算管

30、段的总能量时，应将管段上所有的沿程损失和局部损失算术求和。,三、减小阻力的措施,1、改善固体边界状况,1）增大过流断面几何尺寸、dhf虽然从减小阻力角度采用大管径，但费用会增加，同时有一些用途（如除尘管道），有一个最小风速的要求。因此，管道直径由技术经济比较来确定。,2）减少管长，管越短越好，尽可能采用直管道。,3）减少局部管件，凡是能够不要的管件尽可能不要。,4）提高管壁光滑度。,5）改变局部管件结构-减少局部损失。,进口：光滑喇叭口,突然扩大、缩小：只要布置可能，避免突扩突缩可用渐扩、渐缩或台阶形式。,弯管：采用导叶，使局部从1.00.3。,2、添加剂,流动中加入极少量添加剂，改善结构hf

31、、hm两种。高分子聚合物，如聚氧化已烯等；金属皂，如碱金属皂、铵皂。,第九节 各类管流水力计算,一、简单管道,管径和管壁粗糙度相同的一根管子或这样的数根管子串联在一起的管道系统叫简单管道。,简单管道有三类计算问题：,二、串联管道,由不同直径或粗糙度的数段管子连接在一起的管道叫做串联管道,通过串联管道各管段的流量是相同的。,串联管道的损失应等于各管段损失的总和。,串联管道有两类计算问题:,由连续方程,是由管道尺寸和局部损失系数确定的已知数。,由,已知,试取,计算,查莫迪图,N,Y,计算,三、并联管道,在某处分成几路、在下游某处又汇合成一路的管道叫并联管道。,并联管道的总流量等于个分管道的流量的总

32、和。,并联管道的损失等于个分管道的损失。,按简单管道计算,并联管道的计算问题,具体步骤,例：如图所示的具有并联、串连管路的虹吸管，已知H40m，l1200m，l2100m，l3500m，d10.2m，d20.1m，d30.25m，120.02，30.025，求总流量Q。,【解】管1和管2并联，此并联管路又与管3串连，因此,将 已知数值代入上式，计算得,第十一节 液体的出流,一、孔口的分类,薄壁孔口,一、薄壁孔口定常出流,1.薄壁小孔口定常出流,2.薄壁大孔口定常出流,液体经大孔口的淹没出流，流速和流量仍按上式计算，差别主要在于出流系数不同。,如果壁厚 或侧壁孔口处接一段长的圆管时，此时流动情况

33、为管嘴出流。,二、外伸管嘴定常出流,特点：水流入嘴嘴时如同孔口出流一样，管流股也发生收缩，存在着收缩断面CC，尔后流股逐渐扩张，至出口断面上完全充满管嘴断面流出。,1、管嘴出流流量,以管嘴中心线为基准线，列A-A及B-B断面伯努利方程：,令,则,管嘴流量,由于出口断面B-B流股完全充满（不同于孔口）,则,2、管嘴内的真空度,在C-C断面前后流股与管壁分离，中间形成涡旋区，产生负压，出现了管嘴真空现象。在H0（作用水头）中，压差项 里，pC（绝对压强）小于大气压，从而使H0增大，促使出流流量增大。,三、各种管嘴的出流系数,第七章 流体多维运动学基础,第七章 流体多维运动学基础,本章内容：在许多工

34、程实际问题中，流动参数不仅在流动方向上发生变化，而且在垂直于流动方向的横截面上也要发生变化。要研究此类问题，就要用多维流的分析方法。本章主要讨论理想流体多维流动的基本规律，为解决工程实际中类似的问题提供理论依据，也为进一步研究粘性流体多维流动奠定必要的基础。,第一节 流体运动的连续性方程,当把流体的流动看作是连续介质的流动，它必然遵守质量守恒定律。对于一定的控制体，必须满足输运公式。它表示在控制体内由于流体密度变化所引起的流体质量随时间的变化率等于单位时间内通过控制体的流体质量的净通量。,一、笛卡儿坐标系中微分形式的连续性方程。,设该微元六面体中心点O（x,y,z）上流体质点的速度为、，密度为

35、，于是和 轴垂直的两个平面上的质量流量如图所示。,在 x 方向上，单位时间通过EFGH面流入的流体质量为：,（a）,单位时间通过ABCD面流出的流体质量：,（b）,则在 方向单位时间内通过微元体表面的净通量为（b）-（a），即,（c1）,第八章 理想流体的有旋流动和无旋流动,（c2）,（c3）,因此，单位时间流过微元体控制面的总净通量为:,（c）,微元六面体内由于密度随时间的变化而引起的质量的变化率为：,（d）,同理可得 和 方向单位时间通过微元体表面的净通量分别为:,将式（c），（d）代入下式，取 0，则可得到流场中任一点的连续性方程的一般表达式为:,连续性方程表示了单位时间内控制体内流体质

36、量的增量等于流体在控制体表面上的净通量。,在定常流动中，由于,对于不可压缩流体（=常数）,二、柱坐标系中的连续性方程。,【解】对不可压缩流体连续性方程为：,将已知条件代入上式，有,一、物理模型,第二节 流体微团的运动分析,F点速度:,方程两边加上两个和为0的项，其值不变,组合,注：此关系也称为亥姆霍兹(Helmholtz)速度分解定理，该定理可简述为：在某流场O点邻近的任意点A上的速度可以分成三个部分，与O点相同的平移速度（移运）；绕O点转动在A点引起的速度（旋转运动）；由于变形（包括线变形和角变形）在A点引起的速度（变形运动）。,二、物理意义（以平面流动进行分析）,1.平移运动,向左移动,向

37、上移动,流体力学,2.线变形运动,每秒内单位长度的伸长（或缩短）量称为线应变速度,同理y向线变形速度：,B、C在x方向有速度差,经过dt时间BC边伸长,单位时间单位长度的伸长：,3.角变形运动,B、C在y方向有速度差：,在dt时间内BC线段将旋转：,同理，AB在dt时间线段将旋转：,单位时间内直角ABC变成锐角ABC，变形速度为：,定义XY平面的剪切变形率为：,同理可得：,4.旋转运动,流体微团的旋转角速度的定义为每秒内绕同一转轴的两条互相垂直的微元线段旋转角度的平均值。,规定逆时针旋转角度为正：,BC边旋转的角度为：,BA边旋转的角度为：,轴旋转角速度为两个互相垂直边旋转角速度的一半：,写成

38、矢量形式为：,根据流体微团在流动中是否旋转，可将流体的流动分为两类：有旋流动和无旋流动。,三、有旋运动和无旋运动,在笛卡儿坐标系中：,即当流场速度同时满足：,时流动无旋,需要指出的是，有旋流动和无旋流动仅由流体微团本身是否发生旋转来决定，而与流体微团本身的运动轨迹无关。,流体微团的运动分析的意义：1.将刚体的旋转运动从一般的运动中分离出来，可将流体的 运动分为有旋和无旋两种流动。2.将流体的变形率与内应力联系起来。,【例】某一流动速度场为，其中a是不为零的常数，流线是平行于x轴的直线。试判别该流动是有旋流动还是无旋流动。,【解】由于,所以该流动是有旋运动。,第三节无旋流动,一 速度势函数,对于

39、无旋流场，处处满足：，由矢量分析知，任一标量函数梯度的旋度恒为零，所以速度一定是某个标量函数的梯度，即:,速度势函数 速度势,1.速度势函数的表达式,对于柱面坐标,结论：,无旋条件是速度有势的充要条件。无旋必然有势，有势必须无旋。所以无旋流场又称为有势流场。速度势的存在与流体是否可压缩、流动是否定常无关。,1）速度势函数在任意方向的偏导数等于速度在该方向的投影；,2.速度势的特性,2）流速（流线）与等势面正交，且指向势函数增加的方向；,3）在单连通区域中，任意两点的速度势函数的差值等于该两点间的流速沿任意曲线积分；,有势流动中沿AB曲线的速度线积分等于终点B和起点A的速度势之差。由于速度势是单

40、值的，则该线积分与积分路径无关。,在有势流动中当速度沿封闭周线积分时，周线上的速度环量等于零。,如何根据速度分布求解势函数？,根据无旋条件，速度有势：,代入不可压缩连续性条件可得：,拉普拉斯方程,拉普拉斯算子,4）对于不可压缩流体，速度势是调和函数，满足拉普拉斯方程。,柱坐标系下：,求解不可压缩流体无旋流动问题，便归纳为根据起始条件和边界条件求解拉普拉斯方程问题。,二 流函数,在笛卡儿坐标系中，平面不可压缩流体的连续性方程可写成:,此外，平面流动的流线微分方程为,由数学知识可知，此式是成为某个函数全微分的必要且充分条件,1.流函数的表达式：,若定义某一个函数（流函数）:,即函数永远满足连续方程

41、。很显然，在流线上 0或常数。在每条流线上函数都有它自己的常数值，所以称函数为流函数。,2.流函数的基本性质,1）等流函数线为流线,当 常数时,等流函数线为流线，每条流线有各自的常数值。,2）流体通过两流线间单位高度的体积流量等于两条流线的流函数之差,在xy平面上任取A和B点，AB连线如图所示，则,由不可压缩流体、平面、无旋流动条件有：,将速度和流函数的关系代入上式得,在极坐标系中：,故不可压缩流体的平面无旋流动流函数也满足拉普拉斯方程，也是调和函数。,3）对于不可压缩流体，流函数是调和函数，满足拉普拉斯方程；,4）速度势函数和流函数的关系（流网）,对于不可压缩流体的平面无旋流动，速度势函数和

42、流函数都是调和函数，且具有以下关系：,该数学关系式称为柯西黎曼（CauchyRiemen）条件。由它可得：,上式是等势线簇和流线簇互相垂直的条件，即正交性条件。在平面上它们构成处处正交的网络，称为流网。,例：试证明不可压缩流体平面流动,能满足连续方程，是一个有势流动，并求出速度势。,解：,满足连续方程：,流动为有势流动,本章结束！,第八章 理想流体的有旋流动和无旋流动,第八章 理想不可压缩流体的有旋流动和无旋流动,本章内容：在许多工程实际问题中，流动参数不仅在流动方向上发生变化，而且在垂直于流动方向的横截面上也要发生变化。要研究此类问题，就要用多维流的分析方法。本章主要讨论理想流体多维流动的基

43、本规律，为解决工程实际中类似的问题提供理论依据，也为进一步研究粘性流体多维流动奠定必要的基础。,第一节 微分形式的连续性方程,当把流体的流动看作是连续介质的流动，它必然遵守质量守恒定律。对于一定的控制体，必须满足输运公式。它表示在控制体内由于流体密度变化所引起的流体质量随时间的变化率等于单位时间内通过控制体的流体质量的净通量。,首先推导在笛卡儿坐标系中微分形式的连续性方程。,设该微元六面体中心点O（x,y,z）上流体质点的速度为、，密度为，于是和 轴垂直的两个平面上的质量流量如图所示。,在 x 方向上，单位时间通过EFGH面流入的流体质量为：,（a）,单位时间通过ABCD面流出的流体质量：,（

44、b）,则在x方向单位时间内通过微元体表面的净通量为（b）-（a），即,（c1）,第八章 理想流体的有旋流动和无旋流动,第八章 理想流体的有旋流动和无旋流动,（c2）,（c3）,因此，单位时间流过微元体控制面的总净通量为:,（c）,微元六面体内由于密度随时间的变化而引起的质量的变化率为：,（d）,同理可得 和 方向单位时间通过微元体表面的净通量分别为:,将式（c），（d）代入下式，取 0，则可得到流场中任一点的连续性方程的一般表达式为:,连续性方程表示了单位时间内控制体内流体质量的增量等于流体在控制体表面上的净通量。,在定常流动中，由于,对于不可压缩流体（=常数）,【解】对不可压缩流体连续性方程

45、为：,将已知条件代入上式，有,一、物理模型,第二节 流体微团的运动分析,F点速度:,方程两边加上两个和为0的项，其值不变,组合,注：此关系也称为亥姆霍兹(Helmholtz)速度分解定理，该定理可简述为：在某流场O点邻近的任意点A上的速度可以分成三个部分，与O点相同的平移速度（移运）；绕O点转动在A点引起的速度（旋转运动）；由于变形（包括线变形和角变形）在A点引起的速度（变形运动）。,二、物理意义（以平面流动进行分析）,1.平移运动,向左移动,向上移动,流体力学,2.线变形运动,每秒内单位长度的伸长（或缩短）量称为线应变速度,同理y向线变形速度：,B、C在x方向有速度差,经过dt时间BC边伸长

46、,单位时间单位长度的伸长：,3.角变形运动,B、C在y方向有速度差：,在dt时间内BC线段将旋转：,同理，AB在dt时间线段将旋转：,单位时间内直角ABC变成锐角ABC，变形速度为：,定义XY平面的剪切变形率为：,同理可得：,4.旋转运动,流体微团的旋转角速度的定义为每秒内绕同一转轴的两条互相垂直的微元线段旋转角度的平均值。,规定逆时针旋转角度为正：,BC边旋转的角度为：,BA边旋转的角度为：,轴旋转角速度为两个互相垂直边旋转角速度的一半：,写成矢量形式为：,根据流体微团在流动中是否旋转，可将流体的流动分为两类：有旋流动和无旋流动。,三、有旋运动和无旋运动,在笛卡儿坐标系中：,即当流场速度同时

47、满足：,时流动无旋,需要指出的是，有旋流动和无旋流动仅由流体微团本身是否发生旋转来决定，而与流体微团本身的运动轨迹无关。,【例】某一流动速度场为，其中a是不为零的常数，流线是平行于x轴的直线。试判别该流动是有旋流动还是无旋流动。,【解】由于,所以该流动是有旋运动。,第三节 理想流体的运动方程 定解条件,在x方向：,一、运动微分方程,两端同除以微元体的质量,将加速度的表达式代入,其矢量式为：,公式为理想流体运动微分方程式，物理上表示了作用在单位质量流体上的质量力、表面力和惯性力相平衡。该式推导过程中对流体的压缩性没加限制，故可适用于理想的可压流体和不可压缩流体，适用于有旋流动和无旋流动。,将上式

48、作恒等变形，便可以直接由运动微分方程判定流动是有旋还是无旋流动，在式的第一式右端同时加减、，得:,所以,写成矢量形式,如果流体是在有势的质量力作用下，流场是正压性的，则：,将压强函数对坐标的偏导数有:,将上述关系代入式得:,二、定解条件,方程组的封闭问题,还需要增添一个方程，使方程组封闭，才能求解。,对于不可压缩流体，,对于密度仅是压强的函数的流体,方程组的定解条件,1.初始条件,初始条件是指在起始瞬时t0所给定的流场中每一点的流动参数。,定常流动不需要给定初始条件。,2、边界条件,边界条件是指任一瞬时运动流体所占空间的边界上必须满足的条件。,固体壁面,流体既不能穿透壁面，也不能脱离壁面而形成

49、空隙，即流体与壁面在法线方向的相对分速度应等于零。即：,若固壁是静止的：,在两种不同流体交界面上：,流体交界面,无穷远处,一般给定无穷远处流体的速度、压强和密度。,流道进、出口处,一般给定进、出口处的速度分布、温度、压力等。,一、欧拉积分,当理想正压性流体在有势的质量力作用下作定常无旋流动时，下式右端为零。若在流场中任取一有向微元线段，其在三个坐标轴的投影分别为dx、dy、dz，将它们分别依次乘下式并相加，得:,第四节 理想流体运动方程的积分式,积分,上式为欧拉积分的结果，表明理想正压性流体在有势的质量力作用下作定常无旋流动时，单位质量流体的总机械能在流场中保持不变。,三、伯努利积分,当理想正

50、压性流体在有势的质量力作用下作定常有旋流动时，下式右端第一项等于零。由流线的特性知，此时流线与迹线重合，在流场中沿流线取一有向微元线段，其在三个坐标轴上的投影分别为，将它们的左、右端分别依次乘下式的左、右端，相加有,该积分为伯努利积分。表明理想正压性流体在有势的质量力作用下作定常有旋流动时，单位质量流体的总机械能沿流线保持不变。通常沿不同流线积分常数值有所不同。,第五节 涡线、涡管、涡束和涡通量,一、涡线、涡管、涡束,涡线,涡管,涡束,涡通量,在有旋流动流场的全部或局部区域中连续地充满着绕自身轴线旋转的流体微团，于是形成了一个用角速度 表示的涡量场（或称角速度场）。,1.涡线,涡线是一条曲线，