WKB solutions to the quasi 1-D acoustic wave equation in ducts with non-uniform cross-section and inhomogeneous mean flow properties – Acoustic field and combustion instability

Abstract

Linear modal analysis based on the spatio-temporal harmonic representation of acoustic fluctuations is a well-known reduced order method to predict combustion instability frequencies. However, a limitation of this method is that it is only applicable to domains with homogeneous mean properties. The Wentzel-Kramers-Brillouin (WKB) method facilitates the development of approximate analytical solutions to the one-dimensional (1-D) acoustic wave equation in domains with inhomogeneous mean properties. The “classical” form of the WKB solution to the wave equation is based on the assumption of high frequency (or, short wavelength), and is applicable to uniform cross-section ducts with a mean temperature gradient, but no mean flow. Subsequently, the WKB method was extended to include uniform mean velocity as well. In the current study, we derive a WKB-type solution to the generalized Helmholtz equation in quasi 1-D ducts with non-uniform cross-sectional areas and inhomogeneities in mean flow properties such as the velocity, temperature, density, and pressure. The WKB solution is applicable under the assumptions of high frequency and slowly varying mean properties. In deriving this solution, a novel approach has been developed to relate the first and second spatial derivatives of the axial velocity fluctuations to the pressure fluctuations (and its derivatives).A second WKB-type solution (referred to as the WKB2 solution) is also obtained by relaxing the slowly-varying mean property assumption, but retaining the high frequency approximation. Thirdly, the Helmholtz equation for pressure is solved numerically without invoking the two WKB assumptions. The WKB, WKB2 and numerical solutions are compared for a number of cases involving combinations of non-uniform duct areas and axially varying mean properties. Finally, we apply the WKB solution to predict the longitudinal combustion instabilities in a 1-D combustor characterized by a sudden expansion, with uniform mean flow in each of the two ducts comprising the combustor, and a linear temperature profile downstream of a compact, planar flame.