颗粒绕流问题中几何形体效应的数值模拟.pdfVIP

第 31 卷 第 1 期 北 京 化 工 大 学 学 报 Vol. 31 , No. 1 2004 年 JOURNAL OF BEIJ IN G UN IV ERSIT Y OF CHEMICAL TECHNOLO GY 2004 颗粒绕流问题中几何形体效应的数值模拟 高继贤 　刘 　辉 　李成岳 　杨立英 (北京化工大学可控化学反应科学与技术基础教育部重点实验室 ,北京 　100029) ( ) 摘　要 : 运用计算流体力学 CFD 模拟方法 ,针对球体、正六棱柱、正四面体、圆柱体和环轮体等不同几何形体颗 粒 ,在雷诺数为 10～400 条件下研究了颗粒绕流问题中的几何形体效应。基于模拟结果 ,考察了传统非球形颗粒 几何表征方法的不足 ;据此定义一个新的表征形体效应的参数 ,并在大量模拟计算结果的基础上 ,关联出以上四种 非球形颗粒曳力因数计算公式 ,改进后的曳力因数计算精度大大提高。 关键词 : CFD 模拟 ; 颗粒 ; 曳力因数 ; 两相流 中图分类号: TQ0211 　　颗粒绕流问题是流体力学的传统课题 。早期的 的 5 种不同形体的颗粒, 采用数值模拟方法, 考察单 研究囿于当时计算条件的限制, 往往只能获取近似 颗粒绕流的影响因素及传统等效方法的精度, 寻求 解析解, 如球体绕流的 Stokes 解及其 Ossen 修正 更恰当表征不同形体颗粒几何特性的参数, 提出异 等[ 1] 。随着计算流体力学工具的广泛应用, 使复杂 形颗粒曳力系数的计算方法。 几何形体绕流问题的数值求解成为可能, 这大大增 进了人们对于绕流问题的认识[ 2] 。另一方面, 颗粒 1 　绕流问题的界定和数学模型 绕流问题也是许多工程应用中的基础课题, 如在化 11 　模拟对象 学工程中的吸附床、催化反应器以及流化床中的流 参照相关文献[ 5 6] , 针对 5 种有代表性的催化 固相对运动等。这一问题的重要性在于, 对床层中 剂颗粒形状, 包括球形、圆柱形、正四面体形、正六棱 流动状况的解析, 往往是认识设备内传质、传热和反 柱和环轮体进行了模拟。 应性能的基础。 粒径是表征颗粒的重要几何参数。对非球形粒 近 10 年来, 在催化反应技术领域, 各种几何形 子, 比较常用的有等体积有效直径和等面积有效直 状的颗粒状新型催化剂被不断开发出来。在工程设 径, 统称为相当直径[ 3] 。前者是将非球形颗粒粒径 计中, 作为一种普遍采用的近似或等效方法, 大都采 定义为具有相同体积的球粒子的直径, 设粒子的体 用等体积相当直径和等面积相当直径来处理各种异 积为 Vp , 则体积相当直径 dv 为 ( ) [ 3 4] ( π) 1/ 3 ( ) 形 非球形 颗粒填充床内的流动问题 。对这类 d = 6 V / 1 v p 异形颗粒填充床流动以及压降的实验测定和计算表 若以外表面积 ap 相同的球形粒子的直径来作 明[ 5] , 用前述等体积或等面积相当直径近似处理在