对旋轴流风机内部流场的数值模拟.docVIP

对旋轴流通风机气动性能的数值预估 庞 磊 李孝宽 / 北京市劳动保护科学研究所 李 嵩 王巍雄 / 清华大学航天航空学院 摘要：对某商用对旋轴流通风机设计工况附近多个工况，分别采用单通道定常流和整机非定常流模型进行了气动性能的数值预估，并和实测结果对比，探讨了这两种数值模型用于对旋轴流通风机气动性能预估的可行性。 关键词：对旋轴流通风机；气动性能；数值预估 中图分类号：TH432.1 文献标识码：B 文章编号：1006-8155（2008）04-0020-04 Numerical Prediction of Aerodynamic Performance for a Contra-rotating Axial-flow Fan Abstract: The numerical prediction of aerodynamic performance for a commercial contra-rotating axial-flow fan has been carried out near the design operation using the models of single channel steady flow and whole fan unsteady flow respectively. The comparison between the results from these two models and the results from real test shows that the feasibility of these two models used in prediction of aerodynamic performance for contra-rotating axial-flow fans. Key words: contra-rotating axial-flow fan; aerodynamic performance; numerical prediction 0 引言 对旋轴流通风机具有结构紧凑、流量大、压头大和效率高等优点，广泛用于矿井、隧道等使用空间小又要求强烈通风的场合。由于对旋轴流通风机将两个相对旋转的叶轮、电机和各种支撑都装在一个叶顶间隙很小的筒体中，不仅流场是复杂的非定常流，性能测试也很难准确，所以数值预估和性能实测数据不易吻合。2002年李秋实等用自编程序进行了定常和非定常的数值模拟[1-2]，计算只对设计工况进行，且忽略叶顶间隙的影响。定常计算预估总压比试验结果偏高了1％，而非定常计算却偏高11％，也没有预估效率比较。近年来有人用FLUENT程序对对旋风机又做了一些定常计算[3-6]，所得气动性能与试验结果吻合仍然不佳，有的根本没有试验对比，总的来看，都只能用于定性分析对旋风机的一些工作特点。 1 数值模拟方法 采用商用软件Fluent 6.3对型号为FBD №5/2×7.5kW的一个商用对旋轴流通风机内部三维流场，分别采用单通道定常流和整机非定常流模型进行数值模拟，由此得到风机 气动性能。表1所示为该风机的主要设计参数，叶片采用等厚度薄板，其叶型由叶片上7个截面的不同的双圆弧型线组合而成。网格由Gambit软件自动生成，采用结构和非结构混合型网格填充计算域，计算采用RNG k-ε湍流模式，它比常用的标准k-ε模式可更好地计算分离流动，并采用标准壁函数和SIMPLE算法，动量方程的对流项采用二阶迎风格式，迭代采用多重网格技术加速收敛。建模时除省略电机及轮毂支架外，完全模拟实际试验时的风机整体结构, 特别是两个轮毂间的距离和叶顶间隙，它们对性能影响很大，必须严格确保相同。 1.1 单通道定常数值模拟 由于流动的轴对称性，不同叶片通道内流动几乎相同，于是可以只计算一个叶片通道内的流动来模拟实际试验时的风机整机流场，其建模的计算域和网格如图1所示。为了网格划分的方便，将计算域的几何模型划分为8个子区域：前级动叶区域R4、后级动叶区域R5、进气通道区域R3、出气通道区域R6、进气导流罩区域R2、出气导流罩区域R7，此外还在进气导流罩前添加了长度为5倍叶轮直径的一段延伸区域R1，在出气导流罩外添加了一个半径为5倍叶轮直径，长度为10倍叶轮直径的柱形风室区域R8，上述最后的两个子区域的添加是为了模拟风机在进气法气动性能试验时所具有的真实边界条件。在叶顶间隙区域和叶片的前后缘附近需加密网格，网格总数约70万。 表1 风机设计工况主要设计参数 叶轮直径Dt/ mm 500 转速 N/(r/min) 2900 流量Q/( m3/h) 11000 总压/ Pa 3000 轮毂比ν 0.6 前级叶片数Zf 9 后级叶片数 Zr 7 图