李明辉;李海风
【摘 要】通过对硫铵结晶底入式搅拌器流场的CFD数值模拟,分析了搅拌转速、搅拌桨半径对高效轴流式搅拌器流场的影响.研究表明,搅拌转速对搅拌釜内流场的分布形式影响较小,随着转速的增加,流量和功率随之增加,近似呈线性规律分布,搅拌器内流体具有较好的循环性能;随着搅拌半径的增大,搅拌器内速度分布逐渐均匀,搅拌桨对流体的剪切性能较好. 【期刊名称】《化工机械》
半球型摄像机
无功功率计算【年(卷),期】烧烤箱2014(041)006
【总页数】5页(P794-798)
【关键词】底入式搅拌器;流场;CFD数值模拟
【作 者】李明辉;李海风
【作者单位】中国石油天然气股份有限公司抚顺石化分公司;中国石油天然气股份有限公司抚顺石化分公司
【正文语种】中 文
【中图分类】TQ051.7
底入式搅拌器主要用于化学工业、制药工业、化肥及材料工业等相关领域中，尤其在生物大分子药物高效化、己内酰胺中和与结晶系统和洁净大颗粒硫铵连续结晶系统中应用广泛。传统的径向流搅拌器普遍存在全罐混合差，桨区局部剪切作用强烈及功率因数偏高等缺点，而新型高效轴流式搅拌器与之相比则传质效率高、搅拌功率低、混合效果良好。笔者主要通过对高效轴流式搅拌器流场影响因素的研究，为底入式搅拌设备的优化设计提供理论基础。
1 设计参数及结构设计
根据某石化厂硫铵结晶装置使用工况，制定硫胺结晶搅拌器的设计技术参数如下[1]：
介质 硫铵母液(硫铵、有机物、水)
使用压力 常压
釜内温度 常温
安装方式 底入式
桨叶形式 高效轴向流翼型
搅拌转速 100～150r/min
桨叶回转直径 2 084mm
搅拌轴直径 120～125mm
导流筒内径 2 200 mm
密封形式 双端面机械密封+ 填料密封(API682 Plan 54+32)
硫铵结晶底入式搅拌器总体结构设计如图1所示。
美光隐形眼镜
图1 搅拌器总体结构图1——吊装平台；2——电动机；3——原动机组；4——减速机；5——膜片联轴器； 6——传动装置；7——集装密封；8——安装法兰；9——导流筒； 10——搅拌装置
2 搅拌器流场的数值模拟
以高效轴向流翼型搅拌桨为模拟对象，采用CFD方法来模拟和预测不同搅拌转速及搅拌桨半径对搅拌流场的影响，分析在高效轴向流翼型搅拌桨作用下的流场分布。
搅拌器的主要流动特征是湍流,采用多重参考系法(MRF)，流体为不可压缩的牛顿流体，研究体系为定常体系，其中描述湍流基本方程表达式为[2，3]：
式中 φ——通用变量，可代表u、v、w等求解变量；
Γ——广义扩散系数；
S——广义源项，其具体代表式和湍流模型见文献[2]。
采用FLUENT前处理软件GAMBIT2.2对搅拌釜和搅拌桨分别划分网格，由于搅拌桨区域是
计算中的核心，对其进行加密，网格数约为299 500个；搅拌釜约832 300个网格。两者均采用四面体非结构网格，图2、3显示了轴向截面上的网格。数值方法为SIMPLE算法[4]。计算域为1/3的完整模型，选取区域包括一个完整的搅拌桨叶。
图2 搅拌釜的网格划分图3 桨叶区的网格划分
3 模拟结果与分析
3.1 搅拌器的剪切-循环特性
氧化挂具
搅拌功率是混合设备设计的重要参数,常用功率准数NP和排出流量准数Nqd来表征搅拌器的剪切-循环特性。由于搅拌器对流体搅拌所作的功都消耗在对流体进行剪切和促使流体循环两方面,故用NP/Nqd表征消耗于剪切和循环的能耗比值参数, 其值越大表示叶轮搅拌功率中消耗于剪切的比例越大[5，6]。功率准数NP和排出流量准数Nqd为：
式中 d——叶轮直径；
n——转速；
NP——功率准数；
Nqd——排出流量准数；
P——功率；
Q——桨叶排出流量；
ρ——密度。
3.2 搅拌转速对流场的影响
转速是搅拌桨设计时的一个重要参数。在其他条件不变的情况下，改变桨叶转速进行CFD数值模拟，结果见表1。
表1 不同转速下的模拟数据转速r·min-1流量m3·s-1力矩N·m功率kWNPNqdNP/Nqd1003．781170．4947．800．270．251．051104．371408．7257．530．240．270．901204．761676．7268．480．220．260．831305．091981．6680．930．200．260．781405．562282．0393．200．190．260．711505．952619．47106．980．180
．260．66
图4为搅拌流场的速度分布云图，其流场的分布大致相同，搅拌转速的变化对搅拌釜内流场的分布形式影响较小，导流筒内均匀的分布5个速度区，表明流场中的速度梯度相同，但反应釜内流场的绝对速度随着转速的增加而增加(正相关)。
图4 不同转速下搅拌流场的速度分布云图
由转速与流量、搅拌功率的关系图(图5)可知，随着转速的增加，流量和功率也随之增加，近似呈线性规律分布。图6为转速与NP/Nqd比值的关系曲线，可以看出NP/Nqd值随转速增加而下降，说明当转速增加时，搅拌器内循环性能更好，搅拌桨的轴流性能得到更大的发挥。但是，为控制搅拌器内结晶的二次成核速度，需要在保证结晶器内悬浮液混合良好的基础上减小搅拌转速。
塑料表面电晕处理机图5 不同转速下流量和功率曲线
图6 不同转速下NP/Nqd值曲线
3.3 搅拌桨半径对流场的影响
分别对不同搅拌桨半径条件下的搅拌流场进行CFD数值模拟，结果见表2。
表2 不同搅拌半径下的模拟数据桨半径mm流量m3·s-1力矩N·m功率kWNPNqdNP/Nqd9804．501501．8861．340．210．280．7610004．761721．6670．310．220．270．7910204．961935．2279．040．220．270．8210405．091981．6680．930．200．260．7810605．932535．63103．560．240．290．8310806．002729．32111．470．230．270．85
从图7可以看出，随着搅拌半径的增大，导流筒内的速度梯度从4个逐渐增加到6个，说明搅拌器内速度分布逐渐均匀，导流筒出口和进口处的速度逐渐增强，增加了出口的扩散性，提高了混合效果，此外由于搅拌器底部的速度增加，利于消除底部的搅拌死区。
图7 不同搅拌半径下搅拌流场的速度分布云图
由不同半径下流量和功率的曲线(图8)可以看出，随着搅拌半径的增加，流量和功率曲线大体上是呈增加趋势，说明了釜内循环流速有增加趋势。但是从图9可以看出，NP/Nqd曲线也在缓慢增加，说明随着搅拌半径的增加，剪切性能也在增加。
