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

# Nonlinear Phenomena in Thermoacoustic Systems With Premixed Flames

[+] Author and Article Information
Karthik Kashinath

Engineering Department,
University of Cambridge,
Cambridge CB2 1PZ, UK
e-mail: kk377@cam.ac.uk

Santosh Hemchandra

Institute of Aerodynamics,
RWTH Aachen University,
Aachen 52062, Germany
e-mail: s.hemchandra@aia.rwth-aachen.de

Matthew P. Juniper

Engineering Department,
University of Cambridge,
Cambridge CB2 1PZ, UK
e-mail: mpj1001@cam.ac.uk

1Corresponding author.

2Present address: Department of Aerospace Engineering, Indian Institute of Science, Bangalore 560012, India.

Contributed by the International Gas Turbine Institute (IGTI) of ASME for publication in the Journal of Engineering for Gas Turbines and Power. Manuscript received September 18, 2012; final manuscript received November 30, 2012; published online May 20, 2013. Editor: David Wisler.

J. Eng. Gas Turbines Power 135(6), 061502 (May 20, 2013) (9 pages) Paper No: GTP-12-1365; doi: 10.1115/1.4023305 History: Received September 18, 2012; Revised November 30, 2012

## Abstract

Nonlinear analysis of thermoacoustic instability is essential for the prediction of the frequencies, amplitudes, and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper, an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, however, it can be used to calculate the amplitudes of limit cycles and their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modeled using a nonlinear kinematic model based on the $G$-equation, while the acoustics of planar waves in the tube are governed by linearized momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state for single-mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitude-dependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, along with twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact, resulting in a larger limit cycle amplitude. Multimode simulations show that, in some situations, the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems.

<>
Your Session has timed out. Please sign back in to continue.

## References

Lieuwen, T. C. and Yang, V., 2005, “Combustion Instabilities in Gas Turbine Engines,” Progress in Astronautics and Aeronautics, Vol. 210, AIAA, Reston, Va.
Zinn, B. T. and Lores, M. E., 1972, “Application of the Galerkin Methods in the Solution of Nonlinear Axial Combustion Instability Problems in Liquid Rockets,” Combust. Sci. Technol., 4, pp. 269–278.
Culick, F. E. C., 1976, “Nonlinear Behaviour of Acoustic Waves in Combustion Chambers—Part 1,” Acta Astronaut., 3, pp. 715–734.
Culick, F. E. C., 1976, “Nonlinear Behaviour of Acoustic Waves in Combustion Chambers—Part 2,” Acta Astronaut., 3, pp. 735–757.
Yang, V. and Culick, F. E. C., 1990, “On the Existence and Stability of Limit Cycles for Transverse Acoustic Oscillations in a Cylindrical Combustion Chamber. I: Standing Modes,” Combust. Sci. Technol., 72(1), pp. 37–65.
Yang, V., Kim, S. I., and Culick, F. E. C., 1990, “Triggering of Longitudinal Pressure Oscillations in Combustion Chambers. I: Nonlinear Gas Dynamics,” Combust. Sci. Technol., 72(4), pp. 183–214.
Baum, J. D., Levine, J. N., and Lovine, R. L., 1988, “Pulsed Instability in Rocket Motors: A Comparison Between Predictions and Experiment,” J. Propul. Power, 4(4), pp. 308–316.
Wicker, J. M., Greene, W. D., Kim, S-I., and Yang, V., 1996, “Triggering of Longitudinal Pressure Oscillations in Combustion Chambers. I: Nonlinear Combustion Response,” J. Propul. Power, 12(6), pp. 1148–1158.
Ananthkrishnan, N., Deo, S., and Culick, F. E. C., 2005, “Reduced-Order Modeling and Dynamics of Nonlinear Acoustic Waves in a Combustion Chamber,” Combust. Sci. Technol., 177(2), pp. 221–248.
Poinsot, T., and Candel, S., 1988, “A Nonlinear Model for Ducted Flame Combustion Instabilities,” Combust. Sci. Technol., 61, pp. 121–153.
Dowling, A. P., 1997, “Nonlinear Self-Excited Oscillations of a Ducted Flame,” J. Fluid Mech., 346, pp. 271–290.
Dowling, A. P., 1999, “A Kinematic Model of a Ducted Flame,” J. Fluid Mech., 394, pp. 51–72.
Stow, S. R., and Dowling, A. P., 2004, “Low-Order Modelling of Thermoacoustic Limit Cycles,” ASME Turbo Expo, Vienna, Austria, June 14–17, ASME Paper No. GT2004-54245.
Stow, S. R., and Dowling, A. P., 2008, “A Time-Domain Network Model for Nonlinear Thermoacoustic Oscillations,” ASME Turbo Expo, Berlin, June 9–13, ASME Paper No. GT2008-50770.
Lieuwen, T., 2005, “Nonlinear Kinematic Response of Premixed Flames to Harmonic Velocity Disturbances,” Proc. Combust. Inst., 29, pp. 99–105.
Noiray, N., Durox, D., Schuller, T., and Candel, S. M., 2008, “A Unified Framework for Nonlinear Combustion Instability Analysis Based on the Flame Describing Function,” J. Fluid Mech., 615, pp. 139–167.
Moeck, J. P., Bothien, M. R., Schimek, S., Lacarelle, A., and Paschereit, C. O., 2008, “Subcritical Thermoacoustic Instabilities in a Premixed Combustor,” 14th AIAA/CEAS Aeroacoustics Conference, Vancouver, Canada, May 5–7, AIAA Paper No. 2008-2946.
Subramanian, P., Gupta, V., Tulsyan, B., and Sujith, R. I., 2010, “Can Describing Function Technique Predict Bifurcations in Thermoacoustic Systems?,” 16th AIAA/CEAS Aeroacoustics Conference, Stockholm, Sweden, June 7–9, AIAA Paper No. 2010-3860.
Matveev, I., 2003, “Thermo-Acoustic Instabilities in the Rijke Tube: Experiments and Modeling,” Ph.D. thesis, CalTech, Pasadena, CA.
Juniper, M. P., 2011, “Triggering in the Horizontal Rijke Tube: Non-Normality, Transient Growth and Bypass Transition,” J. Fluid Mech., 667, pp. 272–308.
Schuller, T., Durox, D., and Candel, S., 2003, “A Unified Model for the Prediction of Laminar Flame Transfer Functions: Comparisons Between Conical and v-Flame Dynamics,” Combust. Flame, 134, pp. 21–34.
Preetham, Santosh, H., and Lieuwen, T., 2008, “Dynamics of Premixed Flames Forced by Harmonic Velocity Disturbances,” J. Propul. Power, 24(6), pp. 1390–1402.
Shreekrishna, Hemchandra, S., and Lieuwen, T., 2010, “Premixed Flame Response to Equivalence Ratio Perturbations,” Combust. Theory Model., 14, pp. 681–714.
Abu-Orf, G., and Cant, R. S., 1996, “Reaction Rate Modelling for Premixed Turbulent Methane-Air Flames,” Proceedings of the Joint Meeting of Spanish, Portuguese, Swedish and British Sections of the Combustion Institute, Madeira, Portugal, April 1–4.
Fleifil, M., Annaswamy, A., Ghoneim, Z., and Ghoneim, A., 1996, “Response of a Laminar Premixed Flame to Flow Oscillations: A Kinematic Model and Thermoacoustic Instability Results,” Combust. Flame, 106, pp. 487–510.
Ducruix, S., Durox, D., and Candel, S., 2000, “Theoretical and Experimental Determination of the Transfer Function of a Laminar Premixed Flame,” Proc. Combust. Inst., 28, pp. 765–773.
Baillot, F., Durox, D., and Prud'homme, R., 1992, “Experimental and Theoretical Study of a Premixed Vibrating Flame,” Combust. Flame, 88(2), pp. 149–168.
Baillot, F., Bourehla, A., and Durox, D., 1996, “The Characteristics Method and Cusped Flame Fronts,” Combust. Sci. Technol., 112(1), pp. 327–350.
Cuquel, A., Durox, D., and Schuller, T., 2011, “Theoretical and Experimental Determination of the Flame Transfer Function of Confined Premixed Conical Flames,” 7th Mediterranean Combustion Symposium, Cagliari, Sardinia, Italy, September 11–15.
Michalke, A., 1971, “Instability of a Compressible Circular Free Jet With Consideration of the Influence of the Jet Boundary Thickness,” NASA Report No. TM 75190.
Jiang, G.-S., and Peng, D., 2000, “Weighted ENO Schemes for Hamilton-Jacobi Equations,” SIAM J. Sci. Comput. (USA), 6, pp. 2126–2143.
Gottlieb, S., and Shu, C., 1998, “Total Variation Diminishing Runge-Kutta Schemes,” Math. Comput., 67, pp. 73–85.
Peng, D., Merriman, B., Osher, S., Zhao, H., and Kang, M., 1999, “A PDE-Based Fast Local Level Set Method,” J. Comput. Phys., 155(2), pp. 410–438.
Hemchandra, S., 2009, “Dynamics of Turbulent Premixed Flames in Acoustic Fields,” Ph.D. thesis, Georgia Institute of Technology, Atlanta, GA.
Karimi, N., Brear, M. J., Jin, S.-H., and Monty, J. P., 2009, “Linear and Non-Linear Forced Response of a Conical, Ducted, Laminar Premixed Flame,” Combust. Flame, 156, pp. 2201–2212.
Smereka, P., 2006, “The Numerical Approximation of a Delta Function With Application to Level Set Methods,” J. Comput. Phys., 211, pp. 77–90.
Dowling, A. P., and Morgans, A. S., 2005, “Feedback Control of Combustion Oscillations,” Ann. Rev. Fluid Mech., 37, pp. 151–182.
Rayleigh, J. W. S., 1945, The Theory of Sound, Vol. II, Dover, New York.
Durox, D., Schuller, T., Noiray, N., and Candel, S., 2009, “Experimental Analysis of Nonlinear Flame Transfer Functions for Different Flame Geometries,” Proc. Combust. Inst., 32(1), pp. 1391–1398.
Polifke, W., and Lawn, C., 2007, “On the Low-Frequency Limit of Flame Transfer Functions,” Combust. Flame, 151(3), pp. 437–451.

## Figures

Fig. 1

Schematic of the two-dimensional slot stabilized premixed flame in a duct: L0 is the length of the duct, x˜f is the flame position along the duct, α = 0.02 is the fraction of the duct cross-sectional area occupied by the burner and the flow is from left to right

Fig. 2

Instantaneous images of the flame during one forcing cycle; φ = 1.06, βf = 2.14, ɛ = 0.25, St = 1, and K = 2.5. The thick black line represents the slot burner and the thin black curve is the flame surface. Note the formation of sharp cusps towards the products, which is a distinct characteristic of premixed flames seen in experiments [26,27,35].

Fig. 3

Flame describing function FDF(ω,ɛ) = Q·˜'/Q·˜0/u˜'/u˜0: (a) gain, and (b) phase K = 2.5, βf = 2.14, φ = 1.06

Fig. 4

CIRCE diagram: driving (Ψdriv) – damping (Ψdamp) as a function of amplitude ɛ for thermoacoustic systems with different fundamental frequencies f*: (a) K = 1.0, (b) K = 1.5, and (c) K = 2.5

Fig. 5

CIRCE for thermoacoustic systems with different fundamental frequencies (duct lengths): K = 2.5, φ = 1.06, and βf = 2.14. (a) Slice of Fig. 4(c) at f* = 0.4: A is the stable limit cycle, and (b) slice of Fig. 4(c) at f* = 1.0; A is the unstable limit cycle, and B is the stable limit cycle.

Fig. 6

CIRCE for thermoacoustic systems with different fundamental frequencies (duct lengths): K = 2.5, φ = 1.06, and βf = 2.14. (a) Slice of Fig. 4(c) at f* = 1.2: A and C are the stable limit cycles, B is the unstable limit cycle, and (b) slice of Fig. 4(c) at f* = 1.8: A and C are the unstable limit cycles, B and D are the stable limit cycles.

Fig. 7

Amplitude-dependence of the fundamental of heat release rate oscillations and the phase between the heat release rate and velocity perturbations: K = 2.5, φ = 1.06, and βf = 2.14. (a) q1 at St = 1.0, (b) Δϕ1 at St = 1.0, and (c) sin(Δϕ1) at St = 1.0. The gray shaded area shows regions where driving is negative, while white shows where it is positive. (d) q1 at St = 0.8.

Fig. 8

Time domain calculations for self-excited thermoacoustic systems with one mode and twenty modes: f* = 0.87, K = 1.5, βf = 2.14, φ = 1.06, and ζ = 0.05(c1 = 0.03, c2 = 0.02). (a) Time trace of the acoustic velocity perturbations at the flame location xf, and (b) time trace of the contributions of the first five Galerkin modes to the acoustic velocity at the flame location xf of the 20-mode system.

Fig. 9

Nonlinear heat release rate oscillations in thermoacoustic systems with one mode and twenty modes: f* = 0.87, K = 1.5, βf = 2.14, φ = 1.06, and ζ = 0.05(c1 = 0.03, c2 = 0.02). (a) Time trace of the heat release rate (phase-shifted to avoid overlap of figures), and (b) Fourier transforms of the heat release rate.

Fig. 10

CIRCE diagram: the intersection of driving and damping curves represents the frequency domain prediction for the limit cycle amplitude. Time-domain limit cycle locations are marked by symbols: f* = 0.87, K = 1.5, βf = 2.14, φ = 1.06, and ζ = 0.05(c1 = 0.03, c2 = 0.02).

## Errata

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 Proceedings Articles
Related eBook Content
Topic Collections