| description abstract | The aim of this study is to efficiently calculate parametric thermoacoustic maps of typical combustion chambers. Two configurations are considered: an academic configuration based on a Rijke tube, and an industrial combustion chamber, which is the core of a recently developed microturbine for power generation. Such maps can be understood as the collection of loci of thermoacoustic eigenfrequencies obtained under systematic variations of some defined parameters, while considering the Helmholtz equation as the thermoacoustic model of interest. In this study we consider variations on two parameters: the gain n and time-delay τ associated with a generic flame response model. We also show the feasibility of the proposed approach when considering more realistic flame responses. A straight-forward way to calculate such a thermoacoustic map is by solving the Helmholtz equation, and, thus, the corresponding nonlinear eigenvalue problem (NLEVP), one time per parameter combination. With that approach, the nonlinear eigenvalue problem needs to be solved hundreds or thousands of times if an adequate resolution of the thermoacoustic map is sought. Such a strategy may be computationally unaffordable. In order to overcome this difficulty, this work utilizes an adjoint-based, high-order perturbation method. The actual eigenvalue problem is only solved once at a baseline point. After applying the perturbation equations at that point, a polynomial rational function—the Padé approximant—is obtained to estimate the eigenfrequency drift that results for a small or large perturbation in the flame response. It is demonstrated, for both academic and industrial test cases, that the obtained maps are accurate. Additionally, it is shown that these maps reveal a large variety of thermoacoustic features, such as stability boundaries, intrinsic thermoacoustic modes, and exceptional points. The numerical costs for such calculations are negligible even for the industrial combustion chamber investigated. | |
