Analysis on Flutter Characteristics of Transonic Compressor Blade Row by a Fluid-Structure Coupled Method

Author(s):  
Mingming Zhang ◽  
Anping Hou ◽  
Sheng Zhou ◽  
Xiaodong Yang

A time domain numerical approach is carried out to enhance the understanding of three dimensional blade row aeroelastic characteristics under the parallel computation. The vibration energy of unsteady aerodynamic force on the entire blade row is investigated using numerical solution of 3-D Navier-Stokes equations, coupled with structure finite element models for the blades to identify modal shapes and the structural deformations simultaneously. Interactions between fluid and structure are dealt with in a coupled manner, based on the interface information exchange until convergence in each time step. With this approach good agreement between the numerical results and the experimental data is observed. The flutter mechanism is analyzed according to deformation of the blades. The effect of inter-blade phase angle (IBPA) is included in the analysis by releasing the hypothesis of constant phase angle between adjacent blades in the traveling wave model. The results illustrate fully three dimensional unsteady nonlinear behaviors, such as limit-cycle oscillation. It is shown that all blades flutter at the same mode and frequency, but not at the same amplitude and IBPA. The analysis of the influence of different tip clearance gaps on the flutter characteristics of the blade row is also performed.

1999 ◽  
Vol 121 (4) ◽  
pp. 751-762 ◽  
Author(s):  
G. A. Gerolymos ◽  
I. Vallet

The purpose of this paper is to investigate tip-clearance and secondary flows numerically in a transonic compressor rotor. The computational method used is based on the numerical integration of the Favre-Reynolds-averaged three-dimensional compressible Navier–Stokes equations, using the Launder–Sharma near-wall k–ε turbulence closure. In order to describe the flowfield through the tip and its interaction with the main flow accurately, a fine O-grid is used to discretize the tip-clearance gap. A patched O-grid is used to discretize locally the mixing-layer region created between the jetlike flow through the gap and the main flow. An H–O–H grid is used for the computation of the main flow. In order to substantiate the validity of the results, comparisons with experimental measurements are presented for the NASA_37 rotor near peak efficiency using three grids (of 106, 2 X 106, and 3 X 106 points, with 21, 31, and 41 radial stations within the gap, respectively). The Launder–Sharma k–ε model underestimates the hub corner stall present in this configuration. The computational results are then used to analyze the interblade-passage secondary flows, the flow within the tip-clearance gap, and the mixing downstream of the rotor. The computational results indicate the presence of an important leakage-interaction region where the leakage-vortex after crossing the passage shock-wave mixes with the pressure-side secondary flows. A second trailing-edge tip vortex is also clearly visible.


Author(s):  
Araz Arbabi ◽  
Wahid Ghaly ◽  
Adam Medd

An aerodynamic inverse design method is developed for the simulation of three-dimensional viscous flow over blades, it is implemented into a commercial CFD program, namely ANSYS-CFX, and it is applied to the design of a transonic compressor stage. The implementation is validated for Rotor 37; it is then assessed in the redesign of Stage 67 stator. One set of design choices is to prescribe a target blade pressure loading and blade thickness distributions and a stacking line from hub to tip. The blade walls are assumed to be moving with a virtual velocity that would asymptotically drive the blade to the shape that would correspond to the specified target pressure distribution. This virtual velocity distribution is computed from the difference between the computed and the target pressure distributions. This inverse design approach is fully consistent with the viscous flow assumption and is independent of the CFD approach taken. The Arbitrary Lagrangian-Eulerian formulation of the unsteady Reynolds-Averaged Navier Stokes equations is solved in a time accurate fashion with the blade motion being the source of unsteadiness. At each time step, the blade shape is modified and dynamic meshing is used to remesh the fluid flow domain. To demonstrate the ability of this approach, it is applied to redesign the stator of a transonic axial fan, Stage 67, to improve its performance.


1987 ◽  
Vol 109 (1) ◽  
pp. 83-90 ◽  
Author(s):  
W. N. Dawes

The numerical analysis of highly loaded transonic compressors continues to be of considerable interest. Although much progress has been made with inviscid analyses, viscous effects can be very significant, especially those associated with shock–boundary layer interactions. While inviscid analyses have been enhanced by the interactive inclusion of blade surf ace boundary layer calculations, it may be better in the long term to develop efficient algorithms to solve the full three-dimensional Navier–Stokes equations. Indeed, it seems that many phenomena of key interest, like tip clearance flows, may only be accessible to a Navier–Stokes solver. The present paper describes a computer program developed for solving the three-dimensional viscous compressible flow equations in turbomachine geometries. The code is applied to the study of the flowfield in an axial-flow transonic compressor rotor with an attempt to resolve the tip clearance flow. The predicted flow is compared with laser anemometry measurements and good agreement is found.


Author(s):  
Chunill Hah ◽  
Douglas C. Rabe ◽  
Thomas J. Sullivan ◽  
Aspi R. Wadia

The effects of circumferential distortions in inlet total pressure on the flow field in a low-aspect-ratio, high-speed, high-pressure-ratio, transonic compressor rotor are investigated in this paper. The flow field was studied experimentally and numerically with and without inlet total pressure distortion. Total pressure distortion was created by screens mounted upstream from the rotor inlet. Circumferential distortions of 8 periods per revolution were investigated at two different rotor speeds. The unsteady blade surface pressures were measured with miniature pressure transducers mounted in the blade. The flow fields with and without inlet total pressure distortion were analyzed numerically by solving steady and unsteady forms of the Reynolds-averaged Navier-Stokes equations. Steady three-dimensional viscous flow calculations were performed for the flow without inlet distortion while unsteady three-dimensional viscous flow calculations were used for the flow with inlet distortion. For the time-accurate calculation, circumferential and radial variations of the inlet total pressure were used as a time-dependent inflow boundary condition. A second-order implicit scheme was used for the time integration. The experimental measurements and the numerical analysis are highly complementary for this study because of the extreme complexity of the flow field. The current investigation shows that inlet flow distortions travel through the rotor blade passage and are convected into the following stator. At a high rotor speed where the flow is transonic, the passage shock was found to oscillate by as much as 20% of the blade chord, and very strong interactions between the unsteady passage shock and the blade boundary layer were observed. This interaction increases the effective blockage of the passage, resulting in an increased aerodynamic loss and a reduced stall margin. The strong interaction between the passage shock and the blade boundary layer increases the peak aerodynamic loss by about one percent.


Author(s):  
N. Lymberopoulos ◽  
K. Giannakoglou ◽  
I. Nikolaou ◽  
K. D. Papailiou ◽  
A. Tourlidakis ◽  
...  

Mechanical constraints dictate the existence of tip clearances in rotating cascades, resulting to a flow leakage through this clearance which considerably influences the efficiency and range of operation of the machine. Three-dimensional Navier-Stokes solvers are often used for the numerical study of compressor and turbine stages with tip-clearance. The quality of numerical predictions depends strongly on how accurately the blade tip region is modelled; in this respect the accurate modelling of tip region was one of the main goals of this work. In the present paper, a 3-D Navier-Stokes solver is suitably adapted so that the flat tip surface of a blade and its sharp edges could be accurately modelled, in order to improve the precision of the calculation in the tip region. The adapted code solves the fully elliptic, steady, Navier-Stokes equations through a space-marching algorithm and a pressure correction technique; the H-type topology is retained, even in cases with thick leading edges where a special treatment is introduced herein. The analysis is applied to two different cases, a linear cascade and a compressor rotor, and comparisons with experimental data are provided.


Author(s):  
W. R. Briley ◽  
D. V. Roscoe ◽  
H. J. Gibeling ◽  
R. C. Buggeln ◽  
J. S. Sabnis ◽  
...  

Three-dimensional solutions of the ensemble-averaged Navier-Stokes equations have been computed for a high-turning turbine rotor passage, both with and without tip clearance effects. The geometry is Pratt & Whitney’s preliminary design for the Generic Gas Generator Turbine (GGGT), having an axial chord of 0.5 inch and turning angle of about 160 degrees. The solutions match the design Reynolds number of 3x 106/inch and design inflow/outflow distributions of flow quantities. The grid contains 627,000 points, including 20 radial points in the clearance gap of 0.015 inch, and has a minimum spacing of 10−4 inch adjacent to all surfaces. The solutions account for relative motion of the blade and shroud surfaces and include a backstep on the shroud. Computed results are presented which show the general flow behavior, especially near the tip clearance and backstep regions. The results are generally consistent with experimental observations for other geometries having thinner blades and smaller turning angles. The leakage flow includes some fluid originally in the freestream at 91 percent span. Downstream, the leakage flow behaves as a wall jet directed at 100 degrees to the main stream, with total pressure and temperature higher than the freestream. Radial distributions of circumferentially-averaged flow quantities are compared for solutions with and without tip leakage flow. Two-dimensional solutions are also presented for the mid-span blade geometry for design and off-design inflow angles.


1969 ◽  
Vol 37 (4) ◽  
pp. 727-750 ◽  
Author(s):  
Gareth P. Williams

A method of numerically integrating the Navier-Stokes equations for certain three-dimensional incompressible flows is described. The technique is presented through application to the particular problem of describing thermal convection in a rotating annulus. The equations, in cylindrical polar co-ordinate form, are integrated with respect to time by a marching process, together with the solving of a Poisson equation for the pressure. A suitable form of the finite difference equations gives a computationally-stable long-term integration with reasonably faithful representation of the spatial and temporal characteristics of the flow.Trigonometric interpolation techniques provide accurate (discretely exact) solutions to the Poisson equation. By using an auxiliary algorithm for rapid evaluation of trigonometric transforms, the proportion of computation needed to solve the Poisson equation can be reduced to less than 25% of the total time needed to’ advance one time step. Computing on a UNIVAC 1108 machine, the flow can be advanced one time-step in 2 sec for a 14 × 14 × 14 grid upward to 96 sec for a 60 × 34 × 34 grid.As an example of the method, some features of a solution for steady wave flow in annulus convection are presented. The resemblance of this flow to the classical Eady wave is noted.


Fluids ◽  
2019 ◽  
Vol 4 (2) ◽  
pp. 94 ◽  
Author(s):  
Cornel Marius Murea

A monolithic semi-implicit method is presented for three-dimensional simulation of fluid–structure interaction problems. The updated Lagrangian framework is used for the structure modeled by linear elasticity equation and, for the fluid governed by the Navier–Stokes equations, we employ the Arbitrary Lagrangian Eulerian method. We use a global mesh for the fluid–structure domain where the fluid–structure interface is an interior boundary. The continuity of velocity at the interface is automatically satisfied by using globally continuous finite element for the velocity in the fluid–structure mesh. The method is fast because we solve only a linear system at each time step. Three-dimensional numerical tests are presented.


1991 ◽  
Vol 113 (2) ◽  
pp. 241-250 ◽  
Author(s):  
C. Hah ◽  
A. J. Wennerstrom

The concept of swept blades for a transonic or supersonic compressor was reconsidered by Wennerstrom in the early 1980s. Several transonic rotors designed with swept blades have shown very good aerodynamic efficiency. The improved performance of the rotor is believed to be due to reduced shock strength near the shroud and better distribution of secondary flows. A three-dimensional flowfield inside a transonic rotor with swept blades is analyzed in detail experimentally and numerically. A Reynolds-averaged Navier–Stokes equation is solved for the flow inside the rotor. The numerical solution is based on a high-order upwinding relaxation scheme, and a two-equation turbulence model with a low Reynolds number modification is used for the turbulence modeling. To predict flows near the shroud properly, the tip-clearance flow also must be properly calculated. The numerical results at three different operating conditions agree well with the available experimental data and reveal various interesting aspects of shock structure inside the rotor.


1990 ◽  
Vol 112 (1) ◽  
pp. 109-115 ◽  
Author(s):  
N. M. McDougall

Detailed measurements have been made within an axial compressor operating both at design point and near stall. Rotor tip clearance was found to control the performance of the machine by influencing the flow within the rotor blade passages. This was not found to be the case in the stator blade row, where hub clearance was introduced beneath the blade tips. Although the passage flow was observed to be altered dramatically, no significant changes were apparent in the overall pressure rise or stall point. Small tip clearances in the rotor blade row resulted in the formation of corner separations at the hub, where the blade loading was highest. More representative clearances resulted in blockage at the tip due to the increased tip clearance flow. The effects that have been observed emphasize both the three-dimensional nature of the flow within compressor blade passages, and the importance of the flow in the endwall regions in determining the overall compressor performance.


Sign in / Sign up

Export Citation Format

Share Document