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 ] ...
An extensive investigation of the influence of the leakage flow through a labyrinth seal at supply pressure of 12 bar on the rotordynamics was performed by using numerical calculations and experimental measurements. Toward this end, an experimental rotor setup was established in Shanghai Jiao Tong University. Two labyrinth seals were chosen for comparison, e.g., an interlocking seal and a stepped one. The numerical calculations based on the bulk-flow theory and the perturbation analysis were accomplished. Simultaneous acquisitions of the fluctuating static pressure at the stator wall and the displacement of the whirling rotor were made. The influence of the aerodynamic forcing on the rotor was analyzed in terms of the axial distribution of the mean static pressure, the circumferential distribution of the fluctuating pressure, the fist critical speed and the destabilization rotating speed of the rotor. The experimental results demonstrated that the sinusoidal distribution of the fluctuating static pressure on the stator wall was closely related to the whirling motion of the rotor. The first critical speed of the rotor was reduced by the aerodynamic forcing, resulting in intensified destabilization of the rotor system. Furthermore, the numerical analyses were in good agreement to the experimental measurements. List of symbols AUnsteady cross-sectional area of the cavity (m 2 ) A n Steady annular flow area (m 2 ) [C]Damping matrix C 0Orifice contraction coefficient C 1 Kinetic energy carry-over coefficient C r Steady radial seal clearance (m) C x x , C yy Damping coefficients (N-s/m) C xy , C yx Cross-couple damping coefficients (N-s/m)
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.