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.

1.
Stow
,
S. R.
, and
Dowling
,
A. P.
, 2004, “
Low-Order Modelling of Thermoacoustic Limit Cycles
,” ASME Paper No. GT2004-54245.
2.
Paschereit
,
C. O.
,
Schuermans
,
B.
,
Polifke
,
W.
, and
Mattson
,
O.
, 2002, “
Measurement of Transfer Matrices and Source Terms of Premixed Flames
,”
ASME J. Eng. Gas Turbines Power
0742-4795,
124
, pp.
239
247
.
3.
Krüger
,
U.
,
Hürens
,
J.
,
Hoffmann
,
S.
,
Krebs
,
W.
, and
Bohn
,
D.
, 1999, “
Prediction of Thermoacoustic Instabilities With Focus on the Dynamic Flame Behavior for the 3A-Series Gas Turbine of Siemens KWU
,” ASME Paper No. 99-GT-111.
4.
Krüger
,
U.
,
Hürens
,
J.
,
Hoffmann
,
S.
,
Krebs
,
W.
,
Flohr
,
P.
, and
Bohn
,
D.
, 2000, “
Prediction and Measurement of Thermoacoustic Improvements in Gas Turbines With Annular Combustion Systems
,” ASME Paper No. 2000-GT-0095.
5.
Hsiao
,
G. C.
,
Pandalai
,
R. P.
,
Hura
,
H. S.
, and
Mongia
,
H. C.
, 1998, “
Combustion Dynamic Modeling for Gas Turbine Engines
,” AIAA Paper No. 98-3380.
6.
Hsiao
,
G. C.
,
Pandalai
,
R. P.
,
Hura
,
H. S.
, and
Mongia
,
H. C.
, 1998, “
Investigation of Combustion Dynamics in Dry-Low-Emission (DLE) Gas Turbine Engines
,” AIAA Paper No. 98-3381.
7.
Dowling
,
A. P.
, and
Stow
,
S. R.
, 2003, “
Acoustic Analysis of Gas Turbine Combustors
,”
J. Propul. Power
0748-4658,
19
(
5
), pp.
751
764
.
8.
Dowling
,
A. P.
, and
Stow
,
S. R.
, 2005, “
Acoustic Analysis of Gas Turbine Combustors
,”
Combustion Instabilities in Gas Turbine Engines: Operational Experience, Fundamental Mechanisms, and Modeling
(
Progress in Astronautics and Aeronautics
Vol.
210
),
T. C.
Lieuwen
and
V.
Yang
, eds.,
AIAA
,
New York
, Chap. 13, pp.
369
414
.
9.
Culick
,
F. E. C.
, 1971, “
Non-Linear Growth and Limiting Amplitude of Acoustic Oscillations in Combustion Chambers
,”
Combust. Sci. Technol.
0010-2202,
3
(
1
), pp.
1
16
.
10.
Culick
,
F. E. C.
, 1989, “
Combustion Instabilities in Liquid-Fueled Propulsion Systems: An Overview
,”
AGARD Conference Proceedings No. 450
, Bath, UK, Oct. 6–7, Paper No. 1,
AGARD
,
Neuilly-sur-Seine, France
.
11.
Peracchio
,
A. A.
, and
Proscia
,
W. M.
, 1999, “
Nonlinear Heat-Release/Acoustic Model for Thermoacoustic Instability in Lean Premixed Combustors
,”
ASME J. Eng. Gas Turbines Power
,
121
(
3
), pp.
415
421
. 0742-4795
12.
Pankiewitz
,
C.
, and
Sattelmayer
,
T.
, 2003, “
Time Domain Simulation of Combustion Instabilities in Annular Combustors
,”
ASME J. Eng. Gas Turbines Power
0742-4795,
125
(
3
), pp.
677
685
.
13.
Schuermans
,
B.
,
Bellucci
,
V.
, and
Paschereit
,
C. O.
, 2003, “
Thermoacoustic Modeling and Control of Multi Burner Combustion Systems
,” ASME Paper No. GT2003-38688.
14.
Bellucci
,
V.
,
Schuermans
,
B.
,
Nowak
,
D.
,
Flohr
,
P.
, and
Paschereit
,
C. O.
, 2004, “
Thermoacoustic Modeling of a Gas Turbine Combustor Equipped With Acoustic Dampers
,” ASME Paper No. GT2004-53977.
15.
Richards
,
G. A.
, and
Janus
,
M. C.
, 1998, “
Characterization of Oscillations During Premix Gas Turbine Combustion
,”
ASME J. Eng. Gas Turbines Power
0742-4795,
120
(
2
), pp.
294
302
.
16.
Lieuwen
,
T.
, and
Zinn
,
B.
, 1998, “
The Role of Equivalence Ratio Oscillations in Driving Combustion Instabilities in low NOx Gas Turbines
,”
Proceedings of the 27th International Symposium on Combustion
, Boulder, CO, Aug. 2–7,
The Combustion Institute
,
Pittsburgh, PA
.
17.
Stow
,
S. R.
, and
Dowling
,
A. P.
, 2003, “
Modelling of Circumferential Modal Coupling Due to Helmholtz Resonators
,” ASME Paper No. GT2003-38168.
18.
Stow
,
S. R.
, and
Dowling
,
A. P.
, 2001, “
Thermoacoustic Oscillations in an Annular Combustor
,” ASME Paper No. 2001-GT-0037.
19.
Cheung
,
W. S.
,
Sims
,
G. J. M.
,
Copplestone
,
R. W.
,
Tilston
,
J. R.
,
Wilson
,
C. W.
,
Stow
S. R.
, and
Dowling
,
A. P.
, 2003, “
Measurement and Analysis of Flame Transfer Function in a Sector Combustor Under High Pressure Conditions
,” ASME Paper No. GT2003-38219.
20.
Armitage
,
C. A.
,
Riley
,
A. J.
,
Cant
,
R. S.
,
Dowling
,
A. P.
, and
Stow
,
S. R.
, 2004, “
Flame Transfer Function for Swirled LPP Combustion From Experiment and CFD
,” ASME Paper No. GT2004-52830.
21.
Schuermans
,
B.
, 2003, “
Modeling and Control of Thermoacoustic Instabilities
,” Ph.D. thesis, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
22.
Schuermans
,
B.
,
Paschereit
,
C. O.
, and
Monkewitz
,
P.
, 2006, “
Non-Linear Combustion Instabilities in Annular Gas-Turbine Combustors
,” AIAA Paper No. AIAA-2006-0549.