Research Papers: Gas Turbines: Combustion, Fuels, and Emissions

A Time-Domain Network Model for Nonlinear Thermoacoustic Oscillations

[+] Author and Article Information
Simon R. Stow1

Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, UKsimon.stow@rolls-royce.com

Ann P. Dowling

Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, UK


Present address: Rolls-Royce plc, Derby DE24 8BJ, UK.

J. Eng. Gas Turbines Power 131(3), 031502 (Feb 03, 2009) (10 pages) doi:10.1115/1.2981178 History: Received March 28, 2008; Revised April 25, 2008; Published February 03, 2009

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.

Copyright © 2009 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 1

Schematics of possible application of thermoacoustic network model (top) and modeling approach (bottom)

Grahic Jump Location
Figure 2

Variation of Q̂(m)/Q̂L with β. The solid, dashed, dashed-dotted, and dotted lines denote m=1, 3, 5, and 7, respectively (the value being zero for even m). Note that Q̂(1)/Q̂L is equivalent to the nonlinear flame transfer function relative to the linearized version, Tflame(ω,A)/TflameL(ω).

Grahic Jump Location
Figure 3

Transfer functions, Tn,n(ω), and Green’s functions, Gn,n(t), for original geometry. The solid and dashed lines denote n=±1 (these being equivalent) and n=0, respectively.

Grahic Jump Location
Figure 4

Time-domain solution for original geometry. m′−1(t)/m¯ and m′1(t)/m¯ are shown; in the third plot, the solid line denotes the former while the dashed lines denotes the latter. (m′0(t) is negligible for this solution.)

Grahic Jump Location
Figure 5

Time-domain solution for original geometry, assuming a single-spinning wave, n=1. m′1(t)/m¯ and Q1′(t)/Q¯ are shown; in the third plot, the solid line denotes the former while the dashed lines denotes the latter.

Grahic Jump Location
Figure 6

m′n(t)/m¯ for time-domain solution for long-combustor geometry (Calculation 1)

Grahic Jump Location
Figure 7

m′n(t)/m¯ for time-domain solution for long-combustor geometry (Calculation 2)

Grahic Jump Location
Figure 8

m′n(t)/m¯ for time-domain solution for long-combustor geometry (Calculation 3)

Grahic Jump Location
Figure 9

m′n(t)/m¯ for time-domain solution for long-combustor geometry (Calculation 4)

Grahic Jump Location
Figure 10

Time-domain solution for original geometry, with N=3 and fmax=2 kHz. m′n(t)/m¯ is shown for n=−3, −1, 1, and 3. (m′n(t) is negligible for even n.)



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In