A A A Volume : 44 Part : 2 Proceedings of the Institute of Acoustics Indirect noise is sensitive to the compact-nozzle assumptions Animesh Jain1, University of Cambridge, Cambridge CB2 1PZ, UK Luca Magri2, Imperial College London, London SW7 1AL, UK; University of Cambridge, Cambridge CB2 1PZ, UK ABSTRACT In aircraft engine combustors, incomplete mixing and air cooling give rise to flow inhomogeneities. When accelerated through the nozzle guide vane downstream of the combustor, these inhomogeneities give rise to acoustic waves. This is commonly known as indirect combustion noise. When these sound waves are reflected off the outlet and travel upstream of the combustor, they can lead to thermoacoustic oscillations. To predict indirect noise, models proposed in the literature assume the wavelength of impinging disturbances to be large with respect to the nozzle spatial extent (compact nozzle assumptions). However, in reality, the assumption might not hold; hence, a mismatch between the model and experimental results can be observed. The overarching objective of this study is to show that the acoustic transfer functions are sensitive to the compact-nozzle assumption. We show this on a realistic model of nozzle flow, which contains dissipation (non-isentropic) effects. First, we show that indirect-noise predictions are sensitive to a small change in the Helmholtz number for a nearly compact nozzle, particularly in the subsonic flow regime. Second, we show that modelling dissipation is key to mitigating the unphysical sensitivity of the compact-nozzle assumption. This study highlights the importance of the non-compact and non-isentropic assumptions for the accurate prediction of indirect noise transfer functions and thermoacoustic instability in aeronautical gas turbines. 1. INTRODUCTION Combustion noise is one of the major sources of aero-engine noise. On the one hand, the fluctuations in the unsteady heat release of the flame in the combustor give rise to the direct combustion noise [e.g., 1, 2]. On the other hand, the acceleration of flow inhomogeneities downstream of the combustor generates indirect combustion noise [e.g., 3, 4, 5, 6, 7, 8, 9]. Indirect combustion noise generated by an advected temperature inhomogeneity is referred to as entropy noise [3, 10, 11, 12], whereas, noise generated by an advected compositional inhomogeneity is referred to as compositional noise [9, 13]. These sound waves are perceived as noise when they are emitted downstream of the engine components. Importantly, when these sound waves reflect back to the combustor, they can contribute to the thermoacoustic feedback [e.g., 7, 8, 14, 15]. Hence, research efforts aim at estimating the noise transfer functions accurately with low order models. Marble & Candel (1977) [3] proposed a low-order model to predict entropy noise in a compact nozzle flow, in which the length of the nozzle is negligible as compared to the wavelength of the acoustic waves. Assuming an isentropic flow, the jump conditions were formulated to calculate the transfer functions [3, 16]. The analysis was extended to non-compact nozzles by Leyko et al.(2009) and Duran & Moreau (2013) [12, 17]. The non-compactness of a nozzle is measured using the Helmholtz number, He, which is the between the wavelength of the advected perturbations and that of the acoustic waves (He = ωL/c, where L is the characteristic nozzle length, ω is the perturbation angular frequency and c is the reference speed of sound). The compact nozzle assumption (He = 0) assumes the nozzle as a discontinuity in the flow, hence algebraic jump conditions are imposed across the nozzle to conserve mass, momentum and energy. For non-compact nozzles (He > 0), the equations become differential, which needs to be solved numerically or semi-analytically. All the above-mentioned studies assume the flow to be isentropic. However, factors such as flow recirculation and wall friction make the flow non-isentropic. The effect of non-isentropicity on indirect noise was shown experimentally by De Domenico et al. (2019) [18], who observed a mismatch between experimental observations and the isentropic model for a subsonic flow. Using conservation laws, Jain & Magri (2022) [19] developed a physical model to predict indirect noise in nozzles with dissipation, which was validated with experiments of the literature [18]. In this work, we (i) present the governing equations for entropy noise generation in non-isentropic flows; (ii) analyse the sensitivity of the entropic-acoustic transfer functions with respect to the Helmholtz number for isentropic and non-isentropic flows. The paper is structured as follows. Section 2 describes the mathematical model to predict entropy noise in nozzles with dissipation. Section 3 discusses the sensitivity of the nozzle transfer functions. Conclusions end the paper. 2. MATHEMATICAL MODEL We consider a single component, calorically perfect, ideal and adiabatic flow. The effects of dissipation (non-isentropicity) averaged over the nozzle cross-section are considered using the friction factor [19]. The cut-off frequency is considered sufficiently large to model the acoustics as quasi-one dimensional. The equations of conservation of mass, momentum and energy are [19], respectively where, u is the velocity, M is the Mach number, x is the nozzle axial coordinate (x = 0 at the inlet), p is the pressure, ρ is the density, s is the entropy, T is the temperature, R is the gas constant, γ is the heat capacity ratio, A is the cross-sectional area, D is the nozzle diameter, tan α is the spatial derivative of the nozzle profile (tan α = (1/2)(dD/dx)), and f is the Fanning friction factor. These equations are closed by Gibbs equation, where, cp is the specific heat at constant pressure. Figure 1: Nozzle schematic with nomenclature 2.1. Linearisation The acoustics are modeled as linear perturbations that develop over a steady mean flow. A generic variable is decomposed as (.) = + (.)′, where is the steady mean flow term and (.)′is the first-order perturbation. On linearising Equation 1-Equation 3, we obtaine [19], where, and the linearized Gibbs’ equation is [3] The equations are non-dimensionalized with τ = tfa , η = x / L, = D / L and = / cref , where fa is the frequency of the advected perturbations, L is the nozzle length, and cref is the reference speed of sound. The non-dimensional material derivative is / Dτ = He∂/∂t + ∂/∂η. For zero friction (f → 0) the equations tend to that of the isentropic flow [12]. The equations can be solved with a travelling wave approach, Finally, from here, the variables in the differential equations refer are Fourier transformed as (·)(x, t) → (·)(x) exp(2πiτ). The equations can be solved to calculate the indirect noise transfer functions. 3. SENSITIVITY TO THE HELMHOLTZ NUMBER We study the effect of the Helmholtz number on the entropic-acoustic transfer functions, where SR and ST are the entropic-acoustic reflection and transmission coefficients, respectively. The transfer functions are complex, thus they have both phases and magnitudes. The magnitude and the phase of the reflected wave are important parameters to analyse the effect of indirect noise on thermoacoustic stability. Whereas, the magnitude and phase of the transmitted wave contribute to the noise. We use the Taylor expansion in this section to analyse the sensitivity of the nozzle transfer function with the Helmholtz number, The transfer functions are calculated with an expansion at He = He0. The Taylor expansion of the solution corresponds to the terms of the asymptotic expansion [13]. 3.1. Subsonic Flow We analyse a convergent-divergent nozzle with the inlet area that is twice as large as the exit area. Figure 2 shows the entropic-acoustic transfer functions. The magnitude of the transfer functions in subsonic isentropic nozzle flows in nearly compact nozzles (He < 0.05, in the analysed configuration) is significantly different for a small change in the Helmholtz number. For instance, the reflection coefficient is ≈ 35 times larger at He = 0.1 as compared to that at He = 0. Similar observations can be made for the transmission coefficient. Figure 2 shows the expansion of first, second and third order in a subsonic isentropic (f = 0) nozzle flow. The predictions do not match with an expansion is at Helmholtz number, He = 0 (Figure 2 (a - d)). However, the expansion predictions tend to match the model predictions with an expansion at He ≈ 0.1, especially for higher orders (Figure 2 (e - f)). On the contrary, in a non-isentropic flow (f = 0.08), the third-order expansion matches the model predictions for both expansions at He = 0 (Figure 3 (a - d)) and He = 0.1 (Figure 3 (e - h)). This is because the acoustic waves are generated by the strong acceleration of the mean flow in the nozzle (∂/∂η in Equation 6). In a compact nozzle (He = 0), the flow goes through positive acceleration in the convergent part to a negative acceleration in the divergent part. Hence, the indirect noise generated in the convergent section is cancelled out by the noise generated in the divergent section. However, as the Helmholtz number increases, there is a phase difference associated with the waves in convergent and divergent sections. This does not allow cancellation and a drastic change in the magnitudes of the transfer functions are observed [12]. Hence the compact nozzle analysis is limited to very low frequencies in an isentropic subsonic flow. It should be treated with care. In a flow with dissipation, the indirect noise transfer functions are not as highly sensitive to the change in Helmholtz number. For instance, the reflection coefficient is ≈ 1.2 times larger at He = 0.1 as compared to the compact nozzle (Figure 3). Similar observations can be made for the transmission Figure 2: Entropic-acoustic reflection (a, c, e, g) and transmission (b,d, f, h) magnitudes and phase in an isentropic (f = 0) subsonic flow (Mt = 0.6). The black solid lines are the model predictions of Jain & Magri (2022) [19]. The solid line show the predictions using the Taylor expansion. The cross indicates the expansion point. Figure 3: Entropic-acoustic reflection (a, c, e, g) and transmission (b,d, f, h) magnitudes and phase in a non-isentropic (f = 0.08) subsonic flow (Mt = 0.6) . The black solid lines are the model predictions of Jain & Magri (2022) [19]. The solid line show the predictions using the Taylor expansion. The cross indicates the expansion point coefficient. Moreover, the predictions of the Taylor expansion match the model predictions for higher orders for expansion at He = 0 (Figure 3 (a - d)) and for all orders for an expansion at He = 0.1 (Figure 3 (e - h)). This is because, there is an additional term in the momentum equation (∂/∂η + 4 f / 2 in Equation 6). Physically, the dissipation has the effect of decelerating the mean flow and leads to a change in magnitude and phase of the indirect noise transfer functions. Hence, modelling dissipation can help mitigate of the non-physical sensitivity observed in the case with the isentropic flow assumption. In conclusion, the transfer functions are highly sensitive to the Helmholtz number for small Helmholtz numbers (He < 0.05, in the analysed case) in an isentropic subsonic flow. The sensitivity reduces in a flow with dissipation. This highlights the importance of correct estimation of the Helmholtz number and modelling dissipation, even in the nearly compact nozzle regime. 4. CONCLUSIONS Various indirect noise models in the literature consider the nozzle to be compact. In this work, we highlight the importance of non-compact assumptions. We also analyse the effect of dissipation on this assumption. To do so, first, we show the governing equations for a single component flow with dissipation derived using conservation laws. Second, we analyse the sensitivity of the entropic acoustic transfer functions to changes in the Helmholtz number. The analysis is performed for both isentropic and non-isentropic flows. We observe that the transfer functions are highly sensitive around the Helmholtz number, He = 0 in an isentropic subsonic flow. The sensitivity reduces as the Helmholtz number increases. In a subsonic flow with dissipation, the transfer functions are less sensitive to changes in the Helmholtz number. Therefore, dissipation needs to be modelled for the predictions not to be overly sensitive to the compact-nozzle assumption. In conclusion, the compact nozzle assumptions may result in inaccurate indirect noise predictions, especially in isentropic subsonic flows. This highlights the importance of modelling Helmholtz number effects with dissipation. ACKNOWLEDGEMENTS A. J. is supported by the University of Cambridge Harding Distinguished Postgraduate Scholars Programme. L.M. acknowledges the support from the ERC Starting Grant (PhyCo, n. 94938) and TUM Institute for Advanced Study (German Excellence Initiative and the EU 7th Framework Programme n. 291763). REFERENCES Yasser Mahmoudi, Andrea Giusti, Epaminondas Mastorakos, and Ann P Dowling. Low-order modeling of combustion noise in an aero-engine: the effect of entropy dispersion. Journal of Engineering for Gas Turbines and Power, 140(1), 2018. Matthias Ihme. Combustion and engine-core noise. Annual Review of Fluid Mechanics, 49:277– 310, 2017. FE Marble and SM Candel. Acoustic disturbance from gas non-uniformities convected through a nozzle. Journal of sound and vibration, 55(2):225–243, 1977. Warren C Strahle. Noise produced by fluid inhomogeneities. AIAA Journal, 14(7):985–987, 1976. NA Cumpsty. Jet engine combustion noise: pressure, entropy and vorticity perturbations produced by unsteady combustion or heat addition. Journal of Sound and Vibration, 66(4):527– 544, 1979. JE Ffowcs Williams and MS Howe. The generation of sound by density inhomogeneities in low mach number nozzle flows. Journal of Fluid Mechanics, 70(3):605–622, 1975. Wolfgang Polifke, Christian Oliver Paschereit, and Klaus Döbbeling. Constructive and destructive interference of acoustic and entropy waves in a premixed combustor with a choked exit. Int. J. Acoust. Vib, 6(3):135–146, 2001. Aimee S Morgans and Ignacio Duran. Entropy noise: A review of theory, progress and challenges. International Journal of Spray and Combustion Dynamics, 8(4):285–298, 2016. Luca Magri, Jeff O’Brien, and Matthias Ihme. Compositional inhomogeneities as a source of indirect combustion noise. Journal of Fluid Mechanics, 799, 2016. Elizabeth Cuadra. Acoustic wave generation by entropy discontinuities flowing past an area change. The Journal of the Acoustical Society of America, 42(4):725–732, 1967. Friedrich Bake, Christoph Richter, Bernd Mühlbauer, Nancy Kings, Ingo Röhle, Frank Thiele, and Berthold Noll. The entropy wave generator (EWG): a reference case on entropy noise. Journal of Sound and Vibration, 326(3-5):574–598, 2009. Ignacio Duran and Stephane Moreau. Solution of the quasi-one-dimensional linearized euler equations using flow invariants and the magnus expansion. Journal of Fluid Mechanics, 723:190–231, 2013. Luca Magri. On indirect noise in multicomponent nozzle flows. Journal of Fluid Mechanics, 828, 2017. Chee Su Goh and Aimee S Morgans. The influence of entropy waves on the thermoacoustic stability of a model combustor. Combustion Science and Technology, 185(2):249–268, 2013. Emmanuel Motheau, Franck Nicoud, and Thierry Poinsot. Mixed acoustic–entropy combustion instabilities in gas turbines. Journal of Fluid Mechanics, 749:542–576, 2014. NA Cumpsty and FE Marble. The interaction of entropy fluctuations with turbine blade rows; a mechanism of turbojet engine noise. Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences, 357(1690):323–344, 1977. Matthieu Leyko, Franck Nicoud, and Thierry Poinsot. Comparison of direct and indirect combustion noise mechanisms in a model combustor. AIAA journal, 47(11):2709–2716, 2009. Francesca De Domenico, Erwan O Rolland, and Simone Hochgreb. A generalised model for acoustic and entropic transfer function of nozzles with losses. Journal of Sound and Vibration, 440:212–230, 2019. Animesh Jain and Luca Magri. A physical model for indirect noise in non-isentropic nozzles: transfer functions and stability. Journal of Fluid Mechanics, 935:A33, 2022. 1aj584@cam.ac.uk 2l.magri@imperial.ac.uk Previous Paper 2 of 808 Next