Prediction of sound emission from open turbulent premixed flames

The sound emission from open turbulent flames is dictated by the two-point spatial correlation of rate of change of fluctuating heat release rate and this correlation has not been investigated directly in the past studies. Turbulent premixed flame data from DNS and laser diagnostics are analyzed to study this correlation function and the two-point spatial correlation of the fluctuating heat release rate. This shows that the correlation functions have simple Gaussian forms whose integral length scale is related to the laminar flame thickness and amplitude depends on the spatial distribution of the time-mean rate of heat release. These results and RANS-CFD solution of open turbulent premixed flames are post-processed to obtain the far field SPL, which agrees well with measured values. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Analysis and simulation of combustion oscillations in a laboratory spray combustor

The unstable combustion that can occur in combustion chambers is a major problem for aeroengines and ground-based industrial gas turbines. Nowadays, CFD provides a flexible, low cost tool to supplement direct measurement. This paper presents simulations of combustion oscillations in a liquid-fuelled experimental rig at the University of Cambridge. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone and to develop inlet and outlet boundary conditions just upstream and downstream of the combustion region enabling the CFD calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in good agreement with experimental measurements. More details about the unstable combustion can be obtained from the simulation results. The approach developed here is expected to provide a powerful tool for the design and operation of stable combustion systems. Copyright © 2009 by ASME.

Calculations and predications of combustion oscillations

The combustion oscillations are the phenomena which we may meet in developing the clean, safe and efficient power and propulsion systems. This paper summarizes authors' systematic work on fuel spray combustion oscillations in the recent years. Combining CFD calculations and stability analysis, a new approach of predicting combustion stabilities was developed. With this approach, detailed flow information and unstable modes can be obtained by CFD and solving perturbation equations, respectively. The results provide the guidance on understanding combustion instability mechanisms and developing the control strategies.

Turbomachinery computational fluid dynamics: asymptotes and paradigm shifts.

This paper reviews the development of computational fluid dynamics (CFD) specifically for turbomachinery simulations and with a particular focus on application to problems with complex geometry. The review is structured by considering this development as a series of paradigm shifts, followed by asymptotes. The original S1-S2 blade-blade-throughflow model is briefly described, followed by the development of two-dimensional then three-dimensional blade-blade analysis. This in turn evolved from inviscid to viscous analysis and then from steady to unsteady flow simulations. This development trajectory led over a surprisingly small number of years to an accepted approach-a 'CFD orthodoxy'. A very important current area of intense interest and activity in turbomachinery simulation is in accounting for real geometry effects, not just in the secondary air and turbine cooling systems but also associated with the primary path. The requirements here are threefold: capturing and representing these geometries in a computer model; making rapid design changes to these complex geometries; and managing the very large associated computational models on PC clusters. Accordingly, the challenges in the application of the current CFD orthodoxy to complex geometries are described in some detail. The main aim of this paper is to argue that the current CFD orthodoxy is on a new asymptote and is not in fact suited for application to complex geometries and that a paradigm shift must be sought. In particular, the new paradigm must be geometry centric and inherently parallel without serial bottlenecks. The main contribution of this paper is to describe such a potential paradigm shift, inspired by the animation industry, based on a fundamental shift in perspective from explicit to implicit geometry and then illustrate this with a number of applications to turbomachinery.

Interpretation of measured profile losses in unsteady wake-turbine blade interaction studies

The interaction of wakes shed by a moving bladerow with a downstream bladerow causes unsteady flow. The meaning of the freestream stagnation pressure and stagnation enthalpy in these circumstances has been examined using simple analyses, measurements and CFD. The unsteady flow in question arises from the behaviour of the wakes as so-called negative-jets. The interactions of the negative-jets with the downstream blades lead to fluctuations in static pressure which in turn generate fluctuations in the stagnation pressure and stagnation enthalpy. It is shown that the fluctuations of the stagnation quantities created by unsteady effects within the bladerow are far greater than those within the incoming wake. The time-mean exit profiles of the stagnation pressure and stagnation enthalpy are affected by these large fluctuations. This phenomenon of energy separation is much more significant than the distortion of the time-mean exit profiles that is caused directly by the cross-passage transport associated with the negative-jet, as described by Kerrebrock and Mikolajczak. Finally, it is shown that if only time-averaged values of loss are required across a bladerow, it is nevertheless sufficient to determine the time-mean exit stagnation pressure.

Method for determining erosion resistance of fibre-insulating materials

A method was developed for the estimation of the erosive wear of fiber-insulating materials. The wear increases with increasing impact velocity of the particles, increasing impact angle, particle size and the thermal ageing of the fibre elements. Through CFD simulation of the particle-containing gas flow, the erosion depth can be predicted.

Multi-objective design and optimisation of bypass outlet-guide vanes

This paper describes the development of an automated design optimization system that makes use of a high fidelity Reynolds-Averaged CFD analysis procedure to minimize the fan forcing and fan BOGV (bypass outlet guide vane) losses simultaneously taking into the account the down-stream pylon and RDF (radial drive fairing) distortions. The design space consists of the OGV's stagger angle, trailing-edge recambering, axial and circumferential positions leading to a variable pitch optimum design. An advanced optimization system called SOFT (Smart Optimisation for Turbomachinery) was used to integrate a number of pre-processor, simulation and in-house grid generation codes and postprocessor programs. A number of multi-objective, multi-point optimiztion were carried out by SOFT on a cluster of workstations and are reported herein.

Heat transfer and aerodynamics of turbine blade tips in a linear cascade

Local measurements of the heat transfer coefficient and pressure coefficient were conducted on the tip and near tip region of a generic turbine blade in a five-blade linear cascade. Two tip clearance gaps were used: 1.6% and 2.8% chord. Data was obtained at a Reynolds number of 2.3 × 10 5 based on exit velocity and chord. Three different tip geometries were investigated: a flat (plain) tip, a suction-side squealer, and a cavity squealer. The experiments reveal that the flow through the plain gap is dominated by flow separation at the pressure-side edge and that the highest levels of heat transfer are located where the flow reattaches on the tip surface. High heat transfer is also measured at locations where the tip-leakage vortex has impinged onto the suction surface of the aerofoil. The experiments are supported by flow visualisation computed using the CFX CFD code which has provided insight into the fluid dynamics within the gap. The suction-side and cavity squealers are shown to reduce the heat transfer in the gap but high levels of heat transfer are associated with locations of impingement, identified using the flow visualisation and aerodynamic data. Film cooling is introduced on the plain tip at locations near the pressure-side edge within the separated region and a net heat flux reduction analysis is used to quantify the performance of the successful cooling design. copyright © 2005 by ASME.

Evaluation of wave drag reduction by flow control

An analytical expression is proposed to estimate the wave drag of an aerofoil equipped with shock control. The analysis extends the conventional approach for a single normal shock wave, based on the knowledge that all types of successful shock control on transonic aerofoils cause bifurcated λ-shock structures. The influence of surface curvature on the λ-shock structure has been taken into account. The extended method has been found to produce fairly good agreement with the results obtained by CFD methods while requiring negligible computational effort. This new formulation is expected to be beneficial in the industrial design process of transonic aerofoils and wings where a large number of computational simulations have to be performed.

A new alternative for reduction of secondary flows in low pressure turbines

This paper describes a new flow mechanism for the reduction of secondary flows in Low Pressure Turbines using the benefit of contoured endwalls. The extensive application of contoured endwalls in recent years has provided a deeper understanding of the physical phenomenon that governs the reduction of secondary flows. Based on this understanding, the endwall geometry of a linear cascade of solid-thin profiles typical of Low Pressure Turbines has been redesigned. Experimental data are presented for the validation of this new solution. Based on these data, a reduction of 72% in the SKEH and 20% in the mixed-out endwall losses can be obtained. CFD simulations are also presented to illustrate the effect of the new endwall on the secondary flows. Furthermore, an explanation of the flow mechanism that governs the reduction of the SKEH and the losses is given. Copyright © 2006 by ASME.

Compressor wake/leading-edge interactions at off design incidences

In this paper, the effects of wake/leading-edge interactions were studied at off-design conditions. Measurements were performed on the stator-blade suction surface at midspan. The leading-edge flow-field was investigated using hotwire micro-traverses, hotfilm surface shear-stress sensors and pressure micro-tappings. The trailing-edge flow-field was investigated using hotwire boundary-layer traverses. Unsteady CFD calculations were also performed to aid the interpretation of the results. At low flow coefficients, the time-averaged momentum thickness of the leading-edge boundary layer was found to rise as the flow coefficient was reduced. The time-resolved momentum-thickness rose due to the interaction of the incoming rotor wake. As the flow coefficient was reduced, the incoming wakes increased in pitch-wise extent, velocity deficit and turbulence intensity. This increased both the time-resolved rise in the momentum thickness and the turbulent spot production within the wake affected boundary-layer. Close to stall, a drop in the leading-edge momentum thickness was observed in-between wake events. This was associated with the formation of a leading-edge separation bubble in-between wake events. The wake interaction with the bubble gave rise to a shedding phenomenon, which produced large length scale disturbances in the surface shear stress. Copyright © 2008 by ASME.

Towards topology-free optimization - An application to turbine internal cooling geometries

The application of automated design optimization to real-world, complex geometry problems is a significant challenge - especially if the topology is not known a priori like in turbine internal cooling. The long term goal of our work is to focus on an end-to-end integration of the whole CFD Process, from solid model through meshing, solving and post-processing to enable this type of design optimization to become viable & practical. In recent papers we have reported the integration of a Level Set based geometry kernel with an octree-based cut- Cartesian mesh generator, RANS flow solver, post-processing & geometry editing all within a single piece of software - and all implemented in parallel with commodity PC clusters as the target. The cut-cells which characterize the approach are eliminated by exporting a body-conformal mesh guided by the underpinning Level Set. This paper extends this work still further with a simple scoping study showing how the basic functionality can be scripted & automated and then used as the basis for automated optimization of a generic gas turbine cooling geometry. Copyright © 2008 by W.N.Dawes.

An experimental and numerical study of an oscillating transonic shock wave in a duct

A combined experimental and numerical study of a transonic shock wave in a parallel walled duct subject to downstream pressure perturbations has been conducted. Experiments and simulations have been carried out with a shock strength of M∞ = 1.4 for pressure perturbation frequencies in the range 16-90 Hz. The dynamics of unsteady shock motion and the interaction structure between the unsteady transonic shock wave and the turbulent tunnel floor boundary layer have been investigated. It is found that the (experimentally measured) dynamics of shock motion are generally well predicted by the computational scheme, especially at relatively low (≈ 40 Hz) frequencies. However, at higher frequencies (≈ 90 Hz), some subtle differences between the shock dynamics measured in experiments and those predicted by Computational Fluid Dynamics (CFD) exist. There is evidence from experiments that variations in shock / boundary layer interaction (SBLI) structure caused by shock motion are responsible for a change in the nature of shock dynamics between low and high frequency. In contrast, numerical results at low and high frequencies do not differ significantly and this suggests that the numerical method is not fully capturing the physics of the unsteady flow. Possible reasons for this are considered and a number of areas where CFD is unable to replicate experimental observations are identified. Significantly, CFD predicts changes in SBLI structure due to shock motion that are much too large and this may explain why none of the subtle effects on shock dynamics seen in experiments occur in CFD. Further work developing numerical methods that demonstrate a more realistic sensitivity of SBLI structure to unsteady shock motion is required. Copyright © 2010 by P.J.K. Bruce.

Simulation of solitary wave impact on coastal structures using weakly compressible and incompressible SPH methods

In this paper, we engage a Lagrangian, particle-based CFD method, named Smoothed Particle Hydrodynamic (SPH) to study the solitary wave motion and its impact on coastal structures. Two-dimensional weakly compressible and incompressible SPH models were applied to simulate wave impacting on seawall and schematic coastal house. The results confirmed the accuracy of both models for predicting the wave surface profiles. The incompressible SPH model performed better in predicting the pressure field and impact loadings on coastal structures than the weakly compressible SPH model. The results are in qualitatively agreement with experimental results. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).