将9-16风机模型导入ICEM 进行网格划分，网格划分过程中对离心风机关键部位要进行加密处理，如叶轮、集流器、蜗舌、进气箱的转角处等。对风机的进口与出口适当延长，以保证计算的稳定性。设计原理分析原风机蜗壳内壁型线采用的是传统蜗壳型线设计方法，即不考虑壁面粘性摩擦的影响，气流动量矩保持不变，运用不等边基圆法绘制的近似阿基米德螺旋线。考虑到离心风机结构的复杂且不规
9-16风机
将9-16风机模型导入ICEM 进行网格划分，网格划分过程中对离心风机关键部位要进行加密处理，如叶轮、集流器、蜗舌、进气箱的转角处等。对风机的进口与出口适当延长，以保证计算的稳定性。设计原理分析原风机蜗壳内壁型线采用的是传统蜗壳型线设计方法，即不考虑壁面粘性摩擦的影响，气流动量矩保持不变，运用不等边基圆法绘制的近似阿基米德螺旋线。考虑到离心风机结构的复杂且不规则性，本文采用非结构四面体网格进行划分，其中无进气箱的离心风机网格数量约370万，网格质量为0.3以上；带进气箱的离心风机网格数量为380万，网格质量为0.3以上。
9-16风机采用标准k-模型，壁面函数为Scalable，数值计算方法为高阶求解格式，求解格式为一阶格式。由于通风机转速低，马赫数小，可认为气流为不可压缩定常流动。进口给定质量流量，出口给定静压,壁面条件为无滑移边界，转速为1 480r/min，并将流动区域分为静止域与旋转域，两者通过Interface连接，连接模型为普通连接，坐标变换为转子算法，网格连接方式为GGI。测量时，除地面外无其他的反射条件，测点位置D距地面的高度与风机出口中心持平，水平方向上与出气口轴线成45°，距离出气口中心L=1m。本文所研究的某离心风机叶轮有均布的16 个前向的大小叶片，其内部流场较为复杂，为了揭示9-16风机内的流场特性，对风机进行全三维数值模拟。先单独分析了进气箱内部流场特性，然后对进气箱与风机进行一体化分析，研究进气箱对离心风机性能的影响。
1）9-16风机在进气箱出口与叶轮进口处有涡旋产生，其位置与流量大小相关，涡旋的存在导致叶轮流道发生了堵塞，是离心风机效率降低的原因之一。
2）加进气箱后，风机叶轮尾缘的“尾迹-射流”现象更加的严重，且在小流量区风机内部流场存在偏心现象。
3）加进气箱后9-16风机不仅效率有所降低，其全开流量与压力与无进气箱相比也有所下降，加进气箱后离心风机较优工况点向小流量区偏移，进气箱内部流场的复杂性以及出口速度的不均匀性对风机内部的流场分布产生了影响。
4）相比于无进气箱的情况下，加进气箱后，风机随流量的增加，噪声提升的更快，且在大流量区明显高于不带进气箱的噪声。
5）与实验测试结果对比分析，结果表明采用数值模拟研究风机性能是可行的。
为了提高掘进工作面离心风机导流效果， 提出对9-16风机圆弧形集流器加米字支撑架改造。从图中可以看出，每一种方式都有着不错的降噪效果，其中C型改进风机降噪效果好，在额定工况点附近总A声级能降低约7dB(A)。通过建立离心风机几何模型和数值模型,并施加边界条件，利用Fluent 软件对加米字圆弧集流器和普通圆弧集流器离心风机进行了整机内部流场数值模拟， 采用Tecplot 软件进行后处理，显示同流量下离心风机的压力云图。
将建立好的9-16风机三维模型导入ICEM 软件进行混合网格的划分。其中进出口和叶轮区域采用结构化网格，而蜗壳部分由于其内部结构复杂，尤其是电动机周围结构并非规则模型，故采用适应性较强的非结构化四面体网格，具体网格如图3 所示。综合考虑动静耦合区域对数值模拟预测结果的影响，在进行网格划分时，对边界层进行加密处理，其较低网格质量雅克比[14]在0.3 以上。消声蜗壳为A组合形式时与原风机的出口A声级随流量变化的对比图。为了保证数值计算结果的准确性，避免网格误差对其模拟结果造成影响，对9-16风机进行网格无关性验证，如表1 所示。综合考虑计算精度和计算效率可知，当网格数为25 万左右时预测结果较为合理，终确定整个计算域的网格数为2513558。k-ε 模型作为为普遍有效的湍流模型，能够计算大量的各种回流和薄剪切层流动，被广泛应用于各类风机的数值求解计算中。
由于有梯度扩散项，模型k-ε 方程为椭圆形方程，故其特性同其他椭圆形方程，需要边界条件：9-16风机出口或对称轴处k / n0和/ n0。但上述边界条件只针对高雷诺数而言，在固体壁面附近，流体粘性应力将取代湍流雷诺应力，并在临近固体壁面的粘性底层占主要作用。而多翼离心风机由于结构尺寸小、相对马赫数低，气体黏性力在流体流动过程中起重要作用，因此，在实际运用过程中，标准k-ε 模型由于未充分考虑粘性力的影响，导致计算模型出现偏差。B组合改进风机全压降低了约5．0Pa，9-16风机效率下降了约0．9%。运用Visual C++将上述修正函数编写为UDF代码，并导入Fluent 内置Calculation module。为符合实际运行状态，9-16风机进出口边界条件设置为压力入口和压力出口，出口压降与动能成正比，从而避免在进口和出口定义一致的速度分布[15]。后以CFD 计算的定常结果作为初始条件，进行非定常数值计算。
（作者： 来源：）
