0引言如何以最小的能耗获得所需的搅拌效果一直以来是固液悬浮技术研究的重点.因此，找出固液体系的临界离底转速显得尤为重要［1］.当颗粒全部处于运动时，且颗粒在槽底停留时间不超过1~2 s，即认为达到了完全离底悬浮，而达到此状态所需的最低搅拌转速称为临界转速NJS.在固液悬浮操作中，对完全离底悬浮的研究较多，也发表了不少有关临界转速的关联式.其中，Zwietering提出的完全离底悬浮临界转速关联式应用的最为广泛.目前得出的NJS值大多是在实验室尺度的槽中通过冷模实验估算获得的，再放大到工业反应器中.而对于工业上应用的高温、高压等复杂条件下的大尺度反应器，是无法通过实验获得所需的详尽信息的.因而需要寻找一种新的途径来对工业中应用的搅拌槽内的临界悬浮转速NJS进行预测.计算流体力学（CFD)技术的发展让这种预测成为可能.关于完全离底悬浮临界转速NJS方面的报道比较少，本文的目的是预测搅拌槽内固液悬浮的临界悬浮转速，将模拟结果和文献数据进行对比，验证Fluent软件对搅拌槽内临界离底悬浮转速进行预测的合理性.1计算体系与模拟方法1.1计算体系
1.1.1搅拌槽结构模拟采用的搅拌槽槽体为圆柱形，四块挡板均布，搅拌槽直径T=0.5 m，液面高H=0.46 m，挡板宽W=T20，挡板离底距离为0.01 m.浆型采用45度斜叶桨（搅拌桨示意图如图1所示,参数如表1所示），桨叶离底距离C=T3.图1PBT搅拌桨示意图
Fig.1Configuration of PBT airfoil impeller表1搅拌桨参数表
Table.1parameters of the PBT airfoil impeller
参数D/mmD·T-1b/mmθ/（°）σb搅拌桨2000.4404541.1.2计算物系介质为磷石膏和硫酸，固体相的体积分数为4.65%，颗粒的平均直接为72 μm，密度为2 320 kg/m3.流体操作密度为1 500kg/m3，操作粘度为5 MPa·s.
1.1.3网格划分根据流动的对称性，选取槽体的四分之一为计算域.由于模拟所采用的搅拌桨的结构不规则，网格划分采用非结构化的网格，共1 014 262个网格节点，如图2所示.图2网格划分示意图
Fig.2Meshing for the stirred tank1.2模拟方法计算使用的软件是Fluent 6.3，利用Euleian模型模拟该固—液两相流［23］.Euleian模型分别解每一相动量方程，并通过连续相对分散相的曳力与动量交换项来模拟两相间相互作用.并采用标准的壁面函数法来处理近壁区域流动的计算.而液相湍流模型采用标准kε方程模型.计算中假设固—液两相间无质量传递，由于固相浓度较低，所以相应的阻力系数采用wenyu模型［45］.流动场的计算使用多重参考系法，并设定流动性质为定常流动.相应的压力速度耦合选择Phase Coupled SIMPLE法，并采用一阶迎风的差分格式［67］.2结果与讨论均匀悬浮临界转速的确定最常用的方法是通过测釜内各点的固相浓度，根据釜内固相浓度分布（固体颗粒的体积分数分布情况反映了固体颗粒的的浓度分布状况）的均匀度来判断.衡量搅拌釜内固体颗粒浓度分布均匀性的判据很多.广泛采用的是浓度分布的标准偏差σ.σ=1n∑ni=1(cico-1)2
式中：n为测点数目.cico为测点固相浓度与全釜平均固相浓度之比.σ越小，固体颗粒在釜内分布的均匀程度就越高.随着功率（转速）增加，σ减小并趋于定值，即达到均匀悬浮.此时所需要的最小转速（功率）即为均匀悬浮临界转速.达到均匀悬浮时的均匀度（σ值）与叶轮型式和转速有关［8］.本文研究选取数值模拟的搅拌转速为20 r/min、30 r/min、40 r/min、50 r/min、60 r/min，70 r/min，其对应的雷诺数分别为0.83×105，1.25×105，1.67×105，2.08×105，2.5×105，2.92×105.2.1均匀悬浮状态均匀悬浮状态下的速度、浓度分布如图3所示.第10期董厚生，等：磷酸反应槽内临界离底悬浮转速的CFD模拟
武汉工程大学学报第33卷

图3固液悬浮流场分布
Fig.3Flow field for dilute suspension图3为60 r/min转速下的Fluent数值模拟后的混合流的液相速度矢量图与固相体积分数分布图.图3（a）为经过搅拌轴垂直面内的流场，其中矢量长度代表速度大小，搅拌桨区域的速度值较大，而自由液面区的速度值较小.图3（b）为固相体积分数分布图，表示了槽内的固体分布情况.2.2不同转速下的颗粒浓度分布图4列举了20~70 r/min转速下的浓度分布图.图4不同转速下垂直截面固相体积分数分布
Fig.4Vertical direction volume fraction contours of the
solids with different rotation speeds从图4中固体颗粒的体积分数分布可以看出，20 r/min转速时，液面处的浓度大小跟其他转速下的浓度大小不在一个数量级上.这是由于在液面处，速度过小时，重力起主导作用造成的.另外，随着转速的增加，槽内颗粒浓度更加趋于均匀状态，在搅拌槽的自由液面区颗粒浓度也随之增大.随着转速的不断增加，槽内颗粒浓度分布的标准偏差σ趋于定值时，即认为达到了均匀悬浮.此时所需的最小转速即为均匀悬浮临界转速［9］.2.3临界悬浮转速归纳以上几种转速下槽内的颗粒浓度分布状况，绘制沿图5中测试线L中的几个测试点随着转速的增加其颗粒浓度的标准偏差σ，结果如图6所示.图5测试点示意图
Fig.5The distribution of testpoints图6颗粒浓度偏差图
Fig.6the concentration variance of solids从图6中可以得出转速在20~70 r/min内，槽内上层颗粒浓度偏差均较大，往下则偏差越来越小，当转速达到40 r/min时，槽内的浓度分布标准偏差σ减小并趋于稳定，即认为达到了均匀悬浮状态，由此可推断得到本体系的的颗粒悬浮转速NJS=40 r/min，根据文献经验公式计算得出NJS为38~45 r/min，经计算两者误差为5%~125%［1012］.3结语以上是使用计算流体力学的方法对搅拌槽内的低浓度固液悬浮状况进行了模拟，通过模拟计算，得出不同搅拌转速下的固体颗粒浓度分布图.利用浓度判据的方法推断出临界离底悬浮转速NJS.模拟结果与文献结果比较吻合，不过，此结果仍需用实验来验证，初步验证了CFD技术对NJS预测的可能性，模拟方法与结果可对设备的设计和放大起到一定的指导作用.参考文献：