山东冠熙环保设备有限公司主营：轴流风机,耐高温高湿风机,烘干设备用风机,离心风机,除尘风机

将风机模型导入ICEM 进行网格划分，网格划分过程中对离心风机关键部位要进行加密处理，如叶轮、集流器、蜗舌、进气箱的转角处等。对风机的进口与出口适当延长，以保证计算的稳定性。考虑到离心风机结构的复杂且不规则性，本文采用非结构四面体网格进行划分，其中无进气箱的离心风机网格数量约370万，网格质量为0.3以上；带进气箱的离心风机网格数量为380万，网格质量为0.3以上。

风机采用标准k-ɛ模型，壁面函数为Scalable，数值计算方法为高阶求解格式，求解格式为一阶格式。由于通风机转速低，马赫数小，可认为气流为不可压缩定常流动。进口给定质量流量，出口给定静压,壁面条件为无滑移边界，转速为1 480r/min，并将流动区域分为静止域与旋转域，两者通过Interface连接，连接模型为普通连接，坐标变换为转子算法，网格连接方式为GGI。本文所研究的某离心风机叶轮有均布的16 个前向的大小叶片，其内部流场较为复杂，为了揭示风机内的流场特性，对风机进行全三维数值模拟。先单独分析了进气箱内部流场特性，然后对进气箱与风机进行一体化分析，研究进气箱对离心风机性能的影响。

为改善风机受气体粘性影响导致流动分离加剧的现象，在传统蜗壳型线设计理论的基础上，研究气体粘性力矩对蜗壳壁线分布的影响，并采用动量矩修正方法对其进行改型设计。另外，为真实反映风机内流场分布情况，在标准k-ε 计算模型的扩散项中加入粘性应力作用，使其高计算误差降低至3%。对比分析改型前后风机数值模拟计算和试验测量结果可知，采用修改的k-ε 模型进行计算发现改型后风机内旋涡强度减小，蜗壳出口靠近蜗舌处流动分离得到改善。试验结果表明：改型风机出口静压提升约25Pa，较大全压效率较原型机提升约10%。

同时，由于蜗壳张开度扩大能够抑制流动分离，使蜗舌附近区域的旋涡强度及其影响区域减小，从而有效地降低了多翼离心风机噪声2.5dB。多翼离心风机广泛应用于国民经济的各个领域，是工业生产中主要耗能设备之一，蜗壳作为离心风机中不可或缺的基本元件，其结构的不对称性及内部流动的复杂性会对叶轮出口气流角造成较大影响，使其沿圆周方向呈现出明显的不对称性。而在风机实际运行过程中，风机叶轮出口气流与蜗壳壁面间存在强烈的非定常干涉，使得蜗壳壁面成为风机的主要噪声源。因此提高蜗壳型线设计水平，不仅能改善风机气动性能，还能达到降低噪声的效果。目前国内外学者对离心风机蜗壳型线的研究，主要集中在寻找能真实反映蜗壳内流体流动状态的设计方法。

将建立好的风机三维模型导入ICEM 软件进行混合网格的划分。其中进出口和叶轮区域采用结构化网格，而蜗壳部分由于其内部结构复杂，尤其是电动机周围结构并非规则模型，故采用适应性较强的非结构化四面体网格，具体网格如图3 所示。综合考虑动静耦合区域对数值模拟预测结果的影响，在进行网格划分时，对边界层进行加密处理，其较低网格质量雅克比[14]在0.3 以上。为了保证数值计算结果的准确性，避免网格误差对其模拟结果造成影响，对风机进行网格无关性验证，如表1 所示。综合考虑计算精度和计算效率可知，当网格数为25 万左右时预测结果较为合理，终确定整个计算域的网格数为2513558。k-ε 模型作为为普遍有效的湍流模型，能够计算大量的各种回流和薄剪切层流动，被广泛应用于各类风机的数值求解计算中。

由于有梯度扩散项，模型k-ε 方程为椭圆形方程，故其特性同其他椭圆形方程，需要边界条件：风机出口或对称轴处k / n0和/ n0。但上述边界条件只针对高雷诺数而言，在固体壁面附近，流体粘性应力将取代湍流雷诺应力，并在临近固体壁面的粘性底层占主要作用。而多翼离心风机由于结构尺寸小、相对马赫数低，气体黏性力在流体流动过程中起重要作用，因此，在实际运用过程中，标准k-ε 模型由于未充分考虑粘性力的影响，导致计算模型出现偏差。运用Visual C++将上述修正函数编写为UDF代码，并导入Fluent 内置Calculation module。为符合实际运行状态，风机进出口边界条件设置为压力入口和压力出口，出口压降与动能成正比，从而避免在进口和出口定义一致的速度分布[15]。后以CFD 计算的定常结果作为初始条件，进行非定常数值计算。