A Navier Stokes Analysis of Airfoils in Oscillating Transonic Cascades for the Prediction of Aerodynamic Damping

Author(s):  
Reza S. Abharl ◽  
Michael Giles

An unsteady, compressible, two dimensional, thin shear layer Navier Stokes solver is modified to predict the motion-dependent unsteady flow around oscillating airfoils in a cascade. A quasi 3-D formulation is used to account for the streamwise variation of streamtube height. The code uses Ni’s Lax-Wendroff algorithm in the outer region, an implicit ADI method in the inner region, conservative coupling at the interface, and the Baldwin-Lomax turbulence model. Computational mesh consists of an O-grid around each blade plus an unstructured outer grid of quadrilateral or triangular cells. The unstructured computational grid was adapted to the flow to better resolve shocks and wakes. Motion of each airfoil was simulated at each time step by stretching and compressing the mesh within the O-grid. This imposed motion consists of harmonic solid body translation in two directions and rotation, combined with the correct inter-blade phase angles. Validity of the code is illustrated by comparing its predictions to a number of test cases, including an axially oscillating flat plate in laminar flow, the Aeroelasticity of Turbomachines Symposium Fourth Standard Configuration (a transonic turbine cascade), the Seventh Standard Configuration (a transonic compressor cascade). The overall comparison between the predictions and the test data is reasonably good. A numerical study on a generic transonic compressor rotor was performed in which the impact of varying the amplitude of the airfoil oscillation on the normalized predicted magnitude and phase of the unsteady pressure around the airfoil was studied. It was observed that for this transonic compressor, the non-dimensional aerodynamic damping was influenced by the amplitude of the oscillation.

1997 ◽  
Vol 119 (1) ◽  
pp. 77-84 ◽  
Author(s):  
R. S. Abhari ◽  
M. Giles

An unsteady, compressible, two-dimensional, thin shear layer Navier–Stokes solver is modified to predict the motion-dependent unsteady flow around oscillating airfoils in a cascade. A quasi-three-dimensional formulations is used to account for the stream-wise variation of streamtube height. The code uses Ni’s Lax–Wendroff algorithm in the outer region, an implicit ADI method in the inner region, conservative coupling at the interface, and the Baldwin–Lomax turbulence model. The computational mesh consists of an O-grid around each blade plus an unstructured outer grid of quadrilateral or triangular cells. The unstructured computational grid was adapted to the flow to better resolve shocks and wakes. Motion of each airfoil was simulated at each time step by stretching and compressing the mesh within the O-grid. This imposed motion consists of harmonic solid body translation in two directions and rotation, combined with the correct interblade phase angles. The validity of the code is illustrated by comparing its predictions to a number of test cases, including an axially oscillating flat plate in laminar flow, the Aeroelasticity of Turbomachines Symposium Fourth Standard Configuration (a transonic turbine cascade), and the Seventh Standard Configuration (a transonic compressor cascade). The overall comparison between the predictions and the test data is reasonably good. A numerical study on a generic transonic compressor rotor was performed in which the impact of varying the amplitude of the airfoil oscillation on the normalized predicted magnitude and phase of the unsteady pressure around the airfoil was studied. It was observed that for this transonic compressor, the nondimensional aerodynamic damping was influenced by the amplitude of the oscillation.


2000 ◽  
Vol 123 (2) ◽  
pp. 409-417 ◽  
Author(s):  
Bjo¨rn Gru¨ber ◽  
Volker Carstens

A parametric study which investigates the influence of viscous effects on the damping behavior of vibrating compressor cascades is presented here. To demonstrate the dependence of unsteady aerodynamic forces on the flow viscosity, a computational study was performed for a transonic compressor cascade of which the blades underwent tuned pitching oscillations while the flow conditions extended from fully subsonic to highly transonic flow. Additionally, the reduced frequency and Reynolds number were varied. In order to check the linear behavior of the aerodynamic forces, all calculations were carried out for three different oscillation amplitudes. Comparisons with inviscid Euler results helped identify the influence of viscous effects. The computations were performed with a Navier-Stokes code, the basic features of which are the use of an AUSM upwind scheme, an implicit time integration, and the implementation of the Baldwin-Lomax turbulence model. In order to demonstrate the possibility of this code to correctly predict the unsteady behavior of strong shock-boundary layer interactions, the experiment of Yamamoto and Tanida on a self-induced shock oscillation due to shock-boundary layer interaction was calculated. A significant improvement in the prediction of the shock amplitude was achieved by a slight modification of the Baldwin Lomax turbulence model. An important result of the presented compressor cascade investigations is that viscous effects may cause a significant change in the aerodynamic damping. This behavior is demonstrated by two cases in which an Euler calculation predicts a damped oscillation whereas a Navier-Stokes computation leads to an excited vibration. It was found that the reason for these contrary results are shock-boundary-layer interactions which dramatically change the aerodynamic damping.


Author(s):  
Björn Grüber ◽  
Volker Carstens

A parametric study which investigates the influence of viscous effects on the damping behaviour of vibrating compressor cascades is presented here. To demonstrate the dependence of unsteady aerodynamic forces on the flow viscosity, a computational study was performed for a transonic compressor cascade of which the blades underwent tuned pitching oscillations while the flow conditions extended from fully subsonic to highly transonic flow. Additionally, the reduced frequency and Reynolds number were varied. In order to check the linear behavior of the aerodynamic forces, all calculations were carried out for three different oscillation amplitudes. Comparisons with inviscid Euler results helped identify the influence of viscous effects. The computations were performed with a Navier-Stokes code, the basic features of which are the use of an AUSM upwind scheme, an implicit time integration, and the implementation of the Baldwin-Lomax turbulence model. In order to demonstrate the possibility of this code to correctly predict the unsteady behavior of strong shock-boundary layer interactions, the experiment of Yamamoto and Tanida on a self-induced shock oscillation due to shock-boundary layer interaction was calculated. A significant improvement in the prediction of the shock amplitude was achieved by a slight modification of the Baldwin Lomax turbulence model. An important result of the presented compressor cascade investigations is that viscous effects may cause a significant change in the aerodynamic damping. This behaviour is demonstrated by two cases in which an Euler calculation predicts a damped oscillation whereas a Navier-Stokes computation leads to an excited vibration. It was found that the reason for these contrary results are shock-boundary-layer interactions which dramatically change the aerodynamic damping.


Author(s):  
Paul J. Petrie-Repar ◽  
Andrew McGhee ◽  
Peter A. Jacobs ◽  
Rowan Gollan

In this paper, analytical maps of aerodynamic damping for a two-dimensional compressor cascade (Standard Configuration 10) are presented. The maps are shown as contour plots of the aerodynamic damping as a function of operating condition. The aerodynamic dampings were calculated by a linearized Navier-Stokes flow solver. The flutter boundaries over a wide range of operating conditions are clearly shown on the damping maps and were found to be strongly dependent on the mode frequency and the mode shape. Extremely low values of negative aerodynamic damping were predicted for some off-design operating conditions where flow separation occurred. A damping map was also constructed based on inviscid flow simulations. There were differences in the viscous and inviscid flutter boundaries particularly at off-design inflow angles. The extremely low values of negative aerodynamic damping were only predicted by the viscous simulations and not the inviscid simulations.


Author(s):  
W. N. Dawes

This paper presents a numerical simulation of the unsteady flow in a transonic compressor stage. The simulations were performed using an unstructured mesh, 3D Navier-Stokes solver with a sliding interface between the mesh blocks containing the rotor and stator blades to permit a time-resolved calculation. The focus in the simulation is the endwall flow, its contrast with the mid-span flow and, in particular, the interaction of the rotor overtip leakage flow with the downstream stator. The endwall region of both the rotors and the stators was predicted to be proportionally responsible for much more loss production than the mid-span sections of the blading. Unsteady effects were predicted to be responsible for rather more of the total loss production in the rotors than in the stators but the unsteady loss production in the rotors did not seem to be associated with the endwall flow but rather occurred over the bulk of the blade span caused by unsteady shock motion. By contrast, the rotor overtip leakage flow was shown to cause a considerable degradation of the stator performance near the casing. In the stators nearly a third of the loss production near the endwall was predicted to be directly due to unsteady effects associated with the rotor vortex interaction and a number of individual physical loss generating mechanisms were identified.


1997 ◽  
Vol 119 (1) ◽  
pp. 122-128 ◽  
Author(s):  
S. L. Puterbaugh ◽  
W. W. Copenhaver

An experimental investigation concerning tip flow field unsteadiness was performed for a high-performance, state-of-the-art transonic compressor rotor. Casing-mounted high frequency response pressure transducers were used to indicate both the ensemble averaged and time varying flow structure present in the tip region of the rotor at four different operating points at design speed. The ensemble averaged information revealed the shock structure as it evolved from a dual shock system at open throttle to an attached shock at peak efficiency to a detached orientation at near stall. Steady three-dimensional Navier Stokes analysis reveals the dominant flow structures in the tip region in support of the ensemble averaged measurements. A tip leakage vortex is evident at all operating points as regions of low static pressure and appears in the same location as the vortex found in the numerical solution. An unsteadiness parameter was calculated to quantify the unsteadiness in the tip cascade plane. In general, regions of peak unsteadiness appear near shocks and in the area interpreted as the shock-tip leakage vortex interaction. Local peaks of unsteadiness appear in mid-passage downstream of the shock-vortex interaction. Flow field features not evident in the ensemble averaged data are examined via a Navier-Stokes solution obtained at the near stall operating point.


1993 ◽  
Vol 115 (2) ◽  
pp. 283-295 ◽  
Author(s):  
W. N. Dawes

This paper describes recent developments to a three-dimensional, unstructured mesh, solution-adaptive Navier–Stokes solver. By adopting a simple, pragmatic but systematic approach to mesh generation, the range of simulations that can be attempted is extended toward arbitrary geometries. The combined benefits of the approach result in a powerful analytical ability. Solutions for a wide range of flows are presented, including a transonic compressor rotor, a centrifugal impeller, a steam turbine nozzle guide vane with casing extraction belt, the internal coolant passage of a radial inflow turbine, and a turbine disk cavity flow.


Author(s):  
Minghao Pan ◽  
Paul Petrie-Repar ◽  
Hans Mårtensson ◽  
Tianrui Sun ◽  
Tobias Gezork

In turbomachines, forced response of blades is blade vibrations due to external aerodynamic excitations and it can lead to blade failures which can have fatal or severe economic consequences. The estimation of the level of vibration due to forced response is dependent on the determination of aerodynamic damping. The most critical cases for forced response occur at high reduced frequencies. This paper investigates the determination of aerodynamic damping at high reduced frequencies. The aerodynamic damping was calculated by a linearized Navier-Stokes flow solver with exact 3D non-reflecting boundary conditions. The method was validated using Standard Configuration 8, a two-dimensional flat plate. Good agreement with the reference data at reduced frequency 2.0 was achieved and grid converged solutions with reduced frequency up to 16.0 were obtained. It was concluded that at least 20 cells per wavelength is required. A 3D profile was also investigated: an aeroelastic turbine rig (AETR) which is a subsonic turbine case. In the AETR case, the first bending mode with reduced frequency 2.0 was studied. The 3D acoustic modes were calculated at the far-fields and the propagating amplitude was plotted as a function of circumferential mode index and radial order. This plot identified six acoustic resonance points which included two points corresponding to the first radial modes. The aerodynamic damping as a function of nodal diameter was also calculated and plotted. There were six distinct peaks which occurred in the damping curve and these peaks correspond to the six resonance points. This demonstrates for the first time that acoustic resonances due to higher order radial acoustic modes can affect the aerodynamic damping at high reduced frequencies.


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.


Sign in / Sign up

Export Citation Format

Share Document