0
Research Papers

Thermoacoustics of Can-Annular Combustors

[+] Author and Article Information
G. Ghirardo

Ansaldo Energia Switzerland,
Römerstrasse 36,
Baden 5401, Switzerland
e-mail: giulio.ghirardo@ansaldoenergia.com

C. Di Giovine

Ansaldo Energia Switzerland,
Römerstrasse 36,
Baden 5401, Switzerland

J. P. Moeck

Institut für Strömungsmechanik und
Technische Akustik,
Technische Universität Berlin,
Berlin 10623, Germany;
Department of Energy and
Process Engineering,
Norwegian University of Science and
Technology,
Trondheim 7491, Norway

M. R. Bothien

Ansaldo Energia Switzerland,
Römerstrasse 36,
Baden 5401, Switzerland

1Corresponding author.

Manuscript received June 22, 2018; final manuscript received June 28, 2018; published online September 14, 2018. Editor: Jerzy T. Sawicki.

J. Eng. Gas Turbines Power 141(1), 011007 (Sep 14, 2018) (10 pages) Paper No: GTP-18-1287; doi: 10.1115/1.4040743 History: Received June 22, 2018; Revised June 28, 2018

Can-annular combustors consist of a set of independent cans, connected on the upstream side to the combustor plenum and on the downstream side to the turbine inlet, where a transition duct links the round geometry of each can with the annular segment of the turbine inlet. Each transition duct is open on the sides toward the adjacent transition ducts, so that neighboring cans are acoustically connected through a so-called cross-talk open area. This theoretical, numerical, and experimental work discusses the effect that this communication has on the thermoacoustic frequencies of the combustor. We show how this communication gives rise to axial and azimuthal modes, and that these correspond to particularly synchronized states of axial thermoacoustic oscillations in each individual can. We show that these combustors typically show clusters of thermoacoustic modes with very close frequencies and that a slight loss of rotational symmetry, e.g., a different acoustic response of certain cans, can lead to mode localization. We corroborate the predictions of azimuthal modes, clusters of eigenmodes, and mode localization with experimental evidence.

FIGURES IN THIS ARTICLE
<>
Copyright © 2019 by ASME
Your Session has timed out. Please sign back in to continue.

References

Pennell, D. A. , Bothien, M. R. , Ciani, A. , Granet, V. , Singla, G. , Thorpe, S. , Wickstroem, A. , Oumejjoud, K. , and Yaquinto, M. , 2017, “An Introduction to the Ansaldo GT36 Constant Pressure Sequential Combustor,” ASME Paper No. GT2017-64790.
Noiray, N. , Bothien, M. , and Schuermans, B. , 2011, “Investigation of Azimuthal Staging Concepts in Annular Gas Turbines,” Combust. Theory Modell., 15(5), pp. 585–606. [CrossRef]
Connor, J. O. , Acharya, V. , and Lieuwen, T. , 2015, “Transverse Combustion Instabilities: Acoustic, Fluid Mechanic, and Flame Processes,” Prog. Energy Combust. Sci., 49, pp. 1–39. [CrossRef]
Bothien, M. , 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. 061505. [CrossRef]
Bauerheim, M. , Nicoud, F. , and Poinsot, T. , 2016, “Progress in Analytical Methods to Predict and Control Azimuthal Combustion Instability Modes in Annular Chambers,” Phys. Fluids, 28(2), p. 021303. [CrossRef]
Krebs, W. , Bethke, S. , Lepers, J. , Flohr, P. , Prade, B. , Johnson, C. , and Sattinger, S. , 2005, “Thermoacoustic Design Tools and Passive Control: Siemens Power Generation Approaches,” Combustion Instabilities in Gas Turbine Engines: Operational Experience, Fundamental Mechanisms, and Modeling (Progress in Astronautics and Aeronautics), T. Lieuwen and V. Yang , eds., American Institute of Aeronautics and Astronautics, Arlington, TX, pp. 89–112.
Macquisten, M. A. , Whiteman, M. , and Moran, A. J. , 2004, “Application of Low Order Thermo-Acoustic Network to DLE Staged Combustor,” ASME Paper No. GT2004-54161.
Kaufmann, P. , Krebs, W. , Valdes, R. , and Wever, U. , 2008, “3D Thermoacoustic Properties of Single Can and Multi Can Combustor Configurations,” ASME Paper No. GT2008-50755.
Farisco, F. , Panek, L. , Janus, B. , and Kok, J. B. W. , 2015, “Numerical Investigation of the Thermo-Acoustic Influence of the Turbine on the Combustor,” ASME Paper No. GT2015-42071.
Farisco, F. , Panek, L. , and Kok, J. B. W. , 2016, “Thermo-Acoustic Cross-Talk Between Cans in a Can-Annular Combustor,” International Symposium on Thermoacoustic Instabilities in Gas Turbines and Rocket Engines, Munich, Germany, May 30–June 2, pp. 1–12.
Panek, L. , Huth, M. , and Farisco, F. , 2017, “Thermo-Acoustic Characterization of Can-Can Interaction of a Can-Annular Combustion System Based on Unsteady CFD LES Simulation,” First Global Power and Propulsion Forum, Zurich, Switzerland, Jan. 16–18, Paper No. GPPF-2017-81.
Farisco, F. , Panek, L. , and Kok, J. B. , 2017, “Thermo-Acoustic Cross-Talk Between Cans in a Can-Annular Combustor,” Int. J. Spray Combust. Dyn., 9(4), pp. 452–469. [CrossRef]
Stow, S. R. , and Dowling, A. P. , 2003, “Modelling of Circumferential Modal Coupling Due to Helmholtz Resonators,” ASME Paper No. GT2003-38168.
Mensah, G. A. , and Moeck, J. P. , 2016, “Efficient Computation of Thermoacoustic Modes in Annular Combustion Chambers Based on Bloch-Wave Theory,” ASME J. Eng. Gas Turbines Power, 138(8), p. 081502. [CrossRef]
Bauerheim, M. , Cazalens, M. , and Poinsot, T. , 2015, “A Theoretical Study of Mean Azimuthal Flow and Asymmetry Effects on Thermo-Acoustic Modes in Annular Combustors,” Proc. Combust. Inst., 35(3), pp. 3219–3227. [CrossRef]
Ghirardo, G. , Boudy, F. , and Bothien, M. R. , 2018, “Amplitude Statistics Prediction in Thermoacoustics,” J. Fluid Mech., 844, pp. 216–246. [CrossRef]
Brillouin, L. , 1953, Wave Propagation in Periodic Structures: Electric Filters and Crystal Lattices, 2nd ed., Dover Publications, Mineola, NY, pp. 139–140.
Munjal, M. , 1987, Acoustics of Ducts and Mufflers, Wiley, Hoboken, NJ.
Bender, C. , and Orszag, S. , 1978, Advanced Mathematical Methods for Scientists and Engineers, McGraw-Hill Book, New York.
Miller, P. , 2006, Applied Asymptotic Analysis, American Mathematical Society, Providence, RI.
Stow, S. R. , Dowling, A. P. , and Hynes, T. P. , 2002, “Reflection of Circumferential Modes in a Choked Nozzle,” J. Fluid Mech., 467, pp. 215–239. [CrossRef]
Marble, F. E. , and Candel, S. , 1977, “Acoustic Disturbance From Gas Non-Uniformities Convected Through a Nozzle,” J. Sound Vib., 55(2), pp. 225–243. [CrossRef]
Bauerheim, M. , Duran, I. , Livebardon, T. , Wang, G. , Moreau, S. , and Poinsot, T. , 2016, “Transmission and Reflection of Acoustic and Entropy Waves Through a Stator-Rotor Stage,” J. Sound Vib., 374, pp. 260–278. [CrossRef]
Morgans, A. S. , and Duran, I. , 2016, “Entropy Noise: A Review of Theory, Progress and Challenges,” Int. J. Spray Combust. Dyn., 8(4), pp. 285–298. [CrossRef]
Trefethen, L. N. , 2000, Spectral Methods in MATLAB, Society for Industrial and Applied Mathematics, Philadelphia, PA.
Krebs, W. , Walz, G. , Flohr, P. , and Hoffmann, S. , 2001, “Modal Analysis of Annular Combustors: Effect of Burner Impedance,” ASME Paper No. GT2001-GT-0042.
Morse, P. M. , and Feshback, H. , 1953, Methods of Theoretical Physics, Vol. 1, McGraw-Hill, New York.
Schuermans, B. , Bellucci, V. , and Paschereit, C. O. , 2003, “Thermoacoustic Modeling and Control of Multi Burner Combustion Systems,” ASME Paper No. 2003-GT-38688.
Gray, R. M. , 2006, “Toeplitz and Circulant Matrices: A Review,” Found. Trends Commun. Inf. Theory, 2(3), pp. 155–239. [CrossRef]
Moeck, J. P. , Paul, M. , and Paschereit, C. O. , 2010, “Thermoacoustic Instabilities in an Annular Rijke Tube,” ASME Paper No. GT2010-23577.
Triantafyllou, M. , and Triantafyllou, G. , 1991, “Frequency Coalescence and Mode Localization Phenomena: A Geometric Theory,” J. Sound Vib., 150(3), pp. 485–500. [CrossRef]
Schuermans, B. , Polifke, W. , and Paschereit, C. O. , 1999, “Modeling Transfer Matrices of Premixed Flames and Comparison With Experimental Results,” ASME Paper No. 99-GT-132.
Goldmeer, J. , Vandervort, C. , and Sternberg, J. , 2017, “New Capabilities and Developments in GE's DLN 2.6 Combustion Systems,” Power-Gen International, Las Vegas, NV, Dec. 5–7, pp. 1–15.
Sewell, J. B. , and Sobieski, P. A. , 2005, “Monitoring of Combustion Instabilities: Calpine's Experience,” Combustion Instabilities in Gas Turbine Engines: Operational Experience, Fundamental Mechanisms, and Modeling (Progress in Astronautics and Aeronautics), T. C. Lieuwen and V. Yang , eds., American Institute of Aeronautics and Astronautics, Arlington, TX, pp. 147–162.
Pierre, C. , and Cha, D. , 1989, “Strong Mode Localization in Nearly Periodic Disordered Structures,” AIAA J., 27(2), pp. 227–241. [CrossRef]

Figures

Grahic Jump Location
Fig. 2

(a) Geometry of the 2D model of the transition duct, simplified from the original geometry of Fig. 1. (b) Legend, common to the other four figures. ((c) and (d)) First nontrivial eigenmode of the studied geometry for an axial mode, i.e., m = 0 (λ = 2L). Absolute value of the acoustic pressure on the left and absolute value of the acoustic velocity on the right. The solutions of the axial case match the solutions of the problem with Neumann conditions applied over the whole boundary. ((e) and (f)) same of (c) and (d) but for a push–pull mode, i.e., m = N/2 = 7, for a set of N = 14 cans (λ = 4.26L). In (f), the acoustic velocity has a strong increase where the gap starts and is inhomogeneous over the gap length.

Grahic Jump Location
Fig. 3

First eigenmode with azimuthal order m = 2, for a set of N = 14 communicating cans. Because the solution is complex-valued, we use different colormaps for the absolute value and the phase angle, presented in the two legends in (a). In (b), we present the solution in one can only: we observe that the phase is approximately constant on the cross section at the upstream end of the transition duct, on the left, far from the cross-talk area. The amplitude decreases slightly from left to right. In (c), we present the ensemble of all the 14 cans pulsating together for the mode in (b), reconstructed using Eq. (1). The phase changes mostly at the cross-talk area and changes two times (since the azimuthal order m is 2) the quantity 2π along the annulus. At the upstream end of each transition duct, the phase is approximately constant and the mode presents a certain azimuthal phase pattern. The phase between two cans is the difference of the phase value between two locations. Refer to Fig. 4(a) for experimental evidence of a second-order azimuthal mode like this one in an engine. (a) Legend for (b) and (c). (b) First azimuthal m = 2 mode in one can, pressure field. Legend in (a). (c) Phase in the 14 cans of the same mode of (a). The first can from the left is the same appearing in (b).

Grahic Jump Location
Fig. 4

Experimental evidence of a second azimuthal and of a push–pull mode. The timeseries of 14 pressure sensors in a 14-can can-annular combustor are processed at the frequency of one thermoacoustic instability. We present the amplitude of pulsation at that frequency on the left axis (in arbitrary units) and the phase difference on the right vertical axis at this frequency in each can. The phase difference is between the first can and the nth can, so that the value is zero in the first can. (a) and (b) refer to a different frequency of oscillation and a different operating condition. (a) The phase pattern corresponds to an azimuthal mode of order m = 2 that is rotating, because the phase changes with an approximately constant slope twice the amount of 2π moving along the annulus. This matches the theoretical prediction of Fig. 3(c). (b) The phase difference between most adjacent cans is very close to ±π (m = N/2 = 7). We observe some variation of the amplitude between cans, characterized as mode localization in The Effect of Asymmetry section of the paper.

Grahic Jump Location
Fig. 5

Phase of the equivalent reflection coefficient for the azimuthal modes for N = 14 cans, calculated with reference cross section at the red inlet of the transition duct of Fig. 2(a). The Helmholtz number is defined as He = /c.

Grahic Jump Location
Fig. 6

Validation of the geometrical simplifications for the 2D model for N = 12. The dashed lines are obtained from a fixed complex 3D geometry, while the continuous lines are obtained with the bidimensional model proposed in this paper, with the nondimensional values estimated from the 3D geometry without applying further corrections. The good agreement validates the geometrical approximations made to map the complex 3D geometry to the simpler 2D model.

Grahic Jump Location
Fig. 7

Sensitivity of the first two eigenfrequencies of a can system as a function of the geometry aspect ratio L/H, for all azimuthal wavenumber m, for N = 14. We plot the Helmholtz number of all the eigenmodes with He < 2π for each azimuthal mode m, as a function of the two nondimensional parameters governing can to can acoustic communication. As L/H increases, the Helmholtz number of the azimuthal modes converges slowly to π/2 and 3π/2 for the first two modes. For a fixed value of L/H, the eigenfrequencies for m = 5, 6, 7 are rather close, forming a cluster. These modes are closer together the larger is L/H.

Grahic Jump Location
Fig. 1

Geometry of a typical transition duct and the simplifications leading to the two-dimensional (2D) model, appearing ultimately in Fig. 2(a). The communicating gap between two adjacent transition ducts is colored in gray, with an axial length Lgap. (a) three-dimensional (3D) sketch of a transition duct, (b) front view of the air volume, and (c) side view of the air volume.

Grahic Jump Location
Fig. 8

Spectrogram of the timeseries of the azimuthal modes reconstructed from an engine. We can observe that more than one azimuthal mode is active in a small range of frequencies, forming a cluster. Within the cluster, each mode peaks at a frequency that increases with the azim. order m. The highest amplitude is not always in the same mode, and the system's energy appears to move randomly between the modes, with some statistical preference for the m = 3, 4, 5 modes.

Grahic Jump Location
Fig. 11

Complex eigenvalues for the three cases considered. The eigenvalues have been normalized with the angular eigenfrequency of the quarter-wave mode of an isolated can. The vertical axis is the normalized growth rate of oscillation of each mode. The red numbers refer to the asymmetric case and describe the mth azimuthal order of the respective symmetric eigenmode. The amplitude and phase of each asymmetric mode is presented in Fig. 12 in the plot with the same red number in the top left corner.

Grahic Jump Location
Fig. 10

Sketch of the nth can of the can-annular model. The can communicates with the others via the blue cross-talk area where the f, g waves are sketched.

Grahic Jump Location
Fig. 9

(a) Sketch of the reflection/transmission transfer functions T(d) for an incoming traveling wave f. Each function T(d)(ω) is the transmitted Riemann invariant from one can to a can that is d cans apart. For example, T(0) is the reflection of a wave propagating downstream of one can and being reflected to the same can, and T(3) can be interpreted for example as the transmission of a wave traveling downstream can 4 and propagating upstream in can 7, or any other two cans that are 3 cans apart. ((b) and (c)) Gain and phase of the transfer functions T(d) between all cans.

Grahic Jump Location
Fig. 12

Visualization of the azimuthal mode shapes associated with the eigenvalues of the asymmetric case (Fig. 11). Ordered from left to right and top to bottom with increasing oscillation frequency. The bar height indicates the pressure amplitude and the bar color the phase (same color legend as in Fig. 3), at the interface between each can and each transition duct. The number in the top left corner is the azimuthal order m of the corresponding mode in the symmetric case. Couple of modes with the same m differ very little in frequency and can result in a slowly modulated linear combination with approximately constant amplitude along the annulus. This is not the case for the perturbed push–pull mode in the bottom right corner, where cans 1, 2, 3, 10, 11, and 12 present much larger amplitude than the others.

Tables

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 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