螺旋槽干氣密封微尺度效应的数值分析.docVIP

螺旋槽干气密封微尺度效应的数值分析 王兵 张会强 1 清华大学，航空航天学院，北京100084，中国 2 清华大学，北京3E能源重点实验室，北京100084，中国 2009年12月15日收稿；2010年11月修订；2010年11月录用；2010年11月电子出版 摘要：干气密封系统是在不同的工业应用中被广泛使用的非接触式密封技术。螺旋槽干气密封是利用流体动压效应实现密封和润滑过程，同时由于气体进出导致的减速而在两个密封面之间形成高压气膜。关于在微尺度上凹槽与气压膜对密封性能的影响的深入研究还很少。本文探讨微尺度对影响螺旋槽干气密封性能的修正雷诺方程的数值解。雷诺方程离散通过有限差分法与二阶计划的手段和逐次超松弛（SOR）迭代方法求解。密封气膜流Knudsen数从0.005到0.120，随着气膜深度和密封压力的变化。计算结果表明，由于微尺度的影响，气压膜与密封气体泄漏的平均压力增加，开启力度变大，而气膜刚度显着下降；摩擦扭矩和功率消耗保持不变，即使在密封压力与转速很低的条件下。在本文中，在不同转子面旋转速度的密封性能也被描述。本文研究阐明了微尺度对螺旋槽干气密封和密封性能的影响，预计将有助于改善设计干气密封系统运行中的滑移流动状态。 关键词：微尺度效应，螺旋槽，干气密封，数值分析，雷诺方程 1 引言 干气密封技术已被广泛应用于航空等各个行业，因为它的许多优点，包括低泄漏，低功耗，运行寿命长。运用流体动力学压力效应，干气密封系统仅仅一种非接触式密封技术[1]。非常薄的气膜在干气的转子和定子面之间形成密封系统，并具有润滑功能。一般来说，肩面沟槽的深度从3到10微米。研究人员发现，要获得高性能，螺旋槽是功能最全的肩面配置之一[2-3]。在操作过程中，当转子的转速大于干气密封系统的临界速度时，定子由于凹槽进入流体使压力增加的原因而升高。稳定的气膜，厚度从3到5微米，取决于密封压力，直到静态压力等于动态压力形成于转子与定子面之间。开启力度，摩擦功耗，密封介质泄漏，气膜刚度，除其他因素，在密封系统设计中必须考虑。[4-6] 影响密封性能的主要因素有槽深度，气膜厚度，转速，和密封介质的物理性质，化学性质。众所周知，气膜和凹槽只有几个微米的变化幅度，因此，在一定条件下，微尺度效果是显着的，并可以大大影响流动性，因此，密封性能好。然而，设计者可用的关于微尺度效应的信息却不足，这就产生了在干气密封系统中微效应研究的必要性。最近，阮[7]，采用有限元数值解雷诺方程，研究速度滑移条件对于密封性能的影响。福井等[8]，通过对泊肃叶流的速率校正衍生出雷诺方程并建成了干气密封的理论分析数据库。彭等[9]，提出了校正的润滑方程，适用于纳米流。他们分析了在硬盘磁头表面的摩擦学问题。然而，这些调查提供的关于微尺度对干气密封性能影响的信息有限。 关于干气密封的理论分析，研究人员一直在进行数值研究。通过雷诺方程能解决并获得干气膜的压力分布[10]。三种类型的数值计算方法，可求解雷诺方程：有限元（FE），有限体积（FV），和有限差分（FD）[11-15]。FE和FV方法虽然适合于复杂的几何形状，但在代码开发上存在缺点。相比之下，FD方法的计算机代码更容易开发。此外，FD方法的收敛速度和稳定性可以借助一个直接的解决方法或者逐次超松弛（SOR）迭代方法而增强。因此，FD方法对整个密封面压力场的分析很有效。 在本文中，通过雷诺方程的数值解分析，进行研究微尺度效应对干气密封系统中的流动性和的密封性能的影响。基于FD的快速迭代法可以通过泊肃叶流速率因素校正传统的雷诺方程。干气密封系统中主要表现的开启力度，泄漏，气膜刚度等参数，在压力场的数据的基础上进行计算。通过压力场和性能参数的比较，阐述微尺度效应对螺旋槽干气密封系统的影响。 2 方程和方法 2.1 物理模型和方程 图1所示的是干气密封旋转环物理模型示意图，环面开槽并在绕其轴旋转。在机械平衡条件下，流体填充凹槽的深度为h，填充转子与定子面之间的间隙的深度为δ。转子角速度为ω。轴的半径为r0。肩环的半径为r2.。槽根半径，如图1所示，是为r1。槽的面积比例为β=b1/b，在这里，b1是槽宽，b是槽区域的宽度。间隙比γ=（δ+h）。比密封压力是ε=p0/pi。在这里，p0是外部的压力或密封介质压力，pi是内部压力或环境介质的压力。 图1。物理模型示意图 为方便理论分析，做以下假设： （1）流体在气膜处服从等温理想气体模型。（2）流体的粘度是恒定的。（3）在气膜的流动是层流。（4）环面平稳，呈刚性。（5）转子无轴心差。 干气密封在圆柱坐标系的稳定方程写成 转子角速度ω =2πnr，nr是轴的旋转速度，m为流体的动力粘度。 上述的雷诺方程采用了简化的二维模型，在转子环面上由气膜深度分布函数表示，定义为：