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

Linear Growth Rate Estimation From Dynamics and Statistics of Acoustic Signal Envelope in Turbulent Combustors

[+] Author and Article Information
Nicolas Noiray

CAPS Laboratory,
Mechanical and Process
Engineering Department,
ETH Zürich,
Zürich 8092, Switzerland
e-mail: noirayn@ethz.ch

Contributed by the Combustion and Fuels Committee of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received July 19, 2016; final manuscript received July 21, 2016; published online October 18, 2016. Editor: David Wisler.

J. Eng. Gas Turbines Power 139(4), 041503 (Oct 18, 2016) (11 pages) Paper No: GTP-16-1350; doi: 10.1115/1.4034601 History: Received July 19, 2016; Revised July 21, 2016

Considerable research and development efforts are required to meet the targets of future gas turbine technologies in terms of performance, emissions, and operational flexibility. One of the recurring problems is the constructive coupling between flames and combustor's acoustics. These thermoacoustic interactions can cause high-amplitude dynamic pressure limit cycles, which reduce the lifetime of the hot gas path parts or in the worst-case scenario destroy these mechanical components as a result of a sudden catastrophic event. It is shown in this paper that the dynamics and the statistics of the acoustic signal envelope can be used to identify the linear growth rates hidden behind the observed pulsations, and the results are validated against numerical simulations. This is a major step forward and it will contribute to the development of future gas turbine combustors, because the knowledge of these linear growth rates is essential to develop robust active and passive systems to control these combustion instabilities.

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

References

Lieuwen, T. , 2012, Unsteady Combustor Physics, Cambridge University Press, New York.
Noiray, N. , and Schuermans, B. , 2013, “ Deterministic Quantities Characterizing Noise Driven Hopf Bifurcations in Gas Turbine Combustor,” Int. J. Non-Linear Mech., 50, pp. 152–163. [CrossRef]
Lieuwen, T. C. , 2003, “ Statistical Characteristics of Pressure Oscillations in a Premixed Combustor,” J. Sound Vib., 260(1), pp. 3–17. [CrossRef]
Seume, J. R. , Vortmeyer, N. , Krause, W. , Hermann, J. , Hantschk, C. , Zangl, P. , Gleis, S. , Vortmeyer, D. , and Orthmann, A. , 1998, “ Application of Active Combustion Instability Control to a Heavy Duty Gas Turbine,” ASME J. Eng. Gas Turbines Power, 120(4), pp. 721–726. [CrossRef]
Noiray, N. , and Schuermans, B. , 2012, “ Theoretical and Experimental Investigations on Damper Performance for Suppression of Thermoacoustic Oscillations,” J. Sound Vib., 331(12), pp. 2753–2763. [CrossRef]
Bothien, M. , Moeck, J. , and Paschereit, C. O. , 2010, “ Comparison of Linear Stability Analysis With Experiments by Actively Tuning the Acoustic Boundary Conditions of a Premixed Combustor,” ASME J. Eng. Gas Turbines Power, 132(12), p. 121502. [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. , and Juniper, M. P. , 2013, “ Azimuthal Instabilities in Annular Combustors: Standing and Spinning Modes,” Proc. R. Soc. A, 469(2157), p. 20130232. [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]
Morse, P. M. , and Ingard, K. U. , 1968, Theoretical Acoustics, McGraw-Hill, New York.
Culick, F. E. C. , 2006, “ Unsteady Motions in Combustion Chambers for Propulsion Systems (AGARDograph),” OrganisationNorth Atlantic Treaty OrganisationNeuilly-sur-Seine Cedex, France Report No. NATO/RTO-AG-AVT-039.
Ghirardo, G. , Juniper, M. , and Moeck, J. , 2015, “ Stability Criteria for Standing and Spinning Waves in Annular Combustors,” ASME Paper No. GT2015-43127.
Poinsot, T. , and Veynante, D. , 2001, Theoretical and Numerical Combustion, Edwards, Philadelphia, PA.
Clavin, P. , Kim, J. S. , and Williams, F. A. , 1994, “ Turbulence-Induced Noise Effects on High-Frequency Combustion Instabilities,” Combust. Sci. Technol., 96, pp. 61–84. [CrossRef]
Lieuwen, T. , and Banaszuk, A. , 2005, “ Background Noise Effects on Combustor Stability,” J. Propul. Power, 21(1), pp. 25–31. [CrossRef]
Balachandran, R. , Ayoola, B. O. , Kaminski, C. F. , Dowling, A. P. , and Mastorakos, E. , 2005, “ Experimental Investigation of the Nonlinear Response of Turbulent Premixed Flames to Imposed Inlet Velocity Oscillations,” Combust. Flame, 143, pp. 37–55. [CrossRef]
Krylov, N. , and Bogoliubov, N. , 1949, Introduction to Nonlinear Mechanics, Princeton University Press, London, UK.
Lieuwen, T. , 2005, “ Online Combustor Stability Margin Assessment Using Dynamic Pressure Data,” ASME J. Eng. Gas Turbines Power, 127(3), pp. 478–482. [CrossRef]
Yi, T. , and Gutmark, E. J. , 2008, “ Online Prediction of the Onset of Combustion Instability Based on the Computation of Damping Ratios,” J. Sound Vib., 310, pp. 442–447. [CrossRef]
Schulze, M. , Kathan, R. , and Sattelmayer, T. , 2015, “ Impact of Absorber Ring Position and Cavity Length on Acoustic Damping,” J. Spacecr. Rockets, 52(3), pp. 917–927. [CrossRef]
Hagen, G. , 2011, “ Stochastic Averaging for Identification of Feedback Nonlinearities in Thermoacoustic Systems,” ASME J. Dyn. Syst., Meas., Control, 133(6), p. 061017. [CrossRef]
Meija, D. , Miguel-Brebion, M. , and Selle, L. , 2016, “ On the Experimental Determination of Growth and Damping Rates for Combustion Instabilities,” Combust. Flame, 169, pp. 287–296. [CrossRef]
Lax, M. , 1967, “ Classical Noise. V. Noise in Self-Sustained Oscillators,” Phys. Rev., 160(2), pp. 290–307. [CrossRef]
Lax, M. , Cai, W. , and Xu, M. , 2006, Random Processes in Physics and Finance, Oxford University Press, New York.
Friedrich, R. , Siegert, S. , Peinke, J. , Lück, S. , Siefert, M. , Lindemann, M. , Raethjen, J. , Deuschl, G. , and Pfister, G. , 2000, “ Extracting Model Equations From Experimental Data,” Phys. Lett. A, 271(3), pp. 217–222. [CrossRef]
Gradisek, J. , Siegert, S. , Friedrich, R. , and Grabec, I. , 2000, “ Analysis of Time Series From Stochastic Processes,” Phys. Rev. E, 62(3146–3155), pp. 217–222.
Kolmogorov, A. N. , 1931, “ Über die Analytischen Methoden in der Wahrscheinlichkeitsrechnung,” Math. Ann., 140, pp. 415–458. [CrossRef]
Friedrich, R. , Peinke, J. , Sahimi, M. , and Tabar, M. R. R. , 2011, “ Approaching Complexity by Stochastic Methods: From Biological Systems to Turbulence,” Phys. Rep., 506(5), pp. 87–162. [CrossRef]
Siegert, S. , Friedrich, R. , and Peinke, J. , 1998, “ Analysis of Data Sets of Stochastic Systems,” Phys. Lett. A, 243, pp. 275–280. [CrossRef]
Stratonovich, R. C. , 1963, Topics in the Theory of Random Noise, Vol. 2, Gordon & Breach, New York.
Roberts, J. B. , and Spanos, P. D. , 1986, “ Invited Review No. 1: Stochastic Averaging: An Approximate Method of Solving Random Vibration Problems,” Int. J. Non-Linear Mech., 21(2), pp. 111–134. [CrossRef]
Kim, K. , and Hochgreb, S. , 2012, “ Measurements of Triggering and Transient Growth in a Model Lean-Premixed Gas Turbine Combustor,” Combust. Flame, 159(3), pp. 1215–1227. [CrossRef]
Moeck, J. , Bothien, M. , Schimek, S. , Lacarelle, A. , and Paschereit, C. O. , 2008, “ Subcritical Thermoacoustic Instabilities in a Premixed Combustor,” AIAA Paper No. 2008-2946.
Ushakov, O. V. , Wunsche, H.-J. , Henneberger, F. , Khovanov, I. A. , Schimansky-Geier, L. , and Zaks, M. A. , 2005, “ Coherence Resonance Near a Hopf Bifurcation,” Phys. Rev. Lett., 95(12), p. 123903. [CrossRef] [PubMed]
Zakharova, A. , Vadivasova, T. , Anishchenko, V. , Koseska, A. , and Kurths, J. , 2010, “ Stochastic Bifurcations and Coherencelike Resonance in a Self-Sustained Bistable Noisy Oscillator,” Phys. Rev. E, 81(1), p. 011106. [CrossRef]
Kabiraj, L. , Steinert, R. , Saurahb, A. , and Paschereit, C. O. , 2015, “ Coherence Resonance in a Thermoacoustic System,” Phys. Rev. E, 92(4), p. 042909. [CrossRef]
Waugh, I. C. , and Juniper, M. P. , 2011, “ Triggering in a Thermoacoustic System With Stochastic Noise,” Int. J. Spray Combust. Dyn., 3(3–4), pp. 225–242. [CrossRef]
Gotoda, H. , Shinoda, Y. , Kobayashi, M. , Okuno, Y. , and Tachibana, S. , 2014, “ Detection and Control of Combustion Instability Based on the Concept of Dynamical System Theory,” Phys. Rev. E, 89(2), p. 022910. [CrossRef]
Gopalakrishnan, E. A. , and Sujith, R. I. , 2015, “ Effect of External Noise on the Hysteresis Characteristics of a Thermoacoustic System,” J. Fluid Mech., 776, pp. 334–353. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Pulsation data from a heavy-duty gas turbine combustor at three different conditions: (a) nominal condition, (b) off-design operating points where dynamic pressure amplitude exceeds allowable limits (figure taken from Ref. [2]). Top row: joint probability density functions of modal amplitude and its time derivative. Bottom row: time traces of modal amplitude η, which is proportional to the acoustic pressure p.

Grahic Jump Location
Fig. 2

Nonlinear response of the normalized heat release rate coherent fluctuations qc as function of the amplitude η of the monochromatic acoustic forcing for a turbulent premixed flame. (∘) experimental measurements, adapted from Fig. 13(b) in Ref. [16], where qc is assumed proportional to the CH* chemiluminescence.

Grahic Jump Location
Fig. 3

Postprocessing of the signals from time domain simulations (60 s) of the stochastic Van der Pol oscillator given by Eqs. (9) and (11), for different linear growth rates ν (the signals are normalized by the maximum amplitude reached during the simulation). For these simulations, β = 60 rad/s, κ=3.4, Γ=107, and ωj=2πfj with fj = 120 Hz. Three values of the damping α are considered (72, 60, and 48 rad/s). These simulations, respectively, correspond to a linearly stable thermoacoustic system (top row), a marginally stable one when damping α balances the linear driving from the flames β (middle row), and a stochastic limit-cycle when the growth rate ν=(β−α)/2 is positive (bottom row). On each rows, from left to right: the probability density function of the simulated acoustic signal P(η), a small portion of the time series and its envelope (black curves), the probability density function of this envelope P(A), and the power spectral densities of the envelope SAA and of the acoustic pressure Sηη.

Grahic Jump Location
Fig. 4

(a) Stochastic bifurcation diagram associated with the acoustic amplitude dynamics. The bifurcation parameter is the growth rate ν that has been varied by setting different values of the damping α. A contour plot of the theoretical probability density function (p.d.f. P(A)) given by Eq. (25) is shown together with the deterministic bifurcation curve (solid line). (e)–(g) P(A) at selected linear growth rates νi with a variation from linearly stable (ν<0) to linearly unstable (ν>0) conditions. The shaded areas are the theoretical p.d.f., and the thick blue lines correspond to the p.d.f. from the normalized histograms of the simulated stochastic Van der Pol oscillators (see Figs. 3(c), 3(h), and 3(m)). (b)–(d) Corresponding potential and related contributions.

Grahic Jump Location
Fig. 5

Evolution of half the identified cutoff pulsation ωc/2 (green line), which is associated with the power spectral density SAA of the envelope of the acoustic pressure signal, as a function of the linear growth rate. The cutoff pulsation has been estimated from time-domain simulations of a stochastically forced Van der Pol oscillator. The theoretical value ωc/2=|ν| deduced from the linearized model is also plotted (red line). In the insets, one can see two example of SAA with the estimated and theoretical values of ωc.

Grahic Jump Location
Fig. 6

System identification using Eq. (35) to process the data statistics from the time domain simulations (symbols). This is done in order to estimate F(A), which can be afterward used to deduce the potential V(A). These data correspond to the ones presented in Fig. 4. The dashed lines in the upper diagrams correspond to a parametric system identification using the Van der Pol approximation, for which the identified ν, κ, and Γ give the best curve fit on the extracted amplitude-dependent transition moments (symbols colored by the probability density function of the amplitude). The theoretical curves are shown as black solid lines. The corresponding probability density functions are also given in the second row.

Grahic Jump Location
Fig. 7

Postprocessing of the signals from time domain simulations (240 s) of the self-sustained oscillator given by Eqs. (9) and (B1), for different linear growth rates ν. Same parameters as the Van der Pol case presented in Fig. 3.

Grahic Jump Location
Fig. 8

(a) Stochastic bifurcation diagram of the acoustic amplitude dynamics for an arctangent saturation. The bifurcation parameter is the growth rate ν that has been varied by setting different values of the damping α (same numerical values as for the Van der Pol case). A contour plot of the theoretical p.d.f. P(A) given by Eqs. (25) and (B4) is shown together with the deterministic bifurcation curve (solid line). (e)–(g) P(A) at selected linear growth rates νi. The shaded areas are the theoretical p.d.f., and the thick blue lines correspond to the p.d.f. from the normalized histograms of the simulated model with arctangent saturation (see Fig. 7). (b)–(d) Corresponding potential from Eq. (B4).

Grahic Jump Location
Fig. 9

Nonlinear response of the normalized heat release rate coherent fluctuations qc as function of the amplitude η of the monochromatic acoustic forcing for a turbulent premixed flame. (∘) experimental measurements of the driving, adapted from Fig. 6 in Ref. [32].

Grahic Jump Location
Fig. 10

Postprocessing of the signals from time domain simulations (240 s) of the self-sustained oscillator resulting from a coherent heat release rate feedback containing linear, cubic, and quintic terms (see Fig. 9). The arbitrarily chosen numerical values of the corresponding coefficients are ωj=2π×120 rad/s, β=60 rad/s, κ=9, γ=0.7, Γ=107, and three values for α: from top to bottom, 80 rad/s in the subthreshold region, and α: 72 and 68 rad/s in the bistable region.

Grahic Jump Location
Fig. 11

(a) Stochastic bifurcation diagram of the acoustic amplitude A for the subcritical bifurcation obtained using the nonlinear heat release approximation of Figs. 9 and 10. The bifurcation parameter is the growth rate ν that has been varied by setting different values of the damping α. A contour plot of the theoretical p.d.f. P(A) deduced from the potential (C2) is shown together with the deterministic bifurcation curve (solid line). (e)–(g) P(A) at selected linear growth rates νi corresponding toFig. 10. The shaded areas are the theoretical p.d.f., and the thick blue lines correspond to the p.d.f. from the normalized histograms of the simulations. (b)–(d) Corresponding potential from Eq. (C2).

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