《武汉工程大学学报》  2014年02期 56-62   出版日期:2014-02-28   ISSN:1674-2869   CN:42-1779/TQ
搅拌器内部二维流场数值模拟


0引言搅拌器是使液体、气体介质强迫对流并均匀混合的器件,是以桨叶旋转来混合流体的设备,用于混合不同的物质[14].搅拌器使用非常广泛,在各个领域中都有应用[5].搅拌作是工业过程的基础环节,提高搅拌效率意义重大,高效混合流体是搅拌器的最终目的[6].计算流体动力学(简称CFD)技术具有初步性能预测、内部流动预测、数值试验和流动诊断等作用.FLUENT软件是目前处于世界领先地位的商业CFD软件包之一,是模拟和分析复杂几何区域内流体流动的专用软件.数值模拟是一种寻求提高搅拌叶轮性能的设计思想和设计方案中最快、最经济的方法[78].本文采用CFD技术用FLUENT软件来模拟在不同桨叶直径、转速和桨叶数的情况下,得出搅拌器内流场分布图像,故而可以直接观看到它的内部流动情况,以期为搅拌器的优化设计工作提供参考.1控制方程在定常条件下,搅拌流场的不可压缩流动可用以下方程组描述(1)连续性方程: Pt+(ρu)x+(ρv)y+(ρw)z=0(2)动量方程:uit+xi(uiuj)=Fi-1ρPxi+υ2uixixj(3)湍动能方程:uikxi-xiυ+υtσkkxi=Pr-ε(4)湍动能耗散率方程: uiεxi-xiυ+υtσεεxi=Cε1·ε·Pr-Cε2·ε2k2数值模拟过程 本文研究的是桨叶直径、桨叶转速、桨叶数对搅拌器搅拌性能的影响.采用桨叶直径分别为400 mm(长H=190 mm,宽L=20 mm)、500 mm(长H=240 mm,宽L=20 mm)、600 mm(长H=290 mm,宽L=20 mm)以研究桨叶直径对搅拌器的影响.采用桨叶转速分别为2 rad/s、4 rad/s、6 rad/s以研究桨叶转速对搅拌器的影响.采用四叶桨式搅拌器、六叶桨式搅拌器、八叶桨式搅拌器,浆片都是直列式,以研究桨叶数对搅拌器的影响.搅拌桶的半径R=40 cm,坐标系原点位于搅拌叶轮中心,基于旋转的速度和工质水,可以确定搅拌器内部流动为湍流.以四叶桨式搅拌器为例,如图1所示.图1旋转碗的尺寸示意图Fig.1The size sketch of rotating bowl 2.1建立模型 搅拌器二维模型使用Gambit建立,外圆和内圆组成区域1,内圆和桨叶组成区域2.以Map划分方式划分区域1,以Pave平铺方式划分区域2.网格图如图2~4所示.图2四叶桨式搅拌器网格划分图Fig.2Four vane agitator mesh figure图3六叶桨式搅拌器网格划分图Fig.3Six vane agitator mesh figure图4八叶桨式搅拌器网格划分图Fig.4Eight vane agitator mesh figure2.2边界条件调节条件:将搅拌桶的外壁和搅拌器桨叶的边定义为WALL,环形面和内部的面交界处是两条重合的边,设置它们为INTERFACE边界条件,名称分别叫做“interface_1”和“interface_2”.由于使用的是MRF模型,还需要定义区域,即对不同的子域进行命名.指定外面环面的类型为FLUID,名称为“stationary”;指定内部的面类型为FLUID,名称为“swirl”.2.3求解计算FLUENT软件选择单精度二维求解器;设置计算区域尺寸和单位制尺寸中的网格尺寸按厘米进行生成;使用压力基求解器;湍流模型选择kepsilon双方程模型;定义流体的物理性质中,在FLUENT Database中调出水的物理参数;求解控制参数的设置中,差分格式采用一阶迎风方式,压力速度耦合求解方式中选择SIMPLE算法.第2期徐建民,等:搅拌器内部二维流场数值模拟武汉工程大学学报第36卷3结果及分析3.1桨叶直径的影响以四叶桨式搅拌器为例,2 rad/s的桨叶转速,分别模拟桨叶直径为400 mm、500 mm、600 mm的桨叶,速度矢量图如图5~7所示.结果分析:从图5~7的速度矢量图中可以看到,随着桨叶直径的增大,使整个搅拌设备内的流场分布更均匀,流体混合更充分,搅拌器内速度分布变好,速度在各个区域内都很同步,但桨叶直径的增大,使得材料用量增加,制造费用提高,安装麻烦,还会增加其本身的受力,使其使用寿命降低.图5400 mm 桨叶直径速度矢量图Fig.5Velocity vector diagram of blade diameter d=400 mm图6500 mm桨叶直径速度矢量图Fig.6Velocity vector diagram of blade diameter d=500 mm图7600 mm桨叶直径速度矢量图Fig.7Velocity vector diagram of blade diameter d=600 mm3.2桨叶转速的影响以四叶桨式搅拌器为例,桨叶直径为400 mm,分别对桨叶转速为2 rad/s、4 rad/s、6 rad/s进行数值模拟,速度矢量图如图8~10所示.图82 rad/s桨叶速度的速度矢量图Fig.8Velocity vector diagram of blade angular velocity ω=2 rad/s图94 rad/s桨叶速度的速度矢量图Fig.9Velocity vector diagram of blade angular velocity ω=4 rad/s图106 rad/s桨叶速度的速度矢量图Fig.10Velocity vector diagram of blade angular velocity ω=6 rad/s结果分析:从图8~10的速度矢量图中可以看出,随着转速的提高,搅拌器内流场分布更加均匀,搅拌效果明显增强,但是转速的提高会加大搅拌轴的受力,电机功率也会增大,能耗增加,从而可能导致搅拌器和搅拌轴的结构尺寸需要重新设计.3.3桨叶数的影响本文以桨叶尺寸直径为400 mm转速为2 rad/s为例,分别对四叶桨式搅拌器、六叶桨式搅拌器、八叶桨式搅拌器进行数值模拟,速度矢量图如图11~13所示.图11四叶桨式搅拌器的速度矢量图Fig.11Velocity vector diagram of four vane agitator图12六叶桨式搅拌器的速度矢量图Fig.12Velocity vector diagram of six vane agitator图13八叶桨式搅拌器的速度矢量图Fig.13Velocity vector diagram of eight vane agitator结果分析:从图11~13的速度矢量图可以看出,随着桨叶数的增加,搅拌器内的搅拌效果得到增强,流场分布比较均匀,但效果不是很明显.同时增加桨叶数会带来一些负面效果,它会增加材料用量从而使得成本增加;也会使得材料应力增大,减少使用寿命,因此对于增加桨叶数的方法来提高搅拌效果应该综合各方面考虑,选择合适的桨叶数.4结语搅拌器桨叶直径、桨叶转速和桨叶数的增加有利于搅拌的混合均匀,但是桨叶直径和桨叶数的增加不仅使得搅拌器桨叶加工复杂,技术要求高,而且使得材料用量增加,成本也会提高;桨叶转速也受制于搅拌器和搅拌轴的结构尺寸,不能一味地增大.只有考虑各种因素,才能确定出最经济最合适的桨叶直径、转速和桨叶数.致谢感谢武汉工程大学机电工程学院过程装备与腐蚀际护实验室的全体人员对本研究的支持与帮助.