1、2022-7-221第三章第三章 流体动力学基础流体动力学基础11 描述流体运动的两种方法描述流体运动的两种方法 16 伯努利（伯努利（Bernoulli）方程的应用）方程的应用 18 液体的空化和空蚀现象液体的空化和空蚀现象17 定常流动的动量方程和动量矩方程定常流动的动量方程和动量矩方程12 流体运动的一些基本概念流体运动的一些基本概念14 理想流体的运动微分方程理想流体的运动微分方程13 流体运动的连续性方程流体运动的连续性方程15 理想流体微元流束的伯努力方程理想流体微元流束的伯努力方程2022-7-222 流体运动学研究流体的运动规律，如速度、加速度等运动参数的变化规律，而流体动力学
2、则研究流体在外力作用下的运动规律，即流体的运动参数与所受力之间的关系。本章主要介绍流体运动学和流体动力学的基本知识，推导出流体动力学中的几个重要的基本方程：连续性方程、动量方程和能量方程，这些方程是分析流体流动问题的基础。2022-7-223第一节 描述流体运动的两种方法 连续介质模型的引入，使我们可以把流体看作为由无数个流体质点所组成的连续介质，并且无间隙地充满它所占据的空间。我们把流体质点运动的全部空间称为流场。由于流体是连续介质，所以描述流体运动的各物理量(如速度、加速度等)均应是空间点的坐标和时间的连续函数。根据着眼点的不同，流体力学中研究流体的运动有两种不同的方法，一种是拉格朗日（L
3、agrange）方法，另一种是欧拉（Euler）方法。拉格朗日方法又称随体法，是从分析流场中个别流体质点着手来研究整个流体运动的。这种研究方法，最基本2022-7-224 的参数是流体质点的位移，在某一时刻，任一流体质点的位置可表示为：X=x(a，b，c，t)y=y(a，b，c，t)z=z(a，b，c，t)(3-1)式中a、b、c为初始时刻任意流体质点的坐标，即不同的a、b、c代表不同的流体质点。对于某个确定的流体质点，a、b、c为常数，而t为变量，则得到流体质点的运动规律。对于某个确定的时刻，t为常数，而a、b、c为变量，得到某一时刻不同流体质点的位置分布。通常称a、b、c为拉格朗日变量，它
4、不是空间坐标的函数，而是流体质点标号。2022-7-225 将式（3-1）对时间求一阶和二阶导数，可得任意流体质点的速度和加速度为：(3-2)(3-3),(tcbautxu),(tcbavtyv),(tcbawtzw),(22tcbaatxtuaxx),(22tcbaatytvayy),(22tcbaatztwazz2022-7-226 同样，流体的密度、压强和温度也可写成a、b、c、的函数，即=（a，b，c，），P=P(a，b，c，)，t=t(a，b，c，)。欧拉法，又称局部法，是从分析流场中每一个空间点上的流体质点的运动着手，来研究整个流体的运动的，即研究流体质点在通过某一空间点时流动参数
5、随时间的变化规律。所以流体质点的流动是空间点坐标（x，y，z）和时间t的函数，例如：流体质点的三个速度分量、压强和密度可表示为：u=u(x，y，z，t)v=v(x，y，z，t)（3-4）w=w(x，y，z，t)式中，u，v，w分别表示速度矢量在三个坐标轴上的分量：kwj vi uV2022-7-227 P=p(x，y，z，t)=（x，y，z，t)(3-5)式（3-4）中，当参数x，y，z不变而改变时间t，则表示空间某固定点的速度随时间的变化规律。当参数t不变，而改变x，y，z，则代表某一时刻，空间各点的速度分布。x，y，z有双重意义，一方面它代表流场的空间坐标，另一方面它代表流体质点在空间的位
6、移。根据流体连续介质假设，每一个空间点上都有流体质点所占据。而占据每一个空间点上的流体质点都有自己的速度，有速度必然产生位移。也就是说，空间坐标x，y，z也是流体质点位移的变量，它也是时间t的函数：x=x(t)y=y(t)z=z(t)(3-6)2022-7-228 式（3-6）是流体质点的运动轨迹方程，将上式对时间求导就可得流体质点沿运动轨迹的三个速度分量 （3-7）现在用欧拉法求流体质点的加速度。由于加速度定义为在dt时刻内，流体质点流经某空间点附近运动轨迹上一段微小距离时的速度变化率，于是可按复合函数的求导法则，分别将式（3-4）中三个速度分量对时间取全导数，并将式（3-7）代入，即可得流
7、体质点在某一时刻经过某空间点时的三个加速度分量txuddtvddytwddz2022-7-229 (3-8)用矢量 表示加速度，即 。根据矢量分析的点积公式 （3-9）式中 是矢量微分算子。由式（3-8）可知，用欧拉法求得的流体质点的加速度由两部分组成；第一部分是由于某一空间点上的流体质点zwwywvxwutwazvwyvvxvutvazuwyuvxuutuazyxakajaiaazyxVVtVa)(kzjyix2022-7-2210 的速度随时间的变化而产生的，称为当地加速度，即式（3-8）中等式右端的第一项 、；第二部分是某一瞬时由于流体质点的速度随空间点的变化称为迁移加速度，即式（3-8
8、）中等式右端的后三项 、等；当地加速度和迁移加速度之和称为总加速度。为了加深对当地加速度和迁移加速度的理解，现举例说明这两个加速度的物理意义。如图3-1所示，不可压缩流体流过一个中间有收缩形的变截面管道，截面2比截面1小，则截面2的速度就要比截面1的速度大。所以当流体质点从1点流到2点时，由于截面的收缩引起速度的增加，从而产生了迁移加速度，如果在某一段时间内流进管道的流体输入量有变化（增加或减少），则管道中每一点上流体质点的速tutvtwxuuyuvzuw2022-7-2211图 3-1 中间有收缩形的变截面管道内的流动2022-7-2212 度将相应发生变化（增大或减少），从而产生了当地加速
9、度。应该注意，流体质点和空间点是两个截然不同的概念，空间点指固定在流场中的一些点，流体质点不断流过空间点，空间点上的速度指流体质点正好流过此空间点时的速度。用欧拉法求流体质点其他物理量的时间变化率也可以采用式(3-9)的形式，即 （3-10）式中，括弧内可以代表描述流体运动的任一物理量，如密度、温度、压强，可以是标量，也可以是矢量。称为全导数，称为当地导数，称为迁移导数。)()()(D)D(VtttD)D(t)()()(V2022-7-2213 由上述可知，采用欧拉法描述流体的流动，常常比采用拉格朗日法优越，其原因有三。一是利用欧拉法得到的是场，便于采用场论这一数学工具来研究。二是采用欧拉法，
10、加速度是一阶导数，而拉格朗日法，加速度是二阶导数，所得的运动微分方程分别是一阶偏微分方程和二阶偏微分方程，在数学上一阶偏微分方程比二阶偏微分方程求解容易。三是在工程实际中，并不关心每一质点的来龙去脉。基于上述三点原因，欧拉法在流体力学研究中广泛被采用。当然拉格朗日法在研究爆炸现象以及计算流体力学的某些问题中还是方便的。2022-7-2214 【例例3-1】已知用拉格朗日变量表示得速度分布为 u=(a+2)et-2，v=(b+2)et-2，且t=0时，x=a，y=b。求（1）t=3时质点分布；（2）a=2，b=2质点的运动规律；（3）质点加速度。【解解】根据（3-2）式得 将上式积分，得 上式中
11、c1、c2为积分常数，它仍是拉格朗日变量的函数。利用t=0时，x=a，y=b得c1=-2，c2=-22)2(teatx2)2(tebty12)2(cteaxt22)2(ctebyt2022-7-2215 X=(a+2)et-2t-2 y=(b+2)et-2t-2 （1）将t=3代入上式 得 X=(a+2)e3-8 y=(b+2)e3-8 （2）a=2，b=2时 x=4et-2t-2 y=4et-2t-2 (3)teatu)2(tebtv)2(2022-7-2216 【例例3-2】在任意时刻，流体质点的位置是x=5t2，其迹线为双曲线xy=25。质点速度和加速度在x和y方向的分量为多少？【解解】
12、根据式（3-7）得 由式（3-8）得ttttxu10)5(dddd2txxxttvdd12525ddddy23221010)5(125ttt10tuax430ttvay2022-7-2217第二节 流体运动的一些基本概念 在讨论流体运动的基本规律和基本方程之前，为了便于分析、研究问题，先介绍一些有关流体运动的基本概念。一、定常流动和非定常流动一、定常流动和非定常流动 根据流体的流动参数是否随时间而变化，可将流体的流动分为定常流动和非定常流动，现举例说明如下：如图3-2所示装置，将阀门A和B的开度调节到使水箱中的水位保持不变，则水箱和管道中任一点(如1点、2点和3点等)的流体质点的压强和速度都不
13、随时间而变化，但由于1、2、3各点所处的空间位置不同，故其压强和速度值也就各2022-7-2218图 3-2 流体的出流2022-7-2219图 3-2 流体的出流2022-7-2220 不相同。这时从管道中流出的射流形状也不随时间而变。这种运动流体中任一点的流体质点的流动参数(压强和速度等)均不随时间变化，而只随空间点位置不同而变化的流动，称为定常流动。现将阀门A关小，则流入水箱的水量小于从阀门B流出的水量，水箱中的水位就逐渐下降，于是水箱和管道任一点流体质点的压强和速度都逐渐减小，射流的形状也逐渐向下弯曲。这种运动流体中任一点流体质点的流动参数(压强和速度等)随时间而变化的流动，称为非定常
14、流动。由上可见，定常流动的流场中，流体质点的速度、压强和密度等流动参数仅是空间点坐标x、y、z的函数，而与时间t无关，用表示任一流动参数（即可表示u，v，w，p，等），则 =(x，y，z)(3-11)2022-7-2221 由于是定常流动，故其流动参数对时间的偏导数等于零，即 （3-12）因此，定常流动时流体加速度可简化成 （3-13）由式(3-13)可知，在定常流动中只有迁移加速度。例如图3-2中，当水箱的水位保持不变时，2点到3点流体质点的速度减小，而4点到5点速度增加，都是由于截面变化而引起的迁移加速度。若迁移加速度为零，则为均匀流动，例如流体质点在等截面管道中的流动(3点到4点)。在供
15、水和通风系统中，只要泵和风机的转速不变，运转稳定，则水管和风道中的流体流动都是定常流动。又如0tVVa)(2022-7-2222 火电厂中，当锅炉和汽轮机都稳定在某一工况下运行时，主蒸汽管道和给水管道中的流体流动也都是定常流动。可见研究流体的定常流动有很大的实际意义。二、迹线与流线二、迹线与流线 迹线是流场中某一质点运动的轨迹。例如在流动的水面上撒一片木屑，木屑随水流漂流的途径就是某一水点的运动轨迹，也就是迹线。流场中所有的流体质点都有自己的迹线，迹线是流体运动的一种几何表示，可以用它来直观形象地分析流体的运动，清楚地看出质点的运动情况。迹线的研究是属于拉格朗日法的内容，迹线表示同一流体质点在
16、不同时刻所形成的曲线，其数学表达式为：（3-14）twzvyuxdddd2022-7-2223 式(3-14)就是迹线微分方程，是自变量。流线是某一瞬时在流场中所作的一条曲线，在这条曲线上的各流体质点的速度方向都与该曲线相切，因此流线是同一时刻，不同流体质点所组成的曲线，如图3-3所示。流线可以形象地给出流场的流动状态。通过流线，可以清楚地看出某时刻流场中各点的速度方向，由流线的密集程度，也可以判定出速度的大小。流线的引入是欧拉法的研究特点。例如在流动水面上同时撤一大片木屑，这时可看到这些木屑将连成若干条曲线，每一条曲线表示在同一瞬时各水点的流动方向线就是流线。1、流线的基本特性 (1)在定常
17、流动时，因为流场中各流体质点的速度不随2022-7-2224图 3-3 流线的概念2022-7-2225 时间变化，所以通过同一点的流线形状始终保持不变，因此流线和迹线相重合。而在非定常流动时，一般说来流线要随时间变化，故流线和迹线不相重合。(2)通过某一空间点在给定瞬间只能有一条流线，一般情况流线不能相交和分支。否则在同一空间点上流体质点将同时有几个不同的流动方向。只有在流场中速度为零或无穷大的那些点，流线可以相交，这是因为，在这些点上不会出现在同一点上存在不同流动方向的问题。速度为零的点称驻点，速度为无穷大的点称为奇点。(3)流线不能突然折转，是一条光滑的连续曲线。(4)流线密集的地方，表
18、示流场中该处的流速较大，稀疏的地方，表示该处的流速较小。2022-7-2226 2、流线微分方程 现由矢量分析法导出流线微分方程。设在某一空间点上流体质点的速度矢量 ，通过该点流线上的微元线段 。由流线的定义知，空间点上流体质点的速度与流线相切。根据矢量分析，这两个矢量的矢量积应等于零，即 即 上式又可写成kwjviuVkzj yi xLdddd0d d d k j i dzyx v wu LV0dd0dd0ddzuxwywzvxvyu2022-7-2227 （3-15）式（3-15）就是流线的微分方程，式中时间t是个参变量。【例例3-3】有一流场，其流速分布规律为：u=-ky，v=kx，w=
19、0，试求其流线方程。【解解】由于w=0，所以是二维流动，二维流动的流线方程微分为 将两个分速度代入流线微分方程（3-15），得到 即 xdx+ydy=0 积分上式得到 x2+y2=c 即流线簇是以坐标原点为圆心的同心圆。),(d),(d),(dtzyxwztzyxvytzyxuxvyuxddxyyxkdkd2022-7-2228 三、流管、流束和总流三、流管、流束和总流 在流场中任取一条不是流线的封闭曲线，通过曲线上各点作流线，这些流线组成一个管状表面，称之为流管。如图3-4所示。因为流管是由流线构成的，所以它具有流线的一切特性，流体质点不能穿过流管流入或流出(由于流线不能相交)。流管就像固体
20、管子一样，将流体限制在管内流动。过流管横截面上各点作流线，则得到充满流管的一束流线簇，称为流束。当流束的横截面积趋近于零时，则流束达到它的极限流线。在流束中与各流线相垂直的横截面称为有效截面。流线相互平行时，有效截面是平面。流线不平行时，有效截面是曲面，如图3-5所示。有效截面面积为无限小的流束2022-7-2229 和流管，称为微元流束和微元流管。在每一个微元流束的有效截面上，各点的速度可认为是相同的。无数微元流束的总和称为总流。自然界和工程中所遇到的管流或渠流都是总流。根据总流的边界情况，可以把总流流动分为三类：（1）有压流动 总流的全部边界受固体边界的约束，即流体充满流道，如压力水管中的
21、流动。（2）无压流动 总流边界的一部分受固体边界约束，另一部分与气体接触，形成自由液面，如明渠中的流动。（3）射流 总流的全部边界均无固体边界约束，如喷嘴出口的流动。在总流的有效截面上，流体与固体边界接触的长度称为湿周，用符号表示。2022-7-2230图 3-4 流管和流束2022-7-2231图 3-5 有效截面2022-7-2232 总流的有效截面面积与湿周之比称为水力半径，用符号Rh表示，即 关于湿周和水力半径的概念在非圆截面管道和管束的水力计算中常常用到。四、流量和平均流速四、流量和平均流速 单位时间内通过有效截面的流体体积称为体积流量，以qv表示。其单位为m3/s、m3/h等。单位
22、时间内通过有效截面的流体质量称为质量流量,以qm表示，其单位为kg/s、t/h等。由于微元流束有效截面上各点的流速V是相等的，所以通过微元流束有效截面积为的体积流量dqv和质量流量dqm分别为：dqv=VdA (3-16)dqm=VdA (3-17)ARh2022-7-2233 由于流束是由无限多的微元流束组成的，所以通过流束有效截面面积为的流体体积流量和质量流量分别由式(3-16)和式(3-17)积分求得，即 (3-18)(3-19)以上计算必须先找出微元流束的速度V在整个流束有效截面上的分布规律，这在大部分工程问题中是不能用解析法来确定的。在工程计算中为了方便起见，引入平均流速的概念。平均
23、流速是一个假想的流速，即假定在有效截面上各点都以相同的平均流速流过，这时通过该有效截面上的体积流量仍与各点以真实流速流动时所得到的体积流量相同。AVAVqdAmAVqd2022-7-2234 若以表示平均流速，按其定义可得：(3-20)(3-21)五、一维、二维和三维流动五、一维、二维和三维流动 一般的流动都是在三维空间的流动，流动参数是x、y、z三个坐标的函数，在流体力学中又称这种流动为三维流动。当我们适当地选择坐标或将流动作某些简化，使其流动参数在某些情况下，仅是x、y两个坐标的函数，称这种流动为二维流动。是一个坐标的函数的流动，称为一维流动。如图3-6所示的带锥度的圆管内黏性流体的流动，
24、流体质点运动参数，如速度，即是半径r的函数，又是沿轴AVAVAVqAAVddAqVV2022-7-2235图 3-6 管内流动速度分布2022-7-2236 线距离的函数，即：u=u(r，x)。显然这是二元流动问题。工程上在讨论其速度分布时，常采用其每个截面的平均值u。就将流动参数如速度，简化为仅与一个坐标有关的流动问题，这种流动就叫一维流动，即：u=u(x)。如图3-7所示的绕无限翼展的流动就是二维流动，二维流动的参数以速度为例，可写成：如图3-8所示的绕有限宽翼展的流动就是三维流动，三维流动的参数以速度为例，可写成：六、均匀流和非均匀流六、均匀流和非均匀流 根据流场中同一条流线各空间点上的
25、流速是否相同，可将总流分为均匀流和非均匀流。若相同则称为均匀流，j ),(i ),(xxvyxuVkzyxwzyxvzyxuV),(j ),(i ),(2022-7-2237图 3-7 绕无限翼展的流动2022-7-2238图 3-8绕有限翼展的流动2022-7-2239 否则称为非均匀流。由此定义可知在均匀流中，流线是彼此平行的直线，过水断面（有效截面）是平面。如在等直径的直管道内的水流都是均匀流(图3-9)。注意在均匀流中各流线上的流速大小不定彼此相等在非均匀流中，流线或者是不平行的直线，或者是曲线，如图3-10所示。一般非均匀流的过水断面（有效截面）是曲面。非均匀流按流速的大小和方向沿流
26、线变化的缓、急程度又可分为缓（渐）变流和急变流两种（图3-11）。流速的大小和方向沿流线逐渐改变的非均匀流，称为缓（渐）变流。显然，缓（渐）变流的流线的曲率半径r较大，流线之间的夹角较小。因此，缓（渐）变流是一种流线几乎平行又近似直线的流动，其极限情况就是均匀流。缓（渐）变流的有效截面可看作平面，但是缓（渐）变流各个过水断面的形状和大小是沿程逐渐改变的，各个过水断面上的流速分布图形也是沿程逐渐改变的。流速的大小和2022-7-2240方向沿流线急剧变化的非均匀流，称为急变流。显然其流线之间的夹角较大，或者流线曲率半径较小，或者两者兼而有之。2022-7-2241图 3-9 均匀流2022-7-
27、2242图 3-10 非均匀流2022-7-2243急变流缓变流缓变流缓变流缓变流缓变流急变流急变流急变流急变流图 3-11 缓变流和急变流2022-7-2244第三节 流体流动的连续性方程 连续性方程是质量守恒定律在流体力学中的应用。我们认为流体是连续介质，它在流动时连续地充满整个流场。在这个前提下，当研究流体经过流场中某一任意指定的空间封闭曲面时，可以断定：若在某一定时间内，流出的流体质量和流入的流体质量不相等时，则这封闭曲面内一定会有流体密度的变化，以便使流体仍然充满整个封闭曲面内的空间；如果流体是不可压缩的，则流出的流体质量必然等于流入的流体质量。上述结论可以用数学分析表达成微分方程，
28、称为连续性方程。2022-7-2245 一、直角坐标系下连续性微分方程式一、直角坐标系下连续性微分方程式 设在流场中任取一个微元平行六面体，其边长分别为dx、dy和dz，如图3-12所示。假设微元平行六面体形心的坐标为x、y、z，在某一瞬时t经过形心的流体质点沿各坐标轴的速度分量为u、v、w，流体的密度为。现讨论流体经六面体各面的流动情况。先分析x轴方向，由式(3-4)和式(3-5)可知，u和都是坐标和时间的连续函数，即u=u(x，y，z，t)和=(x，y，z，t)。根据泰勒级数展开式，略去高于一阶的无穷小量，得在d时间内，沿轴方向从左边微元面积dydz流入的流体质量为2022-7-2246图
29、 3-12 流场中的微元平行六面体2022-7-2247 同理可得在dt时间内从右边微元面积dydz流出的流体质量为 (3-22)上述两者之差为在dt时间内沿x轴方向流体质量的变化，即 (3-23)tzytzyxxutzyxxddd,2d,2dtzyxtuuxttzyxtutzyxuxttzyxddd2d2dddd2d),(2d),(tzyxtuuxtddd2d2dtzyxuxtzyxxuxxudddd)(ddddd2022-7-2248 同理可得，在dt时间内沿y轴和z轴方向流体质量的变化分别为：因此，在dt时间内经过微元六面体的流体质量总变化为 (3-24)由于流体是作为连续介质来研究的，
30、所以式(3-24)所表示的六面体内流体质量的总变化，唯一的可能是因为六面体内流体密度的变化而引起的。因此式(3-24)应和由于流体密度的变化而产生的六面体内的流体质量变化相等。设开始瞬时流体的密度为，经过dt时间后的密度为tzyxvydddd)(tzyxwzdddd)(tzyxzwyvxuddddttttzyxd)d,(2022-7-2249 则可求出在dt时间内，六面体内因密度的变化而引起的质量变化为 (3-25)根据连续性条件，式(3-24)和式(3-25)应相等，经简化得到 (3-26)式（3-26）为可压缩流体非定常三维流动的连续性方程。若流体是定常流动，则 ，上式成为 (3-27)式
31、（3-27）为可压缩流体定常三维流动的连续性方程。若流体是不可压缩的，不论是定常或非定常流动均tzyxtzyxzyxttddddddddddd0zwyvxut0t0zwyvxu2022-7-2250 为常数，故式(3-27)成为 (3-28)式（3-28）为不可压缩流体三维流动的连续性的方程。它的物理意义是：在同一时间内通过流场中任一封闭表面的体积流量等于零，也就是说，在同一时间内流入的体积流量与流出的体积流量相等。在流体力学中时常讨论所谓平面（二维）流动，即平行任何一个坐标平面的流动。若这种流动的流动参数（如速度、压强）只沿x、y两个坐标轴方向发生变化，则式（3-28）可以写成 (3-29)
32、由于在推导上述连续性方程时，没有涉及作用力的问题，所以不论是对理想流体还是实际流体都是适用的。0zwyvxu0yvxu2022-7-2251 二、微元流束和总流的连续性方程二、微元流束和总流的连续性方程 在工程上和自然界中，流体流动多数都是在某些周界所限定的空间内沿某一方向流动，即一维流动的问题，所谓一维流动是指流动参数仅在一个方向上有显著的变化，而在其它两个方向上的变化非常微小，可忽略不计。例如在管道中流动的流体就符合这个条件。在流场中取一微元流束(图3-13)。假定流体的运动是连续的、定常的，则微元流管的形状不随时间而改变。又根据流管的特性，流体质点不能穿过流管表面，因此在单位时间内通过微
33、元流管的任一有效截面的流体质量都应相等，即 1V1dA1=2V2dA2=VdA=常数 （3-30）式中 dA1、dA2分别为1、2两个有效截面的面积，m2；2022-7-2252图 3-13 流场中的微元流束2022-7-2253 V1 、V2分别为dA1和dA2上的流速，也称为真实流速，m/s；1 、2分别为和处的流体密度，kg/m3。对于由无限多微元流束所组成的总流（例如流体在管道中的流动），可对式（3-30）进行积分得 （3-31）式中 A1 和A2分别为总流1和2两个有效截面的面积，m2。式（3-31）为一维流动积分形式总流的连续性方程。设 和 是总流两个有效截面l和2上的平均流速，则
34、式（3-31）可写成 (3-32)常数AAAAVAVAVddd212221111V2V222111AVAV2022-7-2254 式中1和2分别代表截面和上的平均密度，kg/m3。式（3-32）表示当流动为可压缩流体定常流体动时，沿流动方向的质量流量为一个常数。对不可压缩均质流体常数，则式（3-32）成为 (3-33)式（3-33）为不可压缩流体一维定常流动的总流连续性方程。该式说明一维总流在定常流动条件下，沿流动方向的体积流量为一个常数，平均流速与有效截面面积成反比，即有效截面面积大的地方平均流速小，有效截面面积小的地方平均流速就大。2211AVAV2022-7-2255 【例例3-4】假设
35、有一不可压缩流体三维流动，其速度分布规律为）U=3（x+y3)，v=4y+z2，w=x+y+2z。试分析该流动是否连续。【解解】根据式（3-28）所以 故此流动不连续。不满足连续性方程的流动是不存在的 3xu4yv2zw09 zwyvxu2022-7-2256 【例例3-5】有一不可压缩流体平面流动，其速度分布规律为u=x2siny，v=2xcosy，试分析该流动是否连续。【解解】根据式（3-29）所以 故此流动是连续的。yxxusin2yxyvsin20)sin2(sin2yxyxyvxu2022-7-2257 【例例3-6】有一输水管道，如图3-14所示。水自截面1-1流向截面2-2。测得
36、截面1-1的水流平均流速 m/s，已知d1=0.5m，d2=1m，试求截面2-2处的平均流速 为多少？【解解】由式（3-33）得 (m/s)2V2V22221144dVdV5.015.02222112ddVV2022-7-2258图 3-14 输水管道2022-7-2259第四节 理想流体的运动微分方程 在流动的理想流体中，取出一个微元平行六面体的微团，它的各边长度分别为dx、dy和dz，如图3-15所示。由于是理想流体，没有黏性，运动时不产生内摩擦力，所以作用在流体微团上的外力只有质量力和压强。该压强与静压强一样，垂直向内，作用在流体微团的表面上。假设六面体形心的坐标为x、y、z，压强为p。
37、先分析x方向的运动，在垂直于x轴的左右两个平面中心点上的压强各等于 由于是微元面积，所以这些压强可以作为各表面上的2dxxpp2dxxpp2022-7-2260图 3-15 推导欧拉运动微分方程用图2022-7-2261 平均压强。设在六面体形心上的单位质量的质量力分量为fx、fy和fz，则作用在微元平行六面体的流体微团上的质量力在轴方向的分量为 fxdxdydz 又流体微团的加速度在x轴上的投影为 ，则根据牛顿第二定律得x轴方向的运动微分方程 将上式各项除以流体微团的流体质量dxdydz，化简后得：同理 (3-34)DtDuDtDuzyxzyxxppzyxxppzyxfxddddd2ddd2
38、ddddDtDuxpfx1DtDvypfy1DtDwzpfz12022-7-2262 这就是理想流体的运动微分方程，早在1755年就为。对于静止的流体u=v=w=0，则由式（3-34）可以直接得出流体平衡微分方程，即欧拉平衡微分方程式（2-3）。因此欧拉平衡微分方程只是欧拉运动微分方程的一个特例。如果把加速度写成展开式，可将欧拉运动微分方程写成如下形式 (3-35)zwwywvxwutwzpfzvwyvvxvutvypfzuwyuvxuutuxpfzyx1112022-7-2263 在一般情况下，作用在流体上的质量力fx、fy和fz 是已知的，对理想不可压缩流体其密度为一常数。在这种情况下，式
39、（3-35）中有四个未知数u、v、w和p，而式（3-35）中有三个方程，再加上不可压缩流体的连续性方程（3-28），就从理论上提供了求解这四个未知数的可能性。2022-7-2264 第五节 理想流体微元流束的伯努利方程 一、理想流体微元流束的伯努利方程一、理想流体微元流束的伯努利方程 理想流体的运动微分方程（3-35）只有在少数特殊情况下才能求解。在下列几个假定条件下：(1)不可压缩理想流体的定常流动；(2)沿同一微元流束（也就是沿流线）积分；(3)质量力只有重力。即可求得理想流体微元流束的伯努利方程。假定流体是定常流动，则有 ，0t0zwyvxu2022-7-2265 因此式(3-35)可写
40、成 (3-36)假如流体微团沿流线的微小位移ds在三个坐标轴上的投影为dx、dy和dz。现用dx、dy和dz分别乘以式（3-36）的第一式、第二式和第三式，则可得到zwwywvxwuzpfzvwyvvxvuypfzuwyuvxuuxpfzyx1112022-7-2266 (3-37)由流线微分方程（3-15）有 udy=vdx ydz=wdy (3-38)wdx=udz 将式（3-38）代入式（3-37）中的对应项，则得zzwwzywvzxwuzzpzfyzvwyyvvyxvuyypyfxzuwxyuvxxuuxxpxfzyxdddd1ddddd1ddddd1d2022-7-2267 (3-3
41、9)将式（3-39）的三个方程相加，得到 (3-40)由于式（3-40）中的dx、dy和dz是流体微团沿流线微小位移ds的三个分量，所以要沿流线（或微元流束）进行积分。wwzzwwyywwxxwwzzpzfvvzzvvyyvvxxvvyypyfuuzzuuyyuuxxuuxxpxfzyxddddd1dddddd1dddddd1dwwvvuuzzpyypxxpzfyfxfzyxdddddd1)ddd(2022-7-2268 式（3-40)中的 假设质量力只有重力，fx=0，fy=0，fz=-g，即z轴垂直向上，oxy为水平面。则式(3-40)可写成 又假设为不可压缩均质流体，即=常数，积分后得
42、或 (3-41)式（3-41）称为理想流体微元流束的伯努利方程。方程右边的常数对不同的流线有不同的值。该方程的适用范围pzzpyypxxpdddd2222d21)(d21dddVwvuwwvvuu0d21d1d2Vpzg常数22Vpgz常数gVgpz222022-7-2269 是：理想不可压缩均质流体在重力作用下作定常流动，并沿同一流线（或微元流束）。若1、2为同一条流线（或微元流束）上的任意两点，则式（3-41）也可写成 (3-42)在特殊情况下，绝对静止流体V=0，由式(3-41)可以得到静力学基本方程 二、方程的物理意义和几何意义二、方程的物理意义和几何意义 为了进一步理解理想流体微元流
43、束的伯努利方程，现来叙述该方程的物理意义和几何意义。1、物理意义 理想流体微元流束的伯努利方程式（3-41）中，左端gVgpzgVgpz2222222111常数gpz2022-7-2270 前两项的物理意义，在静力学中已有阐述，即第一项z表示单位重量流体所具有的位势能；第二项p/(g)表示单位重量流体的压强势能；第三项V2/(2g)理解如下：由物理学可知，质量为m的物体以速度V运动时，所具有的动能为Mv2/2，则单位重量流体所具有的动能为V2/(2g)即(mV2/2)/(mg)=V2/(2g)。所以该项的物理意义为单位重量流体具有的动能。位势能、压强势能和动能之和称为机械能。因此，伯努利方程可
44、叙述为：理想不可压缩流体在重力作用下作定常流动时，沿同一流线（或微元流束）上各点的单位重量流体所具有的位势能、压强势能和动能之和保持不变，即机械能是一常数，但位势能、压强势能和动能三种能量之间可以相互转换，所以伯努利方程是能量守恒定律在流体力学中的一种特殊表现形式。2022-7-2271 2、几何意义图 理想流体微元流束的伯努利方程式（3-41）中，左端前两项的几何意义，同样在静力学中已有阐述，即第一项z表示单位重量流体的位置水头，第二项p/(g)表示单位重量流体的压强水头，第三项V2/(2g)与前两项一样也具有长度的量纲。它表示所研究流体由于具有速度V，在无阻力的情况下，单位重量流体所能垂直
45、上升的最大高度，称之为速度水头。位置水头、压强水头和速度水头之和称为总水头。由于它们都表示某一高度，所以可用几何图形表示它们之间的关系，如图3-16所示。因此伯努利方程也可叙述为：理想不可压缩流体在重力作用下作定常流动时，沿同一流线(或微元流束)上各点的单位重量流体所具有的位置水头、压强水头和速度水头之和保持不变，即总水头是一常数。2022-7-2272图 3-16 总水头线和静水头线2022-7-2273第六节 伯努利（Bernoulli）方程的应用 理想流体微元流束的伯努利方程，在工程中广泛应用于管道中流体的流速、流量的测量和计算，下面以应用最广泛的皮托管和文特里流量计为例，介绍它们的测量
46、原理和伯努利方程的应用。一、皮托管一、皮托管 在工程实际中，常常需要来测量某管道中流体流速的大小，然后求出管道的平均流速，从而得到管道中的流量，要测量管道中流体的速度，可采用皮托管来进行，其测量原理如图3-17所示。在液体管道的某一截面处装有一个测压管和一根两端2022-7-2274VBAZZ图 3-17 皮托管测速原理2022-7-2275 开口弯成直角的玻璃管（称为测速管）。将测速管（又称皮托管）的一端正对着来流方向，另一端垂直向上，这时测速管中上升的液柱比测压管内的液柱高h。这是由于当液流流到测速管入口前的A点处，液流受到阻挡，流速变为零，则在测速管入口形成一个驻点A。驻点A的压强PA称
47、为全压，在入口前同一水平流线未受扰动处（例如B点）的液体压强为 PB，速度为V。应用伯努利方程于同一流线上的、两点，则有 则 (3-43)022gpzgVgpzABgVgpgphBA22ghppvBA222022-7-2276 式（3-43）表明，只要测量出流体的运动全压和静压水头的差值h，就可以确定流体的流动速度。由于流体的特性，以及皮托管本身对流动的干扰，实际流速比用式(3-43)计算出的要小，因此，实际流速为 (3-44)式中 流速修正系数，一般由实验确定，=0.97。如果测定气体的流速，则无法直接用皮托管和静压管测量出气柱差来，必须把两根管子连接到一个形差压计上，从差压计上的液面差来求
48、得流速，如图3-18所示，则 用式(3-43)，则得 （3-45）ghV2)(液液ghppBA122液液液液ghhgV2022-7-2277图 3-18 用皮托管和静压管测量气体流速2022-7-2278 考虑到实际情况，(3-45a)在工程应用中多将静压管和皮托管组合成一件，称为皮托静压管，又称动压管，习惯上常简称它为皮托管，其示意图如图3-19所示。图中1点为总压测点，2点为静压测点，将总静压孔的通路分别连接于差压计的两端，则差压计的指示为总压和静压的差值，从而可由式(3-43)求得测点的流速。皮托-静压管的构造尺寸及使用时的连接方式如图3-20所示。12液液ghV2022-7-2279图
49、 3-19 皮托-静压管2022-7-2280图 3-20 皮托-静压管构造及连接方式2022-7-2281 二、文特里二、文特里(Venturi)流量计流量计 文特里流量计主要用于管道中流体的流量测量，主要是由收缩段、喉部和扩散段三部分组成，如图3-21所示。它是利用收缩段，造成一定的压强差，在收缩段前和喉部用形管差压计测量出压强差，从而求出管道中流体的体积流量。以文特里管的水平轴线所在水平面作为基准面。列截面1-1，2-2的伯努利方程 (3-46)由一维流动连续性方程 (3-47)gVgpgVgp20202222112121VAAV 2022-7-2282图 3-21 文特里流量计原理图2
50、022-7-2283 将式(3-47)代入到式(3-46)，整理得 (3-48)由流体静力学 (3-49)将式(3-49)代入到式(3-48)，则 (3-50)式(3-50)表明，若液，A2，A1已知，只要测量出h液，就可以确定流体的速度。流量为：(3-51)/(1)(2212212AAppV液液ghpp)(21)/(1)(22122AAhgV液液)/(1)(242122222AAhgdVAqV液液2022-7-2284 考虑到实际情况 (3-52)式中Cd为流量系数，通过实验测定。文特里流量计是节流装置中的一种，除此之外还有孔板，喷嘴等，其基本原理与文特里流量计基本相同，不再叙述。三、伯努利
