Acceleration of an unstructured hybrid mesh RANS solver by porting to GPU architectures

The modern CFD process consists of mesh generation, flow solving and post-processing integrated into an automated workflow. During the last several years we have developed and published research aimed at producing a meshing and geometry editing system, implemented in an end-to-end parallel, scalable manner and capable of automatic handling of large scale, real world applications. The particular focus of this paper is the associated unstructured mesh RANS flow solver and the porting of it to GPU architectures. After briefly describing the solver itself, the special issues associated with porting codes using unstructured data structures are discussed - followed by some application examples. Copyright © 2011 by W.N. Dawes.

Passive control of combustion instability in a low emissions aeroderivative gas turbine

This paper describes the conceptual ideas, the theoretical validation, the laboratory testing and the field trials of a recently patented fuel-air mixing device for use in high-pressure ratio, low emissions, gaseous-fueled gas turbines. By making the fuel-air mixing process insensitive to pressure fluctuations in the combustion chamber, it is possible to avoid the common problem of positive feedback between mixture strength and the unsteady combustion process. More specifically, a mixing duct has been designed such that fuel-air ratio fluctuations over a wide range of frequencies can be damped out by passive design means. By scaling the design in such a way that the range of damped frequencies covers the frequency spectrum of the acoustic modes in the combustor, the instability mechanism can be removed. After systematic development, this design philosophy was successfully applied to a 35:1 pressure ratio aeroderivative gas turbine yielding very low noise levels and very competitive NOx and CO measurements. The development of the new premixer is described from conceptual origins through analytic and CFD evaluation to laboratory testing and final field trials. Also included in this paper are comments about the practical issues of mixing, flashback resistance and autoignition.

Robust airfoil design with respect to boundary layer transition

In order to disign an airfoil of which maximum lift coefficient (CL max) is not sensitive to location of forced top boundary layer transition. Taking maximizing mean value of CL max and minimizing standard deviation as biobjective, leading edge radius, manximum thickness and its location, maximum camber and its location as deterministic design variables, location of forced top boundary layer transition as stochastic variable, XFOIL as deterministic CFD solver, non-intrusive polynomial chaos as substitute of Monte Carlo method, we completed a robust airfoil design problem. Results demonstrate performance of initial airfoil is enhanced through reducing standard deviation of CL max. Besides, we also know maximum thickness has the most dominating effect on mean value of CL max, location of maximum thickness has the most dominating effect on standard deviation of CL max, maximum camber has a little effect on both mean value and standard deviation, and maximum camber is the only element of which increase can lead increase of mean value and standard deviation at the same time. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Shallow water equations: Magic and limitations

D Liang from Cambridge University explains the shallow water equations and their applications to the dam-break and other steep-fronted flow modeling. They assume that the horizontal scale of the flow is much greater than the vertical scale, which means the flow is restricted within a thin layer, thus the vertical momentum is insignificant and the pressure distribution is hydrostatic. The left hand sides of the two momentum equations represent the acceleration of the fluid particle in the horizontal plane. If the fluid acceleration is ignored, then the two momentum equations are simplified into the so-called diffusion wave equations. In contrast to the SWEs approach, it is much less convenient to model floods with the Navier-Stokes equations. In conventional computational fluid dynamics (CFD), cumbersome treatments are needed to accurately capture the shape of the free surface. The SWEs are derived using the assumptions of small vertical velocity component, smooth water surface, gradual variation and hydrostatic pressure distribution.

Assessment of computational fluid dynamics and experimental data for shock boundary-layer interactions

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop, numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric, and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and, in general, it was difficult to discern clear trends in the data. For the Reynolds-averaged Navier-Stokes (RANS) methods, the choice of turbulence model appeared to be the largest factor in solution accuracy. Scale-resolving methods, such as large-eddy simulation (LES), hybrid RANS/LES, and direct numerical simulation, produced error levels similar to RANS methods but provided superior predictions of normal stresses. Copyright © 2012 by Daniella E. Raveh and Michael Iovnovich.

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.

Contra-rotating open rotor operation for improved aerodynamics and noise at takeoff

The contra-rotating open rotor is, once again, being considered as an alternative to the advanced turbofan to address the growing pressure to cut aviation fuel consumption and carbon dioxide emissions. One of the key challenges is meeting community noise targets at takeoff. Previous open rotor designs are subject to poor efficiency at takeoff due to the presence of large regions of separated flow on the blades as a result of the high incidence needed to achieve the required thrust. This is a consequence of the fixed rotor rotational speed constraint typical of variable pitch propellers. Within the study described in this paper, an improved operation is proposed to improve performance and reduce rotorrotor interaction noise at takeoff. Three-dimensional computational fluid dynamics (CFD) calculations have been performed on an open rotor rig at a range of takeoff operating conditions. These have been complemented by analytical tone noise predictions to quantify the noise benefits of the approach. The results presented show that for a given thrust, a combination of reduced rotor pitch and increased rotor rotational speed can be used to reduce the incidence onto the front rotor blades. This is shown to eliminate regions of flow separation, reduce the front rotor tip loss and reduce the downstream stream tube contraction. The wakes from the front rotor are also made wider with lower velocity defect, which is found to lead to reduced interaction tone noise. Unfortunately, the necessary increase in blade speed leads to higher relative Mach numbers, which can increase rotor alone noise. In summary, the combined CFD and aero-acoustic analysis in this paper shows how careful operation of an open rotor at takeoff, with moderate levels of re-pitch and speed increase, can lead to improved front rotor efficiency as well as appreciably lower overall noise across all directivities. Copyright © 2011 by ASME.

The influence of turbine stator well coolant flow rate and passage configuration on cooling effectiveness

Market competitiveness for aero engine power plant dictates that improvements in engine performance and reliability are guaranteed a priori by manufacturers. The requirement to accurately predict the life of engine components makes exacting demands of the internal air system, which must provide effective cooling over the engine duty cycle with the minimum consumption of compressor section air. Tests have been conducted at the University of Sussex using a turbine test facility which comprises a two stage turbine with an individual stage pressure ratio of 1.7:1. Main annulus air is supplied by an adapted Rolls-Royce Dart compressor at up to 440 K and 4.8 kg s-1. Cooling flow rates ranging from 0.71 to 1.46 Cw, ent, a disc entrainment parameter, have been used to allow ingress or egress dominated stator well flow conditions. The mechanical design of the test section allows internal cooling geometry to be rapidly re-configured, allowing the effect of jet momentum and coolant trajectory to be investigated. An important facet to this investigation is the use of CFD to model and analyse the flow structures associated with the cavity conditions tested, as well as to inform the design of cooling path geometry. This paper reports on the effectiveness of stator well coolant flow rate and delivery configurations using experimental data and also CFD analysis to better quantify the effect of stator well flow distribution on component temperatures. Copyright © 2011 by Rolls-Royce plc.

Experimental investigation of turbine stator well rim seal, reingestion and interstage seal flows using gas concentration techniques and displacement measurements

Gas turbine engine performance requires effective and reliable internal cooling over the duty cycle of the engine. Life predictions for rotating components subject to the main gas path temperatures are vital. This demands increased precision in the specification of the internal air system flows which provide turbine stator well cooling and sealing. This in turn requires detailed knowledge of the flow rates through rim seals and interstage labyrinth seals. Knowledge of seal movement and clearances at operating temperatures is of great importance when prescribing these flows. A test facility has been developed at the University of Sussex, incorporating a two stage turbine rated at 400 kW with an individual stage pressure ratio of 1.7:1. The mechanical design of the test facility allows internal cooling geometry to be rapidly re-configured, while cooling flow rates of between 0.71 CW, ENT and 1.46 C W, ENT, may be set to allow ingress or egress dominated cavity flows. The main annulus and cavity conditions correspond to in cavity rotational Reynolds numbers of 1.71×106< Reφ <1.93×106. Displacement sensors have been used to establish hot running seal clearances over a range of stator well flow conditions, allowing realistic flow rates to be calculated. Additionally, gas seeding techniques have been developed, where stator well and main annulus flow interactions are evaluated by measuring changes in gas concentration. Experiments have been performed which allow rim seal and re-ingestion flows to be quantified. It will be shown that this work develops the measurement of stator well cooling flows and provides data suitable for the validation of improved thermo-mechanical and CFD codes, beneficial to the engine design process. Copyright © 2011 by Rolls-Royce plc.

Shock control bump robustness enhancement

Robustness enhancement for Shock Control Bumps (SCBs) on transonic wings is an ongoing topic because most designs provide drag savings only in a relatively small band of the airfoil polar. In this paper, different bump shapes are examined with CFD methods which are validated first by comparison with wind tunnel results. An evaluation method is introduced allowing the robustness assessment of a certain design with little computational effort. Shape optimizations are performed to trim SCB designs to maximum performance on the one hand and maximum robustness on the other hand. The results are analysed and different and parameters influencing the robustness are suggested. Copyright © 2012 by Klemens Nuebler.

The design of multi-element airfoils through multi-objective optimization techniques

This paper presents the development and the application of a multi-objective optimization framework for the design of two-dimensional multi-element high-lift airfoils. An innovative and efficient optimization algorithm, namely Multi-Objective Tabu Search (MOTS), has been selected as core of the framework. The flow-field around the multi-element configuration is simulated using the commercial computational fluid dynamics (cfd) suite Ansys cfx. Elements shape and deployment settings have been considered as design variables in the optimization of the Garteur A310 airfoil, as presented here. A validation and verification process of the cfd simulation for the Garteur airfoil is performed using available wind tunnel data. Two design examples are presented in this study: a single-point optimization aiming at concurrently increasing the lift and drag performance of the test case at a fixed angle of attack and a multi-point optimization. The latter aims at introducing operational robustness and off-design performance into the design process. Finally, the performance of the MOTS algorithm is assessed by comparison with the leading NSGA-II (Non-dominated Sorting Genetic Algorithm) optimization strategy. An equivalent framework developed by the authors within the industrial sponsor environment is used for the comparison. To eliminate cfd solver dependencies three optimum solutions from the Pareto optimal set have been cross-validated. As a result of this study MOTS has been demonstrated to be an efficient and effective algorithm for aerodynamic optimizations. Copyright © 2012 Tech Science Press.

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.

A novel MPT noise prediction methodology for highly-integrated propulsion systems with inlet flow distortion

Embedded propulsion systems, such as for example used in advanced hybrid-wing body aircraft, can potentially offer major fuel burn and noise reduction benefits but introduce challenges in the aerodynamic and acoustic integration of the high-bypass ratio fan system. A novel approach is proposed to quantify the effects of non-uniform flow on the generation and propagation of multiple pure tone noise (MPTs). The new method is validated on a conventional inlet geometry first. The ultimate goal is to conduct a parametric study of S-duct inlets in order to quantify the effects of inlet design parameters on the acoustic signature. The key challenge is that the mechanism underlying the 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. 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 MPT noise generation mechanisms while greatly reducing computational cost. A single, 3-D full-wheel unsteady 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 mean flow. Several numerical tools were developed to enable the implementation of this new approach. Parametric studies were conducted to determine appropriate grid and time step sizes for the propagation of acoustic waves. The Ffowcs-Williams and Hawkings integral method is used to propagate the noise to far field receivers. Non-reflecting boundary conditions are implemented through the use of acoustic buffer zones. The body force modeling approach is validated and proof-of-concept studies demonstrate the generation of disturbances at both blade-passing and shaft-order frequencies using the perturbed body force method. The full methodology is currently being validated using NASA's Source Diagnostic Test (SDT) fan and inlet geometry. Copyright © 2009 by Jeff Defoe, Alex Narkaj & Zoltan Spakovszky.

High-order DES and zonal LES of a labyrinth seal

Hybrid numerical large eddy simulation (NLES) and detached eddy simulation (DES) methods are assessed on a labyrinth seal geometry. A high sixth order discretization scheme is used and is validated using a test case of a two dimensional vortex. The hybrid approach adopts a new blending function and along with DES is initially validated using a simple cavity flow. The NLES method is also validated outside of RANS zones. It is found that there is very little resolved turbulence in the cavity for the DES simulation. For the labyrinth seal calculations the DES approach is problematic giving virtually no resolved turbulence content. It is seen that over the tooth tips the extent of the LES region is small and is likely to be a strong contributor to excessive flow damping in these regions. On the other hand the zonal Hamilton-Jacobi approach did not suffer from this trait. In both cases the meshes used are considered to be hybrid RANS-LES adequate. Fortunately (or perhaps unfortunately) the DES profiles are in agreement with the time mean experimental measurements. It is concluded that for an inexperienced CFD practitioner this could have wider implications particularly if transient results such as unsteady loading are desired. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

The Effect of Inlet Guide Vanes on Inlet Flow Distortion Transfer and Transonic Fan Stability

An investigation was carried out into the effects of variable inlet guide vanes (VIGVs) on the performance and stability margin of a transonic fan in the presence of inlet flow distortion. The study was carried out using computational fluid dynamics (CFD) and validated with experimental data. The capability of CFD to predict the changes in performance with or without VIGVs in the presence of an inlet flow distortion is assessed. Results show that the VIGVs improve the performance and stability margin and do so by reducing the amount of swirl at inlet to the rotor component of the fan.