连续性方程和运动方程.PPT

1、1高等传递过程高等传递过程动量传递动量传递张浩勤教授张浩勤教授2主要内容主要内容 第二章第二章 动量传递的基本规律：动量传递的基本规律： 连续性方程连续性方程 运动方程运动方程 奈维奈维-斯托克斯方程（层流）斯托克斯方程（层流） 第三章第三章 奈维奈维-斯托克斯方程应用斯托克斯方程应用 平壁间流平壁间流动动 园管内流动园管内流动 爬流爬流 势流势流 平面流平面流 第四章第四章 边界层流动边界层流动 第五章第五章 湍流理论（雷诺方程（湍流）动湍流理论（雷诺方程（湍流）动量传递的概念量传递的概念3动量传递的应用 了解流场中速度、压力的分布规律； 解决与流体输送相关的单元操作过程中的问题； 为热质传

2、递研究奠定基础。4数学补充数学补充52场论基础知识场论基础知识6（2）场论）场论78本章重点：1 掌握研究问题的两种方法并能够应用； 2 偏导数、全导数、随体导数的差别； 3 连续性方程的推导和应用条件；4 运动方程推导思路和应用条件；5 熟悉场论中的数学表示方法。第二章 连续性方程与运动方程9研究对象 牛顿流体、单一组分、等温流动牛顿流体、单一组分、等温流动系统 10第一节 描述流动问题的两种观点 拉格朗日（Lagrange）观点： 在流体运动的空间中选择某一固定质量的流体微元，观测者随此质点运动（相对坐标系）。观测其特征变化来研究整个流体运动规律。 质量固定，位置和体积可变。 如随船观水，

3、气球探测。11 欧拉观点 ： 流体运动的空间中固定某一位置和体积，分析这点所通过的流体的特性变化来研究整个流体的运动规律 位置和体积固定，质量随时间变化。 如岸上观水，地面观测站。 强调:对同一流场，无论采用哪种观点，其结果都是相同的，只不过采用的观点得当，分析问题方便、简捷一些。12流线与轨线 轨线：流体质点在流场中的运动轨迹流体质点在流场中的运动轨迹，是拉格朗日法考察流体运动所得的结果。 轨线上某一点的切线代表质点的运动方向，轨迹给出了轨迹给出了同一质点在不同时刻的速度方向。（在黑板上画图）1314 流线：流线时这样的空间曲线，对于某一固定时刻，曲线上任一点处的速度方向和曲线在该点的切线方

4、向重合。流线是欧拉法考察的结果。 流线是同一时刻由不同流体质点所占据的空间曲线。它给出该时刻不同质点的运动方向。 （在黑板上画图）1516注意： 流线的性质：在任一时刻通过流场中任何一个空间点都有一条流线，所以流场中的流线是流线簇。在流线簇中，流线越密，代表速度越大。 流线是不能相交的。因为空间上每一点只能有一个速度方向，所以不能有两条流线同时通过一点。即流体不能穿越流线流动。 特例：在速度为零或无限大的空间点上例外，速度为零的点称为驻点，速度无限大的点称为奇点。 流线的形状和位置，在稳态流动中不随时间变化，在非稳态流动中，一般要随时间变化。 对于稳态流动，流场中任何参数均不随时间变化，故流线

5、方程与轨线方程重合。1718流线方程19202122系统与控制体 控制体：对于某个坐标系来说，固定不变的任何体积称之为控制体。控制体的边界称为控制面。控制面总是封闭的。应当指出，占据控制体的诸流体质点可以随时间而变化。 特点：控制面相对于坐标系是固定的。在控制面上可以有质量交换（即有流体进入或流出控制面）、热量交换（能量(热、功、内能等)输入或输出控制面）；控制面上也可以受到控制体以外施加在控制体物质上的力而引起动量交换。控制体内部质量、热量、动量的储存量可能改变。控制体可以是运动的，也可以是固定的。 控制体是欧拉法的结果。23 系统：所取控制体无质量穿越其表面，既没有流体进出，则此固定质量的

6、控制体称为系统。 特点：系统的边界面的形状、位置可以随时间而变化，系统的边界随着流体一起运动。 在系统的边界上，可以有热量交换、动量交换，但是没有质量交换。 系统是拉格朗日观点的结果。 小结：拉格朗日观点 系统 轨线 欧拉观点 控制体 流线 24不同的导数 偏导数： 某固定点处流体密度随时间的变化率。全导数： 流体密度由于位置和时间变化而产生的变化率（观测者在流体中以任意速度运动）。 随体导数： 观测者随流体随波逐流运动，即观测者在流体中与流体流速完全相同的速度运动。此时：ddzzddyyddxxddddzuddyuddxuzyx;DD25随体导数ddzzddyyddxxddzuyuxuddz

7、yxzuyuxuDDzyxztuytuxtutDDtctzyx等也有类似表达式、浓度对温度26随体导数 一般情况，算符 可用下式表示： 算符 所表示的函数称为随体导数或实体导数、拉格朗日导数。zuyuxuDDzyxDDDD),(zyxfddDDddzzddyyddxxdd27第二节 连续性（微分质量衡算）方程 前提条件：单组分等温流动系统前提条件：单组分等温流动系统 分析方法：欧拉法分析方法：欧拉法 控制体：控制体： 流体质点流体质点 28（输出的质量流率）（输入的质量流率）累积的质量速率029第二节 连续性（微分质量衡算）方程在x左侧面：输入微元体积的质量流率输出微元体积的质量流率zxydz

8、dxdy(x,y,z)dydzuxdydzdydzuxdxdxdfffdydzdxxuuxx0)(dxxuuxx)(30微分质量衡算方程 于是得到x方向输出与输入微元体积的质量流率之差： 同理在y方向： Z方向：dydzdxxuuxx)(dxdydzxudydzuxx)(dxdzdyyuuyy)(dxdydzyudxdzuyy)(dxdydzzuuzz)(dxdydzzudxdyuzz)(31微分质量衡算方程（输出的质量流率）（输入的质量流率）累积的质量流率质量衡算： 出入累积0dxdydzzuyuxuzyx)()()(dxdydz)(dxdydzzuyuxuzyx)()()(0)(dxdyd

9、zzuyuxuzyx)()()(0)(32微分质量衡算方程 写成向量形式： 展开：连续方程式一般形式zuyuxuzyx)()()(0)(0)()()()()()(zuyuxuzuyuxuzyxzyx0)()u33微分质量衡算方程的进一步分析zyxuzuyuxDD0)()()()()()(zuyuxuzuyuxuzyxzyx与随体导数定义：得：01)()()(DDzuyuxuzyx01)()()(DDzuyuxuzyx（2-7b）34随体导数的意义对流导数对流导数：由于流体质点运动，从一个点转移到另一个点时发生的变化；所以上述方程式表明：流体微元体积上的一个点在d时间内从进入微元体积的空间位置（

10、x,y,z)移动到微元体积的空间位置（x+dx,y+dy,z+dz)时，流体密度随间的变化率.z(x,y,z)xydzdxdyzyxuzuyuxDD35微分质量衡算方程的进一步分析 v=1，对该式求随体导数，得：（2-9）比较（2-7b)与（2-9）： 体积变形率 速度向量的散度011DDDDvvzuyuxuDDvvzyx136体积变性率和线性变型率x1 x21212xxuuxx12xxuu12xx 37几种特殊情况下连续方程简化稳态流动，密度不随时间变化，即简化为：对于不可压缩流体，于时间与空间无关： 不可压缩流体的连续性方程。 （2-12）0zuyuxuzyx)()()(0)(zuyuxu

11、zyx)()()(00zuyuxuzyx0u38柱坐标和球坐标连续性方程式zxy(x,y,z)或（r,）zxy(x,y,z)或（r,z）39柱坐标和球坐标连续性方程式 柱坐标： 球坐标：0)()(1)(1zruzurrrur0)(sin1)sin(sin1)(12ururrrurr40 思考题：推导球坐标系的连续性方程。 参见浙大教材41第三节 运动方程 衡算基础：动量守恒方程 研究方法：拉格朗日法42用应力表示的运动方程 （一）动量守恒定律在流体微元上的表达式（拉格朗日法） 牛顿第二定律：物体动量随时间的变化率等与该牛顿第二定律：物体动量随时间的变化率等与该物体所受外力的矢量和。物体所受外力

12、的矢量和。 （2-16） F诸外力的向量合；M流体的质量 U流体的速度向量；时间。 DuDMF)(惯性力外力（质量）加速度43 拉格朗日法：拉格朗日法：在流体运动的空间内，选择某一固定质量的流体微元（M为常数为常数），观察者追随此流体微元且一起运动（在相对在相对坐标系下，可以用随体导数的概念来描坐标系下，可以用随体导数的概念来描述述），并根据此运动流体微元的变化状况来研究整个流场中流体运动规律。 固定质量的流体微元： 体积 质量 =常数 所以所以 dxdydzdv dxdydzdvDuDdxdydzFd44微分动量衡算方程 对于微元流体 在x、y、z三个方向： 力：质量力或体积力FB，作用在整

13、个微元流体上； 表面力或机械力，作用在微元流体诸表面上的外力，计为FS.它又可分为法向力和剪应力。DuDdxdydzdFFdiDuDdxdydzdFFdzzizDuDdxdydzdFFdyyiyDuDdxdydzdFFdxxix45(二)、作用于微元体上的外力分析作用于微元体上的外力分析 合外力（惯性力）=质量力+表面力 质量力 作用在所有流体质点上的力，重力 离心力 电场力等。 表面力：作用在控制面上的力,在此即作用在流动流动着的流体微元表面的力着的流体微元表面的力. 机械力(接触力) 这些力是由该控制体毗邻的流体所产生的,由静压力和粘性力提供.对于单位面积而言,称为表面应力.表面应力分为法

14、向应力和切向应力(剪应力). 表面力=法向应力+切向应力46质量力（重力） 在x方向上： dFxB=Xdxdydz X-单位质量流体的质量力在x方向上的分量。 重力Xgcon=Fxg 当X方向为水平方向时， X=Fxg0,90度 当X方向为垂直方向，Xg9.81m/s2 X与重力方向可以相同，也可以不同gx47表面力yzxxxxyxzxy 第一个下表表示应力分量的作用面与x轴垂直。第二个下标y表示应力方向为 y轴方向。 xx 表示法向 应力分量。拉伸方向（向外）为正，压缩方向（向内）为负。 微元流体在运动时，由于法向应力和剪应力的存在，使其发生形变。48表面力每一表面的机械应力均可分解成三个平

15、行于x、y、z三个坐标轴的应力分量3个六个表面，36=18个zxydzdxdy上）(dyyyxyx右）(dxxxxxx后）(dzzzxzx前）(zx左）(xx下）(yx49当小微元体体积缩小为一点时，相对表面上的法向应力与切线应力都是相应地大小相等、方向相反的。故只需采用9个机械应力就可以完全表达：3个法向分量，6个切线分量。 以后将证明该矩阵为对称矩阵50 6个切向应力分量之间的关系 上述6个剪应力可以使小微元旋转且彼此不独立。可以由此关联起来。 这四个剪应力对于旋转轴线产生力矩：力矩质量力矩质量旋转半径旋转半径的平方的平方角加速度角加速度2dyyyxyx2dxxxyxy2dxxxyxy2d

16、yyyxyx dy/2 dx/2 o dx/2 dy/2xy51表面力力矩质量力矩质量旋转半径旋转半径角加速度角加速度（角加速度）旋转半径）2)()2)(22)2)(22dxdydzdydzdxdyydyydxdydzdxxdxxyxyxyxyxxyxyxyxy当小微元体积趋近于0使旋转半径趋近于0 同理： （角加速度）旋转半径）2(yxxyyxxyxzzxzyyz52用应力表示的运动方程zxydzdxdy上）(dyyyxyx右）(dxxxxxx后）(dzzzxzx前）(zx左）(xx下）(yxdxdydxdydzzdxdzdxdzdyydydzdydzdxxdFzxzxzxyxyxyxxxx

17、xxxxS)dxdydzzyxdFzxyxxxxS)(简化后：X方向表面力53X方向总的外力分量dFx 外力分量=质量力分量表面力分量xSxBxdFdFdFzyxXDDuzxyxxxxzyxgconDDuzxyxxxx（2-27a）54以应力项表示的粘性流体运动微分方程zyxXDDuzxyxxxxzyxYDDuzyyyxyyzyxZDDuzzyzxzz的偏微分方程法线应力、剪应力有关单位质量流体质量力、与密度、速度、时间、0),(zyxZYXuuufzxyxxxzyx55问题与讨论系方程。另外还需补充若干个关变量个数，量之间关系，减少独立三个方程，必须找出变个：未知量体积力）。，个：已知量，个

18、：自变量)()(),(,10(3;,4xzzxzyyzyxxyzzyyxxzyxuuuZYXzyxzyxXDDuzxyxxxxzyxYDDuzyyyxyyzyxZDDuzzyzxzz56应力与应变速率的关系 寻找关系如何着手？粘性流体在运动时各层之间会发生相对运动，那么粘性与流体的形变之间必然有一定的联系。 在三维流动中，应力与应变速率之间的关系十分复杂，法向应力的作用难以判断。 思路：将刚体力学中应力与应变的关系，改进后用于流体力学。 参考书：参考书： 王绍亭王绍亭 ，陈涛，动量热量与质量传递，天津科学，陈涛，动量热量与质量传递，天津科学技术出版社，技术出版社，1986年。年。57一维流动

19、剪应力（u联系起来) dydux58596061三维流动, 显然,由于粘性力的作用,流体微元会发生变形.62 对三维流动，如图所示的流体微元，其体积为。由力的分析可知，它在流动过程中会发生体积形变，即由原来的长方形六面微元体变为一个菱形六面微元体。对于x-y平面而言，起作用的切向应力分量有4个，其中xy 和yx分别作用在与x-y平面相垂直的4个平面上。相对应边上的两个应力方向等值反向。经过微分时间d后，原来的长方形变为菱形(图中虚线所示)，相邻两条边线的夹角减小。 636465剪应力)(xuyuyxyxxy)(zuyuyzzyyz)(xuzuzxxzzx（2-34a)（2-34b)（2-34c

20、)与速度关联起来66法向应力 在三维流动中，判断法向应力的作用更为困难。其推导过程较长，我们不打算详细介绍。同学们可参见其他参考书中的推导。在此，我们主要介绍一些基本的概念。 流体静止时，法向应力在数值上等于流体的静压力，但方向相反。 xx =yy =zz =p 或 p=1/3(xx +yy +zz)67 在流动流体中，法向应力由下列两种类型的应力所提供，其一为静压力，它使流体微元承受压缩应力，发生体积形变；其二是由流体流动时的粘性应力的作用产生，其结果是使流体微元在法线方向上承受拉伸或压缩应力，发生线性形变。法向应力与压力和粘性应力的关系，可写为： 表示粘性应力引起的附加法向应力。 6869

21、法向应力)(32)(2zuyuxuxupzyxxxx)(32)(2zuyuxuyupzyxyyy)(32)(2zuyuxuzupzyxzzz（2-35a)（2-35b)（2-35c)与速度关联起来70剪应力和法线应力)(xuyuyxyxxy)(zuyuyzzyyz)(xuzuzxxzzx（2-34a)（2-34b)（2-34c)(32)(2zuyuxuxupzyxxxx)(32)(2zuyuxuyupzyxyyy)(32)(2zuyuxuzupzyxzzz（2-35a)（2-35b)（2-35c)71粘性流体的运动微分方程（Navier-Stokes方程）将（2-34 2-35) 代入上式：z

22、yxXDDuzxyxxxx)()()(32)(2222222222222zxuzuyxuyuzxuyxuxuxuxpXDDuzxyxzyxxx72粘性流体的运动微分方程（Navier-Stokes方程） 其它方向略。见2-36 5个未知数，ux，uy，uz，p加上连续性方程和状态方程f(，p)=0，5个方程，原则上可解。但由于非线性偏微分方程，目前还无法求其通解。 为此，需根据实际加以简化，去掉一些项，使之可解)(31)(222222zuyuxuxzuyuxuxpXDDuzyxxxxx73讨论惯性力 质量力 压力 粘性力)(312uuPFDuDg74方程简化 对于不可压缩流体 见2-37)(2

23、22222zuyuxuxpXDDuxxxx)(222222zuyuxuxpXuzuuyuuxuuDDuxxxxxzxyxxx75柱坐标222222221)(11)(zuururrurrrrpXxuuruururuuurrrrrrxrrrr分量：22222221)(111zuururrurrrprXzuuruuururuuurzrr分量：2222221)(11zururrurrrpXzuuururuuuzzzzzzzzzrz分量：76球坐标ururururururrurrrrpXruuuruururuuurrrrrrrrrrrsin2cot222sin1)(sinsin111sin)(22222

24、22222222分量：urruurururrurrrprXruruuuruururuuurrr2222222222222sincos2sin2sin2)(sinsin1)(111cotsin分量：77球坐标ururruururrurrrprXruuruuuruururuuurr222222222222sincos2sin2sinsin2)(sinsin1)(1sin11cotsin分量：78讨论奈维-斯托克斯方程是用于牛顿流体，层流流动。奈维-斯托克斯方程是单位体积微元力的平衡式 惯性力 重力 压力 粘性力在不同的流动系统中，四种力所起的作用不同，视具体的情况可以简化。例如，对于理想流体，粘度

25、等于0，粘性力项必然等于0。一般而言，对粘性流体管内流动，重力的作用较小。但对于瀑布类的流动，重力的作用是不可低估的。对于静止流体，速度等于0，可以简化为静力学方程式。79奈维-斯托克斯方程原则上是可解方程。 未知数 ux uy uz p 共5个 方程数 奈维-斯托克斯方程 3 连续性方程1 流体状态方程 1 共5个 例如，对于不可压缩均质流体 =常数实际上，非线性方程组的解析解求起来十分困难，只能对简单的定解条件的情况得到。 大量的复杂的情况可借助于数值求解。80定解条件奈维-斯托克斯方程和定解条件一起才构成完整的数学模型。初始条件：传递现象满足的初始状态条件。=0，ux uy uz p 的值 对于稳态流动，没有初始条件。边界条件：传递现象在边界上满足的条件。边界形式多种多样，具体问题要具体分析。81常见的有几种：固体壁面的粘附条件：粘性流体在静止壁面上速度为零。 在运动壁面上 流体速度与运动壁面的速度相等。自由表面 通常自由表面指一个流动的液 体暴露于气体中的那一部分界面。自由表面上压力相等。一般的情况，认为在气液界面上，切应力是连续的。如果， 气液界面的摩擦力可