Thermoacoustic instability is a persistent problem in rocket engines and gas turbines because oscillations are so sensitive to many design parameters. This, however, presents an excellent opportunity for the application of adjoint methods and the combination of experiments and modelling.

The first scientific reports of thermoacoustic oscillations appeared over two centuries ago. Their practical consequences have been evident since liquid rocket engine development in the 1930's: they cause thrust oscillations, structural damage, increased heat transfer, and component or payload failure. Despite decades of research by Germany from the 1930's, by the USA and USSR during the cold war, and recently by the gas turbine industry, these oscillations remain a severe problem today.

In order to achieve high power to weight ratios and high efficiency, rocket and gas turbine engines all have low acoustic damping and high energy densities: up to 50 GW/m3 for liquid rockets, 1 GW/m3 for solid rockets, and 0.1 GW/m3 for jet engines and afterburners. Consequently, large amplitude oscillations can be sustained even if the thermoacoustic mechanism is only slightly (i.e. ~0.1%) efficient for at least one acoustic mode.

Engine development typically consists of component tests, sector tests, full combustor tests, and full engine tests. Thermoacoustic instability tends to recur during the later stages and is rarely predicted reliably by component tests and analysis. This is because the efficiency of the thermoacoustic mechanism is exceedingly sensitive to small changes to the system for reasons described in the following review paper.

Sensitivity and nonlinearity in Thermoacoustics

M. Juniper, R. I. Sujith

*Annual Review of Fluid Mechanics* **50** 661--689 (2018) doi:10.1146/annurev-fluid-122316-045125

M. Juniper, R. I. Sujith

pdf

Open Access

doi: https://doi.org/10.1146/annurev-fluid-122316-045125

Tutorial 1: Obtaining thermoacoustic eigenvalue sensitivities with adjoint methods

Tutorial 1 Matlab files

Tutorial 2: Tools from nonlinear dynamics

Tutorial 2 Matlab files

Nine decades of rocket engine and gas turbine development have shown that thermoacoustic oscillations are difficult to predict but can usually be eliminated with relatively small ad hoc design changes. These changes can, however, be ruinously expensive to devise. This review explains why linear and nonlinear thermoacoustic behaviour is so sensitive to parameters such as operating point, fuel composition, and injector geometry. It shows how non-periodic behaviour arises in experiments and simulations and discusses how fluctuations in thermoacoustic systems with turbulent reacting flow, which are usually filtered or averaged out as noise, can reveal useful information. Finally, it proposes tools to exploit this sensitivity in the future: adjoint-based sensitivity analysis to optimize passive control designs, and complex systems theory to warn of impending thermoacoustic oscillations and to identify the most sensitive elements of a thermoacoustic system.

Open Access

doi: https://doi.org/10.1146/annurev-fluid-122316-045125

Tutorial 1: Obtaining thermoacoustic eigenvalue sensitivities with adjoint methods

Tutorial 1 Matlab files

Tutorial 2: Tools from nonlinear dynamics

Tutorial 2 Matlab files

Nine decades of rocket engine and gas turbine development have shown that thermoacoustic oscillations are difficult to predict but can usually be eliminated with relatively small ad hoc design changes. These changes can, however, be ruinously expensive to devise. This review explains why linear and nonlinear thermoacoustic behaviour is so sensitive to parameters such as operating point, fuel composition, and injector geometry. It shows how non-periodic behaviour arises in experiments and simulations and discusses how fluctuations in thermoacoustic systems with turbulent reacting flow, which are usually filtered or averaged out as noise, can reveal useful information. Finally, it proposes tools to exploit this sensitivity in the future: adjoint-based sensitivity analysis to optimize passive control designs, and complex systems theory to warn of impending thermoacoustic oscillations and to identify the most sensitive elements of a thermoacoustic system.

This extreme sensitivity allows thermoacoustic systems to be stabilized with small design changes, usually at the full engine test stage. These changes can, however, be ruinously expensive to devise. The challenge is to devise and implement these design changes cheaply, quickly, and earlier in the design process.

There are many ways to model a thermoacoustic system. Acoustic network models are widely used in industry. Typically there will be many stable modes and just a handful of unstable modes:

We wish to stabilize the handful of unstable modes by altering the design parameters. There may be several thousand of these and it is impractical to work out the influence of each one by varying each in turn. With adjoint methods, however, we can evaluate the influence of every design parameter on a given eigenvalue with a single calculation. Therefore we need only as many calculations as there are unstable eigenvalues.

Our first application of this technique was to a simple model of a hot wire in a tube and the second was to a Burke-Schumann flame in a tube:

Sensitivity analysis of a time-delayed thermoacoustic system via an adjoint-based approach

L. Magri and M. P. Juniper

*Journal of Fluid Mechanics* **719** 183--202 (2013) doi:10.1017/jfm.2012.639

L. Magri and M. P. Juniper

pdf

doi: https://doi.org/10.1017/jfm.2012.639

We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat-release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release time-delay coefficient, the heat-release parameter, and the hot-wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates.

doi: https://doi.org/10.1017/jfm.2012.639

We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat-release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release time-delay coefficient, the heat-release parameter, and the hot-wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates.

Global modes, receptivity, and sensitivity analysis of diffusion flames coupled with duct acoustics

L. Magri, M. P. Juniper

*Journal of Fluid Mechanics* **752** 237--265 (2014) doi:10.1017/jfm.2014.328

L. Magri, M. P. Juniper

pdf

Open Access

doi: https://doi.org/10.1017/jfm.2014.328

In this theoretical and numerical paper, we derive the adjoint equations for a thermo-acoustic system consisting of an infinite-rate chemistry diffusion flame coupled with duct acoustics. We then calculate the thermo-acoustic system's linear global modes (i.e. the frequency/growth rate of oscillations, together with their mode shapes), and the global modes' receptivity to species injection, sensitivity to base-state perturbations and structural sensitivity to advective-velocity perturbations. Some of these could be found by finite difference calculations but the adjoint analysis is computationally much cheaper. We then compare these with the Rayleigh index. The receptivity analysis shows the regions of the flame where open-loop injection of fuel or oxidizer will have the greatest influence on the thermo-acoustic oscillation. We find that the flame is most receptive at its tip. The base-state sensitivity analysis shows the influence of each parameter on the frequency/growth rate. We find that perturbations to the stoichiometric mixture fraction, the fuel slot width and the heat-release parameter have most influence, while perturbations to the Peclet number have the least influence for most of the operating points considered. These sensitivities oscillate, e.g. positive perturbations to the fuel slot width either stabilizes or destabilizes the system, depending on the operating point. This analysis reveals that, as expected from a simple model, the phase delay between velocity and heat-release fluctuations is the key parameter in determining the sensitivities. It also reveals that this thermo-acoustic system is exceedingly sensitive to changes in the base state. The structural-sensitivity analysis shows the influence of perturbations to the advective flame velocity. The regions of highest sensitivity are around the stoichiometric line close to the inlet, showing where velocity models need to be most accurate. This analysis can be extended to more accurate models and is a promising new tool for the analysis and control of thermo-acoustic oscillations.

Open Access

doi: https://doi.org/10.1017/jfm.2014.328

In this theoretical and numerical paper, we derive the adjoint equations for a thermo-acoustic system consisting of an infinite-rate chemistry diffusion flame coupled with duct acoustics. We then calculate the thermo-acoustic system's linear global modes (i.e. the frequency/growth rate of oscillations, together with their mode shapes), and the global modes' receptivity to species injection, sensitivity to base-state perturbations and structural sensitivity to advective-velocity perturbations. Some of these could be found by finite difference calculations but the adjoint analysis is computationally much cheaper. We then compare these with the Rayleigh index. The receptivity analysis shows the regions of the flame where open-loop injection of fuel or oxidizer will have the greatest influence on the thermo-acoustic oscillation. We find that the flame is most receptive at its tip. The base-state sensitivity analysis shows the influence of each parameter on the frequency/growth rate. We find that perturbations to the stoichiometric mixture fraction, the fuel slot width and the heat-release parameter have most influence, while perturbations to the Peclet number have the least influence for most of the operating points considered. These sensitivities oscillate, e.g. positive perturbations to the fuel slot width either stabilizes or destabilizes the system, depending on the operating point. This analysis reveals that, as expected from a simple model, the phase delay between velocity and heat-release fluctuations is the key parameter in determining the sensitivities. It also reveals that this thermo-acoustic system is exceedingly sensitive to changes in the base state. The structural-sensitivity analysis shows the influence of perturbations to the advective flame velocity. The regions of highest sensitivity are around the stoichiometric line close to the inlet, showing where velocity models need to be most accurate. This analysis can be extended to more accurate models and is a promising new tool for the analysis and control of thermo-acoustic oscillations.

Once the technique had been proven, we were able to move to thermoacoustic network models, which are used in industry. Using adjoint methods, we showed how to stabilize Rama Balachandran's thermoacoustic system by making small design changes:

Adjoint-based sensitivity analysis of low-order thermoacoustic networks using a wave-based approach

J. G. Aguilar, L. Magri, M. P. Juniper

*Journal of Computational Physics* **341** 163--181 (2017) doi:10.1016/j.jcp.2017.04.013

J. G. Aguilar, L. Magri, M. P. Juniper

pdf

Open Access

doi: https://doi.org/10.1016/j.jcp.2017.04.013

Strict pollutant emission regulations are pushing gas turbine manufacturers to develop devices that operate in lean conditions, with the downside that combustion instabilities are more likely to occur. Methods to predict and control unstable modes inside combustion chambers have been developed in the last decades but, in some cases, they are computationally expensive. Sensitivity analysis aided by adjoint methods provides valuable sensitivity information at a low computational cost. This paper introduces adjoint methods and their application in wave-based low order network models, which are used as industrial tools, to predict and control thermoacoustic oscillations. Two thermoacoustic models of interest are analysed. First, in the zero Mach number limit, a nonlinear eigenvalue problem is derived, and continuous and discrete adjoint methods are used to obtain the sensitivities of the system to small modifications. Sensitivities to base-state modification and feedback devices are presented. Second, a more general case with non-zero Mach number, a moving flame front and choked outlet, is presented. The influence of the entropy waves on the computed sensitivities is shown.

Open Access

doi: https://doi.org/10.1016/j.jcp.2017.04.013

Strict pollutant emission regulations are pushing gas turbine manufacturers to develop devices that operate in lean conditions, with the downside that combustion instabilities are more likely to occur. Methods to predict and control unstable modes inside combustion chambers have been developed in the last decades but, in some cases, they are computationally expensive. Sensitivity analysis aided by adjoint methods provides valuable sensitivity information at a low computational cost. This paper introduces adjoint methods and their application in wave-based low order network models, which are used as industrial tools, to predict and control thermoacoustic oscillations. Two thermoacoustic models of interest are analysed. First, in the zero Mach number limit, a nonlinear eigenvalue problem is derived, and continuous and discrete adjoint methods are used to obtain the sensitivities of the system to small modifications. Sensitivities to base-state modification and feedback devices are presented. Second, a more general case with non-zero Mach number, a moving flame front and choked outlet, is presented. The influence of the entropy waves on the computed sensitivities is shown.

The application to annular combustors is technically more complicated because they contain degenerate eigenvalues and a second order analysis is required:

Stability analysis of thermoacoustic nonlinear eigenproblems in annular combustors. Part 1: sensitivity

L. Magri, M. Bauerheim, M. P. Juniper

*Journal of Computational Physics* **325** 395--410 (2016) doi:10.1016/j.jcp.2016.07.032

L. Magri, M. Bauerheim, M. P. Juniper

pdf

Open Access

doi: https://doi.org/10.1016/j.jcp.2016.07.032

We present an adjoint-based method for the calculation of eigenvalue perturbations in nonlinear, degenerate and non self-adjoint eigenproblems. This method is applied to a thermo-acoustic annular combustor network, the stability of which is governed by a nonlinear eigenproblem. We calculate the first- and second-order sensitivities of the growth rate and frequency to geometric, flow and flame parameters. Three different configurations are analysed. The benchmark sensitivities are obtained by finite difference, which involves solving the nonlinear eigenproblem at least as many times as the number of parameters. By solving only one adjoint eigenproblem, we obtain the sensitivities to any thermo-acoustic parameter, which match the finite-difference solutions at much lower computational cost.

Open Access

doi: https://doi.org/10.1016/j.jcp.2016.07.032

We present an adjoint-based method for the calculation of eigenvalue perturbations in nonlinear, degenerate and non self-adjoint eigenproblems. This method is applied to a thermo-acoustic annular combustor network, the stability of which is governed by a nonlinear eigenproblem. We calculate the first- and second-order sensitivities of the growth rate and frequency to geometric, flow and flame parameters. Three different configurations are analysed. The benchmark sensitivities are obtained by finite difference, which involves solving the nonlinear eigenproblem at least as many times as the number of parameters. By solving only one adjoint eigenproblem, we obtain the sensitivities to any thermo-acoustic parameter, which match the finite-difference solutions at much lower computational cost.

Thermoacoustic Helmholtz solvers solve the wave equation in an arbitrarily complex geometry. They are more versatile and more expensive than network models. The following paper shows how to derive and implement adjoint thermoacoustic Helmholtz solvers in 1D and 2D:

Sensitivity Analysis of Thermoacoustic Instability with Adjoint Helmholtz Solvers

M. Juniper

*Physical Review Fluids* **3** 110509 (2018) doi:10.1103/PhysRevFluids.3.110509

M. Juniper

pdf

Open Access

doi: https://doi.org/10.1103/PhysRevFluids.3.110509

Matlab code

Erratum

Gas turbines and rocket engines sometimes suffer from violent oscillations caused by feedback between acoustic waves and flames in the combustion chamber. These are known as thermoacoustic oscillations and they often occur late in the design process. Their elimination usually requires expensive tests and re-design. Full scale tests and laboratory scale experiments show that these oscillations can usually be stabilized by making small changes to the system. The complication is that, while there is often just one unstable natural oscillation (eigenmode), there are very many possible changes to the system. The challenge is to identify the optimal change systematically, cheaply, and accurately. This paper shows how to evaluate the sensitivities of a thermoacoustic eigenmode to all possible system changes with a single calculation by applying adjoint methods to a thermoacoustic Helmholtz solver. These sensitivities are calculated here with finite difference and finite element methods, in the weak form and the strong form, with the discrete adjoint and the continuous adjoint, and with a Newton method applied to a nonlinear eigenvalue problem and an iterative method applied to a linear eigenvalue problem. This is the first detailed comparison of adjoint methods applied to thermoacoustic Helmholtz solvers. Matlab codes are provided for all methods and all figures so that the techniques can be easily applied and tested. This paper explains why the finite difference of the strong form equations with replacement boundary conditions should be avoided and why all of the other methods work well. Of the other methods, the discrete adjoint of the weak form equations is the easiest method to implement; it can use any discretization and the boundary conditions are straightforward. The continuous adjoint is relatively easy to implement but requires careful attention to boundary conditions. The Summation by Parts finite difference of the strong form equations with a Simultaneous Approximation Term for the boundary conditions (SBP--SAT) is more challenging to implement, particularly at high order or on non-uniform grids. Physical interpretation of these results shows that the well-known Rayleigh criterion should be revised for a linear analysis. This criterion states that thermoacoustic oscillations will grow if heat release rate oscillations are sufficiently in phase with pressure oscillations. In fact, the criterion should contain the adjoint pressure rather than the pressure. In self-adjoint systems the two are equivalent. In non-self-adjoint systems, such as all but a special case of thermoacoustic systems, the two are different. Finally, the sensitivities of the growth rate of oscillations to placement of a hot or cold mesh are calculated, simply by multiplying the feedback sensitivities by a number. These sensitivities are compared successfully with experimental results. With the same technique, the influence of the viscous and thermal acoustic boundary layers is found to be negligible, while the influence of a Helmholtz resonator is found, as expected, to be considerable.

Open Access

doi: https://doi.org/10.1103/PhysRevFluids.3.110509

Matlab code

Erratum

Gas turbines and rocket engines sometimes suffer from violent oscillations caused by feedback between acoustic waves and flames in the combustion chamber. These are known as thermoacoustic oscillations and they often occur late in the design process. Their elimination usually requires expensive tests and re-design. Full scale tests and laboratory scale experiments show that these oscillations can usually be stabilized by making small changes to the system. The complication is that, while there is often just one unstable natural oscillation (eigenmode), there are very many possible changes to the system. The challenge is to identify the optimal change systematically, cheaply, and accurately. This paper shows how to evaluate the sensitivities of a thermoacoustic eigenmode to all possible system changes with a single calculation by applying adjoint methods to a thermoacoustic Helmholtz solver. These sensitivities are calculated here with finite difference and finite element methods, in the weak form and the strong form, with the discrete adjoint and the continuous adjoint, and with a Newton method applied to a nonlinear eigenvalue problem and an iterative method applied to a linear eigenvalue problem. This is the first detailed comparison of adjoint methods applied to thermoacoustic Helmholtz solvers. Matlab codes are provided for all methods and all figures so that the techniques can be easily applied and tested. This paper explains why the finite difference of the strong form equations with replacement boundary conditions should be avoided and why all of the other methods work well. Of the other methods, the discrete adjoint of the weak form equations is the easiest method to implement; it can use any discretization and the boundary conditions are straightforward. The continuous adjoint is relatively easy to implement but requires careful attention to boundary conditions. The Summation by Parts finite difference of the strong form equations with a Simultaneous Approximation Term for the boundary conditions (SBP--SAT) is more challenging to implement, particularly at high order or on non-uniform grids. Physical interpretation of these results shows that the well-known Rayleigh criterion should be revised for a linear analysis. This criterion states that thermoacoustic oscillations will grow if heat release rate oscillations are sufficiently in phase with pressure oscillations. In fact, the criterion should contain the adjoint pressure rather than the pressure. In self-adjoint systems the two are equivalent. In non-self-adjoint systems, such as all but a special case of thermoacoustic systems, the two are different. Finally, the sensitivities of the growth rate of oscillations to placement of a hot or cold mesh are calculated, simply by multiplying the feedback sensitivities by a number. These sensitivities are compared successfully with experimental results. With the same technique, the influence of the viscous and thermal acoustic boundary layers is found to be negligible, while the influence of a Helmholtz resonator is found, as expected, to be considerable.

Ensembling geophysical models with Bayesian Neural Networks

U. Sengupta, M. Amos, J. Scott Hosking, C. E. Rasmussen, P. J. Young, M. P. Juniper

*34th Conference on Neural Information Processing Systems (NeurIPS 2020), Vancouver, Canada*, (2020)

U. Sengupta, M. Amos, J. Scott Hosking, C. E. Rasmussen, P. J. Young, M. P. Juniper

pdf

Open Access

doi: https://doi.org/10.17863/CAM.60032

Ensembles of geophysical models improve prediction accuracy and express uncertainties. We develop a novel data-driven ensembling strategy for combining geophysical models using Bayesian Neural Networks, which infers spatiotemporally varying model weights and bias, while accounting for heteroscedastic uncertainties in the observations. This produces more accurate and uncertainty-aware predictions without sacrificing interpretability. Applied to the prediction of total column ozone from an ensemble of 15 chemistry-climate models, we find that the Bayesian neural network ensemble (BayNNE) outperforms existing methods for ensembling physical models, achieving a 49.4% reduction in RMSE for temporal extrapolation, and a 67.4% reduction in RMSE for polar data voids, compared to a weighted mean. Uncertainty is also well-characterized, with 91.9% of the data points in our extrapolation validation dataset lying within 2 standard deviations and 98.9% within 3 standard deviations.

Open Access

doi: https://doi.org/10.17863/CAM.60032

Ensembles of geophysical models improve prediction accuracy and express uncertainties. We develop a novel data-driven ensembling strategy for combining geophysical models using Bayesian Neural Networks, which infers spatiotemporally varying model weights and bias, while accounting for heteroscedastic uncertainties in the observations. This produces more accurate and uncertainty-aware predictions without sacrificing interpretability. Applied to the prediction of total column ozone from an ensemble of 15 chemistry-climate models, we find that the Bayesian neural network ensemble (BayNNE) outperforms existing methods for ensembling physical models, achieving a 49.4% reduction in RMSE for temporal extrapolation, and a 67.4% reduction in RMSE for polar data voids, compared to a weighted mean. Uncertainty is also well-characterized, with 91.9% of the data points in our extrapolation validation dataset lying within 2 standard deviations and 98.9% within 3 standard deviations.

Adjoint methods can also be used in uncertainty quantification (UQ) and greatly reduce the cost of UQ, as long as the uncertainties do not vary very nonlinearly with the parameters.

Stability analysis of thermoacoustic nonlinear eigenproblems in annular combustors. Part 2: Uncertainty Quantification

L. Magri, M. Bauerheim, F. Nicoud, M. P. Juniper

*Journal of Computational Physics* **325** 411--421 (2016) doi:10.1016/j.jcp.2016.08.043

L. Magri, M. Bauerheim, F. Nicoud, M. P. Juniper

pdf

Open Access

doi: https://doi.org/10.1016/j.jcp.2016.08.043

Monte Carlo and Active Subspace Identification methods are combined with first- and second-order adjoint sensitivities to perform (forward) uncertainty quantification analysis of the thermo-acoustic stability of two annular combustor configurations. This method is applied to evaluate the risk factor, i.e., the probability for the system to be unstable. It is shown that the adjoint approach reduces the number of nonlinear-eigenproblem calculations by as much as the Monte Carlo samples.

Open Access

doi: https://doi.org/10.1016/j.jcp.2016.08.043

Monte Carlo and Active Subspace Identification methods are combined with first- and second-order adjoint sensitivities to perform (forward) uncertainty quantification analysis of the thermo-acoustic stability of two annular combustor configurations. This method is applied to evaluate the risk factor, i.e., the probability for the system to be unstable. It is shown that the adjoint approach reduces the number of nonlinear-eigenproblem calculations by as much as the Monte Carlo samples.

Adjoint methods are ideal for identifying the small design changes that can stabilize thermoacoustic modes. They require, however, an accurate model of the system. Devising accurate thermoacoustic models is challenging because extreme sensitivity to parameters introduces considerable systematic error if parameters cannot be estimated accurately. This is arguably the most intractible problem in thermoacoustics.

One solution is to calculate parameters from automated experimental measurements and inverse uncertainty quantification. Parameters learned on laboratory-scale rigs can then be updated as tests are performed on larger scale rigs. Our approach is to measure growth and decay rates extremely accurately using automated experiments:

Experimental sensitivity analysis and control of thermoacoustic systems

G. Rigas, N. Jamieson, L. K. B. Li, M. Juniper

*Journal of Fluid Mechanics* **787** R1 (2016) doi:10.1017/jfm.2015.715

G. Rigas, N. Jamieson, L. K. B. Li, M. Juniper

Experimental sensitivity analysis via a secondary heat source in an oscillating thermoacoustic system

N. Jamieson, G. Rigas, M. P. Juniper

*International Journal of Spray and Combustion Dynamics* **9** (4) 230--240 (2017) doi:10.1177/1756827717696325

N. Jamieson, G. Rigas, M. P. Juniper

pdf

Open Access

doi: https://doi.org/10.1177/1756827717696325

In this paper, we report the results of an experimental sensitivity analysis on a vertical electrically-heated Rijke tube. We examine the stability characteristics of the system due to the introduction of a secondary heat source. The experimental sensitivity analysis is quantified by measuring the shift in linear growth and decay rate as well as the shift in the linear frequency during periods of growth and decay of thermoacoustic oscillations. Linear growth and decay rate measurements agree qualitatively well with the theoretical predictions from adjoint-based methods of Magri & Juniper (*J. Fluid Mech.*, vol. **719**, 2013, pp. 183--202). A discrepancy in the linear frequency measurements highlight deficiencies in the model used for those predictions and shows that the experimental measurement of sensitivities is a stringent test of any thermoacoustic model. The findings suggest that adjoint-based methods are, in principle, capable of providing industry with a cheap and efficient tool for developing optimal control strategies for more complex thermoacoustic systems.

Open Access

doi: https://doi.org/10.1177/1756827717696325

In this paper, we report the results of an experimental sensitivity analysis on a vertical electrically-heated Rijke tube. We examine the stability characteristics of the system due to the introduction of a secondary heat source. The experimental sensitivity analysis is quantified by measuring the shift in linear growth and decay rate as well as the shift in the linear frequency during periods of growth and decay of thermoacoustic oscillations. Linear growth and decay rate measurements agree qualitatively well with the theoretical predictions from adjoint-based methods of Magri & Juniper (

We obtain tens of thousands of datapoints and then infer the uncertainty in parameters. This work is described further in the Physics-based Statistical Learning project.