船舶流体力学课件4.

主讲教师：孙雷宗智船舶工程学院B4.1流体系统的随体导数B4.2积分形式的连续性方程B4.3伯努利方程及其应用B4.4积分形式的动量方程及其应用B4.5积分形式的动量矩方程B4.6积分形式的能量方程B4积分形式的基本方程B4积分形式的基本方程积分形式的流体力学基本方程描述空间有限体积域上的流体运动规律，主要涉及流体质量、动量、动量矩和能量等物理量在有限体积域上的积分值（广延量）随时间和位置的变化规律，它在工程上有广泛应用。主要内容：流体系统的随体导数；积分形式的连续性方程、动量方程、动量矩方程和能量方程及其应用，伯努利方程及其应用等。重点：（1）有限控制体分析，输运公式；（2）有多个一维出入口的控制体上的连续性方程；（3）伯努利方程；（4）有多个一维出入口的控制体上的定常动量方程等。B4.1流体系统的随体导数系统广延量由于为流体系统内物理量的空间分布函数，在系统（system）上积分：称为系统广延量。当取密度、动量、动量矩和能量函数时，分别可得系统质量、系统动量、系统动量矩和系统能量等。),(tr()(,)syssysNtrtd控制体广延量()(,,,)CVCVNtxyztd控制体表面为CS，一流体系统sys（实线包围区域）在t时刻刚好与控制体重合，以后流体系统可以与控制体形状不同。右图为控制体形状变化示意图：B4.1流体系统的随体导数有限控制体分析，输运公式在流场中取一固定不变形的有限控制体CV（图中虚线包围的区域）B4.1流体系统的随体导数t时刻物理量的空间分布函数（单位体积之值），在系统上的积分1452345()()ttCV控制体控制面2345CSAAAA()(,)syssysNtrtd由时间导数的定义，系统广延量的时间导数可表示为00d11lim()()limdddsystttttNNttNtttt由于控制体积分区域可分割成数块，()tt()tCVB4.1流体系统的随体导数23045023400d1limddddddd1limdd11limddlimddsysCVtttttCVtttCVCVttttttttNttttt5tt右端第一项代表控制体广延量对时间的导数01IlimdddCVCVCVttttttB4.1流体系统的随体导数230231IIlim()d()d()d()dAAtttAAAtAttAAvnvnvnvn右端第三项代表单位时间内通过控制面流出控制体的广延量（正值）450451IIIlim()d()d()d()dAAtttAAAtAttAAvnvnvnvnoutd()ddAtvn4501IIIlimddtttt右端第二项代表单位时间内通过控制面流入控制体的广延量（负值）2301IIlimddttttind()ddAtvnB4.1流体系统的随体导数将I，II，III式代入原系统广延量的时间导数公式，并用D/DT代替d/dt()sysCVCSDNdvndADTtrr上式被称为雷诺输运公式，简称输运公式。将II与III相加可得2345II+III()d()dAAAACSAAvnvn上式代表单位时间内通过控制面净流出控制体的广延量。类似于流体质点的随体导数（质点导数）概念，用控制体上的欧拉坐标表示流体系统的随体导数，关系式为：B4.1流体系统的随体导数()sysCVCSDNdvndADttsysDNDtdCVt()CSvndA表示系统与控制体重合时系统广延量对时间的随体导数，又称系统导数；表示控制体广延量随时间的变化率，又称当地变化率，反映流场的不定常性（定常时为零）；表示通过控制面净流出控制体的广延量流量，又称为迁移变化率，反映流场的不均匀性（均匀时为零）。定常流场输运公式上式表明在定常流场中，当系统与控制体重合时，系统广延量的变化只取决于控制面上的流动，与控制体内的流动无关（见下图）。()sysCSDNvndADtB4.1流体系统的随体导数思考题：运输公式：是对固定控制体导出的，若控制体作匀速运动时，下面哪个结论是对的：（A）仍然适用；（B）不再适用；（C）形式不变，但需将迁移项中v改为相对速度vr。B4.1流体系统的随体导数()sysCVCSDNdvndADttB4.2积分形式的连续性方程上式称为积分形式的连续性方程，适用于任何流体的定常和不定常流动。设，系统质量为(,)rtsyssysmd根据质量守恒定律：dd0ddsyssysmdtt()0CVCSdvndAt由输运公式可得：上式表明：通过控制面净流出的质流量等于控制体内流体质量随时间的减少率。B4.2.1固定控制体不可压缩流体实际上，对固定不变形的控制体，上面式子中的当地项中微分和积分运算可变换，迁移项中为绝对速度。v当密度为常数时，式中当地项为零，迁移项中密度项可消去，得上式的物理意义是：对不可压缩流体的流动，从任何固定不变形的控制面净流出的体积流量恒为零。CSA0d)(nv0ddCSCVAtnv对不可压缩流体一维流管流动B4.2.1固定控制体()()0outinoutinvndAvndA12QQ2211VAVA令截面1,2上的流量大小分别为Q1,Q2，由流量公式可得由平均速度公式可得早在16世纪初，达.芬奇就发现了这一规律。outinQQB4.2.1固定控制体若控制面上有多个出入口，设出入口的流量大小为Qout,Qin，由前面的公式可得()()outinVAVA思考题：对于连续性方程：的说法，下列哪个是对的（）（A）仅适用于不可压缩流体的定常流动的；（B）也适用于不可压缩流体的不定常流动；（C）适用于任何流体的定常流动。CSdAnv0)(B4.2.1固定控制体可压缩流体定常运动B4.2.1固定控制体()0CSvndA()()outinVAVA对密度可变流体的定常流动，可得上式的物理意义是：对可压缩流体定常流动，从任何固定不变形的控制面净流出的质流量恒为零。对一维流管流动，设出入口的质量流量大小分别为和，从质量流量公式可得outminminoutmmB4.2.1固定控制体对有多个出入口的控制面上的定常流动，由前面的公式可得inoutmm()()outinVAVA例题B4.2.1：主动脉弓流动：多个一维出入口连续性方程已知：下图是人主动脉弓模型示意图。血液从升主动脉1经主动脉弓流向降主动脉5，方向改变约180°，主动脉弓上分支出头臂干动脉2，左颈总动脉3和左锁骨下动脉4。设所有管截面均为圆形，管直径分别为d1=2.5cm,d2=1.1cm,d3=0.7cm,d4=0.8cm,d5=2.0cm。已知平均流量分别为Q1=6L/min,Q3=0.07Q1,Q4=0.04Q1,Q5=0.78Q1。试求：(1)管2的平均流量Q2；(2)各管的平均速度（用cm/s表示）。解：由取图中虚线所示控制体，有多个出入口。血液按不可压缩流体处理，由式∑Qout=∑QinQ1=Q2+Q3+Q4+Q5例题B4.2.1：主动脉弓流动：多个一维出入口连续性方程（1）管2的流量为Q2=Q1－(Q3+Q4+Q5)=Q1－(0.07+0.04+0.78)Q1=0.11Q1=0.66L/min（2）各管的平均速度为cm/s4.20s/min60cm5.2/lcm1000l/min644232111dQVcm/s6.11s/min60cm.11/lcm1000l/min66.044232222dQVcm/s2.18s/min60cm.70/lcm1000l/min6.07044232333dQVcm/s0.8s/min60cm.80/lcm1000l/min6.04044232444dQVcm/s8.24s/min60cm.02/lcm1000l/min6.78044232555dQVB4.2.2运动控制体无论是惯性系还是非惯性系，只要将迁移项中的速度改为相对于控制体的相对速度，即可得运动控制体形式的连续性方程：0)(dCVCSdAtnvr对具有多个一维出入口的定常流动为()()routrinVAVA上两式常在旋转控制体（如流体机械）中运用。思考题：所谓非惯性系是仅指：（A）做加速运动的控制体；（B）做匀角速度旋转的控制体；（C）做非匀角速度旋转的控制体；（D）包括以上三个答案。B4.2.2运动控制体相对于惯性系（静止或匀速运动的参考系）加速运动的参考系称为非惯性系参考系。地球有自转和公转，我们在地球上所观察到的各种力学现象，实际上是非惯性系中的力学问题。已知：下图为洒水器示意图。臂长R=150mm，喷水管面积A=40mm2，喷口偏转角水从中心转轴底部流入，总流量Q=120mL/s，从两喷口流出。喷管角速度为=500转/分求：（1）管内水流的相对速度Vr。（2）管口水流的绝对速度V。解：取包围喷管，并与喷管一起旋转的控制体，如图中虚线所示。对站在控制体上的观察者而言，水以速度Vr沿两支喷管做定常直线流动。由下式：30()()routrinVAVA111222rrVAVAQ例题B4.2.2：洒水器：运动控制体连续性方程可得水为不可压缩流体，且，由两臂对称方程，上式化为：管内相对速度为：喷口的牵连速度为：由喷口的速度矢量合成，绝对速度为：2112rrrVVVAAA212rVAQ6362120010/15/22(4010)rQmsVmsAmsmmsrRU/85.7)15.0(min/602min)/500(221/2221/2[2cos][(15/)(7.85/)2(15/)(7.85/)cos30]9.1/rrVVUVUmsmsmsmsms例题B4.2.2：洒水器：运动控制体连续性方程B4.3伯努利方程及其应用伯努利方程首次以动能与压强势能相互转换的形式确定了流体运动中速度与压强之间的关系。伯努利方程由伯努利（D.Bernouli，1738）首先提出，后来由欧拉（L.Euler）完善其理论推导过程。B4.3.1沿流线的伯努利方程沿流线的欧拉运动方程在无粘性流体的重力流场中沿流线S取一圆柱形体积元控制体（如图）,控制元长δs,端面面积为δA;两端面上的压强分别为p和p+δp,重力为gδAδs,在流线切线方向（即速度方向）运用牛顿第二定律可得d(,)δδcosδ(δ)δδδdpvatgAspApsAAsst整理后取极限可得：1d(,)cosdpvstgstB4.3.1沿流线的伯努利方程d(,)dvstvvvtts由几何关系coszs将流体元的加速度表达为欧拉形式代回原式得：1d(,)dzpvstvvgvsstts式中s为流线坐标，z为高度坐标，p为圆柱形体积元端面压强，v为圆柱形体积元速度。B4.3.1沿流线的伯努利方程将上式沿流线积分，可得：2dd2vvpsgzt