A mathematical model of calculating rotordynamic coefficients associated with leakage steam flow through labyrinth seals was presented. Particular attention was given to incorporating thermal properties of the steam fluid into prediction of leakage flow and subsequent derivation of rotordynamic coefficients, which quantitatively characterize influence of aerodynamic forcing of the leakage steam flow on the rotordynamics. By using perturbation analysis, we determined periodic and analytic solutions of the continuity and circumferential momentum equations for the time-dependent flow induced by non-axisymmetric rotation of the rotor encompassed by a labyrinth seal. Pressure distributions along labyrinth seal cavities and rotordynamic coefficients were compared at the same condition for air and steam flows. Influence of steam flow through the labyrinth seal cavities on rotordynamic coefficients was analyzed in terms of inlet pressure, inlet swirl velocity and rotor speed.
List of symbolA unsteady cross-sectional area of the cavity [m 2 ] A 0 steady annular flow area [m 2 ] B tooth height [m] C 0 orifice contraction coefficient C 1 kinetic energy carry-over coefficient C r steady radial clearance [m] C x x , C yy damping coefficients [N s/m] C xy , C yx cross-couple damping coefficients [N s/m] Dh 0 , Dh steady and unsteady hydraulic diameter of the cross-sectional area of the cavity [m] F total air reaction force [N] F x , F y annular seal reaction force components [N] h 0 stagnant enthalpy [J/Kg] h iseal enthalpy at the ith orifice [J/Kg] H unsteady labyrinth seal radial clearance [m] H 1 (t, θ) perturbation clearance [m] K x x , K yy direct stiffness coefficients [N/m] K xy , K yx cross-couple stiffness coefficients [N/m] L cavity length [m] N tooth number n specific heat ratio P 0i , P i steady and unsteady pressure at the ith cavity [Pa] P 0iseal , P iseal steady and unsteady pressure at the ith orifice [Pa] P 1i (t, θ) perturbation pressure in the ith cavity [Pa] q 0i ,q i steady and unsteady leakage flow rate per unit length [kg/s m] q 1i (t, θ) perturbation leakage flow rate per unit length [kg/s m] R steam constant [J/kg K] R s rotor radius [m] T i temperature at the ith cavity [K] T iseal temperature at the ith orifice [K] U iseal velocity at the ith orifice [m/s] V 0i , V i steady and unsteady circumferential velocities at the ith cavity [m/s] V 1i (t, θ) perturbation pressure at the ith cavity [m/s] (x, y)x and y coordinates of rotor whirling motionGreek symbols ε eccentricity ratio η imaginary number ρ i , ρ 0i steady and unsteady gas density of the ith cavity [kg/m 3 ] ρ 1 (t, θ) perturbation gas density of the ith cavity [kg/m 3 ] ρ 0iseal , ρ iseal steady and unsteady gas density at the ith clearance [kg/m 3 ] ρ 1iseal (t, θ) perturbation gas density at the ith clearance [kg/m 3 ] τ r 0i , τ ri steady and unsteady shear stresses at the rotor wall [N/m 2 ] τ r 1i (t, θ) perturbation shear stresses at the rotor wall [N/m 2 ] τ s0i , τ si steady and unsteady shear stresses at the stator wall [N/m 2 ] ...