Omnidirectional acoustic cloaking against airborne sound realized by a locally resonant sonic material.

We report that a locally resonant sonic material realizes omnidirectional acoustic invisibility in air. To achieve acoustic cloaking in the low-frequency regime, we axisymmetrically placed elastic rods comprised of silicone rubber and lead around a cloaked object. The radii of the rods are designed to minimize their total scattering cross section for a given frequency. The optimization is performed using an algorithm incorporating multiple scattering theory and gradient-based nonlinear programming. We numerically confirmed that the designed cloaking device suppressed the scattering cross section by almost 92% for all incident directions at the target frequency.

Introduction

Invisibility cloaking remains one of the most exotic phenomena in optics and acoustics. The primary concept behind invisibility cloaking is to cover a scatterer with carefully designed objects such that they suppress or cancel the scattering.

After Pendry et al.[1] reported that a coordinate-transformation technique can determine a material parameter distribution for perfect cloaking, many studies have been devoted to theoretical, numerical, and experimental realization of this transformation-based cloaking in acoustics[2-10]. Because the transformation optics and acoustics require a continuously inhomogeneous and anisotropic material with extreme values, its experimental realization is not easy. The scattering-cancellation technique is another approach to realize optical and acoustic cloaking[11-18]. While the scattering-cancellation approach does not require the continuous and anisotropic material distribution, higher and lower mass densities than the background medium are still required, which is very difficult for airborne sound cloaking.

For acoustic invisibility, a more feasible approach is to design a cloaking object’s geometry rather than manipulate its material properties. For example, Garcia-Chocano et al.[19] reported that aluminum rods distributed around a cylindrical object can significantly suppress its acoustic scattering by optimizing the positions of the rods. This geometry technique is not limited to cylindrical objects, but more complex geometries can be designed by optimization techniques[20-27]. However, the main drawback of geometry-based approaches is that designed devices often lack cylindrical symmetry, i.e., invisibility cloaking only works for specific incident directions.

This challenge was recently tackled by Jo et al.[28] by arranging axisymmetric cylindrical lattices around a core object and optimizing their radii and positions to suppress acoustic scattering. Because they used a rigid material for the cloaking device, the cloaking effect originates from the Bragg scattering among the lattice. This implies that the size of the device must be several times larger than the target wavelength.

In this report, we demonstrate a locally resonant sonic material that realizes omnidirectional acoustic cloaking from airborne sound in the low-frequency regime. Our cloaked object model is similar to those proposed by Garcia-Chocano et al.[19] and Jo et al.[28], i.e., a rigid circular object with surrounding elastic rods. We used silicone rubber and lead to construct the surrounding rods in order to induce local resonance in the low-frequency range[29]. We fixed the positions of the rods and optimized their radii to suppress the total scattering cross section for a given frequency. The calculation is performed using multiple scattering theory[30-32]. After the optimization, we observed that the designed cloaking device reduces the scattering cross section by almost for every direction.

Model and formulation

We define the model to achieve omnidirectional acoustic cloaking of a cylindrical object in air (mass density and bulk modulus ). As shown in Fig. 1, we consider a rigid-material cylindrical core surrounded by rods. Each rod is comprised of silicone rubber[33] (mass density , bulk modulus , and shear modulus ) encapsulating lead[34] (mass density , bulk modulus , and shear modulus ) placed at the center.

Figure 1

Rigid-material core covered with rubber-lead rods. The number of rods in each layer, from innermost to outermost, are 14, 20, 28, and 34. The materials are placed in air and illuminated by an incident acoustic wave. The radii are optimized by the proposed method. The core is placed at the origin. The centers of the rods are located at (1.25R, 0), (1.75R, 0), (2.25R, 0), and (2.75R, 0) (from innermost to outermost), where R is the radius of the core.

For a given operating frequency, we optimized the radii of each rubber-lead rod to minimize their total scattering cross section. The cross section is computed using multiple scattering theory. To this end, we assume that sound pressure p is governed by the Helmholtz equation in air, i.e.,where is a given angular frequency. Since the silicone rubber and lead have non-negligible shear moduli, we solve Navier’s equationunder the plane-strain condition to obtain the displacement u and stress in the elastic rods, where is the mass density in solid. The partial differential equations are coupled by the following interface conditions:where n is the unit outward normal vector.

Acoustic cloaking is achieved by minimizing the scattering cross section . To accomplish this, we set the design parameters as the radii of the rods aligned along each layer as shown in Fig. 1. To ensure that the rods are well-separated, we also impose the constraintswhere R is the radius of the core. The nonlinear optimization problem with linear constraints is solved by sequential least-squares quadratic programming[35, 36], implemented in NLopt[37].

Results

Designed cloaking device

In Table 1, we list the optimized radii of the rubber-lead rods when we gave a plane wave propagating in the positive direction with normalized frequency , where c is the speed of sound in air. The minimized objective value is , where is the scattering cross section of the bare core at the same frequency. From the results, we observed that the designed cloak suppresses scattering by almost 92%.

Table 1

Optimized radii of the rubber-lead rods.

i	1	2	3	4
Rout(i)/R	0.19909	0.23752	0.16915	0.16365
Rin(i)/R	0.02452	0.09350	0.00797	0.00000

To confirm that the designed cloak is omnidirectional, we varied the direction of the incident plane wave and plotted the corresponding scattering cross section in Fig. 2. From the results, we note that the scattering cross section is only minorly sensitive to incident angle, and all the values are far less than unity. This implies that the designed cloak works for any incident angle, thus achieving omnidirectional cloaking.

Figure 2

Scattering cross section of the designed cloak at various incident angles.

The omnidirectional cloak successfully hides the core object from any time-harmonic waves as they can be written as a superposition of plane waves. We placed a monopole source at and defined the incident wave as instead of the plane wave. Figure 3 illustrates the sound fields for the bare and cloaked objects. From the results, we observe that the forward and backward scattering are significantly suppressed outside the designed layers. The layers localize the incident wave around the core and magnify its amplitude. This implies that high-quality resonance contributes to the cloaking effect.

Figure 3

Acoustic field p and scattered field when a monopole source is present. The source is located at . The dots in (b) represent the positions of the rubber-lead rods.

To confirm that the cloaking is associated with the resonance, we calculated the dipolar component of each rubber-lead rod’s scattering matrix S. The results are shown in Fig. 4. At the target frequency , the innermost rod exhibits a sharp dipolar resonance whereas its neighbor has a higher resonant frequency. It is well known that a dipolar resonance contributes to a dramatic shift in macroscopic mass density, indicating that the cloaking is achieved by optimizing layer-dependent effective material parameters.

Figure 4

Dipolar component spectrum of the scattering matrix S for each rubber-lead rod.

We are also interested in how the cloaking performance depends on the operating frequency. In Fig. 5, we plotted the spectrum of the scattering cross section when a plane wave propagates in the positive direction. To check the accuracy of multiple scattering theory, we conducted the same simulation using a mode-matching finite element method[38], implemented in FreeFem++[39]. Both multiple scattering theory and the finite element method yielded consistent values, verifying the calculations in this study. In the spectrum, we observed a sharp dip at the target frequency, meaning that the cloaking only works in a narrow frequency band. However, the spectrum showed that the scattering cross section is also suppressed in another narrow band , though this range is not targeted in the optimization process. This narrow-band cloaking also appears to be associated with the dipolar resonance at in Fig. 4. In the lower frequency range, the cloaked object exhibited larger scattering cross section than the bare one. This is because the resonance effect is no longer available in this spectrum; thus, the designed rods simply increased the geometrical cross section.

Figure 5

Scattering cross section of the designed cloaked material when the upward-propagating plane wave is given for various frequencies.

Viscoelasticity

As the proposed structure utilizes the local resonance of rubber-lead cells, we are interested in whether the cloaking effect is valid even when the materials have non-negligible viscosity, which may break the low-frequency local resonance. Acoustic waves lose their energy mainly due to the viscosity of the solid materials and viscothermal boundary layers at the solid-fluid interfaces[40]. To model the material losses, we replace the elastic tensor C of the silicone rubber with the complex value , where is the original constant, and is the loss tangent.

We performed the same optimization for and list the optimized radii in Table 2. From the results, we observe that even the small viscosity changes the optimized values. In addition, we plot the spectrum of the scattering cross section for the viscoelastic system in Fig. 6. As expected, the introduced viscosity deteriorates the minimum value in the spectrum, while the cloaking bandwidth is slightly broadened. The result indicates that the cloaking effect is still observed even if a small viscosity is introduced in the system.

Table 2

Optimized radii of the rubber-lead rods for various .

i		1	2	3	4
tanδ=0.00	Rout(i)/R	0.19909	0.23752	0.16915	0.16365
	Rin(i)/R	0.02452	0.09350	0.00797	0.00000
tanδ=0.01	Rout(i)/R	0.19207	0.21153	0.19316	0.16445
	Rin(i)/R	0.00000	0.00000	0.02379	0.00000
tanδ=0.04	Rout(i)/R	0.18720	0.20698	0.19073	0.16765
	Rin(i)/R	0.00000	0.00000	0.03251	0.00000
tanδ=0.10	Rout(i)/R	0.19043	0.22659	0.18315	0.17142
	Rin(i)/R	0.00000	0.05137	0.00000	0.00000

Figure 6

Scattering cross section of the designed cloaked materials with loss tangent .

Conclusions

In conclusion, we designed a locally resonant sonic material to realize acoustic cloaking against airborne sound. Multiple scattering theory is formulated to calculate the scattering cross section of multiple elastic rods. The geometry of the elastic rods are optimized such that their cross section is minimized at the target frequency. We successfully demonstrated that the designed cloaking works for every incident direction within a narrow frequency band. Another optimization showed that the proposed cloaking system is still valid even if a small amount of loss is introduced in the solid. Future directions of this work include a comprehensive study of viscothermal losses induced by the boundary-layer effect in air. Another open question is whether the proposed system can realize broadband acoustic cloaking.

Methods

Multiple scattering theory

Multiple scattering theory expresses a solution p to the coupled scattering problem aswhere denotes the angle between x and the axis, is the wavenumber, and is the Bessel function of nth order. Under the time-harmonic assumption with , the Hankel functions represent radiating waves from each rods. The coefficients and are associated with the ith rod (where index denotes the core), whose center is , and are determined by solving a linear system of equations, called a multiple-scattering equation, of the following form:The matrix S is comprised of scattering matrices of the rubber-lead rods, which depend on the configurations (radii, material parameters, and frequency) of the rods, while the matrix T is a function of only the rods’ positions and frequency. The vector b is determined by the incident wave[32].

Mode-matching finite element method

The finite element simulation was conducted by the open source software FreeFem++[39]. To solve the exterior coupled problem, FreeFem++ discretizes the variational equationin air region with elastic counterpartfor test functions and , where is the strain field induced by the displacement , and represents the rubber-lead rods surrounding the core. The two formulations are coupled by the interface conditions (3) and (4). To impose the Sommerfeld radiation condition, we truncated the unbounded air region by a fictitious circle of radius . The radius is large enough to enclose the rods and core. Outside the fictitious disk , we have the following multipole expansionwith given coefficients (incident wave). The unknown coefficients are obtained by matching the sound pressure p and its normal flux with the interior solution on . For details, see the reference[38].