Calculation of Rotordynamic Coefficients and Leakage for Annular Gas Seals by Means of Finite Difference Techniques

1989 ◽  
Vol 111 (3) ◽  
pp. 545-552 ◽  
Author(s):  
R. Nordmann ◽  
F. J. Dietzen ◽  
H. P. Weiser

The compressible flow in a seal can be described by the Navier-Stokes equations in connection with a turbulence model (k–ε model) and an energy equation. By introducing a perturbation analysis in these differential equations we obtain zeroth order equations for the centered position and first order equations for small motions of the shaft about the centered position. These equations are solved by a finite difference technique. The zeroth order equations describe the leakage flow. Integrating the pressure solution of the first order equations yields the fluid forces and the rotordynamic coefficients, respectively.

1987 ◽  
Vol 109 (3) ◽  
pp. 388-394 ◽  
Author(s):  
F. J. Dietzen ◽  
R. Nordmann

For modelling the turbulent flow in a seal the Navier-Stokes equations in connection with a turbulence model (k-ε-model) are solved by a finite-difference method. A motion of the shaft around the centered position is assumed. After calculating the corresponding flow field and the pressure distribution, the rotordynamic coefficients of the seal can be determined. These coefficients are compared with results obtained by using the bulk flow theory of Childs [1] and with experimental results.


1992 ◽  
Vol 114 (3) ◽  
pp. 431-438 ◽  
Author(s):  
F. Simon ◽  
J. Freˆne

An analysis is developed to calculate the static and dynamic characteristics of annular eccentric seals. Effects of inertia forces in the film, tapered geometry and rotor misalignment are taken into account. Derivation of the governing equations for incompressible flow is based on the Navier-Stokes equations, the continuity equation and a turbulence model using the nonlinear analysis developed by Elrod and Ng. The inlet boundary conditions define the initial swirl and the pressure drop due to the fluid acceleration. Perturbation of the flow variables yields a set of zeroth-order and first-order equations. Integration of the zeroth-order equations yields the steady-state solution which defines the seal leakage, the static load and the moment of misalignment. The eccentric and misaligned rotordynamic coefficients are obtained by integration of the first-order pressure equations. Comparisons are made between the stiffness, damping and inertia coefficients derived herein and both theoretical results based on other models and experimental data which were previously published.


1997 ◽  
Vol 119 (3) ◽  
pp. 556-567 ◽  
Author(s):  
Mihai Arghir ◽  
Jean Freˆne

The paper presents a method to calculate the rotordynamic coefficients of circumferentially-grooved liquid seals operating in centered position and turbulent flow regimes. The method is based on the integration of the averaged Navier-Stokes equations and uses a coordinate transformation proposed by Dietzen and Nordmann (1987). The effect of the coordinate transformation on the components of the stress tensor is included in the first order transport equations. To ensure grid independent solutions, numerical boundary conditions for the first-order velocities were formulated using the logarithmic law. The perturbation of the turbulent viscosity was also considered. A pressure recovery effect at the exit section was included in the first order mathematical model. The method is validated by calculations for straight and circumferentially-grooved seals. Comparisons with experimental and theoretical results show a good agreement for straight seals and for seals with few grooves, and a reasonable agreement for severe industrial cases (high Reynolds numbers and large number of grooves).


Author(s):  
Joachim Schettel ◽  
Martin Deckner ◽  
Klaus Kwanka ◽  
Bernd Lu¨neburg ◽  
Rainer Nordmann

The main goal of this paper is to improve identification methods for rotordynamic coefficients of labseals for turbines. This aim was achieved in joint effort of the Technische Universita¨t Mu¨nchen, working on experimental identification methods for rotordynamic coefficients, the University of Technology, Darmstadt, working on prediction methods, and Siemens AG, realizing the results. The paper focuses on a short comb-grooved labyrinth seal. Short labseals, amongst others the above mentioned comb-grooved labyrinth, were examined. by means of a very accurately measuring test rig. The rotor was brought into statically eccentric positions relative to the stator, in order to measure the circumferential pressure distribution as a function of pressure, rotating speed and entrance swirl. The data collected were used to validate results obtained with a numerical method. The theoretical approach is based on a commercial CFD tool, which solves the Navier Stokes equations using numerical methods. As a result, a detailed model of the flow within the test rig is produced. The efforts of computation here are greater than when compared with the likewise wide-spread Bulk flow models, however improved accuracy and flexibility is expected. As the validation of the model is successful, it could then be used to gain further insight in the flow within the seal, and to understand the results better. This showed that rotordynamic coefficients of labseals gained from different test rigs are not necessarily comparable.


1976 ◽  
Vol 78 (2) ◽  
pp. 355-383 ◽  
Author(s):  
H. Fasel

The stability of incompressible boundary-layer flows on a semi-infinite flat plate and the growth of disturbances in such flows are investigated by numerical integration of the complete Navier–;Stokes equations for laminar two-dimensional flows. Forced time-dependent disturbances are introduced into the flow field and the reaction of the flow to such disturbances is studied by directly solving the Navier–Stokes equations using a finite-difference method. An implicit finitedifference scheme was developed for the calculation of the extremely unsteady flow fields which arose from the forced time-dependent disturbances. The problem of the numerical stability of the method called for special attention in order to avoid possible distortions of the results caused by the interaction of unstable numerical oscillations with physically meaningful perturbations. A demonstration of the suitability of the numerical method for the investigation of stability and the initial growth of disturbances is presented for small periodic perturbations. For this particular case the numerical results can be compared with linear stability theory and experimental measurements. In this paper a number of numerical calculations for small periodic disturbances are discussed in detail. The results are generally in fairly close agreement with linear stability theory or experimental measurements.


Sign in / Sign up

Export Citation Format

Share Document