Effects of Boundary Layer Ingestion on the Aero-Acoustics of Transonic Fan Rotors

The use of boundary-layer-ingesting, embedded propulsion systems can result in inlet flow distortions where the interaction of the boundary layer vorticity and the inlet lip causes horseshoe vortex formation and the ingestion of streamwise vortices into the inlet. A previously-developed body-force-based fan modeling approach was used to assess the change in fan rotor shock noise generation and propagation in a boundary-layer-ingesting, serpentine inlet. This approach is employed here in a parametric study to assess the effects of inlet geometry parameters (offset-to-diameter ratio and downstream-to-upstream area ratio) on flow distortion and rotor shock noise. Mechanisms related to the vortical inlet structures were found to govern changes in the rotor shock noise generation and propagation. The vortex whose circulation is in the opposite direction to the fan rotation (counter-swirling vortex) increases incidence angles on the fan blades near the tip, enhancing noise generation. The vortex with circulation in the direction of fan rotation (co-swirling vortex) creates a region of subsonic relative flow near the blade tip radius which decreases the sound power propagated to the far-field. The parametric study revealed that the overall sound power level at the fan leading edge is set by the ingested streamwise circulation, and that for inlet designs in which the streamwise vortices are displaced away from the duct wall, the sound power at the upstream inlet plane increased by as much as 9 dB. By comparing the far-field noise results obtained to those for a conventional inlet, it is deduced that the changes in rotor shock noise are predominantly due to the ingestion of streamwise vorticity.

A body force-based method for prediction of multiple-pure-tone noise: Validation

Emissions, fuel burn, and noise are the main drivers for innovative aircraft design. Embedded propulsion systems, such as for example used in hybrid-wing body aircraft, can offer fuel burn and noise reduction benefits but the impact of inlet flow distortion on the generation and propagation of turbomachinery noise has yet to be assessed. A novel approach is used to quantify the effects of non-uniform flow on the creation and propagation of multiple pure tone (MPT) noise. The ultimate goal is to conduct a parametric study of S-duct inlets to quantify the effects of inlet design parameters on the acoustic signature. The key challenge is that the effects of distortion transfer, noise source generation and propagation through the non-uniform flow field are inherently coupled such that a simultaneous computation of the aerodynamics and acoustics is required to capture the mechanisms at play. The technical approach is based on a body force description of the fan blade row that is able to capture the distortion transfer and the blade-to-blade flow variations that cause the MPT noise while reducing computational cost. A single, 3-D full-wheel CFD simulation, in which the Euler equations are solved to second-order spatial and temporal accuracy, simultaneously computes the MPT noise generation and its propagation in distorted inlet flow. A new method of producing the blade-to-blade variations in the body force field for MPT noise generation has been developed and validated. The numerical dissipation inherent to the solver is quantified and used to correct for non-physical attenuation in the far-field noise spectra. Source generation, acoustic propagation and acoustic energy transfer between modes is examined in detail. The new method is validated on NASA's Source Diagnostic Test fan and inlet, showing good agreement with experimental data for aerodynamic performance, acoustic source generation, and far-field noise spectra. The next steps involve the assessment of MPT noise in serpentine inlet ducts and the development of a reduced order formulation suitable for incorporation into NASA's ANOPP framework. © 2010 by Jeff Defoe, Alex Narkaj & Zoltan Spakovszky.

When will a flame describing function approach to thermoacoustics work well?

In any thermoacoustic analysis, it is important not only to predict linear frequencies and growth rates, but also the amplitude and frequencies of any limit cycles. The Flame Describing Function (FDF) approach is a quasi-linear analysis which allows the prediction of both the linear and nonlinear behaviour of a thermoacoustic system. This means that one can predict linear growth rates and frequencies, and also the amplitudes and frequencies of any limit cycles. The FDF achieves this by assuming that the acoustics are linear and that the flame, which is the only nonlinear element in the thermoacoustic system, can be adequately described by considering only its response at the frequency at which it is forced. Therefore any harmonics generated by the flame's nonlinear response are not considered. This implies that these nonlinear harmonics are small or that they are sufficiently filtered out by the linear dynamics of the system (the low-pass filter assumption). In this paper, a flame model with a simple saturation nonlinearity is coupled to simple duct acoustics, and the success of the FDF in predicting limit cycles is studied over a range of flame positions and acoustic damping parameters. Although these two parameters affect only the linear acoustics and not the nonlinear flame dynamics, they determine the validity of the low-pass filter assumption made in applying the flame describing function approach. Their importance is highlighted by studying the level of success of an FDF-based analysis as they are varied. This is achieved by comparing the FDF's prediction of limit-cycle amplitudes to the amplitudes seen in time domain simulations.

A six sensor method for measuring acoustic properties in ducts.

An accurate description of sound propagation in a duct is important to obtain the sound power radiating from a source in both near and far fields. A technique has been developed and applied to decompose higher-order modes of sound emitted into a duct. Traditional experiments and theory based on two-sensor methods are limited to the plane-wave contribution to the sound field at low frequency. Due to the increase in independent measurements required, a computational method has been developed to simulate sensitivities of real measurements (e.g., noise) and optimize the set-up. An experimental rig has been constructed to decompose the first two modes using six independent measurements from surface, flush-mounted microphones. Experiments were initially performed using a loudspeaker as the source for validation. Subsequently, the sound emitted by a mixed-flow fan has been investigated and compared to measurements made in accordance with the internationally standardized in-duct fan measurement method. This method utilizes large anechoic terminations and a procedure involving averaging over measurements in space and time to account for the contribution from higher-order modes. The new method does not require either of these added complications and gives detail about the underlying modal content of the emitted sound.

Flexible conversion ratio fast reactors: Overview

Conceptual designs of lead-cooled and liquid salt-cooled fast flexible conversion ratio reactors were developed. The performance achievable by the unity conversion ratio cores of these reactors was compared to an existing supercritical carbon dioxide-cooled (S-CO2) fast reactor design and an uprated version of an existing sodium-cooled fast reactor. All concepts have cores rated at 2400 MWt. The cores of the liquid-cooled reactors are placed in a large-pool-type vessel with dual-free level, which also contains four intermediate heat exchangers (IHXs) coupling a primary coolant to a compact and efficient supercritical CO2 Brayton cycle power conversion system. The S-CO2 reactor is directly coupled to the S-CO2 Brayton cycle power conversion system. Decay heat is removed passively using an enhanced reactor vessel auxiliary cooling system (RVACS) and a passive secondary auxiliary cooling system (PSACS). The selection of the water-cooled versus air-cooled heat sink for the PSACS as well as the analysis of the probability that the PSACS may fail to complete its mission was performed using risk-informed methodology. In addition to these features, all reactors were designed to be self-controllable. Further, the liquid-cooled reactors utilized common passive decay heat removal systems whereas the S-CO2 uses reliable battery powered blowers for post-LOCA decay heat removal to provide flow in well defined regimes and to accommodate inadvertent bypass flows. The multiple design limits and challenges which constrained the execution of the four fast reactor concepts are elaborated. These include principally neutronics and materials challenges. The neutronic challenges are the large positive coolant reactivity feedback, small fuel temperature coefficient, small effective delayed neutron fraction, large reactivity swing and the transition between different conversion ratio cores. The burnup, temperature and fluence constraints on fuels, cladding and vessel materials are elaborated for three categories of material - materials currently available, available on a relatively short time scale and available only with significant development effort. The selected fuels are the metallic U-TRU-Zr (10% Zr) for unity conversion ratio and TRU-Zr (75% Zr) for zero conversion ratio. The principal selected cladding and vessel materials are HT-9 and A533 or A508, respectively, for current availability, T-91 and 9Cr-1Mo steel for relatively short-term availability and oxide dispersion strengthened ferritic steel (ODS) available only with significant development. © 2009 Elsevier B.V. All rights reserved.

Plane-strain discrete dislocation plasticity with climb-assisted glide motion of dislocations

A small-strain two-dimensional discrete dislocation plasticity (DDP) framework is developed wherein dislocation motion is caused by climb-assisted glide. The climb motion of the dislocations is assumed to be governed by a drag-type relation similar to the glide-only motion of dislocations: such a relation is valid when vacancy kinetics is either diffusion limited or sink limited. The DDP framework is employed to predict the effect of dislocation climb on the uniaxial tensile and pure bending response of single crystals. Under uniaxial tensile loading conditions, the ability of dislocations to bypass obstacles by climb results in a reduced dislocation density over a wide range of specimen sizes in the climb-assisted glide case compared to when dislocation motion is only by glide. A consequence is that, at least in a single slip situation, size effects due to dislocation starvation are reduced. By contrast, under pure bending loading conditions, the dislocation density is unaffected by dislocation climb as geometrically necessary dislocations (GNDs) dominate. However, climb enables the dislocations to arrange themselves into lower energy configurations which significantly reduces the predicted bending size effect as well as the amount of reverse plasticity observed during unloading. The results indicate that the intrinsic plasticity material length scale associated with GNDs is strongly affected by thermally activated processes and will be a function of temperature. © 2013 IOP Publishing Ltd.

Large eddy simulation of turbine internal cooling ducts

Abstract Large-Eddy Simulation (LES) and hybrid Reynolds-averaged Navier–Stokes–LES (RANS–LES) methods are applied to a turbine blade ribbed internal duct with a 180° bend containing 24 pairs of ribs. Flow and heat transfer predictions are compared with experimental data and found to be in agreement. The choice of LES model is found to be of minor importance as the flow is dominated by large geometric scale structures. This is in contrast to several linear and nonlinear RANS models, which display turbulence model sensitivity. For LES, the influence of inlet turbulence is also tested and has a minor impact due to the strong turbulence generated by the ribs. Large scale turbulent motions destroy any classical boundary layer reducing near wall grid requirements. The wake-type flow structure makes this and similar flows nearly Reynolds number independent, allowing a range of flows to be studied at similar cost. Hence LES is a relatively cheap method for obtaining accurate heat transfer predictions in these types of flows.

Low-order modelling of the response of ducted flames in annular geometries

The adoption of lean premixed prevaporised combustion systems can reduce NOx emissions from gas turbines, but unfortunately also increases their susceptibility to thermoacoustic instabilities. Initially, acoustic waves can produce heat release fluctuations by a variety of mechanisms, often by perturbing the equivalence ratio. If correctly phased, heat release fluctuations can subsequently generate more acoustic waves, which at high amplitude can result in significant structural damage to the combustor. The prediction of this phenomenon is of great industrial interest. In previous work, we have coupled a physics based, kinematic model of the flame with a network model to provide the planar acoustic response necessary to close the feedback loop and predict the onset and amplitude of thermoacoustic instabilities in a lab-scale, axisymmetric single burner combustor. The advantage of a time domain approach is that the modal interaction, the influence of harmonics, and flame saturation can be investigated. This paper extends this approach to more realistic, annular geometries, where both planar and circumferential modes must be considered. In lean premixed prevaporised combustors, fluctuations in equivalence ratio have been shown to be a dominant cause of unsteady combustion. These can occur, for example, due to velocity perturbations in the premix ducts, which can lead to equivalence ratio fluctuations at the fuel injectors, which are subsequently convected downstream to the flame surfaces. Here, they can perturb the heat release by locally altering the flame speed, enthalpy of combustion, and, indirectly, the flame surface area. In many gas turbine designs, particularly aeroengines, the geometries are composed of a ring of premix ducts linking a plenum and an annular combustor. The most unstable modes are often circumferential modes. The network model is used to characterise the flow response of the geometry to heat fluctuations at an appropriate location, such as the fuel injectors. The heat release at each flame holder is determined in the time domain using the kinematic flame model derived, as a function of the flow perturbations in the premix duct. This approach is demonstrated for an annular ring of burners on a in a simple geometry. The approach is then extended to an industrial type gas turbine combustor, and used to predict the limit cycle amplitudes. Copyright © 2012 by ASME.

Acoustic analogy in swirling mean flow applied to predict rotor trailing-edge noise

The present study aims at accounting for swirling mean flow effects on rotor trailing-edge noise. Indeed, the mean flow in between the rotor and the stator of the fan or of a compressor stage is highly swirling. The extension of Ffowcs-Williams & Hawkings' acoustic analogy in a medium at rest with moving surfaces and of Goldstein's acoustic analogy in a circular duct with uniform mean flow to a swirling mean flow in an annular duct is introduced. It is first applied to tonal noise. In most cases, the swirl modifies the pressure distribution downstream of the fan. In several configurations, when the swirl is rather close to a solid body swirl, it is often sufficient to apply a simple Doppler effect correction when predicting the duct modes in uniform mean flow in order to predict accurately the noise radiated with swirl. However, in other realistic configurations, the swirling mean-flow effect cannot be addressed using this simple Doppler effect correction. Second, a rotor trailing-edge noise model accounting for both the effects of the annular duct and the swirling mean flow is developed and applied to a realistic fan rotor with different swirling and sheared mean flows (and as a result different associated blade stagger angles). The benchmark cases are built from the Boeing 18-inch Fan Rig Broadband Noise Test. In all cases the swirling mean flow has an effect. In some cases the a simple Doppler effect may address it, but, in other realistic configurations our acoustic analogy with swirl is needed. © 2012 by the authors. Published by the American Institute of Aeronautics and Astronautics, Inc.

Angle-of-attack effects on counter-rotating propellers at take-off

Due to their potential for significant fuel consumption savings, Counter-Rotating Open Rotors (CRORs) are currently being considered as an alternative to high-bypass turbofans. When CRORs are mounted on an aircraft, several 'installation effects' arise which are not present when the engine is operated in isolation. This paper investigates how flow features arising from one such effect - The angle-of-attack of the engine centre-line relative to the oncoming flow - can influence the design of CROR engines. Three-dimensional full-annulus unsteady CFD simulations are used to predict the time-varying flow field experienced by each rotor and emphasis is put on the interaction of the frontrotor wake and tip vortex with the rear-rotor. A parametric study is presented that quantifies the rotorrotor interaction as a function of the angle-of-attack. It is shown that angle-of-attack operation significantly changes the flow field and the unsteady lift on both rotors. In particular, a frequency analysis shows that the unsteady lift exhibits sidebands around the rotor-rotor interaction frequencies. Further, a non-linear increase in the total rear-rotor tip unsteadiness is observed for moderate and high angles-of-attack. The results presented in this paper demonstrate that common techniques used to mitigate CROR noise, such as modifying the rotor-rotor axial spacing and rear-rotor crop, can not be applied correctly unless angle-of-attack effects are taken into account. Copyright © 2012 by ASME.

Nonlinear thermoacoustics of ducted premixed flames: The influence of perturbation convection speed

When a premixed flame is placed within a duct, acoustic waves induce velocity perturbations at the flame's base. These travel down the flame, distorting its surface and modulating its heat release. This can induce self-sustained thermoacoustic oscillations. Although the phase speed of these perturbations is often assumed to equal the mean flow speed, experiments conducted in other studies and Direct Numerical Simulation (DNS) conducted in this study show that it varies with the acoustic frequency. In this paper, we examine how these variations affect the nonlinear thermoacoustic behaviour. We model the heat release with a nonlinear kinematic G-equation, in which the velocity perturbation is modelled on DNS results. The acoustics are governed by linearised momentum and energy equations. We calculate the flame describing function (FDF) using harmonic forcing at several frequencies and amplitudes. Then we calculate thermoacoustic limit cycles and explain their existence and stability by examining the amplitude-dependence of the gain and phase of the FDF. We find that, when the phase speed equals the mean flow speed, the system has only one stable state. When the phase speed does not equal the mean flow speed, however, the system supports multiple limit cycles because the phase of the FDF changes significantly with oscillation amplitude. This shows that the phase speed of velocity perturbations has a strong influence on the nonlinear thermoacoustic behaviour of ducted premixed flames. © 2013 The Combustion Institute.

Models for acoustic propagation through turbofan exhaust flows

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 previous solutions is that they do not allow for the jet nozzle to protrude a finite distance beyond the end of the nacelle (or in certain configurations being buried a finite distance upstream). In this paper we use the matrix Wiener-Hopf technique, which will allow precisely the finite nacelle-jet nozzle separation to be included. The crucial step in our work is to factorise a certain matrix as a product of terms analytic and invertible in the upper/lower halves of the complex plane. The way we do this matrix factorisation is quite different in the buried and protruding nozzle cases. In the buried case our solution method is the so-called pole-removal technique. In the technically more demanding protruding case, however, we must first use Pade approximants to generate a uniformly-valid, meromorphic representation of a certain function, before the same pole-removal method can be applied. Sample results are presented, investigating in particular the effects of exit plane stagger. © 2007 by B Veitch and N Peake.

Aeroacoustic source measurement methods for characterizing the sound generated by ducted flow devices with higher-order modes.

The International Organization for Standardization (ISO) method 5136 is widely used in industry and academia to determine the sound power radiated into a duct by fans and other flow devices. The method involves placing the device at the center of a long cylindrical duct with anechoic terminations at each end to eliminate reflections. A single off-axis microphone is used on the inlet and outlet sides that can theoretically capture the plane-wave mode amplitudes but this does not provide enough information to fully account for higher-order modes. In this study, the "two-port" source model is formulated to include higher-order modes and applied for the first three modes. This requires six independent surface pressure measurements on each side or "port." The resulting experimental set-up is much shorter than the ISO rig and does not require anechoic terminations. An array of six external loudspeaker sources is used to characterize the passive part of the two-port model and the set-up provides a framework to account for transmission of higher-order modes through a fan. The relative importance of the higher-order modes has been considered and their effect on inaccuracies when using the ISO method to find source sound power has been analyzed.

Aerodynamic design of high end wall angle turbine stages-part I: Methodology development

A design methodology is presented for turbines in an annulus with high end wall angles. Such stages occur where large radial offsets between the stage inlet and stage outlet are required, for example in the first stage of modern low pressure turbines, and are becoming more prevalent as bypass ratios increase. The turbine vanes operate within s-shaped ducts which result in meridional curvature being of a similar magnitude to the bladeto-blade curvature. Through a systematic series of idealized computational cases, the importance of two aspects of vane design are shown. First, the region of peak end wall meridional curvature is best located within the vane row. Second, the vane should be leant so as to minimize spanwise variations in surface pressure-this condition is termed "ideal lean." This design philosophy is applied to the first stage of a low pressure turbine with high end wall angles. © 2014 by ASME.

Measurements of intake separation hysteresis in a model fan and nacelle RIG

A 1/20-scale, low speed model rig representing the fan and nacelle of a high bypass ratio jet engine has been tested under crosswind conditions. The flow conditions under which the intake flow separates and reattaches have been found to exhibit considerable hysteresis. This phenomenon has been examined by a careful test procedure in which the crosswind angle has been slowly increased and then decreased. Measurements of the hysteresis associated with separation and reattachment are presented for independent variations in stream-tube contraction ratio, ground clearance, fan operating point and Reynolds number. The results reveal that particular care must be taken to allow for any hysteresis when testing intakes under crosswind conditions. They also indicate that separation hysteresis is particularly sensitive to fan operating point and the position of the ground plane. These findings suggest that it is important for high Reynolds number intake tests and calculations to include a ground plane and a model of the downstream turbomachinery. © 2002 by the author(s).