伯努利方程的应用.doc

，伯努利方程及其应用伯努利，1738，瑞士。动能与压强势能相互转换。沿流线的伯努利方程将牛顿第二定律应用于控制体内的流体元，沿流线切线方向整理后因为将流体元的加速度转换成欧拉形式的加速度，沿流线的质点导数为则导出此式为一维欧拉方程，使用下述关系将方程沿流线积分。两边乘以ds得：沿流线积分此式为欧拉方程的积分式，适合于可压、无粘不定常运动。对于不可压定常流动，则可简化为此式为伯努利方程，三项分别表示单位质量流体具有的动能、位置势能和压强势能。即总机械能守恒。应用伯努利方程时常采用沿流线上任两点的总机械能值相等的形式。伯努利方程使用的限制条件（1）无粘性流体，（2）不可压流体（3）定常流（4）沿流线。加入能量损失就可适应粘性流体。皮托（pitot）测速管：总压强与动压强皮托测速管又称为皮托-静压管，简称皮托管，为纪念法国人皮托命名。皮托测速管由粗细两根同轴的圆管组成，细管（直径约为1.5 mm）前端开孔（O点），粗管（直径约为mm）在距前端适当长距离处的侧壁上开数个小孔（B点），在孔后足够长距离处两管弯成柄状测速时管轴线沿来流方向放置设正前方的流速保持为，静压强为，流体密度为。粗细两管中的压强被引入形测压计中，形管中液体密度。试求用形管液位差表示流速的关系式。解：设流动符合不可压缩无粘性流体定常流动条件。从皮托管正前方点到端点再到侧壁孔点的线是一条流线，点的速度和压强分别为和，沿流线段按（）式列伯努利方程+在皮托管端点，流体速度降至零，称为驻点，称为驻点压强，形管右支管测到的即是驻点压强由于，由（a）式可得上式中称为动压强，为流体质点的动能全部转化为压强势能时应具有的压强（b） 式表明驻点压强为静压强和动压强之和，故称为总压强由（b）式动压强可表为 由于皮托管较细，流线上的两点的位置差可忽略，伯努利方程为因，由上式，即形管左支通过皮托管侧壁小孔测到的是当地静压强在形管内列压强关系式可得由于实际流体具有粘性及皮托管加工误差等原因，流体动压强转化为形管内液位差读数存在误差，需乘上一个修正系数,由（c）,（d）式可得称为皮托管系数，可通过用标准皮托管作标定测量后确定由（e）式可得小孔出流：托里拆利公式及缩颈效应从一个大容器侧壁下部距离液面为h处开一个小孔，设液面水位不变，求出流速度和出流流量。解：=沿流线法向方向的速度压强关系式由牛顿第二定律：得考虑到几何关系，有整理，得忽略重力，得若密度为常数，则有此式为沿流线法向方向的伯努利方程，应用条件为（1）无粘性流体，（2）不可压流体（3）定常流（4）沿流线法向。如果流线位直线时，曲率半径为无限大，则此式与静压力公式相同。沿总流的伯努利方程将伯努利方程三项机械能在有效截面A上按质量流量积分，总机械能沿流束仍保持守恒，即以截面平均流速V代替不均匀的速度分布，引入动能修正因子。有考虑到质量守恒，得对于一个缓变流的两个截面，有例题：Venturi管：沿总流的伯努利方程应用连续性方程伯努利方程的意义不可压缩粘性流体内流管道入口流动示意图，设管直径为d,管口外均流速度为U。从 开始，流体在壁面上被滞止，形成边界层。边界层外仍保持为均流，称为核心流。由壁面不滑移条件引起壁面附近的流速降低，为满足质量守恒定律，核心流流速增大，速度廓线由平坦逐渐变为凸出。随着边界层厚度不断增长，核心流不断加速，直至 处四周的边界层相遇，核心流消失，整个管腔被边界层流动充满，此后速度廓线不再变化。称 为入口段流动或发展中流动的速度廓线，均可通过求解N-S方程获得。入口段的压强损失，可利用动量方程求解。由例B4.4.1D推导得管道入口段压强损失系数为 式中p0,pL 分别为x=0和x=L 处的压强。 称为达西摩擦因子，它是管道形状，雷诺数 和管壁粗糙度的函数，在充分发展定常流动中 为常数（将在C3.6中详细讨论）。（C3.2.1）式中的项为入口段中相应于充分发展段中的压强损失。K为入口段中特有的附加压强损失，它由两部分组成：将均流加速成充分发展流动所需要的压强系数；由于入口段壁面速度廓线陡峭，壁面切应力大于充分发展段的壁面切应力，为克服这部分阻力差值所需要增加的压强系数。入口段和充分发展段的压强损失系数曲线如图C3.2.1b(分别为实线和虚线)所示。入口段附加压强损失K是入口段长度L，雷诺数及管道形状因子的函数，可运用有限差分法求解N-S方程获得。根据计算的K值可估算入口段的长度L。圆管入口段长度与直径的比值的典型公式为层流, 湍流对层流，最大的入口段长度为LMax=0.062300d=138d, (Re=2300)对湍流，由于边界层厚度增加较快，入口段长度比层流短的多Lt= (2040)d, (Re=104106)在实际的工程长管线中，如口段长度所占的比例往往是微不足道的，因此除特殊要求为其通常不予考虑，全长均按充分发展流动处理，但对一些较短的管道，则应该考虑入口段影响。平行平板间的层流流动是N-S方程具有解析解的典型例子之一，包括固定平板间的压差流，平板间作相对平移运动的剪切流及两种流动同时存在的一般库埃特流。分析库埃特流不仅有理论意义而且有工程背景，如气体或液体在活塞表面与缸壁间的缝隙中的泄漏流动，机床中滑块与导轨面的间隙中的润滑油流动，及滑动轴承的轴颈和轴承的间隙中润滑油流动等等。由于缝隙（b）很小，流动雷诺数不大，属于层流流态，均可用简化的无限大平行平面间的粘性流体定常层流模型来分析。 圆管湍流流动湍流尚无确切和全面的定义。湍流运动是有各种大小和不同涡量的涡旋叠加而形成的流动，在湍流中随机运动和拟序运动并存。圆管流动沿程水头形式的伯努利方程式推广形式为上式中 称为水头损失，量纲是L。圆管流动中水头损失由两部分组成：（1） 沿程损失（）是沿等截面管流动时管壁粘性切应力引起的摩擦损失；（2） 局部损失(), 是由截面积变化，流动分离和二次流等局部因素引起的损失。达西公式在水平直圆管定常流动中只有沿程损失，因，由（B4.6.13b）式中可得用量纲分析法求得 实验表明与成正比关系，习惯上用代替。称为圆管沿程阻力因子或无量纲摩擦因子，因此上式可表为将上式代入（C3.6.1）式可得上式中称为几何因子，V为管内平均速度，为速度水头。局部损失管入出口、管截面变化部位，弯头和三通、各种阀门等。原因：（1）截面变化引起速度重新分布；（2）流体元相互碰撞和增加摩擦；（3）二次流；（4）流动分离成涡。计算式：K为局部损失因子，V为制定部位的平均速度。入口K=0.5，出口K=1管截面突然扩大：局部损失，如图所示，平均速度分别是V1和V2，求局部损失因子。解：取图示虚线所示控制体CV，由连续性方程实验证明在死水区（拐角分离区）的压强 。忽略侧壁上的切应力作用，由动量方程由伯努力方程由(a)(b)(c)式 （d）式中为小管中的速度。含局部损失的管道损失当管道流动中局部损失在总损失中所占比重不能忽略时，管道计算中应将沿程损失和局部损失均考虑在内，全部损失为所有沿程损失和局部损失之和在可压缩流体流动中要考虑的流动参数除速度和压强外还要加上密度和温度。连续性方程不再独立，必须与能量方程和状态方程联合求解，求解的结果显然与不可压缩流体的流动规律不同。例如在一定条件下，可压缩气体在截面积逐渐减小的收缩管道内作减速流动，而在截面积逐渐增大的扩张管道内可作加速流动，均不违背质量守恒定律。通常用压强p、密度和温度T三个物理量表示气体的状态，称为基本状态参数。完全气体的基本状态参数满足如下方程式中R称为气体常数，由下式决定称为通用气体常数，数值是8 314.3J/kgmolK,为气体平均分子量。气体的内能通常指分子热运动所具有的动能。完全气体的内能是温度的单值函数其微分式可表示为称为气体的比定容热容，它也仅是温度的函数单位质量气体的焓称为比焓，记为h(J/kg)，定义为h=e+在热力学中称为流动功，在流体力学中称为压能。焓是内能和压能之和。完全气体的焓的微分式可表示为称为气体的比定压热容，它仅是温度的函数引入比热比与的关系为如空气热力学第一定律表述为：对气体所加的热能等于气体内能的增加和气体对外所做功之和。表达式为对完全气体，可分别表示为单位质量气体的熵称为比熵定义为微分式为热力学第二定律表述为：气体在绝热的可逆过程中熵值保持不变；在不可逆过程中熵值必定增加。对完全气体，绝热而又可逆的过程称为等熵过程，ds=0,气体作无摩擦绝热流动时为等熵流动。对完全气体等熵流动可得完全气体作等熵流动时的状态参数关系式，常用表达式为常数声速来表示流体的可压缩性。声速是弹性介质中微弱扰动传播速度的总称。在图C5.2.2a中有一竖向的微弱压强扰动波在精致的流体介质（V=0）中以声速c向左运动。设某瞬间波前的流体压强和密度分别为p,。波后的流体速度变成dV压强为p+dp,密度为。对地面上的观察者而言，这是一个非定常流动。为了便于考察波前波后流体状态参数的变化关系，在扰动层上取一薄层控制体CV，两边的面积均为A，并将坐标系固定在控制体上与波一起前进，对站在坐标系上的观察者而言，流动是定常的。左边的流体压强为p,密度为，以速度c流入控制体；然后以压强p+dp,密度，及速度为c-dV流出控制体。由一维连续性方程（B4.2.11）式有展开后得略去二阶小量项可得角标s 表示等熵过程。上式是完全气体的理论声速公式当一个强烈的压缩扰动在超声速流场中传播时,在一定条件下将形成强压缩波阵面,称为激波.激波是流动参数的强间断面,流体通过激波后流动参数发生突跃地升高,而速度则突跃地降低。截面变化与Ma关系设管内流动沿管轴x方向，平均速度为V。由一维定常流动欧拉运动方程式式中p为截面上的平均压强，可得由一维可压缩定常流动连续性方程对式取微分，并用除，可得利用（C5.4.2）式及声速公式,由(C5.4.3)式可得因Ma=V/c,上式可化为或由(C5.4.2)式从(C5.4.2)和(C5.4.5)式可看到:Ma1时,号说明(图C5.4.1):亚声速流(Ma1)在收缩管(dA0)和减压（dpo)中，将减