第三章-流体动力学基础

第三章流体动力学基础3-1描述流体运动的两种方法着眼点不同拉格朗日法（Lagrange）：流体质点欧拉法（Euler）：空间跟踪追迹法设立观察站法一、拉格朗日描述法与质点系（a,b,c）为t=t0起始时刻质点所在的空间位置坐标，称为拉格朗日变数。任何质点在空间的位置（x,y,z）都可看作是（a,b,c）和时间t的函数：或r＝r（a,b,c,t）（1）（a,b,c）=const,t为变数,可以得出某个指定质点在任意时刻所处的位置。（2）（a,b,c）为变数,t=const,可以得出某一瞬间不同质点在空间的分布情况。流体质点任一物理量B（如速度、压力、密度等）表示为：B＝B（a,b,c,t）质点系：在t＝0时紧密毗邻的具有不同起始坐标（a,b,c)的无数质点组成一个有确定形状、有确定流动参数的质点系。经过t时间之后，质点系的位置和形状发生变化。二、欧拉描述法与控制体欧拉法不直接追究质点的运动过程，而是以充满运动流体质点的空间——流场为对象。流体质点的物理量B是时空（x,y,z,t）的连续函数：B＝B（x,y,z,t）（x,y,z,）——欧拉变量速度场：u＝u(x,y,z,t),v＝v(x,y,z,t),w＝w(x,y,z,t).控制体：将孤立点上的观察站扩大为一个有适当规模的连续区域。控制体相对于坐标系固定位置，有任意确定的形状，不随时间变化。控制体的表面为控制面，控制面上有流体进出。三、两种描述方法之间的联系如果标号参数为(a,b,c)的流体质点，在t时刻正好到达(x,y,z)这个空间点上，则有B＝B(x,y,z,t)＝B(x(a,b,c,t),y(a,b,c,t),z(a,b,c,t),t)＝B(a,b,c,t)3-2流体运动的几个基本概念一、物理量的质点导数质点导数定义：流体质点的物理量随时间的变化率。随体导数如速度V和加速度a为21、拉格朗日描述中的随体导数V和a在直角坐标系中展开：和以速度在直角坐标系为例：流体质点运动速度在欧拉法中，V＝V(x,y,z,t),由于位置又是时间t的函数，所以流速是t的复合函数，对流速求导可得加速度：写成分量形式2、欧拉描述中随体导数用哈密顿算子表示：局部（当地）加速度：同一空间点上流体速度随时间的变化率。定常流动该项为0。迁移（位变）加速度：同一时刻由于不同空间点的流体速度差异而产生的速度变化率。均匀流场该项为0。对于任一物理量B：局部（当地）导数，表示流场的非定常性。迁移（位变）导数，表示流场的均匀性。质点导数例题：解：二、定常流与非定常流（或恒定流与非恒定流）三、均匀流与非均匀流四、一元流、二元流与三元流按流体运动要素所含空间坐标变量的个数分：（1）一元流一元流(one-dimensionalflow):流体在一个方向流动最为显著，其余两个方向的流动可忽略不计，即流动流体的运动要素是一个空间坐标的函数。若考虑流道（管道或渠道）中实际液体运动要素的断面平均值，则运动要素只是曲线坐标s的函数，这种流动属于一元流动。（2）二元流二元流(two-dimensionalflow):流体主要表现在两个方向的流动，而第三个方向的流动可忽略不计，即流动流体的运动要素是二个空间坐标（不限于直角坐标）函数。（3）三元流三元流（three-dimensionalflow):流动流体的运动要素是三个空间坐标函数。五、迹线与流线迹线流体质点在流场中的运动轨迹线。是拉格朗日法描述流体运动的基础。1、迹线流线是流场中这样一条曲线，曲线上任一点的切线方向与该点的流速方向重合。流线是欧拉法描述流体运动的基础。图为流线谱中显示的流线形状。2、流线流线的作法：在流场中任取一点，绘出某时刻通过该点的流体质点的流速矢量u1，再画出距1点很近的2点在同一时刻通过该处的流体质点的流速矢量u2…，如此继续下去，得一折线1234…，若各点无限接近，其极限就是某时刻的流线。流线方程:设dr为流线上A处的一微元弧长矢量:V为流体质点在A点的流速:根据流线的定义，可以求得流线的微分方程：展开后得到：——流线微分方程dr流线的性质：1.在某一时刻，过某一空间点只有一条流线。流线不能相交，不能突然转折。三种例外：2.对于非定常流动，流线具有瞬时性。3.一般情况下，流线迹线不重合。定常流动中流线形状不随时间变化，而且流体质点的迹线和流线重合驻点相切点奇点脉线在一段时间内，会有不同的流体质点相继经过同一空间固定点，在某一瞬时将这些质点所处的位置点光滑连接而成的曲线。流线、迹线和脉线是本质不同的三种描述流体运动的线，定常时互相重合。六、流管与流束1.流面在流场中作一条任意的空间曲线L（非流线），过此曲线的每一点作流线，这些无数密集的流线所构成的曲面。性质：（与流线相似）（1）在某一时刻，过一条曲线只有一个流面；（2）非定常时，流面形状随时间变化；（3）流体不能穿越流面。2.流管与流束流管定义流管性质：（1）不能相交；（2）形状和位置在非定常时随时间变化；（3）不能在流场内部中断，只能始于或终于流场的边界。如物面，自由面等。流束除了有流管的性质以外，还具有：（1）截面上的速度处处相等；（2）微小截面看成是平面。流束定义：截面面积很小的流管，微元流管。流束的极限是流线。流管截面：以L为周界可以作很多的面，可以是平面或曲面。有效截面（过流断面）：截面上的流速方向处处与该面垂直缓变流动：如果微小流束（流线）间的夹角及流束的曲率都非常小，这种流动称为缓变流动。反之急变流。缓变流的过流断面可看作是平面。急变流的过流断面是曲面缓变流七、流量、净通量1、流量单位时间内通过某一过流断面的流体量。体积流量qv或Q表示，质量流量qm。体积流量（m3/s）：质量流量（kg/s）：如果dA不是过流断面，而是与微元流束相交的任意断面，则体积流量（m3/s）：质量流量（kg/s）：2、净通量流过全部封闭控制面A的流量称为净流量，或净通量。AmAvdAnvρqdAnvqAvρvdAρqAvvdAqAmAvAvdAnvq八、过流断面上的平均速度与动能动量修正系数1、断面平均速度过流断面上各点的流速是不相同的，所以常采用一个平均值来代替各点的实际流速，称断面平均流速。2、动能及动能修正系数动能（kineticenergy）：是指物体由于机械运动而具有的能量。单位时间内通过过流断面的流体动能是：动能修正系数——是实际动能与按断面平均流速计算的动能的比值。AvdAAqvAvdAρvvdmEAk32212113121212233dAvAvAvρdAρvαAA注意：动能修正系数是无量纲数，它的大小取决于总流过水断面上的流速分布，分布越均匀，α值越小，越接近于1.0。层流流速分布湍流流速分布2、动量及动量修正系数动量（momentum)是物体运动的一种量度，是描述物体机械运动状态的一个重要物理量。单位时间内通过过流断面的流体动量是：动量修正系数——是实际动量与按断面平均流速计算的动量的比值。动量修正系数是无量纲数，它的大小取决于总流过水断面的流速分布，分布越均匀，β值越小，越接近于1.0。AdAρvdmvK21112222dAvAvAvρdAρvβAA断面流速分布动能修正系数动量修正系数圆管层流旋转抛物面=2.0β=4/3圆管紊流对数规律=1.05~1.1β=1.02~1.05层流流速分布湍流流速分布§3-3连续方程式一、基本原理特例0AdAnv0VρdVt特例1定常流动则0AdAnvρ特例2不可压缩流动为常数则0VρdVt流管流动的连续性方程的应用：恒定流动时：对于不可压缩流体，则222111AvρAvρ2211AvAv连续性方程的积分形式：由奥－高公式根据控制体与时间的无关性直角坐标系下连续性方程的微分形式即想一想：恒定、不可压情况下，连续性方程的微分形式。二、连续性方程的微分形式0VAρdVtdAnvρdVvρdAnvρVA)(VVdVtρρdVt0)(vρtρ0)()()(zρvyρvxρvtρzyx§3-4流体微团的运动分析一、流体与刚体比较刚体的运动是由平移和绕某瞬时轴的转动两部分组成。流体质点的运动，一般除了平移、转动外，还要发生变形（角变形和线变形）。二、流体微元的速度分解A(x,y,z)点速度为vx,vy,vz，则C点的速度为：xyεdtdαzωdtdβ三、有旋流和无旋流根据流体微团是否绕自身轴旋转，可分为有旋流和无旋流。1.定义：有旋流（vortex）：亦称“涡流”。流体质点（微团）在运动中不仅发生平动（或形变），而且绕着自身的瞬时轴线作旋转运动。如旋风即为空气的涡流。当流体速度变化较大，由于流体粘滞阻力、压强不均匀等因素的影响，就容易形成涡流。无旋流（potentialflow）亦称“势流”、“有势流”。流体在运动中，它的微小单元只有平动或变形，但不发生旋转运动，即流体质点不绕其自身任意轴转动。注意：无旋流和有旋流决定于流体质点本身是否旋转，而与运动轨迹无关。2.有旋流和无旋流的特性（1）若wx=wy=wz=0，即则流动为无旋流，否则，为有旋流。有旋流（涡流）——wx、wy、wz中任一个或全部不等于零的流体运动，绕自身轴有旋转的运动。（与通常的旋转不同）流场内流体质点具有绕质点自身任意轴的角速度。（2）有旋流的特征是存在角速度。角速度是一个矢量，所以可如同用流线描述流动一样，可用涡线描述流动的旋转变化。涡线——在同一瞬时线上各质点的转速矢量都与该曲线相切。无旋流一般存在于无粘性理想流体中。有旋流一般存在于有粘性实际流体中。0,0,0yvxvzvxvzvyvxyxzyz例题已知流体流动的流速场为，判断该流动是无旋流还是有旋流？解：故液体流动是无旋流。§3-5实际流体的运动微分方程式一、作用在流体微元上的应力应力矩阵二、本构方程确定应力与应变的方程式叫本构方程。其中p：在平衡流体，代表一点上的流体静压强；在理想流体，代表一点上的流体动压强；在不可压实际流体，代表一点上的流体动压强的算术平均值。三、纳维－斯托克斯方程式不可压实际流体的运动方程式——N-S方程想一想理想流体、静止情况下的方程。§3-6伯努利方程式及其应用一、流线上的伯努利方程式假设单位质量的流体质点某瞬时的速度为v＝vxi+vyj+vzk，经dt时间，质点沿流线移动一段微小距离ds＝dxi+dyj+dzk＝vxdti+vydtj+vzdtk，为求出单位质量流体移动ds距离与外力作功的能量关系，将ds的三个投影分别与N-S方程的三个式子相乘，然后相加，得下面分别对式中的四类项进行简化1.质量力项，假设质量力有势2.压强项3.粘性摩擦力项4.导数项将结果代回原式，则可得则——适用范围：非定常、质量力有势。——适用范围：定常、质量力有势。——适用范围：定常、重力场、不可压流体。——适用范围：理想、定常、重力场、不可压流体。那么，实际流体在定常、重力场、不可压条件下，在流线上任意两点间可列出伯努利方程为：理想流体在相同条件下，在流线上任意两点间的伯努利方程为：二、粘性总流的伯努利方程式粘性流体在定常、重力场、不可压条件下，在流线上任意两点间可列出伯努利方程为其中用代替，则fhfhgvgpzgvgpz2222222111在实际工程中，我们遇到的往往是过流断面具有有限大小的流动，我们称它们为总流。因此我们应将沿流线的伯努利方程推广到沿总流上去。将上式乘以gdqv，然后对整个总流断面积分，这样就获得总流的能量关系式AvfvAvAρgdqhρgdq)gvzρgp(ρgdq)gvzρgp(2122222221111)为单位时间内通过断面A的势能总和。vAρgdqz)ρgp(假设两个过流断面上的流动为缓变流动,在缓变流动情况下，过流断面可以近似地认为是一个平面。由于过流断面是与流