1、第二章第二章 流体运动学和动力学基础 本章作业:习题本章作业:习题1,3,6,8,9,10,11,13,15 流场流场(流场及其描述方法流场及其描述方法,迹线、流线和流管迹线、流线和流管)流体微团运动的分析流体微团运动的分析(散度、旋度和速度位散度、旋度和速度位)连续方程和流函数连续方程和流函数(连续方程、流函数连续方程、流函数)旋涡运动旋涡运动(涡线、涡管及旋涡强度、环量诱导速度及相关定理涡线、涡管及旋涡强度、环量诱导速度及相关定理)欧拉运动方程及其积分欧拉运动方程及其积分(欧拉运动方程、伯努利方程欧拉运动方程、伯努利方程)流体力学中的动量定理流体力学中的动量定理(一般原理及例子一般原理及例
2、子)2.1.1 拉格朗日方法与欧拉方法拉格朗日方法与欧拉方法1 1、LagrangeLagrange方法(方法(拉格朗日方法,拉格朗日方法,质点法)质点法)着眼于流场中每一个运动着的流体质点,跟踪观察每一个流体质点的运动轨着眼于流场中每一个运动着的流体质点,跟踪观察每一个流体质点的运动轨迹以及运动参数(速度、压强、加速度等)随时间的变化,然后综合所有流体迹以及运动参数(速度、压强、加速度等)随时间的变化,然后综合所有流体质点的运动,得到整个流场的运动规律。质点的运动,得到整个流场的运动规律。2、Euler方法方法 观察者相对于坐标系是固定不动的,着眼于不同流体质点通过空间固定点的观察者相对于坐
3、标系是固定不动的,着眼于不同流体质点通过空间固定点的流动行为,通过记录不同空间点流体质点经过的运动情况,从而获得整个流场流动行为,通过记录不同空间点流体质点经过的运动情况,从而获得整个流场的运动规律。的运动规律。一个速度场2022-10-262沈阳航空工业学院飞行器设计教研室2022-10-263沈阳航空工业学院飞行器设计教研室右边第右边第1项:项:表示速度对时间的偏导数,是由流场的非定常性引起的,称为局部表示速度对时间的偏导数,是由流场的非定常性引起的,称为局部加速度,或当地加速度;加速度,或当地加速度;右边其他项:右边其他项:表示因流体质点位置迁移引起的加速度,称为迁移加速度,位变加表示因
4、流体质点位置迁移引起的加速度,称为迁移加速度,位变加速度,或对流加速度。速度,或对流加速度。二者的合成称为全加速度,或随体加速度。二者的合成称为全加速度,或随体加速度。zwwywvxwutwdtdwzvwyvvxvutvdtdvzuwyuvxuutudtduxxxxxxyzyyyyyxyzzzzzzxyzvvvvavvvtxyzvvvvavvvtxyzvvvvavvvtxyz加速度描述加速度描述2022-10-264沈阳航空工业学院飞行器设计教研室算子表示随流体质点运动的导数,称随体导数。除速度外,对流场中其它变量也成立。如对于压强p,有zwyvxutdtdzpwypvxputpdtdp推广推
5、广2022-10-265沈阳航空工业学院飞行器设计教研室流线流线:流场中的瞬时光滑曲线,在曲线上流体质点的速度方向与各该点流场中的瞬时光滑曲线,在曲线上流体质点的速度方向与各该点的切线方向重合。的切线方向重合。迹线:迹线:流体质点在一定时间内所经过的所有空间点的集合。流体质点在一定时间内所经过的所有空间点的集合。场、定常与非定常场、定常与非定常流管、流面、流量:流管、流面、流量:vsdAnVmA)(dAnVQA)(dAnVgGA)(流量流量是单位时间内穿过指定截面的流体量(体积、质量或重量),例如穿过上述流管中任意截面A的体积流量 、质量流量 和重量流量 可分别表为QmG其中,是局部速度向量,
6、是密度,是微元面积 的法线向量VndAxyzdxdydzvvvVdswdzvdyudx或或2022-10-267沈阳航空工业学院飞行器设计教研室流体微团平动速度:流体微团平动速度:流体微团线变形速率:流体微团线变形速率:流体微团角变形速率(剪切变形速率):流体微团角变形速率(剪切变形速率):流体微团旋转角速度:流体微团旋转角速度:zwyvxuzyx ,yuxvxwzuzvywzyx21,21,21yuxvxwzuzvywzyx21,21,21),(),(),(tzyxwtzyxvtzyxu 2.2.1 流体微团的基本运动形式流体微团的基本运动形式2022-10-268沈阳航空工业学院飞行器设计
7、教研室按速度泰勒级数展开有按速度泰勒级数展开有 zyxxytzyxwzzwyywxxwtzyxwtzzyyxxwzyxzxtzyxvzzvyyvxxvtzyxvtzzyyxxvzyxyztzyxuzzuyyuxxutzyxutzzyyxxuzxyyxxyzxzyzxzy)(),(),(),()(),(),(),()(),(),(),(2.2.2 流体微团速度分解定理流体微团速度分解定理2022-10-269沈阳航空工业学院飞行器设计教研室 2.2.3 散度及其意义散度及其意义散度在流体力学里表示流体微团的相对体积膨胀率相对体积膨胀率(单位时间单位体积的增长量)。三个相互垂直方向的线变形率之和在
8、向量分析中称为速度V的散度,散度,符号为divV,即zwyvxuVVdiv在密度不变的不可压流动里,微团的体积不变,其速度的散度必为零。2022-10-2610沈阳航空工业学院飞行器设计教研室一个流场,如果各处的都等于零,这样的流场称为无旋流场,其流动称为无无旋流旋流。否则为有旋流场,其流动称有旋流有旋流。根据数学上Stokes定律 2.2.4 旋度和位函数旋度和位函数流体微团绕自身轴的旋转角速度的三个分量为x,y,x,合角速度可用矢量表示为这个值在向量分析里记为(1/2)rotV,称为V的旋度。旋度。即有即有旋度为旋转角速度的二倍:旋度为旋转角速度的二倍:VVrotkjizyx2121ALA
9、dVrotrdV如果是无涡流场,那么其旋度为零,由此得到 说明速度场的曲线积分与路径无关,仅是坐标位置的函数。0LrdVyyxxzzvvvvvvcurlyzzxxyi+j+k。22022-10-2611沈阳航空工业学院飞行器设计教研室;2.2.4 旋度和位函数旋度和位函数在数学上表示下列微分代表某个函数的全微分,即L0 LdrdVwdzvdyudxrdVd上式中这个函数称为速度势函数或速度位,其存在的充分必要条件是无涡流动。速度势函数仅是坐标位置和时间的函数。即速度势函数与速度分量的关系为说明速度势函数在某个方向的偏导数等于速度矢量在那个方向的分量。021Vrot),(tzyxzwyvxu x
10、-y2022-10-2612沈阳航空工业学院飞行器设计教研室微分形式的连续方程:000)(0)()()(VdtdzwyvxuzwyvxutVtzwyvxut对于不可压缩流体,连续方程变为0 xu 0 0zwyvVdtd 2.3.1 连续方程连续方程2022-10-2613沈阳航空工业学院飞行器设计教研室结论:流体在运动时,应服从质量守恒定律,这条定理在空气动力学中的流体在运动时,应服从质量守恒定律,这条定理在空气动力学中的数学表达式称为连续方程或质量方程数学表达式称为连续方程或质量方程.矢量表达形式矢量表达形式()0divt对于定常不可压流体的极坐标方程对于定常不可压流体的极坐标方程另一形式另
11、一形式()0yxzvvvddtxyz1()0rzrvvvrrz2022-10-2614沈阳航空工业学院飞行器设计教研室2.3.2 流函数流函数又又故有故有xydv dyv dxddydxyx,xyvvyx(2-23)(2-24)(2-25)速度位与流函数关系:等位线族与流线族正交速度位与流函数关系:等位线族与流线族正交zwyvxu x-y(x,y)称为流函数称为流函数2022-10-2615沈阳航空工业学院飞行器设计教研室欧拉方程微分形式(牛顿第二定律的描述)(教材更直观):上三式即为笛卡儿坐标系下理想流体运动的欧拉方程(1755年)。表明了流体质点的加速度等于质量力减去压力梯度。写成另一种形
12、式,为 2.3.2 Euler运动微分方程组运动微分方程组zpfdtdwypfdtdvxpfdtduzyx111pfdtVd1zpfzwwywvxwutwypfzvwyvvxvutvxpfzuwyuvxuutuzyx1112022-10-2616沈阳航空工业学院飞行器设计教研室如果上式右边项为零,有0sdV)(2 0222sCVVs这样在曲线上,下式成立。这就是Bernoulli积分,或伯努利方程。上式表明,对于理想正压流体的定常上式表明,对于理想正压流体的定常流动,在质量力有势条件下,单位体积流体微团沿着这条特定曲线流动,在质量力有势条件下,单位体积流体微团沿着这条特定曲线s的势能、的势能、
13、压能和动能之和不变,即总机械能不变。(压能和动能之和不变,即总机械能不变。(1738年)年)2.3.3 Bernoulli积分方程及其物理意义积分方程及其物理意义0sdVBernoulli积分成立的条件,是(1)沿着任意一条流线,Bernoulli积分成立。这是因为,在此情况下0sdV VV /Vsd2022-10-2617沈阳航空工业学院飞行器设计教研室 现任取一体积,边界表面积为现任取一体积,边界表面积为S0的确定系统作为考察对象。的确定系统作为考察对象。(1 1)连续方程(质量守恒)连续方程(质量守恒)表示,在系统内不存在源和汇的情况下,系统的质量不随时间变化。表示,在系统内不存在源和汇
14、的情况下,系统的质量不随时间变化。(2)动量方程动量方程 表示:系统的动量对时间的变化率等于外界作用于系统上的所有外力的合表示:系统的动量对时间的变化率等于外界作用于系统上的所有外力的合力。力。(3 3)动量矩方程)动量矩方程 表示:系统对某点的动量矩对时间的变化率等于外界作用于系统上所有外表示:系统对某点的动量矩对时间的变化率等于外界作用于系统上所有外力对同一点力矩之和。力对同一点力矩之和。2.4.2 Lagrange型积分方程型积分方程00ddtddtdM00000SndSpdfFdVdtddtKd000000)(dSprdfrFrdVrdtddtMdSnr2022-10-2618沈阳航空
15、工业学院飞行器设计教研室SnzzSSnyySSnxxSdSVwwdtFdfdSznpdSVvvdtFdfdSynpdSVuudtFdfdSxnp),cos(),cos(),cos(由动量积分方程,可得积分形式动量方程的一个重要方面在于人们往往不需要知道控制体中的流动细节,只需要知道控制面边界处的流动属性来求作用力,这个作用力可以包含摩擦力的影响在内,例如用上述方程来求物体受到的阻力等。上述方程常常用于定常流动的气体中,用于定常流时上式中的当地变化率一项等于零,用于气体则质量力可以忽略。2022-10-2619沈阳航空工业学院飞行器设计教研室2.6 旋涡运动旋涡运动-要点要点2.6.1 涡线、涡
16、管及旋涡强度涡线、涡管及旋涡强度涡线式充满运动流体的旋涡场中的一系列曲线,它具有这样的性质:涡线式充满运动流体的旋涡场中的一系列曲线,它具有这样的性质:在某瞬时该曲线上微团的旋转角速度向量(旋转轴线方向按右手定则)在某瞬时该曲线上微团的旋转角速度向量(旋转轴线方向按右手定则)都和曲线相切,如图都和曲线相切,如图其微分方程为其微分方程为xyzdxdydz2022-10-2620沈阳航空工业学院飞行器设计教研室2.6.3 直线涡的诱导速度及毕奥萨瓦定律直线涡的诱导速度及毕奥萨瓦定律我们把流场中由旋涡存在而产生的速度称为诱导速度。其大小可由毕奥萨我们把流场中由旋涡存在而产生的速度称为诱导速度。其大小
17、可由毕奥萨瓦公式来确定。瓦公式来确定。在不可压流动中,其数学表达式为在不可压流动中,其数学表达式为34dLrdwr或或2sin4dwdLr式中式中dL为涡线上的微段长度;为涡线上的微段长度;r 为流场中任意点至微段的距离;为流场中任意点至微段的距离;为微段为微段dL与与r之间的夹角;之间的夹角;为旋涡强度。为旋涡强度。d的方向垂直于的方向垂直于ONM平面,见图平面,见图2.6 旋涡运动旋涡运动-要点要点2022-10-2621沈阳航空工业学院飞行器设计教研室2.6.4 海姆霍兹旋涡定理海姆霍兹旋涡定理定理一 在同一瞬间沿旋涡线或涡管的旋涡强度不变。定理二 涡线不能在流线中中断;只能在流体边界上
18、中断或形成定理三 在理想流中,涡的强度不随时间变化,既不会增强也不会削弱或消失。闭合圈。2.6 旋涡运动旋涡运动-要点要点2022-10-2622沈阳航空工业学院飞行器设计教研室本章基本要求本章基本要求了解两种描述流场的方法的区别与特点,重点掌握欧拉法下加速度的表达和意义掌握流体微团的几种变形和运动及其数学表达,掌握流体微团的运动分解与刚体运动的异同;了解系统分析方法与控制体分析方法的区别与联系;基本方程是本章重点,积分形式方程要掌握质量方程和动量方程的表达和意义,并会用它们解决实际工程问题;微分形式方程要重点掌握连续方程和欧拉方程的表达和意义;了解微元控制体分析方法;掌握伯努利方程的表达、意义、条件和应用;重点需要掌握的概念:流线、流量、散度、旋度、位函数、流函数、环量与涡的表达、意义及其相互之间的关系;2022-10-2623沈阳航空工业学院飞行器设计教研室