伯努利方程的应用

1、，伯努利方程及其应用伯努利，1738,瑞士。动能与压强势能相互转换。沿流线的伯努利方程将牛顿第二定律应用于控制体内的流体元，沿流线切线方向 整理后 因为将流体元的加速度转换成欧拉形式的加速度，沿流线的质点导数为 则导出此式为一维欧拉方程，使用下述关系将方程沿流线积分。两边乘以ds得： 沿流线积分此式为欧拉方程的积分式，适合于可压、无粘不定常运动。对于不可压定常流动，则可简化为此式为伯努利方程，三项分别表示单位质量流体具有的动能、位 I , .: _ . I置势能和压强势能。即总机械能守恒。应用伯努利方程时常采用沿流线上任两点的总机械能值相等的形式。伯努利方程使用的限制条件（1）无粘性流体，（2

2、）不可压流体（3） X J Z jT Ii II定常流（4）沿流线。加入能量损失就可适应粘性流体。皮托（pitot ）测速管：总压强与动压强皮托测速管又称为皮托-静压管，简称皮托管，为纪念法国人皮托 1 I命名。皮托测速管由粗细两根同轴的圆管组成，细管（直径约为1.5 mm）前端开孔（O点），粗管（直径约为6 mg在距前端适当 长距离处的侧壁上开数个小孔（B点），在孔后足够长距离处两管 弯9。成柄状.测速时管轴线沿来流方向放置.设正前方的流速保 持为V,静压强为p,流体密度为P。粗细两管中的压强被引入U 形测压计中，U形管中液体密度Pm。试求用U形管液位差 的 表 不流速v的关系式。解：设流动

3、符合不可压缩无粘性流体定常流动条件。从皮托管正 前方A点到端点0再到侧壁孔B点的A 0 B线是一条流线，A点 的速度和压强分别为v和p ,沿流线A0段按（B4. 3. 4）式列伯努利方程22VP V。 P。gzAgzo272D在皮托管端点O,流体速度降至零v0=0,称为驻点，P0称为驻点压 强，U形管右支管测到的即是驻点压强.由于zA=zo,由（a）式可 得上式中2pv2称为动压强，为流体质点的动能全部转化为压强势能时应具有的压强.（b）式表明驻点压强为静压强和动压强之和， 故po称为总压强.由（b）式动压强可表为一，一 A ,C由于皮托管较细，流线上的AB两点的位置差可忽略，伯努利方程为因V

4、b = V,由上式Pb =p,即U形管左支通过皮托管侧壁小孔测到的是当地静压强.在U形管内列压强关系式一 X. y r- s._i I i- i可得由于实际流体具有粘性及皮托管加工误差等原因，流体动压强转化为U形管内液位差读数存在误差，需乘上一个修正系数k，由（c） , （d）式可得k称为皮托管系数，可通过用标准皮托管作标定测量后确定.由（e）式可得小孔出流：托里拆利公式及缩颈效应从一个大容器侧壁下部距离液面为h处开一个小孔，设液面水位不变，求出流速度和出流流量。解：沿流线法向方向的速度压强关系式由牛顿第二定律：得2pVg A ncos p A - (p考虑到几何关系，有整理，得p） A= A

5、 -A nnR忽略重力，得 若密度为常数，则有此式为沿流线法向方向的伯努利方程，应用条件为（1）无粘性流体，（2）不可压流体（3）定常流（4）沿流线法向如果流线位直线时，曲率半径为无限大，则此式与静压力公式相同沿总流的伯努利方程将伯努利方程三项机械能在有效截面A上按质量流量积分，总机械能沿流束仍保持守恒，即以截面平均流速 V代替不均匀的速度分布，引入动能 修正因子。有考虑到质量守恒，得例题：Venturi对于一个缓变流的两个截面，有管：沿总流的伯努利方程一1 ;r , J-, . * _(gzi 色. g4gz)i : J %(gz 琢-m h- g马gz4)I j 1| jp,3g h g(

6、z4 Z3)p应用连续(Tm i)g h性方程伯努利方程的意义不可压缩粘性流体内流管道入口流动示意图，设管直径为d,管口外均流速度为 Uo从开始，流体在壁面上被滞止，形成边界层。边界层外仍保持为均流，称为核心流。由壁面不滑移条件引起壁面附近的流速降低，为满足质量守恒定律，核心流流速增大，速度廓线由平坦 逐渐变为凸出。随着边界层厚度不断增长，核心流不断加速，直至处四周的边界层相遇，核心流消失，整个管腔被边界层流动充满， 此后速度廓线不再变化。 称为入口段流动或发展中流动的速度廓线，均可通过求解 N-S方程 获得。入口段的压强损失，可利用动量方程求解。由例 B4.4.1D推导得管道入口段压强损失系

7、数为式中P0，Pl分别为x=0和x=L处的压强。称为达西摩擦因子，它是管道形状，雷诺数和管壁粗糙度的函数，在充分发展定常流动中为常数（将在C3.6中详细讨论）。（C3.2.1 ）式中的项为入口段中相应于充分发展段中的压强损失。K为入口段中特有的附加压强损失，它由两部分组成：将均流加速成充分发展流动所需要的压强系数； 入口段附加压强损失 K是入口段长度L,雷诺数Red及管道形状因子的函数，可运用有限差分法求解 N-S方程获得。根据计算的 K值可估算入口段的长度 L。圆管入口段长度与直径的比值的典型公式为-=0.06Red层流，人=4.4（Red）1/6湍流dd对层流，最大的入口段长度为LMax=

8、0.06 X 2300 x d=138d,（Re=2300）对湍流，由于边界层厚度增加较快，入口段长度比层流短的多Lt=（2040）d,（Re=10 4106）在实际的工程长管线中，如口段长度所占的比例往往是微不足道的，因此除特殊要求为其通常不予考虑，全长均按充分发展流动处理，但对一些较短的管道，则应该考虑入口段影响。平行平板间的层流流动是 N-S方程具有解析解的典型例子之一， 包括固定平板间的压差流， 平板间作相对平 移运动的剪切流及两种流动同时存在的一般库埃特流。分析库埃特流不仅有理论意义而且有工程背景，如气体或液体在活塞表面与缸壁间的缝隙中的泄漏流动，机床中滑块与导轨面的间隙中的润滑油流

9、动， 及滑动轴承的轴颈和轴承的间隙中润滑油流动等等。由于缝隙（b）很小，流动雷诺数不大，属于层流流态，均可用简化的无限大平行平面间的粘性流体定常层 流模型来分析。圆管湍流流动湍流尚无确切和全面的定义。湍流运动是有各种大小和不同涡量的涡旋叠加而形成的流动，在湍流中随机运动和拟序运动并存。圆管流动沿程水头形式的伯努利方程式推广形式为上式中称为水头损失，量纲是 Lo圆管流动中水头损失由两部分组成：沿程损失（hf）是沿等截面管流动时管壁粘性切应力引起的摩擦损失；局部损失（hm）,是由截面积变化，流动分离和二次流等局部因素引起的损失。达西公式在水平直圆管定常流动中只有沿程损失，因V1 =V2,z1 =

10、z2,由（Bb）式中可得hf = RP2 =P:g : g用量纲分析法求得实验表明Ap与l/d成正比关系，习惯上用 MRe,w/d b弋替f（Re/d ）。称为圆管沿程阻力因子或无量纲摩擦因子，因此上式可表为将上式代入（C3.6.1 ）式可得上式中l/d称为几何因子，V为管内平均速度，V 2 / 2g为速度水头。局部损失管入出口、管截面变化部位，弯头和三通、各种阀门等原因：（1）截面变化引起速度重新分布；（2）流体元相互碰撞和增加摩擦；（3）二次流；（4）流动分离成涡计算式:2gK为局部损失因子，V为制定部位的平均速度。入口 K=0.5 ,出口 K=1Vl和V2,求局部损失因子。管截面突然扩大

11、：局部损失，如图所示，平均速度分别是解：取图示虚线所示控制体CV由连续性方程实验证明在死水区（拐角分离区）的压强P =P,。忽略侧壁上的切应力作用，由动量方程由伯努力方程由（a）（b）（c）式（d）式中Vi为小管中的速度。含局部损失的管道损失当管道流动中局部损失在总损失中所占比重不能忽略时，管道计算中应将沿程损失和局部损失均考虑在内，全部损失为所有沿程损失和局部损失之和在可压缩流体流动中要考虑的流动参数除速度和压强外还要加上密度和温度。连续性方程不再独立， 必须与能量方程和状态方程联合求解， 求解的结果显然与不可压缩流体的流动规律不同。例如在一定条件下， 可压缩气体在截面积逐渐减小的收缩管道内

12、作减速流动，而在截面积逐渐增大的扩张管道内可作加速流动，均不违背质量守恒定律。通常用压强p、密度P和温度T三个物理量表示气体的状态，称为基本状态参数。 完全气体的基本状态参数满足如下方程式中R称为气体常数，由下式决定Rm称为通用气体常数，数值是 8314.3J/kg ?mol?K, m为气体平均分子量。气体的内能通常指分子热运动所具有的动能。完全气体的内能是温度的单值函数 其微分式可表示为Cv称为气体的比定容热容，它也仅是温度的函数Cv （T）单位质量气体的始称为比始，记为h(J/kg),定义为h=e+BPp在热力学中称为流动功，在流体力学中称为压能。烂是内能和压能之和。完全气体的燃的微分式可

13、表示P为Cp称为气体的比定压热容，它仅是温度的函数Cp(T)引入比热比与Cp,Cv的关系为如空气热力学第一定律表述为：对气体所加的热能等于气体内能的增加和气体对外所做功之和。表达式为对完全气体，可分别表示为单位质量气体的嫡称为比嫡定义为微分式为热力学第二定律表述为：气体在绝热的可逆过程中嫡值保持不变；在不可逆过程中嫡值必定增加。对完全气体，绝热而又可逆的过程称为等嫡过程，ds=0,气体作无摩擦绝热流动时为等嫡流动。对完全气体等嫡流动可得完全气体作等嫡流动时的状态参数关系式，常用表达式为-p-7 =常数pf声速来表示流体的可压缩性。声速是弹性介质中微弱扰动传播速度的总称。在图C5.2.2a中有一

14、竖向的微弱压强扰动波在精致的流体介质(V=0)中以声速c向左运动。设某瞬间波前的流体压强和密度分别为p, P。波后的流体速度变成 dV压强为p+dp,密度为P + dP。对地面上的观察者而言，这是一个非定常流动。为了便于考察波前波后流体状态参数的变化关系，在扰动层上取一薄层控制体CV,两边的面积均为 A,并将坐标系固定在控制体上与波一起前进，对站在坐标系上的观察者而言，流动是定常的。左边的流体压强为p,密度为P,以速度c流入控制体；然后以压强 p+dp,密度P + dP,及速度为c-dV流出控制体。由一维连续性方程(B4.2.11 )式有展开后得略去二阶小量项可得角标s表不等嫡过程。上式是完全

15、气体的理论声速公式当一个强烈的压缩扰动在超声速流场中传播时，在一定条件下将形成强压缩波阵面，称为?t波.激波是流动参数的强间断面，流体通过激波后流动参数发生突跃地升高，而速度则突跃地降低。截面变化与Ma关系设管内流动沿管轴 x方向，平均速度为 V。由一维定常流动欧拉运动方程式式中p为截面上的平均压强，可得由一维可压缩定常流动连续性方程对式取微分，并用 PVA除，可得利用（C5.4.2或由（C5.4.2）式从（C5.4.2Ma1时，Ma2 _1 1时，Ma2 -10,dA与dV同号，与dp异号说明（图C5.4.1）:亚声速流（Ma1）在收缩管（dA0）和减压（dpo）中，将减速dV0）和增压（dp0）与不可压缩流动相似。超声速流（Ma1）