A类颗粒流化床特性的DEM—CFD模拟.pdfVIP

第 26卷第 5期 化学反应工程与工艺 Vol26．NO5 2010年 1O月 ChemicalReactionEngineeringandTechnology Oct．2010 文章编号 ：1001--7631 (2010)05--0411--07 A类颗粒流化床特性的DEM-CFD模拟 殷 亮 蒋军成 (南京工业大学城市建筑和安全学 院，江苏 南京 210009) 摘要 ： 采用离散元方法一计算流体力学 (DEM—CFD)对耦合边界下 A类颗粒流化床特性进行模拟。讨论 了在床底的表观速度取决于流场情况下 ，不同进 口流速日寸该类颗粒 的流化状态 。分析 了流动与颗粒运动的 相互作用，及其对最小流化速度和最小气泡速度 的影响。发现与均匀床底的表观速度相 比，当床底 的流动 速度取决于流动与颗粒运动耦合时，床体 的稳定性变差，床底的速度分布呈现不均匀 ，局部更早 出现气 泡 ，并沿着该位置的轴向向床顶面发展 ，形成沟流 ，而此时其它部位会产生死 区。 关键词 ； 流化床 ；A类颗粒 ；耦合边界 ；离散元方法 ；计算流体 力学；数值模拟 中图分类号 ： TK229 文献标识码 ；A 近年来 ，在气固粒子流系统的模拟 中，基于离散元方法 (DEM)的模拟技术越来越 多地被研究 人员重视 。该技术以连续体处理流体相 ，通过求解牛顿运动方程追踪粒子轨迹 。DEM 的模型分为软 球模型和硬球模型。Cundall等[1]开发的软球模型在模拟准静态系统时能处理多颗粒间碰撞 ，其处理 的颗粒允许承受微小的变形 ，可采用这些变形来计算颗粒间的弹性力、塑性力及摩擦力 ，并采用牛顿 运动定律描述颗粒运动。而利用硬球模型模拟时，要处理一系列 的碰撞 ，且同时只能是一次碰撞 ，这 样颗粒间的作用力就不是显式计算了。已有不少学者使用 DEM进行流态化 的研究。Rhodes等 讨 论了颗粒间作用力对流态化特性的影响，Ye等L3探讨颗粒和气体属性对流化态的影响，着重比较了 DEM模型与经验公式 。Zhang等L4对上行和下行反应器中颗粒的团聚行为作 了DEM研究，讨论了颗粒 参数对团聚行为的影响。Renzo等 ]研究了液体和气体的流化特性，考察了系统对表观速度的响应。 将 DEM 与计算流体力学 (CFD)耦合来描述颗粒流体的流动，已经成功地模拟和预测颗粒工艺 的性能 ，如混合工艺[6]、颗粒 中密度 的旋流除尘器[7及气泡流化床反应器 中木材 的气化[8]。Zhong 等L9应用 DEM—CFD研究了喷动流化床颗粒运动特性，对颗粒在各 区的变化规律做 了总结。熊源泉 等L1对气固喷射器 内气 固两相流动进行 了三维数值模拟 。目前 国内学者多数 以 B类大颗粒作为 DEM—CFD的研究对象，国外文献主要报道了床底表观速度固定的情况 ，而实际过程中分布板离人 口 有一定高度 ，很难保证床底截面的表观速度是一致的，床底截面的表观速度取决于流动与颗粒运动耦 合的共同作用 。本工作基于三维 A类颗粒床 ，研究床底表观速度为耦合 的情况下，进 口流速不同时 A类颗粒的流化状态 ，并探索动力学机理，以期为DEM 应用于化学反应床提供参考。 1 数学模型和计算方法 1．1 颗粒运动模型 在流化床的DEM—CFD模型中，流体相的控制方程需要描述两相流中流体的流动，固体相要用 不同的方法和不同的尺度求解 。流化床内颗粒的运动包括平动和转动两种方式 ，其运动方程分别为： dvi 一 + F + F + Fg (1) 收稿 日期：2010—04—06；修订 日期 ：2010—07—27 作者简介：殷 亮 (1978一)，男，讲师。E-mail；umoni@sina．corn 412 化学反应工程与工艺 2010年 1O月 J