Large-eddy simulation of complex geometry jets

Computations are made for chevron and coflowing jet nozzles. The latter has a bypass ratio of 6:1. Also, unlike the chevron nozzle, the core flow is heated, making the inlet conditions reminiscent of those for a real engine. A large-eddy resolving approach is used with circa 12 × 10 6 cell meshes. Because the codes being used tend toward being dissipative the subgrid scale model is abandoned, giving what can be termed numerical large-eddy simulation. To overcome near-wall modeling problems a hybrid numerical large-eddy simulation-Reynolds-averaged Navier-Stokes related method is used. For y + ≤ 60 a Reynolds-averaged Navier-Stokes model is used. Blending between the two regions makes use of the differential Hamilton-Jabobi equation, an extension of the eikonal equation. For both nozzles, results show encouraging agreement with measurements of other workers. The eikonal equation is also used for ray tracing to explore the effect of the mean flow on acoustic ray trajectories, thus yielding a coherent solution strategy. © 2011 by Cambridge University.

Large eddy simulation of a hot, high speed coflowing jet

Computations are made of a short cowl coflowing jet nozzle with a bypass ratio 8 : 1. The core flow is heated, making the inlet conditions reminiscent of those for a real engine. A large eddy resolving approach is used with a 12 × 106 cell mesh. Since the code being used tends towards being dissipative the sub-grid scale (SGS) model is abandoned giving what can be termed Numerical Large Eddy Simulation (NLES). To overcome near wall modelling problems a hybrid NLES-RANS (Reynolds Averaged Navier-Stokes) related method is used. For y+ ≤ 60 a κ-l model is used. Blending between the two regions makes use of the differential Hamilton-Jabobi (HJ) equation, an extension of the eikonal equation. Results show encouraging agreement with existing measurements of other workers. The eikonal equation is also used for acoustic ray tracing to explore the effect of the mean flow on acoustic ray trajectories, thus yielding a coherent solution strategy. Copyright © 2011 by ASME.

Large-eddy simulations and measurements of a small-scale high-speed coflowing jet

Measurements and predictions are made of a short-cowl coflowing jet with a bypass ratio of 8:1. The Reynolds number is 300,000, and the inlet Mach numbers are representative of aeroengine conditions. The low Reynolds number of the measurements makes the case well suited to the assessment of large-eddy-simulation-related strategies. The nozzle concentricity is carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both laser Doppler anemometry and particle image velocimetry. The simulations are completed on 6× 106, 12× 106, and 50 × 106 cell meshes. To overcome near-wall modeling problems, a hybrid large-eddy-simulation-Reynolds-averaged-Navier-Stokesrelated method is used. The near-wall Reynolds-averaged-Navier-Stokes layer is helpful in preventing nonphysical separation from the nozzle wall.Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Comparison of les to LDA and PIV measurements of a small scale high speed coflowing jet

Measurements and predictions are made of a short cowl co-flowing jet with a bypass ratio of 8:1. The Reynolds number for computations and measurements are matched at 300,000 and the Mach numbers representative of realistic jet conditions with core and co flow velocities of 240m/s and 216m/s respectively. The low Reynolds number of the measurements makes the case well suited to the assessment of large eddy resolving computational strategies. Also, the nozzle concentricity was carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both two dimensional coincident LDA and PIV systems. The computational simulations are completed on a modest 12×106 mesh. The simulation is now being run on a 50×106 mesh using hybrid RANSNLES (Numerical Large Eddy Simulation). Close to the nozzle wall a k-l RANS model is used. For an axisymmetric jet, comparison is made between simulations which use NLES, RANSNLES and also a simple imposed velocity profile where the nozzle is not modeled. The use of a near wall RANS model is shown to be beneficial. When compared with the measurements the NLES results are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

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.

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.