Sound trapping and waveguiding in locally resonant viscoelastic phononic crystals

We describe the trapping and absorption of audible sound in centimeter-scale claddings of two-dimensional, locally resonant phononic crystals. In a square lattice of local resonators consisting of steel cores and cellulose shells, embedded in a viscous foam, dual acoustic-range band gaps extending from about 200 to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2850 \text { Hz}$$\end{document}2850Hz are achieved. The spectral range consists of a low-frequency, local resonance gap, separated from a higher frequency Bragg resonance gap, by narrow bands of slow-sound modes. We demonstrate that thin claddings of such phononic crystal, of only three unit cells in thickness, can effectively prevent sound transmission, by a combination of reflection and absorption, over much of the audible spectrum. Moreover, frequency-selective sound transmission can be enabled by engineering waveguide channels that transmit sound through the local resonance gap, the Bragg gap, or both. This offers a path to sound-sculpting claddings that can surround a noise-generating source. The viscoelastic foam in our cladding is treated using a fractional Voigt model, capable of describing experimentally observed responses.


Viscoelasticity in the fractional Voigt model
At a given position x = (x 1 , x 2 , x 3 ) , the field u = (u 1 , u 2 , u 3 ) defines the displacement from equilibrium of the infinitesimally small parcel of elastic material.The local deformation is given by the symmetric strain tensor ε ij 18 : Linear elasticity is assumed by dropping the nonlinear term (1/2) k (∂u k /∂x i )(∂u k /∂x j ).
The energy loss in a damped harmonic oscillator arises from a frictional force that is proportional to the magnitude and opposite in direction to the velocity of the oscillator.Similarly, viscoelasticity in a linear continuum is described by a damping term of the stress that is proportional to the rate of strain 19 : where σ ij denotes the stress tensor, C ijkl is a fourth order stiffness tensor, and D ijkl is a fourth order damping tensor.The dynamical field equation in a viscoelastic medium of density ρ satisfies Newton's second law: By forming a scalar product of the dynamical equation and the velocity field, and using the symmetry of the strain tensor, we obtain a continuity equation describing the energy and energy flow densities in the material: where u KE is the kinetic energy density, u PE is the elastic potential energy density, S is the energy flux, and p loss is the power loss density due to viscoelasticity: In oscillatory systems, it is customary to represent viscoelasticity by complex-valued, effective elastic moduli.At a driving frequency ω , all components oscillate with the time dependence of a complex exponential e −iωt in the steady state.By Eq. ( 2), the damping tensor can be represented by the imaginary part of a generalized, complexvalued stiffness tensor Cijkl (ω) ≡ C ijkl − iωD ijkl .The loss tangent η ijkl (ω) ≡ −Im( Cijkl (ω))/Re( Cijkl (ω)) denotes the ratio of the imaginary part to the real part of the complex-valued stiffness tensor.For simplicity, we assume that the viscoelastic dissipation is isotropic, and is the same for volumetric and shear deformations η ijkl = η.
The fractional Voigt model for viscoelastic material, depicted in the insert of Fig. 1, consists of an elastic spring connected in parallel to a generalized dashpot 15 .The strain ε is uniform across the components.The elastic spring of bulk modulus K produces a stress σ e = Kε that is linearly proportional to the applied strain.On the other hand, the generalized dashpot of power exponent α and characteristic time τ produces a stress σ d = K(−iωτ ) α ε .For components connected in parallel, the stresses are additive: where a complex-valued, effective modulus K(ω) and a loss tangent η(ω) are identified: (1) ∂u i ∂x j + ∂u j ∂x i .
(2) For positive values of α , the effective modulus reduces to the elastic modulus of the elastic spring at zero fre- quency K(ω = 0) = K .The mechanical response of the generalized dashpot (σ d = K(−iωτ ) α ε) is purely elastic when α = 0 , and purely viscous when α = 1 .The generalized fractional dashpot can be interpreted as an analytic interpolation between the purely elastic and purely viscous limits.The strain-stress relation of the generalized dashpot can be formally expressed in terms of Riemann-Liouville fractional calculus 20,21 , and is experimentally shown to accurately characterize soft materials 15,22,23 .For simplicity, we employ the simple fractional Voigt model in this article.More realistic models, such as the fractional Zener model, generally require more springs and generalized dashpots of different power exponents 24,25 .
For a weakly damped foam, we choose α = 0.05 , τ = 1 s , K(ω = 0) = 15,000 Pa and mass den- sity ρ f = 25 kg m −3 .The wavelength f = 2π/(ωRe[(ρ f / K(ω)) 1/2 ]) and absorption length scale ) are plotted as a function of frequency in Fig. 1.For a phononic crystal structure to exhibit wave-functional properties, the absorption length is at least an order of magnitude longer than the elastic wavelength.This weakly damped regime is referred to as the metamaterial regime in the literature 26 .In this paper, we consider only the weakly damped regime as depicted in Fig. 1.In contrast, when the absorption length is comparable to or shorter than the elastic wavelength, the oscillation is damped in the viscoelastic medium before wave interference has any observable effects.The strongly damped regime is referred to as the boundary layer material regime 26 .

Locally resonant viscoelastic phononic crystals
Phononic band gaps are observed in locally resonant phononic crystals.The absence of wave modes in a phononic band gap implies that wave propagation is evanescent.The band gaps can be broadly classified into local resonance gaps and Bragg gaps.Local resonance gaps occur over the frequency range where the effective masses and moments of inertia of the acoustically connected local resonators are simultaneously negative.The frequency range of local resonance gaps is insensitive to lattice arrangements, provided that the resonators remain acoustically connected.Bragg gaps are caused by the macroscopic resonances of the array of local resonators, and sensitively depend on the lattice configuration.Frequency-sensitive waveguide modes can be engineered in a desired spectral range, by line defects of local resonance and crystal periodicity.
The supercell method is employed to study point-defect cavity modes and line-defect waveguide modes in periodic systems.In the present work, we study cavity modes and one-dimensional waveguides in two-dimensional locally resonant phononic crystals.Cylindrically symmetric core-shell resonators are arranged in a square lattice with a lattice constant a = 1 cm .Each local resonator is composed of a dense steel core, coupled to a stiff cellulose shell, via a layer of soft, absorbing foam.The core, interstitial foam and shell constitute 20% , 20% and 10% by volume respectively, so that the radius of the core R 1 = a √ 0.2/π ≈ 0.2523 cm , the inner radius of the cellulose shell R 2 = a √ 0.4/π ≈ 0.3568 cm and the outer radius R 3 = a √ 0.5/π ≈ 0.3990 cm .The background and the interstitial foam are assumed to be the same type of isotropic, fractional Voigt viscoelastic material, with power exponent α = 0.05 , characteristic time τ = 1 s , and Lamé constants at zero frequency ˜ (ω = 0) = 7500 Pa and μ(ω = 0) = 3750 Pa .The mass density and Lamé parameters of steel (cellulose) are ρ c = 7940 kg m −3 , c = 107.5 GPa and µ c = 78.15GPa 27 ( ρ s = 1350 kg m −3 , s = 1.21 GPa and µ s = 0.519 GPa 28 ).The unit cell is depicted in Fig. 2.
With a small power exponent, the effective elastic moduli are fairly insensitive to frequency at frequencies above 1/τ .From 50 Hz to 3000 Hz , the real parts of the effective Lamé parameters increase mod- estly by approximately 13% from Re( ˜ )(f = 50 Hz) ≈ 17500 Pa to Re( ˜ )(f = 3000 Hz) ≈ 19700 Pa , and from Re( μ)(f = 50 Hz) = 8730 Pa to Re( μ)(f = 3000 Hz) = 9870 Pa .For simplicity, we choose the val- ues at 600 Hz to calculate the supercell band structures and transmission spectra with lossless foam: Re( ˜ )(f = 600 Hz) ≈ 18800 Pa and Re( μ)(f = 600 Hz) = 9390 Pa .These particular choices of elastic con- stants with zero imaginary parts are referred to as lossless foam in this article.The actual band structures and www.nature.com/scientificreports/transmission peaks are red-shifted at frequencies below 600 Hz and blue-shifted above 600 Hz .As frequency scales inversely with the square root of the elastic moduli, the simplifying assumption introduces about 3.7% error at 50 Hz and 2.4% error at 3000 Hz.
The actual band structure with frequency-dependent elastic moduli can be calculated by the Cutting Surface Method 16,29 , where a generalized dispersion ω(K, p , µ p ) is evaluated parametrically as a function of a hypothetical elastic moduli p and µ p .The resultant dispersion hypersurfaces are cut by the condition ( p , µ p ) = (Re( ˜ )(ω), Re( μ)(ω)).
The in-plane bands of the 1 × 9 supercell are plotted in Fig. 3, complemented with the displacement field profile of representative eigenmodes at the edge of the Brillouin zone M .Each supercell contains nine resona- tors.Each resonator has a rigid core and a rigid shell, individually providing two in-plane translational degrees of freedom and one in-plane rotational degree of freedom.Therefore, the low-frequency response of the 1 × 9 supercell is predominantly recaptured by the 9 × 2 × (2 + 1) = 54 degrees of freedom.The first 27 bands from 0 to 217 Hz are associated with the in-phase translational and rotational oscillations between the cores and shells of individual resonators.These are illustrated by the in-phase translational oscillations at 122 Hz and in-phase rotational oscillations at 200 Hz at M .The 28th to 54th bands from 686 to 1200 Hz are associated with their anti-phase oscillations.At M , the eigenmodes associated with anti-phase rotational and translational oscillations occur at 686 Hz and 1007 Hz respectively.The in-phase and anti-phase bands are separated by a local resonance gap from 217 to 686 Hz .Above the anti-phase bands is a Bragg gap from 1200 to 2850 Hz .The 55th band occurs above the Bragg gap when the wavelength in the background medium becomes comparable to the lattice spacing.A representative eigenmode at 2854 Hz is drawn.The band structure of the supercell is verified to be consistent with that of a square lattice, with each unit cell containing a single resonator, by band folding.
The finite-element method (FEM) is employed to simulate the reflection, transmission and absorption across a finite-thickness locally resonant phononic crystal.The two-dimensional system is assumed to be infinite in the out-of-plane direction, and repeating in the lateral direction.We consider a cladding that is 3 unit cells thick.The phononic crystal cladding unit cell is chosen to be a 3 × 9 array of cylindrically symmetric core-shell resonators embedded in foam.We use this large unit cell as our reference for the introduction of defects in subsequent analysis.The phononic crystal slab is sandwiched between two layers of air.At the air-solid interfaces, the boundary conditions are specified by the continuities of normal stress, tangential stress and normal displacement.At solid-solid interfaces within the phononic crystal, both the normal and tangential components of stress and displacement fields are continuous.The two ends of the simulation domain are terminated separately by a perfectly matched layer (PML) with a thickness of several wavelengths.For normal incidence, periodic boundary conditions are imposed on the wavefield on the lateral sides.The fine mesh has at least 21 grid points per wavelength to ensure numerical accuracy.
The transmission, reflection and absorption coefficients across the 3-cm layer of locally resonant acoustic phononic crystal are plotted in Fig. 3.The results are verified to be consistent with the simpler alternative use of a 3 × 1 array of resonators as the unit cell.The coupling between air and the finite-thickness phononic crystal depends sensitively on the frequency of the incident field, and corresponds to the supercell band diagram in Figure 2. A unit cell of the phononic crystal consists of a dense steel rod of radius R 1 ≈ 0.2523 cm , coupled to a concentric, circular, annular cellulose shell of inner radius R 2 ≈ 0.3568 cm and outer radius R 3 ≈ 0.3990 cm .The interstitial space and the background are filled with viscoelastic foam.The resonators are arranged in a square lattice of lattice constant a = 1 cm.Fig. 3.The reflection troughs, transmission peaks and absorption peaks below 200 Hz are associated with the acoustic bands arising from the in-phase oscillations between the core and the shell.The reflection plateau, flat transmission minimum and absorption minimum from 200 to 800 Hz are related to the local resonance gap.Between 800 to 1200 Hz , there are appreciable reflection troughs, and transmission and absorption peaks, which are associated with the anti-phase bands.Over a broad frequency range of the Bragg gap from 1200 to 2400 Hz , waves are predominately reflected by the three-unit-cell-thick phononic crystal cladding.

Sound trapping by frequency-selective cavity modes
The elastic displacement field is evanescent in the band gap of the phononic crystal.Point defects in the otherwise periodic lattice create spatially localized cavity modes, and enable frequency-selective coupling of sound into the band gaps of the original phononic crystal 30,31 .When a resonator is replaced by the background material, cavity modes are created exclusively in the Bragg gap of the original structure.In contrast, removal and subsequent replacement of a resonator by a defect shell creates cavity modes exclusively in the local resonance gap.In the following, we delineate the two types of cavity modes and illustrate their interactions with an incident pressure wave in FEM simulations.
Cavity modes are introduced in the Bragg gap by the removal of a resonator in the otherwise perfect lattice.We refer to this kind of positional disorder as a Bragg point defect.The transmission spectrum across such a locally resonant phononic crystal cladding, of 3 unit cells in thickness, is plotted in Fig. 4. Compared to the base case in Fig. 3, there are additional transmission peaks at 1517 Hz and 2462 Hz .These are associated with spatially localized dipolar and monopolar deformations in the background material, respectively, depicted in the lower left panel in Fig. 4. In contrast, the broad transmission minimum over the local resonance band gap from 200 to 700 Hz is unaffected by the Bragg point defect.
Cavity modes are introduced exclusively in the local resonance gap by the replacement of a core-shell resonator by a defect shell.In this case, the steel core of the resonator is replaced by foam.It is referred to as a local resonance point defect.The transmission spectrum across a 3-cm locally resonant phononic crystal slab with a local resonance point defect is plotted in Fig. 4. Compared to the base case in Fig. 3, there are additional transmission peaks at 551 Hz and 792 Hz .The elastic displacement fields are depicted in the lower right panel in Fig. 4. The former transmission peak at 551 Hz is associated with the longitudinal translational oscillation of the defect shell.The disturbance is localized near the defect, because wave propagation is evanescent within the local resonance band gap.On the other hand, the latter transmission peak at 792 Hz lies in the frequency range of the anti-phase oscillatory pass bands.The longitudinal translational oscillation of the defect shell is coupled to the anti-phase rotational oscillation of the nearby core-shell resonators.In the original structure, these rotational modes do not couple to the incident pressure field by reflectional and discrete translational symmetries.The defect shell lifts the discrete translational symmetry, enabling excitations of the rotational modes via the translational oscillation of the defect shell.Defect modes can couple effectively to degrees of freedom that are previously inaccessible, unlocking rich and diverse physical responses.

Frequency-selective waveguides through dual phononic band gaps
Wave propagation is evanescent within the band gaps of a phononic crystal.Frequency-selective transmission of sound can be enabled by engineering specific waveguide channels that allow passage of sound through the otherwise impenetrable phononic band gaps 30,[32][33][34][35][36][37] .Different waveguide designs facilitate sound transmission through parts of the local resonance gap, through parts of the Bragg gap, or both gaps.These involve removal and replacement of local resonators along lines connecting the interior and exterior of the phononic crystal cladding.In the following, we describe the properties of the different types of waveguide modes in viscoelastic, locally resonant, phononic crystal slabs.We also numerically simulate their interactions with an incident pressure wave.

Bragg waveguides.
Waveguide modes are introduced in the Bragg gap by the removal of a line of resonators and re-adjustments of the distance between the phononic crystals on either side of the vacancy.We denote by d wg the additional distance between the adjacent resonators across the line defect compared to the original lattice.When a line of resonators is removed, d wg = a and the distance between the centers of the adjacent reso- nators across the waveguide is a + d wg = 2a.
The in-plane bands of a supercell with 8 remaining local resonators and a line defect of width d wg = a is plotted in Fig. 5. Compared to the band structure of the base structure in Fig. 3, numerous bands appear in the original Bragg gap from 1200 to 2850 Hz .These waveguide modes occur at practically the same frequencies using larger supercells, confirming the convergence of our supercell method.Similar supercell calculations are repeated with smaller values of d wg .The supercell band diagrams for d wg = 0.6a and d wg = 0.2a are plotted in Figs. 5 and 6 respectively.A striking common feature is the presence of bands in the Bragg gap of the original structure.The number of waveguide modes and their range of frequencies are reduced as d wg is reduced.On the other hand, the local resonance gap remains unaffected by such line defects.
The waveguide mode elastic displacement fields, for d wg = 0.2a are plotted in Fig. 6, for sound propagation at the center of the irreducible Brillouin zone Ŵ .Arranged in the order of increasing frequency, these additional waveguide modes are associated with the fundamental longitudinal oscillation at 1949 Hz , transverse oscillation at 2574 Hz and rotational oscillation at 2826 Hz , within the line defect.With a wider line defect, the fundamental modes occur at lower frequencies, while more overtone modes appear within the original Bragg gap.Plotted in Fig. 5 are overtone transverse modes at 1924 Hz and 2416 Hz with one and two nodes, respectively, along the center of the waveguide for d wg = 0.6a .In general, with increasing waveguide width, higher harmonics appear and the original Bragg gap is completely filled with propagating sound modes when d wg ≥ 0.6a.
The transmission, reflection and absorption coefficients across the 3-cm cladding with a line defect d wg = 0.2a are plotted in Fig. 6.Consistent with the base case in Fig. 3, reflection troughs, transmission peaks and absorption peaks appear over the frequency ranges of the in-phase and anti-phase pass bands.Moreover, the reflection plateau, flat transmision minimum and absorption minimum from 200 to 800 Hz correspond to the local resonance gap, which is insensitive to small amounts of positional disorder.On the other hand, in the case of lossless foam, there are appreciable transmission peaks from 1900 to 2600 Hz , within the Bragg gap.Very small remnants of these transmission peaks appear with absorbing foam.With the viscoelastic foam, reflection troughs and absorption peaks coincide with these transmission peaks.These indicate the frequency-selective coupling of the incident acoustic wave with the waveguide modes of the cladding.
The coupling between the incident pressure wave and the phononic crystal cladding is restricted by symmetry requirements and boundary conditions.By reflectional symmetry, the pressure wave at normal incidence can excite the longitudinal waveguide modes, but not the transverse and rotational modes.Moreover, shear and compressional waves propagate at different speeds of sound, and are generally decoupled in linear, isotropic solids.The compressional pressure wave in air does not couple to the transverse and rotational shear waves in the elastic or viscoelastic solids.
Local resonance defects and waveguides.Waveguide modes are introduced exclusively in the local resonance gap by a line defect where the local resonators are replaced by non-resonant shells.The supercell band structure with 8 local resonators and a line defect of non-resonant shells is plotted in Fig. 7. Compared to the supercell band structure of the base structure in Fig. 3, three additional bands appear in the local resonance gap.On the other hand, the Bragg gap is unaffected by this line defect.This is because the non-resonant shells occupy the same regions as the local resonators in the original structure.At high frequencies relative to the resonant and normal mode frequencies of the local resonator, the effective mass and effective moment of inertia of the core-shell resonator approach the mass and moment of inertia of the shell alone.The relevant waveguide mode displacement fields are plotted in Fig. 7, at the center of the irreducible Brillouin zone Ŵ .In the order of increas- ing frequency, these eigenmodes are associated with the longitudinal oscillation at 307 Hz , transverse oscillation at 522 Hz and rotational oscillation at 632 Hz , of the defect shell.
The transmission, reflection and absorption coefficients across the 3-cm layer of locally resonant phononic crystal are plotted in Fig. 7. Compared to the base case in Fig. 3, a striking feature is the emergence of reflection troughs, transmission peaks and absorption peaks from 200 to 800 Hz , within the frequency range of the local resonance gap.At normal incidence, the pressure wave in air does not couple to the waveguide modes involving the transverse and rotational oscillations of the defects.
Waveguide with a rigid rod.Waveguide modes extend throughout both the local resonance gap and the Bragg gap by the removal of a line of resonators and concomitant addition of a rigid cellulose rod in the waveguide foam.The supercell band structure with 8 local resonators and a line defect with a rigid rod waveguide is plotted in Fig. 8.The cellulose rod, of width 0.1a, is made of the same material as the shell of the resonator.Compared to the supercell band structure of the base structure in Fig. 3, additional bands appear in both the local resonance gap and the Bragg gap.The removal of a line of resonators results in propagating sound modes in the Bragg gap of the original phononic crystal.At the same time, the rigid cellulose rod introduces translational oscillation modes within the local resonance gap.
The displacement fields of these relevant waveguide modes are plotted in Fig. 8 at the center Ŵ of the irre- ducible Brillouin zone.Within the original local resonance gap are two waveguide modes associated with the longitudinal, in-plane translational oscillation of the rigid cellulose rod at 241 Hz and its transverse, in-plane translational oscillation at 427 Hz .In-plane rotational oscillation of the rod is inhibited, because the phononic crystal is assumed to be infinite in the longitudinal direction and the relevant planar moment of inertia is theoretically infinite.On the other hand, the Bragg gaps are occupied by odd and even symmetry oscillations in the background foam, on either side of the rigid rod.In the order of increasing frequency, these are the odd-parity longitudinal oscillation at 1629 Hz , even-parity longitudinal oscillation at 1662 Hz , odd-parity transverse oscil- lation at 2664 Hz , even-parity transverse oscillation at 2676 Hz , even-parity rotational oscillation at 2799 Hz and odd-parity rotational oscillation at 2804 Hz.
The transmission, reflection and absorption coefficients across the 3-cm layer of locally resonant phononic crystal are plotted in Fig. 8. Compared to the base case in Fig. 3, there are additional reflection dips and absorption peaks in both the local resonance gap and the Bragg gap.Not all waveguide modes, described above, can  Arranged in the order of increasing frequency from bottom to top, the waveguide modes are associated with the longitudinal oscillation of the solid waveguide at 241 Hz , its transverse oscillation at 427 Hz , odd- parity longitudinal oscillation of the background foam in the waveguide at 1629 Hz , even-parity longitudinal oscillation at 1662 Hz , odd-parity transverse oscillation at 2664 Hz , even-parity transverse oscillation at 2676 Hz , even-parity rotational oscillation at 2799 Hz and odd-parity rotational oscillation at 2804 Hz .(c) The transmission T spectrum across a locally resonant, lossless, phononic crystal slab with a rigid rod waveguide is plotted.Compared to the base case, the additional transmission peaks centered at 240 Hz , 1600 Hz , 1900 Hz and 2300 Hz arise from the longitudinal waveguide modes.(d) Similarly, the reflection R , transmission T and absorption A coefficients across the same geometry with the viscoelastic foam are plotted as a function of frequency.Reflection dips and absorption peaks occur at frequencies comparable to the transmission peaks in the case with the lossless foam.These observations agree with the supercell band diagram and waveguide modes.be excited by an incident pressure wave due to symmetry considerations.The decoupling of shear and compressional waves in isotropic solids inhibits the coupling to the transverse and rotational waveguide modes.At normal incidence, reflection symmetry restricts the excited displacement field in the phononic crystal cladding to be of even parity.
In the above illustrations, overall transmission levels remain quite low due to our use of a relatively lossy background foam.A less lossy, environmentally friendly, natural fiber-reinforced composites 38 may function quite effectively within locally resonant phononic crystal claddings to selectively allow more transmission of desired sound frequencies.

Discussion
In summary, we have illustrated the sound trapping and waveguiding properties of various locally resonant phononic crystals.This reveals remarkably broadband control of reflection, transmission and absorption of audible sound in phononic crystal claddings of only a few unit cells in thickness.This is facilitated by the dual band gap nature of the underlying locally resonant medium.Our work reveals important design principles in the construction of acoustic claddings for the mitigation of unwanted noise and the spectral sculpting of sound for the surrounding environment.
Locally resonant phononic crystals exhibit a local resonance gap and a distinct Bragg gap.The former occurs when the effective mass and effective moment of inertia are simultaneously negative and relies on the anti-phase oscillations of acoustically connected resonators.The latter, resulting from the macroscopic resonances of a periodic collection of resonators, scales with the lattice spacing.Spatially localized cavity modes are created in the local resonance gap by replacing a core-shell resonator by a defect shell.Likewise, waveguides in the local resonance gap can be engineered by replacing a line of local resonators by non-resonant shells.The resulting cavity and waveguide modes involve longitudinal, transverse and rotational oscillations of the non-resonant shells.On the other hand, cavity modes appear in the Bragg gap when a core-shell resonator is replaced by the background foam.The frequency extent of the waveguide modes can be controlled by the width of the positional line defect.For a narrow waveguide with d wg = 0.2a , the waveguide states fill the upper half of the Bragg gap.These cor- respond to the fundamental longitudinal, transverse and rotational oscillations within the waveguide foam.For a wider waveguide with d wg ≥ 0.6a , the higher harmonics of foam distortion fill the entirety of the Bragg gap.
The acoustic response to an incident pressure wave of the locally resonant phononic crystal that we have considered is predominated by the translational and rotational oscillations of the rigid components.However, the rotational degrees of freedom do not couple effectively when the pressure wave is normally incident.Defects in the phononic crystal lift the discrete translational symmetry of the lattice, enabling richer and more diverse acoustic responses.
Reflection symmetry can be lifted through oblique incidence of the pressure wave.The FEM simulations with the various phononic crystal slabs were repeated at oblique angles of incidence.While the resultant transmission and reflection spectra show small frequency shifts of the transmission peaks, excitations of the rotational modes are insignificant.A plausible explanation is the stark contrast in the speed of sound between air and the phononic crystal constituents.The speed of sound in air is c air ≈ 343 m s −1 .The speed of longitudinal disturbance in foam is c f = 1/Re[(ρ f / K) 1/2 ] ≈ 40 m s −1 at 600 Hz .The order of magnitude contrast in the speeds of sound implies that the wave is strongly refracted towards the normal direction, when entering the phononic crystal from air.
In our present work, all numerical simulations are performed in the frequency domain.In principle, considerable computational resources can be saved using time-domain methods.When the finite-thickness structure is excited by a Gaussian wave pulse with a finite spectral width, a broad frequency range of acoustic response can be simulated simultaneously in one trial.The reflection and transmission coefficients over the desired spectral range can be calculated by numerically efficient fast Fourier transform (FFT) of the reflected and transmitted pulses respectively.However, the conventional implementation of finite-difference time-domain (FDTD) method is problematic for structures containing materials with large density contrasts 39 .Future work may focus on numerically stable and accurate implementation of time-domain simulations for elastic and viscoelastic composites involving large contrasts in material parameters.One possibility is to model resonators as rigid bodies with frequency-dependent, effective masses and moments of inertia 13,14,16,17 .Another possibility is employing an adaptive mesh, which is denser in foam with low speeds of sound, and coarse in stiff solids with high speeds of sound.Nevertheless, the mesh in the regions of the stiff solids cannot be too coarse to faithfully represent their geometries.
Future analysis would benefit from more realistic models of environmentally friendly viscoelastic materials.For simplicity, we assumed a simple isotropic fractional Voigt model with identical loss factors for compressional and shear deformations.More realistic models of viscoelastic materials may involve multiple dashpots of different power exponents and characteristic times.The actual loss factors for shear and volumetric deformations may differ, depending on the relative motion between the air in foam and the elastic skeleton 40,41 .Biodegradable "green" materials with less sound absorption are advocated to replace non-recyclable melamine and polyurethane for acoustic panels.These may enable enhanced sound transmission at higher audible frequencies.
An important generalization of our analysis is three-dimensional viscoelastic phononic crystals.It is anticipated that three-dimensional, locally resonant phononic crystals may exhibit an even broader bandwidth of audible sound control.For example, a body-centered cubic crystal of local resonators exhibits dual audiblerange band gaps over a frequency range from 600 to 3700 Hz 13,14 .The design of waveguide channels in a three- dimensional phononic crystal cladding may provide unprecedented, real-world, opportunities for audible sound control at will.

Figure 1 .
Figure 1.A microscopic fractional Voigt model consists of a spring of spring constant K, connected in parallel to a generalized fractional dashpot of power exponent α and characteristic time τ (see insert).The wavelength f (blue) and the absorption length scale ξ f (red) are plotted in logarithmic scale as a function of frequency.The relevant material constants of the viscoelastic foam are ρ f = 25 kg m −3 , K(ω = 0) = 15,000 Pa , α = 0.05 and τ = 1 s.

Figure 3 .
Figure 3. (a) The in-plane bands of the 1 × 9 supercell of a square lattice of cylindrically symmetric core-shell resonators is plotted along a high symmetry path in the irreducible Brillouin zone M → Ŵ → X → M .A local resonance gap spans from 217 to 686 Hz .The Bragg gap ranges from 1200 to 2850 Hz .(b) Representative eigenmodes at M are drawn, which are associated with in-phase translational oscillations at 122 Hz , in-phase rotational oscillations at 200 Hz , anti-phase rotational oscillations at 686 Hz , anti-phase translational oscillations at 1007 Hz and coherent Bragg scattering at 2854 Hz .(c) The transmission T spectrum across the lossless, locally resonant, phononic crystal of 3 unit cells in thickness are plotted.(d) Similarly, the transmission T , reflection R and absorption A coefficients across the locally resonant, viscoelastic, phononic crystal of 3 unit cells in thickness are plotted as a function of frequency.The transmission minima, reflection plateau and absorption minima from 200 to 800 Hz , and from to 1300 to 2500 Hz are consistent with the local resonance gap and Bragg gap in the band diagram.

Figure 4 .
Figure 4. (a) The transmission coefficient T across the locally resonant phononic crystal of 3 unit cells in thickness, with the Bragg point defect, is plotted as a function of frequency.Compared to the base case in Fig. 3, there are additional transmission peaks at 1517 Hz and 2462 Hz , associated with dipolar and monopolar cavity modes localized at the positional defect, respectively.(b) The transmission coefficient across the locally resonant phononic crystal, with the local resonance point defect, is plotted.The additional transmission peak at 551 Hz in the local resonance gap of the original phononic crystal, is associated with the longitudinal translational oscillation of the defect shell.On the other hand, the transmission peak at 792 Hz is caused by the coupling of the translational oscillation of the defect shell with the anti-phase rotational modes of the adjacent core shell resonators, lying in the anti-phase pass bands of the original phononic crystal.(c, d) The displacement field profiles of the relevant point defect modes are plotted.

Figure 5 .
Figure 5. (a) The in-plane band structure of the supercell with Bragg waveguide of width d wg = a is plotted along a high symmetry path in the irreducible Brillouin zone M → Ŵ → X → M .(b) Similarly, the band structure of the supercell with Bragg waveguide of width d wg = 0.6a is plotted.Multiple waveguide bands appear in the Bragg gap of the original structure.The local resonance gap remains almost unaltered by positional disorder.(c) For the latter case, the elastic displacement fields of selected eigenmodes, at the Brillouin zone center Ŵ , are plotted.The direction of the displacement field is indicated by blue arrows.With a wide waveguide, overtone modes, localized in the line defect, may exist in the Bragg gap.

Figure 6 .
Figure 6.(a) The in-plane band structure of the supercell with a Bragg waveguide of width d wg = 0.2a is plotted along a high symmetry path in the irreducible Brillouin zone M → Ŵ → X → M .Waveguide bands occupy the upper half of the Bragg gap of the original structure.The local resonance gap is insensitive to the slight positional disorder.(b) The displacement fields of the waveguide modes, within the Bragg gap, at the Brillouin zone center Ŵ , are plotted.The direction of the displacement field is indicated by blue arrows.The waveguide modes are associated with longitudinal oscillation at 1949 Hz (bottom), transverse oscillation at 2574 Hz (middle), and rotational oscillation at 2826 Hz (top).(c) The transmission T spectrum across the lossless, locally resonant, phononic crystal claddings, with a line defect d wg = 0.2a , of 3 unit cells in thickness, are plotted.The transmission peaks 1900 Hz to 2600 Hz are associated with the waveguide states within the original Bragg gap.(d) Similarly, the transmission T , reflection R and absorption A coefficients across the same geometry with the lossless foam replaced by a viscoelastic fractional Voigt material are plotted as a function of frequency.Very small remnants of the transmission peaks in the Bragg gap are still discernible.The reflection plateau and absorption minimum from 200 to 800 Hz arise from the local resonance gap which is insensitive to positional disorder.

Figure 7 .
Figure 7. (a) The in-plane band structure of the supercell with a local resonance defect waveguide is plotted along a high symmetry path in the irreducible Brillouin zone M → Ŵ → X → M .Waveguide bands occupy the local resonance gap of the original structure.The Bragg gap is unaffected by the local resonance defects.(b) The elastic displacement fields of the waveguide modes, within the local resonance gap, at the Brillouin zone center Ŵ , are plotted.The direction of the displacement field is indicated by blue arrows.The waveguide modes are associated with the longitudinal oscillation of the non-resonant shell at 307 Hz (bottom), its transverse oscillation at 522 Hz (middle), and its rotational oscillation at 632 Hz (top).(c) The transmission spectrum T across a locally resonant, phononic crystal slab with a waveguide line defect of non-resonant shells is plotted.The additional transmission peaks centered at 310 Hz and 460 Hz , within the local resonance gap of the original phononic crystal, indicate the coupling of the incident pressure wave with the local resonance waveguide modes.The broad transmission minimum from 1200 to 2800 Hz are consistent with the Bragg gap.(d) Similarly, the reflection R , transmission T and absorption A coefficients across the viscoelastic phononic crystal slab of the same geometry are plotted as a function of frequency.The reflection dips and absorption peaks from 200 to 800 Hz arise from the waveguide modes along the line defect.These observations agree with the supercell band diagram and waveguide modes.

Figure 8 .
Figure 8.(a) The in-plane band structure of the supercell with a rigid rod waveguide is plotted along a high symmetry path in the irreducible Brillouin zone M → Ŵ → X → M .Waveguide bands occupy both the local resonance gap and the Bragg gap of the original structure.(b)The displacement fields of the waveguide modes at the Brillouin zone center Ŵ , are plotted.The direction of the displacement field is indicated by blue arrows.Arranged in the order of increasing frequency from bottom to top, the waveguide modes are associated with the longitudinal oscillation of the solid waveguide at 241 Hz , its transverse oscillation at 427 Hz , odd- parity longitudinal oscillation of the background foam in the waveguide at 1629 Hz , even-parity longitudinal oscillation at 1662 Hz , odd-parity transverse oscillation at 2664 Hz , even-parity transverse oscillation at 2676 Hz , even-parity rotational oscillation at 2799 Hz and odd-parity rotational oscillation at 2804 Hz .(c) The transmission T spectrum across a locally resonant, lossless, phononic crystal slab with a rigid rod waveguide is plotted.Compared to the base case, the additional transmission peaks centered at 240 Hz , 1600 Hz , 1900 Hz and 2300 Hz arise from the longitudinal waveguide modes.(d) Similarly, the reflection R , transmission T and absorption A coefficients across the same geometry with the viscoelastic foam are plotted as a function of frequency.Reflection dips and absorption peaks occur at frequencies comparable to the transmission peaks in the case with the lossless foam.These observations agree with the supercell band diagram and waveguide modes.