A Study of Fan-Distortion Interaction Within the NASA Rotor 67 Transonic Stage

The performance of a transonic fan operating within nonuniform inlet flow remains a key concern for the design and operability of a turbofan engine. This paper applies computational methods to improve the understanding of the interaction between a transonic fan and an inlet total pressure distortion. The test case studied is the NASA rotor 67 stage operating with a total pressure distortion covering a 120-deg sector of the inlet flow field. Full-annulus, unsteady, three-dimensional CFD has been used to simulate the test rig installation and the full fan assembly operating with inlet distortion. Novel post-processing methods have been applied to extract the fan performance and features of the interaction between the fan and the nonuniform inflow. The results of the unsteady computations agree well with the measurement data. The local operating condition of the fan at different positions around the annulus has been tracked and analyzed, and this is shown to be highly dependent on the swirl and mass flow redistribution that the rotor induces ahead of it due to the incoming distortion. The upstream flow effects lead to a variation in work input that determines the distortion pattern seen downstream of the fan stage. In addition, the unsteady computations also reveal more complex flow features downstream of the fan stage, which arise due to the three dimensionality of the flow and unsteadiness. © 2012 American Society of Mechanical Engineers.

A validation of commercial flow solver performance on hybridised octree meshes

In recent years we have been developing a meshing system which is aimed at eliminating the bottleneck represented by building meshes for real-world, complex turbomachinery configurations. This system is based on a rapid octree meshing technology which is then made conformal to the bodies present. The objective of this paper is to demonstrate that this class of mesh is not only very fast to produce but also fit-for-purpose in the sense that simulations generated with third-party commercial flow solvers like Fluent have the same accuracy as those performed on more conventional meshes. A range of standard examples and test cases will be presented. Copyright © 2011 by ASME.

A Physical Interpretation of Stagnation Pressure and Enthalpy Changes in Unsteady Flow

This paper provides a physical interpretation of the mechanism of stagnation enthalpy and stagnation pressure changes in turbomachines due to unsteady flow, the agency for all work transfer between a turbomachine and an inviscid fluid. Examples are first given to illustrate the direct link between the time variation of static pressure seen by a given fluid particle and the rate of change of stagnation enthalpy for that particle. These include absolute stagnation temperature rises in turbine rotor tip leakage flow, wake transport through downstream blade rows, and effects of wake phasing on compressor work input. Fluid dynamic situations are then constructed to explain the effect of unsteadiness, including a physical interpretation of how stagnation pressure variations are created by temporal variations in static pressure; in this it is shown that the unsteady static pressure plays the role of a time-dependent body force potential. It is further shown that when the unsteadiness is due to a spatial nonuniformity translating at constant speed, as in a turbomachine, the unsteady pressure variation can be viewed as a local power input per unit mass from this body force to the fluid particle instantaneously at that point. © 2012 American Society of Mechanical Engineers.

Stall Warning by Blade Pressure Signature Analysis

At low mass flow rates, axial compressors suffer from flow instabilities leading to stall and surge. The inception process of these instabilities has been widely researched in the past---primarily with the aim of predicting or averting stall onset. In recent times, attention has shifted to conditions well before stall and has focused on the level of irregularity in the blade passing signature in the rotor tip region. In general, the irregularity increases in intensity as the flow rate through the compressor is reduced. Attempts have been made to develop stall warning/avoidance procedures based on the level of flow irregularity, but little effort has been made to characterize the irregularity itself, or to understand its underlying cause. Work on this project has revealed for the first time that the increase in irregularity in the blade passing signature is highly dependent on both tip-clearance size and eccentricity. In a compressor with small, uniform, tip-clearance, the increase in blade passing irregularity that accompanies a reduction in flow rate will be modest. If the tip-clearance is enlarged, however, there will be a sharp rise in irregularity at all circumferential locations. In a compressor with eccentric tip-clearance, the increase in irregularity will only occur in the part of the annulus where the tip-clearance is largest, regardless of the average clearance level. In this paper, some attention is also given to the question of whether the irregularity observed in the prestall flow field is due to random turbulence or to some form of coherent flow structure. Detailed flow measurements reveal that the latter is the case. From these findings, it is clear that a stall warning system based on blade passing signature irregularity would be difficult to implement in an aero-engine where tip-clearance size and eccentricity change during each flight cycle and over the life of the compressor. © 2013 American Society of Mechanical Engineers.

Shock Propagation and MPT Noise From a Transonic Rotor in Nonuniform Flow

One of the major challenges in high-speed fan stages used in compact, embedded propulsion systems is inlet distortion noise. A body-force-based approach for the prediction of multiple-pure-tone (MPT) noise was previously introduced and validated. In this paper, it is employed with the objective of quantifying the effects of nonuniform flow on the generation and propagation of MPT noise. First-of-their-kind back-to-back coupled aero-acoustic computations were carried out using the new approach for conventional and serpentine inlets. Both inlets delivered flow to the same NASA/GE R4 fan rotor at equal corrected mass flow rates. Although the source strength at the fan is increased by 38 dB in sound power level due to the nonuniform inflow, far-field noise for the serpentine inlet duct is increased on average by only 3.1 dBA overall sound pressure level in the forward arc. This is due to the redistribution of acoustic energy to frequencies below 11 times the shaft frequency and the apparent cut-off of tones at higher frequencies including blade-passing tones. The circumferential extent of the inlet swirl distortion at the fan was found to be two blade pitches, or 1/11th of the circumference, suggesting a relationship between the circumferential extent of the inlet distortion and the apparent cut-off frequency perceived in the far field. A first-principles-based model of the generation of shock waves from a transonic rotor in nonuniform flow showed that the effects of nonuniform flow on acoustic wave propagation, which cannot be captured by the simplified model, are more dominant than those of inlet flow distortion on source noise. It demonstrated that nonlinear, coupled aerodynamic and aero-acoustic computations, such as those presented in this paper, are necessary to assess the propagation through nonuniform mean flow. A parametric study of serpentine inlet designs is underway to quantify these propagation effects. © 2013 American Society of Mechanical Engineers.

Shroud heat transfer measurements from a rotating cavity with an axial throughflow of air

The paper discusses measurements of heat transfer obtained from the inside surface of the peripheral shroud. The experiments were carried out on a rotating cavity, comprising two 0.985-m-dia disks, separated by an axial gap of 0.065 m and bounded at the circumference by a carbon fiber shroud. Tests were conducted with a heated shroud and either unheated or heated disks. When heated, the disks had the same temperature level and surface temperature distribution. Two different temperature distributions were tested; the surface temperature either increased, or decreased with radius. The effects of disk, shroud, and air temperature levels were also studied. Tests were carried out for the range of axial throughflow rates and speeds: 0.0025 ≤ m ≤ 0.2 kg/s and 12.5 ≤ Ω ≤ 125 rad/s, respectively. Measurements were also made of the temperature of the air inside the cavity. The shroud Nusselt numbers are found to depend on a Grashof number, which is defined using the centripetal acceleration. Providing the correct reference temperature is used, the measured Nusselt numbers also show similarity to those predicted by an established correlation for a horizontal plate in air. The heat transfer from the shroud is only weakly affected by the disk surface temperature distribution and temperature level. The heat transfer from the shroud appears to be affected by the Rossby number. A significant enhancement to the rotationally induced free convection occurs in the regions 2 ≤ Ro ≤ 4 and Ro ≥ 20. The first of these corresponds to a region where vortex breakdown has been observed. In the second region, the Rossby number may be sufficiently large for the central throughflow to affect the shroud heat transfer directly. Heating the shroud does not appear to affect the heat transfer from the disks significantly.

Blade Loading and Its Application in the Mean-Line Design of Low Pressure Turbines

In order to minimize the number of iterations to a turbine design, reasonable choices of the key parameters must be made at the preliminary design stage. The choice of blade loading is of particular concern in the low pressure (LP) turbine of civil aero engines, where the use of high-lift blades is widespread. This paper considers how blade loading should be measured, compares the performance of various loss correlations, and explores the impact of blade lift on performance and lapse rates. To these ends, an analytical design study is presented for a repeating-stage, axial-flow LP turbine. It is demonstrated that the long-established Zweifel lift coefficient (Zweifel, 1945, "The Spacing of Turbomachine Blading, Especially with Large Angular Deflection" Brown Boveri Rev., 32(1), pp. 436-444) is flawed because it does not account for the blade camber. As a result the Zweifel coefficient is only meaningful for a fixed set of flow angles and cannot be used as an absolute measure of blade loading. A lift coefficient based on circulation is instead proposed that accounts for the blade curvature and is independent of the flow angles. Various existing profile and secondary loss correlations are examined for their suitability to preliminary design. A largely qualitative comparison demonstrates that the loss correlations based on Ainley and Mathieson (Ainley and Mathieson, 1957, "A Method of Performance Estimation for Axial-Flow Turbines," ARC Reports and Memoranda No. 2974; Dunham and Came, 1970, "Improvements to the Ainley-Mathieson Method of Turbine Performance Prediction," Trans. ASME: J. Eng. Gas Turbines Power, July, pp. 252-256; Kacker and Okapuu, 1982, "A Mean Line Performance Method for Axial Flow Turbine Efficiency," J. Eng. Power, 104, pp. 111-119). are not realistic, while the profile loss model of Coull and Hodson (Coull and Hodson, 2011, "Predicting the Profile Loss of High-Lift Low Pressure Turbines," J. Turbomach., 134(2), pp. 021002) and the secondary loss model of (Traupel, W, 1977, Thermische Turbomaschinen, Springer-Verlag, Berlin) are arguably the most reasonable. A quantitative comparison with multistage rig data indicates that, together, these methods over-predict lapse rates by around 30%, highlighting the need for improved loss models and a better understanding of the multistage environment. By examining the influence of blade lift across the Smith efficiency chart, the analysis demonstrates that designs with higher flow turning will tend to be less sensitive to increases in blade loading. © 2013 American Society of Mechanical Engineers.

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.

Large Eddy Simulation for Turbines: Methodologies, Cost and Future Outlooks

Flows throughout different zones of turbines have been investigated using large eddy simulation (LES) and hybrid Reynolds-averaged Navier–Stokes-LES (RANS-LES) methods and contrasted with RANS modeling, which is more typically used in the design environment. The studied cases include low and high-pressure turbine cascades, real surface roughness effects, internal cooling ducts, trailing edge cut-backs, and labyrinth and rim seals. Evidence is presented that shows that LES and hybrid RANS-LES produces higher quality data than RANS/URANS for a wide range of flows. The higher level of physics that is resolved allows for greater flow physics insight, which is valuable for improving designs and refining lower order models. Turbine zones are categorized by flow type to assist in choosing the appropriate eddy resolving method and to estimate the computational cost.

Application of scalar dissipation rate model to siemens DLE combustors

The standard design process for the Siemens Industrial Turbomachinery, Lincoln, Dry Low Emissions combustion systems has adopted the Eddy Dissipation Model with Finite Rate Chemistry for reacting computational fluid dynamics simulations. The major drawbacks of this model have been the over-prediction of temperature and lack of species data limiting the applicability of the model. A novel combustion model referred to as the Scalar Dissipation Rate Model has been developed recently based on a flamelet type assumption. Previous attempts to adopt the flamelet philosophy with alternative closure models have failed, with the prediction of unphysical phenomenon. The Scalar Dissipation Rate Model (SDRM) was developed from a physical understanding of scalar dissipation rate, signifying the rate of mixing of hot and cold fluids at scales relevant to sustain combustion, in flames and was validated using direct numerical simulations data and experimental measurements. This paper reports on the first industrial application of the SDRM to SITL DLE combustion system. Previous applications have considered ideally premixed laboratory scale flames. The industrial application differs significantly in the complexity of the geometry, unmixedness and operating pressures. The model was implemented into ANSYS-CFX using their inbuilt command language. Simulations were run transiently using Scale Adaptive Simulation turbulence model, which switches between Large Eddy Simulation and Unsteady Reynolds Averaged Navier Stokes using a blending function. The model was validated in a research SITL DLE combustion system prior to being applied to the actual industrial geometry at real operating conditions. This system consists of the SGT-100 burner with a glass square-sectioned combustor allowing for detailed diagnostics. This paper shows the successful validation of the SDRM against time averaged temperature and velocity within measurement errors. The successful validation allowed application of the SDRM to the SGT-100 twin shaft at the relevant full load conditions. Limited validation data was available due to the complexity of measurement in the real geometry. Comparison of surface temperatures and combustor exit temperature profiles showed an improvement compared to EDM/FRC model. Furthermore, no unphysical phenomena were predicted. This paper presents the successful application of the SDRM to the industrial combustion system. The model shows a marked improvement in the prediction of temperature over the EDM/FRC model previously used. This is of significant importance in the future applications of combustion CFD for understanding of hardware mechanical integrity, combustion emissions and dynamics of the flame. Copyright © 2012 by ASME.