1. Introduction
A major problem in today’s industrialized world is noise and its effects on people’s health. Various solutions have been proposed to solve this problem using noise reduction techniques, which are applied at three different stages: near the source, at the transmission, and at the receiver. The most common solution, and the most widely used because of its simplicity, is the use of devices at the transmission level or, more precisely, the installation of noise barriers (NBs) between the source and the receiver [
1]. The insulation performance of NB can be improved by a suitable choice of shape: the geometry of the edge or the curvature of the main flat surface can be designed to modify the diffraction of the upper part [
2]. Incorporating porous materials into them also helps to improve their acoustic perfor-mance, especially with regard to the reflection of sound waves between the source and body of the NB [
3].
In this field, a new method has emerged in the last decades, which has new properties in contrast to conventional flat walls for noise reduction: sonic crystal noise barriers (SCNB) [
4,
5,
6,
7,
8]. SCNBs are made of periodic materials called sonic crystals, which are arrays of scatterers embedded in air [
9]. While conventional barriers do not allow air or water to pass through, SCNBs are permeable to fluids. Moreover, they can be arranged to attenuate noise in specific target frequency ranges, which is one of their advantages over conventional barriers that cannot distinguish between different noises. SCNBs exhibit the property of Bragg bandgaps (Bragg-BGs), with frequency ranges for which acoustic waves are not able to propagate due to the interference generated by the multiple scattering of waves [
10,
11,
12]. The central frequency of the first Bragg-BG is f
c = c/2a, where c is the speed of sound in the fluid (m/s), and the lattice constant (a) is the distance between scatterers and is characteristic of the unit cell of a crystal lattice [
13]. Bragg-BGs have improved the ability to control wave propagation through many additional applications, such as waveguides [
14] or diffusers [
15,
16]. Some approaches have been developed to improve the insulation level and range of Bragg-BGs, e.g., by using fractal geometries [
17] or different radii of the scatterers that make up the SCs [
18].
In addition to Bragg scattering, BGs are also produced by local resonances. In particular, they are used to induce low-frequency BGs in SCs and have been studied in one- and two-dimensional periodic arrays [
19,
20,
21,
22,
23]. In acoustic systems, these resonances are usually generated by Helmholtz Resonators (HRs) [
24]. The BGs based on local resonances of HRs (HR-BGs) depend on the alignment of the mouth of the resonator with respect to the incident sound wave. Some studies have shown that an opposite arrangement of the HRs mouths can help to achieve sound insulation and air ventilation simultaneously [
25], and that the maximum of the absorption coefficient is extremely angle dependent [
26]. The HR-BG has a smaller frequency range than the Bragg-BG, but several authors have proposed to extend the BGs of SCs by multiresonant scatterers [
27,
28,
29,
30,
31,
32]. HR-BG and Bragg-BG influence each other [
33,
34] and sometimes the interaction between the two wave propagation phenomena does not improve the insulation capacity of the SCNB. [
35]. As far as we know, there are no studies on the effect of this interaction on the transmission of a SCNB and on the specific influence of the resonator on the Bragg-BG. The main purpose of this work is to study the effect of the interaction between the local resonance and the Bragg interference in order to increase the noise insulation produced by SCNBs.
2. Materials and Methods
The development and validation of innovative solutions for noise reduction usually requires an enormous amount of time, specific laboratory resources and high costs with little scope for change until a new prototype is manufactured. To circumvent these drawbacks, the scientific community has developed and validated several numerical algorithms to evaluate the acoustic performance of new devices, as it is crucial to accurately predict their acoustic properties even before they are designed. These simulation techniques have contributed to the advancement of acoustic technology, especially to the development of novel periodic materials, such as SCs. These numerical simulation techniques, in combination with optimization strategies, are efficient design tools as they allow the development of devices with improved acoustic performance and new features [
15,
16,
36].
The periodicity and symmetries of SCs allow the calculation of the model to be simplified and to achieve an approximation to the real case of long domains with high accuracy, e.g., the evaluation of the insulation solution for an extended path, such as highways or railways. In recent years, a variety of simulation methods have been used to evaluate the effectiveness of periodic structures in acoustics. In the present work, domain-discretized techniques are used for numerical simulations.
2.1. The Finite-Difference Time-Domain (FDTD) Method
The FDTD method uses finite differences as approximations for the spatial and temporal derivatives of the fundamental equations. The method was first proposed by K. Yee for electromagnetics [
37]. Maloney and Cummings [
38] adapted the method to the field of acoustics by converting the momentum conservation and continuity equations into central difference equations, thus obtaining updated formulas for acoustic pressure and particle velocity. The seminal work of Cao et al. showed how efficient FDTD works for band structure calculations [
39], and Miyashita focuses on the study of waveguides based on SCs [
14]. FDTD has also been used in recent work on sonic crystals [
15,
16,
40,
41,
42]. According to the work of Peiró-Torres et al. [
43], FDTD is not the best method to determine the insertion loss of a SC due to the ratio between computational cost (memory) and accuracy. The high memory requirement is the main disadvantage of volumetric methods, such as FEM and FDTD. However, in terms of computation time, FDTD stands out from other methods, such as BEM and FEM, because it is a time-domain method that can cover a large frequency range with a single simulation. To ensure the stability of the method, the Courant–Friedrichs–Lewy number must not be greater than 1 [
44]. Moreover, this number should be as high as possible to minimize numerical dispersion, i.e., numerical error in phase velocity.
As with other numerical methods, in FDTD there is a limit to the maximum frequency that can be correctly simulated, which is related to the size of the mesh elements, also known as the spatial step. The typical rule of thumb is 8 elements per wavelength. However, since the numerical model, in this case, contains circular geometries (cylindrical scatterers), it is preferable to define finer meshes to avoid the need to apply conformal methods to minimize the staircase effect. In this work a mesh with square structure has been used.
2.2. Numerical Model
We consider a two-dimensional SC, which consists of a set of HRs distributed in an infinite square periodic lattice.
Figure 1a shows the unit cell of the two-dimensional periodic system under investigation, which is defined by a scatterer with the outer radius R and a HR inside. Periodic boundary conditions are set at the X and Y boundaries so that the infinitely extended SC is defined by a periodic structure with a square lattice and a lattice constant (a) defining the unit cell. The irreducible Brillouin zone is also shown and represents the central cell of the reciprocal lattice. In a two-dimensional system, the filling fraction is the ratio between the area occupied by the scatterer and the area of the entire unit cell. The resonator is characterized by the length of the neck (
L), the width of the mouth (
w), and the inner surface area (S =
), where
r is the inner radius of the surface of the HR. The first resonance frequency of the resonator is [
24]:
where ∆ is the correction factor for neck length added account for the effects of the neck-cavity junction caused by a sudden change in the circular cross-sections (usually between 1.6 and 1.8) [
45]. For simplicity, the geometry of the resonator for all the simulations in this work was chosen to neglect the viscothermal losses in the medium (see
Appendix A).
The band structure describes the dispersion relation of waves in the SC. FDTD simulations were used to calculate the band structure. In the band structure, propagation bands and BGs can be identified. In this work, we are interested in two propagation directions: 0° and 90°. The band structure is calculated in the two corresponding directions in reciprocal space or k-space: 0° (ΓX) and 90° (ΓX’), where G (k
x = 0, k
y = 0), X (k
x = π/a, k
y = 0) and X’ (k
x = 0, k
y = π/a) denote the boundaries of the first Brillouin zone. Numerical simulations using the FDTD method are performed to evaluate the transmission of two SCNBs with a design based on the unit cell defined in
Figure 1a. Periodic Boundary Conditions (PBC) are considered to reproduce an infinite SC in the Y-direction but truncated to three layers of scatterers in the X-direction (see
Figure 1b,c). Since our SC has square symmetry, the periodicity is the same in the
x- and
y-axis, and in the two transmission models, both the distance between the scatterers and the transverse dimension (
y-axis) of the models are a = 17 cm, with a filling fraction ff = π (r/a)2 = 0.4. This means that the external radius (R) is 6.86 cm and the distance between the scatterers is a/2-R = 1.64 cm.
Thus, the geometric model of the SCNBs and the unit cell have the same topology, namely a two-dimensional square-periodic structure. An incident plane wave with frequency f and acoustic pressure P
i propagates along the
x-axis from left to right in such a way that it is aligned with the mouth of the HR in
Figure 1b and perpendicular to it in
Figure 1c (hereafter referred to as the 0° and 90° alignments). Perfectly Matched Layers (PMLs) are applied to both
x-axis boundaries in the models to absorb waves propagating outward from the domain.
The transmission of the system is analyzed in the evaluation points of a line at a distance of 1 m from the SC by evaluating the transmission coefficient (
τ) defined by:
where <·> stands for the averaging of both, the squared incident acoustic pressure and the squared transmitted acoustic pressure (
Pt) at the evaluation points. From
the equation for the Insertion Loss (IL) in decibels can be derived as follows:
and indicates how the sound pressure level at the evaluation points changes when the SCNB is positioned in the acoustic path of the incident wave.
3. Bragg and Helmholtz Resonances Interaction
The effect of multiple scattering interference significantly reduces the transmission of a SCNB in the BGs and increases its insulation. In general, this phenomenon applies not only to the Bragg-BG, but also to the HR-BG [
20]. In the model proposed here, the scatterers have a local resonator. Depending on the central frequencies and the width of their BGs, both Bragg-BG and HR-BG can overlap, leading to a different type of interaction between the phenomena [
33,
34].
Next, the transmission of an incident plane wave from a SCNB is investigated numerically using the models shown in
Figure 1b for the 0° alignment and
Figure 1c for the 90° alignment. The aim is to analyze the physical interaction between the first Bragg resonance and the first local resonance and their mutual interaction. The value of the specific central frequency of the HR-BG influences the width and the insulation of the Bragg-BG. Two SCNBs are studied with their HR-BG below and above the Bragg-BG with a fixed lattice constant (a = 0.17 m), i.e., a central Bragg frequency of fc = 1000 Hz. According to (1), the HR resonance frequency depends on the neck, the mouth, and the surface of the two-dimensional resonator. The resonant frequency of the HR could be changed by suitably varying each of these parameters according to (1). However, due to geometric constraints in the design of the resonator inside the scatterer, it is appropriate to vary the surface of the resonator to cover a wide range of HR resonant frequencies with realistic geometries.
Figure 2 shows results of FDTD simulations of the transmission coefficient of a SCNB in a solid line based on the models shown in
Figure 1b,c together with the band structure of the SC, with circles with 0° alignment in red and 90° alignment in blue for ΓX and ΓX’, respectively. For reference, the SCNB is shown with Closed Cylinders, SCNB-CC, i.e., without local resonance, in black. The scatterers of the SCNB have a HR with a mouth w = 1 cm, a neck L = 1 cm, and an r = 4.5 cm, corresponding to a resonant frequency f
HR = 400 Hz < fc. The frequency range includes the first Bragg-BG with the central frequency f
c and the first HR-BG, with resonance frequency f
HR. Both phenomena have higher resonances, but they are far from interacting with the first resonances and their mutual interactions are beyond the scope of this work.
Although the band structure was obtained from an infinite periodic SC and the transmission loss from the SCNB with an array of three layers of resonators, the agreement between the presence of BGs and transmission minima is significant.
Figure 2-left shows the case of a HR with a small surface and, thus, a resonant frequency lower than the central frequency of the Bragg-BG (f
HR < fc). At low frequencies, below and around the local resonance, oscillations are observed in the transmission, which are due to Fabry–Perot interference phenomena [
46] due to the multiple reflections of the waves at the truncated ends of the SC. This is thus a phenomenon caused by the finite length of the crystal and cannot be explained by the band structure of the infinite SC. The interaction between the two BGs is visible in the frequency range of the HR-BG below the Bragg-BG, where the greatest difference from the case of the SCNB-CC is observed. A significant reduction in transmission is observed for both 0° alignment and 90° alignment due to the local resonance.
Figure 2-right shows the case where the resonance frequency of the HR has been tuned above the Bragg-BG (higher surface). The scatterers of the SCNB have a HR with a mouth w = 1 cm, a neck L = 1 cm, and r = 1.2 cm, corresponding to f
HR = 1400 Hz > fc. Again, the reduction of the transmission coefficient can be observed in the frequency range where the local resonance and the Bragg resonance interact. In
Supplementary Materials an animation shows the interaction of Bragg diffraction and the local resonance with the f
HR in the range from 400 Hz to 1400 Hz.
Figure 3 shows the numerical results of IL for several SCNBs consisting of arrays of HRs tuned to f
c = 1000 Hz (a = 0.17 m) in a color scale. Each sub-figure contains the results of IL for 60 numerical simulations of arrays of HRs with different resonant frequencies. The resonant frequency of the HR is plotted on the
y-axis and the frequency of the propagating waves on the
x-axis. For illustration, the lower and upper limits of the Bragg-BG of the SCNB-CC are drawn with vertical white dashed lines at f = 713 Hz and f = 1092 Hz, respectively (criterion IL = 18 dB). Since the phenomenon of local resonance does not occur in an SCNB-CC, the IL is identical for all SCNBs along the entire
y-axis. The white dotted line corresponds to the value of the theoretical resonance frequency of the HR, calculated from (1).
Figure 3 left shows the case of the SCNB with 0° alignment. The IL depends strongly on the resonance frequency of the HRs. The nature of the interaction between Bragg-BG and HR-BG is determined by their relative frequency resonances. When f
HR is far from f
c (f
HR < 600 Hz and f
HR > 1200 Hz), both resonances are well separated. When the resonances are closed (600 Hz < f
HR < 1200 Hz), the interaction is significant and important changes in the IL are observed with respect SCNB-CC. The IL has an asymmetric shape because the transmission is much higher at frequencies f < f
HR than at frequencies f > f
HR. For example, at f
HR = 1400 Hz, the Bragg-BG is similar to that of the SCNB-CC at the same frequency, which agrees quite well with the white vertical dashed lines. The interaction between the resonances of the IL of the SCNB with 90° alignment is shown in
Figure 3 right. It is also asymmetric with respect to the local resonant frequency, but the insulation increase occurs at frequencies below the Bragg frequency. The increase of IL due to the interaction between the Bragg resonance and the local resonance thus strongly depends on the specific orientation of the HRs of the SCNB with respect to the incident plane wave.
Figure 4 shows the IL as a function of frequency for specific cases with f
HR = 400 Hz (f
HR < f
c) and f
HR = 1400 Hz (f
HR > f
c) in red and blue respectively.
Figure 4 left shows the 0° alignment. The SCNB-CC is shown in dashed black. The colored area visually shows the change of the IL of the SCNB with HRs with respect to the SCNB-CC. It can be observed that in the case where f
HR is above the Bragg-BG, the IL increases compared to the SCNB-CC. If the resonance frequency is below, the opposite is true: the insulation of the SCNB decreases. Instead, the IL decreases compared to the SCNB-CC when the f
HR is above the Bragg-BG, as can be seen in
Figure 4 right for a 90° alignment. The IL of the SCNB, on the other hand, increases when the resonance is below. The above results prove that the orientation of the HR mouth in a SCNB is directly related to its insulating properties. Moreover, the resonant frequency of the resonator can be tuned to increase the IL of the SCNB at the Bragg-BG by choosing the alignment of the HRs accordingly. It should be noted that the coexistence of Bragg resonances and HR resonances not only leads to an increase or decrease of IL at the Bragg-BG, but also to a frequency shift of f
c. This is due to the phase shift produced by the HR, which can be explained by a change in the lattice constant experienced by that the wave during propagation [
35].
4. SCNB Scatterer with Two HRs
As a result of the interaction between the Bragg resonance and the local resonances, the increase in insulation of the SCNB depends strongly on the alignment of the mouth of the HR. The relative value of the local resonance (f
HR) with respect to the central frequency of the Bragg-BG (f
c) is an important factor as it determines the interaction between the two phenomena. The use of resonators in SCNBs does not necessarily increase their sound insulation. Therefore, careful attention should be paid to the alignment of the mouths and the resonant frequencies of the HRs in the scatterers of SCNBs.
Figure 5 shows a scatterer for a SCNB consisting of two HR aligned at 0° and 90°, with variable and coupled surfaces. The inner moving wall can be rotated by an angle θ, and simultaneously changes the resonances of the two HRs. Both HR have identical mouths w = 3 cm and L = 1.5 cm. As in
Section 3, IL is computed numerically in a SCNB consisting of scatterers arranged in a two-dimensional square lattice with a filling fraction ff = 0.4 and the central frequency of the Bragg-BG is f
c = 1000 Hz.
Figure 6 shows the behavior of a SCNB with cylinders with two HRs and the influence on the IL due to the movement of the inner moving wall.
Figure 6 left shows the IL in color scale for the SCNB with scatterers with coupled HRs as a function of the rotation angle θ of the wall analyzed with the same numerical model as in
Figure 1b,c. The change in SCNB insulation behavior results from θ-driven surfaces changes in the HRs. Thus, if the HR with 0° alignment has a smaller surface (θ = 0°), the HR with 90° alignment has the largest surface and the IL in the BG-Bragg is high. However, the surfaces ratio of the HRs is changed, IL varies. When the surfaces of the HRs are inverted (θ = 180°), the IL in the BG-Bragg decreases drastically and only the resonances of the respective HRs are observed. HRs do not behave the same with identical surfaces due to the uneven alignment of their mouths in relation to the incident wave. This asymmetry favors the increase of the insulation of the structure when frequencies of the HRs are properly tuned for each alignment.
Figure 6 right shows the IL as a function of frequency for the cases θ = 0°, θ = 180° and the SCNB-CC. As seen in
Section 3, the IL of SCNB at the Bragg-BD increases in 0° alignment for f
HR < f
c and decreases in 90° alignment for f
HR > f
c. This condition is fulfilled for both coupled-surface HRs when θ = 0°. In this case, the surface of the HR with 0° alignment is actually minimum and the surface of the 90° alignment is maximum. It can be seen that the frequency range corresponding to the Bragg-BG (between 500 Hz and 1300 Hz), the IL increases by about 20 dB compared to the SCNB-CC due to the interaction between Bragg and local resonances. The central frequencies of HR-BG of both HRs appear as local maxima in the IL at f = 610 Hz and f = 1165 Hz for the HR aligned at 90°. In this configuration, the resonators have resonance frequencies at the edges of the Bragg-BG. The interaction between the local resonances and the Bragg resonance leads to a reduction in the transmission of the SCNB or, in other words, an increase in the insertion loss. When the inner wall is rotated to θ = 180°, the surfaces of the HRs are reversed and the acoustic behavior of the SCNB changes drastically. The IL of the Bragg-BG changes by more than 30 dB when the wall is rotated from θ = 0° to θ = 180° and the resonance frequencies now shift to f = 534 Hz and f = 1520 Hz. Although both configurations are geometrically equivalent because their HRs have identical surfaces, the central frequencies of the HR resonances are different because of the different alignment of the cylinders with respect to the incident wave.
As
Figure 6 right shows, the IL reaches SPL = 40 dB in the Bragg-BG frequency range. In [
27], the combination of concentric closed-cylinders and cylinders with HR gives numerical and experimental results with less than 20 dB in the Bragg-BG frequency range. Several settings of HRs scatterers have been tested numerically and experimentally in [
23], including porous materials inside the cylinder, but none of the combinations brought more than 25 dB in the whole Bragg-BG frequency range. Thus, our proposal of a barrier with resonators tuned at the edges of Bragg-BG offers a very significant increase in IL compared to other proposals in the literature.
Since SC is not an isotropic material, the acoustic behavior of SCNB with other incidences may be considered. This approach is made considering that the noise source is far away enough, so when source is closer, it would be interesting analyze how the SCNB performs. Instead of rotating the separator of the two HR cavities, rotate the whole cylinder can be also a good choice. By doing this, each scatter will have two different degrees of freedom. Implementing automatic control according to different stimulus could also be an extremely complex and versatile system, with the adaptative ability to perform several cases. Another interesting prospect to study would be now explore periodical configurations and see how the resonant scatterer performs.