液-固脉冲流化床中浓度波传播与衰减.pdfVIP

第 2 卷 第 5 期 过 程 工 程 学 报 Vol.2 No.5 2 0 0 2 年 10 月 The Chinese Journal of Process Engineering Oct. 20 02 液– 固脉冲流化床中浓度波传播与衰减 1 1 2 晋国栋 , 刘大有 , 舒其望 (1. 中国科学院力学研究所非线性力学国家重点实验室，北京 100080; 2. Brown University, Providence, RI 02912, USA) 摘 要：采用数值计算和实验验证相结合的方法研究液– 固脉冲流化床中浓度波的传播和衰减. 当 脉冲开半周期T 和闭半周期T 都远大于颗粒弛豫时间τ 时，两相的惯性力之差在一个周期的绝大 2 1 p 部分时间中相对于重力很小，可忽略，这时由双流体模型的动量方程可推导出推广的 Richardson–Zaki 公式，双流体模型简化为局部平衡模型. 采用五阶精度WENO 格式求解浓度波传 播方程，得到了脉冲流化过程中浓度波传播与衰减的规律，与实验结果符合良好. 关键词：液-固脉冲流化床；浓度波；数值模拟；WENO 格式 中图分类号：TQ051.1 文献标识码：A 文章编号：1009–606X(2002)05–0392–08 1 前 言 [1,2] 脉冲流化作为改善流态化质量的一种方法，已有许多实验研究和实际应用 ，但对其物理机 制和动态过程还不太清楚. 为弄清脉冲流化过程的物理机制，近年来开展了一系列的理论与实验 研究[3–5] ，并采用双流体模型对气– 固、液– 固脉冲流化过程进行了模拟，计算结果与实验结果符合 良好. 但由于基本形式的双流体模型在数学上对初值问题是不适定的[6,7] ，流化速度U(t) 的脉冲幅度 受到一定限制. 以方波形式变化的U(t)为例，通过对双流体模型的分析表明，除了在U(t)突变后十 分短暂的时间(与单颗粒弛豫时间τ 同量级，约0.05 s) 内，一般说来，两相惯性力之差远小于重力 p 和相间阻力，可以忽略. 这时，由动量方程可推导出推广的Richardson–Zaki关系式，即局部平衡模 型( LEM , Local equilibrium model ). 该模型形式简单，反映了物理过程的最主要因素，避免了因基 本形式的双流体模型对初值问题的不适定性所带来的模拟困难, 同时可提高计算效率，节约计算 时间. 由Jiang 等[8]发展的WENO(Weighted essentially non-oscillatory scheme)格式是用来求解双曲型 微分方程的一种高精度有限差分格式. 由于该方法具有精度高、无振荡特性，在计算流体力学和 其他领域得到了广泛的应用. 本文采用五阶精度WENO 格式，求解 U(t) 以方波和正弦波形式变化 的脉冲流化过程，并对方波形式变化的脉冲流化过程进行了实验验证. 2 数学模型及数值方法 本文基于如下假设：一维等截面流动、两相不可压缩、忽略壁面阻力和由于颗粒碰撞和脉动 引起的固相压强，则双流体模型中相k(k=p, f 分别表示颗粒相和流体相) 的连续性方程和动量方程 分别为[9] ∂α ∂(α u ) k + k k 0, ( 1)