第二章 连续性方程与运动方程.pptVIP

. 牛顿型流体的本构方程 对于牛顿型流体的一维流动，当速度梯度与y轴方向相同时，剪应力与剪切速率(或形变速率)成正比，即 剪应力 τ与速度关联起来 对于三维流动．情况要复杂得多，每一剪应力与其相应两方向的形变速率有关。经分析推导，其关系为 法向应力 τ与速度关联起来 流体静止时，法向应力在数值上即为流体的静压力。当流体流动时，这一关系并不成立。它是由两部分组成的：其一是流体的压力，它使流体微元承受压缩，发生体积形变；其二由流体的粘性作用引起，它使流体微元在法线方向上承受拉伸或压缩发生线性形变。 下面六个剪应力和法线应力方程是直角坐标下牛顿型流体的本构方程 A 粘性流体的运动微分方程（Navier-Stokes方程） 将本构方程代入上式动量衡算方程，得： 整理，得 粘性流体的运动微分方程（Navier-Stokes方程） 5个未知数，ux，uy，uz，ρ，p加上连续性方程和状态方程f(ρ，p)=0，5个方程，原则上可解。但由于非线性偏微分方程，目前还无法求其通解。为此，需根据实际加以简化，去掉一些项，使之可解。 运动方程的最终形式为 柱坐标 第三章 连续性方程与运动方程 Euler观点和Lagrange观点 Euler观点 ： 流体运动的空间中固定某一位置和体积，分析这点所通过的流体的特性变化来研究整个流体的运动规律。位置和体积固定，质量随时间变化。如岸上观水，地面观测站。 Lagrange观点： 在流体运动的空间中选择某一固定质量的流体微元，观测者随此质点运动。观测其特征变化来研究整个流体运动规律。质量固定，位置和体积可不固定。如随船观水，气球探测。 物理量的时间导数 在动量、热量与质量传递过程中．众多物理量如密度、速度、温度等随时间的变化率，是传递过程速率大小的量度。物理量的时间导数有三种：偏导数、全导数和随体导数。下面以测量大气的温度t随时间θ的变化为例说明之。气温随空间位置和时间变化，可表为t＝t(x, y, z,θ)，t为空间对时间的连续函数。 偏导数：某固定点处物理参数随时间的变化率。为了测定大气的温度，可以将测温计装于观测站的某个空间位置，观测者记录下不同时刻的空气温度，此时得到的温度随时间的变化以 表示之，称为温度t的偏导数。 全导数：物理参数由于位置和时间变化而产生的变化率（观测者在流体中以任意速度运动）。测量大气温度也可采用下述方法：将测温计装在飞机上。飞机以一定的速度v在空间 飞行。观察者记录下不同时刻的空气温度*此时得到的温度随时间的变化以dt／dθ表示之，称为温度t的全导数。 全导数的表达式可由对t进行全微分得到 随体导数：观测者随流体随波逐流运动，即观测者在流体中与流体流速完全相同的速度运动。第三种测量大气温度的方法是将测温计装于探空气球上。此时探空气球随空气一起漂动，其速度与周围大气的速度相同。观察者记录下不同时刻的大气温度。如此获得的温度t随时间θ的变化称为随体导数(Substantial derivatives)亦称拉格朗日导数(Largrangian derivatives)，以Dt/Dθ表示。 连续性方程的推导 单组份系统： （输出的质量流率）—（输入的质量流率） ＋累积的质量速率＝0 在x左侧面： 输入微元体积的质量流率 输出微元体积的质量流率 z x y dz dx dy (x,y,z) dy dz ρux dy dz 连续性方程的推导 于是得到x方向输出与输入微元体积的质量流率之差： 同理在y方向： Z方向： 连续性方程的推导 （输出的质量流率）—（输入的质量流率）＝ 累积的质量流率＝ 质量衡算： 出—入＋累积＝0 写成向量形式 连续性方程的进一步分析 由于密度ρ是空间（x，y，z）和时间θ的连续函数，即ρ＝f(x,y,z,θ），那么ρ的随体导数 将连续性方程 展开，得到 于是 即 为连续性方程的另一表达形式。 随体导数 的意义 局部导数： 表示ρ在空间的一个固定点处随时间的变化； 对流导数： 表示密度由一点移动到另一点时所发生的变化。； 的物理意义为：当流体质点在dθ时间内由空间的一点（x, y, z)移动到另一点（x+dx, y+dy, z+dz)时，流体密度ρ随时间的变化率。 由于 将上式对时间求随体导数，亦即 上式的左侧 表示流体微元的体积膨胀速率或形变速率，右侧是速度向量的散度 几种特殊情况下连续方程简化 稳态流动，密度不随时间变化，即 上式可简化为： 对于不可压缩流体，ρ＝常数，则无论稳态还是非稳态：