流体动力学基础 流体力学

1、第四章第四章 流体动力学基础流体动力学基础 第一节第一节 流体的运动微分方程流体的运动微分方程 第二节第二节 元流的伯努利方程元流的伯努利方程 第三节第三节 总流的伯努利方程总流的伯努利方程 第四节第四节 总流的动量方程总流的动量方程 第五节第五节 理想流体的无旋流动理想流体的无旋流动 第一节第一节 流体的运动微分方程流体的运动微分方程 一、理想流体运动微分方程一、理想流体运动微分方程 在运动的理想流体中，取微小平行六面体(质点)， 正交的三个边长dx,dy,dz,分别平行于x,y,z坐标轴(图4 1)。设六面体的中心点o，速度压强，分析该微 小六面体方向的受力和运动情况。 1. 1.表面力：

2、表面力：理想流体内不存在切应力只有压强方 向受压面(abcd面和abcd面)形心点 图 图4141连续性微分方程连续性微分方程 的压强为： 受压面上的压力为： 质量力： 由牛顿第二定律 得： ( ) -( ) + dt du dxdydz x dxpp x p M 2 1 dxpp x p N 2 1 dydzpP MM dydzpP NN dxdydzXFBx dt du x x mF dxp x p 2 1 dxp x p 2 1 dydz dxdydzX 化简得： （41） 将加速度项展成欧拉法表达式 : （42） dt du z p dt du y p dt du x p z y x

3、Z Y X 1 1 1 z u zy u yx u xt u z p z u zy u yx u xt u y p z u zy u yx u xt u x p zzzz yyyy xxxx uuuZ uuuY uuuX 1 1 1 uupf t u 1 用矢量表示为： （43） 上式即理想流体运动微分方程式，又称欧拉运动 微分方程式。该式是牛顿第二定律的表达式，因此是 控制理想流体运动的基本方程式。 1755年欧拉在所著的流体运动的基本原理中 建立了欧拉运动微分方程式，及上一节所述的连续性 微分方程式。对于理想流体的运动，含 有和 四个未知量，由式(330)和式(336)组成的基本方 程组，

4、满足未知量和方程式数目一致，流动可以求解。 因此说，欧拉运动微分方程和连续性微分方程奠定了 理想流体动力学的理论基础。 zyxuuu, 二、粘性流体运动微分方程二、粘性流体运动微分方程 1 1粘性流体的动压强粘性流体的动压强 理想流体因无粘滞性，运动时不出现切应力，只 有法向应力，即动压强。用类似分析流体静压强 特性的方法，便可证明任一点动压强的大小与作用 面的方位无关，是空间坐标和时间变量的函数， 即 （，）。 粘性流体的应力状态和理想流体不同，由于粘性 作用，运动时出现切应力，使任一点的法向应力的 大小与作用面的方位有关。如以应力符号的第个 下角标表示作用面的方位， 第二个角标表示应力的方

5、向，则法向应力 进步研究证明，任一点任意三个正交面上的法向应力之和 都不变，即 据此，在粘性流体中，把某点三个正文面上的法向应力的平 均值定义为该点的动压强以p表示： （44） 如此定义，粘性流体的动压强也是空间坐标和时间变量的函 数 。 zzyyxx ppp pppppp zzyyxx tzyxpp, zzyyxx pppp 3 1 2. 2. 应力和变形速度的关系应力和变形速度的关系 粘性流体的应力与变形速度有关，其中法向应力与线变形 速度有关，切应力则与角变形速度有关 流动中某点的动压强是过该点三个相互正交平面上法向应力 的平均值，同某一平面上的法向应力有一定差值，称为附加法向应 力，以

6、表示，它是流体微团在法线方向上发生线变形 (伸长或缩短)引起的。 （45） 切应力与角变形速度的关系，在简单剪切流动中符合牛顿内 摩擦定律 zzyy ppp xx , z u zzzz y u yyyy x u xxxx z y x pppp pppp pppp 2 2 2 dy du u 将牛顿内摩擦定律推广到一般空间流动,得出 （46） 3 3粘性流体运动微分方程粘性流体运动微分方程 采用类似于推导理想流体运动微分方程式（41）的方法， 取微小平行 六面体，根据牛顿第二定律建立以应力(包括切 应力)表示的运动微分方程式，并以式(45)、式(46)代人 整理，使得到粘性液体运动微分方程： y

7、 u x u yxxy x u z u xzzx z u y u zyyz x y zx y z 用矢量表示为 式中： 拉普拉斯算子。 粘性流体运动微分方程，又称为纳维 斯托克斯方程(简写为NS方程)。 z u zy u yx u xt u zz p z u zy u yx u xt u yy p z u zy u yx u xt u xx p zzzz yyyy z y xx uuuuZ uuuuY uuuuX 2 1 2 1 2 1 uuupf t u 2 1 2 2 2 2 2 2 2 zyx （48） NS方程表示作用在单位质量流体上的质量力、表面力(压力 和粘性力) 的相平衡。由NS

8、方程式和连续性微分方程式组成的基 本方程组，原则上可以求解速度场和压强场p,可以说粘性流体的运 动分析，归结为对NS方程的研究。 例41 理想流体速度场为 为常数。试求：（1）流动是否可能实现；（2）流线方程；（3）等 压面方程(质量力忽略不计) 解 (1)由连续性微分方程 满足连续性条件，流动是可能实现的。 (2)由流线方程 得 ： baubxuayuzyx,0, 0 z u y u x uzyx yxu dy u dx bx dy ay dx aydybxdx 积分得流线方程 a,b同号，流线是双曲线a,b异号,流线是圆。 (3)由欧拉运动微分方程式，不计质量力： 将方程组化为全微分形式:

9、 aby x u u abx y u u y xy p x yx p 1 1 )( 1 )()( 1 ydyxdxabdp ydyxdxabdy y p dx x p caybx 22 积分，得 令p=常数 即得等压面方程 等压面是以坐标原点为中心的圆。 2 22 c yx abp cyx 22 第二节第二节 元流的伯努利方程元流的伯努利方程 一、理想流体运动微分方程的伯努利积分一、理想流体运动微分方程的伯努利积分 理想流体运动微分方程式是非线性偏微分方程组，只有特 定条件下的积分，其中最为著名的是伯努利(Daniel Bernoull， 17001782，瑞士科学家)积分。 （410） z

10、u zy u yx u xz p z u zy u yx u xy p z u zy u yx u xx p zzz yyy z y x uuuZ uuuY uuuX 1 1 1 由理想流体运动微分方程式 各式分别乘以沿流线的坐标增量dx，dy，dz，然后相加得： 、 1.引人限定条件： 作用在流体上的质量力只有重力：X=Y=0,Z=-g; dzdydxZdzYdyXdx z p y p x p 1 )( dz dt du dy dt du dx dt du z y x gdzZdzYdyXdx)( dt du z p dt du y p dt du x p z y x Z Y X 1 1 1

11、 .不可压缩，恒定流：,Czyxpp, p z p y p x p ddpdzdydx 11 .恒定流流线与迹线重合：dx=uxdt，dy=uydt，dz=uzdt 则 2 222 zyx z y x uuu ddz dt du dy dt du dx dt du 相加带入后得： dzdydxZdzYdyXdx z p y p x p 1 )( dz dt du dy dt du dx dt du z y x g up z g up z 22 2 22 2 2 11 1 Cgz u g p 2 2 Cz g u p 2 2 理想流体运动微分方程沿流线的积分称为伯努利积分， 由于元流的过流断面积

12、无限小，所以沿流线的伯努利方 程就 是元流的伯努利方程。推导该方程引入的限定条件，就 是理想流体元流伯努利方程的应用条件，归纳起来有：理想 流体；恒定流动；质量力中只有重力；沿元流(流线)；不可压 缩流体。 1. 1.物理意义式物理意义式 式子中的前两项 、 和 的物理意义分 别是单位重量流体具有的比位能压能或比势能；单位重量 流体具有的动能。 g p z g u 2 2 g p z , 三项之和 是单位重量流体具有的机械能，沿 同一无流(沿同一流线)。单位重量流体的机械能守恒。伯努利 方程又称为能量方程。 g u g p z 2 2 2. 2.流体意义流体意义 各项的流体力学意义为：z是位置

13、水头， 压强水头； 两项之和 是测压管水头， 是流速水头，三项之和 称为总水头 g p g u 2 2 g p zHp g u g p zH 2 2 图42水头线 表示理想流体的恒定流动，沿同一元流(沿同一流线)各断式 (423)则面的总水头相等理想流体的水头线是水平线 3. 3.几何意义 几何意义 各项的几何意义是不同的几何高度：z是位置高度，测压管高 度。总结如下： p 例例42 应用皮托(Pito，H.)管测量点流速 前文指出，流速水头可直接量测，现以均匀管流为例加以说明。 设均匀管流，欲量测过流断面上某点A的流速(图43)。在该点放置 一根两端开口，前端弯转90的细管，使前端管口正对来

14、流方向， 另一端垂直向上，此管称为测速管。来流在A点受测速管的阻滞速 度为零，动能全部转化为压能测速管中液面升高。 另在A点上游的同一流线取相距很近的o点，因这两点相距很 近，o点的压强p实际上等于放置测速管以前A点的压强 应用理想流 体元流伯努利方程： （425） （426） 0 2 2 h g p g p g u g p g u g p 2 2 图43点流速的测量 式中o点的压强水头，由另根测压管量测， 于是测速管和测压管中液面的高度差,就是A 点的流速水头，该点的流速： （427） 根据上述原理，将测速管和测风管组合 成测量点流速的仪器，图44所示，与迎流 孔(测速孔)相通的是测速管，与

15、侧面顺流孔 （测压孔或环形窄缝）相通的是测压管。考 虑到粘性流体从迎流孔至顺流孔存在粘性效 应，以及皮托管队员流场的干扰等影响，引 用修正系数C： 02 2ghC g pp gCu 图44 毕托管构造 02 2gh g pp gu 式中C是修正系数数值接近于1.0，由实验测定。 【例4-3】 有一贮水装置如图（4-5）所示，水池足够大，当阀 门关闭时，压强计读数为2.8个大气压强。而当将阀门全开，水从管 中流出时，压强计读数是0.6个大气压强，试求当水管直径d=12cm 时，通过出口的体积流量(不计流动损失)。 【解解】 当阀门全开时列1-l、2-2截面的伯努利方程 当阀门关闭时，根据压强计的

16、读数， 应用流体静力学基本方程 ， 求出值： g V g pp g p H aaa 2 6.0 00 2 2 图45 aaa ppgHp8 . 2 OmH g p H a 2 28 9806 980608 .28 .2 所以管内流量： sm g p HgV a /78.20 9806 980606 .0 8 .2806.92 6 .0 2 2 smVdq V /235.078.20785.0 4 3 2 2 三、粘性流体元流的伯努利方程三、粘性流体元流的伯努利方程 能量守恒原理：能量可以从一种形式转换成另一种形式， 既不能创造、也不能消灭，总能量是恒定的 粘性流体流动时，单位重量流体具有的机械

17、能沿程减少， 总水头线是沿程下降。 因此，设为粘性流体元流单位重量流体由过流断面11 运动至过流断面22的机械能损失，称为元流的水头损失， 根据能量守恒原理，便可得到粘性流体元流的伯努利方程 水头损失 也具有长度的量纲。 w h g up z 21 2 11 22 2 22 wg up hz 第三节第三节 总流的伯努利方程总流的伯努利方程 上一节的最后得到了粘性液体元流的伯努利方程式，为了 解决实际问题，还需要将其推广到总流中去。 一、渐变流及其性质一、渐变流及其性质 渐变流： 凡质点的迁移加速度(位变加速度)很小，的流动， 或者说流线近于平行直线的流动渐变流是均匀流的宽延，所以 均匀流的性质

18、，对于渐变流都近似成立，主要是： 1渐变流的过流断面近于平面。面上各点的速度方向近于 平行； 2恒定渐变流过流断面上的动压强按静压强的规律分布，即： cz p 二、总流的伯努利方程二、总流的伯努利方程 设恒定总流，过流断面11、22为渐变流断面，面积为 A1,A2由元流的伯努利方程： 图47急变流和渐变流 g up z 21 2 11 22 2 22 wg up hz 单位时间通过元流两过流断面的能量关系 总流是由无数元流构成的，上式对总流过流断面积分 便得到单位时间通过总流两过流断面的总能量关系 分别确定三种类型的积分 第一类积分： 因所取过流断面是渐变流断面 dQz g up 21 2 1

19、1 dQz g up 22 2 22 dQhw 1 1 111 A p dAuz 1 2 1 112 A g u dAu 2 2 222 A p dAuz 2 2 2 222 A g u dAu Q w aQh A p udAz cz p 第二类积分： 动能校正系数， 式中 是为校正以断面平均速度计算的动能与实际功 能的差异而引入的校正系数， 值取决于过流断面上 的流速分布情况，分布均匀的流动 。 通常取 A g u udA 2 2 A g u udA 2 2 A g u dA 2 3 Q g v 2 2 Av dAu dA dA A A g v A g u 3 3 2 3 2 3 10. 1

20、05. 11 Q p zudAz A p 第三类积分： 为总流单位重量流体由11至22断面的平均机械 能损失，称总流的水头损失 两断面间无分流及汇流，Q1Q2Q，并以 除 上式，得 Q w dQh w h QhdQh w Q w QhQ g v Q p z Q g v Q p z w 2 2 2 22 2 2 11 1 2gQ wg vp g vp hzz 2221 2 222 2 111 2. 2. 伯努利方程的适用条件伯努利方程的适用条件 式式(437)(437)即粘性流体总流的伯努利方程。将元流的伯努利方即粘性流体总流的伯努利方程。将元流的伯努利方 程推广为总流的伯努利方程，引入了某些限

21、制条件，也就是程推广为总流的伯努利方程，引入了某些限制条件，也就是 总流伯努利方程的适用条件包括：总流伯努利方程的适用条件包括： . .不可压缩流体恒定流；不可压缩流体恒定流； . .质量力只有重力；质量力只有重力； 不可压缩流体不可压缩流体( (以上引自粘性流体元流的伯努利方程以上引自粘性流体元流的伯努利方程) )； . .所取过流断面为渐变流断面；所取过流断面为渐变流断面； . .两断面间无分流和汇流；两断面间无分流和汇流； . .两断面间无能量的输入或支出；两断面间无能量的输入或支出； . .不存在相对运动。不存在相对运动。 3. 伯努利方程的方法步骤伯努利方程的方法步骤 . .断面选择

22、断面选择 通常选择未知量所在的断面和已知量最多的断面， 它们 都必须是渐 变流断面； . .代表点选择代表点选择 无压流一般选择自由液面，有压流一般选在管 道中心； . .位置基准面选择位置基准面选择 习惯选择在过各代表点最低者的水平面。 位置准面选择对结果无影响； . .压强基准面选择压强基准面选择 液体一般选取相对压强；气体一般选取绝 对压强。压强准面选择对结果无影响； . .列伯努利方程列伯努利方程 对于初学者，应该分项列出，哪怕是零，也 应该写 出。但一般只用符号代替，而不代入具体数值，以 便推导出未知量的计算公式； . .解伯努利方程解伯努利方程 求解出题目中所要求的未知量； . .

23、给出答案给出答案 给出正确的答案 例例43 用直径d100mm的水管从水箱引水(图49)。水箱水面与 管道出口断面中心的高差H4m保持恒定，水头损失 3m水柱。 试求管道的流量。 解解 这是一道简单的总流问题，应用伯努利方程： 图49管道出流 w h g vp z 21 2 11 wg vp hz 22 2 22 求解的关键是“三选”：选基准面、计算断面和计算点。为便 于计算，选通过管道出口断面中心的水平面为基准面00(图4 9)。计算断面应选在渐变流断面，并使其中一个已知量最多，另 一个含待求量。技以上原则本题选水箱水面为11断面，计算点 在自由水面上、运动参数z1=H,p1=0 (相对压强

24、), v1=0 。选管道出 口断面为22断面，以出H断面的中心运动参数z2=0,p2=0, v2待求。 将各量代人总流伯努利方程： 取 得： wg v hH 2 2 2 0 . 12 smhHgvw/43. 4)(2 2 四、总流伯努利方程应用的修正四、总流伯努利方程应用的修正 伯努利方程是古典水动力学应用最广的基本方程。应用伯努 利方程要重视方程的应用条件，切忌不顾应用条件，随意套用公 式，要对实际问题做具体分析，灵活运用。下面结合三种情况加 以讨论。 1.气体的伯努利方程 总流的伯努利方程式(436)是对不可 压缩流体导出的，气体是可压缩流体，但 是对流速不很大(60ms)，压强变化不 大

25、的统，如工业通风管道、烟道等，气流 在运动过程中密度的变化很小，在这样的 条件下，伯努利方程仍可用于气流。由于 气流的密度同外部空气的密度是相同的数量级，在用相对压强进 行计算时，需要考虑外部大气压在不同高度的差值。 设恒定气流(图410)、气流的密度为 外部空气的密度为 ， 过流断面上计算点的绝对压强 。 列11和22断面的伯努利方程式： a absabsPP21, 图410恒定气流 wg vp hz 22 2 22 g vp z 21 2 11 （438） 进行气流计算，通常把上式表示为压强的形式 （439） 式中pw为压强损失 （440） 将式(439)中的压强用相对压强p1,p2表示，

26、则： （441） （442） 式中 为 处的大气压， 为高程 处的大 压，代人式(437)，整理得： （443） 1 21 2 2 1 11 v pz abs wabs p v pz 2 2 2 22 wwghp aabs ppp 11 1222 zzppp aaabs 12 zzp aa a p 1 z 2 z 1221 2 1 zzp a v w v pp 22 2 2 这里 称为静压； 称为动压。 为单位体积气体所受有效浮力， 为气体沿 浮力方向升高的距离，乘积 为11断面相对于22 断面单位体积气体的位能，称为位压。 式（442）就是以相对压强计算的气流伯努利方程。 当气流的密度和外界

27、空气的密度相同 ，或两计算点的高 度相同 时，位压为零，式（442）化简为： （444） 式中静压与动压之和称为全压。 当气流的密度远大于外界空气的密度( )，此时相当 于液体总流，式(443)中 可忽略不计，认为各点的当地大气 压相同，式(443)化简为： 21, p p 2 , 2 2 2 2 1 vv g a 12 zz 12 zzg a a 21 zz 21 2 1 v p w v pp 22 2 2 a a （445） 除以 ，即 （446） 由此可见，对于液体总流来说，压强 不论是绝对压强， 还是相对压强，伯努利方程的形式不变。 2.有能量输入或输出有能量输入或输出 总流伯努利方程

28、式(437)是在两过流断面问除水头损失之外， 在无能量输入或输出的条件下导出的。当面过流断面间有水泵、 风机(图411)或水轮机(图412)等流体机械时，存在能量的输入 或输出。 此种情况，根据能量守恒原理，计入单位重量流体经流体机械 获得或失去的机械能， 1221 2 1 zzp v wg v pp 22 2 2 g g v z p 2 2 1 1 1 w p h g v z 2 2 2 2 2 21, pp 式(429)便扩展为有能量输入或输出的伯努利方程式： （447） 式中：+H表示单位重量流体通过流体机械(如水泵)获得的机械 能，对于水泵称为水泵的扬程； -H 表示单位重量流体给流体

29、机械(如水轮机)的机械 能，又称为水轮机的设计水头。 H p z g v 2 1 1 2 11 w p h v z 2 2 2 2 2 图412有能量输出的总流 图411有能量输入的总流 3.两断面间有分流或汇流两断面间有分流或汇流 总流的伯努利方程式(436)，是在 两过流断面间无分流和汇流的条件下导 出的。而实际的供水供气管道沿程多有 分流和汇流这种情况式(436)是否还 能用呢?对于两断面间有分流的流动(图 413)，设想11断面的来流，分为两 股(以虚线划分)分别通过22、33断面。 对 (11断面中的一部分)和22断 面列伯努利方程，其间无分流： （448） 图413沿程分流 1 1

30、 g v g p z 2 1 1 2 1 2 1 2 22 2 2 w h g v g p z 因所取11断面为渐变流断面。面上各点的势能相等，则： （449） 如11断面流速分布较为均匀，则： （450） 故 （451） 近似成立。同理可得： （452） 由以上分析，对于实际I程中沿程分流的总流，当所取过流断面为渐 变流断面，断面上流速分布较为均匀，并计人相应断面之间的水头 损失。 g P Z g P Z 1 11 g v g P Z g v g P Z 22 2 1 1 1 2 11 1 21 2 2 2 1 2 1 22 21 w P h g v g Z g v g p Z 31 2 3

31、 3 3 2 1 1 1 22 wh g v g P Z g v g P Z g v g v 22 2 1 2 1 第四节第四节 总流的动量方程总流的动量方程 总流的动量方程是继连续性方程式、伯努利方程式(436)之后 的第三个积分形式基本方程，它们在流体力学及水力学中习惯地 被称为三大方程，下面由动量原理，推导总流的动量方程。 一、总流的动量方程一、总流的动量方程 设恒定总流，取过流断面、为渐变流断面，面积 为以过流断面及总流的例表面围成的空间为控制体(图314)。控 制体内的流体，经dt时间，由运动到位置。 在流过控制体的总流内，任取元流12,断面面积dA1,dA2，点 流速为 ，dt时间

32、，元流动量的增量 （453） （454） 21, uu 2121 KKKd dtt KK 22212111 KK 1122 KKKd 2221KK 2111 KK dt时间，总流动量的增量，因为过流断面为渐变流断面，各点的 流速平行，按平行矢量和的法则，定义为方向的基 本单位向量，为方向的基本单位向量 （455） 对于不可压缩液体，并引入校正系数，以断面 平均流速v代替点流速 积分得： （456） 式中 是为校正以断面平均速度计算的动量与实际动量的差 异而引入的校正系数，称为流速分布不均匀动量校正系数： （457） 2i2u 1i1u Kd 2 2222 2 iudtdAu A 1 1111

33、1 iudtdAu A Kd 21 2 222 vAvdt 1 111 vAvdt 1122 vvdtQ dtF Av dAu A 2 2 值取决于过流断面上的速度分布，速度分布较均匀的流动， 1.021. 05，通常取 1.0 由动量原理，质点系动员的增量等于作用于该质点系上的外 力的冲量： (458) 投影式： (459) 式(458)、式(459)就是恒定总流的动量方程。方程表 明，作用于控制体内流体上的外力，等于单位时间控制体流出动 量与流人动量之差。综合推导式(447)规定的条件，总流动量方 程的应用条件有：恒定流；过流断面为渐变流断面，不可压缩流 体。 dtF 1122 vvdtQ

34、 F 1122 vvQ zzz yyy xxx vvQF vvQF vvQF 1122 1122 1122 二、动量方程应用举例二、动量方程应用举例 【例例39】 水平放置在混凝土支座上的变直径弯管，弯管两端 与等直径管相连接处的断面11上压力表读数p1=17.6104Pa ， 管中流量qv=0.1m3/s，若直径d1=300，d2=200，转角=60， 如图414所示。求水对弯管作用力F的大小。 【解解】 水流经弯管，动量 发生变化，必然产生作用力F。而 F与管壁对水的反作用力R平衡。 管道水平放置在xoy面上，将R分 解成Rx和Ry两个分力。 取管道进、出两个截面和管内壁 为控制面，如图所

35、示，坐标按图示方向设置。 图414 .根据连续性方程可求得： .列管道进出口的伯努利方程 ，则： .所取控制体受力分析，进、出口控制面上得总压力： sm d q v v /42. 1 3 . 0 41 . 0 4 2 2 1 1 sm d q v v /18. 3 2 . 0 41 . 0 4 2 2 2 2 g v g p g v g p 22 2 22 2 11 2 2 2 2 112 vvpp218. 342. 11000106 .17 223 Pa 3 102 .17 43.123 .0 4 106 .17 23 111 ApP 40.52 .0 4 106 .17 23 222 Ap

36、P （kN） （kN） 壁面对控制体内水的反力Rx、Ry，其方向先假定如图(414)所示。 .写出动量方程 选定坐标系后，凡是作用力（包括其分力）与坐标轴方向一致 的，在方程中取正值；反之，为负值。 沿x轴方向 沿y轴方向 coscos 1 221 vvqRPP Vx coscos 122 1 PPvvqR Vx 568. 060cos43.1240. 560cos42. 118. 31 . 0 (KN) sin0sin 11 vqRP Vx sinsin 11 vqPR Vy 88.1060sin42. 11 . 060sin43.12 (KN) 管壁对水的反作用力 水流对弯管的作用力F与R

37、大小相等，方向相反。总流动量方程 是动量原理的总流表达式，方程给出了总流动量变化与作用力之 间的关系。根据这一特点，求总流与边界面之间的相互作用力问 题，以及因水头损失难以确定运用伯努利方程受到限制的问题， 适于用动量方程求解。 三、动量矩方程三、动量矩方程 上面对动量定理的推导过程中所用之方法、步骤，对动量矩定 理也完全适用，而所得结果与动量定理完全相似，只要在以上的 相应式个，将动量换成动量短就成为动量矩定理；这里不作重复 的推演。 恒定流动的动量矩定理为： 89.1088.10568. 0 2 2 22 YX RRR (KN) 上式表明，在流出面上的流出动量 矩与流入面上的流入动量矩之差

38、等于 外力矩之和。 常见的流体机械中，离心式水泵、 风机都是将其机械能转换为流体的动 能和压能的。水轮机则是利用流体的 动能使叶片机械转动向外输出功率， 其工作原理都是相同的。 图415表示水轮机叶轮的两个叶 片所形成的槽道，流体自叶轮外径 的圆周面流入槽道， 经叶轮内 径的 圆周面流出槽道，进入叶轮 中心区域的导管沿轴向流出；叶轮叶 片就是在流体流动时获得力矩而转动 向外作功的。 iin A n A FrdAVVrdAVVr INou (460) 图415 1r 2r 假定叶片数目足够多，则叶片间的槽道可近似为一元流动，各 截面上的速度是均匀的。还假定叶轮作等角速 的旋转，则叶 轮个流场虽为

39、不定常，但叶轮中的总体动量矩不随时间变化，可 适用定常的动量矩公式，下面我们来导出水轮机(也称涡轮机)的动 量矩公式。 先选取控制面：半径的 进口圆周团和半径 的出口圆周团 之间的流体表面，其中包括各叶片与流体的接触面； 现在分析控制面上的运动情况及受力情况。设流体以相对速度 经半径 的圆周团流入叶片槽道，由于半径 的圆周速度即牵 连速度 ，则流体流入槽道的绝对速度为 (461) 设绝对速度为 与圆周切向夹角为 则其径向分量 和周 向分量 的大小分别为： (462) 1r 2 r rv 1 r 1 r 1 1rVe 111erVVV 111sinVVn 111cosVVt 1V1 1Vn 1V

40、t (463) 同理，流体在流出半径 圆周面上的相对速度 ，牵连速 度 ，则绝对速度为 (464) 设绝对速度为 与圆周切向夹角为 ，则其径向分量 和周 向分量 的大小分别为： (465) (466) 在流量为Q的情况下，流出控制面的动量矩为其切向动量 与半径 的乘积，即： (467) 同理，流入控制团的动量矩为其切向动量与半径之乘积，即 ： (468) 假定无粘性力作用，则控制面中的两圆周面上的压力合力不产 生力矩，只有叶片对流体的作用力矩。 2r 2Vr 22rVe 222erVVV 2V2 2Vn 2Vt 222sinVVn 222cosVVt 2tQV 2r 22222cosrQVrQ

41、Vt 11111cosrQVrQVt 则根据动量矩定理，(464)式减(465)式等于外力矩： (469) 根据作用反作用原理，叶片上获得流体所给的作用力矩力 （470） 这就是欧拉涡轮方程式，是涡轮机械的基本方程式。叶轮所获得 的功率为 （471） 当流出叶片槽道的绝对速度 的方向取半径方向，即 时， 则叶轮获得的力矩公式变为 （472） 相应地，叶轮所获得的功率公式为 （473） 11221112220coscosrVtrVtQrVrVQM 22112221110coscosrVtrVtQrVrVQM 22112211110coscosVeVtVeVtQVeVVeVQMp 2V 90 11

42、1110cosrQVrQVMt 11111coseteVQVVQVP 第五节第五节 理想流体的无旋流动理想流体的无旋流动 在第三章中，在微团运动分析的基础上，见流体的运动分为有旋流动和无旋 流动。理论研究证明只有不可压缩理想流体，运动初始无旋。严格地说，粘性流 体的运动都是有旋流动，但在实际流动中，多有粘性的影响很小，从静止转入流 动（初始无旋）的情况，诸如通风车间，用吸风装置抽气，工作区内形成的气流; 水库中的静水，因闸门开启形成的闸孔出流或堰流；以及空气或水绕物体流动 时，在边界层外面，广阔区域的流动等，都可视为无旋流动。 一、势函数一、势函数： 根据曲线积分定理，无旋流的条件式（550）

43、是表达式 成为某一函数的全微分的必要和充分条件 （474） （475） 得： ， ， （ 476） zzyyxxdududu dzudyudxud zyx dzdydxd zyx xx u yy u zz u gradu 函数 仿照应力场势函数，静电场势函数的定义，称为 速度势函数。由此得出，无旋流是有速势的流动，简称势流; 反之，有速势的流动是无旋流，两者含义相同。 将式（4-37）不可压缩流体连续性微分方程 ： （478） 即： （479） 式中 拉普拉斯算子 式（478）是著名的拉普拉斯方程，满足拉普拉斯方程的函 数是调和函数。所以，调和函数的一切性质，也是速度势函数 拥有的性质。 ),

44、(zyx x u x y u y z u z 0 2 2 2 2 2 2 zyx 0 2 2 2 2 2 2 2 2 zyx 由以上分析可知，不可压缩流体无旋流动的问题，归结为在 给定的边界条件下，求解拉普拉斯方程，一旦求得速度势 ， 就可由式（476）求得流速 ，解得压强，问题得 到解决。 二、流函数二、流函数 对于平面运动，有连续性微分方程 ，移项得 根据曲线积分定理，前式是表达式 成为某一函数 的全微分的必要和充分条件 （480） 比较 （481） 得 （482） ),(zyxuuuu x u x y u y 0 x u x y u y dxudyu yx yx, dxudyud yx

45、dydxd yx xx u yy u 函数 称为流函数。由流函数的引出条件可知，凡是不可 压缩流体的平面的流动，连续性微分方程成立，不论无旋流动 或有旋流动，都存在流函数，而只有无旋流动才有流速势，可 见流函数比流速势更具有普遍意义。 1.流函数具有以下性质： 流函数的等值线是流线 证明： 流函数值相等 ，由式得流函数等值线方 程 则 上式即平面流动的流线方程，故 流函数的等值线是流线，给流线以不同值，便得到族流线给流 函数以不同值，便得到流线族。 .两条流线的流函数的差值，等于通过该两流线间的单 宽流量： yx, 0,dc 0dxudyu yx yx u dy u dx （483） 这一性质

46、也可表述为：平面流动中，通过任一曲线的单宽 流量，等于该曲线两端流函数的差值。 .平面无旋流动的等流函数线（流线）与等势线正交。 证明：对于平面无旋流动，同时存 在流速势函数和流函数，由等流函数线方程 某一点的斜率 由等势线方程 dlynuxnudludq yxn ,cos,cos dluu dl dx ydl dy x dxudyu yx d 12 2 1 2 1 ddqq q q 0dxudyud yx x y u u dx dy m 1 0dyudxud yx 图416流函数 同一点等势线斜率 （484） 等流函数线与等势线正交，故等势线也就是过流断面线。 .平面无旋流动，流函数是调和函

47、数。 证明：因为平面无旋流动 则 得 带入上式，得 （485） x y u u dx dy m 2 1 21 y x x y u u u u mm 0 2 1 y u x u x y z x u y u yx , 0 y u x u x y 0 2 2 2 2 yx 即： 平面无旋流动的流函数满足拉普拉斯方程，是调和函数。 式中 拉朴拉斯算子 （486） 式即柯西黎曼条件。满足拉普拉斯方程和柯西黎曼条 件，是一对共轭调和函数。 三、几种常见的基本平面势流：三、几种常见的基本平面势流： 拉普拉斯方程在复杂的边界条件下，虽然难以求解，一 些简单的平面势流，其流速势和流函数却不难求得。研究这些 简单

48、的平面势流的意义在于通过简单势流的叠加，往往能组合 成符合某些给定边界条件的复杂流场。 0 2 2 2 2 2 2 yx xy yx 1均匀直线流均匀直线流 均匀直线流是流场中各点速度大小相等，方向相同的流动，是 一种最简单的平面势。速度场 ， ; 速度势 （487） （488） 若均匀直线流流速平行于轴 （489） 若均匀直线流流速平行于轴 （490） au x bu y dyudxu yx byax dxudyu yx bxay ayaxuy, 0 bxbyux, 0 图417均匀直线流 2源流源流 如图418所示，在平面势流中，源流就是流体从潭点均匀地向各个 方向出流的流动。组成这种流型的线，就是源点所在平面势流中i面 上，从源点0出发的一族射线。 速度场 速度 流函数 等势线方程 等势线是以o点为圆心的同心圆。 流线方程 ，流线是由o点引出的射线以直角坐标 系表示。 r q r u 2 0 u rd udrur rdr q r q ln 22 drurdur 22 q r qr d crc, 图418平面源流 cc, （491） （492） 3汇流汇流 流体从四周沿经向均匀的流入的流动称之为汇流。流入汇 点单位厚度流量称为汇流强度-q。汇流的速度势和流函数的表 达式与源流的相同，符号相反