第四章——流体动力学分析基础课件

1、11. 1. 基本物理定律基本物理定律 u质量守恒定律质量守恒定律u能量守恒定律（热力学第一定律）能量守恒定律（热力学第一定律）u动量守恒定律（牛顿运动定律）动量守恒定律（牛顿运动定律）连续性方程连续性方程动量方程和动量矩方程动量方程和动量矩方程内维尔内维尔斯托克斯方程斯托克斯方程伯努利方程伯努利方程能量方程能量方程22. 2. 微分方法和积分方法微分方法和积分方法 微分方法：微分方法：将基本物理定律应用到流体微元或将基本物理定律应用到流体微元或微元控制体微元控制体上上, 可得到可得到微分形式的基本方程微分形式的基本方程，求解方程可得到物，求解方程可得到物理理 量的空间分布规律。量的空间分布规

2、律。 积分方法：积分方法：将基本物理定律应用到有限体积将基本物理定律应用到有限体积控制体控制体上，可得上，可得 积分形式的基本方程积分形式的基本方程。求解方程可得到物理量在。求解方程可得到物理量在 有限体积区域上的总体量的变化规律。有限体积区域上的总体量的变化规律。3系统：系统：一定质量的流体质点的集合一定质量的流体质点的集合, , 相当于热力学中的相当于热力学中的闭口系闭口系3. 3. 系统和控制体系统和控制体控制体：控制体：流场中确定的空间区域流场中确定的空间区域 。相当于热力学中的。相当于热力学中的开口系开口系, , 其其边界面称为边界面称为控制面控制面. 41. 1. 雷诺输运方程雷诺

3、输运方程描述了描述了系统内系统内流体参数变化与流体参数变化与控制体内控制体内流体参数的变化之流体参数的变化之间的关系。间的关系。定义定义B B为系统内任一物理量，为系统内任一物理量，为单位质量的该物理量，则为单位质量的该物理量，则dBdm或者或者BdmdVtuCVSttu50limssttttsBBdBdtt 系统内物理量随时间变化可以表示为：系统内物理量随时间变化可以表示为：uCVSIIIIII其中其中,sIIIIIc vIIIIttttttBBBBBB,sc vttBB6,0limc vIIIIc vttttsBBBBdBdtt 所以：所以：uCVS展开后有：展开后有：,000limlim

4、limc vc vIIIIttttttttttsBBBBdBdtttt 控制体内控制体内B B的变化率的变化率B B通过控制面通过控制面净流出率净流出率B B通过控制面通过控制面净流入率净流入率IIIIII7,0,limc vc vtttc vtc vBBdBddVtdtdt 控制体内物理量控制体内物理量B B随时间变化率：随时间变化率：由于控制体相对静止且固定不变由于控制体相对静止且固定不变,c vc vc vddVdVdVdtttuCVSIIIIII8单位时间内通过微元控制面的流出的体积通量为：单位时间内通过微元控制面的流出的体积通量为：IIIIIIinonidAodAtoViVooood

5、VVndA dt单位时间内通过微元控制面的流入的体积通量为：单位时间内通过微元控制面的流入的体积通量为：iiiidVV ndAdt物理量物理量B B的净流出率：的净流出率：*, c sBdVV n dA9得到雷诺输运方程得到雷诺输运方程.c vc ssdBdV n dAdtt 4-8雷诺输运方程雷诺输运方程。表示了系统内物理量。表示了系统内物理量B B随时间的变化率，随时间的变化率，等于控制体内该物理量随时间的变化率加上通过控制面该物等于控制体内该物理量随时间的变化率加上通过控制面该物理量的静流出率。理量的静流出率。式中式中B B是系统内任一物理量，是系统内任一物理量，是单位质量的该物理量，即

6、是单位质量的该物理量，即 =dB/dm；s系统，系统，c.v 控制体控制体，c.s 控制面，控制面，V 速度速度，n 控制面的外法线方向控制面的外法线方向102. 2. 雷诺输运方程的物理意义雷诺输运方程的物理意义流体流体质点质点参数参数B B的随体导数的随体导数 = = 当地导数当地导数 + + 迁移导数迁移导数流体流体系统系统参数参数B B的随体导数的随体导数 = B= Bc.vc.v对时间的导数对时间的导数 + B+ B的净流出率的净流出率以流体质点和空间坐标点为研究对象，适用于微分分析以流体质点和空间坐标点为研究对象，适用于微分分析以系统和控制体为研究对象，适用于控制体分析以系统和控制

7、体为研究对象，适用于控制体分析定常条件下定常条件下. c ssdBV n dAdt4-9定常条件下系统内物理量定常条件下系统内物理量B B的变化仅与通过控制面的流动有关，的变化仅与通过控制面的流动有关，与其内部状态无关，可以得到积分形式的控制方程与其内部状态无关，可以得到积分形式的控制方程111. 1. 连续性方程连续性方程质量守恒定律质量守恒定律系统内的流体在流动过程中质量不发生变化。系统内的流体在流动过程中质量不发生变化。 令雷诺输运方程中的物理量令雷诺输运方程中的物理量B B为系统的为系统的质量质量m m, , 则单位则单位质量的物理量质量的物理量=dB/dm=1=dB/dm=1，雷诺输

8、运方程为：，雷诺输运方程为：.c vc ssdmdV n dAdtt 4-10由质量守恒定律知，系统内的质量不变，（由质量守恒定律知，系统内的质量不变，（dm/dtdm/dt）s s=0=0，所以，所以.0c vc sdV n dAt 4-11积分形式的连续性方程积分形式的连续性方程，表示通过控制面的，表示通过控制面的净质量流出率净质量流出率 等于控制体内部等于控制体内部质量的减少率。质量的减少率。 适用于任何流体的定常和不定常流动。适用于任何流体的定常和不定常流动。 122. 2. 不可压缩流体的连续性方程不可压缩流体的连续性方程对于定常流动或者不可压缩流体，式（对于定常流动或者不可压缩流体

9、，式（4-114-11）可以简化为：）可以简化为：.0c sV n dA 4-12考虑图考虑图4-44-4所示的微元所示的微元流管流管内内不可压缩流体的流动，在流管不可压缩流体的流动，在流管壁面上壁面上（Vn）=0, ,截面截面1 1上上（V1n1）=-V1，截面，截面2上上（V2n2）=V2，因此，因此22110V dAV dA4-13对不可压缩流体的流动，通过控制面净流出的对不可压缩流体的流动，通过控制面净流出的体积流量体积流量 恒为零恒为零 13 对于任意有限截面的流管，如果对于任意有限截面的流管，如果 和和 为为A A1 1、A A2 2两个两个有效截面上的平均流速，则有有效截面上的平

10、均流速，则有1V2V1122V AV A4-13不可压缩流体不可压缩流体一维流动的连续性方程一维流动的连续性方程。【例【例4-14-1】已知油的密度为】已知油的密度为850kg/m850kg/m3 3，在内径为，在内径为0.2m 0.2m 的输油管的输油管道截面上的流速为道截面上的流速为2m/s2m/s，求另一内径为，求另一内径为0.05m 0.05m 的截面上的流速的截面上的流速及管道内的质量流量。及管道内的质量流量。【解】由不可压缩流体连续性方程【解】由不可压缩流体连续性方程 有有1122V AV A222112/20.2/0.0532(/ )VV ddm s其质量流量其质量流量211(

11、/4)53.4(/ )GVAd Vkg s 143. 3. 可压缩流体定常流动的连续性方程可压缩流体定常流动的连续性方程.0c sV n dA4-14对于可压缩的定常流动，由式对于可压缩的定常流动，由式4-11: 4-11: 对可压缩流体定常流动，通过控制面净流出的对可压缩流体定常流动，通过控制面净流出的质量流量质量流量 恒为零恒为零 对于任意有限截面的流管，如果对于任意有限截面的流管，如果 和和 为为A A1 1、A A2 2两个有两个有效截面上的平均流速，效截面上的平均流速， 1 1 2 2为两个有效截面上的密度，为两个有效截面上的密度，有有1V2V111222V AV A4-14a12m

12、mQQ4-14b或者或者15能量守恒定律能量守恒定律热力学第一定律热力学第一定律：系统内的：系统内的能量变化率能量变化率等于等于 单位时间内外界对系统所做的单位时间内外界对系统所做的功功加上单位时加上单位时 间内外界传递给系统的间内外界传递给系统的热量热量dEQWdt4-15dEdt系统能量对时间的变化率，是空间和时间的函数（系统能量对时间的变化率，是空间和时间的函数（状态量状态量）Q系统热量随时间的变化率，是时间的函数（系统热量随时间的变化率，是时间的函数（过程量过程量）系统吸热，系统吸热，Q Q为正值；系统放热，为正值；系统放热，Q Q为负值。为负值。 W系统与外界系统与外界作功作功随时间

13、的变化率，是时间的函数（随时间的变化率，是时间的函数（过程量过程量）环境对系统作功，环境对系统作功，W W为正值；系统对环境作功，为正值；系统对环境作功，W W为负值。为负值。 16对式中对式中4-154-15系统能量变化率应用雷诺输运方程，则系统能量变化率应用雷诺输运方程，则B=EB=E，=dE/dm=e=dE/dm=e.c vc ssdEe deV n dAdtt 联立式（联立式（4-154-15）得）得.c vc se deV n dAQWt 单位时间内输入系统的单位时间内输入系统的热量热量与环境对系统与环境对系统作功之和作功之和，等于控，等于控制体内制体内能量对时间变化率能量对时间变化

14、率加上通过控制体表面的加上通过控制体表面的能量流率能量流率。通常条件下，不考虑系统与外界的热交换，即认为通常条件下，不考虑系统与外界的热交换，即认为Q=0Q=04-164-17通过控制面的作功通过控制面的作功.nc sWV dA17其中表面应力其中表面应力 ，对于理想无粘流体，切应力为，对于理想无粘流体，切应力为零，且法向应力零，且法向应力 ，所以，所以nnpnppn . c sWp V n dA代入代入4-174-17，对于定常流动有，对于定常流动有.0c sepV n dA重力场中重力场中22uVeegz所以：所以：2.02unc sVpegzV dA重力场中重力场中理想理想流体流体定常绝

15、热定常绝热流动的能量方程流动的能量方程4-24181. 1. 伯努利方程伯努利方程将重力场中理想流体定常绝热流动的能量方程（将重力场中理想流体定常绝热流动的能量方程（4-244-24）应用到微）应用到微元流管中，在微元流管壁面上元流管中，在微元流管壁面上V Vn n=0=0， 在流入截面在流入截面A A1 1上，上，V Vn n=-V=-V1 1，在流出截面上在流出截面上V Vn n=V=V2 2，则，则21222211222111022uuAAVpVpVegzdAVegzdA在在被积函数在微元面上被积函数在微元面上积分，积分，2222112211022uuVpVpegzegz不可压理想流体与

16、外界无热交换的条件，内能、密度为常数不可压理想流体与外界无热交换的条件，内能、密度为常数4-262222112122VpVpgzgz4-271922Vpgzconst伯努利方程伯努利方程。表示了不可压缩理想流体在重力场中作定常。表示了不可压缩理想流体在重力场中作定常流动时，沿流线流动时，沿流线单位质量单位质量流体的动能、位置势能和压强势能之流体的动能、位置势能和压强势能之和守恒的规律。式中各项的单位为和守恒的规律。式中各项的单位为J/kg J/kg 4-28即即伯努利方程的伯努利方程的限制条件限制条件为：为：定常流动；定常流动；无粘流体（忽略粘性影响）；无粘流体（忽略粘性影响）；不可压缩流体；

17、不可压缩流体；沿流线沿流线 伯努利方程的条件虽然苛刻，但应用广泛伯努利方程的条件虽然苛刻，但应用广泛 20对对4-284-28式除以重力加速度式除以重力加速度g g，可以得到，可以得到单位重量单位重量流体的伯努利方程流体的伯努利方程22Vpzconstgg4-29式式4-29中，各项的单位为中，各项的单位为J/N，即米（，即米（m）。）。2. 2. 伯努利方程的物理意义伯努利方程的物理意义 V2/(2g) 为单位重量流体具有的动能为单位重量流体具有的动能速度水头速度水头；z 为单位重量流体所具有的位势能为单位重量流体所具有的位势能位置水头位置水头； p/(g)为单位重量流体的压强势能为单位重量

18、流体的压强势能压强水头压强水头。 因此伯努利方程描述了理想不可压缩流体在重力作用下作定常因此伯努利方程描述了理想不可压缩流体在重力作用下作定常流动时，沿同一流线上各点的单位重量流体所具有的位势能、压强流动时，沿同一流线上各点的单位重量流体所具有的位势能、压强势能和动能之和守恒。势能和动能之和守恒。21位置水头、压强水头和速度水头之和称为位置水头、压强水头和速度水头之和称为总水头总水头H。因此伯努。因此伯努利方程也可叙述为：理想不可压缩流体在重力作用下作定常流利方程也可叙述为：理想不可压缩流体在重力作用下作定常流动时，沿同一流线动时，沿同一流线(或微元流束或微元流束)上各点的单位重量流体所具有上

19、各点的单位重量流体所具有的位置水头、压强水头和速度水头之和为常数。的位置水头、压强水头和速度水头之和为常数。223. 3. 总流伯努利方程总流伯努利方程 式式4-28描述单位质量流体沿描述单位质量流体沿流线流线流动时总机械能守恒。在由流动时总机械能守恒。在由无数流线组成的流束中，将伯努利方程中三项机械能在无数流线组成的流束中，将伯努利方程中三项机械能在有效截面有效截面A上按质量流量积分，总机械能沿流束仍保持守恒，即上按质量流量积分，总机械能沿流束仍保持守恒，即 22AVpgzdQconst利用利用有效截面有效截面上的压强分布满足静力学规律，得到上的压强分布满足静力学规律，得到22ApVzQdQ

20、constgg 用总流有效截面上的平均速度用总流有效截面上的平均速度 代替不均匀的速度分布，为代替不均匀的速度分布，为此引入动能修正因子此引入动能修正因子， V2323322322AAAAAVdQV dAV dAgVV QV dAdQg 所以所以总流总流的伯努利方程为：的伯努利方程为： 22Vpzconstgg方程成立的限制条件是：方程成立的限制条件是： （1）忽略粘性摩擦；）忽略粘性摩擦； （2）不可压缩流体；）不可压缩流体； （3）定常流动；）定常流动； （4）无能量的输入和输出）无能量的输入和输出4-30244. 4. 伯努利方程的应用伯努利方程的应用（）小孔出流（）小孔出流 如图所示一

21、敞口水箱，侧壁下部开一小孔，假定水箱内水如图所示一敞口水箱，侧壁下部开一小孔，假定水箱内水位保持不变，求水从小孔流出的速度与孔口到液面的垂直距离位保持不变，求水从小孔流出的速度与孔口到液面的垂直距离之间的关系。之间的关系。 从自由液面上从自由液面上1和小孔和小孔2之间找一根流之间找一根流线，建立线，建立1-2两点之间的伯努利方程两点之间的伯努利方程 2222112122VpVpgzgz 由边界条件由边界条件V1=0，Z1=h, P1=Pa，Z2=0, P2=Pa,得到得到22Vgh4-30托里拆利公式托里拆利公式，液面上流体质点的位能全部转化为小孔出，液面上流体质点的位能全部转化为小孔出 流的

22、动能。流的动能。 12zh25（2）皮托管）皮托管 皮托管皮托管通过测量总压与静压之差来测量流体速度的一种装置通过测量总压与静压之差来测量流体速度的一种装置26 当液流流到测速管入口前的当液流流到测速管入口前的0点处，液流受到阻挡，流速变点处，液流受到阻挡，流速变为零，形成驻点。驻点的压强称为为零，形成驻点。驻点的压强称为总压总压，选择过，选择过0点的一个流点的一个流线上的另一点线上的另一点1，列伯努利方程，列伯努利方程2200110122VpVpzzgggg其中，其中，V0=0，Z0=Z1，p0= gh0，p1= gh1，可以得到，可以得到10122Vg hhg h4-31 工程上由于流体的

23、粘性使流动产生能量损失，以及皮托管本工程上由于流体的粘性使流动产生能量损失，以及皮托管本身对流动的干扰，实际流速比用上式计算出的要小，因此引入校身对流动的干扰，实际流速比用上式计算出的要小，因此引入校正系数正系数12Vg h4-31b27 如果测定气体的流速，则无法直接用皮托管和静压管测量如果测定气体的流速，则无法直接用皮托管和静压管测量出气柱差来，必须把两根管子连接到一个形差压计上，从差压出气柱差来，必须把两根管子连接到一个形差压计上，从差压计上的液面差来求得流速，此时计上的液面差来求得流速，此时121Vg h液 在工程应用中多将静压管和皮托管组合成一件，称为皮托在工程应用中多将静压管和皮托

24、管组合成一件，称为皮托静压管，又称动压管，习惯上常简称它为皮托管静压管，又称动压管，习惯上常简称它为皮托管28比托管的应用比托管的应用空速管空速管 29（3）文丘里管流量计）文丘里管流量计 文丘里流量计用于管道中流体的流量测量，主要是由收缩段、文丘里流量计用于管道中流体的流量测量，主要是由收缩段、喉部和扩散段三部分组成，如图所示。其原理是管道收缩，流速喉部和扩散段三部分组成，如图所示。其原理是管道收缩，流速增加，压强降低；通过测量压强的变化来求出管道中流体的体积增加，压强降低；通过测量压强的变化来求出管道中流体的体积流量。流量。30 在中心流线上列在中心流线上列1-2两点间的伯两点间的伯努利方

25、程努利方程2211221222VpVpzzgggg 当管道水平放置时当管道水平放置时z1=z2，根据不可压流体的连续性方程，根据不可压流体的连续性方程V1A1=V2A2，即，即V1=V2A2/A1，可得，可得1212224421212()2()1 (/) 1/ppppVAAdd所以，流量为所以，流量为221244212()41/dppQdd4-3231得到流量表达式的另一种形式得到流量表达式的另一种形式12ppg h由流体静力学压强分布关系式由流体静力学压强分布关系式224421241/dg hQdd4-32a工程上为了便于使用通常将两测压管做成工程上为了便于使用通常将两测压管做成U型管的形式

26、，同样考虑粘性引起的能量型管的形式，同样考虑粘性引起的能量损失后，损失后，2244212 ()41/ghQdddm32【例【例4-2】如图】如图4-9所示，水沿渐缩管道垂直向上流动。已知所示，水沿渐缩管道垂直向上流动。已知d1=0.3m, d2=0.2m，压力表显示相对压强为，压力表显示相对压强为p1=196kPa, p2=98.1kPa, h=2m，不计摩擦损失，试计算流量。，不计摩擦损失，试计算流量。【解】在渐缩管道的入口和出口建立伯努利方程【解】在渐缩管道的入口和出口建立伯努利方程2211221222VpVpzzgggg利用利用z1=0，z2=2m，以及连续性方程，以及连续性方程2441

27、222211/2ppVgzdd1/2122244212 ()/13.97(/ )1/ppgzVm sdd2322/40.439(/ )Qd Vms331. 1. 定常流动的动量方程定常流动的动量方程 用于求解流体边界上流体与固体的相互作用，是动量定用于求解流体边界上流体与固体的相互作用，是动量定理在流体流动问题上的应用。理在流体流动问题上的应用。动量定理动量定理 ：系统内流体动量对时间变化率等于作用在系统上的系统内流体动量对时间变化率等于作用在系统上的 外力矢量和。外力矢量和。sdmVFdt通过雷诺输运方程通过雷诺输运方程B=mV，则，则 =V，带入雷诺输运方程，带入雷诺输运方程.c vc s

28、sdmVVdV V n dAdtt 对于定常流动有对于定常流动有 . c sFV V n dA 4-3534式式4-35表示了作用在控制体上的合力等于流出、流入控制体的表示了作用在控制体上的合力等于流出、流入控制体的净动量流率。上式是净动量流率。上式是 一个矢量方程，通常进行分解后求解一个矢量方程，通常进行分解后求解 .xc sFu V n dA.yc sFv V n dA.zc sFw V n dA 式式4-36表明定常流条件下，作用于控制体上表明定常流条件下，作用于控制体上合力合力沿三个坐沿三个坐标轴的标轴的分量分量与流出、流入控制面的与流出、流入控制面的净动量流率净动量流率在三个坐标轴的

29、在三个坐标轴的分量分量相等。相等。 4-36注注：（：（1）Fx, Fy, Fz, u, v, w 可正可负，取决于坐标轴的方向可正可负，取决于坐标轴的方向 （2）（）（Vn） 当流进控制体为负，流出控制体为正当流进控制体为负，流出控制体为正35（1）合力项）合力项 作用在系统上的作用在系统上的合力合力包括作用在控制体上的包括作用在控制体上的质量力质量力和控制面上和控制面上的的面积力面积力 msFFF若流体仅处于重力场中，则质量力为若流体仅处于重力场中，则质量力为mFmgk面积力包括（面积力包括（1）与固体接触的控制面受到固体的作用力）与固体接触的控制面受到固体的作用力 （2）周围流体接触作用

30、产生的面积力）周围流体接触作用产生的面积力其中质量力其中质量力,mc vFfd,sc sFpn dA其中其中-n表示流体在控制面上受到的作用力指向控制面内部表示流体在控制面上受到的作用力指向控制面内部36（2）净动量流率项）净动量流率项 当选择合适的控制体使当选择合适的控制体使V和和在控制面上均匀分布时，净动量流率为在控制面上均匀分布时，净动量流率为.nnc soutinV V n dAVV AVV A 其中其中Vn= Vn是控制面上的法向速度的是控制面上的法向速度的值值 ， Vn A表示通过控制面表示通过控制面的体积流量的体积流量Q，所以，所以. c soutinV V n dAQVQV 分

31、解到三个坐标轴上为分解到三个坐标轴上为 4-38a.xoutinc sFu V n dAQuQu.youtinc sFv V n dAQvQv.zoutinc sFw V n dAQwQw4-38b4-38372. 2. 动量方程的应用动量方程的应用 动量方程应用要点：动量方程应用要点：（1）选择合适的控制体，包含尽可能多的已知条件）选择合适的控制体，包含尽可能多的已知条件（2）选择确定的坐标系以简化条件）选择确定的坐标系以简化条件（3）将外力、动量流率向确定的坐标轴投影，求该坐标方向上）将外力、动量流率向确定的坐标轴投影，求该坐标方向上 的分量的分量（4）假定待求力的方向与所选坐标方向一致，

32、若结果为负，则）假定待求力的方向与所选坐标方向一致，若结果为负，则 说明力的实际方向与假定方向相反说明力的实际方向与假定方向相反38【例【例4-3】如图】如图4-10所示，水从固定喷嘴定常流出，垂直冲击一所示，水从固定喷嘴定常流出，垂直冲击一平板。水离开喷嘴的速度平板。水离开喷嘴的速度V1=20m/s, 喷嘴出口面积喷嘴出口面积A1=0.005m2。假定水冲击平板后沿平板流动。水的密度为假定水冲击平板后沿平板流动。水的密度为1000kg/m3,试确定试确定支撑这块平板所需的水平力。支撑这块平板所需的水平力。【解】选择包围平板并切割【解】选择包围平板并切割水流和支撑架的矩形为控制水流和支撑架的矩

33、形为控制体，沿水射流方向建立体，沿水射流方向建立x轴，轴，垂直方向为垂直方向为y轴。支撑架对轴。支撑架对控制体有作用力控制体有作用力Rx，假定其，假定其沿沿x正向，列正向，列x方向上的动量方向上的动量方程：方程：39支撑该平板的水平力大小为支撑该平板的水平力大小为2000（N) ，方向与，方向与x轴反向轴反向 11 12000( )xRVVAN所以所以.mxsxc sFFu V n dA重力场中重力场中0mxF sxxFR而面积力而面积力动量流率项有一个流入项和两个流出项，其中流出项与动量流率项有一个流入项和两个流出项，其中流出项与x轴轴垂直，在垂直，在x轴分量为零轴分量为零 11 1.out

34、inc su V n dAQuQuVVA40【例】【例】 水平放置在混凝土支座上的变直径弯管，弯管两端与等水平放置在混凝土支座上的变直径弯管，弯管两端与等直径管相连接处的断面直径管相连接处的断面1-1上压力表读数上压力表读数p1=17.6104Pa，管中，管中流量流量Qv=0.1m3/s，若直径，若直径d1=300，d2=200，转角，转角=600，如图所示。求水对弯管作用力如图所示。求水对弯管作用力F的大小的大小41【解】【解】 水流经弯管，动量发生变化，必然产生作用力水流经弯管，动量发生变化，必然产生作用力F。而。而F与与管壁对水的反作用力管壁对水的反作用力R平衡。管道水平放置在平衡。管道

35、水平放置在xoy面上，将面上，将R分解分解成成Rx和和Ry两个分力。取管道进、出两个截面和管内壁为控制面，两个分力。取管道进、出两个截面和管内壁为控制面，如图所示，坐标按图示方向设置。如图所示，坐标按图示方向设置。（1）根据连续性方程可求得：）根据连续性方程可求得：12210.1 41.42(/ )/40.3VQvm sd22220.1 43.18(/ )/40.2VQvm sd42得得2-2面上压强面上压强 ：2232112()/217.2 10 ()ppvvPa（3）写出动量方程）写出动量方程假定壁面对控制体内水的作用力为假定壁面对控制体内水的作用力为Rx、Ry，其方向如图，其方向如图，沿

36、沿x轴方向列动量方程轴方向列动量方程1122cosxVoutVinp Ap ARQ uQ u212211(cos )cos568( )xVRQvvp Ap AN （2） 列管道进、出口的伯努利方程列管道进、出口的伯努利方程gvgpgvgp2222221143沿沿y轴方向列动量方程轴方向列动量方程111sin(0sin )yVp ARQv111sinsin10880()yVRp AQ vN管壁对水的作用力合力管壁对水的作用力合力2222( 0.568)10.8810.89()xyRRRkN水流对弯管的作用力水流对弯管的作用力F与与R大小相等，方向相反大小相等，方向相反441. 1. 定常流动的角

37、动量方程（动量矩方程）定常流动的角动量方程（动量矩方程） 角动量方程角动量方程描述了作用于流体系统上的描述了作用于流体系统上的力矩力矩与与角动量角动量随时间的随时间的变化关系。可以确定流体和外界之间作用力的变化关系。可以确定流体和外界之间作用力的位置位置。动量矩定理：动量矩定理：单位时间内流体单位时间内流体系统对转动轴的动量矩（系统对转动轴的动量矩（角动角动量量）的变化率，等于作用于系）的变化率，等于作用于系统上所有外力对同一轴的力矩统上所有外力对同一轴的力矩之和。之和。ooodMrFrmVdt4-4445.oooc sMrFrVV n dA4-47 表示定常流动时作用在控制体（系统）上所有力

38、的表示定常流动时作用在控制体（系统）上所有力的力矩矢力矩矢量和量和，等于流入、流出控制面的，等于流入、流出控制面的净角动量流率。净角动量流率。BrmV令雷诺输运方程中令雷诺输运方程中 ，则，则 rV.oooc vc ssdrmVrVdrVV n dAdtt 将上式代入将上式代入4-44得得在定常条件下在定常条件下.ooooc vc sMrFrVdrVV n dAt 462. 2. 角动量方程的应用角动量方程的应用 【例【例4-6】一草坪洒水器在水平】一草坪洒水器在水平面（面（xy平面）内绕平面）内绕z轴等角速度轴等角速度旋转，转速为旋转，转速为120r/min。如图。如图所示，水从中心垂直管进

39、入，所示，水从中心垂直管进入，经过转臂两端的喷嘴喷出，进经过转臂两端的喷嘴喷出，进水流量水流量Qi=0.006m3/s，喷嘴出，喷嘴出口截面积口截面积A0=0.001m2，撒水臂，撒水臂长长R=0.2m。（。（1）为使洒水器）为使洒水器维持该等角速度旋转，外界需维持该等角速度旋转，外界需加的阻力矩为多少？（加的阻力矩为多少？（2）如果）如果阻力矩为零，则洒水器的旋转阻力矩为零，则洒水器的旋转角速度将为多少？角速度将为多少？47【解】建立如图所示坐标系，【解】建立如图所示坐标系，z轴为轴为旋转轴，选择洒水器旋转臂与旋转轴，选择洒水器旋转臂与x轴重轴重合时围绕洒水器建立控制体。合时围绕洒水器建立控

40、制体。（1）设维持）设维持120r/min需要的力矩为需要的力矩为ToosomMrFrFT其中控制体周围大气压产生的压力对其中控制体周围大气压产生的压力对O点力矩为零，控制体内点力矩为零，控制体内质量力对质量力对O点对称，其力矩也为零。假设外力矩点对称，其力矩也为零。假设外力矩T沿沿z轴正向，轴正向，所以所以.ooc sMTrVV n dV而而1122.()()()ooooooooiic srVV n dArVQrVQrVQ48其中其中o1，o2表示两个出口，表示两个出口，i表示进口。对于进口处表示进口。对于进口处ri=0(速度平速度平行于行于)z轴，两个出口的动量矩大小相等，方向相同，所以有

41、轴，两个出口的动量矩大小相等，方向相同，所以有2()()oooooiTkrVQrVQVo为洒水器出口的为洒水器出口的绝对速度绝对速度（相对于静止的坐标系的速度），（相对于静止的坐标系的速度），等于出口的等于出口的相对速度相对速度减去喷嘴运动的减去喷嘴运动的牵连速度。牵连速度。()()(3.02.51)0.492iooutoQVVRjRjjjA其中其中0.2orRii所以所以()(0.20.49 )iooiTkQ rVQij31000 6 10 (0.2 0.49)0.588kk 49（2）阻力矩为零时的转速）阻力矩为零时的转速()0iooTkQ rV()02ioQRiRjA所以所以15(/ )

42、2ioQrad sA R转速转速60143.24( /min)2nr501. 1. 连续性方程连续性方程 方程的导出：方程的导出：积分形式的连续性方程的一般形式积分形式的连续性方程的一般形式.0c vc sdVV n dAt.()c sc vV n dAV dV高斯定理高斯定理: 物理量对控制面的面积分物理量对控制面的面积分,等于该物理量的等于该物理量的散度散度在在 控制面所包围的控制体内的体积分控制面所包围的控制体内的体积分dSAdvAn所以所以.()0c vVdVt 对积分形式的连续性方程进行数学变换对积分形式的连续性方程进行数学变换51由于控制体选择的任意性，可知被积函数必为零由于控制体

43、选择的任意性，可知被积函数必为零()0Vt4-54微分形式的微分形式的连续性方程连续性方程一般形式一般形式在直角坐标系中可以分解为在直角坐标系中可以分解为()()()0uvwtxyz4-56展开后并利用随体导数的概念，可得展开后并利用随体导数的概念，可得0DVDt4-57方程的适用条件：满足连续介质假设的方程的适用条件：满足连续介质假设的任何流动任何流动52 微元六面体的质量守恒分析微元六面体的质量守恒分析 如图所示，设流体流过以如图所示，设流体流过以M(x,y,z)为基点，以为基点，以dx,dy,dz为为边长的微元控制体。边长的微元控制体。 在在 t 时间内沿时间内沿x方向净流出控制体（流出

44、质量减去流入质量）方向净流出控制体（流出质量减去流入质量）的质量为的质量为()udxdydz tx同理在同理在 t 时间内沿时间内沿y方向和方向和z方向净流出控制体的质量为方向净流出控制体的质量为()vdxdydz ty()wdxdydz tz53 在在t时间内，微元六面体流体的质量变化为时间内，微元六面体流体的质量变化为t dxdydzdxdydzdxdydz ttt 按质量守恒定律，在按质量守恒定律，在t 时间内沿三个方向净时间内沿三个方向净流出流出控制体的控制体的总质量应等于控制体内总质量应等于控制体内减少减少的质量：的质量： ()()()uvwdxdydz tdxdydz txyzt化

45、简后可化简后可得得0zwyvxut或者或者()0Vt4-5454对于定常流动，由式对于定常流动，由式4-54可得连续性方程变为：可得连续性方程变为：()0V4-58 对于不可压缩流体的流动（定常或者非定常），由式对于不可压缩流体的流动（定常或者非定常），由式4-57可得微分形式的连续性方程为：可得微分形式的连续性方程为：0V4-60 表示了不可压缩流体流动时，速度的表示了不可压缩流体流动时，速度的散度散度为零。在直角为零。在直角坐标系中坐标系中0uvwxyz4-610 xxyyzz 或者或者55【例【例4-8】不可压缩二维平面流动，】不可压缩二维平面流动，y方向上的速度分量为方向上的速度分量为

46、v=y2-y-x，求，求x方向的速度分量方向的速度分量u，假定，假定x=0时时u=0【解】将不可压缩流体的连续性微分方程应用到二维流动【解】将不可压缩流体的连续性微分方程应用到二维流动0uvxy210uyx 将将y方向上的速度分量方向上的速度分量v代入得代入得积分后有积分后有(12 )( )uy xf y利用边界条件利用边界条件x=0时时u=0，得得f(y)=0所以所以u=(1-2y)x=x-2xy562. 2. 内维尔内维尔斯托克斯方程斯托克斯方程 对图示微元六面体应用对图示微元六面体应用牛顿第二定律，有牛顿第二定律，有DVFdxdydzDt其中作用在微元六面体上的质量力其中作用在微元六面体

47、上的质量力()mxxdFf dxdydz()yxxxzxsxpdFdxdydzxyz作用在微元六面体上的面积力包括法向力和切向力作用在微元六面体上的面积力包括法向力和切向力考虑考虑x方向上的动量平衡方向上的动量平衡()()xmxsxDuFdFdFdxdydzDt4-6457将质量力和面积力代入到将质量力和面积力代入到4-64并化简后得到运动微分方程：并化简后得到运动微分方程：yxxxzxxpDufDtxyz 其中应力的第一个下标表示应力其中应力的第一个下标表示应力作用面的法线方向，作用面的法线方向，第二第二个下标表示个下标表示应力的方向。应力的方向。（注：这里的法向应力不等于静压强）（注：这里

48、的法向应力不等于静压强）4-68同理可以得到同理可以得到y方向和方向和z方向上的运动微分方程方向上的运动微分方程xyyyzyypDvfDtxyzyzxzzzzpDwfDtxyz4-694-7058其中法向应力其中法向应力xxxxpp yyyypp zzzzpp 上式中上式中p为流体的静压强，为流体的静压强，xxxx、yyyy、zzzz，为流体的粘性，为流体的粘性变形引起的法向应力。因此式（变形引起的法向应力。因此式（4-684-68）式（式（4-704-70）变为）变为yxxxzxxDupfDtxxyzxyyyzyyDvpfDtyxyzyzxzzzzDwpfDtzxyz4-68a4-70a4-

49、69a59写成矢量形式为：写成矢量形式为：而而ijij为作用在微元六面体上的为作用在微元六面体上的黏性应力张量黏性应力张量xxxyxzijyxyyyzzxzyzzijijDVffpDt 4-71以应力形式表示的粘性流体的以应力形式表示的粘性流体的运动微分方程运动微分方程其中其中ijij为作用在微元六面体上的为作用在微元六面体上的应力张量应力张量xxxyxzxxxyxzijyxyyyzyxyyyzzxzyzzzxzyzzpppppp60（1）理想流体的欧拉运动微分方程）理想流体的欧拉运动微分方程对于理想流体对于理想流体ijij=0=0，式，式4-714-71简化为简化为DVfpDt在笛卡尔坐标系

50、中可以分解为在笛卡尔坐标系中可以分解为4-72xDupfDtxyDvpfDtyzDwpfDtz4-73理想流体的理想流体的欧拉运动微分方程欧拉运动微分方程61（2）黏性流体的内维）黏性流体的内维-斯托克斯方程斯托克斯方程 对于流体的三维流动，斯托克斯提出了对于流体的三维流动，斯托克斯提出了广义牛顿内摩擦广义牛顿内摩擦定律，定律，给出了给出了应力应力和和应变应变之间的关系之间的关系22()3ijijijijVp1()2jiijjiuuxx10ijijij对于切向应力对于切向应力2()2xyxyyxyxuvyx2()2yzyzzyzywvyz4-742()2xzzxxzzxuwzx62 上式表明流

51、体的法向应力不仅与流体的静压强相关，而上式表明流体的法向应力不仅与流体的静压强相关，而且与流体的且与流体的线变形速率线变形速率和和速度的散度速度的散度有关，而三项和为有关，而三项和为22()3xxupVpx22()3yyvpVpy22()3zzwpVpz4-75对于法向应力对于法向应力1()3xxyyzzpppp 式式4-76表明粘性流体三个相互垂直方向上的法向应力的表明粘性流体三个相互垂直方向上的法向应力的平均值等于流体静压强的值平均值等于流体静压强的值4-7663 将以应变形式表示出来的应力（式将以应变形式表示出来的应力（式4-74式式4-75）代入）代入粘性流体的运动微分方程（式粘性流体

52、的运动微分方程（式4-68式式4-70）并整理后得：）并整理后得：22222211()3xDupuuufVDtxxyzx22222211()3yDvpvvvfVDtyxyzy22222211()3zDwpwwwfVDtzxyzz内维内维-斯托克斯斯托克斯（Navier-Stokes）方程的一般形式，期）方程的一般形式，期中中为为拉普拉斯拉普拉斯算子算子写成矢量形式为：写成矢量形式为：11()3dVfpVVdt 4-7964对于对于不可压缩流体不可压缩流体N-S方程简化为方程简化为1dVfpVdt 对于对于理想无粘流体理想无粘流体N-S方程简化为方程简化为欧拉运动微分方程欧拉运动微分方程4-72

53、，对，对于于静止流体静止流体N-S方程简化为静止流体的方程简化为静止流体的欧拉平衡微分方程欧拉平衡微分方程。4-80注：注： 理论上理论上N-S方程包含方程包含u,v,w和和p四个未知量，结合流体的连四个未知量，结合流体的连续性方程可以构成封闭的方程组，完成对流场的积分求解续性方程可以构成封闭的方程组，完成对流场的积分求解 实践中仅能对某些简单的流动在适当的假设下进行求解实践中仅能对某些简单的流动在适当的假设下进行求解 对于可压缩流体，还要结合对于可压缩流体，还要结合状态方程状态方程求解求解 当涉及到作功、传热和内能变化时需要结合当涉及到作功、传热和内能变化时需要结合能量方程能量方程求解求解N

54、-S方程在方程在柱坐标柱坐标和和球坐标球坐标下的形式见教材和参考书下的形式见教材和参考书653. 3. 能量方程（简介）能量方程（简介） 类似于连续性方程和运动方程的导出过程，将热力学第类似于连续性方程和运动方程的导出过程，将热力学第一定律应用到微元控制体一定律应用到微元控制体DEDeQ WdxdydzDtDt 对该微元体进行传热、与外界作功和内能变化进行分析后对该微元体进行传热、与外界作功和内能变化进行分析后可以得到可以得到与外界的与外界的热交换热交换内热源内热源表面力表面力作功作功系统内能系统内能变化率变化率()()ijjiiiDeTf VuqDtxxx 质量力质量力作功作功664. 4.

55、 基本微分方程组的定解条件基本微分方程组的定解条件 对于工程流体力学常见的流动问题利用连续性方程和运动方对于工程流体力学常见的流动问题利用连续性方程和运动方程即可联立积分求解，为了给出程即可联立积分求解，为了给出特解特解需要给出需要给出定解条件定解条件以确定积以确定积分常数。分常数。（1）初始条件）初始条件 非定常流动，待求变量在非定常流动，待求变量在t=t0时刻的空间分布时刻的空间分布00( , ,)( , )V x y z tVx y z00( , ,)( , )p x y z tpx y z00( , ,)( , )x y z tx y z00( , ,)( , )T x y z tTx

56、 y z定常流动不需要初始条件定常流动不需要初始条件67（2 2）边界条件）边界条件 包围流场的边界上的流动参量值。三种典型的边界条件：包围流场的边界上的流动参量值。三种典型的边界条件：壁面无滑移条壁面无滑移条件件 粘性流体与固体壁面接触时，在壁面上流体速度等于固体粘性流体与固体壁面接触时，在壁面上流体速度等于固体壁面的速度壁面的速度wVV对无粘性流体，无需满足无滑移条件，但法向速度仍应连续对无粘性流体，无需满足无滑移条件，但法向速度仍应连续68 进口与出口条件进口与出口条件对对内流内流问题，进口和出口的速度和压强分布通常由实验给出；问题，进口和出口的速度和压强分布通常由实验给出；对对外流外流问题，通常用无穷远处的速