传输原理 4 动量微分方程

冶金传输原理北京科技大学冶金与生态工程学院冶金传输原理第一部分流体力学第四章：动量传输的微分方程吴铿2009.03.15第四章流体动力学•4.1连续性微分方程•4.2理想流体运动方程•4.3伯努利方程•4.4伯努利方程的应用•4.5实际流体运动方程4.1连续性微分方程•微元控制体(流体的密度为ρ)dxxvABCEFGHDdydzzvyv流进ABCD面的质量流量流出EFGH面的质量流量xxxdydzdxdydzxX轴方向的净流入的质量(流入和流出之差)。xdxdydzxxyzoxdydz()xxdxdydzx4.1连续性微分方程•同理可以给出y和z两个方向上的表达式。•流体考虑三个方向上净流入质量为：yxzdxdydzxyz净流入量导致微元控制体内流体密度的变化，微元体内单位时间质量的累积(或累积的质量)()dxdydzdxdydzdxdydztt质量守衡定律：净流入质量=累积的质量(质量的累积)yxzdxdyddxdydzyztxz0yxzztxy由前面的方程可以得到连续性方程式中未涉及力的问题，是运动学方程，对理想流体和粘性流体均适用，实际存在的流体都必满足连续性方程xyzDDttxyz0yxzxyz()0yxzDDtxyz将上式展开，考虑可得：不可压缩流，ρ为常数4.1连续性微分方程0yxzztxy0yxzxyz对于稳定流动0t例题1：已知某二维不可压缩流的速度分布(或速度场)如下，试分析此流场是否存在？解：因为流场是否存在必须满足连续性方程，所以可用上面的直角坐标公式来判断。由于：26()xxy32yyz4zxyz对于园柱坐标，不可压缩流：10rzrrrzr4.1连续性微分方程例题2试判断下列平面流场是否连续？6xx2yy4zz624120yxzxyz22sincos,2cosrrr0zz2sincosrr4sincosr2sincosrrr12sincos4sincos2sincos0rrrrr解：因为代入前面的园柱坐标公式可得说明此流场是不连续的，亦即不存在。说明此流场是连续的。4.1连续性微分方程•取一个微分控制体如图所示。由于动量是矢量，首先考察x轴方向上的动量。这一方向上动量输入、输出控制体的方式有三种，即沿着x轴，以及沿着y轴和z轴的输入、输出。•沿x轴输入和输出的动量速率(流密度)分别为：4.2理想流体运动方程(欧拉方程),xxdydzxxxxdxdydzx净动量速率为xxdxdydzx4.2理想流体运动方程(欧拉方程)•同理，在沿着y轴和z轴方向上，x轴方向上的净动量速率分别为xydxdydzyxzdxdydzz•三个方向之和即为轴方向总的净动量速率：xyxxxzdxdydzxyz•设x轴方向上单位质量的质量力为，则整个微元控制体的x轴方向的质量力为：Xdxdydz4.2理想流体运动方程(欧拉方程)•在x轴方向上，控制体内累积的动量速率为xdxdydzt•将上述公式代入动量守恒定律中，于是有xxyxxxzPXxxyzt•展开上式右侧各项，考虑连续性方程，并假定密度不变，则有：xxxxxyzPXxtxyzD1DxPXxt4.2理想流体运动方程(欧拉方程)•关于y、z方向上的力与动量速率，也应该有类似的关系，它们的表达式如下：yyyyxyzPYytxyzzzzzxyzPZztxyzD1DyPYytD1DzPZzt•将前面的动量速率与力的分量式(共3个)写成矢量式，可得：1DDPtF•这就是理想流体运动方程，是1775年欧拉(Euler)首先提出的，亦称欧拉方程。4.3伯努利方程•因为只考虑定常流动，所以欧拉方程中的P、υx、υy和υz都只是坐标x、y、z的函数，而与时间t无关。也就是说：•因此在定常流动下，欧拉方法的速度表达式：0yxzttt,,,xxxyzt,,xxxyz,,xxxyzxxxxxxxxyzddxdydzdtxdtydtzdtxyz•而对时间的全微分并不为零，其表达式如下：4.3伯努利方程•同理：yyyxyzydytzdxzzzxyzzdytzdx•在流体力学中，从上面的等对时间的全微分不为零看出，流体力学中有了新的加速度概念，称为位置变化造成的加速度，简称位变加速度，而等不为零时对应的加速度称为时变加速度。,,xxxyzxt•不难看出，位变加速度是专属于流体力学的新概念。4.3伯努利方程•根据对于欧拉方程，考虑以下特殊条件：1.理想流体；2.稳定流动；3.不可压缩流体；4.质量力只有重力；5.质点沿一条特定流线运动。1xDPXxDtxxxxyzxytxz两边乘以dx1xxxxyzPXdxdxdxdxdxxxyz沿流线移动，流线微分方程式,,zxyxyzdydzddzxxdydxyzdxdydz4.3伯努利方程类似的三式相加2222xyz（）2()()2Pgdzdd21()()()2PPPXdxYdyZdzdxdydzdxyz22PzCg1()xxxxPXdxdxdxdydzxxyz2()2xxxdd21()2yPYdydydy21()2zPZdzdzdzPPPdPdxdydzxyz稳定流时4.3伯利例方程•理想流体沿流线的伯努利方程,其物理意义可以分为：①能量意义：方程中每一项表示单位重力流体所具有的能量。gzi和Pi/ρ分别代表单位重力流体所具有的位能和压力能；而υi2/2代表单位重力流体所具有的动能。•它说明理想流体沿流线流动时，单位重力的流体所具有的位能、压力能和动能三者之间尽管可以转化，但三者之和必为常数。这显然是机械能守恒定律的推广。2211221222PPgzgz2211221222PPzzgggg4.3伯利例方程②几何意义：方程中每一项的量纲与长度相同，表示单位重力流体所具有的水头。zi表示所研究点相对某一基准面的几何高度，称为位置水头；Pi/ρg表示所研究点处压强大小的高度，有与该压强相当的液柱高度，称为测压管高度，或称为测压管水头；υi2/2g表示所研究点处速度大小的高度，称为测速高度，或称为速度水头。伯努利方程说明在重力作用下的理想流体定常流动中，几何高度、测压管高度和测速高度之和为一个常数，称为水力高度和总水头。•柏努利方程都是在一定条件下积分得到的，应用时必须注意下列限制条件：理想流体；不可压缩；定常流动；作用于流体上的力仅有重力；不考虑流体旋转。4.3伯利例方程•意义：反映在重力作用下理想不可压缩流体稳定流动中，沿同一流线上，单位重量流体具有的位能、压能和动能的相互转换和守恒关系。•例题5在一均匀的平行于x方向的流动中，放置一个半径为R的静止圆柱体，已知在圆柱体表面处的流速υθ，试求其压强分布。解：假设距圆柱体很远处的压强为P∞，速度为υ∞，在忽略位能项的影响，有：22022PPP式中P0为滞止压强。4.3伯利例方程•以为在无旋流场中滞止压强是一常数，在圆柱体表面处的速度为：sin2•因此圆柱体表面处的压强是：2222(2sin)2222PPP22[14sin]2P2214sin2PP•图给出了绕圆柱体有势流动的压强分布曲线。4.3伯利例方程4.4伯努利方程的应用•伯努利方程应用——文丘里管•用途：文丘里管用于测量管路中的流速或流量。•原型：在管路中加接一段截面收缩的管子，并与压力计相连如右所示。选取沿管轴的一条流线及•流线上的两点1和2，对应的管子截面分别为A1、A2，流速为υ1、υ2。设γ1、γ2分别为管路中流体与U形管中流体的重度。若测量出压力计的高度差值，即可得到管路中的流速或流量值。•由流管的流量公式，在流量不变的情况下，可以得到：A1υ1=A2υ24.4伯努利方程的应用•根据重力场中不可压缩定常流动的伯努利方程得：2211221222PPgg•因为P1-P2＝γ2h，同时考虑到γ1γ2，故P1P2，上式因此变为：2222112hg222111212AhgA•由此得流速和流量：1212122.1ghAA211122212.11ghQAAA4.4伯努利方程的应用•伯努利方程应用——毕托管•用途：测量流场内某点流速的仪器。•原型:直角管两端开口，一端面向来流，另一端向上，管内液面高出水面H。•A端形成一驻点(速度为0)，驻点处的压力称为总压力。•B点在A点的上游，与A点位于同一水平流线，不受侧管影响。BAHH000()ABPHHPH4.4伯努利方程的应用2BAB2PPg2BAU2PHgB＝-PUB2gH式中：PA－总压，PB－静压动压，N/m2应用伯努利方程于A、B两点2BBA0002PPg4.4伯努利方程的应用例题6一毕托管安装在某烟道内，与毕托管连接的酒精压差计读数为h=5mm，酒精的相对密度为d=0.8，若烟气温度400℃时其重度为γg=5.13N/m3，求测点处烟气的流速。UB2gH30.8109.81229.810.00512.3/5.13UgHms4.4伯努利方程的应用例题7如图所示为一虹吸管，水从一个大容器经虹吸管流入大气中，若出口截面上流速均匀分布，试求出口处的流速及A点处流体的压强。设液面高度保持不变。解题思路：两个未知数，两个方程，伯努利和质量守恒。2112221222PPgzzgAA22114.5实际流体运动方程•理想流体运动的微分方程——欧拉方程•实际的流体有粘性——粘性运动微分方程增加了一个粘性项•TX、TY、TZ是作用在流体上的粘性力在x、y、z轴上的投影；•在下面来求三个坐标轴上粘性力的投影。DDDDDDxxyyzzPXTxtPYTytPZTztDDDDDDxyzPXxtPYytPZzt4.5实际流体运动方程•由牛顿内摩擦定律(或牛顿粘性定律)，考虑流体流动时流体的内摩擦力(又称粘性力)可