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

Impact of Heat Release Distribution on the Spinning Modes of an Annular Combustor With Multiple Matrix Burners

[+] Author and Article Information
Davide Laera

DMMM,
Sez. Macchine ed Energetica,
Politecnico di Bari,
Via Re David 200,
Bari 70125, Italy
e-mail: davide.laera@poliba.it

Kevin Prieur

Laboratoire EM2C,
CNRS, CentraleSupélec,
Université Paris-Saclay,
Grande Voie des Vignes,
Chatenay-Malabry cedex 92295, France;
Safran Tech, E&P,
Châteaufort, CS 80112,
Magny-Les-Hameaux 78772, France

Daniel Durox, Thierry Schuller, Sébastien Candel

Laboratoire EM2C,
CNRS, CentraleSupélec,
Université Paris-Saclay,
Grande Voie des Vignes,
Chatenay-Malabry cedex 92295, France

Sergio M. Camporeale

DMMM,
Sez. Macchine ed Energetica,
Politecnico di Bari,
Via Re David 200,
Bari 70125, Italy

1Corresponding author.

Contributed by the Combustion and Fuels Committee of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received August 7, 2016; final manuscript received September 12, 2016; published online January 10, 2017. Assoc. Editor: Song-Charng Kong.

J. Eng. Gas Turbines Power 139(5), 051505 (Jan 10, 2017) (10 pages) Paper No: GTP-16-1393; doi: 10.1115/1.4035207 History: Received August 07, 2016; Revised September 12, 2016

The present article reports a numerical analysis of instability coupled by a spinning mode in an annular combustor. This corresponds to experiments carried out on the MICCA test facility equipped with 16 matrix burners. Each burner response is represented by means of a global experimental flame describing function (FDF). A harmonic balance nonlinear stability analysis is carried out by combining the FDF with a Helmholtz solver to determine the system dynamics trajectories in a frequency-growth rate plane. The influence of the distribution of the volumetric heat release corresponding to each burner is investigated in a first stage. Even though each of the 16 burners is compact with respect to the transverse mode wavelength, and the 16 flames occupy the same volume, this distribution of heat release is not compact in the azimuthal direction and simulations reveal an influence of this volumetric distribution on frequencies and growth rates. This study emphasizes the importance of providing a suitable description of the flame zone geometrical extension and correspondingly an adequate representation of the level of heat release rate fluctuation per unit volume. It is found that these two items can be deduced from a knowledge of the heat release distribution under steady-state operating conditions. Once the distribution of the heat release fluctuations is unequivocally defined, limit cycle simulations are performed. For the conditions explored, simulations retrieve the spinning nature of the self-sustained mode that was identified in the experiments both in the plenum and in the combustion chamber.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Lieuwen, T. , and Yang, V. , eds., 2005, Combustion Instabilities in Gas Turbine Engines: Operational Experience, Fundamental Mechanisms and Modeling (Progress in Astronautics and Aeronautics), American Institute of Aeronautics and Astronautics, Reston, VA.
Candel, S. M. , 1992, “ Combustion Instabilities Coupled by Pressure Waves and Their Active Control,” Symp. Combust., 24(1), pp. 1277–1296. [CrossRef]
Bourgouin, J. F. , Durox, D. , Moeck, J. P. , Schuller, T. , and Candel, S. , 2015, “ A New Pattern of Instability Observed in an Annular Combustor: The Slanted Mode,” Proc. Combust. Inst., 35(3), pp. 3237–3244. [CrossRef]
Prieur, K. , Durox, D. , Schuller, T. , and Candel, S. , “ A Hysteresis Phenomenon Leading to Spinning or Standing Azimuthal Instabilities in an Annular Combustor,” Combust. Flame in press. [CrossRef]
Durox, D. , Bourgouin, J. , Moeck, J. , Philip, M. , Schuller, T. , and Candel, S. , 2013, “ Nonlinear Interactions in Combustion Instabilities Coupled by Azimuthal Acoustic Modes,” International Workshop on Non-Normal and Nonlinear Effects in Aero-and Thermo-Acoustics, Munich, Germany, June 18–21, pp. 18–21. https://mediatum.ub.tum.de/doc/1253524/1253524.pdf
Bourgouin, J.-F. , Durox, D. , Moeck, J. P. , Schuller, T. , and Candel, S. , 2014, “ Characterization and Modeling of a Spinning Thermoacoustic Instability in an Annular Combustor Equipped With Multiple Matrix Injectors,” ASME J. Eng. Gas Turbines Power, 137(2), p. 21503. [CrossRef]
Worth, N. A. , and Dawson, J. R. , 2013, “ Self-Excited Circumferential Instabilities in a Model Annular Gas Turbine Combustor: Global Flame Dynamics,” Proc. Combust. Inst., 34(2), pp. 3127–3134. [CrossRef]
Worth, N. A. , and Dawson, J. R. , 2013, “ Modal Dynamics of Self-Excited Azimuthal Instabilities in an Annular Combustion Chamber,” Combust. Flame, 160(11), pp. 2476–2489. [CrossRef]
Krebs, W. , Flohr, P. , Prade, B. , and Hoffmann, S. , 2002, “ Thermoacoustic Stability Chart for High-Intensity Gas Turbine Combustion Systems,” Combust. Sci. Technol., 174(7), pp. 99–128. [CrossRef]
Bothien, M. R. , Noiray, N. , and Schuermans, B. , 2015, “ Analysis of Azimuthal Thermo-Acoustic Modes in Annular Gas Turbine Combustion Chambers,” ASME J. Eng. Gas Turbines Power, 137(6), p. 61505. [CrossRef]
Mastrovito, M. , Camporeale, S. , Forte, A. , and Fortunato, B. , 2005, “ Analysis of Pressure Oscillations Data in Gas Turbine Annular Combustion Chamber Equipped With Passive Damper,” ASME Paper No. GT2005-69056.
Schuermans, B. , Paschereit, C. , and Monkewitz, P. , 2006, “ Non-Linear Combustion Instabilities in Annular Gas-Turbine Combustors,” AIAA Paper No. 2006-549.
Noiray, N. , Bothien, M. , and Schuermans, B. , 2011, “ Investigation of Azimuthal Staging Concepts in Annular Gas Turbines,” Combust. Theory Model., 15(5), pp. 585–606. [CrossRef]
Ghirardo, G. , and Juniper, M. P. , 2013, “ Azimuthal Instabilities in Annular Combustors: Standing and Spinning Modes,” Proc. R. Soc. A, 469(2157), p. 20130232. [CrossRef]
Noiray, N. , and Schuermans, B. , 2013, “ On the Dynamic Nature of Azimuthal Thermoacoustic Modes in Annular Gas Turbine Combustion Chambers,” Proc. R. Soc. A, 469(2151), p. 20120535. [CrossRef]
Ghirardo, G. , Ćosić, B. , Juniper, M. P. , and Moeck, J. P. , 2015, “ State-Space Realization of a Describing Function,” Nonlinear Dyn., 82(1), pp. 9–28. [CrossRef]
Ghirardo, G. , Juniper, M. P. , and Moeck, J. P. , 2015, “ Stability Criteria for Standing and Spinning Waves in Annular Combustors,” ASME Paper No. GT2015-43127.
Selle, L. , Lartigue, G. , Poinsot, T. , Koch, R. , Schildmacher, K. U. , Krebs, W. , Prade, B. , Kaufmann, P. , and Veynante, D. , 2004, “ Compressible Large Eddy Simulation of Turbulent Combustion in Complex Geometry on Unstructured Meshes,” Combust. Flame, 137(4), pp. 489–505. [CrossRef]
Wolf, P. , Balakrishnan, R. , Staffelbach, G. , Gicquel, L. Y. M. , and Poinsot, T. , 2012, “ Using LES to Study Reacting Flows and Instabilities in Annular Combustion Chambers,” Flow, Turbul. Combust., 88(1), pp. 191–206. [CrossRef]
Parmentier, J.-F. , Salas, P. , Wolf, P. , Staffelbach, G. , Nicoud, F. , and Poinsot, T. , 2012, “ A Simple Analytical Model to Study and Control Azimuthal Instabilities in Annular Combustion Chambers,” Combust. Flame, 159(7), pp. 2374–2387. [CrossRef]
Bauerheim, M. , Parmentier, J. F. , Salas, P. , Nicoud, F. , and Poinsot, T. , 2014. “ An Analytical Model for Azimuthal Thermoacoustic Modes in an Annular Chamber Fed by an Annular Plenum,” Combust. Flame, 161(5), pp. 1374–1389. [CrossRef]
Pankiewitz, C. , and Sattelmayer, T. , 2003, “ Time Domain Simulation of Combustion Instabilities in Annular Combustors,” ASME J. Eng. Gas Turbines Power, 125(3), pp. 677–685. [CrossRef]
Stow, S. R. , and Dowling, A. P. , 2009, “ A Time-Domain Network Model for Nonlinear Thermoacoustic Oscillations,” ASME J. Eng. Gas Turbines Power, 131(3), p. 031502. [CrossRef]
Morgans, A. S. , and Stow, S. R. , 2007, “ Model-Based Control of Combustion Instabilities in Annular Combustors,” Combust. Flame, 150(4), pp. 380–399. [CrossRef]
Nicoud, F. , Benoit, L. , Sensiau, C. , and Poinsot, T. , 2007, “ Acoustic Modes in Combustors With Complex Impedances and Multidimensional Active Flames,” AIAA J., 45(2), pp. 426–441. [CrossRef]
Campa, G. , and Camporeale, S. M. , 2014, “ Prediction of the Thermoacoustic Combustion Instabilities in Practical Annular Combustors,” ASME J. Eng. Gas Turbines Power, 136(9), p. 91504. [CrossRef]
Silva, C. F. , Nicoud, F. , Schuller, T. , Durox, D. , and Candel, S. , 2013, “ Combining a Helmholtz Solver With the Flame Describing Function to Assess Combustion Instability in a Premixed Swirled Combustor,” Combust. Flame, 160(9), pp. 1743–1754. [CrossRef]
Camporeale, S. M. , Fortunato, B. , and Campa, G. , 2010, “ A Finite Element Method for Three-Dimensional Analysis of Thermo-Acoustic Combustion Instability,” ASME J. Eng. Gas Turbines Power, 133(1), p. 11506. [CrossRef]
Bourgouin, J.-F. , Durox, D. , Moeck, J. P. , Schuller, T. , and Candel, S. , 2013, “ Self-Sustained Instabilities in an Annular Combustor Coupled by Azimuthal and Longitudinal Acoustic Modes,” ASME Paper No. GT2013-95010.
Durox, D. , Prieur, K. , Schuller, T. , and Candel, S. , 2015, “ Different Flame Patterns Linked With Swirling Injector Interactions in an Annular Combustor,” ASME Paper No. GT2015-42034.
Pierce, A. D. , 1991, Acoustics: An Introduction to Its Physical Principles and Applications, Acoustical Society of America, Melville, NY.
Munjal, M. L. , 1987, Acoustics of Ducts and Mufflers With Application to Exhaust and Ventilation System Design, Wiley, New York.
Noiray, N. , Durox, D. , Schuller, T. , and Candel, S. , 2008, “ A Unified Framework for Nonlinear Combustion Instability Analysis Based on the Flame Describing Function,” J. Fluid Mech., 615, pp. 139–167. [CrossRef]
Boudy, F. , Durox, D. , Schuller, T. , Jomaas, G. , and Candel, S. , 2011, “ Describing Function Analysis of Limit Cycles in a Multiple Flame Combustor,” ASME J. Eng. Gas Turbines Power, 133(6), p. 061502. [CrossRef]
Palies, P. , Durox, D. , Schuller, T. , and Candel, S. , 2010, “ The Combined Dynamics of Swirler and Turbulent Premixed Swirling Flames,” Combust. Flame, 157(9), pp. 1698–1717. [CrossRef]
Palies, P. , Durox, D. , Schuller, T. , and Candel, S. , 2011, “ Experimental Study on the Effect of Swirler Geometry and Swirl Number on Flame Describing Functions,” Combust. Sci. Technol., 183(7), pp. 704–717. [CrossRef]
Schuller, T. , Durox, D. , and Candel, S. , 2003, “ A Unified Model for the Prediction of Laminar Flame Transfer Functions,” Combust. Flame, 134(1–2), pp. 21–34. [CrossRef]
Sensiau, C. , Nicoud, F. , and Poinsot, T. , 2008, “ Thermoacoustic Analysis of an Helicopter Combustion Chamber,” AIAA Paper No. AIAA-2008-2947. http://arc.aiaa.org/doi/pdf/10.2514/6.2008-2947
Cuquel, A. , Silva, C. , Nicoud, F. , Durox, D. , and Schuller, T. , 2013, “ Prediction of the Nonlinear Dynamics of a Multiple Flame Combustor by Coupling the Describing Function Methodology With a Helmholtz Solver,” ASME Paper No. GT2013-95659.
Comsol, A. , 2005, COMSOL Multiphysics User's Guide, The Comsol Group, Stockholm, Sweden.
Lehoucq, R. B. , Sorensen, D. C. , and Yang, C. , 1998, ARPACK Users’ Guide: Solution of Large-Scale Eigenvalue Problems With Implicitly Restarted Arnoldi Methods, Vol. 6., SIAM, Philadelphia, PA.
Kinsler, L. E. , Frey, A. R. , Coppens, A. B. , and Sanders, J. V. , 1999, “ Fundamentals of Acoustics,” Fundamentals of Acoustics, 4th ed., Lawrence E. Kinsler , A. R. Frey , A. B. Coppens , and J. V. Sanders , eds., Wiley-VCH, Weinheim, Germany, p. 560.
Palies, P. , Durox, D. , Schuller, T. , and Candel, S. , 2011, “ Nonlinear Combustion Instability Analysis Based on the Flame Describing Function Applied to Turbulent Premixed Swirling Flames,” Combust. Flame, 158(10), pp. 1980–1991. [CrossRef]
Vagelopoulos, C. M. , Egolfopoulos, F. N. , and Law, C. K. , 1994, “ Further Considerations on the Determination of Laminar Flame Speeds With the Counterflow Twin-Flame Technique,” Symp. Combust., 25(1), pp. 1341–1347. [CrossRef]
Bauerheim, M. , Salas, P. , Nicoud, F. , and Poinsot, T. , 2014, “ Symmetry Breaking of Azimuthal Thermo-Acoustic Modes in Annular Cavities: A Theoretical Study,” J. Fluid Mech., 760, pp. 431–465. [CrossRef]

Figures

Grahic Jump Location
Fig. 7

Details of the flame region domain and of the reference point for the velocity fluctuations in one sector of the annular combustor

Grahic Jump Location
Fig. 6

Interpolated flame describing function (FDF). The experimental data are displayed as white dots: (a) gain and (b) phaseφ.

Grahic Jump Location
Fig. 5

Pressure measurements outside the combustion chamber used to determine the end correction for the first 1A–1L azimuthal mode

Grahic Jump Location
Fig. 4

(a) Three-dimensional model of the MICCA chamber with the temperature field in Kelvin and (b) longitudinal cut showing geometrical details of the matrix injector representation

Grahic Jump Location
Fig. 3

Pressure signals recorded by microphones in the plenum (a) and the combustion chamber (b) for the spinning mode analyzed in this article obtained at an equivalence ratio ϕ = 0.98 and a bulk velocity ub = 0.49 ms−1

Grahic Jump Location
Fig. 2

Top view of the MICCA chamber with microphone locations indicated

Grahic Jump Location
Fig. 1

(a) Photograph of the MICCA combustor with a zoomed view of a matrix injector and the waveguide outlet. (b) Schematic representation of the experimental setup with details of the matrix injector position (c) and of the mixture feeding line (d).

Grahic Jump Location
Fig. 14

Trajectory of the MICCA combustor colored with the growth rate projected over the contour plot of the gain of the FDF. The zoom frame highlights the limit cycle condition for a vanishing damping rate (triangular mark) and for a finite damping rate (circular mark).

Grahic Jump Location
Fig. 8

Three different flame domains characterized by the same flame volume Vf  = 16.6 cm−3. Azimuthal and longitudinal extensions: (a) lx = 6.8 cm, lz = 0.4 cm; (b) lx = 4.2 cm, lz = 1.2 cm; and (c) lx = 3.6 cm, lz = 1.6 cm.

Grahic Jump Location
Fig. 9

Dynamical trajectories in the frequency (f) growth rate (α) plane colored by (a) the velocity fluctuation level |u′/u¯| and (b) the gain of the FDF. Circular marks indicate the flame volume shown in Fig. 8(a), triangular marks refer to the case shown in Fig. 8(b), and square marks to the model of Fig. 8(c).

Grahic Jump Location
Fig. 10

Four different flame volumes characterized by the same longitudinal dimension but different extension in the azimuthal direction. Azimuthal dimension (lx): (a) lx = 6.8 cm; (b) lx = 4.2 cm; (c) lx = 3.6 cm; and (d) lx = 0.8 cm. Axial extension lz = 0.8 cm.

Grahic Jump Location
Fig. 11

Trajectories in the frequency (f) growth rate (α) plane colored by the velocity fluctuation level |u′/u¯|. Circle marks indicate the flame volume shown in Fig. 10(a), triangular marks refer to the case shown in Fig. 10(b), square marks to the model of Fig. 10(c) and diamonds indicate the last one shown in Fig.10(d).

Grahic Jump Location
Fig. 12

(a) Theoretical conical flame of diameter a and height h over a cylindrical volume. (b) Image of the flame recorded in the longitudinal combustor equipped with the perforated plate injector in a steady-state configuration.

Grahic Jump Location
Fig. 13

Acoustic response of the MICCA combustor from 300 Hz to 500 Hz measured by microphones located in the (a) plenum and (b) combustion chamber. The frequency bandwidth Δf determined at half maximum provides the damping rate in both volumes.

Grahic Jump Location
Fig. 15

Pressure mode magnitude |p̂| with pressure contour lines plotted on a cylindrical surface equidistantly located from the lateral walls

Grahic Jump Location
Fig. 16

(a) Pressure distribution along the azimuthal direction in the plenum (continuous line) and in the combustion chamber (dashed line). (b) Pressure phase evolution along the azimuthal direction in the plenum (continuous line) and in the combustion chamber (dashed line).

Grahic Jump Location
Fig. 17

Comparison between the plenum PM1 signal (continuous line) and chamber MC1 signal (dot-dashed line). When the angular shift due to the microphones position is taken into account, the two signals manifest a phase shift of approx. 0.18 rad.

Tables

Errata

Discussions

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