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

Thermoacoustic Damping Rate Determination From Combustion Noise Using Bayesian Statistics

[+] Author and Article Information
Nicolai V. Stadlmair

Lehrstuhl für Thermodynamik,
Technische Universität München,
Garching 85748, Germany
e-mail: stadlmair@td.mw.tum.de

Tobias Hummel

Lehrstuhl für Thermodynamik,
Technische Universität München,
Garching 85748, Germany
e-mail: hummel@td.mw.tum.de

Thomas Sattelmayer

Lehrstuhl für Thermodynamik,
Technische Universität München,
Garching 85748, Germany
e-mail: sattelmayer@td.mw.tum.de

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 5, 2017; final manuscript received September 12, 2017; published online June 27, 2018. Editor: David Wisler.

J. Eng. Gas Turbines Power 140(11), 111501 (Jun 27, 2018) (7 pages) Paper No: GTP-17-1435; doi: 10.1115/1.4038475 History: Received August 05, 2017; Revised September 12, 2017

In this paper, we present a method to determine the quantitative stability level of a lean-premixed combustor from dynamic pressure data. Specifically, we make use of the autocorrelation function of the dynamic pressure signal acquired in a combustor where a turbulent flame acts as a thermoacoustic driver. In the proposed approach, the unfiltered pressure signal including several modes is analyzed by an algorithm based on Bayesian statistics. For this purpose, a Gibbs sampler is used to calculate parameters like damping rates and eigenfrequencies in the form of probability density functions (PDF) by a Markov-chain Monte Carlo (MCMC) method. The method provides a robust solution algorithm for fitting problems without requiring initial values. A further advantage lies in the nature of the statistical approach since the results can be assessed regarding its quality by means of the PDF and its standard deviation for each of the obtained parameters. First, a simulation of a stochastically forced van-der-Pol oscillator with preset input values is carried out to demonstrate accuracy and robustness of the method. In this context, it is shown that, despite a large amount of uncorrelated background noise, the identified damping rates are in a good agreement with the simulated parameters. Second, this technique is applied to measured pressure data. By doing so, the combustor is initially operated under stable conditions before the thermal power is gradually increased by adjusting the fuel mass flow rate until a limit-cycle oscillation is established. It is found that the obtained damping rates are qualitatively in line with the amplitude levels observed during operation of the combustor.

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


Rayleigh, J. L., 1878, “ The Explanation of Certain Acoustical Phenomena,” Nature, 18(4), pp. 319–321.
Lieuwen, T. C. , 2012, Unsteady Combustor Physics, Cambridge University Press, Cambridge, UK. [CrossRef]
Wieber, P. R. , 1966, “ Acoustic Decay Coefficients of Simulated Rocket Combustors,” NASA Lewis Research Center, Cleveland, OH, Technical Report No. NASA TN D-3425. https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19660014565.pdf
Phillips, B. , and Morgan, C. J. , 1967, “ Mechanical Absorption of Acoustic Oscillations in Simulated Rocket Combustion Chambers,” NASA Lewis Research Center, Cleveland, OH, Technical Report No. NASA TN D-3792. https://ntrs.nasa.gov/search.jsp?R=19670007361
Sattelmayer, T. , Kathan, R. , Köglmeier, S. , Kaess, R. , and Nicole, A. , 2015, “ Validation of Transverse Instability Damping Computations for Rocket Engines,” J. Propul. Power, 31(4), pp. 1148–1158. [CrossRef]
Kathan, R. , 2013, “ Verlustmechanismen in Raketenbrennkammern,” Ph.D. thesis, Technische Universität München, Munich, Germany.
Fiala, T. , Kathan, R. , and Sattelmayer, T. , 2011, “ Effective Stability Analysis of Liquid Rocket Combustion Chambers: Experimental Investigation of Damped Admittances,” 62nd International Astronautical Congress, Cape Town, South Africa, Oct. 3–7, Paper No. IAC11-C4.3.11. https://www.researchgate.net/publication/264706972_Effective_Stability_Analysis_of_Liquid_Rocket_Combustion_Chambers_Experimental_Investigation_of_Damped_Admittances
Bodisco, T. , Reeves, R. , Situ, R. , and Brown, R. , 2012, “ Bayesian Models for the Determination of Resonant Frequencies in a DI Diesel Engine,” Mech. Syst. Signal Process., 26, pp. 305–314. [CrossRef]
Wagner, M. , Jörg, C. , and Sattelmayer, T. , 2013, “ Comparison of the Accuracy of Time-Domain Measurement Methods for Combustor Damping,” ASME Paper No. GT2013-94844.
Stadlmair, N. V. , Wagner, M. , and Hirsch, C. , 2015, “ Experimentally Determining the Acoustic Damping Rates of a Combustor With a Swirl Stabilized Lean Premixed Flame,” ASME Paper No. GT2015-42683.
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]
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]
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(1–2), pp. 442–447. [CrossRef]
Nair, V. , and Sujith, R. I. , 2014, “ Multifractality in Combustion Noise: Predicting an Impending Combustion Instability,” J. Fluid Mech., 747, pp. 635–655. [CrossRef]
Lieuwen, T. , 2005, “ Online Combustor Stability Margin Assessment Using Dynamic Pressure Data,” ASME J. Eng. Gas Turbines Power, 127(3), pp. 478–482. [CrossRef]
Lieuwen, T. C. , 2005, “ Systems and Methods for Detection of Combustor Stability Margin,” U.S. Patent No. US20050247064 A1. https://www.google.com/patents/US20050247064
Noiray, N. , and Schuermans, B. , 2013, “ Deterministic Quantities Characterizing Noise Driven Hopf Bifurcations in Gas Turbine Combustors,” Int. J. Non-Linear Mech., 50, pp. 152–163. [CrossRef]
Noiray, N. , and Schuermans, B. A. , 2013, “ On the Dynamic Nature of Azimuthal Thermoacoustic Modes in Annular Gas Turbine Combustion Chambers,” Proc. R. Soc. A, 469(2151), p. 20120535.
Mejia, 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]
Geman, S. , and Geman, D. , 1984, “ Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images,” IEEE Trans. Pattern Anal. Mach. Intell., 6(6), pp. 721–741. [CrossRef] [PubMed]
Metropolis, N. , Rosenbluth, A. W. , Rosenbluth, M. N. , Teller, A. H. , and Teller, E. , 1953, “ Equation of State Calculations by Fast Computing Machines,” J. Chem. Phys., 21(6), pp. 1087–1092. [CrossRef]
Hastings, W. K. , 1970, “ Monte Carlo Sampling Methods Using Markov Chains and Their Applications,” Biometrika, 57(1), pp. 97–109. https://academic.oup.com/biomet/article-abstract/57/1/97/2721936
Culick, F. E. , 2006, “ Unsteady Motions in Combustion Chambers for Propulsion Systems,” NATO Research and Technology Organization, Neuilly-sur-Seine, France, Technical Report No. AG-AVT-039. http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADA466461
Gardiner, C. , 2009, Stochastic Methods, Springer, Berlin.
Ahlfors, L. V. , 1966, Complex Analysis: An Introduction to the Theory of Analytic Functions of One Complex Variable, McGraw-Hill, New York.
Elster, C. , and Wübbeler, G. , 2016, “ Bayesian Regression versus Application of Least Squares—An Example,” Metrologia, 53(1), p. S10.
Robert, C. P. , and Casella, G. , 2004, Monte Carlo Statistical Methods (Springer Texts in Statistics), Vol. 95, Springer, New York. [CrossRef]
Rubinstein, R. Y. , and Kroese, D. P. , 2011, Simulation and the Monte Carlo Method, Wiley, Hoboken, NJ.
Plummer, M. , 2003, “ JAGS: A Program for Analysis of Bayesian Graphical Models Using Gibbs Sampling,” Third International Workshop on Distributed Statistical Computing (DSC), Vienna, Austria, Mar. 20–22, p. 125. https://www.r-project.org/conferences/DSC-2003/Proceedings/Plummer.pdf
Noiray, N. , 2016, “ Linear Growth Rate Estimation From Dynamics and Statistics of Acoustic Signal Envelope in Turbulent Combustors,” ASME Paper No. GT2016-58169.
Lieuwen, T. C. , 2003, “ Statistical Characteristics of Pressure Oscillations in a Premixed Combustor,” J. Sound Vib., 260(1), pp. 3–17. [CrossRef]
Mordechai, S. , 2011, Applications of Monte Carlo Method in Science and Engineering, InTech, Rijeka, Croatia. [CrossRef]
Noiray, N. , and Denisov, A. , 2017, “ A Method to Identify Thermoacoustic Growth Rates in Combustion Chambers From Dynamic Pressure Time Series,” Proc. Combust. Inst., 36(3), pp. 3843–3850. [CrossRef]
Stadlmair, N. V. , Mohammadzadeh Keleshtery, P. , Zahn, M. , and Sattelmayer, T. , 2017, “ Impact of Water Injection on Thermoacoustic Modes in a Lean Premixed Combustor Under Atmospheric Conditions,” ASME Paper No. GT2017-63342.
Betz, M. , Zahn, M. , Wagner, M. , Stadlmair, N. V. , Schulze, M. , Hirsch, C. , and Sattelmayer, T. , 2017, “ Impact of Damper Parameters on the Stability Margin of an Annular Combustor Test Rig,” ASME Paper No. GT2017-64239.
Hummel, T. , Berger, F. M. , Stadlmair, N. V. , and Schuermans, B. , 2017, “ Extraction of Linear Growth and Damping Rates of High-Frequency Thermoacoustic Oscillations From Time Domain Data,” ASME Paper No. GT2017-64233.


Grahic Jump Location
Fig. 1

Block diagram of the second-order oscillator system to generate synthetic time series

Grahic Jump Location
Fig. 2

Numerical validation test cases for single mode system (a) and multimode system (b): time series of raw signal, probability density of pressure signal p′, normalized PSD of p′, autocorrelation and envelope of the autocorrelation

Grahic Jump Location
Fig. 3

Numerical validation test case A: probability densities of the fitting parameters ν1 and f=ω/2π and its corresponding standard deviations σ(…), autocorrelation function of the signal, fitted autocorrelation

Grahic Jump Location
Fig. 4

Comparison of different simulation test cases A and B: autocorrelation of the signal, fitted autocorrelation, identified parameters ν1, ν2, and simulation parameters

Grahic Jump Location
Fig. 5

Evolution of the error for the damping rates identifying one mode I = 1 and two modes I = 2

Grahic Jump Location
Fig. 6

Atmospheric single burner test-rig with lean-premixed, swirl-stabilized flame

Grahic Jump Location
Fig. 7

Signal characteristics of the measurement data for ϕ=0.63, Pth = 41.80 kW: time series of raw signal, probability density of pressure signal p′, normalized PSD of p′, autocorrelation of the signal and envelope of the autocorrelation function

Grahic Jump Location
Fig. 8

Evolution of the damping rate (top) ν1 and ν2 and normalized modal amplitudes (bottom) η1 and η2 in the combustor for varying equivalence ratios



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