化学工程课件-连续性方程与运动方程.pptVIP

剪应力 （2-34a) （2-34b) （2-34c) τ与速度关联起来 法向应力 （2-35a) （2-35b) （2-35c) τ与速度关联起来 剪应力和法线应力 （2-34a) （2-34b) （2-34c) （2-35a) （2-35b) （2-35c) 粘性流体的运动微分方程（Navier-Stokes方程） 将（2-35) 代入上式： 粘性流体的运动微分方程（Navier-Stokes方程） 5个未知数，ux，uy，uz，ρ，p加上连续性方程和状态方程f(ρ，p)=0，5个方程，原则上可解。但由于非线性偏微分方程，目前还无法求其通解。 为此，需根据实际加以简化，去掉一些项，使之可解 柱坐标 球坐标 球坐标 动量传递 第二章 连续性方程与运动方程 连续性方程（微分质量） 微分能量方程 运动方程（微分动量） 微分质量衡算方程 单组份系统： （输出的质量流率）—（输入的质量流率） ＋累积的质量速率＝0 在x左侧面： 输入微元体积的质量流率 输出微元体积的质量流率 z x y dz dx dy (x,y,z) dy dz ρux dy dz 微分质量衡算方程 于是得到x方向输出与输入微元体积的质量流率之差： 同理在y方向： Z方向： 微分质量衡算方程 （输出的质量流率）—（输入的质量流率）＝ 累积的质量流率＝ 质量衡算： 出—入＋累积＝0 微分质量衡算方程 写成向量形式： 展开： 连续方程式一般形式 几种算法符号及意义谢树艺，《工程数学—矢量分析与场论》，人民教育出版社，1978年，北京 哈米尔顿（Hamilton)算子： 梯度 散度： 微分质量衡算方程的进一步分析 由于密度ρ是空间（x，y，z）和时间的连续函数，及： ρ＝f(x,y,z,θ） 将密度ρ进行全微分： 写成全导形式 不同的导数 偏导数： 某固定点处流体密度ρ随时间的变化率。 全导数： 流体密度由于位置和时间变化而产生的变化率（观测者在流体中以任意速度运动）。 随体导数： 观测者随流体随波逐流运动，即观测者在流体中与流体流速完全相同的速度运动。此时： 随体导数 随体导数 一般情况，算符 可用下式表示： 算符 所表示的函数称为随体导数或实体导数、拉格朗日导数。 欧拉观点和Lagrange观点 欧拉观点 ： 流体运动的空间中固定某一位置和体积，分析这点所通过的流体的特性变化来研究整个流体的运动规律。位置和体积固定，质量随时间变化。如岸上观水，地面观测站。 Lagrange观点： 在流体运动的空间中选择某一固定质量的流体微元，观测者随此质点运动。观测其特征变化来研究整个流体运动规律。质量固定，位置和体积可不固定。如随船观水，气球探测。 微分质量衡算方程的进一步分析 与随体导数定义： 得： 随体导数的意义 局部导数： 在一个固定 点（x,y,z)该量ρ随时间的变化； 对流导数： 由于流体质点运动，从一个 点转移到另一个点时发生的变化； 所以上述方程式表明：流体微元体积上的一个点 在dθ时间内从进入微元体积的空间位置（x,y,z)移动到微元体积的空间位置（x+dx,y+dy,z+dz)时，流体密度ρ随间的变化率. z (x,y,z) x y dz dx dy 微分质量衡算方程的进一步分析 由 ∵ ρv=1，对该式求随体导数，得： 可得： ∴ （2-7b） （2-9） 比较（2-7b)与（2-9）： 体积变形率 速度向量的散度 体积变性率和线性变型率 x1 x2 体积变形率 速度向量的散度 几种特殊情况下连续方程简化 稳态流动，密度不随时间变化，即 简化为： 对于不可压缩流体，ρ于时间与空间无关： （2-12） （2-12）不可压缩流体的连续性方程。 柱坐标和球坐标连续性方程式 z x y (x,y,z)或 （r,Φ,θ） z x y (x,y,z)或 （r,θ,z） θ Φ θ 柱坐标和球坐标连续性方程式 柱坐标： 球坐标： 微分动量衡算方程 用应力表示的运动方程 （2-16） F—诸外力的向量合；M—流体的质量 U—流体的速度向量；θ—时间。 对于固定质量 且运动的流体而言 速度u对时间θ的随体导数，表示流体的加速度 惯性力＝外力＝（质量）＊加速度 微分动量衡算方程 对于微元流体 在x、y、z三个方向： 力：质