《两相流与热物理-刘婕》第三章第三节.pptVIP

3.3.4. 能量方程 每个颗粒有自己的内能 ，自己的动能 。单位质量固体微粒流能量，是大量颗粒个体能量统计平均。 内能： 随机运动动能（不考虑颗粒的内部自由度：振动，转动等）： 有些学者例似于气体分子运动论那样，将 称为颗粒流的动力学温度[4] [4].L. Huilin, D. Gidaspow, E. Manger, “Kinetic theory of fluidized binary granular mixtures”, Physical Review E, 64: 61301-61309, 2001. 对于  因此我们得到 的方程: (3.3.20) 对于 于是我们得到关于 的方程： (3.3.21) 将(3.3.21)与(3.3.20)相加，我们可以得到固体微粒流的总能量方程： (3.3.22) 将同类量合并得到固体微粒流用运动论推导得到的能量守恒方程 (3.3.22) 第三节、流体——固体两相流基本方程 §3. 1固体微粒拟流体模型及其基本方程 我们假设固体微粒为一种拟流体的连续介质，和流体一起充满整个流场。考察一个由固定不动的封闭光滑曲面A所围的微元体 内，其中包含了质量为 的体积为 的固体颗粒，质量为 体积为 的流体。这个微元体内混合物的总质量 ，总体积 。 A 固体微粒的相密度定义为： (3.1.1) 固体微粒拟流体的分密度定义为： (3.1.2) 固体微粒拟流体的体积分数（容积率）定义为： (3.1.3) 流体相的体积分数为 (3.1.4) 显然有 流体的相密度定义为： (3.1.5) 流体的分密度定义为： (3.1.6) 显然，有关系 3．1．1 状态方程 若流体是理想气体，则根据完全气体状态方程，气相分压为： 混合物压力p近似定义为： (3.1.7) 气相分压可改写为： (3.1.8) 注意到它又该近似地为两相分压之和： 故可得固相分压是： (3.1.9) 固体微粒拟流体的状态方程可简单表示为： (3.1.10) 3．1．2 连续性方程 根据Leibniz公式： 由于封闭曲面是固定的，即 ，我们于是得到： 因此，在这个微元体中相k质量的变化为： 单位时间相k由控制体净流出的质量为： 如果在控制体内有相变或化学反应发生，单位时间相k由此产生的质量减小为： 根据控制体内质量守恒有： 由于控制体是任意取的，因此下述微分型成立： (3.1.11) 在笛卡尔坐标系下，方程可以进一步分为： 固体微粒流： (3.1.12) 流体： (3.1.13) 用同样的方法，我们可以推导得到两相混合物的连续性方程： (3.1.14) 当有下关系时，则有(3.1.14)＝(3.1.12)＋(3.1.13)。 (3.1.15) (3.1.16) 3．1．3 动量方程 相k在控制体内单位时间动量的变化为： 单位时间由控制体净流出的动量为： 单位时间控制体内由于相变或化学反应产生质量变化带来的动量变化为： K相在控制体内所受到的体积力（如重力等）为： K相在控制面上所受到的面积力为： 两相之间相互作用力，既出现在控制体内，也存在于控制面上，它包括阻力、压差力、附加质量力、Besset力、Saffman力、Magnus力等等，在两相系单位体积内，这个力大小表示为 。 于是，根据牛顿第二定律我们得到： 由于控制体是任意取的，故在笛卡尔坐标系下我们有：  (3.1.17)‘ 应用连续性方程（3.1.11），我们得到：