006_旋转机械空化磨损形成机理的CFD研究.pdf

旋转机械空化磨损形成机理的 CFD 研究 曾小康 周慧辉 昝元峰 熊万玉 （中国核动力研究设计院中核集团核反应堆热工水力技术重点实验室，成都，610041 ） 摘要：针对田湾 1 号和 2 号机组主循环泵运行一年后轴瓦表面上存在着类似于空化磨损的凹凸不平区 域，本文采用 CFD 方法研究了主循环泵径向止推轴承室内的压力场，指出了轴承室内形成低压区域的因素 为流动分离效应和流动加速效应。本文重点分析了静轴瓦楔形槽和动轴瓦引水槽重合时的空化现象及其形 成机理，确定了为避免空化现象的发生而提高轴承室进口静压的数值。本文为采用 CFD 方法对复杂转动机 械进行水力分析提供了很好的例子，其研究成果为田湾主循环泵轴承冷却回路头箱改造提供了依据。 关键词：旋转机械 空化磨损 CFD 机理分析 1 引 言 田湾核电站的一号机组和二号机组的主泵均采用了水润滑轴承，以水为润滑剂的水润滑 轴承具有以下诸多优点：无污染，绿色环保；来源广泛，价格低廉；水不燃烧，没有安全隐 患；较之油润滑轴承成本低，易维护保养；能够降低摩擦、振动和噪声；简化润滑轴承系统 的结构等。 图 1.静轴瓦瓦片及其磨损图 图2.动轴瓦瓦片及其磨损图 Fig.1 The abrasion of the actionless axletree Fig.2 The abrasion of the action axletree 但水润滑方式也带来了一些相应的技术难点和问题，比如：在田湾核电站的主泵运行约 一年后，检修人员发现主泵轴瓦上存在着一些类似于因汽蚀磨损而产生的、凹凸不平的非光 滑磨损面（见图 1 和图 2 ）。 在田湾主泵的轴承室内，因高速流动而形成低压区的情况可分为三种： （1）在动、静轴瓦间极窄的区域中，紧靠静轴瓦的流体处于低速流动状态，紧靠动轴 瓦的流体由于“拖曳力” 的作用处于高速流动状态，流体由静轴瓦流向动轴瓦时就会发生压力 突降； （2 ）在动轴瓦的导水槽中，由于动轴瓦的高速旋转使导水槽及其附近的流体产生相当 大的离心力，在离心力的驱动下，迫使流体从环形内腔加速进入导水槽中，在导水槽进口处 产生“流动分离效应”，引起压力突降； （3 ）在动轴瓦导水槽和静轴瓦楔形槽重合时，由于静轴瓦楔形槽沿半径方向、而动轴 瓦导水槽与半径方向存在角度，从而前者对离心力的阻力远小于后者，流体在重合部分产生 “加速效应”和“流动分离效应”，引起压力突降。 对于第一种情况，可近似地采用伯努利方程进行计算分析；但对于第二种情况和第三种 情况，由于几何结构和流场情况复杂，各种相关因素较多，必须采用 CFD 方法计算得到压 力场和速度场，才能分析空化发生的原因及空化的强度等问题，在此基础上给出具体的优化 建议。 2 计算思路及方法 本次计算的目的是获得主泵轴瓦流道几何结构、进口条件和旋转运动对空化形成的影 响，分析空化是否发生、发生的位置和大小。那么，为保证计算的可实现性，计算的技术难 点就集中在流道几何结构的简化和网格划分、旋转计算模型的选取和边界条件的设置。值得 注意的是，相对于一般的转动问题，本次计算的原型存在着厚度仅为几微米的液膜，几微米 的厚度超出了 CFD 方法能解析的最小尺度，但由于液膜为动静轴瓦间的极小缝隙而在计算 中不能略掉，所以必须对液膜进行处理，使其能满足 CFD 计算的要求又不能对液膜以外的 区域的流场压力场产生改变，液膜的处理是本次计算的关键技术难点。以上的四个技术难点 不是孤立存在，而是交叉在一起的，难点产生的根源主要来源于超薄液膜和复杂几何结构。 目前尚没有相似的算例可供参考，导致整个计算难度大大增加，因此本计算可作为一个带有 探索性的成功算例。 为了解决计算中的诸多技术难点，从而实现计算目的，我们制定了以下总体技术思路： （1）在几何模型简化上，对远离磨损区的流道结构进行简化，包括去除倒角和死水区 域、把不规则流道简化为规则流道、把轴瓦上下游的复杂流道都简化为圆环形流道等,静轴 瓦和动轴瓦的简化如图 3 和图 4 。 （2 ）在网格划分上，考虑到非结构化网格若要达到计算要求的精度，其网格数将远超 过工作站电脑的最大允许网格数，所以必须采用结构化网格，结构化网格容易控制数目和网 格精度，并对液膜也更容易处理。对全周期的网格，采用周期性网格（静轴瓦为 1/10 周期， 动轴