六直叶Rushton涡轮搅拌桨三维流场研究的论文.pdfVIP

中国工程热物理学会 流体机械 学术会议论文 编号：097007 六直叶Rushton 涡轮搅拌桨三维流场研究 1,2 1 1 栾德玉 周慎杰 陈颂英 （1 山东大学机械工程学院，济南 250062） （2 青岛科技大学机电工程学院，青岛 266061） Tel E-mail: qddy05@163.com 摘 要：利用计算流体力学(CFD)的滑移网格法计算了六直叶 Rushton 涡轮搅拌设备内的流场，模拟 了搅拌桨叶端及附近区域流场的速度分布及湍动行为，并与实验数据作对比。结果表明：CFD 计算的 时均径向和切向速度与已有的实验结果吻合很好，搅拌桨近桨区的速度分布并非对叶片高度的中心位 置严格对称，搅拌设备的流场结构并非完全由搅拌桨的行为决定；采用标准k −ε模型对湍流动能的 预测仍然偏低，要进一步改进对桨叶附近流动场湍动能的预报，应采用能够反映各向异性的模型算法。 关键词：滑移网格；Rushton 涡轮搅拌桨；数值模拟 0 前言 搅拌槽广泛应用于各种工业操作中,如液体混合、气体分散、固体悬浮及化学反应 等，如何准确地描述和模拟搅拌槽内的混合过程及流动状况,为混合设备的设计优化提 供理论指导，是混合技术的一个重要发展方向。近年来,利用计算流体力学(CFD)方法研 究搅拌槽内的流动场已经成为搅拌反应器研究的新热点，人们采用实验研究和数值模拟 相结合的手段，在三叶桨搅拌设备内部流场的研究方面作了大量工作，并取得了许多有 益的成果[1,2]。本文采用滑移网格法对六直叶 Rushton 涡轮搅拌设备内的全流场进行模 拟,同时考虑由于搅拌桨与挡板的相对运动所引起的流场周期性变化，并输出瞬时流场 及时均流场。 1 计算流体力学模型 非稳态状态下，场变量ϕ的输运过程应满足守恒定律，控制方程用(1)式表示： ∂ ρϕ +div ρϕu div =Γgradϕ +S ⑴ ( ) ( ) ( ) ϕ ∂t 方程左端第一项为场变量ϕ随时间的变化率，第二项为对流项，方程右端第一项为 扩散项,第二项为源项。上式对应于质量、动量、湍流动能和耗散率的守恒方程时，分 别为1、u 、k 和ε 。 j 考虑到搅拌槽为圆柱形平底结构，采用柱坐标形式： ∂ ∂ 1 ∂ 1 ∂ ρϕ + ρu ϕ + rρu ϕ + ρu ϕ ( ) ( Z ) ( r ) ( θ ) ∂t ∂z r ∂r r ∂θ ⑵ ∂ ∂ ϕ  1 ∂  ∂ϕ  1 ∂ Γ∂ ϕ  =Γ + Γr +  +S ∂z ∂z r ∂r ∂r r ∂θ r ∂θ       其展开形式如表1 所示。 1 表1 控制方程中的通用变量ϕ ϕ Γ S 1 0