An improved multimodal method for sound propagation in nonuniform lined ducts

An efficient method is proposed for modeling time harmonic acoustic propagation in a nonuniform lined duct without flow. The lining impedance is axially segmented uniform, but varies circumferentially. The sound pressure is expanded in term of rigid duct modes and an additional function that carries the information about the impedance boundary. The rigid duct modes and the additional function are known a priori so that calculations of the true liner modes, which are difficult, are avoided. By matching the pressure and axial velocity at the interface between different uniform segments, scattering matrices are obtained for each individual segment; these are then combined to construct a global scattering matrix for multiple segments. The present method is an improvement of the multimodal propagation method, developed in a previous paper [Bi et al., J. Sound Vib. 289, 1091-1111 (2006)]. The radial rate of convergence is improved from O(n(-2)), where n is the radial mode indices, to O(n(-4)). It is numerically shown that using the present method, acoustic propagation in the nonuniform lined intake of an aeroengine can be calculated by a personal computer for dimensionless frequency K up to 80, approaching the third blade passing frequency of turbofan noise.     

Comparison of measured broadband noise attenuation spectra from circular flow ducts and from lined engine intakes with predictions

Sound propagation in lined circular ducts is investigated in the presence of uniform and sheared flow. The modal solutions are obtained by solving an eigenvalue equation which, in the case of sheared flow, is derived by using finite differences and by matching the pressure and the radial component of the particle velocity at the interface of the regions of uniform and sheared flow. For the uniform flow region, standard Bessel function solutions are used. The attenuation of acoustic energy at a given frequency and for a given liner length is computed on the assumption that at the inlet to the lined duct, the acoustic energy is equally distributed among the propagating modes. The total number of propagating modes is determined from the hard wall “cut off” condition. The failure to find some of the modal solutions on the attenuation computed in this way is discussed. It is shown that the reliability of this method of computing liner attenuation depends on the ability to successfully compute most of the modal solutions over a large range of frequencies, flow conditions and duct wall impedance values. A numerical technique is developed which uses a fraction of the total number of solutions to compute the total attenuations without appreciable loss of accuracy. Measured attenuation spectra from a flow duct facility and from lined intake ducts of the RB.211 engine are compared with predictions. In general very good agreement between predictions and measurements is obtained.     

Broadband noise reduction by circular multi-cavity mufflers operating in multimodal propagation conditions

In this study, sound propagation through a circular duct with non-locally lining is investigated both numerically and experimentally. The liner concept is based on perforated screens backed by air cavities. Dimensions of the cavity are chosen to be of the order or bigger than the wavelength so acoustic waves within the liner can propagate parallel to the duct surface. This gives rise to complex scattering mechanisms among duct modes which renders the muffler more effective over a broader frequency range. This work emanates from the Cleansky European HEXENOR project which aim is to identify the best multi-cavity muffler configuration for reduction of exhaust noise from helicopter turboshaft engines. Here, design parameters are the cavity dimensions in both longitudinal and azimuthal directions. The best cavity configuration must in addition fit weight specifications which implies that the number of walls separating each cavity should be chosen as small as possible. To achieve these objectives, the scattering matrix of the lined duct section is obtained experimentally for two specific muffler configurations operating in multimodal propagation conditions. The good agreement with numerical predictions serves to validate the perforate plate impedance model used in our calculation. Finally, given an incident acoustic pressure which is representative of typical combustion noise spectrum, the best cavity configuration achieving the maximum overall acoustic Transmission Loss is selected numerically. The study also illustrates how the acoustic performances are dependent on the nature of the incident field.     

流管实验装置中声传播计算的模态方法

流管实验装置是测量有流动情况下航空发动机消声短舱内声衬声阻抗的主要装置。本文发展了一种解析的模态匹配方法进行在平均流有声衬条件下矩形流管中声传播的计算。用同伦方法求解特征值问题,并与用环绕积分求解的结果进行比较。声场通过轴向阻抗间断面的声压和声质点速度积分相等计算。第一个算例是无流动、硬壁、有限长、考虑端口反射的情况,并与北航流管实验台测量数据进行了对比;第二个算例为有流动情况下有限长声衬管道不考虑端口反射的声场计算,它与文献中NASA流管实验结果和CAA计算结果符合得很好。     

On the acoustic modes in a cylindrical duct with an arbitrary wall impedance distribution

The present paper considers the propagation of sound in a cylindrical duct, with a wall section of finite length covered by an acoustic liner whose impedance is an arbitrary function of position. The cases of (i) uniform wall impedance, and wall impedance varying along the (ii) circumference or (iii) axis of the duct, or (iv) both simultaneously, are explicitly considered. It is shown that a nonuniform wall impedance couples modes with distinct azimuthal l or axial m wave numbers, so that their radial wave numbers k can no longer be calculated separately for each pair (m,l). The radial wave numbers are the roots of an infinite determinant, in the case when the wall impedance varies either (i) circumferentially or (ii) radially. If the wall impedance varies (iv) both radially and circumferentially, then the radial wave numbers are the roots of a doubly infinite determinant, i.e., an infinite determinant in which each term is an infinite determinant. The infinite determinants specifying the radial wave numbers are written explicitly for sound in a cylindrical nozzle with a uniform axial flow, in which case the radial eigenfunctions are Bessel functions; the method of calculation of the radial wave numbers applies equally well to a cylindrical nozzle with shear flow and/or swirling flows, with the Bessel functions replaced by other eigenfunctions. The radial wave numbers are calculated by truncation of the infinite determinants, for several values of the aspect ratio, defined as the ratio of length to diameter. It is shown that a nonuniform wall impedance will give rise to additional modes compared with a uniform wall impedance. The radial wave numbers specify the eigenfrequencies for the acoustic modes in the duct; the imaginary parts of the eigenfrequencies specify the decay of the sound field with time, and thus the effectiveness of the acoustic liner.     

Application of the Wiener–Hopf method for describing the propagation of sound in cylindrical and rectangular channels with an impedance jump in the presence of a flow

Exact solutions to problems of the propagation of acoustic modes in lined channels with an impedance jump in the presence of a uniform flow are constructed. Two problems that can be solved by the Wiener- Hopf method—the propagation of acoustic modes in an infinite cylindrical channel with a transverse impedance jump and the propagation of acoustic modes in a rectangular channel with an impedance jump on one of its walls—are considered. On the channel walls, the Ingard–Myers boundary conditions are imposed and, as an additional boundary condition in the vicinity of the junction of the linings, the condition expressing the finiteness of the acoustic energy. Analytical expressions for the amplitudes of the transmitted and reflected fields are obtained.     

Experimental validation of numerical simulations for an acoustic liner in grazing flow: Self-noise and added drag

A coordinated experimental and numerical simulation effort is carried out to improve our understanding of the physics of acoustic liners in a grazing flow as well our computational aeroacoustics (CAA) method prediction capability. A numerical simulation code based on advanced CAA methods is developed. In a parallel effort, experiments are performed using the Grazing Flow Impedance Tube at the NASA Langley Research Center. In the experiment, a liner is installed in the upper wall of a rectangular flow duct with a 2 in. by 2.5 in. cross section. Spatial distribution of sound pressure levels and relative phases are measured on the wall opposite the liner in the presence of a Mach 0.3 grazing flow. The computer code is validated by comparing computed results with experimental measurements. Good agreements are found. The numerical simulation code is then used to investigate the physical properties of the acoustic liner. It is shown that an acoustic liner can produce self-noise in the presence of a grazing flow and that a feedback acoustic resonance mechanism is responsible for the generation of this liner self-noise. In addition, the same mechanism also creates additional liner drag. An estimate, based on numerical simulation data, indicates that for a resonant liner with a 10 percent open area ratio, the drag increase would be about 4 percent of the turbulent boundary layer drag over a flat wall.     

Effect of flow on quasi-one-dimensional acoustic wave propagation in a variable area duct of finite length

The general equation for the velocity potential of quasi-one-dimensional acoustic wave motion in a variable area, finite duct with one-dimensional flow is derived by using a perturbation technique. The non-linear second-order partial differential equation is linearized and then solved, by either a power series expansion method or the Runge-Kutta fourth-order method, for harmonic time dependence. The boundary condition taken at the duct mouth is that of matching the impedance of the duct sound field to that of the radiation field at the duct opening. Three axial Mach number variations along the duct axis are considered and the results obtained are compared with those for the case of constant Mach number, to determine the influence of the axial velocity gradient on sound propagation. The effect of flow on the radiation impedance is also considered.     

Acoustic propagation in partially choked converging-diverging ducts

A computer model based on the wave-envelope technique is used to study acoustic propagation in converging-diverging hard walled and lined circular ducts carrying near sonic mean flows. The influences of the liner admittance, boundary layer thickness, spinning mode number, and mean Mach number are considered. The numerical results indicate that the diverging portion of the duct can have a strong reflective effect for partially choked flows.     

The use of Herschel–Quincke tubes to improve the efficiency of lined ducts

The scattering matrix results of Herschel–Quincke (HQ) resonators installed in combination with an acoustic liner (HQ-Liners) are presented in this paper. This approach aims at controlling both tonal and broadband noise to improve the liner efficiency. It uses circumferential arrays of Herschel–Quincke tubes on a main duct in a serial association with a locally reacting liner of known impedance. Results for the scattering matrix of this system are deduced from an analytical model and compared with experimental and numerical data showing a good agreement. Analysis of the scattering matrix coefficients points out the modal conversion properties that depend on the number of HQ tubes along the circumference. Results of the transmission loss show that the choice of an optimal HQ configuration with respect to the liner properties can substantially improve the liner efficiency.     

Side-branch resonators modelling with Green׳s function methods

This paper deals with strategies for computing efficiently the propagation of sound waves in ducts containing passive components. In many cases of practical interest, these components are acoustic cavities which are connected to the duct. Though standard Finite Element software could be used for the numerical prediction of sound transmission through such a system, the method is known to be extremely demanding, both in terms of data preparation and computation, especially in the mid-frequency range. To alleviate this, a numerical technique that exploits the benefit of the FEM and the BEM approach has been devised. First, a set of eigenmodes is computed in the cavity to produce a numerical impedance matrix connecting the pressure and the acoustic velocity on the duct wall interface. Then an integral representation for the acoustic pressure in the main duct is used. By choosing an appropriate Green?s function for the duct, the integration procedure is limited to the duct–cavity interface only. This allows an accurate computation of the scattering matrix of such an acoustic system with a numerical complexity that grows very mildly with the frequency. Typical applications involving Helmholtz and Herschel–Quincke resonators are presented.     

Finite difference computation of acoustic scattering by small surface inhomogeneities and discontinuities

The use of finite difference schemes to compute the scattering of acoustic waves by surfaces made up of different materials with sharp surface discontinuities at the joints would, invariably, result in the generations of spurious reflected waves of numerical origin. Spurious scattered waves are produced even if a high-order scheme capable of resolving and supporting the propagation of the incident wave is used. This problem is of practical importance in jet engine duct acoustic computation. In this work, the basic reason for the generation of spurious numerical waves is first examined. It is known that when the governing partial differential equations of acoustics are discretized, one should only use the long waves of the computational scheme to represent or simulate the physical waves. The short waves of the computational scheme have entirely different propagation characteristics. They are the spurious numerical waves. A method by which high wave number components (short waves) in the wave scattering process is intentionally removed so as to minimize the scattering of spurious numerical waves is proposed. This method is implemented in several examples from computational aeroacoustics to illustrate its effectiveness, accuracy and efficiency. This method is also employed to compute the scattering of acoustic waves by scatterers, such as rigid wall acoustic liner splices, with width smaller than the computational mesh size. Good results are obtained when comparing with computed results using much smaller mesh size. The method is further extended for applications to computations of acoustic wave reflection and scattering by very small surface inhomogeneities with simple geometries.     

Application of the equivalent surface source method to the acoustics of duct systems with non-uniform wall impedance

The paper outlines the application of the equivalent surface source method to the analysis of the acoustic field in a partially lined duct with arbitrarily non-uniform wall impedance. Lined sections of the duct wall are represented by unsteady mass source singularities, the strengths of which are determined by solving integral equations. The method is applicable to lined walls of impedance which is non-uniform in the streamwise and/or circumferential direction. Numerical examples are given to show the effects of various design parameters on sound attenuation. Some advantageous features of circumferentially non-uniform wall impedance are demonstrated.     

A straightforward method for wall impedance eduction in a flow duct

The development of the advanced liner technology for aeroengine noise control necessitates the impedance measurement method under realistic flow conditions. Currently, the methods for this need are mainly based on the inverse impedance eduction principle, confronting with the problems of initial guess, high computation cost, and low convergence. In view of this, a new strategy is developed that straightforwardly educes the impedance from the sound pressure information measured on the duct wall opposing to the test acoustic liner embedded in a flow duct. Here, the key insight is that the modal nature of the duct acoustic field renders a summed-exponential representation of the measured sound pressure; thus, the characterizing axial wave number can be readily extracted by means of Prony's method, and further the unknown impedance is calculated from the eigenvalue and dispersion relations based on the classical mode-decomposition analysis. This straightforward method is simple in its basic principle but remarkably has the advantages of ultimately overcoming the drawbacks inherent to the inverse methods, incorporating the realistic multimode nonprogressive wave effects, high computational efficiency, possibly reducing the measurement points, and even avoiding the necessity of the duct exit impedance that bothers perhaps all the existing waveguide methods.     

Time-domain characterization of the acoustic damping of a perforated liner with bias flow

Combustion instabilities are caused by the interaction of unsteady heat releases and acoustic waves. To mitigate combustion instabilities, perforated liners, typically subjected to a low Mach number bias flow (a cooling flow through perforated holes), are fitted along the bounding walls of a combustor. They dissipate the acoustic waves by generating vorticity at the rims of perforated apertures. To investigate the absorption of plane waves by a perforated liner with bias flow, a time-domain numerical model of a cylindrical lined duct is developed. The liners' damping mechanism is characterized by using a time-domain "compliance." The development of such time-domain compliance is based on simplified or unsimplified Rayleigh conductivity. Numerical simulations of two different configurations of lined duct systems are performed by combining a 1D acoustic wave model with the compliance model. Comparison is then made between the results from the present models, and those from the experiment and the frequency-domain model of previous investigation [Eldredge and Dowling, J. Fluid Mech. 485, 307-335(2003)]. Good agreement is observed. This confirms that the present model can be used to simulate the propagation and dissipation of acoustic plane waves in a lined duct in real-time.     

Bulk reaction modeling of ducts with and without mean flow

A general formulation for analysis of sound field in a uniform flow duct lined with bulk-reacting sound-absorbing material is presented here. Presented theoretical model predicts the rate of attenuation for symmetric as well as asymmetric modes in rectangular duct lined with loosely bound (bulk-reacting) sound-absorbing material, which allows acoustic propagation through the lining. The nature of attenuation in rectangular ducts lined on two and four sides with and without mean flow is discussed. Computed results are compared with published theoretical and experimental results. The presented model can be used as guidelines for the acoustic design of silencers, air-conditioning ducts, industrial fans, and other similar applications.     

On the prediction of “buzz-saw” noise in acoustically lined aero-engine inlet ducts

Aero-engines operating with supersonic fan tip speeds generate an acoustic signature containing energy spread over a range of harmonics of the engine shaft rotation frequency. These harmonics are commonly known as the “buzz-saw” tones. The pressure signature attached to a supersonic ducted fan will be a sawtooth waveform. The non-linear propagation of a high-amplitude irregular sawtooth upstream inside the inlet duct redistributes the energy amongst the buzz-saw tones. In most modern aero-engines the inlet duct contains an acoustic lining, whose properties will be dependent on the mode number and frequency of the sound, and the speed of the oncoming flow. Such effects may not easily be incorporated into a time-domain approach; hence the non-linear propagation of an irregular sawtooth is calculated in the frequency domain, which enables liner damping to be included in the numerical model. Results are presented comparing noise predictions in hard-walled and acoustically lined inlet ducts. These show the effect of an acoustic liner on the buzz-saw tones. These predictions compare favourably with previous experimental measurements of liner insertion loss (at blade passing frequency), and provide a plausible explanation for the observed reduction in this insertion loss at high fan operating speeds.     

Modal decomposition method for acoustic impedance testing in square ducts

Accurate duct acoustic propagation models are required to predict and reduce aircraft engine noise. These models ultimately rely on measurements of the acoustic impedance to characterize candidate engine nacelle liners. This research effort increases the frequency range of normal-incidence acoustic impedance testing in square ducts by extending the standard two-microphone method (TMM), which is limited to plane wave propagation, to include higher-order modes. The modal decomposition method (MDM) presented includes four normal modes in the model of the sound field, thus increasing the bandwidth from 6.7 to 13.5 kHz for a 25.4 mm square waveguide. The MDM characterizes the test specimen for normal- and oblique-incident acoustic impedance and mode scattering coefficients. The MDM is first formulated and then applied to the measurement of the reflection coefficient matrix for a ceramic tubular specimen. The experimental results are consistent with results from the TMM for the same specimen to within the 95% confidence intervals for the TMM. The MDM results show a series of resonances for the ceramic tubular material exhibiting a monotonic decrease in the resonant peaks of the acoustic resistance with increasing frequency, resembling a rigidly-terminated viscous tube, and also evidence of mode scattering is visible at the higher frequencies.     

Predicting insertion loss of large duct systems above the plane wave cutoff frequency

If the dimensions of a silencer or muffler component are small compared to an acoustic wavelength, plane wave propagation can be assumed. This is not the case for HVAC (heating, ventilation, and air conditioning) duct systems, and large diesel engine mufflers commonly used in ship and generator sets. For such applications, the wave behavior in the inlet and outlet ducts is three-dimensional. In this paper, the finite element method is utilized to simulate large duct systems with an aim to predict the insertion loss. The boundary condition on the source side is a diffuse field applied by determining a suitable cross-spectral force matrix of the excitation. At the termination, the radiation impedance is calculated utilizing a wavelet algorithm. Simulation results are compared to published measurement results for HVAC plenums and demonstrate good agreement.