涡轮桨搅拌槽内混合过程的数值模拟-北京化工大学学报.pdf

第 31 卷 第 6 期 北 京 化 工 大 学 学 报 Vol. 31 , No. 6 2004 年 JOURNAL OF BEIJ IN G UN IV ERSIT Y OF CHEMICAL TECHNOLO GY 2004 涡轮桨搅拌槽内混合过程的数值模拟 张国娟 　闵　健 　高正明 　牛国瑞 　施力田 (北京化工大学化学工程学院 ,北京 　100029) 摘 　要 : 文中采用 FLU EN T 软件对六直叶涡轮桨搅拌槽 内的混合过程进行了数值模拟 , 选用多重参考系法 ( ) ε MRF 及标准 k 模型 ,将速度场与浓度场方程分开进行求解 ,所得的混合时间的模拟结果与实验值相吻合 。同 ( ) 时用计算机流体力学 CFD 方法研究了不同的加料点、监测点位置及操作条件对混合时间的影响规律 ,模拟结果 表明:混合过程主要由搅拌槽内的流体流动所控制 ,混合时间与加料点及监测点位置密切相关 。研究结果对于工 业搅拌反应器的优化具有一定的参考意义 。 ( ) 关键词 : 混合时间; 计算流体力学 CFD ; 数值模拟 ; 搅拌槽 ; 涡轮桨 中图分类号 : TQ052 ( ) ( 　　随着计算流体力学 CFD 技术及软件的迅速发 流场均处于近似稳态 的完全湍流状态 Re 1 × 展 ,国外在搅拌混合过程的数值模拟方面做了大量 104 ) ;所有方程采用隐式格式进行求解 。 研究工作 ,尤其对单层标准六直叶涡轮桨搅拌槽内 2 　模拟策略 混合过程的数值模拟 ,从文献[ 13 ] 报道来看 ,模拟 用的网格较稀 ,所设浓度收敛残差较大 ,所得模拟结 2 1 　搅拌槽结构及网格划分 果与实验值或文献值有较大误差 ,一般约 20 % 。国 以圆柱形平底有机玻璃槽为模拟对象 ,搅拌槽 内在这方面刚起步 ,周国忠等[4 ] 对单层和双层六直 直径为 0476 m ,液位高度与槽径相同 ,槽内有四块 叶涡轮桨的混合过程进行了数值模拟 ,其结果没有 挡板 ,宽度为槽径的 1/ 10 ,离槽壁 0008 m 。六直叶 实验数据验证 ,研究主要侧重模拟方法的比较 。与 涡轮搅拌桨 ,桨叶直径与槽径比为 04 ,桨的离底距 已有的实验研究[5 ] 相对应 ,本文对单层六直叶涡轮 离为槽径的 1/ 3 ,以水为工作介质 ,实验中在槽一侧 桨搅拌槽内的混合过程进行了 CFD 数值模拟 ,研究 的液面处加示踪剂饱和 KCl 溶液 ,在槽另一侧的底 了搅拌转速 、不同的示踪剂加料点及监测点位置对 部用电导率仪监测示踪剂浓度随时间的变化 ,从而 混合时间的影响规律 ,并与实验结果进行了比较 。 确定混合时间。 因示踪剂加入后搅拌槽内浓度场不对称 ,所以 1 　模拟方法 选取整个槽体作为计算域 , 网格创建采用结构化与 本工作采用 CFD 商用软件 FL U EN T 6 0 进行 非结构化相结合技术 ,桨叶区采用四面体非结构化 模拟计算 。从数值模拟的角度来看 ,模拟搅拌槽的 网格 ,且对该处