Wave propagation in a slowly-varying duct with mean vortical swirling flow

The propagation of unsteady disturbances in a slowlyvarying cylindrical duct carrying mean swirling flow is investigated using a multiple-scales technique. This is applicable to turbomachinery flow behind a rotor stage when the swirl and axial velocities are of the same order. The presence of mean vorticity couples acoustic and vorticity equations which produces an eigenvalue problem that is not self-adjoint unlike that for irrotational mean flow. In order to determine the amplitude variation along the duct, an adjoint solution for the coupled system of equations is derived. The solution breaks down where a mode changes from cut on to cut off. In this region the amplitude is governed by a form of Airy's equation, and the effect of swirl is to introduce a small shift in the origin of the Airy function away from the turning-point location. The variation of axial wavenumber and amplitude along the duct is calculated. In hard-walled ducts mean swirl is shown to produce much larger amplitude variation along the duct compared with a nonswirling flow. Mean swirl also has a large effect in ducts with finite-impedance walls which differs depending on whether modes are co-rotating with the swirl or counter rotating. © 2001 by A.J. Cooper, Published by the American Institute of Aeronautics and Astronautics, Inc.

Models for acoustic propagation through turbofan exhaust flows-lined ducts

Turbomachinery noise radiating into the rearward arc is an important problem. This noise is scattered by the trailing edges of the nacelle and the jet exhaust, and interacts with the shear layers between the external flow, bypass stream and jet, en route to the far field. In the past a range of relevant model problems involving semi-infinite cylinders have been solved. However, one limitation of these previous solutions is that they do not allow for the jet nozzle protruding a finite distance beyond the end of the nacelle (or in certain configurations being buried a finite distance upstream). With this in mind, we have used the matrix Wiener-Hopf technique to allow precisely this finite nacelle-jet nozzle separation to be included. We have previously reported results for the case of hard-walled ducts, which requires factorisation of a 2 × 2 matrix. In this paper we extend this work by allowing one of the duct walls, in this case the outer wall of the jet pipe, to be acoustically lined. This results in the need to factorise a 3 × 3 matrix, which is completed by use of a combination of pole-removal and Pad́e approximant techniques. Sample results are presented, investigating in particular the effects of exit plane stagger and liner impedance. Here we take the mean flow to be zero, but extension to nonzero Mach numbers in the core and bypass flow has also been completed. Copyright © 2009 by Nigel Peake & Ben Veitch.

Sound transmission in strongly-curved slowly-varying cylindrical and annular lined ducts with flow

In this paper we consider the propagation of acoustic waves along a curved hollow or annular duct with lined walls. The curvature of the duct centreline and the wall radii vary slowly along the duct, allowing application of an asymptotic multiple scales analysis. This generalises Rienstra's analysis of a straight duct of varying cross-sectional radius. The result of the analysis is that the modal wavenumbers and mode shapes are determined locally as modes of a torus with the same local curvature, while the amplitude of the modes evolves as the mode propagates along the duct. The duct modes are found numerically at each axial location using a pseudo-spectral method. Unlike the case of a straight duct, there is a fundamental asymmetry between upstream and downstream propagating modes, with some mode shapes tending to be concentrated on either the inside or outside of the bend depending on the direction of propagation. The interaction between the presence of wall lining and curvature is investigated in particular; for instance, in a representative case it is found that the curvature causes the first few acoustic modes to be more heavily damped by the duct boundary than would be expected for a straight duct. Analytical progress can be made in the limit of very high mode order, in which case well-known 'whispering gallery' modes, localised close to the wall, can be identified.

Generation and scattering of sound in vortical flow

We present solutions to scattering problems for unsteady disturbances to a mean swirling flow in an annular duct with a rigid 'splitter'. This situation has application to rotor-stator interaction noise in aeroengines, where the flow downstream of the fan is swirling and bifurcates into the by-pass duct and the engine core. We also consider the trailing edge extension of this problem. Inviscid mean flow in a cylindrical annulus is considered, with both axial and swirling (azimuthal) velocity components. The presence of vorticity in the mean flow couples the acoustic and vorticity modes of irrotational flow. Instead we have one combined spectrum of acoustic-vorticity waves in which the 'sonic' and 'nearly-convected' modes are fully coupled. In addition to the aeroacoustics application the results offer insight into the behaviour of these acoustic-vorticity waves, and the precise nature of the coupling between the two types of mode. Two regimes are discussed in which progress has been made, one for a specialised mean flow, uniform axial flow and rigid body swirl, and a second regime in which the frequency is assumed large, valid for any axisymmetric mean flow. The Wiener-Hopf technique is used to solve the scattering problems mathematically, and we present numerical evaluations of these solutions. Several new effects are seen to arise due to the mean vorticity, in particular the generation of sound at a trailing edge due to the scattering of a nearly convected disturbance, in contrast to the way a convected gust silently passes a trailing edge in uniform mean flow.

Wake, shock and potential field interactions in a 1.5 stage turbine: Part I: Vane-rotor and rotor-vane interaction

The composition of the time-resolved surface pressure field around a high-pressure rotor blade caused by the presence of neighboring blade rows was studied, with the individual effects of wake, shock and potential field interaction being determined. Two test geometries were considered: first, a high-pressure turbine stage coupled with a swan-necked diffuser exit duct; secondly, the same high-pressure stage but with a vane located in the downstream duct. Both tests were carried out at engine-representative Mach and Reynolds numbers. By comparing the results to time-resolved computational predictions of the flowfield, the accuracy with which the computation predicts blade interaction was determined. It was found that in addition to upstream vane-rotor and rotor-downstream vane interactions, a new interaction mechanism was found resulting from the interaction between the downstream vane's potential field and the upstream vane's trailing edge potential field and shock.

Wake, shock and potential field interactions in a 1.5 stage turbine: Part II: Vane-vane interaction and discussion of results

The composition of the time-resolved surface pressure field around a high-pressure rotor blade caused by the presence of neighboring blade rows was studied, with the individual effects of wake, shock and potential field interaction being determined. Two test geometries were considered: first, a high-pressure turbine stage coupled with a swan-necked diffuser exit duct; secondly, the same high-pressure stage but with a vane located in the downstream duct. Both tests were carried out at engine-representative Mach and Reynolds numbers. By comparing the results to time-resolved computational predictions of the flowfield, the accuracy with which the computation predicts blade interaction was determined. Evidence was obtained that for a large downstream vane, the flow conditions in the rotor passage, at any instant in time, are close to being steady state.

The development of turbine exit flow in a swan-necked inter-stage diffuser

This paper describes both the migration and dissipation of flow phenomena downstream of a transonic high-pressure turbine stage. The geometry of the HP stage exit duct considered is a swan-necked diffuser similar to those likely to be used in future engine designs. The paper contains results both from an experimental programme in a turbine test facility and from numerical predictions. Experimental data was acquired using three fast-response aerodynamic probes capable of measuring Mach number, whirl angle, pitch angle, total pressure and static pressure. The probes were used to make time-resolved area traverses at two axial locations downstream of the rotor trailing edge. A 3D time-unsteady viscous Navier-Stokes solver was used for the numerical predictions. The unsteady exit flow from a turbine stage is formed from rotordependent phenomena (such as the rotor wake, the rotor trailing edge recompression shock, the tip-leakage flow and the hub secondary flow) and vane-rotor interaction dependant phenomena. This paper describes the time-resolved behaviour and three-dimensional migration paths of both of these phenomena as they convect downstream. It is shown that the inlet flow to a downstream vane is dominated by two corotating vortices, the first caused by the rotor tip-leakage flow and the second by the rotor hub secondary flow. At the inlet plane of the downstream vane the wake is extremely weak and the radial pressure gradient is shown to have caused the majority of the high loss wake fluid to be located between the mid-height of the passage and the casing wall. The structure of the flow indicates that between a high pressure stage and a downstream vane simple two-dimensional blade row interaction does not occur. The results presented in this paper indicate that the presence of an upstream stage is likely to significantly alter the structure of the secondary flow within a downstream vane. The paper also shows that vane-rotor interaction within the upstream stage causes a 10° circumferential variation in the inlet flow angle of the 2nd stage vane.

Axial compressor response to inlet flow distortions by a CFD analysis

The usual approach to compressor design considers uniform inlet flow characteristics. Especially in aircraft applications, the inlet flow is quite often non uniform, and this can result in severe performance degradation. The magnitude of this phenomenon is amplified in military engines due to the complexity of inlet duct configurations and the extreme flight conditions. CFD simulation is an innovative and powerful tool for studying inlet distortions and can bring this inside the very early phases of the design process. This project attempts to study the effects of inlet flow distortions in an axial flow compressor trying to minimize the use computer resources and computational time. The first stage of a low bypass ratio compressor has been analyzed and its clean and distorted performance compared outlining the principal changes due to uneven flow distribution: drop in mass flow, increase in pressure and temperature ratios, decrease in surge margin. Three different studies have then been conducted to better understand the effects of the level, the type and the frequency of the distortion.

Deposition of small particles from turbulent flows

This paper deals with particle deposition onto solid walls from turbulent flows. The aim of the study is to model particle deposition in industrial flows, such as the one in gas turbines. The numerical study has been carried out with a two fluid approach. The possible contribution to the deposition from Brownian diffusion, turbulent diffusion and shear-induced lift force are considered in the study. Three types of turbulent two-phase flows have been studied: turbulent channel flow, turbulent flow in a bent duct and turbulent flow in a turbine blade cascade. In the turbulent channel flow case, the numerical results from a two-dimensional code show good agreement with numerical and experimental results from other resources. Deposition problem in a bent duct flow is introduced to study the effect of curvature. Finally, the deposition of small particles on a cascade of turbine blades is simulated. The results show that the current two fluid models are capable of predicting particle deposition rates in complex industrial flows.

Performance of a Boundary Layer Ingesting (BLI) propulsion system

This paper presents an assessment of the performance of an embedded propulsion system in the presence of distortion associated with boundary layer ingestion. For fan pressure ratios of interest for civil transports, the benefits of boundary layer ingestion are shown to be very sensitive to the magnitude of fan and duct losses. The distortion transfer across the fan, basically the comparison of the stagnation pressure non-uniformity downstream of the fan to that upstream of the fan, has a major role in determining the impact of boundary layer ingestion on overall fuel burn. This, in turn, puts requirements on the fidelity with which one needs to assess the distortion transfer, and thus the type of models that need to be used in such assessment. For the three-dimensional distortions associated with fuselage boundary layers ingested into a subsonic diffusing inlet, it is found that boundary layer ingestion can provide decreases in fuel burn of several per cent. It is also shown that a promising avenue for mitigating the risks (aerodynamic as well as aeromechanical) in boundary layer ingestion is to mix out the flow before it reaches the engine face.

An experimental study of unsteady separated shock wave boundary layer interactions

An experimental investigation of the unsteady interaction between a turbulent boundary layer and a normal shock wave of strength M∞ = 1.4 subject to periodic forcing in a parallel walled duct has been conducted. Emphasis has been placed on the mechanism by which changes in the global flow field influence the local interaction structure. Static pressure measurements and high speed Schlieren images of the unsteady interaction have been obtained. The pressure rise across the interaction and the appearance of the local SBLI structure have been observed to vary during the cycle of periodic shock wave motion. The magnitude of the pressure rise across the interaction is found to be related to the relative Mach number of the unsteady shock wave as it undergoes periodic motion. Variations in the upstream Influence of the interaction are sensitive to the magnitude and direction of shock wave velocity and acceleration and it is proposed that a viscous lag exists between the point of boundary layer separation and the shock wave position. Further work exploring the implications of these findings is proposed, including studies of the variation in position of the points of boundary layer separation and reattachment.

The prediction of hot streak migration in a high-pressure turbine

Accurate predictions of combustor hot streak migration enable the turbine designer to identify high-temperature regions that can limit component life. It is therefore important that these predictions are achieved within the short time scales of a design process. This article compares temperature measurements of a circular hot streak through a turning duct and a research turbine with predictions using a three-dimensional Reynolds-averaged Navier-Stokes solver. It was found that the mixing length turbulence model did not predict the hot streak dissipation accurately. However, implementation of a very simple model of the free stream turbulence (FST) significantly improved the exit temperature predictions on both the duct and research turbine. One advantage of the simple FST model described over more complex alternatives is that no additional equations are solved. This makes the method attractive for design purposes, as it is not associated with any increase in computational time.

High frequency multimode radiation from ducts with flow

Multimode sound radiation from hard-walled semi-infinite ducts with uniform subsonic flow is investigated theoretically. An analytic expression, valid in the high frequency limit, is derived for the multimode directivity function in the forward arc for a general family of mode distribution functions. The multimode directivity depends on the amplitude and directivity function of each individual mode. The amplitude of each mode is expressed as a function of cut-off ratio for a uniform distribution of incoherent monopoles, a uniform distribution of incoherent axial dipoles and for equal power per mode. The modes' directivity functions are obtained analytically by applying a Lorentz transformation to the zero flow solution. The analytic formula for the multimode directivity with flow is derived assuming total transmission of power at the open-end of the duct. This formula is compared to the exact numerical result for an unflanged duct, computed utilizing a Wiener-Hopf solution. The agreement is shown to be excellent. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Behaviour of unsteady transonic shock / boundary layer interactions with three-dimensional effects

An experimental investigation into the response of transonic SBLIs to periodic down-stream pressure perturbations in a parallel walled duct has been conducted. Tests have been carried out with a shock strength of M ∞ = 1.5 for pressure perturbation frequencies in the range 16-90 Hz. Analysis of the steady interaction at M∞ = 1.5 has also been made. The principle measurement techniques were high speed schlieren photography and laser Doppler anemometry. The structure of the steady SBLI was found to be highly three-dimensional, with large corner flows and sidewall SBLIs. These aspects are thought to influence the upstream transmission of pressure information through the interaction by affecting the post-shock flow field, including the extent of regions of secondary supersonic flow. At low frequency, the dynamics of shock motion can be predicted using an inviscid analytical model. At increased frequencies, viscous effects become significant and the shock exhibits unexpected dynamic behaviour, due to a phase lag between the upstream transmission of pressure information in the core flow and in the viscous boundary layers. Flow control in the form of micro-vane vortex generators was found to have a small impact on shock dynamics, due to the effect it had on the post-shock flow field outside the viscous boundary layer region. The relationship between inviscid and viscous effects is developed and potential destabilising mechanisms for SBLIs in practical applications are suggested. Copyright © 2009 by Paul Bruce and Holger Babinsky.

Experimental investigation into parameters governing corner interactions for transonic shock wave/ boundary layer interactions

Experiments were conducted investigating the interaction between a normal shock wave and a corner boundary layer in a constant area rectangular duct. Active corner suction and passive blowing were applied to manipulate the natural corner flows developing in the working section of the Cambridge University supersonic wind tunnel. In addition robust vane micro-vortex generators were applied to the corners of the working section. Experiments were conducted at Mach numbers of M∞=1.4 and 1.5. Flow visualisation was carried out through schlieren and surface oil flow, while static pressures were recorded via floor tappings. The results indicate that an interplay occurs between the corner flow and the centre line flow. It is believed that corner flow separation acts to induce a shock bifurcation, which in turn leads to a smearing of the adverse pressure gradient elsewhere. In addition the blockage effect from the corners was seen to result in a reacceleration of the subsonic post-shock flow. As a result manipulation of the corner regions allows a separated or attached centre line flow to be observed at the same Mach number. Copyright © 2010 by Babinsky, Burton, Bruce.