Modelling nonlinear thermoacoustic instability in an electrically heated Rijke tube

Abstract

An analysis of thermoacoustic instability is performed for a horizontal Rijke tube with an electrical resistance heater as the heat source. The governing equations for this fluid flow become stiff and are difficult to solve by the computational fluid dynamics (CFD) technique, as the Mach number of the steady flow and the thickness of the heat source (compared to the acoustic wavelength) are small. Therefore, an asymptotic analysis is performed in the limit of small Mach number and compact heat source to eliminate the above stiffness problem. The unknown variables are expanded in powers of Mach number. Two systems of governing equations are obtained: one for the acoustic field and the other for the unsteady flow field in the hydrodynamic zone around the heater. In this analysis, the coupling between the acoustic field and the unsteady heat release rate from the heater appears from the asymptotic analysis. Furthermore, a non-trivial additional term, referred to as the global-acceleration term, appears in the momentum equation of the hydrodynamic zone, which has serious consequences for the stability of the system. This term can be interpreted as a pressure gradient applied from the acoustic onto the hydrodynamic zone. The asymptotic stability of the system with the variation of system parameters is presented using the bifurcation diagram. Numerical simulations are performed using the Galerkin technique for the acoustic zone and CFD techniques for the hydrodynamic zone. The results confirm the importance of the global-acceleration term. Bifurcation diagrams obtained from the simulations with and without the above term are different. Acoustic streaming is shown to occur during the limit cycle and its effect on the unsteady heat release rate is discussed.