干式罗茨真空泵吸气级内流动的瞬态模拟

　　基于动网格方法建立了干式真空泵罗茨型吸气级的三维瞬态数值计算模型。模拟结果与抽速曲线对比，入口压力为1000Pa时误差为11.5%，100Pa时误差为34.2%，表明计算流体力学(CFD)方法不适用于入口压力较低及极限真空时真空泵内的流动研究，但在入口压力较高时具有较好的数值精度。

　　由于干式真空泵的主要设计问题多集中于入口压力较高，负荷较大的运行工况，应用CFD方法研究干式真空泵的流动特性具有实用价值。文中计算了真空泵的性能参数，分析了泵腔内的流动现象和流场的主要特征。

　　干式真空泵作为核心装备广泛应用于半导体产业、太阳能电池制造、化工、制药等重要行业，其研发和生产越来越受到重视。干式真空泵有多种结构形式，其中主要有多级罗茨式、多级爪式、组合式(罗茨+爪型)和螺杆式等，抽气过程中气体的温度和压力存在复杂的变化。目前，除制造装配要求精度高、加工难度大等制约因素外，设计理论也亟待提高，如转子热变形及对间隙的影响、空气动性噪声的产生与控制、颗粒物的抽除及在腔内的运动等问题都值得深入分析，因此，对泵腔的流场进行研究具有重要意义。

　　近年来，计算流体力学(CFD，Computational Fluid Dynamics)发展迅速并被应用于流体机械内部流场的研究。在真空技术领域也不断有CFD的应用成果被报道，Boulon和Cheng等研究了复合分子泵Holweck排气级和涡轮增压泵内的流动模拟问题；Wang等对不同运转条件下蒸汽喷射泵的性能进行了数值计算，发布此是文以多级干式真空泵的罗茨型吸气级泵腔为研究对象，基于CFD通用软件ANSYS-CFX并采用动网格方法对泵内的瞬态流动进行了数值计算，讨论了CFD方法的适用性，分析了气体的流动现象和流场分布，为进一步研究干式真空泵内的物理过程和理论问题进行了有益的尝试。

1、物理模型和数值模拟

　　1.1、物理模型

　　多级干式真空泵内的流动状态比较复杂，除粘滞流外，在直排大气的排气级可能出现湍流，而在压力较低的吸气级泵腔的狭小间隙内，可能存在滑移流、过渡流甚至分子流态。气体流态判别依据由分子平均自由程和几何特征尺寸决定的克努森数Kn，通常CFD方法适用于分析湍流、粘滞流，在滑移流的模拟中也有成功的算例，而在Kn>011的过渡流或分子流时，由于Navier-Stokes方程基于的连续性假设不再成立，理论上不适用。在高真空稀薄气体的模拟中主要采用直接蒙特卡罗(DSMC，Directsimulation Monte Carlo)方法，但随着压力升高，分子间碰撞频率迅速增加需要巨大的计算资源支撑，计算效率较低。值得注意的是根据克努森数Kn选择CFD或DSMC方法并没有明确的界限。基于以下事实:

　　(1)多级干式真空泵中克努森数较大的分子流和过渡流zui可能出现在压力较低的吸气级泵腔间隙处，泵内其它大部分流动空间Kn较小。

　　(2)由于吸气通道流导较大，多级干式真空泵的吸气级选择罗茨型泵腔。

　　(3)多级干式真空泵泵腔几何空间及容积变化过程复杂，整体建模存在一定困难。

　　本文采用CFD方法并以一实际干式罗茨真空泵的吸气级泵腔为对象展开研究，其主要设计参数如表1所示。

表1 主要设计参数

　　根据以上数据，转子与泵腔内壁的设计间隙为0.25mm，转子与转子在轴心连线上的啮合间隙为0.32mm。

　　1.2、动网格方法

　　罗茨转子型线复杂，转动时泵腔的容积和形状不断变化，需要采用动网格方法建立模型。CFX中直角坐标系下动网格的积分形式守恒方程为

　　式中，Q为流体密度；U为流体的速度；W为控制体边界移动速度；Leff为包含湍流粘度的粘度系数；„eff为包含湍流扩散的扩散系数；U为标量，如单位质量工质的焓值或内能；SU、SU为源项；dn为表面的外法向量；i，j代表直角坐标系中的坐标轴。动网格求解应满足几何守恒律