Labyrinth Seal Rotordynamic Forces Using a Three-Dimensional Navier-Stokes Code

1992 ◽  
Vol 114 (4) ◽  
pp. 683-689 ◽  
Author(s):  
D. L. Rhode ◽  
S. J. Hensel ◽  
M. J. Guidry

A finite difference method for determining rotordynamic forces on an eccentric whirling labyrinth cavity has been developed. A coordinate-transformation was applied to the Reynolds time-averaged Navier-Stokes equations in order to use the modified bipolar coordinate system. The SIMPLER algorithm with QUICK differencing and the high Reynolds number k–ε turbulence model are used to compute the complex turbulent flowfield. A circular whirl orbit about the geometric center of the housing was specified for simplicity. The new model was tested against the rotordynamic force measurements, and close agreement was found. For the cases considered, the radial and tangential force components become rotor dynamically less desirable with increasing inlet swirl. Also, circumferential pressure variations are included for enhanced insight into the flowfield.

Author(s):  
J. Jeffrey Moore

Abstract Labyrinth seals are utilized inside turbomachinery to provide non-contacting control of internal leakage. These seals can also play an important role in determining the rotordynamic stability of the machine. Traditional labyrinth seal models are based on bulk-flow assumptions where the fluid is assumed to behave as a rigid body affected by shear stress at the interfaces. To model the labyrinth seal cavity, a single, driven vortex is assumed and relationships for the shear stress and divergence angle of the through flow jet are developed. These models, while efficient to compute, typically show poor prediction for seals with small clearances, high running speed, and high pressure (Childs, 1993). In an effort to improve the prediction of these components, this work utilizes three-dimensional computational fluid dynamics (CFD) to model the labyrinth seal flow path by solving the Reynolds Averaged Navier Stokes equations. Unlike bulk-flow techniques, CFD makes no fundamental assumptions on geometry, shear stress at the walls, as well as internal flow structure. The method allows modeling of any arbitrarily shaped domain including stepped and interlocking labyrinths with straight or angled teeth. When only leakage prediction is required, an axisymmetric model is created. To calculate rotordynamic forces, a full 3D, eccentric model is solved. The results demonstrate improved leakage and rotordynamic prediction over bulk-flow approaches compared to experimental measurements.


2003 ◽  
Vol 125 (4) ◽  
pp. 427-433 ◽  
Author(s):  
J. Jeffrey Moore

Labyrinth seals are utilized inside turbomachinery to provide noncontacting control of internal leakage. These seals can also play an important role in determining the rotordynamic stability of the machine. Traditional labyrinth seal models are based on bulk-flow assumptions where the fluid is assumed to behave as a rigid body affected by shear stress at the interfaces. To model the labyrinth seal cavity, a single, driven vortex is assumed and relationships for the shear stress and divergence angle of the through flow jet are developed. These models, while efficient to compute, typically show poor prediction for seals with small clearances, high running speed, and high pressure.* In an effort to improve the prediction of these components, this work utilizes three-dimensional computational fluid dynamics (CFD) to model the labyrinth seal flow path by solving the Reynolds Averaged Navier Stokes equations. Unlike bulk-flow techniques, CFD makes no fundamental assumptions on geometry, shear stress at the walls, as well as internal flow structure. The method allows modeling of any arbitrarily shaped domain including stepped and interlocking labyrinths with straight or angled teeth. When only leakage prediction is required, an axisymmetric model is created. To calculate rotordynamic forces, a full 3D, eccentric model is solved. The results demonstrate improved leakage and rotordynamic prediction over bulk-flow approaches compared to experimental measurements.


1997 ◽  
Vol 119 (4) ◽  
pp. 900-905 ◽  
Author(s):  
X. Zheng ◽  
C. Liao ◽  
C. Liu ◽  
C. H. Sung ◽  
T. T. Huang

In this paper, computational results are presented for three-dimensional high-Reynolds number turbulent flows over a simplified submarine model. The simulation is based on the solution of Reynolds-Averaged Navier-Stokes equations and two-equation turbulence models by using a preconditioned time-stepping approach. A multiblock method, in which the block loop is placed in the inner cycle of a multi-grid algorithm, is used to obtain versatility and efficiency. It was found that the calculated body drag, lift, side force coefficients and moments at various angles of attack or angles of drift are in excellent agreement with experimental data. Fast convergence has been achieved for all the cases with large angles of attack and with modest drift angles.


2003 ◽  
Vol 478 ◽  
pp. 227-235 ◽  
Author(s):  
J. D. GIBBON ◽  
Charles R. DOERING

Dissipation-range intermittency was first observed by Batchelor & Townsend (1949) in high Reynolds number turbulent flows. It typically manifests itself in spatio-temporal binary behaviour which is characterized by long, quiescent periods in the signal which are interrupted by short, active ‘events’ during which there are large excursions away from the average. It is shown that Leray's weak solutions of the three-dimensional incompressible Navier–Stokes equations can have this binary character in time. An estimate is given for the widths of the short, active time intervals, which decreases with the Reynolds number. In these ‘bad’ intervals singularities are still possible. However, the average width of a ‘good’ interval, where no singularities are possible, increases with the Reynolds number relative to the average width of a bad interval.


Author(s):  
V. Tamimi ◽  
M. Zeinoddini ◽  
A. Bakhtiari ◽  
M. Golestani

In this paper results from simulating the vortex shedding phenomena behind a fixed tapered circular cylinder, at relatively high Reynolds numbers, are reported. Ansys-CFX computational fluid dynamics model, based on solving three-dimensional (3D) incompressible transient Navier Stokes equations, is employed for this purpose. The geometries applied in the models resemble those used in wind tunnel experiments by other researchers. The taper slope along the cylinder span is uniform with a tangent of 24:1. The diameter at mid-span of the cylinder equals to 0.0389 m. The Reynolds number (based on the mid-span diameter) is around 29,000. The computational model has first been calibrated against experiments for uniform 3D cylinders as well as results from a Direct Numerical Simulation of turbulent wake with vortex shedding past a uniform circular cylinder, as obtained by other researchers. The main flow characteristics for tapered cylinders such as vortex dislocations and splitting, cellular vortex shedding, oblique vortex shedding and the variation of the vorticity patterns along the tapered cylinder could be obtained from the simulations.


2009 ◽  
Vol 628 ◽  
pp. 357-369 ◽  
Author(s):  
J. D. CROUCH ◽  
A. GARBARUK ◽  
D. MAGIDOV ◽  
A. TRAVIN

Buffeting flow on transonic aerofoils serves as a model problem for the more complex three-dimensional flows responsible for aeroplane buffet. The origins of transonic aerofoil buffet are linked to a global instability, which leads to shock oscillations and dramatic lift fluctuations. The problem is analysed using the Reynolds-averaged Navier–Stokes equations, which for the foreseeable future are a necessary approximation to cover the high Reynolds numbers at which transonic buffet occurs. These equations have been shown to reproduce the key physics of transonic aerofoil flows. Results from global-stability analysis are shown to be in good agreement with experiments and numerical simulations. The stability boundary, as a function of the Mach number and angle of attack, consists of an upper and a lower branch – the lower branch shows features consistent with a supercritical bifurcation. The unstable modes provide insight into the basic character of buffeting flow at near-critical conditions and are consistent with fully nonlinear simulations. The results provide further evidence linking the transonic buffet onset to a global instability.


2020 ◽  
Vol 14 (4) ◽  
pp. 7369-7378
Author(s):  
Ky-Quang Pham ◽  
Xuan-Truong Le ◽  
Cong-Truong Dinh

Splitter blades located between stator blades in a single-stage axial compressor were proposed and investigated in this work to find their effects on aerodynamic performance and operating stability. Aerodynamic performance of the compressor was evaluated using three-dimensional Reynolds-averaged Navier-Stokes equations using the k-e turbulence model with a scalable wall function. The numerical results for the typical performance parameters without stator splitter blades were validated in comparison with experimental data. The numerical results of a parametric study using four geometric parameters (chord length, coverage angle, height and position) of the stator splitter blades showed that the operational stability of the single-stage axial compressor enhances remarkably using the stator splitter blades. The splitters were effective in suppressing flow separation in the stator domain of the compressor at near-stall condition which affects considerably the aerodynamic performance of the compressor.


Processes ◽  
2021 ◽  
Vol 9 (5) ◽  
pp. 841
Author(s):  
Yuzhen Jin ◽  
Huang Zhou ◽  
Linhang Zhu ◽  
Zeqing Li

A three-dimensional numerical study of a single droplet splashing vertically on a liquid film is presented. The numerical method is based on the finite volume method (FVM) of Navier–Stokes equations coupled with the volume of fluid (VOF) method, and the adaptive local mesh refinement technology is adopted. It enables the liquid–gas interface to be tracked more accurately, and to be less computationally expensive. The relationship between the diameter of the free rim, the height of the crown with different numbers of collision Weber, and the thickness of the liquid film is explored. The results indicate that the crown height increases as the Weber number increases, and the diameter of the crown rim is inversely proportional to the collision Weber number. It can also be concluded that the dimensionless height of the crown decreases with the increase in the thickness of the dimensionless liquid film, which has little effect on the diameter of the crown rim during its growth.


Sign in / Sign up

Export Citation Format

Share Document