汽轮机高低齿迷宫式隔板汽封内部流动的数值研究.docVIP

2005 Fluent 中国用户大会论文集 汽轮机高低齿迷宫式隔板汽封内部流动的数值研究 李 军，丰镇平 （西安交通大学能源与动力工程学院?叶轮机械研究所，西安 710049） 摘 要：针对汽轮机高低齿迷宫式隔板汽封的复杂结构和轴的旋转效应，采用CFD商用软件FLUENT数值求解三维Navier-Stokes方程技术研究其内部流动。采用有限体积方法离散Navier-Stokes方程，标准k?ε两方程紊流模型封闭方程组。数值求解方法采用SIMPLE算法，迎风2阶格式求解对流项，迎风1阶格式求解扩散项。研究分析整体加工直齿和镶片式尖齿两种高低齿迷宫式隔板汽封的流动特性和泄漏量。数值模拟结果证明了齿间距增大，泄漏量相应增加，凸台的位置也影响泄漏量。数值模拟得到的流场揭示了高低齿迷宫式隔板汽封内的复杂流动特性。本研究工作为工程设计高低齿迷宫式隔板汽封提供了理论依据和技术支撑。 关键词：汽轮机 高低齿 迷宫式隔板汽封 数值模拟 1． 前言 现代火力发电厂汽轮机汽封一般采用梳齿迷宫式汽封，迷宫式汽封的密封原理是齿与转子间形成的一系列节流间隙与膨胀空腔，使通过的工作蒸汽产生节流与热力学效应而达到密封效果。汽封按其装置位置的不同可分为轴端汽封（轴封），隔板汽封和通流部分汽封（叶顶汽封）三类。对于隔板汽封，在汽轮机的高压端，缸内蒸汽压力高，为减少蒸汽的泄漏量，一般采用高低齿式汽封。在低压端，常采用光轴汽封，以适应转子和汽缸较大差胀的需要[1]。 随着节能和现代火电技术的发展，超临界汽轮机的应用范围越来越广泛，蒸汽压力高，级间密度大，造成漏气损失增加，而为了防止或减缓汽封片的磨损，又不能选择过小的间隙，因此合适的高低齿迷宫式隔板汽封的设计，既可以减少漏汽损失，又可以避免机组运行过程中发生转子低频振动，提高轴系稳定性[2,3]。 1由于高低齿隔板汽封的复杂性和转子的高速旋转造成其内部的三维紊流，使得实验研究其内部流动特性变得非常困难。随着计算机技术和计算流体动力学水平的发展，数值研究汽封内部流动，提高汽封设计水平，逐渐在汽轮机制造业受到重视而得到相应的发展。蔡虎等人采用数值方法研究了叶顶汽封的流动对叶片主流流场影响情况[4]。Schlienger等人采用实验方法研究了两种不同形状的叶顶汽封内部流动特性，主要研究分析了叶顶汽封内的泄漏流对叶片流道内二次流影响的特性[5]。Denecke等人采用实验方法研究了动叶叶顶汽封齿与沟槽的形状和相对位置对泄漏量的影响[6]。但是对于具有高速旋转的高低齿迷宫式隔板汽封内部流场的数值研究的报道却较少。 本文采用CFD商用软件FLUENT数值方法研究了汽轮机高压整体加工直齿和镶片式尖齿两种不同结构高低齿迷宫式隔板汽封的三维内部流场。采用有限体积方法离散Navier-Stokes方程，k?ε两方程紊流模型封闭方程组。数值求解方法采用SIMPLE算法，迎风2阶格式求解对流项，迎风1阶格式求解扩散项。数值模拟结果证明了齿间距增大，泄漏量相应增加，凸台的位置对泄漏量也有影响。数值模拟得到的流场揭示了高低齿迷宫式隔板汽封内的复杂流动特性。本研究工作为工程设计高低齿迷宫式隔板汽封提供了理论依据和技术支撑。 作者简介：李军（1971年3月），男，安徽省六安市生，副教授，工学博士，研究方向是叶轮机械气动热力学 基金项目：国家自然科学基金资助项127 2005 Fluent 中国用户大会论文集 2． 计算模型和数值方法 图1给出了某高低齿迷宫式隔板汽封一段弧段的三维结构示意图。其中s是间隙，即齿尖与转轴之间的径向距离，t是齿间距，w是凸台高。高低齿隔板汽封内的泄漏流动是由于隔板前后压力差造成的。因为转轴的高速旋转和迷宫式汽封结构的复杂性，使得高低齿隔板汽封内的流场是典型的三维紊流。为了正确分析隔板汽封内部流场特性和计算泄漏量，数值求解三维 Navier-Stokes方程是需要的。 图1 三维迷宫式高低齿隔板汽封结构示意图和计算网格图 vdivρUφ=divΓφgradφ+Sφ 三维隔板汽封内部流动的控制方程如下描述 ()() 式中，ρ为流体密度；U为流体速度矢量；φ为通用变量，可以代表速度u，v，T以及k和ε 等求解变量；Γφ为广义扩散系数；Sφ为广义源项。 为了节省计算时间，我们采用计算一段弧段的方法分析整个隔板汽封内部流场的办法，弧段的两个断面是周期性边界。流体是过热蒸汽。由于计算区域的复杂性，采用非结构化计算网格（如图1所示）。采用有限体积方法离散控制方程，对流项采用2阶迎风格式进行求解，扩散项采用1阶迎风格式进行数值求解。k?ε两方程紊流模型用于封闭求解的方程组，采用SIMPLE算法对离散方程进行数值求解[7]。计算区域中转轴的转速是3000rpm，模拟中将转轴壁面设计成旋转壁面。 v 3． 结果分