A Low-Order Modelling of Ducted Flames With Temporally Varying Equivalence Ratio in Realistic Geometries

2011 ◽  
Author(s):  
Owen S. Graham ◽  
Ann P. Dowling
Keyword(s):  
Time Delay ◽  
Heat Release ◽  
Limit Cycle ◽  
Network Model ◽  
Time Domain ◽  
Gas Turbines ◽  
Structural Damage ◽  
Equivalence Ratio ◽  
Function Analysis ◽  
Kinematic Evolution

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.

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

2012 ◽  
Author(s):  
Owen S. Graham ◽  
Ann P. Dowling
Keyword(s):  
Heat Release ◽  
Gas Turbine ◽  
Network Model ◽  
Time Domain ◽  
Acoustic Waves ◽  
Gas Turbines ◽  
Equivalence Ratio ◽  
Fuel Injectors ◽  
Thermoacoustic Instabilities ◽  
Lean Premixed

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.

A Time-Domain Network Model for Nonlinear Thermoacoustic Oscillations

2008 ◽  
Author(s):  
Simon R. Stow ◽  
Ann P. Dowling
Keyword(s):  
Heat Release ◽  
Frequency Domain ◽  
Limit Cycle ◽  
Network Model ◽  
Time Domain ◽  
Acoustic Waves ◽  
Mean Flow ◽  
The Mean ◽  
The Time Domain ◽  
Thermoacoustic Oscillations

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.

A Time-Domain Network Model for Nonlinear Thermoacoustic Oscillations

2009 ◽  
Vol 131 (3) ◽  
Author(s):  
Simon R. Stow ◽  
Ann P. Dowling
Keyword(s):  
Heat Release ◽  
Frequency Domain ◽  
Limit Cycle ◽  
Network Model ◽  
Time Domain ◽  
Acoustic Waves ◽  
Mean Flow ◽  
The Mean ◽  
The Time Domain ◽  
Thermoacoustic Oscillations

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.

Nonlinear Heat-Release/Acoustic Model for Thermoacoustic Instability in Lean Premixed Combustors

1999 ◽  
Vol 121 (3) ◽  
pp. 415-421 ◽  
Author(s):  
A. A. Peracchio ◽  
W. M. Proscia
Keyword(s):  
Heat Release ◽  
Limit Cycle ◽  
Gas Turbines ◽  
Pressure Fluctuations ◽  
Parametric Model ◽  
Operating Conditions ◽  
Cycle Frequency ◽  
Thermoacoustic Instability ◽  
Lean Premixed ◽  
Industrial Gas Turbines

Lean premixed combustors, such as those used in industrial gas turbines to achieve low emissions, are often susceptible to the thermoacoustic combustion instabilities, which manifest themselves as pressure and heat release oscillations in the combustor. These oscillations can result in increased noise and decreased durability due to vibration and flame motion. A physically based nonlinear parametric model has been developed that captures this instability. It describes the coupling of combustor acoustics with the rate of heat release. The model represents this coupling by accounting for the effect of acoustic pressure fluctuations on the varying fuel/air ratio being delivered to the flame, causing a fluctuating heat release due to both fuel air ratio variations and flame front oscillations. If the phasing of the fluctuating heat release and pressure are proper, an instability results that grows into a limit cycle. The nonlinear nature of the model predicts the onset of the instability and additionally captures the resulting limit cycle. Tests of a lean premixed nozzle run at engine scale and engine operating conditions in the UTRC single nozzle rig, conducted under DARPA contract, exhibited instabilities. Parameters from the model were adjusted so that analytical results were consistent with relevant experimental data from this test. The parametric model captures the limit cycle behavior over a range of mean fuel air ratios, showing the instability amplitude (pressure and heat release) to increase and limit cycle frequency to decrease as mean fuel air ratio is reduced.

Adjoint Methods for Elimination of Thermoacoustic Oscillations in a Model Annular Combustor via Small Geometry Modifications

2018 ◽  
Author(s):  
José G. Aguilar ◽  
Matthew P. Juniper
Keyword(s):  
Sensitivity Analysis ◽  
Time Delay ◽  
Heat Release ◽  
Combustion Chamber ◽  
Heat Release Rate ◽  
Release Rate ◽  
Gas Turbines ◽  
Acoustic Pressure ◽  
Annular Combustor ◽  
Thermoacoustic Oscillations

In gas turbines, thermoacoustic oscillations grow if moments of high fluctuating heat release rate coincide with moments of high acoustic pressure. The phase between the heat release rate and the acoustic pressure depends strongly on the flame behaviour (specifically the time delay) and on the acoustic period. This makes the growth rate of thermoacoustic oscillations exceedingly sensitive to small changes in the acoustic boundary conditions, geometry changes, and the flame time delay. In this paper, adjoint-based sensitivity analysis is applied to a thermoacoustic network model of an annular combustor. This reveals how each eigenvalue is affected by every parameter of the system. This information is combined with an optimization algorithm in order to stabilize all thermoacoustic modes of the combustor by making only small changes to the geometry. The final configuration has a larger plenum area, a smaller premix duct area and a larger combustion chamber volume. All changes are less than 6% of the original values. The technique is readily scalable to more complex models and geometries and the inclusion of further constraints, such that the combustion chamber itself should not change. This demonstrates why adjoint-based sensitivity analysis and optimization could become an indispensible tool for the design of thermoacoustically-stable combustors.

A Weakly Nonlinear Approach Based on a Distributed Flame Describing Function to Study the Combustion Dynamics of a Full-Scale Lean-Premixed Swirled Burner

2017 ◽  
Vol 139 (9) ◽  
Author(s):  
Davide Laera ◽  
Sergio M. Camporeale
Keyword(s):  
Heat Release ◽  
Limit Cycle ◽  
Heat Release Rate ◽  
Release Rate ◽  
Gas Turbines ◽  
Full Scale ◽  
Describing Function ◽  
Weakly Nonlinear ◽  
Lean Premixed ◽  
Flame Describing Function

Modern combustion chambers of gas turbines for power generation and aero-engines suffer of thermo-acoustic combustion instabilities generated by the coupling of heat release rate fluctuations with pressure oscillations. The present article reports a numerical analysis of limit cycles arising in a longitudinal combustor. This corresponds to experiments carried out on the longitudinal rig for instability analysis (LRIA) test facility equipped with a full-scale lean-premixed burner. Heat release rate fluctuations are modeled considering a distributed flame describing function (DFDF), since the flame under analysis is not compact with respect to the wavelengths of the unstable modes recorded experimentally. For each point of the flame, a saturation model is assumed for the gain and the phase of the DFDF with increasing amplitude of velocity fluctuations. A weakly nonlinear stability analysis is performed by combining the DFDF with a Helmholtz solver to determine the limit cycle condition. The numerical approach is used to study two configurations of the rig characterized by different lengths of the combustion chamber. In each configuration, a good match has been found between numerical predictions and experiments in terms of frequency and wave shape of the unstable mode. Time-resolved pressure fluctuations in the system plenum and chamber are reconstructed and compared with measurements. A suitable estimate of the limit cycle oscillation is found.

Dynamic and Flashback Characteristics of the Syngas Premixed Swirling Combustors

2008 ◽  
Author(s):  
Quanbin Song ◽  
Aibin Fang ◽  
Gang Xu ◽  
Yanji Xu ◽  
Weiguang Huang
Keyword(s):  
Heat Release ◽  
Gas Turbines ◽  
Equivalence Ratio ◽  
Dynamic Pressure ◽  
Experimental System ◽  
Frequency Characteristics ◽  
Chemiluminescence Intensity ◽  
Power Spectral ◽  
Density Analysis ◽  
Power Spectral Density Analysis

Flashback and combustion thermoacoustic oscillation are two of the major problems in the lean premixed combustors of gas turbine, especially for hydrogen enriched syngas fuels. This paper introduced the experimental results regarding the characteristics research of thermoacoustic oscillation, dynamic heat release and flashback based on a premixed combustion experimental system. The critical flashback equivalence ratio and three dynamic parameters were tested in four types of combustor using eleven kinds of simulated syngas which have the same LHV. Partial flashback and full flashback are observed. With the increase of the hydrogen concentration, the critical equivalence ratio is decreasing and full flashback occurs more easily. The change of nozzle’s geometry structure has a consequent effect on the flashback characteristics. Flashback occurs more easily when the thermoacoustic oscillation is greater, and the corresponding critical equivalence ratio is smaller. Both flashback status and normal status are relative stable combustion within a certain range of equivalence ratio. As the generation and disappearance of the thermoacoustic oscillation and its amplitude-frequency characteristics somehow could be in some randomicity, the critical equivalence ratio is correspondingly not a constant. The dynamic chemiluminescence intensity which indicates the dynamic heat release has a good corresponding relationship in terms of amplitude and frequency characteristics with the dynamic pressure which indicates the thermoacoustic oscillation. Through power spectral density analysis, dynamic pressure and chemiluminescence intensity signals both reflect the dynamic characteristics of the combustor and can be applied in the research of the combustion instability, including the thermoacoustic oscillation characteristics, fluctuation of the chemical heat release rate and flashback characteristics and so on. The results and methods could give some reference to the development of hydrogen enriched gas turbines for future IGCC and poly-generation systems.

scholarly journals Nonlinear Heat-Release/Acoustic Model for Thermoacoustic Instability in Lean Premixed Combustors

1998 ◽  
Author(s):  
A. A. Peracchio ◽  
W. M. Proscia
Keyword(s):  
Heat Release ◽  
Limit Cycle ◽  
Gas Turbines ◽  
Pressure Fluctuations ◽  
Parametric Model ◽  
Operating Conditions ◽  
Cycle Frequency ◽  
Thermoacoustic Instability ◽  
Lean Premixed ◽  
Industrial Gas Turbines

Lean premixed combustors, such as those used in industrial gas turbines to achieve low emissions, are often susceptible to thermoacoustic combustion instabilities, which manifest themselves as pressure and heat release oscillations in the combustor. These oscillations can result in increased noise and decreased durability due to vibration and flame motion. A physically based nonlinear parametric model has been developed that captures this instability. It describes the coupling of combustor acoustics with the rate of heat release. The model represents this coupling by accounting for the effect of acoustic pressure fluctuations on the varying fuel/air ratio being delivered to the flame, causing a fluctuating heat release due to both fuel air ratio variations and flame front oscillations. If the phasing of the fluctuating heat release and pressure are proper, an instability results that grows into a limit cycle. The nonlinear nature of the model predicts the onset of the instability and additionally captures the resulting limit cycle. Tests of a lean premixed nozzle run at engine scale and engine operating conditions in the UTRC Single Nozzle Rig, conducted under DARPA contract, exhibited instabilities. Parameters from the model were adjusted so that analytical results were consistent with relevant experimental data from this test. The parametric model captures the limit cycle behavior over a range of mean fuel air ratios, showing the instability amplitude (pressure and heat release) to increase and limit cycle frequency to decrease as mean fuel air ratio is reduced.

Time-delay-induced instabilities and Hopf bifurcation analysis in 2-neuron network model with reaction–diffusion term

2018 ◽  
Vol 313 ◽  
pp. 306-315 ◽  
Author(s):  
Swati Tyagi ◽  
Subit K Jain ◽  
Syed Abbas ◽  
Shahlar Meherrem ◽  
Rajendra K Ray
Keyword(s):  
Hopf Bifurcation ◽  
Time Delay ◽  
Network Model ◽  
Bifurcation Analysis ◽  
Reaction Diffusion ◽  
Neuron Network ◽  
Diffusion Term

Optical Measurement of Local and Global Transfer Functions for Equivalence Ratio Fluctuations in a Turbulent Swirl Flame

2013 ◽  
Author(s):  
Bernhard C. Bobusch ◽  
Bernhard Ćosić ◽  
Jonas P. Moeck ◽  
Christian Oliver Paschereit
Keyword(s):  
Transfer Function ◽  
Heat Release ◽  
Heat Release Rate ◽  
Release Rate ◽  
Equivalence Ratio ◽  
Transfer Functions ◽  
Optical Measurement ◽  
Low Frequencies ◽  
Fuel Flow ◽  
Calibration Measurement

Equivalence ratio fluctuations are known to be one of the key factors controlling thermoacoustic stability in lean premixed gas turbine combustors. The mixing and thus the spatio-temporal evolution of these perturbations in the combustor flow is, however, difficult to account for in present low-order modeling approaches. To investigate this mechanism, experiments in an atmospheric combustion test rig are conducted. To assess the importance of equivalence ratio fluctuations in the present case, flame transfer functions for different injection positions are measured. By adding known perturbations in the fuel flow using a solenoid valve, the influence of equivalence ratio oscillations on the heat release rate is investigated. The spatially and temporally resolved equivalence ratio fluctuations in the reaction zone are measured using two optical chemiluminescence signals, captured with an intensified camera. A steady calibration measurement allows for the quantitative assessment of the equivalence ratio fluctuations in the flame. This information is used to obtain a mixing transfer function, which relates fluctuations in the fuel flow to corresponding fluctuations in the equivalence ratio of the flame. The current study focuses on the measurement of the global, spatially integrated, transfer function for equivalence ratio fluctuations and the corresponding modeling. In addition, the spatially resolved mixing transfer function is shown and discussed. The global mixing transfer function reveals that despite the good spatial mixing quality of the investigated generic burner, the ability to damp temporal fluctuations at low frequencies is rather poor. It is shown that the equivalence ratio fluctuations are the governing heat release rate oscillation response mechanism for this burner in the low-frequency regime. The global transfer function for equivalence ratio fluctuations derived from the measurements is characterized by a pronounced low-pass characteristic, which is in good agreement with the presented convection–diffusion mixing model.

Sign in / Sign up

Export Citation Format

Share Document