Low-Order Modelling of the Response of Ducted Flames in Annular Geometries

Author(s):  
Owen S. Graham ◽  
Ann P. Dowling

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.

Author(s):  
Owen S. Graham ◽  
Ann P. Dowling

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the so-called G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies.


Author(s):  
W. S. Cheung ◽  
G. J. M. Sims ◽  
R. W. Copplestone ◽  
J. R. Tilston ◽  
C. W. Wilson ◽  
...  

Lean premixed prevaporised (LPP) combustion can reduce NOx emissions from gas turbines, but often leads to combustion instability. A flame transfer function describes the change in the rate of heat release in response to perturbations in the inlet flow as a function of frequency. It is a quantitative assessment of the susceptibility of combustion to disturbances. The resulting fluctuations will in turn generate more acoustic waves and in some situations self-sustained oscillations can result. Flame transfer functions for LPP combustion are poorly understood at present but are crucial for predicting combustion oscillations. This paper describes an experiment designed to measure the flame transfer function of a simple combustor incorporating realistic components. Tests were conducted initially on this combustor at atmospheric pressure (1.2 bar and 550 K) to make an early demonstration of the combustion system. The test rig consisted of a plenum chamber with an inline siren, followed by a single LPP premixer/duct and a combustion chamber with a silencer to prevent natural instabilities. The siren was used to induce variable frequency pressure/acoustic signals into the air approaching the combustor. Both unsteady pressure and heat release measurements were undertaken. There was good coherence between the pressure and heat release signals. At each test frequency, two unsteady pressure measurements in the plenum were used to calculate the acoustic waves in this chamber and hence estimate the mass-flow perturbation at the fuel injection point inside the LPP duct. The flame transfer function relating the heat release perturbation to this mass flow was found as a function of frequency. The same combustor hardware and associated instrumentation were then used for the high pressure (15 bar and 800 K) tests. Flame transfer function measurements were taken at three combustion conditions that simulated the staging point conditions (Idle, Approach and Take-off) of a large turbofan gas turbine. There was good coherence between pressure and heat release signals at Idle, indicating a close relationship between acoustic and heat release processes. Problems were encountered at high frequencies for the Approach and Take-off conditions, but the flame transfer function for the Idle case had very good qualitative agreement with the atmospheric-pressure tests. The flame transfer functions calculated here could be used directly for predicting combustion oscillations in gas turbine using the same LPP duct at the same operating conditions. More importantly they can guide work to produce a general analytical model.


Author(s):  
Fangyan Li ◽  
Xiaotao Tian ◽  
Ming-long Du ◽  
Lei Shi ◽  
Jiashan Cui

Abstract Thermoacoustic instabilities are commonly encountered in the development of aeroengines and rocket motors. Research on the fundamental mechanism of thermoacoustic instabilities is beneficial for the optimal design of these engine systems. In the present study, a thermoacoustic instability model based on the lean premixed gas turbines (LPGT) combustion system was established. The longitudinal distribution of heat release caused by the intrinsic instability of flame front is considered in this model. Effects of different heat release distributions and characteristics parameters of the premixed gas (Lewis number Le, Zeldovich Number and Prandtl number Pr) on thermoacoustic instability behaviors of the LPGT system are investigated based on this model. Results show that the LPGT system features with two kinds of unstable thermoacoustic modes. The first one corresponds to the natural acoustic mode of the plenum and the second one corresponds to that of the combustion chamber. The characteristic parameters of premixed gases have a large impact on the stability of the system and even can change the system from stable to unstable state.


Author(s):  
Jialin Su ◽  
Andrew Garmory ◽  
Jon Carrotte

Modern low emission combustion systems are more prone to combustion instabilities due to operation at lean conditions. The response of the airflow passing through the injector to incident acoustic waves is therefore of interest. Airflow fluctuations can initiate, for example, perturbations in stoichiometry and velocity that are subsequently delivered into the heat release region. In the case of liquid fuelled gas turbines the atomisation process will also be affected. Such effects can lead to further unsteady heat release and the generation of acoustic waves, thereby leading to combustion instability. This paper describes experimental measurements and the development of a numerical methodology by which the unsteady airflow response of complex, modern, low emission fuel injectors can be characterised. Single and two passage injector configurations have been investigated which broadly capture many of the features associated with modern fuel injectors. Although targeted at low emission (lean burn) liquid fuelled injector geometries, the methodology developed is thought applicable to a wide range of injector configurations. Initially experimental measurements were used to characterise the overall acoustic impedance of each injector design over a range of frequencies. Such information is also required for the low order thermo-acoustic network models, as typically used in the design process, to predict the stability of the combustion system. In addition to the experimental measurements a methodology was developed using unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations in which acoustic boundary conditions were implemented to reproduce the experimental scenarios. Interrogation of the pressure field enabled similar data analysis techniques to be applied to the numerical data for determining the injector acoustic characteristics. Fidelity of the numerical simulations is confirmed by the excellent agreement between the experimental data and numerical simulations. Furthermore, the unsteady flow field within the passages is difficult to access experimentally, but can be examined in more detail from the simulation results. In this way an improved understanding of the passage flows and their individual responses to the incident acoustic pressure waves can be obtained. The numerical approach is aimed at providing a computationally efficient and economic tool for predicting the acoustic characteristics of the complex geometries typical of modern fuel injector designs. Using this tool injector designs with different acoustic response characteristics can be developed relatively quickly.


Author(s):  
Simon R. Stow ◽  
Ann P. Dowling

Lean premixed prevaporised (LPP) combustion can reduce NOx emissions from gas turbines, but often leads to combustion instability. Acoustic waves produce fluctuations in heat release, for instance by perturbing the fuel-air ratio. These heat fluctuations will in turn generate more acoustic waves and in some situations linear oscillations grow into large amplitude self-sustained oscillations. The resulting limit cycles can cause structural damage. Thermoacoustic oscillations will have a low amplitude initially. Thus linear models can describe the initial growth and hence give stability predictions. An unstable linear mode will grow in amplitude until nonlinear effects become sufficiently important to achieve a limit cycle. While the frequency of the linear mode can often provide a good approximation to that of the resulting limit cycle, linear theories give no prediction of its resulting amplitude. In previous work, we developed a low-order frequency-domain method to model thermoacoustic limit cycles in LPP combustors. This was based on a ‘describing function’ approach and is only applicable when there is a dominant mode and the main nonlinearity is in the combustion response to flow perturbations. In this paper that method is extended into the time domain. The main advantage of the time-domain approach is that limit-cycle stability, the influence of harmonics, and the interaction between different modes can be simulated. In LPP combustion, fluctuations in the inlet fuel-air ratio have been shown to be the dominant cause of unsteady combustion: these occur because velocity perturbations in the premix ducts cause a time-varying fuel-air ratio, which then convects downstream. If the velocity perturbation becomes comparable to the mean flow, there will be an amplitude-dependent effect on the equivalence ratio fluctuations entering the combustor and hence on the rate of heat release. Since the Mach number is low, the velocity perturbation can be comparable to the mean flow, with even reverse flow occurring, while the disturbances are still acoustically linear in that the pressure perturbation is still much smaller than the mean. Hence while the combustion response to flow velocity and equivalence ratio fluctuations must be modelled nonlinearly, the flow perturbations generated as a result of the unsteady combustion can be treated as linear. In developing a time-domain network model for nonlinear thermoacoustic oscillations an initial frequency-domain calculation is performed. The linear network model, LOTAN, is used to categorise the combustor geometry by finding the transfer function for the response of flow perturbations (at the fuel injectors, say) to heat-release oscillations. This transfer function is then converted into the time domain through an inverse Fourier transform to obtain the Green’s function, which thus relates unsteady flow to heat release at previous times. By combining this with a nonlinear flame model (relating heat release to unsteady flow at previous times) a complete time-domain solution can be found by stepping forward in time. If an unstable mode is present, its amplitude will initially grow exponentially (in accordance with linear theory) until saturation effects in the flame model become significant, and eventually a stable limit cycle will be attained. The time-domain approach enables determination of the limit-cycle. In addition, the influence of harmonics and the interaction and exchange of energy between different modes can be simulated. These effects are investigated for longitudinal and circumferential instabilities in an example combustor system and results are compared to frequency-domain limit-cycle predictions.


Author(s):  
Hyung Ju Lee ◽  
Kyu Tae Kim ◽  
Jong Guen Lee ◽  
Bryan D. Quay ◽  
Domenic A. Santavicca

An experimental study was conducted to estimate and confirm equivalence ratio fluctuations at the inlet of a lean premixed gas turbine combustor. Fuel injectors were placed at several locations in the mixing section of the combustor, in order to produce different instability characteristics due to the equivalence ratio fluctuations. An IR absorption technique was used to measure the equivalence ratio fluctuations at the inlet of the dump combustor. The measured IR signals were processed in two different ways and the results were compared to confirm the two calibrated equivalence ratio signals. The processed data showed that the two processing methods gave very similar results, and the phase of the measured equivalence ratio fluctuations at the combustor inlet by the IR absorption technique agreed well with that of equivalence ratio fluctuations predicted by time lags in the mixing section. It was, however, not possible to accurately predict the magnitude of the equivalence ratio fluctuations at the combustor inlet by the time lag analysis because the equivalence ratio fluctuations generataed at the fuel injection location is changed by mixing and diffusion as the fuel is convected through the combustor.


Author(s):  
Simon R. Stow ◽  
Ann P. Dowling

Lean premixed prevaporized (LPP) combustion can reduce NOx emissions from gas turbines but often leads to combustion instability. Acoustic waves produce fluctuations in heat release, for instance, by perturbing the fuel-air ratio. These heat fluctuations will in turn generate more acoustic waves and in some situations linear oscillations grow into large-amplitude self-sustained oscillations. The resulting limit cycles can cause structural damage. Thermoacoustic oscillations will have a low amplitude initially. Thus linear models can describe the initial growth and hence give stability predictions. An unstable linear mode will grow in amplitude until nonlinear effects become sufficiently important to achieve a limit cycle. While the frequency of the linear mode can often provide a good approximation to that of the resulting limit cycle, linear theories give no prediction of its resulting amplitude. In previous work, we developed a low-order frequency-domain method to model thermoacoustic limit cycles in LPP combustors. This was based on a “describing-function” approach and is only applicable when there is a dominant mode and the main nonlinearity is in the combustion response to flow perturbations. In this paper that method is extended into the time domain. The main advantage of the time-domain approach is that limit-cycle stability, the influence of harmonics, and the interaction between different modes can be simulated. In LPP combustion, fluctuations in the inlet fuel-air ratio have been shown to be the dominant cause of unsteady combustion: These occur because velocity perturbations in the premix ducts cause a time-varying fuel-air ratio, which then convects downstream. If the velocity perturbation becomes comparable to the mean flow, there will be an amplitude-dependent effect on the equivalence ratio fluctuations entering the combustor and hence on the rate of heat release. Since the Mach number is low, the velocity perturbation can be comparable to the mean flow, with even reverse flow occurring, while the disturbances are still acoustically linear in that the pressure perturbation is still much smaller than the mean. Hence while the combustion response to flow velocity and equivalence ratio fluctuations must be modeled nonlinearly, the flow perturbations generated as a result of the unsteady combustion can be treated as linear. In developing a time-domain network model for nonlinear thermoacoustic oscillations an initial frequency-domain calculation is performed. The linear network model, LOTAN, is used to categorize the combustor geometry by finding the transfer function for the response of flow perturbations (at the fuel injectors, say) to heat-release oscillations. This transfer function is then converted into the time domain through an inverse Fourier transform to obtain Green’s function, which thus relates unsteady flow to heat release at previous times. By combining this with a nonlinear flame model (relating heat release to unsteady flow at previous times) a complete time-domain solution can be found by stepping forward in time. If an unstable mode is present, its amplitude will initially grow exponentially (in accordance with linear theory) until saturation effects in the flame model become significant, and eventually a stable limit cycle will be attained. The time-domain approach enables determination of the limit cycle. In addition, the influence of harmonics and the interaction and exchange of energy between different modes can be simulated. These effects are investigated for longitudinal and circumferential instabilities in an example combustor system and the results are compared with frequency-domain limit-cycle predictions.


Author(s):  
A. Andreini ◽  
B. Facchini ◽  
L. Mangani ◽  
F. Simonetti

In the last years, the more restrictive environmental legislations have constrained gas turbine manufacturers to the development of new low-emission combustors. Lean Premixed technology has become a necessary standard to meet emissions requirements and allowing an heavy reduction of nitrogen oxides emission. This kind of technology, due to the use of lean premixed mixtures, is severely affected by thermoacoustic phenomena which cause damages to combustor components and consequently reduce the overall gas turbine life of a factor of two or more. Specifically, premixed flames pose the threat of pressure oscillations. This phenomenon is the effect of the strong interaction between combustion heat-release and fluid dynamics aspects. In order to investigate thermoacoustic instabilities, a mono-dimensional code was developed and validated. It takes into account only longitudinal frequencies and it is thought to be highly modular to modify or add blocks, corresponding to different thermoacoustic models. The tool is based on a lumped-parameter approach, which consists in considering constant mean flow quantities over each fundamental straight duct element and a nodal point at each duct interface. For each interface, where an acoustic impedance could be present, the linearized fluctuating mass, momentum and energy equations are solved including entropic waves. To validate such tool, several tests, referring to actual test rigs and experimental gas turbine combustor geometries, were performed. The results show a general agreement with empirical data and other numerical results reported in literature in terms of resonance frequencies, stability and modal shapes, both for no flame and fluctuating heat release cases.


Author(s):  
Wyatt Culler ◽  
Janith Samarasinghe ◽  
Bryan D. Quay ◽  
Domenic A. Santavicca ◽  
Jacqueline O’Connor

Combustion instability in gas turbines can be mitigated using active techniques or passive techniques, but passive techniques are almost exclusively used in industrial settings. While fuel staging, a common passive technique, is effective in reducing the amplitude of self-excited instabilities in gas turbine combustors at steady-state conditions, the effect of transients in fuel staging on self-excited instabilities is not well understood. This paper examines the effect of fuel staging transients on a laboratory-scale five-nozzle can combustor undergoing self-excited instabilities. The five nozzles are arranged in a four-around-one configuration and fuel staging is accomplished by increasing the center nozzle equivalence ratio. When the global equivalence ratio is φ = 0.70 and all nozzles are fueled equally, the combustor undergoes self-excited oscillations. These oscillations are suppressed when the center nozzle equivalence ratio is increased to φ = 0.80 or φ = 0.85. Two transient staging schedules are used, resulting in transitions from unstable to stable operation, and vice-versa. It is found that the characteristic instability decay times are dependent on the amount of fuel staging in the center nozzle. It is also found that the decay time constants differ from the growth time constants, indicating hysteresis in stability transition points. High speed CH* chemiluminescence images in combination with dynamic pressure measurements are used to determine the instantaneous phase difference between the heat release rate fluctuation and the combustor pressure fluctuation throughout the combustor. This analysis shows that the instability onset process is different from the instability decay process.


Author(s):  
P. Griebel ◽  
R. Bombach ◽  
A. Inauen ◽  
R. Scha¨ren ◽  
S. Schenker ◽  
...  

The present experimental study focuses on flame characteristics and turbulent flame speeds of lean premixed flames typical for stationary gas turbines. Measurements were performed in a generic combustor at a preheating temperature of 673 K, pressures up to 14.4 bars (absolute), a bulk velocity of 40 m/s, and an equivalence ratio in the range of 0.43–0.56. Turbulence intensities and integral length scales were measured in an isothermal flow field with Particle Image Velocimetry (PIV). The turbulence intensity (u′) and the integral length scale (LT) at the combustor inlet were varied using turbulence grids with different blockage ratios and different hole diameters. The position, shape, and fluctuation of the flame front were characterized by a statistical analysis of Planar Laser Induced Fluorescence images of the OH radical (OH-PLIF). Turbulent flame speeds were calculated and their dependence on operating conditions (p, φ) and turbulence quantities (u′, LT) are discussed and compared to correlations from literature. No influence of pressure on the most probable flame front position or on the turbulent flame speed was observed. As expected, the equivalence ratio had a strong influence on the most probable flame front position, the spatial flame front fluctuation, and the turbulent flame speed. Decreasing the equivalence ratio results in a shift of the flame front position farther downstream due to the lower fuel concentration and the lower adiabatic flame temperature and subsequently lower turbulent flame speed. Flames operated at leaner equivalence ratios show a broader spatial fluctuation as the lean blow-out limit is approached and therefore are more susceptible to flow disturbances. In addition, because of a lower turbulent flame speed these flames stabilize farther downstream in a region with higher velocity fluctuations. This increases the fluctuation of the flame front. Flames with higher turbulence quantities (u′, LT) in the vicinity of the combustor inlet exhibited a shorter length and a higher calculated flame speed. An enhanced turbulent heat and mass transport from the recirculation zone to the flame root location due to an intensified mixing which might increase the preheating temperature or the radical concentration is believed to be the reason for that.


Sign in / Sign up

Export Citation Format

Share Document