






针对离心风机有无进气箱两种结构形式,建立了两种计算模型,利用CFX 软件对两种模型进行数值模拟,研究其内部三维流场特性,基于数值模拟结果分析了进气箱对离心风机的性能影响。数值模拟结果表明:加进气箱后,离心风机的全开流量与压力有所降低,缩短了有效工作区域;在离心风机内部叶轮进口处产生涡旋现象,堵塞了叶轮流道,使风机的效率和压力降低。数值模拟结果与实验测试值对比是比较吻合。进气箱是离心风机重要的组成部分,主要应用于大型离心风机与双吸离心风机。进气箱在其出口处气体发生近90°转弯,内部流场十分复杂,并造成很大的流动损失。其出口速度的不均匀性对离心风机性能影响明显,有必要对其特性进行研究。A.G.Sheard通过研究加进气箱的通风机,在离心风机叶轮进口加导流板控制叶轮进口的非均匀气流,结果表明在叶轮进口加导流板能够提高风机的全压,并得出了叶片根部断裂的原因。使用三维粒子动态分析仪(3D-PDA)对大型风机进气箱内部三维气体流场进行测量,揭示了其内部流动的基本特征,为了解进气箱流场结构和流动机理提供了依据。
离心风机叶片吸力侧形成的低能流积聚的“尾迹区”,形成“射流-尾流”结构。加进气箱后,风机叶轮尾缘处的“尾迹-射流”更加的严重,风机模型尾迹区占了比较大的空间,减少了风机流道有效面积。在小流量区,风机内部的流场分布发生偏心现象(C 处),叶轮流道E 侧,气体比较充实,叶轮流道F 侧气体分布较差,与原始风机内部流场分布相比,其离心风机叶轮流道的充盈性差。离心风机的效率曲线如图6,无进气箱情况下在流量为2.82kg/s,压力为3 106.23Pa 时,达到较率68.64%;加进气箱后在流量为1.68kg/s,压力为2 775.54Pa,达到较率59.45%,通过与原始风机对比可知,加进气箱后其较率降低8.19%。同样由图6 效率曲线对比图可知,加进气箱后风机整体效率降低,与原始离心风机相比其区域比较窄,缩短了工作区域,且加进气箱后较优工况点向小流量区偏移。加进气箱后,离心风机的全开流量降低,与无进气箱相比,流量降低了16.9%。由图7 可知,加进气箱不仅降低了风机的全开流量,其全压也有所减少。风机性能测试采用C 型试验装置对带进气箱的离心风机进行了性能测试,测试标准按GB/T 1236-2017《工业通风机用标准化风道进行性能实验》执行。
将建立好的离心风机三维模型导入ICEM 软件进行混合网格的划分。其中进出口和叶轮区域采用结构化网格,而蜗壳部分由于其内部结构复杂,尤其是电动机周围结构并非规则模型,故采用适应性较强的非结构化四面体网格,具体网格如图3 所示。综合考虑动静耦合区域对数值模拟预测结果的影响,在进行网格划分时,对边界层进行加密处理,其较低网格质量雅克比[14]在0.3 以上。为了保证数值计算结果的准确性,避免网格误差对其模拟结果造成影响,对离心风机进行网格无关性验证,如表1 所示。综合考虑计算精度和计算效率可知,当网格数为25 万左右时预测结果较为合理,终确定整个计算域的网格数为2513558。k-ε 模型作为为普遍有效的湍流模型,能够计算大量的各种回流和薄剪切层流动,被广泛应用于各类风机的数值求解计算中。
由于有梯度扩散项,模型k-ε 方程为椭圆形方程,故其特性同其他椭圆形方程,需要边界条件:离心风机出口或对称轴处k / n0和/ n0。但上述边界条件只针对高雷诺数而言,在固体壁面附近,流体粘性应力将取代湍流雷诺应力,并在临近固体壁面的粘性底层占主要作用。而多翼离心风机由于结构尺寸小、相对马赫数低,气体黏性力在流体流动过程中起重要作用,因此,在实际运用过程中,标准k-ε 模型由于未充分考虑粘性力的影响,导致计算模型出现偏差。运用Visual C 将上述修正函数编写为UDF代码,并导入Fluent 内置Calculation module。为符合实际运行状态,离心风机进出口边界条件设置为压力入口和压力出口,出口压降与动能成正比,从而避免在进口和出口定义一致的速度分布[15]。后以CFD 计算的定常结果作为初始条件,进行非定常数值计算。
