Abstract

Ventilation noise control devices often involve a trade-off between their size and ventilating performance, which limits the ability to reduce low-frequency sound in high-ventilation conditions. To address this challenge, the present study explores the use of Hilbert fractal-based design in ventilated metamaterials for improved acoustic performance. The sound transmission loss (STL) of these metamaterials is compared to that of a simple expansion chamber, which serves as the base case. Various parameters, including Hilbert order (O), channel width (K), ventilated space (l), unit cell thickness (H), and the number of unit cells (N) are investigated. Initially, the transfer matrix method evaluates STL without considering thermoviscous effects, which are later incorporated in numerical simulations and impedance tube experiments. The parametric study reveals that increasing the Hilbert curve order decreases the fundamental frequency, while a higher K value increases it. Additionally, more unit cells enhance STL but reduce its broadband nature. Through the finite element method, band diagrams and eigenmodes of Hilbert and base configurations indicate that increased Hilbert orders result in more bands and correspondence between transmission loss spectra and band gaps. The study also identifies dipole resonance modes in the Hilbert structure, which induce a negative effective bulk modulus that contributes to STL. Real-time performance testing in a twin reverberation chamber demonstrates that the Hilbert structure achieves a 5-dB improvement in STL compared to the base configuration across the 700- to 1400-Hz range. These findings are essential for achieving broadband low-frequency noise reduction while allowing airflow.

1 Introduction

The growth of the human population and the continuous expansion of urban areas have contributed significantly to the increasing levels of environmental noise worldwide [1]. This has led to the emergence of noise pollution as a critical urban sustainability issue that substantially impacts health and quality of life [2,3]. Ventilated windows are vital in maintaining optimal thermal comfort, air quality, and natural lighting in domestic buildings. Besides, they are essential to facilitate heat exchange in heavy machinery in industrial environments. However, these ventilated windows are the potential gateways for the noise to enter, which is a potential drawback. Thus, a well-ventilated system with higher noise control is always challenging as there is a trade-off between ventilation and noise cancellation. Therefore, the development of these well-ventilated systems with proper noise control is utmost needed.

In urban buildings, various traditional methods have been suggested in the literature to address this challenge. These methods include different types of facades, protrusions, plenum windows, active noise control windows, balconies, arrays of quarter-wavelength resonators, and Helmholtz resonators [410]. Additionally, acoustic enclosures, barriers, and curtains are utilized to mitigate machinery noise while allowing ventilation [11]. These techniques are effective for noise control, but they can be bulky for low-frequency noise control, as dictated by the mass frequency law. Furthermore, their acoustic performance is significantly influenced by the incident angle.

Acoustic metamaterial has significantly influenced low-frequency noise control applications in recent years. Acoustic metamaterial is an artificial structure exhibiting unique properties, such as negative effective bulk modulus and negative effective density. Liu et al. [12] experimentally proposed the first metamaterial, which unit cell consists of a lead ball covered with silicone rubber. The arrangement of these unit cells in a 3D lattice structure provides the spectral gaps on two orders smaller than the corresponding wavelength of the structure. Further continuation of research leads to broadening the application of acoustic metamaterials such as low-frequency sound absorption [13,14], noise cancellation [15], acoustic superscatterers [16,17], acoustic holography [18], acoustic black holes [19], and acoustic cloaking [20,21].

In recent years, the sonic crystal-based plenum window performance has been investigated by Lee et al. [22]. They found it is capable of traffic noise attenuation of 4.2 dBA at 1000 Hz. Similar to the sonic crystal arrangement, Shen et al. [23] experimentally demonstrated the flow noise control using a circular arrangement of a shunted Helmholtz resonator. This study helps to achieve omnidirectional acoustic shielding without disturbing the airflow. Melnikov et al. [24] proposed the arrangement of C-shaped meta-atoms for machinery noise control. The C-shaped local resonance structure is equivalent to the Helmholtz resonator, allowing the airflow between the meta-atoms and attenuating the noise simultaneously.

The literature discussed above regarding local resonance-based metamaterials plays an essential role in noise attenuation. This type of metamaterials is primarily based on the Helmholtz or quarter-wavelength resonator. The usual bottle-type Helmholtz resonator does not allow airflow. To address this, Kim and Lee [25] proposed the diffraction-type resonator. This design features an airflow passage, a diameter lesser than the wavelength of the acoustic wave to attain strong diffraction. Moreover, they followed another noise control condition by achieving a negative effective bulk modulus. Wang et al. [10] effectively demonstrated the combined effect of utilizing the quarter-wavelength resonator with absorption material in ventilated windows. Due to this combined effect, they achieved a transmission loss of 10–22 dB in the 500–4000 Hz. Similarly, Yu et al. [26] conducted a theoretical and numerical study of the ventilated duct with a periodic array of subwavelength resonators. The results show a 30 dB noise reduction from 600 to 1600 Hz. Kumar et al. [27] achieved low to mid-frequency noise attenuation by applying Helmholtz resonator behavior to the ventilated metamaterial. Wang et al. [28] investigated the perforated and constrained metamaterial for low-frequency (<500 Hz) noise control. Recently, Xiao et al. [29] studied the effect of ventilated openings on sound transmission loss and found that changing the inclination causes broadband transmission loss.

In the sonic crystal arrangement of the local resonators, a major constraint is their large size requirement for low-frequency noise control. Space-coiling acoustic metamaterials have evolved to overcome this disadvantage of the local resonance phenomenon in sonic crystal arrangement. This structure helps slow down the propagation speed within the structure itself. Due to this advantage, research on space-coiling acoustic metamaterials has grown significantly. Additionally, self-similar fractal structures are proposed over arbitrary space-coiling structures for better sound attenuation [30,31]. Hilbert curve is one of the best examples of a self-similar fractal structure. The use of self-similar fractal structures and Hilbert curves have been applied in various areas of sound manipulation, such as wave filtering [32], bandgap formation [3335], acoustic lens [36,37], sound cloaking, tunneling [38,39] and sound absorption [40]. To improve ventilation with noise attenuation, the study of space-coiling ventilated acoustic metamaterials has been undertaken by providing additional ventilation space. However, there is limited literature on this topic, and existing studies do not compare with local resonance or diffraction-based ventilated metamaterial [4143]. Additionally, a detailed summary of recent acoustic ventilated metamaterials is shown in Table 1. Zhao et al. [44], Man et al. [38], and Comandini et al. [40] focused on different applications of Hilbert fractal structures such as transformer noise control, noise control, quarter bending, sound cloaking, sound tunneling, and sound absorption. Still, there is a lack of literature on the theoretical and real-time experimental application of Hilbert fractal structures in ventilated metamaterials, which this study needs to address.

Table 1

Summary of recent trend in acoustic ventilated metamaterial

LiteratureUnit cell size [mm]Structure configurationTarget FrequencyUnit cell arrangementAnalyticalSimulationexperimentVentilation provisionPerformance
Shen et al. [23]R150, H45Helmholtz Resonator2500 HzCircularNot studied10 dB (STL)
Kumar et al. [27]L60, W60Helmholtz Resonator700–1600 HzNot specifiedNot studiedSound absorption and STL studyImpedance tube study18 dB (STL)
Wang et al. [28]R3, H1.575Perforated and constrain metamaterial430 HzNot specifiedTMMSTL studyImpedance tube study20 dB
Melnikov et al. [24]R12, H310Helmholtz Resonator1500 Hz3D MetacagaeNot studiedSTL study13 dB (STL)
Kumar et al. [43]R50, H15Labyrinthine structure530–1230 HzCircularNot studiedSound absorption studyImpedance tube study>0.5 (SAC)
Kumar et al. [49]R12.5, H154Helmholtz Resonator2500 HzCircular metacageNot studiedSIL study20 dB (SIL)
Krasikova et al. [50]R53, H60Helmholtz Resonator1500–16,000 HzCircularNot studiedBand diagram and STL study19 dB (STL)
Xiang et al. [42]L80,H50Space-coiling channel400–1000 HzNot specifiedTMMSRC studyImpedance tube study0.9 (SRC)
Xiao et al. [29]L200,H20Diffraction-type resonator400–2000 Hz1 D periodic duct arrangementNot studiedSTL studyNot studied30 dB (STL)
Ma et al. [41]L100,H10Space-coiling channel577–860 Hz & 1730–1873 Hz1 D & 2 D periodic arrangementFor band structureBand diagram and STL study10 dB (STL)
LiteratureUnit cell size [mm]Structure configurationTarget FrequencyUnit cell arrangementAnalyticalSimulationexperimentVentilation provisionPerformance
Shen et al. [23]R150, H45Helmholtz Resonator2500 HzCircularNot studied10 dB (STL)
Kumar et al. [27]L60, W60Helmholtz Resonator700–1600 HzNot specifiedNot studiedSound absorption and STL studyImpedance tube study18 dB (STL)
Wang et al. [28]R3, H1.575Perforated and constrain metamaterial430 HzNot specifiedTMMSTL studyImpedance tube study20 dB
Melnikov et al. [24]R12, H310Helmholtz Resonator1500 Hz3D MetacagaeNot studiedSTL study13 dB (STL)
Kumar et al. [43]R50, H15Labyrinthine structure530–1230 HzCircularNot studiedSound absorption studyImpedance tube study>0.5 (SAC)
Kumar et al. [49]R12.5, H154Helmholtz Resonator2500 HzCircular metacageNot studiedSIL study20 dB (SIL)
Krasikova et al. [50]R53, H60Helmholtz Resonator1500–16,000 HzCircularNot studiedBand diagram and STL study19 dB (STL)
Xiang et al. [42]L80,H50Space-coiling channel400–1000 HzNot specifiedTMMSRC studyImpedance tube study0.9 (SRC)
Xiao et al. [29]L200,H20Diffraction-type resonator400–2000 Hz1 D periodic duct arrangementNot studiedSTL studyNot studied30 dB (STL)
Ma et al. [41]L100,H10Space-coiling channel577–860 Hz & 1730–1873 Hz1 D & 2 D periodic arrangementFor band structureBand diagram and STL study10 dB (STL)

Note: R, H, L, and W imply the radius, height, length, and width.

STL, sound transmission loss; SIL, sound insertion loss; SAC, sound absorption coefficient; SRC, sound reduction coefficient; TMM, transfer matrix method.

However, these studies primarily focused on 2D lattice configurations and did not address low-frequency broadband noise reduction in ventilated systems. In contrast, this research highlights the use of a 1D periodic duct arrangement for reducing noise in ventilated systems, particularly focusing on the low-frequency range, specifically below 1600 Hz. This frequency range is crucial for controlling outdoor urban noise in buildings, especially traffic noise. Studies [4548] have demonstrated that traffic noise, a significant contributor to urban sound pollution, has peak sound pressure levels and high spectral energy content within this range.

The present study addresses these gaps by investigating the sound attenuation performance of an innovative Hilbert curve–based acoustic ventilated metamaterials. The evaluation encompasses a comprehensive approach using the transfer matrix method, rigorous computational analysis, and experimental validation with an impedance tube. The impact of various geometric parameters on sound transmission loss is explored by analyzing first-, second-, and third-order Hilbert curves, as well as a base model. Furthermore, the finite element method is employed to construct band diagrams and eigenmodes, elucidating the fundamental mechanisms driving sound transmission loss. Effective properties are examined through phase plots to enhance understanding. The study is grounded in controlled impedance tube experiments and realistic assessments using twin reverberation chamber tests. This multifaceted approach advances the theoretical understanding and paves the way for practical applications of ventilated metamaterials in sound attenuation.

2 Material Design and Methodology

This section provides a detailed approach for developing and analyzing Hilbert curve ventilation metamaterials. It starts by explaining the design principles and construction process of the metamaterial structure. Then, it describes the methods used to assess its acoustic performance, including the transfer matrix method, experimental setup and procedures, and computational studies. Finally, the section concludes with an analysis of the effective properties, focusing on its bulk modulus and density to offer insight into its acoustic behavior.

2.1 Design and Construction of Hilbert Curve–Based Acoustic Ventilated Metamaterials.

The base case geometry features a simple expansion chamber arrangement, which aligns with prior literature [25,26,29,51], highlighting the importance of a simple expansion chamber design in acoustic applications. This study introduces the Hilbert structure in the base case cavity space. The Hilbert curve ventilated metamaterials are designed by dividing the base case cavity into four sections. Further, each section is divided into Here, N = 1, 2, and 3 for first, second, and third orders. The center of 2N×2N section is connected from lower left to lower right clockwise with the Hilbert channel width K, as shown in Fig. 1(a). Additionally, the sound wave propagation length in the 1/4th of the Hilbert structure space is Lp=L(2N11/2N+1), N = 1, 2, and 3. Here, L is the length of the unit cell.

Fig. 1
Schematic details of Hilbert curve–based ventilated metamaterials and base case. (a) Cross-sectional view of air domain in Hilbert curve and base case, (b) unit cell of Hilbert curve and base case with nomenclature, and (c) the base duct and Hilbert structure space are represented as equivalent quarter-wavelength assumption for transfer matrix method.
Fig. 1
Schematic details of Hilbert curve–based ventilated metamaterials and base case. (a) Cross-sectional view of air domain in Hilbert curve and base case, (b) unit cell of Hilbert curve and base case with nomenclature, and (c) the base duct and Hilbert structure space are represented as equivalent quarter-wavelength assumption for transfer matrix method.
Close modal

Figure 1(b) illustrates the 3D unit cell of the first-order (O1), second-order (O2), and third-order (O3) Hilbert curve ventilated metamaterial and base case. In the Hilbert configurations, the blue area depicted in Figs. 1(a) and 1(b) represents the air domain through which acoustic waves propagate, with the Hilbert channel width K ranging from 3 to 6 mm. The overall length of this structure is denoted as L and measures 70 mm. Additionally, H is the total thickness of the unit cell varying by 24, 26, and 28 mm, and l×l defines the spatial dimension for airflow entry varying between 15 × 15, 20 × 20, and 25 × 25 mm2. For the O3 configuration, the K value is fixed at 3 mm due to geometry constraints and no involvement of K in the base case. Therefore, the total combinations of O1, O2, O3, and base case are 270 (108 + 108 + 27 + 27). To simplify, Table 2 shows a comprehensive overview of the various parameters utilized in this study, based on which the acoustic characteristics of the Hilbert structure are studied. Throughout the study, each sample is named by using the parameter values; for example, the sample O2_K3_l20_H24_N2 represents the second-order Hilbert curve, 3-mm Hilbert channel width, 20-mm ventilated space opening, 24 mm is the thickness of the unit cell, and 2 number of the unit cells, respectively.

Table 2

Sample parameters

Order of the Hilbert curve (O)Hilbert channel width (K) in mmVentilated space area (l × l) in mm2Thickness of the unit cell (H) in mmNumber of unit cell (N)Geometry name
First order315 × 15241O1_K3_l15_H24_N1
320 × 20241O1_K3_l20_H24_N1
325 × 25241O1_K3_l25_H24_N1
315 × 15242O1_K3_l15_H24_N2
315 × 15243O1_K3_l15_H24_N3
615 × 15241O1_K6_l15_H24_N1
Second order315 × 15241O2_K3_l15_H24_N1
315 × 15261O2_K3_l15_H26_N1
315 × 15281O2_K3_l15_H28_N1
320 × 20241O2_K3_l20_H24_N1
325 × 25241O2_K3_l25_H24_N1
415 × 15241O2_K4_l15_H24_N1
515 × 15241O2_K5_l15_H24_N1
315 × 15242O2_K3_l15_H24_N2
615 × 15243O2_K6_l15_H24_N3
515 × 15241O2_K5_l15_H24_N1
Third order315 × 15241O3_K3_l15_H24_N1
320 × 20241O3_K3_l20_H24_N1
325 × 25241O3_K3_l25_H24_N1
315 × 15242O3_K3_l15_H24_N2
325 × 25243O3_K3_l25_H24_N3
Base15 × 15241Base_ l15_H24_N1
20 × 20241Base_l20_H24_N1
25 × 25241Base_l25_H24_N1
15 × 15242Base_l15_H24_N2
15 × 15243Base_l15_H24_N3
Order of the Hilbert curve (O)Hilbert channel width (K) in mmVentilated space area (l × l) in mm2Thickness of the unit cell (H) in mmNumber of unit cell (N)Geometry name
First order315 × 15241O1_K3_l15_H24_N1
320 × 20241O1_K3_l20_H24_N1
325 × 25241O1_K3_l25_H24_N1
315 × 15242O1_K3_l15_H24_N2
315 × 15243O1_K3_l15_H24_N3
615 × 15241O1_K6_l15_H24_N1
Second order315 × 15241O2_K3_l15_H24_N1
315 × 15261O2_K3_l15_H26_N1
315 × 15281O2_K3_l15_H28_N1
320 × 20241O2_K3_l20_H24_N1
325 × 25241O2_K3_l25_H24_N1
415 × 15241O2_K4_l15_H24_N1
515 × 15241O2_K5_l15_H24_N1
315 × 15242O2_K3_l15_H24_N2
615 × 15243O2_K6_l15_H24_N3
515 × 15241O2_K5_l15_H24_N1
Third order315 × 15241O3_K3_l15_H24_N1
320 × 20241O3_K3_l20_H24_N1
325 × 25241O3_K3_l25_H24_N1
315 × 15242O3_K3_l15_H24_N2
325 × 25243O3_K3_l25_H24_N3
Base15 × 15241Base_ l15_H24_N1
20 × 20241Base_l20_H24_N1
25 × 25241Base_l25_H24_N1
15 × 15242Base_l15_H24_N2
15 × 15243Base_l15_H24_N3

2.2 Theory of Transfer Matrix Approach for Base and Hilbert Curve Ventilated Metamaterials.

The transmission loss characteristics of the base geometry and designed Hilbert metamaterials are obtained through the transfer matrix method using Eq. (1) [51,52].
(1)
The acoustic pressure (P) and particle velocity (V) on both ends of the unit cell are correlated by using the transfer matrix of the base case (Tβunit) and Hilbert metamaterial (THunit) as shown in Eqs. (2) and (3) [51].
(2)
(3)
Here, Tf and Tb are the transfer matrix of the front and back cover.
(4)
where k0=2πf/C0 is the acoustic wave number in the air is medium, C0 is the speed of sound in air, and A0 is the ratio of the opening area to the total area of unit cell (A0=l2/L2). The effective length of the airflow opening is Lc is defined as [5254]:
(5)
where th refers to the thickness of the Hilbert structure.
To determine the transfer matrix of the base case duct (Tβ) and Hilbert structure (TH), an assumption of the equivalent quarter-wavelength resonator is considered, as shown in Fig. 1(c). These matrices are given below [51]:
(6)
where Ae_b=(L2l2)tb/Le is the area of base duct space, Ae_h=Lp×K is the area of the single effective medium, Le=(L2l2)/2 is the effective length (same for the base, and Hilbert case), Lp=L(2N11/2N+1) is wave propagation length, N = 1, 2, and 3 is the order of Hilbert curve, and K is the thickness of the Hilbert channel.
The acoustic impedance for the base case duct [51] (Zβ) and Hilbert structured space [55] (ZH) is:
(7)
where L is the length of the unit cell, l is the airflow opening length, α is the normalized effective length [51], tb is the thickness of the base case duct space, and th is the thickness of the Hilbert structured space. The Ze represents the effective impedance of the single effective medium equation as follows [54]:
(8)
where ke=2πf/Ce is the wavenumber, Le length of the effective medium, and ζ is the ratio of the effective medium cross-section are (Ae_h) a to the opening cross-section area (l2) of the unit cell [54]. The effective properties depend on the geometry parameter and length of the sound wave travelling path in the Hilbert structure. According to Zhao et al. [56], the effective speed of sound (Ce=(Le×C0)/Lp) is derived from the conservation of sound wave traveling time in the Hilbert channel, and effective density (ρe=(C0×ρ0)/Ce) can be derived from the impedance mismatch condition. Therefore, the transmission loss is calculated using the components of the unit cell transfer matrix (Tβ_unit&TH_unit) as follows:
(9)

The proposed transfer matrix model helps to analyze the acoustic characteristics of the Hilbert fractal structure without considering the thermoviscous effect. This approach will be used in Sec. 3.1 to predict the acoustic performance of the base and Hilbert curve ventilated metamaterials. For a more robust analysis, conducting impedance tube experiments and including the thermoviscous effect in the computational model would be beneficial. These will be studied in the following sections.

2.3 Experimental Setup and Procedure.

This section provides a detailed explanation of the methodologies and equipment utilized in the acoustic experiments. It begins with an overview of the sample preparation process, including the materials used and their specifications. Subsequently, it describes the impedance tube setup and the method for measuring sound transmission loss. Finally, the section concludes with an outline of the reverberation test procedure, covering the test chamber setup, measurement equipment, and the process for conducting the experiments.

2.3.1 Impedance Tube Test.

To realize the experimental setup, the sample is fabricated using Fused Deposition Modeling (FDM) 3D printing technology. The material chosen for this construction is poly-lactic acid (PLA), with a density of 1250kg/m3, a Poisson's ratio of 0.35, and Young's modulus of 3.5 GPa. The 3D-printed samples, as displayed in Figs. 2(a) and 2(b), were manufactured with a diameter of 100 mm and included a face cover. A unit cell size of 70 mm is embedded in the 100-mm diameter structure to support the sample, facilitating attachment to an impedance tube. The face cover is carefully attached to the Hilbert structure model using cyanoacrylate glue to ensure airtightness and prevent potential leakage. The 3D printing process utilized the Flash Forge Creator 3 Pro 3D printing machine, noted for its printing precision of ±0.2 mm, a layer height of 0.1 mm, and a print speed of 90 mm/s. These specific printing parameters are crucial in achieving the desired metamaterial structure, ensuring precision and accuracy throughout the manufacturing process.

Fig. 2
Sample and experimental setup: (a) photograph of 3D-printed samples and face cover, (b) assembled view of unit cell, and (c) schematic details of impedance tube setup
Fig. 2
Sample and experimental setup: (a) photograph of 3D-printed samples and face cover, (b) assembled view of unit cell, and (c) schematic details of impedance tube setup
Close modal

The acoustics performance of the Hilbert curve–based metamaterials is evaluated through experimental tests conducted in the BSWA SW 422 impedance tube system at the Acoustics and Dynamics Lab in IIT Madras. The test is performed based on the transfer matrix method as per ASTM E1050-08 and ASTM E2611-09 Standards. The experiments are performed in the frequency range of 250–1600 Hz. As shown in Fig. 2(c), this experimental test setup incorporated additional hardware components such as a 50 W power amplifier, four 1/4″ microphones, a four-channel data acquisition card—MC3242, and VA-LAB 4 software—for data processing.

A calibration process is performed before commencing the experiments. The four microphones are calibrated using a pistonphone calibrator to ensure accurate measurements. The sample is meticulously secured with Teflon tape along its external diameter to avoid sound leakage and provide a tight fit with an impedance tube. Each sample is tested five times to ensure repeatability and consistency, and the average of five measurements is used to calculate the final results. The impedance tube test results will be used to validate the computational study that incorporates thermoviscous effects, as presented in Sec. 3.1.

2.3.2 Reverberation Chamber Test.

Figure 3 illustrates the test facility comprising two closed pentagon-shaped reverberation chambers. The difference between the volumes of the chambers is kept at more than 10%, following the ISO 10140 standards [57]. The chambers have a cut-off frequency of 315 Hz. The B&K sound source type 4224 generates white noise in the larger chamber. As shown in Fig. 3, the 3D-printed sample is fixed between the two chambers. This study uses two one-fourth-in. condenser microphones (Microtech Gefell-M370, Gefell, Germany). One microphone is placed in the sound source side, and the other is placed in the receiver chamber to record the acoustic data. The sensitivity of each microphone is 12.1 and 12.2 mV, calibrated using a B&K type 4228 pistonphone calibrator. Both microphones have a flat frequency response between 50 and 20,000 Hz within ±1 dB. The microphone signal is sampled at 12,800 per s using the National Instruments data acquisition board (NI-9234). The reverberation time of the larger and smaller chambers is measured and found to be 1.47 and 1.34 s, respectively. The sound transmission loss (TL) is calculated using the following equation:
(10)
Fig. 3
Schematic of a twin reverberation chamber. The S1, S2, and S3 are the microphone positions on the sound source side, and R1, R2, and R3 are the microphone positions on the sound receiver side, and the test sample consists of a 3×3 arrangement of the unit cells
Fig. 3
Schematic of a twin reverberation chamber. The S1, S2, and S3 are the microphone positions on the sound source side, and R1, R2, and R3 are the microphone positions on the sound receiver side, and the test sample consists of a 3×3 arrangement of the unit cells
Close modal

where L1andL2 is the average sound pressure level in the source and receiver room, S is the area of the test sample partition (m2), A=0.612V/T, V is the volume of the sound-receiving reverberation chamber (m3), and T is the reverberation time (s) of the receiver chamber.

The reverberation chamber test method described here will be used to evaluate the acoustic performance of Hilbert-based structures in a more realistic, multi-directional sound field. This approach complements the plane wave conditions of the impedance tube experiments by simulating reflection-rich environments. The results from these tests, presented in a subsequent section, will provide insights into the materials' behavior under diffuse sound field conditions, offering a clearer understanding of their acoustic characteristics in practical applications.

2.4 Computational Study.

The acoustic behavior of the samples is analyzed using COMSOL Multiphysics 6.1 [58] software for eigenfrequency and sound transmission loss studies. The eigenfrequency study helps to obtain the band diagram. Ideally, wherever the bandgap forms there, the sound wave will not propagate. Furthermore, to understand the behavior of sound transmission loss (STL), the studies are conducted using pressure acoustics and thermoviscous acoustic modules. The pressure acoustics study does not account for losses, while the thermoviscous acoustic module includes viscous and thermal losses. The solid domain has a higher acoustic impedance than the air medium, so sound wave propagation is restricted within the air domain. Therefore, the rigid wall approximation is used for the fluid-solid interface to simplify the computational model. The sound waves are propagated in the x-direction of the domain, as shown in Fig. 4. The Perfectly Matched Layer (PML) is employed at both ends of the domain, which helps avoid reflection to the propagation field. In order to accurately represent the physics, the thermoviscous model and boundary layer mesh are utilized in this study. The maximum element size is fixed by L/10, where L is the length of the sample. The boundary layer mesh is composed of five layers, and the thickness of each layer is defined as dvisc(mm)=0.22(100/f0). Here, f0 is 1600Hz. The transmission loss is calculated as
(11)
where τ denotes the power transmission coefficient. Power transmission coefficient is the ratio of transmitted to incident acoustic power. The τ can also be expressed in terms of acoustic pressures as τ=|pt|2|pi|2, where pt is the transmitted acoustic pressure measured at the outlet of the sample and pi is the incident acoustic pressure measured at the inlet of the sample. The proposed computational method with a thermoviscous model for sound transmission loss described in this section will be used to analyze the acoustic performance of the designed structure through a parametric study, as described in Secs. 3.2, 3.3 and 3.4. Additionally, the method for calculating band diagrams outlined here will be employed to explain the dispersive nature of the structure, with results presented in Sec. 3.5.
Fig. 4
Schematic details of the computational domain for sound transmission loss study
Fig. 4
Schematic details of the computational domain for sound transmission loss study
Close modal

2.5 Effective Property Analysis.

Determining the acoustic properties of a heterogeneous medium is a challenging task. To address this issue, the effective medium theory can approximate a medium as homogeneous and isotropic [59]. Furthermore, negative effective parameters play a crucial role in acoustic metamaterial and can help to explain negative refraction, acoustic cloaking, and acoustic lensing. The retrieval method is commonly used for complex metamaterial structures to determine the effective parameters using complex reflection and transmission coefficients. This method involves finding effective parameters using complex transmission and reflection coefficients. The effective impedance ξ and refractive index (n) are obtained using Eqs. (12) and (13) [51,59].

(12)
(13)
where k is the wave number and H is the thickness of the structure.
(14)
(15)
The effective bulk modulus and density can be expressed in Eqs. (16) and (17).
(16)
(17)

The effective bulk modulus and density are calculated using computational results of reflection and transmission coefficients. The results from this method are presented in Sec. 3.6, demonstrating how these effective properties contribute to the overall acoustic performance of the structure.

3 Results and Discussions

This section discusses the performance of the Hilbert curve acoustic metamaterials and a base case using analytical, experimental, and computational methods. It also explores the impact of thermoviscous properties. The study delves into parametric studies to optimize the structure for broadband transmission loss. Additionally, it discusses the results of band diagrams and effective properties to understand the reasons for sound transmission loss.

3.1 Effect of Thermoviscous Property on Acoustic Ventilated Structures.

Figure 5 shows the transmission loss results from the transfer matrix method (theory) and simulation (lossless) results. Three different Hilbert orders and base cases are analyzed. The geometry features are a 3 mm Hilbert channel width (K3), 15 × 15 mm2 ventilated space (l15), a structure thickness of 24 mm (H24), and a single unit cell (N1). Among the results, the Hilbert-based structure performs better than the base case due to the longer effective path length leading to the higher interaction time between the acoustic waves and the structure, which enhances the sound attenuation performance in the Hilbert structure. Notably, the second-order configuration provides broadband sound transmission loss between 700 Hz and 1600 Hz. Additionally, the results show that the resonance frequencies obtained by the transfer matrix method are in acceptable agreement with the lossless computational result.

Fig. 5
Comparison of sound transmission loss performance between theory and simulation for ventilated acoustic metamaterials and base case without considering thermo viscous loss
Fig. 5
Comparison of sound transmission loss performance between theory and simulation for ventilated acoustic metamaterials and base case without considering thermo viscous loss
Close modal

However, in real situations, thermoviscous effects play a crucial role in propagating acoustic waves in complex and narrow geometries. This effect involves the dissipation of sound energy due to thermal conduction and the viscous nature of the medium through which the acoustic waves propagate. In Fig. 6, the transmission loss of the COMSOL simulation result is compared with the impedance tube experiment. Six different samples are used for validation. The computational model is well adapted to the experimental data and agrees well with the resonance frequency of the structures.

Fig. 6
Thermo viscous simulations implemented to validate the impedance tube sound transmission loss experiment for Hilbert space ventilated material and the base case
Fig. 6
Thermo viscous simulations implemented to validate the impedance tube sound transmission loss experiment for Hilbert space ventilated material and the base case
Close modal

All Hilbert configurations followed the distinct asymmetric Fano-like resonance [60,61], except for the base case. Fano-like resonance is an asymmetric profile formed by interference between the sound wave passing through the resonating and non-resonating channels. The asymmetric profile has the peak portion formed due to destructive interference. In the first-order geometry (O1), there is a maximum transmission loss of 40 and 50 dB for 3 (K3) and 6 (K6) mm Hilbert channel width, respectively. This comparison shows that the maximum Hilbert channel width (K) results in the maximum transmission loss, the physical insights of which will be discussed in Sec. 3.3. Among the tests, the O2_K3_l15_H24_N1 and O2_K6_l15_H24_N1 configurations performed in the broad range of transmission loss, which is almost one-octave band, i.e., 695–1600 Hz and 715–1600 Hz, respectively. The computational analysis of the third-order (O3) configuration indicates the presence of two distinct Fano-like resonances, one at around 400 Hz and the other at around 1100 Hz. However, the experimental observations slightly mismatch the computational study. The reason for this is the lower wall thickness of the O3 geometry, which creates unwanted sound leakage. Overall, the computational study agrees well with the experiments. Additionally, widening the (K) from 3 to 6 mm also increases peak transmission loss observed for O1 and O2 configurations.

While comparing the sound transmission loss with and without consideration of thermoviscous loss (Figs. 5, 6, and Fig. 15 in the  Appendix), it is found that considering thermoviscous loss gives a more realistic result. The lossless model predicts resonant behavior but overestimates sound transmission loss by neglecting the viscous damping effect. For example, the peak transmission loss for model O1_k3_l15_H24_N1 considering thermoviscous effects is 38.6 dB, whereas the lossless case is 53.5 dB. For model O2_k3_l15_H24_N1, the transmission loss is 28 dB when considering losses, compared to 43.2 dB in the lossless case. Similarly, for model O3_k3_l15_H24_N1, the lossless case overestimates the peak transmission loss by 15 dB at both the resonance frequencies. On the other hand, considering the thermoviscous impact in the simulation aligns well with the experimental results. However, a curved narrow channel is not involved in the base case, so the thermoviscous effect is not dominant. These results emphasize the necessity of a thermoviscous model to design the fractal narrow structure metamaterial. The effect of various parameters as mentioned in Table 2 is discussed in the following sections.

3.2 Effect of Ventilated Space Opening (l).

Figure 7 shows the transmission loss performance of O1, O2, O3, and base configuration for three differently sized ventilated openings (l×l=15×15,20×20and25×25mm2). In this analysis, other parameters are fixed, including a Hilbert channel width of 3 mm, a thickness of 24 mm, and a single unit cell. The results are derived from COMSOL transmission loss simulations using the thermoviscous model. Among the models, the O2_K3_l15_H24_N1 model provides the best average transmission loss of 22 dB in the frequency range of 695–1600 Hz. However, for O2 configuration, when the ventilated space is increased to 25×25mm2, the average transmission loss reduces to 15 dB, and the resonant frequency shifts to 800 Hz. This trend is observed in other configurations as well. It is clear that increasing the ventilated space results in a resonance frequency shift to the higher range. Besides, the maximum transmission loss decreases due to the reduction in the reflection of the acoustic wave. Despite using the Hilbert channel, there is still a trade-off between airflow and control of sound wave transmission. However, based on observations, the Hilbert-based configuration is better equipped to control sound wave propagation over a broader frequency range than the base case.

Fig. 7
Transmission loss simulation (thermoviscous) results showing the dependence on variation of ventilation space (l) for Hilbert space ventilated metamaterials
Fig. 7
Transmission loss simulation (thermoviscous) results showing the dependence on variation of ventilation space (l) for Hilbert space ventilated metamaterials
Close modal

3.3 Effect of Hilbert Channel Width (K) and Thickness of Unit Cell (H).

Figure 8(a) illustrates the effect of varying Hilbert channel width (K) for O2 configuration. The K value varies from 3 to 6 mm with fixed parameters of l15, H24, and N1. The results are obtained through COMSOL transmission loss simulation with the thermoviscous model. When the K value is increased in this configuration, there is an increase in transmission loss along with a frequency shift. The same trend is observed for the O1 configuration (see Fig. 16 in the  Appendix). To enhance the observed result, the inset plot in Fig. 8(a) illustrates the relationship between the real part of the impedance and the variation in K. It shows that K = 6 mm results in a higher impedance value than K = 3 mm, indicating greater energy dissipation through viscous losses in the fluid, leading to higher transmission loss. Additionally, the resonance frequency plot inset in Fig. 8(a) shows an increase in resonance frequency with respect to K, resulting in the shift of transmission loss to a higher frequency range. In summary, the increase in frequency and transmission loss with increasing K is mainly due to changes in resonance frequencies and impedance characteristics of the structure, respectively.

Fig. 8
Transmission loss simulation (thermo viscous) result showing the dependence on (a) variation of Hilbert space width (K), and (b) variation of the thickness of unit cell (H)
Fig. 8
Transmission loss simulation (thermo viscous) result showing the dependence on (a) variation of Hilbert space width (K), and (b) variation of the thickness of unit cell (H)
Close modal

Figure 8(b) shows the variation of the unit-cell thickness (H) of the O2 configuration with fixed parameters of K3, l15, and N1. It is observed that the resonance frequency does not change significantly with varying H values. However, the transmission loss has an average variation of 2 dB with low to high H value, and a similar trend is observed for O1 and O3 configurations (Fig. 17 in the  Appendix).

3.4 Effect of Number of Unit Cell (N).

This section demonstrates the effect of the number of unit cells in the periodic arrangement (in x-direction as inset in Fig. 9) on the transmission loss capabilities. The analysis considers the samples of O1, O2, O3, and base cases, with varying unit cell numbers ranging from 1 to 3 by using thermoviscous module in COMSOL simulation. Figure 9 reveals no significant difference in the peak frequency when increasing the number of unit cells. However, there is a slight decrease in the broadband range due to an unsmooth transition of impedance mismatch. Additionally, the results show that as the number of unit cells increases, the transmission loss spectrum tends to align more closely with the bandgap in the band diagram (Fig. 18 in the  Appendix). This trend is consistent with the work of Xiao et al. [29], which shows that increasing the number of unit cells leads to alignment with the bandgap in the band diagram.

Fig. 9
Transmission loss simulation (thermo viscous) result showing the dependence on variation of number of unit cells (N)
Fig. 9
Transmission loss simulation (thermo viscous) result showing the dependence on variation of number of unit cells (N)
Close modal

The band diagram is studied in the following section to comprehend the broadband range of sound transmission loss. Moreover, the amplitude of sound transmission loss increases due to enhanced wave reflection within the periodic structure. The results from various geometric configurations suggest that the parameters K and l are more important than H. Therefore, for further analysis, H and l are set at their lower limit of 24 mm and 15 mm, respectively, while K at 3 and 6 mm.

3.5 Eigen Frequency and Band Diagram Analysis.

Figure 10, known as a dispersion curve or band diagram, depicts the relationship between the normalized wavenumber and frequency. This study defines the Brillouin zone [62] π/Htoπ/H from the unit cell. The bandgaps are in band diagram highlights in cyan color. It is observed that increasing the order of the Hilbert structure results in the bandgap shift to a lower frequency. The flat bands observed in O1, O2, and O3 are responsible for the fundamental resonance of their structures. Additionally, increasing the Hilbert channel width (K) shifted the fundamental frequency to a slightly higher range. For instance, in the second-order configuration (O2), increasing K shifted the fundamental frequency from 695 to 732 Hz. Moreover, the transmission loss spectra in Fig. 6 align well with the corresponding bandgap.

Fig. 10
Band diagram for different configurations with acoustic pressure distribution for fundamental resonance
Fig. 10
Band diagram for different configurations with acoustic pressure distribution for fundamental resonance
Close modal

As the tortuosity increases, the fundamental frequency shifts toward the lower frequency, leading to increased bands. In the base case, no flat band is observed due to the absence of a fractal structure. The presence of a flat band is closely connected to the asymmetrical Fano profile in the transmission loss spectra. Fano-like resonance is an asymmetrical profile caused by the interaction of a sound wave moving through resonating and non-resonating channels. The Hilbert structure introduces localized discrete forms of resonant modes at a specific frequency (resonant channel). This resonant mode interacts with the continuous propagating wave mode in the periodic structure (non-resonant channel), leading to Fano resonance. This phenomenon means that the horizontal line represents the resonant frequency, showing that these resonant modes are localized and prevent the acoustic wave from propagating through the structure.

Furthermore, in Fig. 10, the eigenfield or pressure distribution of the fundamental resonance frequency is also included. The bandgap exhibits dipolar characteristics attributed to Fano resonance in all cases except the base case. Besides, the third-order (O3) geometry has an additional bandgap in the higher frequency range due to the increased wave propagation length and belongs to higher order modes (1118 Hz). Among all the configurations, the O2 exhibits a wide bandgap range primarily for higher Hilbert channel width (K = 6 mm), ranging from approximately 700 to 1400 Hz.

Figure 11 illustrates the relationship between bandgap width and the variation of K, H, and l values, as determined from the band diagram using the eigenfrequency module in COMSOL [63]. In Fig. 11(a), a clear trend shows for constant H and l; the band gap width expands with an increase in K, such as in the O2 configuration, the band gap width increases from 525 Hz (695–1220 Hz) at K = 3 mm to 668 Hz (732–1400 Hz) at K = 6 mm. Notably, the base case does not involve K, and there is no variation of K for O3 due to geometry constraints. Figure 11(b) reveals that for the Hilbert configuration, there is no significant variation in band gap width as H increases. However, the base case band gap widens from 296 Hz (1304–1600 Hz) to 434 Hz (1166–1600 Hz) with increasing H. This band gap variation aligns with the transmission loss variation observed in Figs. 8(a) and 8(b). In contrast, Fig. 11(c) demonstrates a decreasing trend in band gap width as l increases, consistent with the results presented in Fig. 7. Overall, these findings suggest that the K value has the most pronounced positive impact on band gap width, while l exhibits an inverse trend. The H value primarily influences the base case configuration, with minimal effect on the Hilbert structure.

Fig. 11
Band gap variation for Hilbert acoustic metamaterials and base case. (a) Variation of Hilbert channel width (K), (b) variation of unit cell thickness (H), and (c) variation of ventilation space (l).
Fig. 11
Band gap variation for Hilbert acoustic metamaterials and base case. (a) Variation of Hilbert channel width (K), (b) variation of unit cell thickness (H), and (c) variation of ventilation space (l).
Close modal
Figure 12 illustrates the phase velocity related to normalized frequency. The phase velocity [64] (cp) is calculated using the following equation:
(18)
where f is the frequency, and k is the wavenumber, which is obtained from a band diagram performed using COMSOL Multiphysics. Normalized frequency is the ratio of the actual frequency (f) to the fundamental frequency (fR) of the structure. As seen in Fig. 12, the phase velocity decreases significantly for different structures compared to the speed of sound in the air. It is worth noting that the phase velocity approaches zero just before the bandgap starts. Therefore, the resonance of the structures occurs in the deep subwavelength region [64], as shown in Table 3.
Fig. 12
Variation of Phase velocity of the band diagram for Hilbert space ventilated material and the base case
Fig. 12
Variation of Phase velocity of the band diagram for Hilbert space ventilated material and the base case
Close modal
Table 3

Subwavelength thickness of the geometry

GeometrySubwavelength thickness
O1_K3_l15_H24_N1λ/10.3
O1_K6_l15_H24_N1λ/10.0
O2_K3_l15_H24_N1λ/20.5
O2_K6_l15_H24_N1λ/19.5
O3_K3_l15_H24_N1λ/38.0
Base_l15_H24_N1λ/11.0
GeometrySubwavelength thickness
O1_K3_l15_H24_N1λ/10.3
O1_K6_l15_H24_N1λ/10.0
O2_K3_l15_H24_N1λ/20.5
O2_K6_l15_H24_N1λ/19.5
O3_K3_l15_H24_N1λ/38.0
Base_l15_H24_N1λ/11.0

3.6 Analysis of Effective Properties.

The results of the effective properties obtained through the retrieval method (Sec. 2.5) are displayed in Fig. 13. It is observed that only the effective bulk modulus has a negative value for all configurations, which is highlighted in cyan color. The negative effective bulk modulus is due to the inward and outward movement of the structure's diagonal (insert in Fig. 13), resulting in the external force reversing the rate of volumetric change. These negative effective properties are responsible for bending the sound wave and leading to high acoustic wave reflection, ultimately causing transmission loss.

Fig. 13
Effective bulk modulus and density for first-, second-, and third-order Hilbert structure and base model. The negative bulk modulus ranges are highlighted in cyan.
Fig. 13
Effective bulk modulus and density for first-, second-, and third-order Hilbert structure and base model. The negative bulk modulus ranges are highlighted in cyan.
Close modal

The effective bulk modulus for O1_K3_l15_H24_N1 and O1_K6_l15_H24_N1 is negative and ranges from 1385 to 1600 Hz and 1411 to 1600 Hz, respectively. In the negative bulk modulus regime, the transmission loss and bandgap are well-aligned due to the decrease in acoustic wave propagation velocity. A similar trend is observed for O2_K3_l15_H24_N1 and O2_K6_l15_H24_N1 configurations with a broader range of transmission loss. For O3, the negative bulk modulus occurs within two specific frequency ranges (372–950 Hz and 1118–1600 Hz), as observed in the bandgap and sound transmission loss study. Moreover, the base structure has the most negligible impact on negative bulk modulus among the different configurations. Evidently, for all Hilbert configurations, the resonance modes are dominated mainly by dipole resonance and, for the base case, by monopole resonance.

3.7 Analysis of Sound Transmission Loss in Twin Reverberation.

In the impedance tube experiment, the O2_K6_l15_H24_N1 configuration exhibited broadband transmission loss effectively. It is worth noting that the impedance tube experiment operates under plane wave conditions only, whereas real-life situations present a diffuse field environment. Consequently, the O2_K6_l15_H24_N1 and base configuration are chosen for the reverberation chamber test to achieve accurate real-time experimental results. Figure 14 compares the transmission loss of O2_K6_l15_H24_N1 with the base configuration. This O2 configuration demonstrates a significant increase in transmission loss across frequencies, ranging from approximately 700–1400 Hz. It consistently outperforms the base case with an average transmission loss of 5 dB. These results confirm the expected formation of a bandgap, which is further supported by impedance tube measurements. This improvement in sound attenuation with the O2 configuration emphasizes the use of Hilbert curve–based ventilated metamaterials to reduce sound propagation, which has implications for various noise control and acoustic engineering applications.

Fig. 14
Experimental comparison of sound transmission loss performance between O2_K6_l15_H24_N1, and Base_l15_H24_N1 using reverberation chamber testing
Fig. 14
Experimental comparison of sound transmission loss performance between O2_K6_l15_H24_N1, and Base_l15_H24_N1 using reverberation chamber testing
Close modal

4 Conclusion

This study evaluated the effectiveness of using the Hilbert curved acoustic metamaterials over the simple expansion chamber for sound attenuation. The transfer matrix method, computational analysis, and experimental results are used to study the transmission loss for Hilbert and base configurations. The study delved into the impact of Hilbert orders and geometry parameters, noting that increasing the Hilbert order reduces the fundamental frequency and affects broadband transmission loss spectra due to resonance frequencies and impedance characteristics of the structure, respectively. Besides, increasing the K value raised the fundamental frequency of the Hilbert configuration. The results showed that the Hilbert-based configuration had a Fano resonance in the transmission loss spectra. Specifically, the Hilbert second-order 3-mm channel width, 15-mm ventilation space, and 24-mm thickness (O2_K3_l15_H24_N1) configuration exhibited a broadband transmission loss of 22 dB within the frequency range of 695–1600 Hz, while the Hilbert second-order 6-mm channel width, 15-mm ventilation space, and 24-mm thickness (O2_K6_l15_H24_N1) configuration achieved a 24-dB transmission loss within the frequency range of 732–1600 Hz.

Furthermore, the study also investigated the band diagram of Hilbert and base configurations to find the relation between broadband transmission loss spectra and bandgaps. The broadband transmission loss spectra matched well with the bandgap in the band diagram. The results showed that increasing the number of orders led to an increase in the number of bands and the O2_K3_l15_H24_N1 and O2_K6_l15_H24_N1 configurations had wider bandgaps in the subwavelength thickness of λ/20. Analysis of effective parameters shows that the negative bulk modulus is due to dipole resonance in the Hilbert configuration and monopole for the base.

Real-time testing in the twin reverberation chamber demonstrated that the Hilbert structure, particularly the O2_K6_l15_H24_N1 configuration, achieved a 5-dB transmission loss compared to the base case within the frequency range of approximately 700–1400 Hz. These results highlight the promising role of Hilbert curve–based ventilated metamaterials in mitigating sound propagation and have implications for various noise control and acoustic engineering applications.

Acknowledgment

The authors would like to thank Professor Chandramouli Padmanabhan, Dr. Jhalu Gorain, and Mr. K. Suresh for their kind help with the impedance tube and twin reverberation chamber facility. The authors are also thankful to the central skill training and fabrication facility at IIT Madras and Mr. R Vasantha Kumar for fabricating 3D printing samples.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Appendix

Additional supporting data of the cases presented in the article are given in Figs. 1518.

Fig. 15
Comparison of sound transmission loss performance between thermoviscous losses and lossless for ventilated acoustic metamaterials and base case by COMSOL simulation
Fig. 15
Comparison of sound transmission loss performance between thermoviscous losses and lossless for ventilated acoustic metamaterials and base case by COMSOL simulation
Close modal
Fig. 16
Variation of Hilbert space width (K) for first-order (a) ventilated area 15×15, (b) ventilated area 20×20, and (c) ventilated area 25×25
Fig. 16
Variation of Hilbert space width (K) for first-order (a) ventilated area 15×15, (b) ventilated area 20×20, and (c) ventilated area 25×25
Close modal
Fig. 17
Variation of thickness of unit cell (H). (a) first order, and (b) third order.
Fig. 17
Variation of thickness of unit cell (H). (a) first order, and (b) third order.
Close modal
Fig. 18
Comparison of band gap width between the band diagram and transmission loss study for second-order geometry with different numbers of unit cells
Fig. 18
Comparison of band gap width between the band diagram and transmission loss study for second-order geometry with different numbers of unit cells
Close modal

References

1.
Morillas
,
J. M. B.
,
Gozalo
,
G. R.
,
González
,
D. M.
,
Moraga
,
P. A.
, and
Vílchez-Gómez
,
R.
,
2018
, “
Noise Pollution and Urban Planning
,”
Curr. Pollut. Rep.
,
4
(
3
), pp.
208
219
.
2.
Lusk
,
S. L.
,
McCullagh
,
M.
,
Dickson
,
V. V.
, and
Xu
,
J.
,
2016
, “
Position Statement: Harmful Effects of Environmental Noise Exposures
,”
Nurs. Outlook
,
64
(
4
), pp.
395
396
.
3.
Héroux
,
M. E.
,
Babisch
,
W.
,
Belojevic
,
G.
,
Brink
,
M.
,
Janssen
,
S.
,
Lercher
,
P.
,
Paviotti
,
M.
, et al
,
2015
, “
WHO Environmental Noise Guidelines for the European Region
,”
Euronoise
,
2015
, pp.
2589
2593
.
4.
Tang
,
S. K.
,
2017
, “
A Review on Natural Ventilation-Enabling Façade Noise Control Devices for Congested High-Rise Cities
,”
Appl. Sci.
,
7
(
2
), p.
175
.
5.
Tong
,
Y. G.
, and
Tang
,
S. K.
,
2013
, “
Plenum Window Insertion Loss in the Presence of a Line Source – A Scale Model Study
,”
J. Acoust. Soc. Am.
,
133
(
3
), pp.
1458
1467
.
6.
Huang
,
H.
,
Qiu
,
X.
, and
Kang
,
J.
,
2011
, “
Active Noise Attenuation in Ventilation Windows
,”
J. Acoust. Soc. Am.
,
130
(
1
), pp.
176
188
.
7.
Mohsen
,
E. A.
, and
Oldham
,
D. J.
,
1977
, “
Traffic Noise Reduction Due to the Screening Effect of Balconies on a Building Façade
,”
Appl. Acoust.
,
10
(
4
), pp.
243
257
.
8.
Badino
,
E.
,
Manca
,
R.
,
Shtrepi
,
L.
,
Calleri
,
C.
, and
Astolfi
,
A.
,
2019
, “
Effect of Façade Shape and Acoustic Cladding on Reduction of Leisure Noise Levels in a Street Canyon
,”
Build. Environ.
,
157
, pp.
242
256
.
9.
Tang
,
S. K.
,
2012
, “
Narrow Sidebranch Arrays for Low Frequency Duct Noise Control
,”
J. Acoust. Soc. Am.
,
132
(
5
), pp.
3086
3097
.
10.
Wang
,
Z. H.
,
Hui
,
C. K.
, and
Ng
,
C. F.
,
2014
, “
The Acoustic Performance of Ventilated Window With Quarter-Wave Resonators and Membrane Absorber
,”
Appl. Acoust.
,
78
, pp.
1
6
.
11.
Murphy
,
E.
, and
King
,
E. A.
,
2014
, “Chapter 7—Noise Mitigation Approaches,”
Environmental Noise Pollution
,
E.
Murphy
and
E. A.
King
, eds.,
Elsevier
,
Boston, MA
, pp.
203
245
.
12.
Liu
,
Z.
,
Zhang
,
X.
,
Mao
,
Y.
,
Zhu
,
Y. Y.
,
Yang
,
Z.
,
Chan
,
C. T.
, and
Sheng
,
P.
, “
Locally Resonant Sonic Materials
,”
Science
,
289
(
5485
), pp.
1734
1736
.
13.
Xie
,
S.
,
Wang
,
D.
,
Feng
,
Z.
, and
Yang
,
S.
,
2020
, “
Sound Absorption Performance of Microperforated Honeycomb Metasurface Panels With a Combination of Multiple Orifice Diameters
,”
Appl. Acoust.
,
158
, p.
107046
.
14.
Long
,
H.
,
Shao
,
C.
,
Liu
,
C.
,
Cheng
,
Y.
, and
Liu
,
X.
,
2019
, “
Broadband Near-Perfect Absorption of Low-Frequency Sound by Subwavelength Metasurface
,”
Appl. Phys. Lett.
,
115
(
10
), p.
103503
.
15.
Wang
,
W.
,
2020
, “
Analysis and Optimization of a Metamaterial for Noise Reduction
,” Master of Science Thesis, Technical University of Munich, Munich, Bavaria, Germany.
16.
Bai
,
Y.
,
Wang
,
X.
,
Luo
,
X.
, and
Huang
,
Z.
,
2022
, “
Acoustic Superscatterer Enables Remote Mitigation of Underwater Source Radiation
,”
J. Appl. Phys.
,
131
(
24
), p.
244502
.
17.
Ramachandran
,
V. P.
, and
Rajagopal
,
P.
,
2023
, “
Bandwidth-Limited Passive Suppression of Cylindrical Source Radiation Using Metamaterial Based Acoustic Superscatterers
,”
J. Sound Vib.
,
560
, p.
117767
.
18.
Li
,
W.
,
Lu
,
G.
, and
Huang
,
X.
,
2022
, “
Acoustic Hologram of the Metasurface With Phased Arrays via Optimality Criteria
,”
Mech. Syst. Signal Process.
,
180
, p.
109420
.
19.
He
,
M. X.
, and
Ding
,
Q.
,
2023
, “
Dynamic Analysis and Design of Metamaterial Plates With Crossed Acoustic Black Holes for Vibration Control
,”
ASME J. Vib. Acoust.
,
145
(
1
), p.
011013
.
20.
Zheng
,
L. Y.
,
Wu
,
Y.
,
Ni
,
X.
,
Chen
,
Z. G.
,
Lu
,
M. H.
, and
Chen
,
Y. F.
,
2014
, “
Acoustic Cloaking by a Near-Zero-Index Phononic Crystal
,”
Appl. Phys. Lett.
,
104
(
16
), p.
161904
.
21.
Li
,
Q.
, and
Vipperman
,
J. S.
,
2019
, “
Two-Dimensional Arbitrarily Shaped Acoustic Cloaks With Triangular Patterns of Homogeneous Properties
,”
ASME J. Vib. Acoust.
,
141
(
2
), p.
021014
.
22.
Lee
,
H. M.
,
Wang
,
Z.
,
Lim
,
K. M.
,
Xie
,
J.
, and
Lee
,
H. P.
,
2020
, “
Novel Plenum Window With Sonic Crystals for Indoor Noise Control
,”
Appl. Acoust.
,
167
, p.
107390
.
23.
Shen
,
C.
,
Xie
,
Y.
,
Li
,
J.
,
Cummer
,
S. A.
, and
Jing
,
Y.
,
2018
, “
Acoustic Metacages for Sound Shielding With Steady Air Flow
,”
J. Appl. Phys.
,
123
(
12
), p.
124501
.
24.
Melnikov
,
A.
,
Maeder
,
M.
,
Friedrich
,
N.
,
Pozhanka
,
Y.
,
Wollmann
,
A.
,
Scheffler
,
M.
,
Oberst
,
S.
,
Powell
,
D.
, and
Marburg
,
S.
,
2020
, “
Acoustic Metamaterial Capsule for Reduction of Stage Machinery Noise
,”
J. Acoust. Soc. Am.
,
147
(
3
), pp.
1491
1503
.
25.
Kim
,
S. H.
, and
Lee
,
S. H.
,
2014
, “
Air Transparent Soundproof Window
,”
AIP Adv.
,
4
(
11
), p.
117123
.
26.
Yu
,
X.
,
Lu
,
Z.
,
Cheng
,
L.
, and
Cui
,
F.
,
2017
, “
On the Sound Insulation of Acoustic Metasurface Using a Sub-Structuring Approach
,”
J. Sound Vib.
,
401
, pp.
190
203
.
27.
Kumar
,
S.
,
Xiang
,
T. B.
, and
Lee
,
H. P.
,
2020
, “
Ventilated Acoustic Metamaterial Window Panels for Simultaneous Noise Shielding and Air Circulation
,”
Appl. Acoust.
,
159
, p.
107088
.
28.
Wang
,
X.
,
Luo
,
X.
,
Yang
,
B.
, and
Huang
,
Z.
,
2019
, “
Ultrathin and Durable Open Metamaterials for Simultaneous Ventilation and Sound Reduction
,”
Appl. Phys. Lett.
,
115
(
17
), pp.
2
5
.
29.
Xiao
,
Z.
,
Gao
,
P.
,
He
,
X.
,
Qu
,
Y.
, and
Wu
,
L.
,
2023
, “
Multifunctional Acoustic Metamaterial for Air Ventilation, Broadband Sound Insulation and Switchable Transmission
,”
J. Phys. D: Appl. Phys.
,
56
(
4
), p.
044006
.
30.
Xin
,
Y. J.
,
Huang
,
R. N.
,
Li
,
P.
,
Qian
,
Q.
,
Yan
,
Q.
,
Sun
,
Y. T.
,
Ding
,
Q.
,
Huang
,
Y. J.
, and
Cheng
,
S. L.
,
2023
, “
Labyrinthine Acoustic Metamaterials With Triangular Self-Similarity for Low-Frequency Sound Insulation at Deep Subwavelength Dimensions
,”
Results Phys.
,
54
, p.
107151
.
31.
Zhao
,
P.
,
Zhang
,
K.
,
Zhao
,
C.
, and
Deng
,
Z.
,
2022
, “
Mechanism of Band Gaps in Self-Similar Triangular Lattice with Koch Fractal
,”
ASME J. Vib. Acoust.
,
144
(
3
), p.
031009
.
32.
Chen
,
Y.
, and
Wang
,
L.
,
2015
, “
Multiband Wave Filtering and Waveguiding in Bio-Inspired Hierarchical Composites
,”
Extreme Mech. Lett.
,
5
, pp.
18
24
.
33.
Wen
,
W.
,
Zhou
,
L.
,
Li
,
J.
,
Ge
,
W.
,
Chan
,
C. T.
, and
Sheng
,
P.
,
2002
, “
Subwavelength Photonic Band Gaps from Planar Fractals
,”
Phys. Rev. Lett.
,
89
(
22
), pp.
1
4
.
34.
Xia
,
B.
,
Li
,
L.
,
Liu
,
J.
, and
Yu
,
D.
,
2018
, “
Acoustic Metamaterial With Fractal Coiling Up Space for Sound Blocking in a Deep Subwavelength Scale
,”
ASME J. Vib. Acoust.
,
140
(
1
), p.
011011
.
35.
Man
,
X.
,
Liu
,
T.
,
Xia
,
B.
,
Luo
,
Z.
,
Xie
,
L.
, and
Liu
,
J.
,
2018
, “
Space-Coiling Fractal Metamaterial With Multi-Bandgaps on Subwavelength Scale
,”
J. Sound Vib.
,
423
, pp.
322
339
.
36.
Song
,
G. Y.
,
Huang
,
B.
,
Dong
,
H. Y.
,
Cheng
,
Q.
, and
Cui
,
T. J.
,
2016
, “
Broadband Focusing Acoustic Lens Based on Fractal Metamaterials
,”
Sci. Rep.
,
6
(
1
), pp.
1
7
.
37.
Cai
,
X.
,
Huang
,
Z.
,
Wang
,
C.
,
Jia
,
P.
,
Yang
,
J.
, and
Zhang
,
L.
,
2023
, “
Acoustic Wave Manipulation by Phase Conjugate Metasurface
,”
ASME J. Vib. Acoust.
,
145
(
2
), p.
021006
.
38.
Man
,
X.
,
Luo
,
Z.
,
Liu
,
J.
, and
Xia
,
B.
,
2019
, “
Hilbert Fractal Acoustic Metamaterials With Negative Mass Density and Bulk Modulus on Subwavelength Scale
,”
Mater. Des.
,
180
, p.
107911
.
39.
Xin
,
Y. J.
,
Huang
,
R. N.
,
Li
,
P.
,
Yan
,
H.
,
Dong
,
X. J.
,
Yan
,
Q.
,
Sun
,
Y. T.
,
Cheng
,
S. L.
, and
Zhao
,
Q. X.
,
2023
, “
Labyrinth Acoustic Metamaterials With Fractal Structure Based on Hilbert Curve
,”
Physica B
,
667
, p.
415150
.
40.
Comandini
,
G.
,
Khodr
,
C.
,
Ting
,
V. P.
,
Azarpeyvand
,
M.
, and
Scarpa
,
F.
,
2022
, “
Sound Absorption in Hilbert Fractal and Coiled Acoustic Metamaterials
,”
Appl. Phys. Lett.
,
120
(
6
), p.
061902
.
41.
Ma
,
P. S.
,
Seo
,
Y. H.
, and
Lee
,
H.
,
2023
, “
Multiband Ventilation Barriers Enabled by Space-Coiling Acoustic Metamaterials
,”
Appl. Acoust.
,
211
, p.
109565
.
42.
Xiang
,
L.
,
Wang
,
G.
,
Zhu
,
C.
,
Shi
,
M.
,
Hu
,
J.
, and
Luo
,
G.
,
2022
, “
Ventilation Barrier With Space-Coiling Channels of Varying Cross-Section for Broadband Sound Insulation
,”
Appl. Acoust.
,
201
, p.
109110
.
43.
Kumar
,
S.
, and
Lee
,
H. P.
,
2020
, “
Labyrinthine Acoustic Metastructures Enabling Broadband Sound Absorption and Ventilation
,”
Appl. Phys. Lett.
,
116
(
13
), p.
134103
.
44.
Zhao
,
X.
,
Liu
,
G.
,
Zhang
,
C.
,
Xia
,
D.
, and
Lu
,
Z.
,
2018
, “
Fractal Acoustic Metamaterials for Transformer Noise Reduction
,”
Appl. Phys. Lett.
,
113
(
7
), p.
7
.
45.
Li
,
Z.
,
2018
, “
Spectral Comparison of Pass-by Traffic Noise
,”
Proceedings of INTER-NOISE 2018 – 47th International Congress and Exposition on Noise Control Engineering Impact Noise Control Engineering
,
Chicago, IL
,
Aug. 26–29
, ASME, pp.
1
8
.
46.
Buratti
,
C.
,
Belloni
,
E.
, and
Moretti
,
E.
,
2014
, “
Façade Noise Abatement Prediction: New Spectrum Adaptation Terms Measured in Field in Different Road and Railway Traffic Conditions
,”
Appl. Acoust.
,
76
, pp.
238
248
.
47.
Mesihovic
,
M.
, and
Milford
,
I.
,
2016
, “
The Need for Updated Traffic Noise Spectra, Used for Calculation of Sound Insulation of Windows and Facades
,”
Proceedings of INTER-NOISE 2016 – 45th International Congress and Exposition on Noise Control Engineering
, vol. 253, no. 1, pp.
7679
7686
.
48.
Yang
,
W.
,
Cai
,
M.
, and
Luo
,
P.
,
2020
, “
The Calculation of Road Traffic Noise Spectrum Based on the Noise Spectral Characteristics of Single Vehicles
,”
Appl. Acoust.
,
160
, p.
107128
.
49.
Kumar
,
S.
, and
Lee
,
H. P.
,
2022
, “
Acoustic Performance of Sonic Metacage Consisting of Helmholtz’s Resonator Columns With Internal Partitions
,”
Appl. Acoust.
,
196
, p.
108887
.
50.
Krasikova
,
M.
,
Krasikov
,
S.
,
Melnikov
,
A.
,
Baloshin
,
Y.
,
Marburg
,
S.
,
Powell
,
D. A.
, and
Bogdanov
,
A.
,
2023
, “
Metahouse: Noise-Insulating Chamber Based on Periodic Structures
,”
Adv. Mater. Technol.
,
8
(
1
), pp.
1
11
.
51.
Jung
,
J. W.
,
Kim
,
J. E.
, and
Lee
,
J. W.
,
2018
, “
Acoustic Metamaterial Panel for Both Fluid Passage and Broadband Soundproofing in the Audible Frequency Range
,”
Appl. Phys. Lett.
,
112
(
4
), p.
041903
.
52.
Ingard
,
U.
,
1953
, “
On the Theory and Design of Acoustic Resonators
,”
J. Acoust. Soc. Am.
,
25
(
6
), pp.
1037
1061
.
53.
Huang
,
S.
,
Fang
,
X.
,
Wang
,
X.
,
Assouar
,
B.
,
Cheng
,
Q.
, and
Li
,
Y.
,
2018
, “
Acoustic Perfect Absorbers via Spiral Metasurfaces With Embedded Apertures
,”
Appl. Phys. Lett.
,
113
(
23
), p.
233501
.
54.
Huang
,
S.
,
Zhou
,
Z.
,
Li
,
D.
,
Liu
,
T.
,
Wang
,
X.
,
Zhu
,
J.
, and
Li
,
Y.
,
2020
, “
Compact Broadband Acoustic Sink With Coherently Coupled Weak Resonances
,”
Sci. Bull.
,
65
(
5
), pp.
373
379
.
55.
Zhou
,
X.
,
Wang
,
X.
, and
Xin
,
F.
,
2023
, “
Ultrathin Acoustic Metamaterial as Super Absorber for Broadband Low-Frequency Underwater Sound
,”
Sci. Rep.
,
13
(
1
), pp.
1
10
.
56.
Zhao
,
J.
,
Zhang
,
L.
, and
Wu
,
Y.
,
2017
, “
Enhancing Monochromatic Multipole Emission by a Subwavelength Enclosure of Degenerate Mie Resonances
,”
J. Acoust. Soc. Am.
,
142
(
1
), pp.
EL24
EL29
.
57.
ISO (International Organization for Standardization)
, “ISO 10140-1:2021: Acoustics – Laboratory Measurement of Sound Insulation of Building Elements,” https://www.iso.org/obp/ui/en/#iso:std:iso:10140:-1:ed-3:v1:en
58.
COMSOL Multiphysics User Guide
,
2010
, “Acoustics Module,” Interfaces, p.
214
.
59.
Fokin
,
V.
,
Ambati
,
M.
,
Sun
,
C.
, and
Zhang
,
X.
,
2007
, “
Method for Retrieving Effective Properties of Locally Resonant Acoustic Metamaterials
,”
Phys. Rev. B
,
76
(
14
), pp.
1
5
.
60.
Nguyen
,
H. Q.
,
Wu
,
Q.
,
Chen
,
H.
,
Chen
,
J. J.
,
Yu
,
Y. K.
,
Tracy
,
S.
, and
Huang
,
G. L.
,
2021
, “
A Fano-Based Acoustic Metamaterial for Ultra-Broadband Sound Barriers
,”
Proc. R. Soc. A Math. Phys. Eng. Sci.
,
477
(
2248
), p.
20210024
.
61.
Cheng
,
Y.
,
Zhou
,
C.
,
Yuan
,
B. G.
,
Wu
,
D. J.
,
Wei
,
Q.
, and
Liu
,
X. J.
,
2015
, “
Ultra-Sparse Metasurface for High Reflection of Low-Frequency Sound Based on Artificial Mie Resonances
,”
Nat. Mater.
,
14
(
10
), pp.
1013
1019
.
62.
Yao
,
Y.
,
Hou
,
Z.
,
Wu
,
F.
, and
Zhang
,
X.
,
2011
, “
Low-Frequency Band Gaps in One-Dimensional Thin Phononic Crystal Plate With Periodic Stubbed Surface
,”
Physica B
,
406
(
11
), pp.
2249
2253
.
63.
Karthik
,
R.
,
Barathula
,
S.
, and
Srinivasan
,
K.
,
2024
, “
Data-Driven Approach for Prediction of Band Diagram of Locally Resonant Sonic Crystal
,”
J. Phys. D: Appl. Phys.
,
57
(
13
), p.
135307
.
64.
Almeida
,
G. N.
,
Vergara
,
E. F.
, and
Lenzi
,
A.
,
2023
, “
Metamaterial Based on Labyrinthine Tubes Loaded Periodically and Ultra-Low Phase Velocity to the Low-Frequency Energy Control
,”
Appl. Phys. Lett.
,
123
(
8
), p.
081706
.