Acoustic cloaking theory

Andrew N Norris


An acoustic cloak is a compact region enclosing an object, such that sound incident from all directions passes through and around the cloak as though the object was not present. A theory of acoustic cloaking is developed using the transformation or change-of-variables method for mapping the cloaked region to a point with vanishing scattering strength. We show that the acoustical parameters in the cloak must be anisotropic: either the mass density or the mechanical stiffness or both. If the stiffness is isotropic, corresponding to a fluid with a single bulk modulus, then the inertial density must be infinite at the inner surface of the cloak. This requires an infinitely massive cloak. We show that perfect cloaking can be achieved with finite mass through the use of anisotropic stiffness. The generic class of anisotropic material required is known as a pentamode material (PM). If the transformation deformation gradient is symmetric then the PM parameters are explicit, otherwise its properties depend on a stress-like tensor that satisfies a static equilibrium equation. For a given transformation mapping, the material composition of the cloak is not uniquely defined, but the phase speed and wave velocity of the pseudo-acoustic waves in the cloak are unique. Examples are given from two and three dimensions.


1. Introduction

The observation that the electromagnetic equations remain invariant under spatial transformations is not new. Ward & Pendry (1996) used it for numerical purposes, but the result was known to Post (1962) who discussed it in his book, and it was probably known far earlier. The recent interest in passive cloaking and invisibility is due to the fundamental result of Greenleaf et al. (2003a,b) that singular transformations could lead to cloaking for conductivity. Not long after this important discovery Leonhardt (2006) and Pendry et al. (2006) made the key observation that singular transformations could be used to achieve cloaking of electromagnetic waves. These results and others have generated significant interest in the possibility of passive acoustic cloaking.

Acoustic cloaking is considered here in the context of the so-called transformation or change-of-variables method. The transformation deforms a region in such a way that the mapping is one-to-one everywhere except at a single point, which is mapped into the cloak inner boundary (figure 1). The acoustic problem is for the infinitesimal pressure p(x, t) that satisfies the scalar wave equation in the surrounding fluid,Embedded Image(1.1)The basic idea is to alter the cloak's acoustic properties (density and modulus) so that the modified wave equation in ω mimics the exterior equation (1.1) in the entire region Ω. This is achieved if the spatial mapping of the simply connected region Ω to the multiply connected cloak ω has the property that the modified equation in ω when expressed in Ω coordinates has exactly the form of (1.1) at every point in Ω.

Figure 1

The undeformed simply connected region Ω is transformed by the mapping Χ into the multiply connected cloak ω. Essentially, a single point O is transformed into a hole (the invisible region) surrounded by the cloak ω. The outer boundary ∂ω+ is coincident with ∂Ω+(=∂Ω) and the inner boundary ∂ω is the image of the point O. Apart from O and ∂ω the mapping is everywhere one-to-one and differentiable.

The objective here is to answer the question: what type of material is required to realize these unusual properties that make an acoustic cloak? While cloaking cannot occur if the bulk modulus and density are simultaneously scalar quantities (see below), it is possible to obtain acoustical cloaks by assuming that the mass density is anisotropic (Chen & Chan 2007; Cummer & Schurig 2007; Cummer et al. 2008). A tensorial density is not ruled out on fundamental grounds (Milton et al. 2006) and in fact there is a strong physical basis for anisotropic inertia. For instance, Schoenberg & Sen (1983) showed that the inertia tensor in a medium comprising alternating fluid constituents is transversely isotropic (TI) with elements 〈ρ〉 in the direction normal to the layering, and 〈ρ−1−1 in the transverse direction, where 〈.〉 is the spatial average. Anisotropic effective density can arise from other microstructures, as discussed by Mei et al. (2007) and Torrent & Sánchez-Dehesa (2008). The general context for anisotropic inertia is the Willis equations of elastodynamics (Milton & Willis 2007), which Milton et al. (2006) showed are the natural counterparts of the electromagnetic (EM) equations that remain invariant under spatial transformation. The acoustic cloaking has been demonstrated, theoretically at least, in both two and three dimensions: a spherically symmetric cloak was discussed by Chen & Chan (2007) and Cummer et al. (2008), while Cummer & Schurig (2007) described a two-dimensional cylindrically symmetric acoustic cloak. These papers use a linear transformation based on prior EM results in two dimensions (Schurig et al. 2006).

The cloaks based on anisotropic density in combination with the inviscid acoustic pressure constitutive relation (bulk modulus) will be called inertial cloaks (ICs). The fundamental mathematical identity behind the ICs is the observation of Greenleaf et al. (2007) that the scalar wave equation is mapped into the following form in the deformed cloak region:Embedded Image(1.2)Here g=(gij) is the Riemannian metric with Embedded Image and Embedded Image. The reader familiar with differential geometry will recognize the first term in equation (1.2) as the Laplacian in curvilinear coordinates. Comparison of the transformed wave equation (1.2) with the IC wave equation provides explicit expressions for the IC density tensor and the bulk modulus (Greenleaf et al. 2008).

We will derive an identity equivalent to (1.2) in §2 using an alternative formulation adapted from the theory of finite elasticity. A close examination of the anisotropic density of the ICs shows that its volumetric integral, the total mass, must be infinite for perfect cloaking. This raises grave questions about the usefulness of the ICs. The rest of this paper provides a solution to this quandary. The main result is that the IC is a special case of a more general class of the acoustic cloaks, defined by anisotropic inertia combined with anisotropic stiffness. The latter is obtained through the use of the pentamode materials (PMs; Milton & Cherkaev 1995). In the same way that an ideal acoustic fluid can be defined as the limit of an isotropic elastic solid as the shear modulus tends to zero, there is a class of limiting anisotropic solids with five (hence penta) easy modes of deformation analogous to shear, and one non-trivial mode of stress and strain. The general cloak comprising PM and IC is called the PM-IC model. The additional degrees of freedom provided by the PM-IC allow us to avoid the infinite mass dilemma of the IC.

We begin in §2 with a new derivation of the IC model, and a discussion of the infinite mass dilemma. The PMs are introduced in §3 where it is shown that they display simple wave properties, such as an ellipsoidal slowness surface. The intimate connection between the PM and the acoustic cloaking follows from theorem 4.2 in §4. The properties of the generalized PM-IC model for cloaking are developed in §4 through the use of an example cloak that can be either pure IC or pure PM as a parameter is varied. Further examples are given in §5, with a concluding summary of the generalized acoustic cloaking theory in §6.

2. The IC

The transformation from Ω to ω is described by the point-wise deformation from XΩ to Embedded Image. In the language of finite elasticity, X describes a particle position in the Lagrangian or undeformed configuration and x is particle location in the Eulerian or deformed physical state. The transformation or mapping defined by Χ is one-to-one and invertible except at the single point X=O (figure 1). We use ∇, ∇X and div, Div to indicate the gradient and divergence operators in x and X, respectively. The component form of div A is ∂Ai/∂xi or ∂Aij/∂xi when A is a vector or a second-order tensor-like quantity, respectively. The deformation gradient is defined as F=∇Xx with inverse F−1=∇X, or in component form FiI=∂xi/∂XI and Embedded Image. The Jacobian of the deformation is Embedded Image or, in terms of volume elements in the two configurations, J=dv/dV. The polar decomposition implies F=VR, where R is proper orthogonal (RRt=RtR=I, det R=1) and the left-stretch tensor V∈Sym+ is the positive definite solution of V2=FFt. The analysis is as far as possible independent of the spatial dimension d, although applications are restricted to d=2 or 3.

The principal result for the IC is given in lemma 2.1.

Embedded Image(2.1)

The r.h.s. can be expressed asEmbedded Image(2.2)Using the chain rule in the form Embedded Image or Embedded Image implies that Embedded Image, which is Embedded Image. The proof follows from the identity (see problems 2.2.1 and 2.2.3 in Ogden 1997):Embedded Image(2.3)

(a) Cloak acoustic parameters

The connection with acoustics is made by identifying the field variable p in lemma 2.1 as the acoustic pressure. The cloak comprises an inviscid fluid with bulk modulus K(x), such that the pressure satisfies the standard relationEmbedded Image(2.4)where v(x, t) is the particle velocity. The IC is defined by the assumption that the momentum balance involves a symmetric second-order inertia tensor ρ according toEmbedded Image(2.5)Although this is a significant departure from classical acoustical theory in assuming an anisotropic mass density, it is by no means unprecedented. Based on the analysis of Schoenberg & Sen (1983), a spatially varying tensor ρ could possibly be achieved by small pockets of layered fluid separated by massless impermeable membranes.

Eliminating the velocity between equations (2.4) and (2.5) gives a single equation for the pressure,Embedded Image(2.6)Consider the uniform wave equation in Ω,Embedded Image(2.7)Using lemma 2.1, we can express this in the deformed physical description as equation (2.6), where the bulk modulus and inertia tensor areEmbedded Image(2.8)For a given deformation F, the identities (2.8) define the unique cloak with spatially varying material parameters K and ρ each defined by the deformation gradient. We note the following identity that is independent of F:Embedded Image(2.9)

Could the cloak possibly have isotropic density? That is, could the cloak be described by a standard acoustic fluid with two scalar parameters, density and bulk moduli? The identity ρ=JV−2 means that ρ=ρI can occur only if V is a multiple of the identity, V=wI for some scalar w=w(x). The deformation of Ω into the smaller region ω could certainly be accomplished at some but not all points by this deformation, which corresponds to a uniform contraction or expansion, with rotation. However, the deformation near the inner surface of the cloak cannot be of this form. In fact, the deformation in the neighbourhood of X=O must be extremely non-uniform and anisotropic. We will discuss this below when we examine a fundamental and severe deficiency of the IC model.

(b) Continuity between the cloak and the acoustic fluid

Let ds, n and dS, N denote the area element and unit normal to the outer boundary ∂ω+ and ∂Ω+(=∂ω+), respectively. These are related by the deformation through Nanson's formula (Ogden 1997), Embedded Image. The nature of the cloak requires that the outer surface is identical in either description, since both must match with the exterior fluid. We, therefore, require that ds=dS at every point on the outer surface, orEmbedded Image(2.10)and equation (2.8) then implies thatEmbedded Image(2.11)Equation (2.11) is a purely kinematic condition.

The interior of the cloak mimics the wave equation in the exterior fluid. The final requirement that the cloak will be acoustically ‘invisible’ is that the pressure and normal velocity match across the outer surface separating the fluid and cloak. These two continuity conditions arise from the balance of force (normal traction) per unit area and the constraint of particle continuity. The condition for pressure is simply that p is continuous across the outer surface, whether one uses the wave equation in physical space, (2.6), or its counterpart in the undeformed simply connect region (2.7). As for the kinematic condition consider its equivalent, the continuity of normal acceleration. This is Embedded Image in physical space, and using equation (2.5) it becomes Embedded Image, which must match with Embedded Image in the fluid. Alternatively, equation (2.11) and the relation Ft∇=∇X imply, as expected, thatEmbedded Image(2.12)The final term is simply the normal acceleration in the undeformed description.

In summary, the continuity conditions at the outer surface in the physical description areEmbedded Image(2.13)

(c) Example: a rotationally symmetric IC

Consider the inverse deformationEmbedded Image(2.14)where Embedded Image and r=|x|. Using F−1=∇X implies thatEmbedded Image(2.15)where Embedded Image and the second-order tensors are Embedded Image and Embedded Image. The bulk modulus and mass density in the cloak follow from equation (2.8) as:Embedded Image(2.16)The anisotropic inertia has the formEmbedded Image(2.17)where the radial and azimuthal principal values ρr and ρ can be read off from equation (2.16) as functions of f.

Introducing the radial and azimuthal phase speeds, Embedded Image and Embedded Image, the mass density tensor can then be expressed as Embedded Image. The quantity r is the square of the radial acoustic impedance, Embedded Image. Equation (2.9) implies that the identity Embedded Image is required for cloaking. The three equations (2.16) for K, ρr and ρ in terms of f can be replaced by the universal relation (2.9), i.e.Embedded Image(2.18)along with simple expressions for the wave speeds in terms of f,Embedded Image(2.19)We will see later that the phase and the wave (group velocity) speeds in the principal directions are identical. Note that f′ is required to be positive. The original quantities can be expressed in terms of the phase speeds asEmbedded Image(2.20)One could, for instance, eliminate f as the fundamental variable defining the cloak in favour of c(r), from which all other quantities can be determined from the differential equation relating the speeds, Embedded Image.

We assume that the cloak occupies Embedded Image with uniform acoustical properties K=1 and ρ=I in the exterior. The areal matching condition (2.11) with Embedded Image is satisfied by F and ρ of equations (2.15) and (2.16) if f is continuous across the boundary, which is accomplished by requiring f(b)=b. The pressure and velocity continuity conditions (2.13) becomeEmbedded Image(2.21)

Note that the cloak density is isotropic if cr=c, which requires that f′=f/r. Thus f=γr with γ constant, but the outer boundary condition f(b)=b implies that γ=1, which is the trivial undeformed configuration.

Perfect cloaking requires that f vanish at r=a. It is clear that c blows up as ra, as does the product r. In order to examine the individual behaviour of K and ρr, consider f∝(ra)α near a for α constant and non-negative. No value of α>0 will keep the radial density ρr bounded, although the unique choice α=1/d ensures that the bulk modulus K(a) remains finite and non-zero. Note that the azimuthal density ρ has a finite limit in two dimensions for power law decay f∝(ra)α, while ρ remains finite in three dimensions if α≤1, otherwise it blows up. Similarly, the radial phase speed scales as cr∝(ra)1−α, which remains finite for α≤1, blowing up otherwise. These results are summarized in table 1.

View this table:
Table 1

Behaviour of quantities near the inner surface r=a for the scaling fξα as ξ=ra↓0. (The total radial mass mr is defined in equation (2.22).)

We use a non-dimensional measure of the total mass in the cloak, Embedded Image. The total mass is isotropic for the symmetric deformation and configuration considered here, Embedded Image, where Embedded Image. Assuming for the moment that f(a) is non-zero, i.e. a near-cloak (Kohn et al. 2008), thenEmbedded Image(2.22)These forms indicate not only that mr→∞ as f(a)→0 but also the form of the blow-up. To leading order, Embedded Image and Embedded Image in two and three dimensions, respectively. The blow-up of mr occurs no matter how f tends to zero. The infinite mass is an unavoidable singularity.

(d) A massive problem with inertial cloaking

Table 1 and the example above illustrate a potentially grievous issue: infinite mass is required for perfect cloaking in the IC model. We now show that the problem is not specific to the rotationally symmetric cloak but is common to all ICs. Consider a ball of radius ϵ around X=O. Its volume dV=O(ϵd) is mapped to a volume with inner surface defined by the finite cloak inner boundary ∂ω and outer surface by a distance O(ϵβ) further out, where β>0 is a local scaling parameter, assumed constant (in terms of the example above and table 1, β=1/α). The mapped current volume is then dv=O(ϵβ) so that J=dv/dV=O(ϵβd). The eigenvalues of V are λ1=O(ϵβ−1), λ2, …, d=O(ϵ−1). The bulk modulus and the principal values of the density matrix are thereforeEmbedded Image(2.23)The principal value ρ1 blows up whether d=2 or d=3. Furthermore, the total mass associated with ρ1 in the mapped volume is m1=O(ϵ2−d). This blows up in three dimensions, and a more careful analysis for two dimensions similar to that for the rotationally symmetric case shows Embedded Image.

In summary, the IC theory, while consistent and formally sound, reveals an underlying and ‘massive’ problem. We will show how this can be circumvented by using a more general cloaking theory that allows for anisotropic stiffness (elasticity) in addition to, or instead of, the anisotropic inertia. The anisotropic elastic material required is of a special type, called a PM (Milton & Cherkaev 1995), which is introduced next.

3. Pentamode materials

We consider Hooke's law in three dimensions in the form Embedded Image, where the six-vectors of stress and strain, and the associated 6×6 matrix of moduli areEmbedded ImageThe Embedded Image terms ensure that products and norms are preserved, e.g. Embedded Image.

The PM is rank one, or in other words, five of the six eigenvalues of Embedded Image vanish (Milton & Cherkaev 1995). The one remaining positive eigenvalue is thereforeEmbedded Image(3.1)Accordingly, the moduli can be defined by the stiffness Embedded Image and a normalized six-vector Embedded Image,Embedded Image(3.2)The stress is described by a single scalar, Embedded Image with Embedded Image, and Embedded Image. Thus,Embedded Image(3.3)The PM (Milton et al. 2006) is so named because there are five easy ways to deform it, associated with the eigenvectors of the five zero eigenvalues of the elasticity stiffness. Pentamodes obviously include isotropic acoustic fluids, for which the only stress–strain eigenmode is a hydrostatic stress, or pure pressure, and the five easy modes are all pure shear. Milton & Cherkaev (1995) describe how PMs can be realized from specific microstructures.

(a) Example: an orthotropic PM

An elastic material with orthotropic symmetry has nine non-zero elements in general: the six Cij=Cji, i, j=1, 2, 3, plus C44, C55 and C66. We set these last three (shear) moduli to zero. The stress σ must then be diagonal in the Cartesian coordinate system, implying that Embedded Image, and therefore,Embedded Imagewith the following relations holding: Embedded Image, Embedded Image, Embedded Image.

(b) Compatibility condition for PMs

The notation Embedded Image and Embedded Image is used to signify the fact that the Embedded Image tensors are normalized by Embedded Image and therefore Embedded Image is given by equation (3.1). We will not follow this normalization in general, but write:Embedded Image(3.4)In other words, the products in (3.4) are the important physical quantities, not K and S individually. The stress in the PM is always proportional to the tensor S and only one strain element is significant, S : ϵ. The rank deficiency of the moduli, which is apparent from (3.2) or (3.4), means that there is no inverse strain–stress relation for the elements of ϵ in terms of the elements of σ.

Static equilibrium of a PM under an applied load leads to a constraint on the spatial variability of the PM stiffness. Consider an inhomogeneous PM with smoothly varying C(x)=K0(x)S0(x)⊗S0(x). Under an applied static load the strain will also be spatially inhomogeneous, but the only part of the strain that is important is the component along the PM eigenvector. With no loss in generality, we may put Embedded Image for some scalar function w. The stress is then σ=qS0, where Embedded Image. Let S=qS0, then the static equilibrium condition Embedded Image becomes Embedded Image. Finally, the PM stiffness is C=KSS where Embedded Image.

The fourth-order stiffness of a smoothly varying PM can always be expressed as C=KSS, where K(x)>0 and S(x)∈Sym satisfies the static equilibrium condition,Embedded Image(3.5)

This identity also arises in a completely different manner later when we consider transformed wave equations. We say that the PM is of canonical form when equation (3.5) applies. The decomposition of lemma 3.1 is unique up to a multiplicative constant. Thus, if a static load is applied to a PM expressed in canonical form, then the stress and strain are σ(x)=c0S and Embedded Image, respectively, for constant c0.

In summary, stability under static loading places a constraint on the PM moduli, which will turn out to be useful when we return to the cloaking problem. The constraint means that the moduli can in general be expressed in canonical form.

(c) Dynamic equations of motion in a PM

The equations for small amplitude disturbances in a PM with anisotropic mass density areEmbedded Image(3.6)andEmbedded Image(3.7)These are, respectively, the specific form of Hooke's law for a PM and the momentum balance incorporating the inertia tensor. In order to make the equations look similar to those for an acoustic fluid, we identify the ‘pseudo-pressure’ p with the negative single stress, Embedded Image. The stress tensor then becomesEmbedded Image(3.8)and the linear constitutive relation can be written asEmbedded Image(3.9)Equations (3.7) and (3.9) imply that the pseudo-pressure satisfies the generalized acoustic wave equation,Embedded Image(3.10)This reduces to the acoustic equation (2.6) with anisotropic inertia and isotropic stiffness when S=I. Finally, assuming that the PM is in canonical form, so that S satisfies the equilibrium condition (3.5), we haveEmbedded Image(3.11)

(d) Wave motion in a PM

The wave properties of PMs are of interest since we will show that they can be used to make the acoustic cloak. Consider plane wave solutions for displacement of the form Embedded Image, for |n|=1 and constant q, k and v, and uniform PM properties. Non-trivial solutions of the equations of motion (3.6) and (3.7) must satisfyEmbedded Image(3.12)The acoustical or Christoffel (Musgrave 2003) tensor K(Sn)⊗(Sn) is rank one and it follows that of the three possible solutions for v2, only one is not zero, the quasi-longitudinal solution,Embedded Image(3.13)The slowness surface is therefore an ellipsoid. Standard arguments for waves in anisotropic solids (Musgrave 2003) show that the energy flux velocity (or wave velocity or ray direction) isEmbedded Image(3.14)Note that this is in the direction Sq, and satisfies Embedded Image, a well-known relation for generally anisotropic solids with isotropic density.

As an example, consider the orthotropic PM with a density tensor of the same symmetry and coincident principal axes. ThenEmbedded Image(3.15a)Embedded Image(3.15b)Embedded Image(3.15c)where Embedded Image, Embedded Image and Embedded Image, and ρ1, ρ2 and ρ3 are the principal inertias.

4. The general acoustic cloaking theory

We now show that the IC is but a special case of a much more general type of acoustic cloak. While the IC depends upon the anisotropic inertia, the general cloaking model can have both anisotropic inertia and stiffness. The additional degree of freedom is obtained by replacing the pressure field with the scalar stress of a PM. The general cloaking model is called PM-IC.

(a) The fundamental identity

Let P∈Sym be non-singular and F is the deformation gradient for the mapping Xx with Embedded Image and V2=FFt. ThenEmbedded Image(4.1)if P satisfiesEmbedded Image(4.2)

The proof is given in appendix A. This clearly generalizes lemma 2.1, and in the context of PMs it implies theorem 4.2.

The pressure p satisfies a uniform wave equation in Ω. Under the transformation Ωω with Embedded Image and V2=FFt, p satisfies the equation for the pseudo-pressure of a PM with stiffness C and anisotropic inertia ρ,Embedded Image(4.3a)whereEmbedded Image(4.3b)and S satisfiesEmbedded Image(4.3c)

Note that the stress tensor S is not uniquely defined, although it must satisfy the equilibrium condition (4.3c). The associated density depends only on the left stretch tensor of F, viz. V. The IC corresponds to the special case of S=I, which is a trivial solution of equation (4.3c). The importance of theorem 4.2 is that the cloaks may simultaneously comprise PM stiffness and anisotropic inertia, which provides a vastly richer potential set of material parameters, not limited to the model of equation (2.6).

Theorem 4.2 implies that the phase speed, wave velocity vector and polarization (not normalized) for plane waves with phase direction n are, from equations (3.13) and (3.14),Embedded Image(4.4)The phase speed and wave velocity are independent of whether the cloak is an IC or the generalized PM-IC. These important wave properties are functions of the deformation only. They can be expressed in revealing forms using the deformation gradient as v=|Ftn| and c=FN, where N=Ftn/|Ftn|. Note that the polarization q does in general depend upon the PM properties through the stress S.

(i) Continuity between the cloak and the acoustic fluid

Continuity conditions at the cloak outer surface in the physical description follow in the same manner as (2.13). The main difference is that the stress in the cloak is not isotropic, and therefore the condition that the shear tractions on the boundary vanish must be explicitly stated. The conditions for the pseudo-pressure which satisfies equation (3.11) areEmbedded Image(4.5)These follow from equations (3.8) and (3.7).

(ii) Rays in the cloak are straight lines in the undeformed space

Although theorem 4.2 implies that the simple wave equation (2.7) in Ω is exactly mapped to equation (3.11) in ω and hence all wave motion properties transform accordingly, including rays, it is instructive to deduce the ray transformation separately. We now demonstrate explicitly that rays in the cloak ω, which are curves that minimize travel time, are just straight lines in Ω. Consider the straight line X(t)=X0+τN, where N is a unit vector in Ω. The associated curve in ω is x(τ)=x(X0+τN). Differentiation yields dx/dτ=FN, which is the same as V2s, where the vector sFtN. Differentiating s(τ), keeping in mind that N is fixed, givesEmbedded Image(4.6)where the compatibility identity Embedded Image has been used. We therefore deduce that straight lines in Ω are mapped to solutions of the coupled ordinary differential equations,Embedded Image(4.7)But these are identically the ray equations in the cloak (see appendix B). They are also the geodesic equations for the metric V−2. The ray equations conserve the quantity s.V2s that is equal to unity, reflecting the fact that s is the slowness vector, s=n/v (see equations (4.4) and (B 4)). An illustration of rays inside the physical cloak is presented in §5.

(iii) Relation to the Milton, Briane and Willis transformations

Milton et al. (2006) examined how the elastodynamic equations transform under general curvilinear transformations. They showed, in particular, that if the deformation is harmonic then the constitutive relation (2.4) and momentum balance (2.5) for a compressible inviscid fluid with isotropic density transform into the equations for a PM with anisotropic inertia, equations (3.6) and (3.7), respectively. The deformation is harmonic if Embedded Image, which realistically limits the transformation to the identity (Milton et al. 2006). This would appear to indicate that acoustic cloaking using the transformation method is impossible, in contradiction to the present result. In fact, as we show next, the Milton, Briane and Willis (MBW) result is a special case of the more general theory embodied in theorem 4.2, one that corresponds to the choice S=J−1V2.

The PM stiffness and inertia tensor found by Milton et al. (2006) are C=J−1V2V2 and ρ=J−1V2 (their eqns (2.12) and (2.13)). These are of the general form required by equation (4.3b) if we identify S as S=J−1V2. Does this satisfy the equilibrium condition (4.3c)? Using equation (2.3) Embedded Image and this vanishes if the deformation is harmonic. The MBW transformation therefore falls under the requirements of theorem 4.2 for the specific choice of S=J−1V2 that satisfies the equilibrium equation (4.3c) only if the deformation is harmonic.

Having shown that the MBW transformation result is a special case of the present theory, it is clear that the transformation as considered here is different from theirs. Milton et al. (2006) demand that all of the equations transform isomorphically, whereas the present theory requires only that the scalar acoustic wave equation is mapped to the scalar wave equation for the PM (see equations (4.3a)). The mapping contains an arbitrary but divergence-free tensor S that defines the particular but non-unique constitutive relation (2.4) and momentum balance (2.5). Consider, for instance the displacement fields u(X) and u in Ω and ω, respectively. Under the transformation of (Milton et al. 2006) Embedded Image (eqn (2.2) of Milton et al. 2006). There is no analogous constraint in the present theory. In other words, we do not require an isomorphism between the equations for all of the field variables. Instead, the scalar wave equation for the acoustic pressure is isomorphic to the scalar equation for the pseudo-pressure of the PM.

(b) Cloaks with isotropic inertia

Theorem 4.2 opens up a vast range of potential material properties. It means that there is no unique cloak associated with a given transformation Ωω and its deformation gradient F. We now take advantage of this non-uniqueness to consider the possibility of isotropic inertia. Equation (4.3b) indicates that the density is isotropic if S is proportional to V. Hence, we deduce lemma 4.3.

A necessary and sufficient condition that the density is isotropic, ρ=ρI, is that there is a scalar function h(x), such thatEmbedded Image(4.8)in which case,Embedded Image(4.9)and the Laplacian is Embedded Image.

There is a general circumstance for which a solution can be found for h. It takes advantage of the second-order differential equality,Embedded Image(4.10)Although F is generally unsymmetric, F=Ft in the special case that the deformation gradient is a pure stretch with no rotation (R=I). We therefore surmise lemma 4.4.

If the deformation gradient is a pure stretch (R=I and hence F coincides with V) then the density is isotropic,Embedded Image(4.11)and the Laplacian becomes Embedded Image.

The infinite mass problem of the IC can be avoided if the material near the inner boundary ∂ω has integrable mass. This could be achieved, for instance, by requiring that the deformation near ∂ω is symmetric (pure stretch). Lemma 4.4 and the scaling arguments of §2d imply that the isotropic density scales as ρ=O(ϵdβ), which is integrable as long as β<d+1 (α>1/(d+1)).

(c) Example: the rotationally symmetric cloak

We again consider the deformation of equation (2.14) for the cloak Embedded Image and assume that the symmetric tensor S has the form Embedded Image. Differentiation yields Embedded Image, and the ‘equilibrium’ condition (4.3c) is satisfied if w(r) and γ(r) are related by Embedded Image. It is convenient to introduce a new function g(r), such that γ=rg′/g and w=(g/r)d−1, which automatically makes Embedded Image. The cloak parameters, therefore, have general rotationally symmetric formsEmbedded Image(4.12)

The functions f and g are independent of one another, and together define a two-degree of freedom class of PM-IC model. The general solution has both anisotropic stiffness and anisotropic inertia. The previous example of the pure IC corresponds to the special case of g=r, for which equation (4.12) gives S=I and K, and ρ agree with equation (2.16).

The form of the stress S indicates the PM-IC has TI symmetry. This is a special case of the orthotropic PM considered earlier. A normal TI solid with axis of symmetry in the x3-direction has five independent elastic moduli: C11, C33, C12, C13 and C44. The last is a shear modulus, the other shear modulus is C66=(C11C12)/2. We set all shear moduli to zero, implying C44=0 and C12=C11, and the remaining independent moduli C11, C33 and C13 satisfy Embedded Image. The PM, therefore, has two independent elastic moduli. Let C33Kr(r), C11K(r) and Embedded Image, then the fourth-order elasticity tensor defined by (4.12) isEmbedded Image(4.13)where the stiffnesses Kr and K, and the principal values of the inertia tensor given by equation (2.17), areEmbedded Image(4.14)The phase speeds cr and c in the principal directions are again given by equation (2.19). This might seem amazing at first sight, but recall that it is predicted from the general theory. That is, the phase speed and wave velocity are independent of how we interpret the cloak material, as an IC or the more general PM-IC. In this example, it means that the phase speed and wave velocity are independent of g.

(i) Pure PM cloak with isotropic density

The inertia is isotropic when ρr=ρ, which occurs if g(r)=f(r). In that case ρ=ρI, and equation (4.14) reduces toEmbedded Image(4.15)We observe that the parameters of equation (4.15) are obtained from the IC parameters in equations (2.16) and (2.17) under the substitutions Embedded Image. Thus, the universal relation analogous to equation (2.18) is nowEmbedded Image(4.16)and by analogy with equation (2.20) the three original material parameters can be expressed using the phases speeds only, asEmbedded Image(4.17)In summary, there is a one-to-one correspondence between the two sets of three material parameters for the limiting cases of the pure IC on the one hand, and the pure PM cloak on the other. Of course, as discussed before, the density and stiffness cannot be simultaneously isotropic. The PM-IC model with material properties (4.12) includes both limiting cases when g=r and f, respectively.

Table 2 summarizes the scaling of the physical quantities for isotropic inertia, similar to the scalings in table 1 for the pure IC. Note that the wave speeds cr, c and the intermediate (C13) modulus Embedded Image have limiting behaviour that is independent of the dimensionality, while the density ρ and the moduli Kr and K depend upon whether the cloak is in two or three dimensions.

View this table:
Table 2

Behaviour of quantities near the vanishing point r=a for the scaling fξα as ξ=ra↓0 with isotropic inertia.

5. Further examples

(a) A non-radially symmetric cloak with finite mass

The examples considered above are rotationally symmetric and rather special in that they can be made using uniformly pure IC, or pure PM, or hybrid PM-IC. The pure IC model is always achievable as lemma 2.1 showed, but it suffers from the infinite mass catastrophe. The pure PM model requires that lemma 4.3 hold at all points, which is not realistic. However, we can always obtain a cloak comprising partly pure PM by requiring the deformation to be locally a pure stretch (lemma 4.4). In particular, by constraining the deformation near the inner surface ∂ω in this manner, the density can be made both isotropic and integrable. We now demonstrate this for a non-rotationally symmetric cloak.

For A∈Sym+, h(ζ), h′(ζ)>0 for ζ∈[0,1], consider the deformation,Embedded Image(5.1)This generalizes the deformation of equation (2.14) (A=I) and has the important property that the deformation gradient is symmetric,Embedded Image(5.2)The inner surface is an ellipse (two-dimensional) or ellipsoid (three-dimensional),Embedded Image(5.3)

The mapping Embedded Image must be the identity on the outer surface of the cloak ∂ω+=∂Ω. This eliminates the transformation (5.1) as a possible deformation in the vicinity of ∂ω+ but it does not rule it out elsewhere. In particular, it can be used on the inner surface ∂ω and for a finite surrounding volume. Then it could be patched to a different mapping closer to the outer boundary of the cloak, one which reduces to the identity on ∂ω+. For instance,Embedded Image(5.4)where ν(x)=1 for all x between ∂ω and some surface Embedded Image, beyond which ν decreases smoothly to zero as x approaches ∂ω+, which is assumed to be a level surface of ζ, i.e. an ellipsoid or an ellipse. We assume that ζ=1 on the outer surface, so thatEmbedded Image(5.5)Let Embedded Image be the level surface ζ0 for constant ζ0∈(0,1). The surface separating the pure PM inner region from the PM-IC outer part of the cloak is thereforeEmbedded Image(5.6)Based on lemma 4.4, the inner part of the cloak between ∂ω and Embedded Image can be constructed from pure PM with isotropic density Embedded Image. The remaining part of the cloak is PM-IC and the mass of the entire cloak will be finite.

For instance, in figure 2, h(ζ)=(1/2)(1+ζ) for ζ∈[0,1], ν=1 for Embedded Image and ν=4(1−ζ) for Embedded Image, and the principal values of A are 0.6 and 1.0. It also shows each ray following a continuous path through the cloak with collinear incident and emergent ray paths. There is a unique ray separating the rays traversing the cloak in opposite senses, and which defines a ‘stagnation point’ at the cloak inner surface. The separation ray is the one that would intersect the singular point in the undeformed space, O in figure 1. This is the origin in figure 2 and since the rays are incident horizontally, the separation ray is defined by x2=0 outside the cloak, and it intersects ∂ω at Embedded Image. The wavefront in effect splits or tears apart at the incident intersect and it reforms at the emergent intersect. The time delay between these two events is infinitesimal since the tearing/rejoining is associated with the instant at which the wavefront would traverse O in the undeformed space. A time-lapse movie illustrating this more vividly may be seen in the electronic supplementary material (20 s long). Another movie showing the ray paths for different directions of incidence can be found in the electronic supplementary material.

Figure 2

Ray paths through a non-radially symmetric cloak. The solid curves are the inner and outer surfaces of the cloak. The dashed curve delineates the inner region in which the deformation gradient is symmetric everywhere and the cloak is pure PM with finite isotropic mass. Two movies of the rays and the wavefronts in this cloak may be viewed in the electronic supplementary material.

(b) Scattering from near-cloaks

A near-cloak or almost perfect cloak is defined here as one with inner surface ∂ω that does not correspond to the single point X=O. We illustrate the issue using the radially symmetric deformation (2.14) with f(a) small but non-zero, and assuming time harmonic motion, with the factor e−ikt understood but omitted. Since the inner surface is not the image of a point, it is necessary to prescribe a boundary condition on the interior surface, which we take as zero pressure on r=a. The specific nature of the boundary condition should be irrelevant as f(a) shrinks to zero.

As before, the cloak occupies Embedded Image, but now f(a)>0. The total response for plane wave incidence is Embedded Image,Embedded Imagein two and three dimensions, respectively, and p0 is a constant. A near-cloak can be defined in many ways: for instance, a power law f(r)=b((rδ)/(bδ))α with 0<δ<a is considered in Norris (2008). Here we assume a linear near-cloak mapping similar to the one examined by Kohn et al. (2008),Embedded Image(5.7)where 0<δ<a. Hence, f(δ)(a)=δ and the radius at which the mapping is zero, r=aδ(ba)/(bδ), defines the size of a smaller but perfect cloak.

Some representative results are shown in figure 3, which illustrates clearly a disparity between the cylindrical and spherical cloakings, even when the physical optics cross-sections are identical. Thus, for f(a)=0.01a, the three-dimensional cross-section is negligible (figure 3d) but the two-dimensional cross-section is two orders of magnitude larger (figure 3c). Ruan et al. (2007) found that the perfect cylindrical EM cloak is sensitive to perturbation. This sensitivity is evident from the present analysis through the dependence on the length δ that measures the departure from perfect cloaking δ=0.

Figure 3

A plane wave is incident from the left with frequency k=10 on the cloak defined by equation (5.7) with a=1 and b=21/(d−1). The outer cloak radius b is chosen so that the geometrical cross-section of the cloak is twice that of the cloaked region in both two dimensions ((a) δ=0.3, Σ=1.48 and (c) δ=0.01, Σ=0.12) and three dimensions ((b) δ=0.3, Σ=0.79 and (d) δ=0.01, Σ=1×10−3). The circular core in the plots is the cloaked region of radius a. The virtual inner radius f(a)=δ is 0.3 or 0.01, and Σ is the total scattering cross-section.

The ineffectiveness of the same cloak in two dimensions when compared with three dimensions can be understood in terms of the scattering cross-section. The leading order far-field is of the form Embedded Image. The optical theorem implies that the total scattering cross-section, and hence the total energy scattered, is determined by the forward scattering amplitude, Embedded Image. Thus,Embedded Image(5.8)The cross-section is dominated in the small kf(a) limit by the n=0 term, with leading order approximations,Embedded Image(5.9)This explains the greater efficacy in three dimensions, and suggests that all things being equal, cylindrical cloaking is more difficult to achieve than its spherical counterpart.

6. Discussion and conclusion

Starting from the idea of an acoustic cloak defined by a finite deformation, we have shown that the acoustic wave equation in the undeformed region is mapped into a variety of possible equations in the physical cloak. Theorem 4.2 implies that the general form of the wave equation in the cloak isEmbedded Image(6.1)where the stress-like symmetric tensor S is divergence free and the inertia tensor is ρ=JSV−2S. The non-unique nature of S for a given fixed deformation opens many possibilities for interpreting the cloak in terms of material properties.

If S is constant (S=I with no loss in generality) then the cloak material corresponds to an acoustic fluid with pressure p defined by a single bulk modulus but with a mass density ρ that is anisotropic, which we call the IC. The IC model is mathematically consistent but physically impossible because it requires a cloak of infinite total mass. There appears to be no way to avoid this if one restricts the cloak material properties to the IC model. If one is willing to use an imperfect cloak with finite mass, and is concerned with fixed frequency waves, then the scattering examples show that significant cloaking can be obtained by shrinking the effective visible radius to be sub-wavelength. The two- and three-dimensional responses for imperfect cloaking are quite distinct, with far better results found in three dimensions.

A cloak of finite mass is achievable by allowing S to be spatially varying and divergence free. The general material associated with equation (6.1), called PM-IC, has both anisotropic inertia and anisotropic elastic properties. The elastic stiffness tensor has the form of a PM characterized by the symmetric tensor S and a single modulus K. Under certain circumstances, characterized in lemmas 4.3 and 4.4, the density becomes isotropic and the material is pure pentamode. More importantly, the total mass can be made finite.

The finite mass problem arises from how we interpret the cloak material in the neighbourhood of its inner surface. It is therefore not necessary to totally abandon the pure IC model, but it does mean that the alternative PM-IC is required at the inner surface. From the examples considered here it appears that one can always use a pure PM model near the inner cloak surface, and thereby achieve finite mass. One method is to force the deformation near the inner surface to be a pure stretch, then lemma 4.4 implies that the density is locally Embedded Image. The total mass remains finite as long as ρ is locally integrable, which is easily achieved.

The theory and simulations of PM-IC and PMs presented here illustrate the wealth of possible material properties that are opened up through the general PM-IC model of acoustic cloaking. The physical implementation is in principle feasible: for instance, anisotropic inertia can be achieved by microlayers of inviscid acoustic fluid (Schoenberg & Sen 1983), while the microstructure required for PMs has been described (Milton & Cherkaev 1995). Fabrication of practical PM-IC materials remains as a challenging but worthwhile goal.


Constructive suggestions from the anonymous reviewers are appreciated.



View Abstract