# A Novel Approach to Density Near-Zero Acoustic Metamaterials

Authors: Milan Sečujski; Norbert Cselyuszka; Vesna Crnojević-Bengin
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

The study demonstrates the possibility of achieving near-zero propagation of sound waves in acoustic metamaterials based on a membrane-based metamaterial unit cell which exhibits effective mass density of Lorentzian type. The unit cell, which represents the acoustic counterpart of the split ring resonator, was previously used as a building block of left-handed metamaterials, as it exhibits negative density at certain frequencies. In this study we show that its application can be extended to achieving propagation of sound waves at a frequency where its effective density equals zero. This effect can be exploited in a range of applications where extremely low phase variation over long physical distances is required, such as energy tunneling or tailoring the acoustic radiation phase pattern in arbitrary ways. After discussing the dependence of the frequency response of the unit cell on the properties of the host, we show that it can be used to design near-zero acoustic filters with low insertion loss and steep roll-off. Finally, we show that it can be used to achieve simultaneous near-zero propagation at multiple, independently chosen frequencies.


---

## Body

## 1. Introduction

Acoustic metamaterials are artificially designed mechanical structures used to control the propagation of sound waves in ways generally not possible in conventional media. They typically include periodic inhomogeneities whose dimensions are sufficiently small in comparison to the wavelength, allowing the material to be treated as homogenous and described by its effective parameters. As opposed to electromagnetic (EM) metamaterials, where specific design is adopted to manipulate the values of effective electric permittivity or magnetic permeability, the research in the field of acoustic metamaterials has principally focused on obtaining metamaterials with unconventional values of compressibilityβ or mass density ρ [1–7].In our previous research [8], we proposed a novel acoustic metamaterial unit cell based on the analogy with the split ring resonator (SRR), an element regularly used to obtain negative magnetic permeability in EM metamaterials. The proposed acoustic metamaterial unit cell (acoustic split ring resonator, aSRR) consisted of a thin elastic membrane placed in a short tube coaxial with an acoustic duct, as shown in Figure 1. Its mechanical layout was defined so as to match its lumped-element equivalent electrical circuit to the one of a typical resonant-permeability EM unit cell consisting of an SRR-loaded microstrip, given in [9]. The effective density of the proposed unit cell was shown to be equal to(1)ρeff=Ae+Ai1/ρ0Ae+1/ρ′-8πτ/ω2AidAi,where Ai=r2π is the cross-sectional area of the tube (internal duct), Ae=(R2-r2)π is the remaining cross-sectional area of the duct external to the tube, d is tube length, τ is membrane tension, ρ0 is air density, and ρ′=ρ01+M/ρ0Aid is the average density in the tube, where M is the mass of the membrane [8]. Unlike other membrane-based approaches to achieving effective mass density which is not always positive [4–7], where effective density exhibits Drude type behaviour, the proposed approach resulted in an effective density of the Lorentz type. In this paper we show that the same unit cell, with an appropriate choice of host medium, can be used to achieve near-zero (NZ) acoustic wave propagation at either of the two frequencies where its effective density, given by (1), is equal to zero.Figure 1
Acoustic split ring resonator (aSRR). (a) Layout and equivalent circuit. (b) Cross-section.
(a)
(b)


## 2. Effective Density from Transmission Line Model and Simulation

In order to validate the analytically derived expression for effective density given by (1), the same parameter has been extracted from a transmission line (TL) model as well as a finite element method (FEM) simulation in COMSOL Multiphysics using the approach analogous to the one developed for EM metamaterials [10, 11]. The values of the parameters used in both approaches were M=0.88mg,  d=18mm,  R=4mm,  r=3mm, and τ=5N/m, and the resulting resonant frequency was fr=1175Hz.In the approach based on TL modeling, each acoustically short section of the aSRR is modeled with series impedanceZ and shunt admittance Y or, alternatively, as a two-port network with the transfer matrix given by(2)pivi=cosh⁡γlZ0sinh⁡γlZ0-1sinh⁡γlcosh⁡γlpovo,where pi, po, vi, and vo are input and output pressures and velocities, l is the length of the section, Z0=Z/Y is its characteristic impedance, and γ=ZY is the propagation constant. The losses due to air viscosity and thermal loss are negligible, but the narrow pipe effect has to be taken into account [12–14]. The series impedance and the shunt admittance in this case are given by(3)Z=jωρ0A1-Fv-1,Y=jAωρ0c21+γ-1Ft,where A is the relevant cross-section area, ω is the angular frequency, ρ0 is the density of air, c is the sound velocity in air, while Fv and Ft are given by (4)Fv=2rv-jJ1rv-jJ0rv-j,Ft=2rt-jJ1rt-jJ0rt-j. In (4) Jm is the Bessel function of the first kind, order m, and rv and rt are given by(5)rv=rρωμ,rt=rvμCpλ=σrv,where Cp is specific thermal capacity at a constant pressure, λ is thermal conductivity, μ is shear viscosity, and σ is the square root of the Prandtl number. The thin membrane is modelled as the impedance:(6)ZM=-jωmS2I1kmrJ0kmr+J1kmrI0kmrI1kmrJ2kmr-J1kmrI2kmr,where r is the membrane radius, Jn and In are the regular and the modified Bessel function of nth order, m is the mass of the membrane,S its surface area, while km is the wave number obtained from flexural wave equation:(7)km2=ωρm′′D.In (7) ρm′′ is the surface mass density of the membrane and D is its flexural rigidity, equal to(8)D=Eh3121-v2,where E is Young’s modulus, h is the thickness of the membrane, and v is Poisson’s coefficient. In our TL model and FEM simulations we adopted the values given for DuPont Kapton FPC in COMSOL Multiphysics (E=2.758GPa, ρm=1420kg/m3, v=0.34), while the thickness of the membrane was set to h=22μm, yielding the previously specified value of membrane mass M.Figure2(a) shows a comparison between the effective density obtained from (1) and the same parameter obtained from TL model and FEM simulation. A good match between the results obtained from analytical expression, the model and the simulation can be seen, having in mind that the difference in the magnitudes can be attributed to the fact that (1) refers to the lossless case.Figure 2
(a) Effective mass density of the aSRR obtained from the analytical expression, TL model, and FEM simulation. (b) Transmission coefficient obtained from TL model and FEM simulation.
(a)
(b)


## 3. Near-Zero Wave Propagation

Effective density as a function of frequency exhibits resonant behaviour, and there are two frequencies,f01<fr and f02>fr, at which its value is equal to zero. From the simulated transmission coefficient, shown in Figure 2(b), it can be seen that the frequency response of the aSRR exhibits a notch around fr; that is, the structure behaves as a single-negative metamaterial and does not support propagating waves. For that reason, near-zero propagation at f01 is characterized by very high attenuation. At f02 attenuation is significantly lower, but the near-zero character of the propagation at f02 is not obvious because it is not easily distinguishable from conventional propagation at frequencies immediately above. It can be concluded that, for the purpose of clearly demonstrating near-zero propagation, an all-pass acoustic duct is not a convenient host medium. A convenient host should exhibit a stopband at the frequency range of interest and it should not support double negative (left-handed) wave propagation. In this specific case, where effective densitycan have both positive and negative values, this implies that the effective compressibility βeff should be strictly positive. An example of such a host is a simple acoustic low-pass filter, either CLC or LCL, provided that the aSRR operates above its cut-off frequency.In order to demonstrate density near-zero (DNZ) wave propagation we opted for a CLC filter as a host and loaded its inductive section with one aSRR, as shown in Figure3(a) (not to scale). A TL model was created, similar to the one described in the previous section, and a corresponding FEM simulation was also carried out in COMSOL Multiphysics. The resulting transmission coefficient is shown in Figure 3(b), where a very good match between the TL model and FEM simulation can be seen. The dashed dotted green line indicates the frequency response of the host low-pass filter. The transmission peak in the stopband of the host corresponds to one of the frequencies where effective density is equal to zero, namely, to f02. The frequency response exhibits a transmission zero as well, and its frequency corresponds to the notch in the frequency response of the aSRR itself. It should be noted that in this case the frequency of the transmission peak is higher than the frequency of the transmission zero.Figure 3
Low-pass acoustic CLC filter (C sections ⌀52 × 20 mm, L section ⌀10 × 80 mm) loaded with one aSRR in the L section. (a) Layout. (b) Frequency response obtained through TL model and FEM simulation.
(a)
(b)It is also possible to demonstrate DNZ wave propagation by loading each of the two capacitive sections of the CLC host with an aSRR, as shown in Figure4(a). A transmission peak and a transmission zero appear in the stopband of the host filter again, their frequencies corresponding to the zero value of the effective density and to the notch in the frequency response of the aSRR, respectively. However, in this case, the frequency of the transmission peak corresponds to f01, and it is lower than the frequency of the transmission zero, as shown in Figure 4(b).Figure 4
(a) Low-pass acoustic CLC filter loaded with two aSRRs in its capacitive sections. (a) Layout. (b) Frequency response.
(a)
(b)It should be noted that when the unit cell is placed in a host structure, its effective density changes because the space between the unit cell and the wall of the host acts like an additional external duct. This includes an additional parallel inductance in the equivalent circuit and eventually results in a shift of the resonant frequency of the aSRR, visible in the insets of both Figures3(b) and 4(b). Consequently, the actual values of f01 and f02 are also different than those indicated in Figure 2(b). The shift in the case of the CLC host is more noticeable because the inductance of the additional external duct is significantly smaller than the inductance of the actual external duct of the aSRR.It can also be shown that the situation does not change significantly if both inductive sections of an LCL low-pass host or just the single inductive section of a CLC low-pass host are loaded with an aSRR. Whether the transmission peak will be at a higher frequency than the transmission zero or vice versa depends only on whether the aSRR is located in the inductive or the capacitive section of the host filter.


## 4. Near-Zero Bandpass Filtering

The versatility of the acoustic split ring resonator will be illustrated with an example of its use in NZ bandpass filtering. Namely, it has been shown that in case aSRR is located in the inductive section of the host filter, the transmission peak is at a higher frequency than the transmission zero, while in case it is located in the capacitive section of the host filter, the opposite holds. Let us suppose that all three sections of a CLC host are loaded with one aSRR each, as shown in Figure5(a), but that the dimensions of the aSRRs are slightly different, chosen so as to set their transmission peaks at the same frequency. The aSRR in the inductive section is expected to introduce a transmission zero at a frequency immediately below the transmission peak, while the aSRRs in the capacitive sections are expected to introduce transmission zeros at a frequency immediately above. The frequency response of the entire structure will thus correspond to a bandpass filter with very small fractional bandwidth and a steep roll-off. This expectation has been confirmed by a FEM simulation, whose results are shown in Figure 5(b). The bandwidth of the obtained filter is equal to 30 Hz (corresponding to a fractional bandwidth of 2.4%), while the insertion loss is as low as 0.2 dB.Figure 5
Near-zero bandpass filter. (a) Layout (dimensions are the same as in the previous case except ford1=18.95mm for the aSRR in the inductive section of the host, and d2=17.4mm for the aSRR in its capacitive sections). (b) Frequency response (FEM simulation).
(a)
(b)It should be noted that, unlike left-handed structures, NZ filters do not suffer from inherently high passband insertion losses. Namely, from Figure2(b), it can be seen that the value of the transmission coefficient at the near-zero frequency is relatively high (i.e., the losses are relatively low). On the other hand, it exhibits significantly lower values in the frequency range where the effective density is negative, that is, at frequencies where double negative devices operate, which results in high insertion losses.Previous research [9] suggests that the mechanisms for reducing insertion loss, increasing the fractional bandwidth or achieving multiband propagation, can be easily extended to the case of DNZ propagation as well. For example, a dual-band aSRR, shown in Figure 6(a), has been shown to exhibit zero density at two independently engineered frequencies, and if placed in an otherwise all-pass acoustic duct, two notches in the frequency response can be observed. Here we show that the same unit cell can support simultaneous near-zero propagation at two frequencies if placed in a convenient host, such as the L section of the CLC filter used in the previous simulations. To that aim, a FEM simulation was carried out, with results shown in Figure 6(b), and it can be seen that two narrow passbands emerge in the stopband of the host. The dimensions used in the simulation were R=5mm, r1=3.2mm, r2=3mm, d1=32mm, d2=4mm, M1=1mg, and M2=0.88mg, setting the two near-zero frequencies at f1=856Hz and f2=1566Hz. At these two frequencies the effective density is equal to zero, as shown in the inset of Figure 6(b), which confirms the NZ character of acoustic wave propagation.Figure 6
(a) Dual-band aSRR. (b) Near-zero wave propagation at two independently chosen frequencies, obtained by loading the inductive section of a CLC host with a dual-band aSRR (FEM simulation).
(a)
(b)


## 5. Conclusions

The paper proposes a novel approach to achieving density near-zero propagation of acoustic waves. A resonant unit cell previously employed for obtaining a double negative metamaterial is used to achieve propagation of acoustic waves with virtually no change in phase over physically long distances. In this way the acoustic energy does not propagate but rather tunnels through a NZ metamaterial. Having in mind the wide range of application of NZ metamaterial structures, including energy tunneling, energy splitting, and perfect transmission through waveguides with sharp bends as well as various types of wave front manipulation, this research undoubtedly offers a valuable contribution in this field. Our future work will include the fabrication of the proposed unit cell and experimental verification of density near-zero propagation, as well as investigation of more complex structures based on the demonstrated phenomenon.



---

*Source: 626593-2015-09-30.xml*