Performance of choked unsteady ejector-nozzles for use in pressure-gain combustors

If the conventional steady flow combustor of a gas turbine is replaced with a device which achieves a pressure gain during the combustion process then the thermal efficiency of the cycle is raised. All such 'Pressure Gain Combustors' (e.g. PDEs, pulse combustors or wave rotors) are inherently unsteady flow devices. For such a device to be practically installed in a gas turbine it is necessary to design a downstream row of turbine vanes which will both accept the combustors unsteady exit flow and deliver a flow which the turbine rotor can accept. The design requirements of such a vane are that its exit flow both retains the maximum time-mean stagnation pressure gain (the pressure gain produced by the combustor is not lost) and minimises the amplitude of unsteadiness (reduces unsteadiness entering the downstream rotor). In this paper the exit of the pressure gain combustor is simulated with a cold unsteady jet. The first stage vane is simulated by a one-dimensional choked ejector nozzle with no turning. The time-mean and rms stagnation pressure at nozzle exit is measured. A number of geometric configurations are investigated and it is shown that the optimal geometry both maximizes time mean stagnation pressure gain (75% of that in the exit of the unsteady jet) and minimizes the amplitude of unsteadiness (1/3 of that in the primary jet). The structure of the unsteady flow within the ejector nozzle is determined computationally. Copyright © 2009 by J Heffer and R Miller.

Airfoil flow and noise computation using monotonically integrated large eddy simulation and acoustic analogy: Effect of the grid resolution.

A new scalable Monotonically Integrated Large Eddy Simulation (MILES) method based on the Compact Accurately Boundary-Adjusting high-REsolution Technique (CABARET) has been applied for the simulation of unsteady flow around NACA0012 airfoil at Re = 400,000 and M = 0.058. The flow solution is coupled with the Ffowcs Williams-Hawkings formulation for far-field noise prediction. The computational modeling results are presented for several computational grid resolutions: 8, 16, and 32 million grid cells and compared with the experimental data available.

Grid-refined LES predictions for Fan-OGV broadband noise

This work forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. In this paper, we focus on LES calculations of noise sources on and close to a blade trailing edge. We consider two test cases; one an isolated NACA0012 airfoil in flow, and the other an industry-standard rotating fan. In the first case, turbulent mean and RMS velocities and energy spectra at different locations are compared with those from experiment. 1,2The sound generated by the unsteady pressure fluctuations on the airfoil surface and by the flow turbulence will be predicted using a Ffowcs Williams Hawkings (FW-H) surface. In the second case, unsteady flow and acoustic fields around the blade passage 3 are presented for a refined mesh, and the rotor-stator tonal noise will be predicted by using the rotor-wake mean velocity profile and the methodology described in Lloyd & Peake 4. Copyright © 2009 by Qinling Li, Nigel Peake & Mark Savill.

Turbomachinery computational fluid dynamics: asymptotes and paradigm shifts.

This paper reviews the development of computational fluid dynamics (CFD) specifically for turbomachinery simulations and with a particular focus on application to problems with complex geometry. The review is structured by considering this development as a series of paradigm shifts, followed by asymptotes. The original S1-S2 blade-blade-throughflow model is briefly described, followed by the development of two-dimensional then three-dimensional blade-blade analysis. This in turn evolved from inviscid to viscous analysis and then from steady to unsteady flow simulations. This development trajectory led over a surprisingly small number of years to an accepted approach-a 'CFD orthodoxy'. A very important current area of intense interest and activity in turbomachinery simulation is in accounting for real geometry effects, not just in the secondary air and turbine cooling systems but also associated with the primary path. The requirements here are threefold: capturing and representing these geometries in a computer model; making rapid design changes to these complex geometries; and managing the very large associated computational models on PC clusters. Accordingly, the challenges in the application of the current CFD orthodoxy to complex geometries are described in some detail. The main aim of this paper is to argue that the current CFD orthodoxy is on a new asymptote and is not in fact suited for application to complex geometries and that a paradigm shift must be sought. In particular, the new paradigm must be geometry centric and inherently parallel without serial bottlenecks. The main contribution of this paper is to describe such a potential paradigm shift, inspired by the animation industry, based on a fundamental shift in perspective from explicit to implicit geometry and then illustrate this with a number of applications to turbomachinery.

Effects of rotating inlet distortion on multistage compressor stability

In multi-spool engines, rotating stall in an upstream compressor will impose a rotating distortion on the downstream compressor, thereby affecting its stability margin. In this paper experiments are described in which this effect was simulated by a rotating screen upstream of several multistage low-speed compressors. The measurements are complemented by, and compared with, a theoretical model of multistage compressor response to speed and direction of rotation of an inlet distortion. For co-rotating distortions (i.e., distortions rotating in the same direction as rotor rotation), experiments show that the compressors exhibited significant loss in stability margin and that they could be divided into two groups according to their response. The first group exhibited a single peak in stall margin degradation when the distortion speed corresponded to roughly 50% of rotor speed. The second group showed two peaks in stall margin degradation corresponding to distortion speeds of approximately 25-35% and 70-75% of rotor speed. These new results demonstrate that multistage compressors can have more than a single resonant response. Detailed measurements suggest that the two types of behavior are linked to differences between the stall inception processes observed for the two groups of compressors and that a direct connection thus exists between the observed forced response and the unsteady flow phenomena at stall onset. For counter-rotational distortions, all the compressors tested showed minimal loss of stability margin. The results imply that counter-rotation of the fan and core compressor, or LP and HP compressors, could be a worthwhile design choice. Calculations based on the two-dimensional theoretical model show excellent agreement for the compressors which had a single peak for stall margin degradation. We take this first-of-a-kind comparison as showing that the model, though simplified, captures the essential fluid dynamic features of the phenomena. Agreement is not good for compressors which had two peaks in the curve of stall margin shift versus distortion rotation speed. The discrepancy is attributed to the three-dimensional and short length scale nature of the stall inception process in these machines; this includes phenomena that have not yet been addressed in any model.

The development of aerodynamic performance measurements in a transient test facility

Transient test facilities offer the potential for the simultaneous study of turbine aerodynamic performance, unsteady flow phenomena and the heat transfer characteristics of a turbine stage. This paper describes the development of aerodynamic performance measurement techniques in the Oxford Rotor Facility (ORF). The solutions to the technological issues involved with transient testing presented in this paper are expected to achieve levels of precision uncertainty comparable with traditional steady flow test rigs. The theoretical background to the measurement of aerodynamic performance is presented together with a comprehensive pre-test uncertainty analysis. The instrumentation scheme for the measurement of stage mass flow rate is discussed in detail, the measurements of shaft power, total inlet enthalpy, and stage pressure ratio are also outlined. The current working section features a 62% scale, 1-1/2 stage, high-pressure shroudless transonic turbine. The required inlet flow conditions are provided by an Isentropic Light Piston Tunnel (ILPT) with a quasi-steady state run time of approximately 70ms. The testing is conducted at engine representative specific speed, pressure ratio, gas-to-wall temperature ratio, Mach number and Reynolds number.

Vortex Transport and Blade Interactions in High Pressure Turbines

This paper presents a study of the three-dimensional flow field within the blade rows of a high-pressure axial flow steam turbine stage. Half-delta wings were fixed to a rotating hub to simulate an upstream rotor passage vortex. The flow field is investigated in a Low-Speed Research Turbine using pneumatic and hot-wire probes downstream of the blade row. The paper examines the impact of the delta wing vortex transport on the performance of the downstream blade row. Steady and unsteady numerical simulations were performed using structured 3D Navier-Stokes solver to further understand the flow field. The loss measurements at the exit of the stator blade showed an increase in stagnation pressure loss due to the delta wing vortex transport. The increase in loss was 21% of the datum stator loss, demonstrating the importance of this vortex interaction. The transport of the stator viscous flow through the rotor blade row is also described. The rotor exit flow was affected by the interaction between the enhanced stator passage vortex and the rotor blade row. Flow underturning near the hub and overturning towards the mid-span was observed, contrary to the classical model of overturning near the hub and underturning towards the mid-span. The unsteady numerical simulation results were further analysed to identify the entropy producing regions in the unsteady flow field.

The effects of a trip wire and unsteadiness on a high speed highly loaded low-pressure turbine blade

This paper presents the effect of a single spanwise 2D wire upon the downstream position of boundary layer transition under steady and unsteady inflow conditions. The study is carried out on a high turning, high-speed, low pressure turbine (LPT) profile designed to take account of the unsteady flow conditions. The experiments were carried out in a transonic cascade wind tunnel to which a rotating bar system had been added. The range of Reynolds and Mach numbers studied includes realistic LPT engine conditions and extends up to the transonic regime. Losses are measured to quantify the influence of the roughness with and without wake passing. Time resolved measurements such as hot wire boundary layer surveys and surface unsteady pressure are used to explain the state of the boundary layer. The results suggest that the effect of roughness on boundary layer transition is a stability governed phenomena, even at high Mach numbers. The combination of the effect of the roughness elements with the inviscid Kelvin-Helmholtz instability responsible for the rolling up of the separated shear layer (Stieger [1]) is also examined. Wake traverses using pneumatic probes downstream of the cascade reveal that the use of roughness elements reduces the profile losses up to exit Mach numbers of 0.8. This occurs with both steady and unsteady inflow conditions.

An experimental and numerical study of an oscillating transonic shock wave in a duct

A combined experimental and numerical study of a transonic shock wave in a parallel walled duct subject to downstream pressure perturbations has been conducted. Experiments and simulations have been carried out with a shock strength of M∞ = 1.4 for pressure perturbation frequencies in the range 16-90 Hz. The dynamics of unsteady shock motion and the interaction structure between the unsteady transonic shock wave and the turbulent tunnel floor boundary layer have been investigated. It is found that the (experimentally measured) dynamics of shock motion are generally well predicted by the computational scheme, especially at relatively low (≈ 40 Hz) frequencies. However, at higher frequencies (≈ 90 Hz), some subtle differences between the shock dynamics measured in experiments and those predicted by Computational Fluid Dynamics (CFD) exist. There is evidence from experiments that variations in shock / boundary layer interaction (SBLI) structure caused by shock motion are responsible for a change in the nature of shock dynamics between low and high frequency. In contrast, numerical results at low and high frequencies do not differ significantly and this suggests that the numerical method is not fully capturing the physics of the unsteady flow. Possible reasons for this are considered and a number of areas where CFD is unable to replicate experimental observations are identified. Significantly, CFD predicts changes in SBLI structure due to shock motion that are much too large and this may explain why none of the subtle effects on shock dynamics seen in experiments occur in CFD. Further work developing numerical methods that demonstrate a more realistic sensitivity of SBLI structure to unsteady shock motion is required. Copyright © 2010 by P.J.K. Bruce.

The nonlinear heat release response of stratified lean-premixed flames to acoustic velocity oscillations

This paper analyzes the forced response of swirl-stabilized lean-premixed flames to high-amplitude acoustic forcing in a laboratory-scale stratified burner operated with CH4 and air at atmospheric pressure. The double-swirler, double-channel annular burner was specially designed to generate high-amplitude acoustic velocity oscillations and a radial equivalence ratio gradient at the inlet of the combustion chamber. Temporal oscillations of equivalence ratio along the axial direction are dissipated over a long distance, and therefore the effects of time-varying fuel/air ratio on the response are not considered in the present investigation. Simultaneous measurements of inlet velocity and heat release rate oscillations were made using a constant temperature anemometer and photomultiplier tubes with narrow-band OH*/CH* interference filters. Time-averaged and phase-synchronized CH* chemiluminescence intensities were measured using an intensified CCD camera. The measurements show that flame stabilization mechanisms vary depending on equivalence ratio gradients for a constant global equivalence ratio (φg=0.60). Under uniformly premixed conditions, an enveloped M-shaped flame is observed. In contrast, under stratified conditions, a dihedral V-flame and a toroidal detached flame develop in the outer stream and inner stream fuel enrichment cases, respectively. The modification of the stabilization mechanism has a significant impact on the nonlinear response of stratified flames to high-amplitude acoustic forcing (u'/U∼0.45 and f=60, 160Hz). Outer stream enrichment tends to improve the flame's stiffness with respect to incident acoustic/vortical disturbances, whereas inner stream stratification tends to enhance the nonlinear flame dynamics, as manifested by the complex interaction between the swirl flame and large-scale coherent vortices with different length scales and shedding points. It was found that the behavior of the measured flame describing functions (FDF), which depend on radial fuel stratification, are well correlated with previous measurements of the intensity of self-excited combustion instabilities in the stratified swirl burner. The results presented in this paper provide insight into the impact of nonuniform reactant stoichiometry on combustion instabilities, its effect on flame location and the interaction with unsteady flow structures. © 2011 The Combustion Institute.

Modern challenges facing turbomachinery aeroacoustics

In this review we describe current scientific and technological issues in the quest to reduce aeroengine noise, in the face of predicted rapid increases in the volume of air traffic, on the one hand, and increasingly strict environmental regulation, on the other. Alongside conventional ducted turbofan designs, new open-rotor contra-rotating power plants are currently under development, which present their own noise challenges. The key sources of tonal and broadband noise, and the way in which noise propagates away from the source, are surveyed in both cases. We also consider in detail two key aspects underpinning the flow physics that continue to receive considerable attention, namely the acoustics of swirling flow and unsteady flow-blade interactions. Finally, we describe possible innovations in open-rotor engine design for low noise.

Lift and the leading-edge vortex

Flapping wings often feature a leading-edge vortex (LEV) that is thought to enhance the lift generated by the wing. Here the lift on a wing featuring a leading-edge vortex is considered by performing experiments on a translating flat-plate aerofoil that is accelerated from rest in a water towing tank at a fixed angle of attack of 15°. The unsteady flow is investigated with dye flow visualization, particle image velocimetry (PIV) and force measurements. Leading-and trailing-edge vortex circulation and position are calculated directly from the velocity vectors obtained using PIV. In order to determine the most appropriate value of bound circulation, a two-dimensional potential flow model is employed and flow fields are calculated for a range of values of bound circulation. In this way, the value of bound circulation is selected to give the best fit between the experimental velocity field and the potential flow field. Early in the trajectory, the value of bound circulation calculated using this potential flow method is in accordance with Kelvin's circulation theorem, but differs from the values predicted by Wagner's growth of bound circulation and the Kutta condition. Later the Kutta condition is established but the bound circulation remains small; most of the circulation is contained instead in the LEVs. The growth of wake circulation can be approximated by Wagner's circulation curve. Superimposing the non-circulatory lift, approximated from the potential flow model, and Wagner's lift curve gives a first-order approximation of the measured lift. Lift is generated by inertial effects and the slow buildup of circulation, which is contained in shed vortices rather than bound circulation. © 2013 Cambridge University Press.

On the time delay in continuous in-cylinder sampling from IC engines

When gas sample is continuously drawn from the cylinder of an internal combustion engine, the sample that appears at the end of the sampling system corresponds to the in-cylinder content sometime ago because of the finite transit time which is a function of the cylinder pressure history. This variable delay causes a dispersion of the sample signal and makes the interpretation of the signal difficult An unsteady flow analysis of a typical sampling system was carried out for selected engine loads and speeds. For typical engine operation, a window in which the delay is approximately constant may be found. This window gets smaller with increase in engine speed, with decrease in load, and with the increase in exit pressure of the sampling system.