Case studies on aeroacoustics damping performances of Coupled Helmholtz Resonators over low frequency ranges

Abstract

Helmholtz resonators are widely applied in industries to reduce the transmission of unwanted noise and passively attenuate tonal noise in pipes/ducts. For example, internal combustion engines take advantage of Helmholtz resonators by combining air box, subwoofers and perforated pipes/plates to reduce engine noise emission nowadays. However, the main drawback of the conventional Helmholtz resonators includes a very narrow bandwidth. To broaden the effective noise damping frequency range of Helmholtz resonators, a number of different improved designs are reported and tested recently in the literature such as applying an array of separated resonators or implementing mechanical-coupled Helmholtz resonators. In this work, we conduct a brief review on the aeroacoustics damping performance of coupled Helmholtz resonators over low frequency range. The noise damping performances could be characterized and quantified in terms of transmission loss and sound absorption coefficient. Both definitions are discussed. Comparison is then made between the separately working resonators and the mechanical-coupled ones. This is achieved by replacing the sharable flexible sidewall with a rigid sidewall. The effect of the mean grazing flow is also discussed and examined on the noise damping performances. An experimental case study of the present authors is included. Additionally, to improve the noise damping performance further, an overview of the resonators neck with different noise damping materials inserted is conducted. Finally, 3D numerical approaches in simulating the aeroacoustics damping performances of coupled Helmholtz resonators are reviewed. This case study and brief overview provides a guidance on improving the design of coupled Helmholtz resonators and an array of Helmholtz resonators.

Keywords

Helmholtz resonator aeroacoustics transmission loss acoustic damper sound absorption coefficient

Introduction

Undesirable engine noises^1–3 could be attenuated by applying Helmholtz resonators.^4–6 A Helmholtz resonator^7–10 has a shape resembling a beer bottle, containing an enclosed volume of air,^11,12 which communicates with the engine pipe via a small opening/neck.^13–15 The cavity volume air may resonate at its natural resonant frequency.^16–18 It depends on the cavity volume and the geometry area of the neck/opening. At resonance, the combustion noises^19–23 could be significantly attenuated by the thermos-viscous effect²⁴ and vortex shedding²⁵ across the neck of the resonator. The resonators’ damping performance is dramatically deteriorated at a far off-resonant frequency.^26–28 To improve the resonators’ noise damping performances over a broad frequency range, actively passive control of acoustic damper such as Helmholtz resonators or perforated liners^29,30 could be applied to tune the resonant frequency of the acoustic dampers in real-time. In this way, it enables acoustic damper corresponding to the operating condition changes of the engines.^31–34

Helmholtz resonators^4,5,8 have a wide application in air-breathing engine industries, such as gas turbine engines¹¹ and internal combustion engines.¹² Extensive research^16–18 has been conducted to optimize the design and the acoustic noise damping performance of such resonators at a low frequency range, typically below 800 Hz.^3,12,18 It is possible that a Helmholtz resonator^35–37 is mechanically coupled with a flexible membrane additionally, which can be used as a coupling sidewall/element to improve the transmission loss^38,39 and stabilize combustion systems.^40,41 Previous experimental and numerical studies^42–45 have shown that Helmholtz resonators could be applicable for sound absorption^46–49 over large bandwidths provided multiple resonators that tuned in different frequencies;^50–55 Helmholtz resonator allows narrow frequency bandwidth with high transmission loss.^56–62 Besides, estimated transmission loss^63–66 can be improved by inserting acoustic porous materials.⁴ Fibreglass is a porous material with good sound absorption and low transmission losses ⁵. Thus, it is imperative to analyse the fundamental approaches in three parameters: (1) acoustic materials properties, (2) coupling mechanism in experimental set-up and (3) frequency variation applied on Helmholtz resonators.^67–73 In this work, we conduct experimental case studies on investigating acoustic damping performances of coupled Helmholtz resonators in different shapes over the low frequency ranges. Here we also consider the effects (1) acoustic porous materials, (2) the grazing flow Mach number. This is aimed to identify the optimum design for combustion engine applications ^74–85. Before the case study is presented, it would be informative to overview of theoretically predicting the resonant frequency of Helmholtz resonators and characterizing the resonator’s acoustic damping performances in various ways.

Overview of Helmholtz resonators and characterizing resonator’s aeroacoustics damping

Overview of the noise damping mechanism and resonant frequency of Helmholtz resonators

Review of noise damping mechanism of conventional Helmholtz resonators

The acoustic damping mechanism^86–88 of Helmholtz resonators consists of the conversion of acoustical energy into vortical energy. They are also known as linear and nonlinear damping mechanism. Incident large-amplitude sound waves interact with the shear layers formed at the resonator neck rims. It is subjected to acoustic pressure fluctuations. The shear layers destabilized to form vortex rings. The vortices are then dissipated with insignificant broadband sound production. As long as certain frequencies, depending on the operating and geometrical characteristics, the incident acoustic waves could be almost completely attenuated by the resonator compared by a resonant cavity and a neck.

Helmholtz resonator^89–91 is a specially built damper, which is enabled to observe resonant response at specific frequencies when subject to external vibration. The volume of the resonator V, neck area S and length L determine the frequency subject to resonation in constructing the Helmholtz resonator

λ ≫ L, λ ≫ V^{\frac{1}{3}}, λ ≫ S^{\frac{1}{2}}

(1)

The necks of Helmholtz resonators allow air to radiate, which contributes to resistance flow against the downstream flow, when sound is propagated from the source. This is defined as the radiation resistance, $R_{r}$ , and the radiation of mass, $R_{m}$ will result in the air pushing out from the neck. The corrected effective length of the neck and the mass of the air contained in the neck, are corresponded at lower frequencies. Assuming that when the air contained in the resonator’s neck moves a unit, this movement provides an additional unit of mass element. Besides, thermos-viscous losses,^92–96 $R_{w}$ , at the side walls and the stiffness result in air mass moving against the walls and due to its compression of air. The total effective length of the neck and the mass are written as $L^{'} = L + 1.2 R .$ Here $L^{'}$ = the effective length of the neck, L= the length of the neck. The mass is given as $m = ρ_{₀} S L^{'}$ .

The stiffness of the system is obtained by modelling the neck to be a vibrating piston. When the piston is compressed inward by x-distance, the volume of the resonator is changed, and the pressure is varied as

Δ V = - S x a n d \frac{∆ ρ}{ρ} = - \frac{∆ V}{V} = \frac{S x}{V}

(2)

P = \frac{ρ_{0} c^{2} ∆ ρ}{ρ} = \frac{ρ_{0} c^{2} S x}{V}

(3)

To maintain the displacement, x, the force, f is needed, where f =SP=sx; s is the effective stiffness, that is, given as

s = \frac{ρ_{0} c^{2} S^{2}}{V}

(4)

The wall resistance and the radiation resistance provide the total resistance against downstream propagation. These downstream waves have an instantaneous complex component, which acts on the resonator opening. This is modelled by

F = S P e^{i ω t}

(5)

The differential equation resulted from manipulating equations (11)–(13) is given as

m \frac{d^{2} x}{d t^{2}} + (R_{r} + R_{w}) \frac{d x}{d t} + s x = S P e^{i ω t}

(6)

Since the resonator^97–101 input mechanical impedance is

Z_{m} = (R_{r} + R_{w}) + j (ω m - \frac{s}{ω})

(7)

where resonance occurs, when reactance diminishes to zero. Then, we derive an equation where the resonant frequency is independent of the cavity shape

ω_{0} = {c (\frac{S}{L^{'} V})}^{\frac{1}{2}}

(8)

Equation (8) is valid for any shape, as long as the length constraints are much smaller than a single wavelength, and the neck is not too large.

Review of coupled Helmholtz resonators

Conventional Helmholtz resonators are highly effective over a narrow frequency band only to reduce the transmission of wanted noise.⁴ This limits its application in practice due to the narrow bandwidth and the ineffectiveness on attenuating acoustic disturbances of varying frequencies. In order to broaden the frequency bandwidth of the noise attenuation, and to be adaptive to dissipate acoustic perturbations of varying frequency, coupled Helmholtz resonators were designed and proposed. The first mechanical-coupled Helmholtz resonators were proposed by Griffin and his co-workers⁴ about 2 decades ago. Other designs of coupled Helmholtz resonators and the tests^63–69 have shown that flexible membrane-coupled Helmholtz resonators could be effective over a wider frequency bandwidth. And the resonant frequencies are altered⁷³ due to the dynamic variation of the cavity volume resulting from the incident large-amplitude sound waves. As revealed in equation (8), the conventional resonators cavity volume V=constant. However, the flexible membrane will lead to the volume fluctuating V≠const. This is periodic pressure fluctuations in the cavities of the coupled Helmholtz resonators. This is confirmed in the experimental test.⁷³ When a thin elastic membrane is stacked between two resonators connected in parallel, it leads to more attenuation of the unwanted sound over a wider frequency bandwidth. Some noise damping materials are filled in the necks of the resonators, enabling to evaluate the transmission loss in the coupled Helmholtz resonators. The model illustrated in Figure 1 shows how the coupled resonators are connected with a closed-end and open-end air space along the air duct downstream.

Figure 1.

Mechanically modelled coupled resonator structure.

In this duct, the problem is modelled in a classical mechanical style, dealing with the actual physical displacement of air within the duct and resonators as well as the membrane. Since the complex pressures of the propagating and reflected waves have already been modelled in Appendix equation (49) and the respective quantities of the resonators have been indicated, we shall use the same terminology in this description as well. However, since the experimental tests have been conducted in an anechoic chamber, we shall assume that the reflected wave coefficient, $C_{-}$ = 0. At the lower frequencies, we can assume that there is continuity in acoustic pressure along the duct.

The system deals with the actual movement of air in the resonator, the rectangular-shaped duct and the membrane. Equation (47) demonstrates the pressure acted by the reflected waves, propagating waves as well as the resonators; thus, the same terminology has been applied. Nevertheless, the reflected wave coefficient, X₋ is assumed to be zero due to the echo-free chamber, in which either sound or electromagnetic waves are absorbed to insulate from exterior noise sources. Another assumption has been made at low frequencies; the pressure in the duct is expressed as such, where the quantities are defined in Appendix 1; L₂ is the distance between two resonators

T ₊ + T ₋ = W ₊ + W ₋ = P_{1}

(9a)

W_{+} e^{- i K L_{2}} + W_{-} e^{i K L_{2}} = X_{+} e^{- j K L_{2}} + X_{-} e^{j K L_{2}} = P_{2}

(9b)

The conservation of the volume of velocity is applied below at any values of the resonant frequencies at Helmholtz resonators as

{V = - S x a n d U}_{T} = U_{W} - S ẋ ₁

(10a)

where

U_{T}

is the volume velocity before the neck of the resonator 1,

U_{W}

is the volume velocity after the resonator’s neck and Sẋ₁ is the volume velocity into the neck of the HR1. Therefore, the surface area and the pressure amplitude in each volume velocity component are written as follows

S ν [T^{+} - T^{-}] = S ν [W^{+} - W^{-}] - [ρ c S^{1} ẋ^{1}] f o r H R 1

(10b)

U_{W} = U_{X} + S ẋ ₂

(11a)

S γ [X ₊ - X ₋] = S γ [W ₊ e ˉ i k l ₂ - W ₋ e i k l ₂] - [ρ c S ₂ ẋ ₂] f o r H R 2

(11b)

where Sγ is the area of the air duct in the neck and

β = \frac{ρ c ₂}{A}

is for each resonator.

The motions of each individual resonator can be mathematically described as follows

m ₁ ẍ ₁ + β ₁ S ₁ 2 x ₁ + β ₁ S ₁ S ₃ x ₃ = - P ₁ S ₁

(12a)

m_{2} ẍ_{2} - β_{2} S_{2} S_{3} x_{3} + β_{2} S_{2}^{2} x_{2} = - P_{2} S_{2}

(12b)

m ₃ ẍ ₃ + β ₁ S ₃ 2 x ₃ + β ₁ S ₁ S ₃ x ₁ + β ₂ S ₃ 2 x ₃ - β ₂ S ₂ S ₃ x ₂ + k ₃ x ₃ + ς_{3} ẋ ₃ = 0

(12c)

where k and

ς

are the spring constant of the membrane and its damping properties, respectively. The Helmholtz resonator’s system can be solved by using Equations (9)–(12).

Characterizing Aeroacoustic Damping Performances of Helmholtz resonators

Definition of transmission loss

Sound waves travel spherically from a point source of the sound, and when it does, it becomes less intense. This is because with a source emitting constant power, the area of the surface of the sphere increases; thus, the sound intensity at any point on the sphere must decrease. Transmission loss occurred, when there is a loss in sound intensity due to the opposing wave propagation from the source to the downstream. Sound pressure is inversely proportional to the distance from the source, as long as the sound does not encounter any obstacles.

The experimental tests mainly focus on the sound transmission loss, geometry differences and the damping materials’ behaviour affecting the Helmholtz resonators. Different methods are introduced to find out the acoustic properties and their behaviours in the system. In the previous researches,^54,76 one used the standing-wave-ratio (SWR) method to measure the maximum and the minimum values in the standing wave tube by mobilizing the microphone at different positions. The applied method requires a minimum distance of one wavelength design test tube, whereas the lower frequencies are not suitable for this application according to λ= v/f ⁵³. The accuracy is unlikely to be improved if the uncertainty of the first minimum position is unable to determine, and errors would be introduced. Additionally, the aspect of the diameter and displacement of the microphone is limited in the SWR method, where a large diameter is advisable.

Incident and reflected acoustic impulses can be separated to determine the acoustic properties of the system in a long tube. A microphone is assembled at the origin, where the tube has to be significantly long enough so that the total duration of the time is sufficient for the sound waves to propagate from upstream to downstream and to reflect back to the origin. The properties of the system would be able to identify through the values obtained from the sound magnitude and phase. In order to have the whole system work correctly, the tube can’t be excessively long because the dissipation will contribute to the major issue during the experimental tests.⁵³ However, the experimental tests were conducted using six microphones. Three of them are located upstream, and the remaining 3 microphones are placed downstream to determine the impact of the acoustic properties filled with noise damping materials in Helmholtz resonators. Upstream microphones are to collect data from the incident and reflected waves to determine the reflection coefficient, the sound absorption coefficient and transmission loss.

Theory of transmission loss

Following equation (47), the pressure of a harmonic wave propagating in an air duct is expressed as a function of distance, x. The equation can be described further with cross- and auto-spectral densities as

p (x_{1}, t) = p_{i} (x_{1}, t) + p_{r} (x_{1}, t) = e^{i ω t} [a (t) e^{- i K x_{1}} + b (t) e^{i K x_{1}}]

(13)

where

x_{1}

represents the position of the first microphone counted from the loudspeaker to the downstream. Similarly, the remaining microphones are numbered accordingly, starting from the loudspeaker. When the sound wave propagates, the cross- and auto-spectral densities of the incident and reflected waves shall be known, before mathematical manipulation is further applied. Three microphones are set up in the upstream duct, where the cross-spectrums^38,39 of the microphones are listed by taking microphone 1 as a reference. The corresponding pressure spectrums measured from the microphones can be expressed as

{\hat{S}}_{11} (f) = [\frac{1}{T}] {{\hat{P}}_{1} (f, T) {\hat{P}}_{1}^{*} (f, T)}

(14)

{\hat{S}}_{22} (f) = [\frac{1}{T}] {{\hat{P}}_{2} (f, T) {\hat{P}}_{2}^{*} (f, T)}

(15)

{\hat{S}}_{33} (f) = [\frac{1}{T}] {{\hat{P}}_{3} (f, T) {\hat{P}}_{3}^{*} (f, T)}

(16)

{\hat{S}}_{12} (f) = [\frac{1}{T}] {{\hat{P}}_{1} (f, T) {\hat{P}}_{2}^{*} (f, T)}

(17)

{\hat{S}}_{13} (f) = [\frac{1}{T}] {{\hat{P}}_{1} (f, T) {\hat{P}}_{3}^{*} (f, T)} = {{\hat{C}}_{13} (f) + i {\hat{Q}}_{13} (f)}

(18)

where,

{\hat{S}}_{12} (f)

and

{\hat{S}}_{13} (f)

are the cross-spectrums of the microphone 2 and microphone 3 with respect to microphone 1.

{\hat{C}}_{1, n} (f)

and

{\hat{Q}}_{1, n} (f)

are the real and imaginary parts of the cross-spectrum of any microphones with respect to microphone 1.

{\hat{P}}_{n} (f)

is the Fourier Fast Transform (FFT)³⁸ of the pressure amplitudes corresponding to that microphone.

By merging equations (9a), (9b) and (10), we can get equations that can be solved analytically via matrix manipulation for any two microphones. In our experimental test, the third microphone adds an additional data point that can be further used for averaging purposes for more reliable and repeatable results

{\hat{S}}_{11} (f) = {\hat{S}}_{A A} (f) + {\hat{S}}_{B B} (f) + 2 [{\hat{C}}_{A B} (f) \cos (k_{i} + k_{r})] x_{1} + {\hat{Q}}_{A B} (f) \sin (k_{i} + k_{r}) x_{1}

(19)

{\hat{S}}_{22} (f) = {\hat{S}}_{A A} (f) + {\hat{S}}_{B B} (f) + 2 [{\hat{C}}_{A B} (f) \cos (k_{i} + k_{r})] x_{2} + {\hat{Q}}_{A B} (f) \sin (k_{i} + k_{r}) x_{2}

(20)

{\hat{C}}_{12} (f) = {\hat{S}}_{A A} (f) \cos k_{i} (x_{1} - x_{2}) + {\hat{S}}_{B B} (f) \cos k_{r} (x_{1} - x_{2}) + {\hat{C}}_{A B} (f) [\cos (k_{r} x_{1} + k_{i} x_{2}) + \cos (k_{i} x_{1} + k_{r} x_{2})] + {\hat{Q}}_{A B} (f) [\sin (k_{r} x_{1} + k_{i} x_{2}) + \sin (k_{i} x_{1} + k_{r} x_{2})]

(21)

{\hat{Q}}_{12} (f) = - {\hat{S}}_{A A} (f) \sin k_{i} (x_{1} - x_{2}) + {\hat{S}}_{B B} (f) \sin k_{r} (x_{1} - x_{2}) + {\hat{C}}_{A B} (f) [- \sin (k_{i} x_{1} + k_{r} x_{2}) + \sin (k_{r} x_{1} + k_{i} x_{2})] + {\hat{Q}}_{A B} (f) [\cos (k_{i} x_{1} + k_{r} x_{2}) - \cos (k_{r} x_{1} + k_{i} x_{2})]

(22)

where

{\hat{S}}_{A A}

{\hat{S}}_{B B}

are incident and reflected plane waves measured upstream;

{\hat{C}}_{A B}

and

{\hat{Q}}_{A B}

are the real and imaginary parts of the cross-spectrum of the incident and reflected plane waves. Theoretically, any pressure differential within the duct would give rise to the presence of a mean flow speed of air. However, since the plane sound waves create differential pressure induced by air flow is negligible compared to the speed of sound. Therefore, we can assume the mean flow to be zero. Thus, equations (19)–(22) can be reduced to

{\hat{S}}_{11} (f) = {\hat{S}}_{A A} (f) + {\hat{S}}_{B B} (f) + 2 [{\hat{C}}_{A B} (f) \cos 2 k x_{1} + {\hat{Q}}_{A B} (f) \sin 2 k x_{1}]

(23)

{\hat{S}}_{22} (f) = {\hat{S}}_{A A} (f) + {\hat{S}}_{B B} (f) + 2 [{\hat{C}}_{A B} (f) \cos 2 k x_{2} + {\hat{Q}}_{A B} (f) \sin 2 k x_{2}]

(24)

{\hat{C}}_{12} (f) = [{\hat{S}}_{A A} (f) + {\hat{S}}_{B B} (f)] c o s k (x_{1} - x_{2}) + 2 [{\hat{C}}_{A B} (f) c o s k (x_{1} + x_{2}) + {\hat{Q}}_{A B} (f) s i n k (x_{1} + x_{2})]

(25)

{\hat{Q}}_{12} (f) = [- {\hat{S}}_{A A} (f) + {\hat{S}}_{B B} (f)] \sin k (x_{1} - x_{2})

(26)

The above equations can also be applied to microphones 1 and 3 by switching out microphone 2 with 3. All the quantities above are obtained from the recorded pressure data, and can be solved for ${\hat{S}}_{A A}$ , ${\hat{S}}_{B B}$ , ${\hat{C}}_{A B}$ and ${\hat{Q}}_{A B}$ via matrix multiplication, as given below

[A] X = b^{T}

(27)

{[A]}^{- 1} A * X = {[A]}^{- 1} * b^{T}

(28)

X = {[A]}^{- 1} * b^{T}

(29)

[A] = [\begin{array}{c} 1 & 1 & \begin{array}{c} 2 * [\cos (k_{i} + k_{r}) x_{1} & 2 * [\sin (k_{i} + k_{r}) x_{1}] \end{array} \\ 1 & 1 & \begin{array}{c} 2 * [\cos (k_{i} + k_{r}) x_{x}] & 2 * [\sin (k_{i} + k_{r}) x_{2}] \end{array} \\ \begin{array}{c} \cos [k_{i} * (x_{1} - x_{2})] \\ - \sin [k_{i} * (x_{1} - x_{2})] \end{array} & \begin{array}{c} \cos [k_{r} * (x_{1} - x_{2})] \\ \sin [k_{r} * (x_{1} - x_{2})] \end{array} & \begin{array}{c} \begin{array}{c} \cos (k_{r} x_{1} + k_{i} x_{2}) + \cos (k_{i} x_{1} + k_{r} x_{2}) \\ \sin (k_{r} x_{1} + k_{i} x_{2}) - \sin (k_{i} x_{1} + k_{r} x_{2}) \end{array} & \begin{array}{c} \sin (k_{r} x_{1} + k_{i} x_{2}) + \sin (k_{i} x_{1} + k_{r} x_{2}) \\ \cos (k_{i} x_{1} + k_{r} x_{2}) - \cos (k_{r} x_{1} + k_{i} x_{2}) \end{array} \end{array} \end{array}]

b^{T} = [{\hat{S}}_{11}, {\hat{S}}_{22}, {\hat{C}}_{12}, {\hat{Q}}_{12}] a n d X^{T} = [{\hat{S}}_{A A}, {\hat{S}}_{B B}, {\hat{C}}_{A B,} {\hat{Q}}_{A B}]

(30)

By combining equations (14)–(30), the pressure spectra of the upstream and downstream waves can be measured to determine the transmission loss. The experimental test was carried out in an anechoic chamber, assuming that there are no downstream reflected waves. Nevertheless, the transmission loss values of the sound wave propagating in the air duct in any geometrical form can be calculated using the following equations

T L = 10 \log_{10} [\frac{{\hat{S}}_{A A} (f)}{{\hat{S}}_{C C} (f)}]

(31)

Δ = 1 - [\frac{{\hat{S}}_{B B} (f) + {\hat{S}}_{C C} (f)}{{\hat{S}}_{A A} (f) + {\hat{S}}_{D D} (f)}]

(32)

where

{\hat{S}}_{C C}

{\hat{S}}_{D D}

represents auto-spectrums of the incident waves and reflected waves; whilst,

{\hat{S}}_{A A}

is the incident and reflected waves measured upstream;

{\hat{C}}_{C D}

and

{\hat{Q}}_{C D}

are defined as real and imaginary parts of the cross-spectrum of the incident and reflected waves. Sound absorption coefficient Δ is also typically used to characterize the acoustical energy being absorbed by the acoustic damper. It is defined in equation (32).

Experimental case studies

Pressure sensors calibration

To demonstrate the procedures and quantify the aeroacoustics damping performances of coupled Helmholtz resonators, an experimental study⁶¹ is conducted as a case study. Here, six microphones were distributed, three upstream and three downstream, and utilized in Helmholtz resonators attached to a pipe with a grazing flow present. The frequency response of each microphone (see Figure 3) would be slightly different depending upon the tolerances and manufacturability. Hence, the microphones are required to calibrate as an additional process prior to the experiments for reliability and repeatability. However, the difference of 4 dB is acceptable, as observed in the previous experimental studies.^55,56,63

A calibration chamber (see Figure 2) was constructed with six holes circularly drilled inside where the microphones are flush-mounted for calibration. The box was designed and built in such a way that it is rigid and heavy enough to prevent the noise generated from structure vibration. Furthermore, the mounting position of the microphones and the loudspeaker are maintained at a constant distance for repeatability. Sound level metres were placed in the calibration chamber between the point of origin and the end of the test tube to measure the estimated loss of sound pressure level in decibels (dB). In addition, a power amplifier function generator, sound level metres, and a velocity test probe were used in the present set-up.

Figure 2.

Calibration box.

Figure 3.

Calibrated amplitude and phase ratio of the microphones used in experimental tests.

Theory of the classical two-microphone technique for determining TL

The distance between the microphones and the sound source is the most influential factor in the calibration process. Therefore, it is enabled to minimize any distance error with respect to the variation of the frequencies and wavelength as well as to ensure no distinct phase shift among microphones is observed. Hence, the magnitude difference in phase shift is concluded as the result of the mounting process. To avoid such scenario, the transfer function of the reference and the target microphones was kept as the execution of relative calibration. Each microphone will have a range of frequency response or transfer function, and the frequency responses for the second and third microphones with respect to the first microphone are written as

\hat{H} ₁ ₂ (f) = \frac{\hat{S} ₁ ₂ (f)}{\hat{S} ₁ ₁ (f)}

(33a)

\hat{H} ₁ ₃ (f) = \frac{\hat{S} ₁ ₃ (f)}{\hat{S} ₁ ₁ (f)}

(33b)

Mathematical computation of the transfer function between the second and third microphones is expressed as

{\frac{\hat{S} ₂ ₃ (f)}{\hat{S} ₂ ₂ (f)}} = {\frac{\hat{S} ₁ ₃ (f)}{\hat{S} ₁ ₁ (f)}}

(34a)

\underset{c a l i b r a t e d}{\underset{⏟}{\hat{S} ₂ ₃ (f)}} = {\frac{\hat{S} ₁ ₃ (f) * \hat{H} ₁ ₂ (f)}{\hat{H} ₁ ₃ (f) * S ₁ ₂ (f)}} . {\frac{\hat{S} ₁ ₃ (f)}{{[\hat{H} ₁ ₂ (f)]}^{2}}}

(34b)

{\frac{\hat{S} ₂ ₃ (f)}{\hat{H} ₂ ₃ (f)}} = {\frac{\hat{S} ₁ ₃ (f) * \hat{S} ₂ ₂ (f)}{S ₁ ₂ (f)}} . {\frac{\hat{H} ₁ ₂ (f)}{\hat{H} ₁ ₃ (f) * {[\hat{H} ₁ ₂ (f)]}^{2}}}

(34c)

Therefore, the microphones can be calibrated using the frequency response function, with the first microphone being the reference. Raw data are collected, then normalized the ratio between the auto- and cross-spectral densities of the microphones to the first microphone at each frequency. From the equations below, n represents non-calibrated quantities and $\hat{S} ₁ ₁ (f)$ does not require correcting because the first microphone normalized the calibration. It is highlighted that the final calibration of the microphones is completed after processing the experimental data

{\hat{S}}_{11} (f) = {[{\hat{S}}_{23} (f)]}_{n}

(35a)

\hat{S} ₂ ₂ (f) = {\frac{{[{\hat{S}}_{22} (f)]}_{n}}{{[\hat{H} ₁ ₂ (f)]}^{2}}}

(35b)

\hat{S} ₃ ₃ (f) = {\frac{{[{\hat{S}}_{33} (f)]}_{n}}{{[\hat{H} ₁ ₃ (f)]}^{2}}}

(35c)

\hat{S} ₁ ₂ (f) = {\frac{{[{\hat{S}}_{12} (f)]}_{n}}{\hat{H} ₁ ₂ (f)}}

(35d)

\hat{S} ₂ ₃ (f) = {\frac{{[{\hat{S}}_{23} (f)]}_{n}}{\hat{H} ₂ ₃ (f)}}

(35e)

Calibration results of the microphones

Acoustic pressure amplitude data was logged using LabView software, and the experiments were conducted in a non-echoic environment to avoid influences from external factors that could impact data collection. Once data were collected, the data processing was done in ‘MATLAB’ software, which is particularly implemented for calibration. Recorded pressure amplitude is reasonably acceptable, although there is a slight difference in magnitude, that is, $\pm$ 2 dB reference to microphone 1, was noticed. There is still the possibility of the minimal phase shift of microphones; however, the recorded phase shift is considered negligible in the real-world scenario. Thus, the phase calibration is not as needed as the amplitude calibration.

The minimal phase shift difference is seen between $\pm$ 0.1 for microphones 1–6. However, further calibration is not required, since the collected phase shift is considered negligible and acceptable to perform the experiments.

Experimental set-up and tests

Overview of the experimental set-up and key test components

Coupled Helmholtz Resonators, shown schematically in Figure 4, are set up with the rectangular-shaped duct. There are two necking attached to the resonators and a membrane-embedded inside. Six microphone holding tubes are drilled, each 10 mm apart from the centre. The space for membrane mounting is available between two connecting resonators, and their necks are held tightly by bolts and nuts to prevent leakage. The chamber is connected to the boundary box, where the loudspeaker and air flow supply are mounted, as shown in Figure 4. The speaker is fixed at a place closest to the inlet duct, and the mean airflow is supplied into the chamber with a top acrylic lip, sealed with scotch tape to ensure that there is no air leakage during the experiments.

Figure 4.

Experimental set-up for Helmholtz resonators.

Loudspeaker and microphones

A loudspeaker is used as a source to transmit sound waves over input frequencies during the experiment. The speaker size of 6 inches is selected, that is, sufficient enough to generate power of 105 dB waveform, which covers the required frequency ranges. The loudspeaker used for calibration and assembled at the chamber is chosen approximately the same size to obtain the closest measurement. However, errors could be introduced if sound wave propagation along the air duct takes longer to develop the acoustic profile fully.

Six microphones are used for these experiments, and the sizes are selected to fit for the accuracy so that the minimum calibration would be acquired and the consideration of phase constraints in tolerances to the tubes drilled at the air duct’s wall. The axial locations of those six microphones should be mounted accurately, and careful handling of those microphones is necessary during the pre- and post-calibration. Figure 5(b) shows the microphones are slotted inside the tube sealed with blue tabs for experiments.

Figure 5.

Photos of a Loudspeaker (left) and a microphone (right)

Helmholtz Resonators

In this experiment as shown schematically in Figure 6, the resonators are $\frac{1}{2}$ shortening in neck length, and the individual resonator is calculated as

ω_{0} = c {[\frac{s}{L^{'} V}]}^{\frac{1}{2}}

(36)

where c = 343 m/s, S = 2.026x

10^{- 3} m^{2}

f_{1} = 140.2

Hz,

f_{2} = 115.7

ω_{o} = 2 π f_{1} = 880.9 \approx 881

ω_{o} = 2 π f_{2} = 726.96 \approx 727

Figure 6.

Geometrical Display of Square-Shaped Helmholtz Resonators.

The resonator’s $f_{1}$ is labelled as $H_{1}$ for the first resonator and $H_{2}$ is labelled as $f_{2}$ for the second resonator. Also, the resonator can be modelled as a single resonator with no partition in between, giving the resonant frequency of $f_{3}$ =89.2 Hz.

Acoustic foam and acoustic wool

Eggcrate foam as shown in Figure 7 is a type of acoustic foam formulated as soundproofing by attenuating the incident sound waves. The purpose of using acoustic foam is to enhance sound quality by reducing the residual sound in the environment. It is a lightweight material made of polyurethane foam and eliminates the noises by controlling the reverberation that sound can make by bouncing off from the walls. Reverberation is defined as several reflections, when a sound signal is reflected and decayed once the objects absorb the signal in the environment. Resulting from the sound reflection, the sound intensity starts to reduce depending on the frequency.

Figure 7.

Pictures of (a) acoustic foam and (b) acoustic wool.

Acoustic wool serves as sound-attenuating insulation. It has the properties of glass wool, so-called fibreglass, which primary function is reducing airborne sound transmission loss and sound reverberation in the environment. It is widely used in industry where lightweight, high tensile strength and exceptional resilience properties are required. It is made of fine and long fibre slabs with a diameter of 6–7 $μ m,$ each stacked by a binder at high temperature.

Flexible plate and rigid plate

The circular rig is fabricated, and an elastic sheet is pasted on it. The complete flexible plate as shown in Figure 8(a) is sandwiched between the two connecting rods of resonators. A similar approach has been carried out for rigid plates (see Figure 8(b)), where the material is made of acrylic embedded between resonators as another parameter for this experiment.

Figure 8.

Pictures of (a) flexible plate and (b) rigid plate as the sharable side wall of the coupled Helmholtz resonators.

Description of the experimental method and theoretical prediction of resonant frequencies

The experiment is conducted in the anechoic chamber with three different parameters to test the transmission loss of resonators. The experimental set-up of cylindrical-shaped Helmholtz resonators is shown in Figure 9. The first parameter is the rigid circular plate, 100 mm in diameter, connected in parallel between the two Helmholtz resonators to evaluate transmission loss. Secondly, the flexible circular plate of the same diameter with rubber attached, a damping mechanism, was placed to conduct the test. Lastly, the resonators’ necks are filled with three different types of the noise damping materials, such as acoustic foam and wool, for conducting the test. Two different geometries: square- and cylinder-shaped resonators, are experimentally tested. The corresponding results are discussed in section 5.

Figure 9.

Experimental set-up of cylindrical-shaped Helmholtz resonators in the anechoic chamber.

The resonator’s parameters and geometry of the overall structure are described in Appendixes 1 and 2. The loudspeaker is connected to the amplifier to produce a range of frequencies, 150–600 Hz, by increasing 5 Hz in each parameter. Six microphones method for transmission loss is adopted; measurement readings are recorded via three microphones located upstream and another three downstream of the resonators. Each microphone is fixed at 10 mm equal intervals away from each other, and the reference is chosen to microphone 1. Mean flow velocity is measured using a hot-wire anemometer in this experiment.

Four different parameters are tested, the air pump supply is turned on to the maximum supply voltage, and the maximum mean flow velocity is measured using a hot-wire anemometer. The subsequent velocities [ $0, \frac{2 u_{\max}}{3}, \frac{1 u_{\max}}{3}, u_{\max}$ ] are calculated, and the frequencies are collected using LabView and processed in MATLAB. Further dimensionless quantity of Mach number, $= \frac{u}{c}$ , which presents the ratio of the measured flow velocity, u, within the boundary to the local speed of the sound, c.

• Theoretical prediction of the resonant frequencies of HRs

The theoretical calculation is carried out for both square-shaped (see Figure 6) and cylindrical-shaped (see Figure 4) Helmholtz resonators to compare with the result from the experiment. The individual resonators frequency is calculated using the equation below

f = \frac{\sqrt{\frac{C^{2} S}{V L_{e f f e c t i v e}}}}{2 π}

(37)

where,

L_{e f f e c t i v e} = L + 1.2 R

is defined as effective length, equivalent to the corrective length. In order to capture the end effect, an additional length is employed at each end of the neck, since the acoustic length of the neck is longer than the physical neck length. An extra length with the value of 0.6 R is added to the length L, where the end side opens to the atmosphere, and has no flange. Two flanges are involved in each HR1 and HR2; therefore,

∆ L

=1.2 R is added to the physical length in this case.L= Physical Length of the Resonator’s Neck (m)R= Radius of the Resonator’s Neck (m)C= Speed of the Sound

(\frac{m}{s^{2}})

S= Surface Area of the Neck (

m^{2})

V= Volume of the Resonator m3

()

• Prediction of the resonant frequency of the cylindrical-shaped Resonators with a sharable Rigid Plate as a sidewall

f_{1} = \frac{\sqrt{\frac{C^{2} S}{V_{1} L_{e f f e c t i v e}}}}{2 π} = \frac{\sqrt{\frac{{(343)}^{2} * 1.964 * 10^{- 3}}{6.847 * 10^{- 4} * 0.11}}}{2 π} = 278.76 H z

{H z f}_{2} = \frac{\sqrt{\frac{C^{2} S}{V_{2} L_{e f f e c t i v e}}}}{2 π} = \frac{\sqrt{\frac{{(343)}^{2} * 1.964 * 10^{- 3}}{1.126 * 10^{- 3} * 0.11}}}{2 π} = 217.38

• Prediction of the resonant frequency of the square-shaped Resonators with a sharable Rigid Plate as a sidewall

f_{1} = \frac{\sqrt{\frac{C^{2} S}{V_{1} L_{e f f e c t i v e}}}}{2 π} = \frac{\sqrt{\frac{{(343)}^{2} * 1.964 * 10^{- 3}}{7.595 * 10^{- 4} * 0.10}}}{2 π} = 277.60 H z

f_{2} = \frac{\sqrt{\frac{C^{2} S}{V_{2} L_{e f f e c t i v e}}}}{2 π} = \frac{\sqrt{\frac{{(343)}^{2} * 1.964 * 10^{- 3}}{1.25 * 10^{- 3} * 0.10}}}{2 π} = 216.38 H z

Experimental results coupled cylindrical- and square-shaped Helmholtz Resonators with a sharable rigid Plate in the presence of a low mach number grazing flow

The experiment was carried out with a rigid plate made of acrylic with a 1.5 mm thickness circular in size and applied between the two resonators to act as a sharable sidewall. Ideally, two prominent peaks of the resonant frequencies could be observed. The experimental tests are conducted for a frequency range between 150 Hz and 600 Hz with an incremental step of 5 Hz.

Figures 10 and 11 show the variation of the sound absorption coefficient, Δ, of squared- and cylindrical-shaped coupled Helmholtz resonators varied with the forcing frequency, as the grazing flow Mach number is set to 4 different values. It can be seen that the maximum sound absorption coefficient could be as high as 80%, which means that the resonators absorb approximately 80% of the incident sound energy. While comparing the performance of square- and cylindrical-shaped resonators, it is observed that the cylindrical-shaped resonators perform better at the tested frequency range, that is, between 150 Hz and 600 Hz, with the presence of local peaks, Δ_max. This confirms that the cylindrical shape is a better design choice.

Figure 10.

Variation of sound absorption coefficient of square-shaped resonators with a sharable rigid plate.

Figure 11.

Variation of sound absorption coefficient of cylindrical-shaped resonators with a sharable rigid plate.

Observing both Figures 10 and 11 also reveals that the low Mach number grazing flow does not dramatically affect the sound absorption coefficient. This is probably because the grazing flow speed is negligible. Finally, the presence of the grazing flow does influence the resonant frequencies. The experimental measured resonant frequencies corresponding to the peaks of Δ are 20% different from the theoretically predicted value, as discussed in Section 3.5.

Coupled Helmholtz Resonators with no sharable plate applied

It would be interesting to know that when no sharable plate is applied to the coupled Helmholtz resonator, what would the aeroacoustics damping performance be changed? Figures 12 and 13 illustrate the variation of the sound absorption coefficient, Δ, of squared- and cylindrical-shaped coupled Helmholtz resonators varied with the forcing frequency. There is no sharable sidewall between the coupled resonators since the grazing flow Mach number is set to 4 different values. It can be seen from Figure 13 that the square-shaped coupled resonators have a large sound absorbing capacity near 500 Hz. In addition, the square-shaped resonators perform better than the cylindrical-shaped counterparts, as shown in Figure 12.

Figure 12.

Sound absorption coefficient of cylindrical-shaped resonators with no sharable plate.

Figure 13.

Variation of the sound absorption coefficient with frequency, as the square-shaped resonators are applied in the absence of any sharable plate.

Coupled Helmholtz resonators with a sharable flexible plate

The sharable rigid plate is replaced with a flexible plate so that the two resonators are mechanically coupled. The rubber membrane is secured on the circular metal plate with adhesive. To observe the acoustic damping performance, the flexible and rigid plates are tested to compare the transmission loss when there is no air velocity flow in the air duct. Figures 14 and 15 show the transmission loss variation with the frequency, as the grazing flow Mach number is set to 0 and 0.03, respectively, and the square-shaped Helmholtz resonators are concerned. A comparison is then made between the flexible sharable plate and the sharable rigid plate.

Figure 14.

Comparison of the transmission losses of square resonators in presence of the sharable rigid and the flexible plate at Ma = 0.0.

Figure 15.

Variation of the transmission loss comparison of square resonators with frequency, as the flexible plate is applied at Ma = 0.003.

It can be seen that the maximum transmission loss for the flexible rubber plate is 7.7 dB at the frequency of 200.1 Hz in the absence of the grazing flow, that is, Ma=0.00 (see Figure 14). When the grazing flow velocity increased to Ma=0.003 (see Figure 15), the peak transmission loss value of 16.1 dB is obtained at 264.9 Hz. It means that the flexible plate leads to a much more significant transmission loss compared to that of the rigid plate.

Coupled resonators filled with acoustic foam and wool

The experiment is carried out with noise damping materials such as acoustic foam and wool filled into the resonators’ necks. Since the necks are filled up, tests were conducted for a frequency range between 150 Hz and 1000 Hz under velocity flow and no-flow conditions. Figures 16 and 17 show the results of the cylindrical-shaped Helmholtz resonators, as the acoustic form and acoustic wool are inserted, and the grazing flow Mach number is set to 2 different values.

Figure 16.

Variation of the sound absorption coefficient of coupled cylindrical-shaped resonators filled with noise damping materials at Ma= 0.0 and 0.003.

Figure 17.

Transmission loss of the cylindrical-shaped coupled Helmholtz resonators filled with Noise Damping Materials at Ma= 0.0 and 0.003.

As seen below in Figure 16, acoustic foam and wool would give the transmission loss at 800.1 Hz, at the peak pressure amplitude of 11.4 dB. The plotted graph for absorption coefficient had some trade-offs between 780 Hz and 830 Hz in the test. Keeping the same neck’s diameter with a slight variation in the length of the neck up to 10 mm, the overall results for the no velocity flow condition would be the same. The sound transmission losses in the resonators are quite consistent for different damping materials.

Discussion of the case studies results

After conducting a series of experimental tests in our case study, there are two peaks of transmission loss for both rigid plates, the cylindrical- and square-shaped Helmholtz resonators, which agree well with the theoretical prediction. Additionally, the experimental results are found to be matched for both the cylindrical- and square-shaped resonators, when there is no plate shared between the two resonators. When there is no plate attached between the coupled resonators, it has been found that the distinct transmission loss peaks are located at higher frequencies than the natural resonant frequency of the resonators. For the acoustic damping materials, only one visible peak was seen before 1 kHz, that is, most likely due to the sound waves being reduced in the amplitude of pressure oscillations as the result of the dissipation in resonators. In the absence of the grazing flow, the flexible plate with the cylindrical-shaped resonators configuration resulted in a larger transmission loss than the square-shaped resonators. Further comparison of the rigid plate case of the cylindrical- and square-shaped coupled resonator in the absence of the grazing flow shows that the cylindrical-shaped resonators are performing better than square-shaped resonators in terms of the aeroacoustic damping performances.

Overview of numerical modelling of coupled Helmholtz resonators

Governing equations of the numerical model

The system equations of a pipe with an aeroacoustic damper applied consist of the mass, momentum and energy conservation.^1,8 Air is the working medium, and it is assumed to be a perfect gas. The systems equations of Navier–Stokes in the time-domain can then be described as

\frac{\partial ρ}{\partial t} + \frac{\partial}{\partial x_{i}} (ρ u_{i}) = 0

(38a)

\frac{\partial}{\partial t} (ρ u_{i}) + \frac{\partial}{\partial x_{j}} (ρ u_{i} u_{j}) + \frac{\partial p}{\partial x_{i}} = \frac{\partial τ_{i j} (u_{k})}{\partial x_{j}}

(38b)

\frac{\partial p}{\partial t} + u_{i} \frac{\partial p}{\partial x_{i}} + c^{2} ρ \frac{\partial u_{i}}{\partial x_{i}} = (κ - 1) τ_{i j} (u_{k}) \frac{\partial u_{i}}{\partial x_{j}} .

(38c)

Where u_i,

ρ

, and p mean the instantaneous velocity vector, density and pressure variables; c denotes the average sound speed;

τ_{i j} (u_{k}) = μ (\frac{\partial u_{i}}{\partial x_{j}} + \frac{\partial u_{j}}{\partial x_{i}} - \frac{2}{3} \frac{\partial u_{k}}{\partial x_{k}} δ_{i j})

is the viscous stress.

κ

is the air’s specific heat ratio. Here we assume that the heat and gravity force-related molecular diffusion is negligible (i.e. λ = 0 and g_i = 0).

The thermodynamic state equation also holds as $p = ρ R T$ . The compressible Navier–Stokes equations, that is, Equations (38a)-(38c) can be linearized to obtain the acoustic governing equations. The instantaneous variables are assumed to consist of a mean and a fluctuating part. In physics, it means that small-amplitude fluctuations (represented by a prime) are superimposed on a mean flow field (denoted by an overbar). In mathematics, it is written as follows

ρ = \bar{ρ} + ρ^{'}, u_{i} = {\bar{u}}_{i} + u_{i}^{'}, p = \bar{p} + p^{'}

(39a)

The fluctuating parts of the air density, flow velocity and acoustic pressure can be expressed in terms of Fourier series expansion as

ρ^{'} = \hat{ρ} {(ω) e}^{i ω t}, u_{i}^{'} = {\hat{u}}_{i} {(ω) e}^{i ω t}, p^{'} = \hat{p} (ω) e^{i ω t}

(39b)

These quantities, $ρ^{'} / \bar{ρ}$ , $| u_{i}^{'} | / \bar{c}$ , $p^{'} / \bar{p}$ and $T^{'} / \bar{T}$ , are assumed to be a small order ε, where ε<< 1 and ${\bar{c}}^{2} = κ R \bar{T}$ is the mean speed of sound. Note: the mean part of the flow variables satisfies the continuity of mass, momentum and energy. An overhat in equation (39b) denotes the flow parameters in the frequency domain.

Substituting Equations (39a) into (38a)–(38c) gives

\frac{\partial ρ^{'}}{\partial t} + \frac{\partial}{\partial x_{i}} ({\bar{u}}_{i} ρ^{'} + \bar{ρ} u_{i}^{'}) = 0

(40a)

\frac{\partial \bar{ρ} u_{i}^{'}}{\partial t} + \frac{\partial}{\partial x_{j}} (\bar{ρ} {\bar{u}}_{j} u_{i}^{'}) + (\bar{ρ} u_{j}^{'} + {\bar{u}}_{j} ρ^{'}) \frac{\partial {\bar{u}}_{i}}{\partial x_{j}} + \frac{\partial p^{'}}{\partial x_{i}} = \frac{\partial τ_{i j} (u_{k}^{'})}{\partial x_{j}}

(40b)

\frac{\partial p^{'}}{\partial t} + \frac{\partial}{\partial x_{i}} ({\bar{u}}_{i} p^{'} + κ \bar{p} u_{i}^{'}) + (κ - 1) (p^{'} \frac{\partial {\bar{u}}_{i}}{\partial x_{i}} + u_{i}^{'} \frac{\partial \bar{p}}{\partial x_{i}}) = (κ - 1) [{(τ_{i j} (u_{k}) \frac{\partial u_{i}}{\partial x_{j}})}^{'}]

(40c)

Equations (40a)–(40c) describe the spatiotemporal evolution of the fluctuating quantities: $ρ^{'}$ , $u_{i}^{'}$ and $p^{'}$ . Note that the pressure fluctuation is selected as a primitive variable in the energy equation equation (41a). Furthermore, the source term of equation (40c) is attributed to the molecular stresses. The linearized mass conservation equation (40a) is needed to further simplify the linearized momentum equation (40b). By conducting linearization, any second-order terms are negligible. These include the term $u_{i}^{'} u_{j}^{'}$ characterizing turbulence. Substituting equation (39b) into equations (40a)–(40c), and eliminating the common factor of $e^{i ω t}$ leads to the linearized Navier–Stokes equations in frequency domain as

i ω \hat{ρ} + \frac{\partial}{\partial x_{i}} ({\bar{u}}_{i} \hat{ρ} + \bar{ρ} {\hat{u}}_{i}) = 0

(41a)

i ω \bar{ρ} {\hat{u}}_{i} + \frac{\partial}{\partial x_{j}} (\bar{ρ} {\bar{u}}_{j} {\hat{u}}_{i}) + (\bar{ρ} {\hat{u}}_{j} + {\bar{u}}_{j} \hat{ρ}) \frac{\partial {\bar{u}}_{i}}{\partial x_{j}} + \frac{\partial \hat{p}}{\partial x_{i}} = \frac{\partial τ_{i j} ({\hat{u}}_{k})}{\partial x_{j}}

(41b)

i ω \hat{p} + \frac{\partial}{\partial x_{i}} ({\bar{u}}_{i} \hat{p} + κ \bar{p} {\hat{u}}_{i}) + (κ - 1) (\hat{p} \frac{\partial {\bar{u}}_{i}}{\partial x_{i}} + {\hat{u}}_{i} \frac{\partial \bar{p}}{\partial x_{i}}) = (κ - 1) [(τ_{i j} ({\hat{u}}_{k}) \frac{\partial {\hat{u}}_{i}}{\partial x_{j}})]

(41c)

These frequency-domain–governing equations (58)–(60) are iteratively determined via COMSOL 5.4. Here, UMFPACK (a set of routines for solving unsymmetric sparse systems) direct solver and Menter’s Shear Stress Transport (SST) turbulence model are applicable due to the non-negligible duct/grazing flow since it introduces much less numerical diffusion and leads to a better prediction of the eddy viscosity. The SST model is a combination of the k-ε model in the free stream and the k-ω model near the solid surface of the duct to capture the turbulence attenuation in details.

Model validation

The unstructured meshes are illustrated in Figure 18. A mesh independence investigation is conducted first. In order to better capture the vortices and the shear layer at the duct-neck joint sections, a finer mesh with 1,260,360 cells is chosen. The selected mesh is enough by benchmarking the present predictions with the available experimental data in Refs. [26–28]. The frequency step is set to 5 Hz and remains the same in all simulated cases. It will be confirmed that the frequency step is fine enough to produce good solutions. Note that the step could be decreased; however, a small frequency step gives rise to a significantly increased computational time and cost. Meanwhile, a smaller step is found to have a negligible effect on the transmission loss variation.^21,22 Therefore the frequency step is set to 5 Hz.

Figure 18.

Computational unstructured mesh of the 2D modelled duct with the modified Helmholtz.

In order to benchmark the developed model, comparison investigations are conducted on a 2D rectangular duct, as a grazing mean flow is present and a conventional Helmholtz resonator is implemented. Our numerical predictions are then compared with the experimental data available in Ref. [1], and the 3D simulation results of the cylindrical duct and the cylindrical resonator. Note that all variables in the current numerical studies are chosen to be the same as the experimental investigations, for example, the geometry and physical dimensions of the HR and the duct and the grazing flow conditions. To be consistent with the experimental measurements, the studied frequency range is set to be between 50 and 250 Hz. Figure 19 shows the comparison of transmission loss varied with the acoustic frequency between the present numerical, 3D time-domain simulation (STAR-CD Ref. [1]) and the experimental results, as the grazing flow Mach number Mu is set to 3 different values. It can be seen that for a given Mu, a good agreement between the numerical, 3D simulation and the experimental results is obtained in general. As the grazing flow Mach number is increased, the maximum transmission loss is found to be reduced dramatically. The corresponding resonant frequency is found to be shifted from 88 Hz to ∼100 Hz, as Mu is increased from 0 to 0.1. This investigation reveals that the developed numerical tool is applicable to evaluate the acoustic attenuation performances of a Helmholtz resonator,^34,35 as there is a grazing mean flow.

Figure 19.

Comparing the current predicted transmission losses with those obtained from the present model, 3D numerical results¹ and the experimental results from Ref. [1].

To simulate the experimental tests as reported in the case study, a 3D numerical model is developed and meshed, as shown in Figure 20. There are three boundary conditions specified to use in this model. Firstly, at the solid boundaries, the outer walls of the resonators and air duct are defined as the hard (wall) boundary conditions.

Figure 20.

Representation of mesh in cylindrical and square resonators.

Superposition of incoming imposed plane waves and exiting plane waves. This condition is validated, if the frequency is maintained below the cut off frequency for the second propagating mode. Also, it is important that enough meshes must be sized correctly to allow the proper resolution of pressure waves, since the mesh frequency cut off is related to the mesh size. The figure below shows the mesh for the cylindrical resonators comprises 36,907 elements in domain, 4456 elements in boundary, and 463 edge elements. In contrast, square resonators consist of 47,568 elements in domain, 5090 elements in boundary, and 445 edge elements. It is seen that the mesh in the cylindrical resonators is finer than of the square resonators. Mesh can be calculated by using [ $\frac{\frac{c}{f}}{10}$ ], where $λ = 10 \frac{c}{f}$ meshes per wavelength.

Other boundary conditions are named plane wave radiation at the inlet and outlet. The wave implies consistent pressure across the air duct in which it propagates. The plane wave can be seen as a rigid piston at the outlet and considered in a 2D situation, where boundary condition implies termination in specific acoustic impedance of the medium. Thirdly, the incident pressure field at the inlet with |p’|_max= 20 Pa plane waves.

After the previous discussions on three boundary conditions, the surface acoustic pressure field simulation is carried out, as shown in Figure 21, at 250 Hz; the sound levels for HR1 are quite close to one another. When the frequency is set to 500 Hz (see Figure 22), acoustic pressures appear to be similar for HR2 as well. Besides, the frequencies are selected to approximate values, and the results are then compared with the mathematically results.

Figure 21.

Representations of acoustic pressure field level for resonators at 250 Hz.

In general, with a proper validation, the frequency-domain model coupled be applied as a prediction tool to determine the acoustics damping performances of coupled Helmholtz resonators. It has great potential to be implemented in combustion-involved engines.^103–114 The acoustic damping performance of such coupled Helmholtz resonators is also related to the structures.^115–119 The conventional theoretical model^120,121 may not be applicable like the experimental^122,123 and numerical approaches^124,125 as discussed in the present work.

Conclusions and future works

In practice, coupled Helmholtz resonators are applicable to broaden the effective noise damping range.^126,130 In this work, we conduct a brief overview of the design and the evolution of coupled Helmholtz resonators. Helmholtz resonator design and optimization can be achieved by two typical approaches. One is the experimental test, and the other is numerical modelling. For this, the recent research and development of the experimental and numerical prediction methods are reviewed. With the proper calibration of acoustic pressure sensors, the experiments could be performed to determine the transmission loss profiles with respect to noise damping materials, the geometry of the resonators and dimensional variations in Helmholtz resonators. The numerical approach could also be applied in either time- or frequency domain. Comparison between the experimental and numerical predictions can then be conducted to validate the model and to optimize the design of coupled Helmholtz resonators.

In this paper, a case study of the transmission loss of geometrical differences in resonators, reduction of neck length, and partition formation between HR1 and HR2 with rigid, flexible plate and noise damping materials have been studied experimentally. It is observed that the mean flow velocity has an insignificant impact in operational ranges of 150–600 Hz regardless of the square- and cylindrical-shaped resonators. More transmission loss peaks are observed in the cylindrical type than the square type, as well as the geometry, which affects the noise damping performances. Besides, the volume of the resonator is inversely proportional to the resonant frequencies. Increasing the cavity volume capacity results in a decrease in the resonant frequency in this case study. The unexpected observation is a shift of the resonant frequency with no rigid plate to the higher value relative to that of the flexible plate and rigid plate when Ma=0.0.

Further experimental and numerical researches could be conducted. It has been shown in our case studies that the cylindrical resonator is experimentally confirmed to perform better in terms of transmission loss than the square resonator. The set-up is done with the air supply pump, which can go up to the maximum mean velocity of 3.2 $\frac{m}{s^{2}}$ . It is better to conduct details studies on the effects of (1) higher grazing flow Mach number and (2) bias flow Mach number on the aeroacoustics damping performances of the coupled Helmholtz resonators. It is also interesting to investigate the amount of the damping materials filled in resonators affecting the sound transmission loss over a broader frequency range. In addition, further numerical investigations especially by using 3D numerical model could be conducted to improve the design and to achieve the optimization of the coupled Helmholtz resonators, with detailed and convincing validation of the numerical model.

Footnotes

Acknowledgements

The first and corresponding authors would like to thank Zhejiang University of Science and Technology and Tongji University for providing the financial support. (The National Natural Science Foundation of China, 51776188).

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: the National Natural Science Foundation of China, 51776188; Natural Science Foundation of Zhejiang Province; LQ21E06001 and Natural Science Foundation of Zhejiang Province; LQ21E06003.

ORCID iD

Chenzhen Ji

Appendix 1

Appendix 2

References

Wang

Zhao

. Numerical analysis on combustion instabilities in end-burning-grain solid rocket motors utilizing pressure-coupled response functions. Aerospace Sci Tech 2020; 98: 105701.

Wang

Rubini

Qin

, et al. A model to predict acoustic resonant frequencies of distributed Helmholtz resonators on gas turbine engines. Appl Sci 2019; 9: 1419.

Selamet

Lee

. Helmholtz resonator with extended neck. The J Acoust Soc America 2003; 113: 1975–1985.

Griffin

Lane

Huybrechts

. Coupled Helmholtz resonators for acoustic attenuation. J Vibration Acoust 2001; 123(1): 11–17.

Cai

Zhao

Chan

Shahsavari

Tailoring reduced mechanisms for predicting flame propagation and ignition characteristics in ammonia and ammonia/hydrogen mixtures. Energy 2022; 260: 125090.

Guan

Zhao

Low

. Experimental evaluation on acoustic impedance and sound absorption performances of porous foams with additives with Helmholtz number. Aerospace Sci Tech 2021; 119: 107120.

Guan

Zhao

Ren

. Aeroacoustic attenuation performance of a Helmholtz resonator with a rigid baffle implemented in the presence of a grazing flow. Int J Aerosp Eng 2020; 2020: 1916239.

Dastourani

Bahman-Jahromi

. Evaluation of aeroacoustic performance of a Helmholtz resonator system with different resonator cavity shapes in presence of a grazing flow. J Aerospa Eng 2021; 34(5): 0402106.

Francés

Navarro-Fuster

Marini

, et al. Estimation of the speed of sound waves using a modular 3D printed Helmholtz resonator. Phys Edu 2021; 56: 055039.

10.

Zhao

Morgans

. Tuned passive control of combustion instabilities using multiple Helmholtz resonators. J Sound Vib 2009; 320: 744–757.

11.

Zhao

A’Barrow

Morgans

Carrotte

. Acoustic damping of a Helmholtz resonator with an oscillating volume. AIAA J 2009; 47: 1672–1679.

12.

Selamet

. Circular asymmetric Helmholtz resonators. J Acoust Soc America 2000; 107(5): 2360–2369.

13.

Zhao

Wang

. Geometric shapes effect of in-duct perforated orifices on aeroacoustics damping performances at low Helmholtz and Strouhal number. J Acoust Soc Am 2019; 145(4): 2126–2136.

14.

Guan

Zhao

, et al. Evaluations of acoustic damping performances of double-layer in-duct perforated plates at low Mach and Helmholtz number. J Acoust Soc America 2019; 146(5): 3512–3523.

15.

Guan

Zhao

, et al. Aeroacoustic damping performance studies on off-axial double-layer in-duct orifices at low Mach and Helmholtz number. Appl Acoust 2019; 156: 46–55.

16.

Jungowski

Grabitz

. Self-sustained oscillation of a jet impinging upon a Helmholtz resonator. J Fluid Mech 1987; 179: 77–103.

17.

Lee

Kang

Kim

, et al. Fabrication and performance evaluation of the Helmholtz resonator inspired acoustic absorber using various materials. Micromachines 2020; 11(11): 983.

18.

Selamet

Iqbal

, et al. Effect of Flow on Helmholtz Resonator Acoustics: A Three-Dimensional Computational Study vs. Experiments, 2011. SAE Technical Paper 2011-01-1521. DOI: 10.4271/2011-01-1521

19.

Chen

Zhao

. Numerical study of non-reacting flowfields of a swirling trapped vortex ramjet combustor. Aerospace Sci Tech 2018; 74: 81–92. DOI: 10.1016/j.ast.2018.01.006

20.

Sun

Zhao

Zhu

Rao

Wang

Large-eddy simulations of self-excited thermoacoustic instability in a premixed swirling combustor with an outlet nozzle. Physics of Fluids 2022; 34: 044112.

21.

Lavrentjev

Rämmal

. Experimental study of noise barrier boards with increased acoustic performance by utilizing Helmholtz resonator effects. Mater Today Proc 2020; 28(4): 2566–2571.

22.

Zhao

. Transient growth of flow disturbances in triggering a Rijke tube combustion instability. Combustion and Flame 2012; 159: 2126–2137.

23.

Zhang

Zhao

Dobriyal

, et al. Theoretical and experimental investigation of thermoacoustics transfer function. Appl Energ 2015; 154: 131–142.

24.

Benouhiba

Rougeot

Andreff

, et al. Origami-based auxetic tunable Helmholtz resonator for noise control. Smart Mater Structures 2021; 30: 035029.

25.

Zhao

. A review of acoustic dampers applied to combustion chambers in aerospace industry. Prog Aerosp Sci 2015; 74: 114–130.

26.

Zhao

, et al. A review of active control approaches in stabilizing combustion systems in aerospace industry. Prog Aerosp Sci 2018; 97(1): 35–60.

27.

Zhao

Yin

Experimental and investigation of geometric shape effect of coupled Helmholtz resonators on aeroacoustics damping performances in presence of low grazing flow. Aerospace Sci Tech 2022; 128: 107799.

28.

Zhao

Becker

Entropy production and enhanced thermal performance studies on counter-flow double-channel hydrogen/ammonia-fuelled micro-combustors with different shaped internal threads. Int. J. Hydrogen Energy 2022; 47(85): 36306–36322.

29.

Cai

Zhao

. Temperature dependence of laminar burning velocity in ammonia/dimethyl ether-air premixed flames. J. Thermal Science 2022; 31: 189–197.

30.

Zhao

Chow

. Thermoacoustic instability of a laminar premixed flame in Rijke tube with a hydrodynamic region. J Sound Vib 2013; 332: 3419–3437.

31.

Zhao

, et al. NOx Emission Reduction Reaction of Ammonia-Hydrogen with Self-sustained Pulsating Oscillations. Therm Sci Eng Prog 2020; 19: 100615.

32.

Dowling

Hughes

. Sound absorption by a screen with a regular array of slits. J Sound Vib 1992; 156: 387–405.

33.

Song

Zhu

Pan

Wang

Zhao

Numerical investigations on the beating behavior of self-excited combustion instability in a hydrogen-fueled Rijke type combustor. Aerospace Sci Tech 2022; 126: 107624.

34.

Mahesh

Ranjith

Mini

. Inverse design of Helmholtz resonator based low-frequency acoustic absorber using deep neural network. J Appl Phys 2021; 129: 174901.

35.

Euve

Pham

Petigjeans

, et al. Time domain modelling of a Helmholtz resonator analogue for water waves. J Fluid Mech 2021; 920: A22.

36.

Mahesh

Mini

. Theoretical investigation on the acoustic performance of Helmholtz resonator integrated microperforated panel absorber. Appl Acoust 2021; 178(1): 108012.

37.

Jun

Nespesny

Pencik

, et al. Optimized method for Helmholtz resonator design formed by perforated boards. Appl Acoust 2021; 184(1): 108341.

38.

Seybert

Ross

. Experimental determination of acoustic properties using a two-microphone random-excitation technique. J Acoust Soc Am 2010; 61: 1362–1370.

39.

Zhao

. Two-dimensional lattice Boltzmann investigation of sound absorption of perforated orifices with different geometric shapes. Aerosp Sci.Technol 2014; 39: 40–47.

40.

Rao

Zhao

Wang

Zhao

Xie

. Similarity phenomena of lean swirling flames at different bulk velocities with acoustic disturbances. Chinese Journal of Aeronautics. DOI: 10.1016/j.cja.2022.07.001

41.

Zhao

Guan

Reinecke

. Characterizing hydrogen-fuelled pulsating combustion on thermodynamic properties of a combustor. Commun Phys 2019; 2: 44.

42.

Zhao

. Feedback control of combustion instabilities using a Helmholtz resonator with an oscillating volume. Combust Sci Technol 2012; 184(5): 694–716.

43.

Singh

Rienstra

. Nonlinear asympotoic impedance model for a Helmholtz resonator liner. J Sound Vibration 2014; 333(15): 3536–3549.

44.

Zhao

. Transmission loss analysis of a parallel-coupled Helmholtz resonator network. AIAA J 2012; 50(6): 1339–1348.

45.

Tam

CKW

Walker

. Numerical simulation of a slit resonator in a grazing flow under acoustic excitation. J Sound Vibration 2008; 313(3–5): 449–471.

46.

Cummings

Eversman

. High amplitude acoustic transmission through duct terminations: theory. J Sound Vib 1983; 91: 503–518.

47.

Rossing

. Springer Handbook of Acoustics. 2nd Edition. New York: Springer-Verlag, 2014.

48.

Mendez

Eldredge

. Acoustic modeling of perforated plates with bias flow for large-eddy simulations. J Comput Phys 2009; 228: 4757–4772.

49.

Tonon

Moers

EMT

Hirschberg

. Quasi-steady acoustic response of wall perforations subject to a grazing-bias flow combination. J Sound Vib 2013; 332: 1654–1673.

50.

Burak

Billson

Eriksson

, et al. Validation of a time- and frequency-domain grazing flow acoustic liner model. AIAA J 2009; 47(8): 1841–1851.

51.

Sun

Cai

Shahsavari

Sun

Zhao

Wang

. RANS simulations on combustion and emission characteristics of a premixed NH3/H2 swirling flame with reduced chemical kinetic model. Chinese Journal of Aeronautics 2021; 34(12): 17–27.

52.

Zhang

Bodony

. Numerical simulation of two dimensional acoustic liners with high-speed grazing flow. AIAA J 2011; 49(2): 365–375.

53.

Davies

Grove

Bell

Rea

Furkert

Sellier

. Preliminary design and test of high altitude two-stage rockets in new zealand. Aerospace Sci Techno 2022; 128: 107741.

54.

Jing

Sun

, et al. Effect of grazing flow on the acoustic impedance of an orifice. AIAA J 2001; 39(8): 1478–1578.

55.

Zhao

Ang

. Numerical and experimental investigation of the acoustic damping effect of single-layer perforated liners with joint bias-grazing flow. J Sound Vibration 2015; 342: 152–167.

56.

Zhao

Ang

. Experimental comparison of noise dissipation effects of single- and double-layer acoustic liners. Appl Acoust 2018; 141: 281–292.

57.

Youcef

Saim

. Numerical Analysis of the Baffles Inclination on Fluid Behavior in a Shell and Tube Heat Exchanger. J Appl Comput Mech 2021; 7(1): 312–320. DOI: 10.22055/jacm.2020.32925.2103

58.

Zhao

NOx emission reduction in ammonia-powered micro-combustors by partially inserting porous medium under fuel-rich condition. Chemical Engineering Journal 2022; 434: 134680.

59.

Tian

Ain

Anjum

, et al. Fractal N/MEMS: From pull-in instability to pull-in stability. Fractals 2021; 29(2): 2150030.

60.

Anjum

Skrzypacz

. A variational principle for a nonlinear oscillator arising in microelectromechanical system. J Appl Comput Mech 2021; 7(1): 78–83.

61.

. Transmission Loss Measurement of Coupled Helmholtz Resoantos Filled with Noise Damping Materials. Nanyang Technological University, 2015.

62.

Kierkegaard

Boji

Efraimsson

. A frequency domain linearized Navier-Stokes equations approach to acoustic propagation in flow ducts with sharp edges. J Acoust Soc Am 2010; 127(2): 710–719.

63.

Zhao

Han

, et al. Parametric measurements of the effect of in-duct orifice edge shape on its noise damping performance. J Sound Vibration 2016; 384: 130–145.

64.

Zhang

Zhao

Han

. Control of combustion instability with a tunable Helmholtz resonator. Aerospace Sci Tech 2015; 41: 55–62.

65.

Zhao

. Lattice Boltzmann investigation of acoustic damping mechanism and performance of an in-duct circular orifice. J Acoust Soc America 2014; 135(6): 3243–3251.

66.

Zhao

. Two-dimensional Lattice Boltzmann Investigation of Sound Absorption of Perforated Orifices with Different Geometric Shapes. Aerospace Sci Tech 2014; 39: 40–47.

67.

Zhao

. Mitigation of premixed flame-sustained thermoacoustic oscillations using an electrical heater. Int J Heat Mass Transfer 2015; 86: 309–318.

68.

Zhao

. Combustion Instabilities in a Bifurcating Tube: Open- and Closed-loop Measurements. AIAA J 2014; 52(11): 2513–2523.

69.

Zhao

Morgans

Dowling

. Tuned Passive Control of the Acoustic Damping of Perforated Liners. AIAA J 2011; 49(4): 725–734.

70.

Zhao

. Feedback control of combustion instabilities using a Helmholtz resonator with an oscillating volume. Combustion Sci Tech 2012; 184(5): 694–716.

71.

Zhong

Zhao

. Time-domain characterization of acoustic damping of a perforated liner with bias flow. J Acoust America Soc 2012; 312(1): 271–281.

72.

Guo

Zhao

Becker

. Temperature difference and stack plate spacing effects on thermodynamic performances of standing-wave thermoacoustic engines driven by cryogenic liquids and waste heat. J. Thermal Science 2022; 31(5): 1434–1451.

73.

Zhao

. Transmission loss analysis of a parallel-coupled Helmholtz resonator network. AIAA J 2012; 50(6): 1339–1346.

74.

Sun

Zhao

Cai

. Dilution effect on thermal performances and NOx emission characteristics of an ammonia-oxygen micro-thermophotovoltaic system. Thermal Science and Engineering Progress 2022; 34: 101401.

75.

Liang

Zhao

. Identifying the skeleton mechanism for oscillatory combustion with functional weight analysis. Combustion and Flame 2022; 244: 112243.

76.

Zhao

. A Real-time Plane-wave Decomposition Algorithm for Characterizing Perforated Liners Damping at Multiple Mode Frequencies. J Acoust Soc America 2011; 129(3): 1084–1192.

77.

Zhao

. Proper orthogonal and dynamic model decomposition analyses of nonlinear combustion instabilities in a solid-fuel Ramjet combustor. Therm Sci Eng Prog; 27(2022): 101147.

78.

Zhao

. Numerical investigations on solid-fueled Ramjet inlet thermodynamic properties effects on generating self-sustained combustion instability. Aerospace Sci Tech; 119(2021): 107079.

79.

Dandsena

Mohapatra

Satapathy

, et al. Noise Control of outdoor unit of split type air-conditioner using periodic scatters made with array of Helmholtz resoantors. Appl Acoust 2021; 179: 108054.

80.

Rajendran

Piascek

Echenagucia

. Design of broadband Helmholtz resonator arrays using the radiation impedance method. J Acoust Soc Am 2022; 457: 151.

81.

Gautam

Celik

Meek

, et al. On the acoustic performance of double degree of freedom Helmholtz resonator based acoustic liners. Appl Acoust 2022; 191: 108661.

82.

Ismail

Kim

Chang

, et al. Sound transmission loss of a Helmholtz resonator-based acoustic metasurface. Appl Acoust 2022; 188: 108569.

83.

Thannickal

Tharakan

Chakravarthy

. Active control of combustion noise by a twin resonator trim adjustment system. Combust Sci Tech 2021; 194: 16. DOI: 10.1080/00102202.2021.1929195

84.

Asakura

. Numerical investigation of the sound insulation effect of a suspended ceiling structure with arrayed Helmholtz resonators by the finite-difference time-domain method. Appl Acoust 2021; 172: 107601.

85.

Yang

Choy

Wang

, et al. Sound radiation and suppression of an unbaffled long enclosure using Helmholtz resonators. Mech Sys Sig Proc 2022; 165: 108408.

86.

Mariia

Anton

Sergei

, et al. Broadband Noise Mitigation Using Coupled Helmholtz Resonators: A Numberical Study. Washington DC, USA: InterNoise2021, pp. 3944–4918.

87.

Pishvar

Harne

R. L.

Nonlinear behaviro of Helmholtz resonator with a compliant wall for low-frequency broad band noise control. J Vib Acoust 2022; 144(3): 031008.

88.

Yang

. Investigation of thermoacoustic oscillation attenuation by modified Helmholtz dampers with dual frequency bands. Appl Acoust 2022; 185: 108433.

89.

Guo

Zhang

Fang

, et al. Wideband low-frequency sound absorption by inhomogeneous multi-layer resonators with extended necks. Compos Structures 2021; 260: 113538.

90.

Gautam

Celik

Azarpeyvand

. An experimental and numerical study on the effect of spacing between two Helmholtz resonators. Acoustics 2021; 3(1): 97–117.

91.

Zhao

Yang

Zhang

. Numerical investigation of the effect of distributed heat sources on heat-to-sound conversion in a T-shaped thermoacoustic system. Energy 2015; 144: 204–213.

92.

Shi

Mak

. Helmholtz resonator with a spiral neck. Appl Acoust 2015; 99: 68–71.

93.

Langfeldt

Hoppen

Gleine

. Broadband low-frequency sound transmission loss improvement of double walls with Helmholtz resonators. J Sound Vib 2020; 476: 115309.

94.

Tang

. On Helmholtz resonators with tapered necks. J Sound Vib 2005; 279: 1085–1096.

95.

Komkin

Mironov

Bykov

. Sound absorption by a Helmholtz resonator. Acoust Phys 2017; 63: 385–392.

96.

Chanaud

. Effects of geometry on the resonance frequency of Helmholtz resonators. J Sound Vib 1994; 178: 337–348.

97.

Zhang

Mak

, et al. Hybrid noise control using multiple Helmholtz resonator arrays. Appl Acoust 2019; 143: 31–37.

98.

Nagaya

Hano

Suda

. Silencer consisting of two-stage Helmholtz resonator with auto-tuning control. J Acoust Soc Am 2001; 110: 289–295.

99.

Langfeldt

Hoppen

Gleine

. Resonance frequencies and sound absorption of Helmholtz resonators with multiple necks. Appl Acoust 2019; 145: 314–319.

100.

de Bedout

Franchek

Bernhard

, et al. Adaptive-passive noise control with self-tuning Helmholtz resonators. J Sound Vib 1997; 202: 109–123.

101.

Coulon

Atalla

Desrochers

. Optimization of concentric array resonators for wide band noise reduction. Appl Acoust 2016; 113: 109–115.

102.

Dowling

Ffowcs Williams

. Sound and Sources of Sound. New York: John Wiley & Sons, pp. 13–31.

103.

Reyhanoglu

. Feedback control of acoustic disturbances transient growth in triggering thermoacoustic instability. J Sound Vibration 2014; 333(16): 3639–3656.

104.

Yang

Sun

. Effect of distributed heat sournces on low frequency thermoacoustic instabilities. J Sound Vibration 2013; 332(12): 3098–3111.

105.

Zhao

. Prediction of stability behaviors of longitudinal and circumferential eigenmodes in a choked thermoacoustic combustor. Aerospace Sci Tech 2015; 46: 12–21.

106.

Zhao

, et al. Nonorthogonality analysis of a thermoacoustic system with a premixed V-shaped flame. Energ Convers Manag 2014; 85: 102–111.

107.

Zhao

Gutmark

Reinecke

. Mitigating self-excite flame pulsating and thermoacoustic oscillations using perforated liners. Sci Bull 2019; 64(13): 941–952.

108.

Cai

Zhao

Wang

. Evaluation of NOX emissions characteristics in a CO2-free micro-power system by implementing a perforated plate. Renew Sust Energ Rev 2021; 145: 111150.

109.

Huang

Chen

, et al. Effects of heat recirculation on combustion characteristics of n-heptane in micro combustors. Appl Therm Eng 2016; 109: 697–708.

110.

Duan

Zhang

, et al. Experimental study on the energy flow of a gasoline-powered vehicle under the NEDC of cold-starting. Appl Thermal Eng 2017; 115: 1173–1186.

111.

Zhao

Teo

. Perforemance of small scale bladeless electromagnetic energy harvesters driven by water or air. Energy 2014; 74: 99–108.

112.

Huang

Yan

, et al. Experimental study of n-heptane/air combustion in meso-scale burners with porous media. Experiemtnal Therm Fluid Sci 2014; 52: 47–58.

113.

Zhao

. Minimizing transient energy growth of nonlinear thermoacoustic oscillations. Int J Heat Mass Trans 2015; 81: 188–197.

114.

Cai

Zhao

Shi

, et al. Mitigating NOX emissions from an ammonia-fuled micro-power system with a perforated plate implemented. J Hazard Mater 2021; 401: 123848.

115.

Zhao

. Entropy-involved energy measure study of intrinsic thermoacoustic oscillations. Appl Energ 2016; 177: 570–578.

116.

Dezvarch

. Upgrading the seismic capacity of pile-supported wharfs using semi-active liquid column gas damper. J Appl Comput Mech 2020; 6(1): 112–124.

117.

Emami

Dezvareh

Mousavi

. Effectiveness of viscous damper on vibrations control of full-height lattice tower under wave and wind excitations in time domain. J Oceanograph 2022; 13(51): 1–13.

118.

Wang

Zhao

. Experimental study of transition to instability in a Rijke tube with axially distributed heat source. Int. J. Heat Mass Transfer 2022; 183(B): 122157.

119.

Tian

. Fractal Pull-in Stability Theory for Microelectromechanical Systems. Front Phys 2021; 9: 606011.

120.

Zhao

. Effect of heat source on transient energy growth analysis of a thermoacoustic system. Energ Convers Manag 2015; 89: 309–317.

121.

Zhao

. Feedback control of self-sustained nonlinear combustion oscillations. J Eng Gas Turbines POWER 2016; 138(6).

122.

Sun

Zhao

Rao

Wang

. Characterizing nonlinear dynamic features of self-sustained thermoacoustic oscillations in a premixed swirling combustor. Appl Energ 2020; 264: 114698.

123.

Zhang

Zhao

Dobriyal

. Theoretical and experimental investigation of thermoacoustic transfer function. Appl Energ 2015; 154: 131–142.

124.

Han

Zhao

. Effect analysis on energy conversion enhancement and NOx emission reduction of ammonia/hydrogen fuelled wavy micro-combustor for microthermophotovoltaic application. Fuel 2020; 289: 119755.

125.

Gao

Wang

Cheng

. Deep auto-encoder network in predictive design of Helmholtz resonator: on-demand prediction of sound absorption peak. Appl Acoust 2022; 191: 108680.

126.

Zhu

Gao

, et al. Nonlinear sound absorption of Helmholtz resonators with serrated necks under high-amplitude sound wave excitation. J Sound Vibration 2022; 537: 117197.

127.

Liu

Xin

. Sound absorption of a perforated planel backed with perforated porous material: Energy dissipation of Helmholtz resonator cavity. Mechancial Syst Signal Process 2023; 185: 109762.

128.

Sabat

Pennec

Leveque

, et al. Single and coupled Helmholtz resoantors for low frequency sound manipulation. J Appl Phys 2022; 132: 064505.

129.

Song

Huang

Zhou

, et al. Perfect acoustic absorption of Helmholtz resoantors via tapered necks. Appl Phys Express 2022; 15(8): 084006.

130.

Smith

MJA

Abrahams

. Tailored acoustic metamaterials Part I: Thin- and thick-walled Helmholtz resonator arrays. Proc R Soc A 2022; 478: 20220124.