Thermal Details in a Rotor–Stator Cavity at Engine Conditions With a Mainstream

1992 ◽  
Vol 114 (2) ◽  
pp. 446-453 ◽  
Author(s):  
S. H. Ko ◽  
D. L. Rhode

This investigation involves a numerical study of enclosed Rotor–Stator cavities of gas turbine engines. The complete elliptic form of the 2-D, axisymmetric Navier–Stokes equations for compressible turbulent flow were solved. Included are the complete fluid and thermal effects of the hot mainstream gas interacting with the cooling cavity purge flow at actual engine flow conditions for generalized geometries. Additional flow conditions above and below those for engine nominal conditions are also considered. The relationships among the important flow parameters are investigated by examining the entire set of computations. The predictions reveal that a small recirculation zone in the stator shroud axial gap region is the primary mechanism for the considerable thermal transport from the mainstream to the turbine blade root/retainer region of the rotor.

Author(s):  
S. H. Ko ◽  
D. L. Rhode

This investigation involves a numerical study of enclosed rotor-stator cavities of gas turbine engines. The complete elliptic form of the 2-D, axisymmetric Navier-Stokes equations for compressible turbulent flow were solved. Included are the complete fluid and thermal effects of the hot mainstream gas interacting with the cooling cavity purge flow at actual engine flow conditions for generalized geometries. Additional flow conditions above and below that for engine nominal conditions are also considered. The relationships among the important flow parameters are investigated by examining the entire set of computations. The predictions reveal that a small recirculation zone in the stator shroud axial gap region is the primary mechanism for the considerable thermal transport from the mainstream to the turbine blade root/retainer region of the rotor.


2005 ◽  
Vol 127 (4) ◽  
pp. 806-815 ◽  
Author(s):  
J. U. Schlüter ◽  
X. Wu ◽  
S. Kim ◽  
S. Shankaran ◽  
J. J. Alonso ◽  
...  

Full-scale numerical prediction of the aerothermal flow in gas turbine engines are currently limited by high computational costs. The approach presented here intends the use of different specialized flow solvers based on the Reynolds-averaged Navier-Stokes equations as well as large-eddy simulations for different parts of the flow domain, running simultaneously and exchanging information at the interfaces. This study documents the development of the interface and proves its accuracy and efficiency with simple test cases. Furthermore, its application to a turbomachinery application is demonstrated.


Author(s):  
Martin P. King ◽  
Michael Wilson ◽  
J. Michael Owen

Buoyancy effects can be significant in the rotating annular cavities found between compressor discs in gas-turbine engines, where Rayleigh numbers above 1012 are common. In some engines, the cavity is ‘closed’, so that the air is confined between four rotating surfaces: two discs and inner and outer cylinders. In most engines, however, the cavity is ‘open’, and there is an axial throughflow of cooling air at the centre. For open rotating cavities, a review of the published evidence suggests a Rayleigh-Be´nard type of flow structure, in which, at the larger radii, there are pairs of cyclonic and anti-cyclonic vortices. The toroidal circulation created by the axial throughflow is usually restricted to the smaller radii in the cavity. For a closed rotating annulus, solution of the unsteady Navier-Stokes equations, for Rayleigh numbers up to 109, show Rayleigh-Be´nard convection similar to that found in stationary enclosures. The computed streamlines in the r-θ plane show pairs of cyclonic and anti-cyclonic vortices; but, at the larger Rayleigh numbers, the computed isotherms suggest that the flow in the annulus is thermally mixed. At the higher Rayleigh numbers, the computed instantaneous Nusselt numbers are unsteady and tend to oscillate with time. The computed time-average Nusselt numbers are in good agreement with the correlations for Rayleigh-Be´nard convection in a stationary enclosure, but they are significantly higher than the published empirical correlations for a closed rotating annulus.


2003 ◽  
Vol 3 ◽  
pp. 195-207
Author(s):  
A.M. Ilyasov ◽  
V.N. Kireev ◽  
S.F. Urmancheev ◽  
I.Sh. Akhatov

The work is devoted to the analysis of the flow of immiscible liquid in a flat channel and the creation of calculation schemes for determining the flow parameters. A critical analysis of the well-known Two Fluids Model was carried out and a new scheme for the determination of wall and interfacial friction, called the hydraulic approximation in the theory of stratified flows, was proposed. Verification of the proposed approximate model was carried out on the basis of a direct numerical solution of the Navier–Stokes equations for each fluid by a finite-difference method with phase-boundary tracking by the VOF (Volume of Fluid) method. The graphical dependencies illustrating the change in the interfase boundaries of liquids and the averaged over the occupied area of the phase velocities along the flat channel are presented. The results of comparative calculations for two-fluid models are also given, according to the developed model in the hydraulic approximation and direct modeling. It is shown that the calculations in accordance with the hydraulic approximation are more consistent with the simulation results. Thus, the model of hydraulic approximation is the most preferred method for calculating stratified flows, especially in cases of variable volumetric content of liquids.


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.


1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.


1996 ◽  
Vol 118 (2) ◽  
pp. 335-340 ◽  
Author(s):  
Mohamed Selmi

This paper is concerned with the solution of the 3-D-Navier-Stokes equations describing the steady motion of a viscous fluid inside a partially filled spinning and coning cylinder. The cylinder contains either a single fluid of volume less than that of the cylinder or a central rod and a single fluid of combined volume (volume of the rod plus volume of the fluid) equal to that of the cylinder. The cylinder rotates about its axis at the spin rate ω and rotates about an axis that passes through its center of mass at the coning rate Ω. In practical applications, as in the analysis and design of liquid-filled projectiles, the parameter ε = τ sin θ, where τ = Ω/ω and θ is the angle between spin axis and coning axis, is small. As a result, linearization of the Navier-Stokes equations with this parameter is possible. Here, the full and linearized Navier-Stokes equations are solved by a spectral collocation method to investigate the nonlinear effects on the moments caused by the motion of the fluid inside the cylinder. In this regard, it has been found that nonlinear effects are negligible for τ ≈ 0.1, which is of practical interest to the design of liquid-filled projectiles, and the solution of the linearized Navier-Stokes equations is adequate for such a case. However, as τ increases, nonlinear effects increase, and become significant as ε surpasses about 0.1. In such a case, the nonlinear problem must be solved. Complete details on how to solve such a problem is presented.


Author(s):  
David Gross ◽  
Yann Roux ◽  
Benjamin Rousse ◽  
François Pétrié ◽  
Ludovic Assier ◽  
...  

The problem of Vortex-Induced Vibrations (VIV) on spool and jumper geometries is known to present several drawbacks when approached with conventional engineering tools used in the study of VIV on risers. Current recommended practices can lead to over-conservatism that the industry needs to quantify and minimize within notably cost reduction objectives. Within this purpose, the paper will present a brief critical review of the Industry standards and more particularly focus on both experimental and Computational Fluid Dynamic (CFD) approaches. Both qualitative and quantitative comparisons between basin tests and CFD results for a 2D ‘M-shape’ spool model will be detailed. The results presented here are part of a larger experimental and numerical campaign which considered a number of current velocities, heading and geometry configurations. The vibratory response of the model will be investigated for one of the current velocities and compared with the results obtained through recommended practices (e.g. Shear7 and DNV guidelines). The strategy used by the software K-FSI to solve the fluid-structure interaction (FSI) problem is a partitioned coupling solver between fluid solver (FINE™/Marine) and structural solvers (ARA). FINE™/Marine solves the Reynolds-Averaged Navier-Stokes Equations in a conservative way via the finite volume method and can work on structured or unstructured meshes with arbitrary polyhedrons, while ARA is a nonlinear finite element solver with a large displacement formulation. The experiments were conducted in the BGO FIRST facility located in La Seyne sur Mer, France. Particular attention was paid towards the model design, fabrication, instrumentation and characterization, to ensure an excellent agreement between the structural numerical model and the actual physical model. This included the use of a material with low structural damping, the performance of stiffness and decay tests in air and in still water, plus the rationalization of the instrumentation to be able to capture the response with the minimum flow perturbation or interaction due to instrumentation.


1998 ◽  
Vol 14 (1) ◽  
pp. 23-29
Author(s):  
Robert R. Hwang ◽  
Sheng-Yuh Jaw

ABSTRACTThis paper presents a numerical study on turbulent vortex shedding flows past a square cylinder. The 2D unsteady periodic shedding motion was resolved in the calculation and the superimposed turbulent fluctuations were simulated with a second-order Reynolds-stress closure model. The calculations were carried out by solving numerically the fully elliptic ensemble-averaged Navier-Stokes equations coupled with the turbulence model equations together with the two-layer approach in the treatment of the near-wall region. The performance of the computations was evaluated by comparing the numerical results with data from available experiments. Results indicate that the present study gives good agreement in the shedding frequency and mean drag as well as in some phase profiles of the mean velocity.


Sign in / Sign up

Export Citation Format

Share Document