About this Journal Submit a Manuscript Table of Contents
Advances in Acoustics and Vibration
Volume 2009 (2009), Article ID 317890, 7 pages
http://dx.doi.org/10.1155/2009/317890
Research Article

Designing 2D Phononic Crystal Slabs with Transmission Gaps for Solid Angle as well as Frequency Variation

Department of Underwater Research, Swedish Defence Research Agency, FOI Kista, 16490 Stockholm, Sweden

Received 29 January 2009; Accepted 15 June 2009

Academic Editor: Mohammad Tawfik

Copyright © 2009 Sven M. Ivansson. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Phononic crystals (PCs) can be used as acoustic frequency selective insulators and filters. In a two-dimensional (2D) PC, cylindrical scatterers with a common axis direction are located periodically in a host medium. In the present paper, the layer multiple-scattering (LMS) computational method for wave propagation through 2D PC slabs is formulated and implemented for general 3D incident-wave directions and polarizations. Extensions are made to slabs with cylindrical scatterers of different types within each layer. As an application, the problem is considered to design such a slab with small sound transmittance within a given frequency band and solid angle region for the direction of the incident plane wave. The design problem, with variable parameters characterizing the scatterer geometry and material, is solved by differential evolution, a global optimization algorithm for efficiently navigating parameter landscapes. The efficacy of the procedure is illustrated by comparison to a direct Monte Carlo method.

1. Introduction

Just as photonic crystals can be used to manipulate light, phononic crystals (PCs) with inclusions in a lattice with single, double, or triple periodicity can be used to manipulate sound [1]. When a sound wave of a certain frequency penetrates the PC, the energy is scattered by the inclusions. According to Bragg’s law, constructive and destructive interference appears in certain directions. It follows that transmission and/or reflection for certain frequencies can be absent (band gaps), even for all angles of incidence (absolute band gaps). The band gaps are ideal and fully developed only in space-filling PCs. In a PC slab with finite thickness, the “band-gap wavefields” are reduced significantly but they do not vanish throughout the gap. Acoustic frequency selective insulators and filters are possible applications.

Several computational methods have been adapted and developed to study wave propagation through PCs. Two such methods are the purely numerical finite-difference time domain (FDTD) method and the semianalytical layer multiple-scattering (LMS) method, which is developed from Korringa-Kohn-Rostoker theory [2]. These methods complement each other. Advantages with the LMS method are its computational speed, which makes it useful for forward modeling in connection with extensive optimization computations, and the physical insight it provides. It appears that the LMS method was first developed for the 3D case with spherical scatterers [3, 4], and recent review papers include Sigalas et al. [5] and Sainidou et al. [6].

The LMS method has also been applied to the 2D case with infinite cylindrical scatterers. It is mainly the in-plane propagation case that has been considered [711], for which one space dimension can be eliminated in the wave equations. Out-of-plane propagation has been treated by Mei et al. [12] and, for band structure calculations, somewhat earlier by Wilm et al. [13] (using the plane-wave method) and by Guenneau et al. [14].

In the present paper, basic LMS equations for propagation of plane waves of any direction through a 2D PC slab are first provided in Section 2. The Poisson summation formula and the Graf addition theorem are utilized. As shown in Section 3, different types of scatterers at the same interface can be allowed, which represents an extension as compared to the treatment in [12]. Coupled equation systems are derived for the different scatterer types. For 3D PC slabs, a corresponding extension has recently proved useful for applications to design of anechoic coatings [15].

Applications to design of 2D PC slabs with broad transmission gaps for incident waves of in-plane as well as out-of-plane directions are presented in Section 4. The design problem is equivalent to a nonlinear optimization problem, and differential evolution [16], a global optimization technique from inverse theory, is used. A 2D PC slab example from Mei et al. [12], with lead cylinders in an epoxy host, is revisited. Geometrical and material parameter values are varied to position and widen its band gap.

2. Basic 2D Layer Multiple-Scattering Computational Method

As in Figure 1, a right-hand Cartesian 𝑥𝑦𝑧 coordinate system is introduced in a fluid-solid medium surrounded by homogeneous half-spaces. The horizontal directions are 𝑥 and 𝑦. The medium is periodic with period 𝑑 in the 𝑥 direction and uniform in the y direction.

317890.fig.001
Figure 1: Horizontal 𝑥𝑦 coordinates and a 𝑧 depth coordinate axis are introduced. There are three scatterer interfaces in this illustration. The medium is periodic with period 𝑑 in the 𝑥 direction, although only finitely many cylindrical scatterers are shown at each interface, and the cylinder axes are parallel to the y-axis.

Sound waves with time dependence exp(𝑖𝜔𝑡), to be suppressed in the formulas, are considered, where 𝜔 is the angular frequency. It follows that an incident plane wave with horizontal wavenumber vector 𝐤=(𝑘,𝜅,0) will give rise to a linear combination of reflected and transmitted plane waves with displacement vectors

𝐮(𝐫)=exp𝑖𝐊𝑠𝐠𝑗𝐫𝐞𝑗.(1) Here, 𝐫=(𝑥,𝑦,𝑧), 𝑗=1,2,3 for a wave of type P,SV,SH, respectively, 𝑠=+() for a wave in the positive (negative) 𝑧 direction, and

𝐊±𝐠𝑗=𝐤𝜔+𝐠±𝑐𝑗2||𝐤||+𝐠21/2=𝜔(0,0,1)𝑐𝑗(sin𝜃cos𝜑,sin𝜃sin𝜑,cos𝜃),(2) where g belongs to the reciprocal lattice

𝑘𝐠=𝑥,𝑘𝑦=,02𝜋𝑚𝑑,0,0(3) with 𝑚 running over the integers. Furthermore, 𝑐1 is the compressional-wave velocity 𝛼 and 𝑐2=𝑐3 are the shear-wave velocity 𝛽. The angular variables 𝜃,𝜑 of 𝐊±𝐠𝑗 are defined by (2), with a possibly complex cos𝜃. The vectors 𝐞𝑗=𝐞𝑗(𝐊±𝐠𝑗) are defined by 𝐞1=(sin𝜃cos𝜑,sin𝜃sin𝜑,cos𝜃), 𝐞2=(cos𝜃cos𝜑,cos𝜃sin𝜑,sin𝜃),  𝐞3=(sin𝜑,cos𝜑,0). It is convenient also to introduce the compressional- and shear-wave wavenumbers 𝑘𝑝=𝜔/𝛼 and 𝑘𝑠=𝜔/𝛽.

As detailed in [17], for example, and references therein, reflection and transmission matrices 𝑅𝐵,𝑇𝐵 and 𝑅𝐴,𝑇𝐴 can now be introduced, for the discrete set of waves specified by (1)–(3). The mentioned reference concerns the doubly periodic case for a 3D PC, with periodicity in the 𝑦 direction as well, but the R/T matrix formalism is the same. There are three scatterer interfaces in the illustration of Figure 1. Individual R/T matrices can be combined recursively [18, 19]. Layer thicknesses, as well as translations of individual scatterer interfaces in the x direction, are conveniently accounted for by phase shifts of the complex amplitudes of the plane-wave components.

2.1. Interface with Periodically Distributed Scatterers of a Common Type

Explicit expressions for the R/T matrices are well known for an interface between two homogeneous half-spaces [19]. To handle an interface with periodically distributed scatterers, as one of the three in Figure 1, the following cylindrical vector solutions to the wave equations can be used [12]:

𝐮𝐿𝑙𝜅𝑖(𝐫)=𝑘𝑝𝑓𝑙𝐮(𝑝𝑟)exp(𝑖𝑙𝜂)exp(𝑖𝜅𝑦),(4)𝑀𝑙𝜅1(𝐫)=𝑘𝑠×𝐮𝑁𝑙𝜅𝐮(𝐫),(5)𝑁𝑙𝜅1(𝐫)=𝑘𝑠𝑓×𝑙(𝑞𝑟)exp(𝑖𝑙𝜂)exp(𝑖𝜅𝑦)𝐞𝑦,(6) where cylindrical coordinates 𝑟,𝜂,𝑦 are used according to 𝐫=(𝑟sin𝜂,𝑦,𝑟cos𝜂), and 𝐞𝑦=(0,1,0). The index 𝑙=0,±1,± , and 𝜅 is a real number. Residual wavenumbers 𝑝, 𝑞 are defined by

𝑝=(𝑘2𝑝𝜅2)1/2𝑘,𝑞=2𝑠𝜅21/2.(7) The notation 𝐮0𝐿𝑙𝜅(𝐫), 𝐮0𝑀𝑙𝜅(𝐫), 𝐮0𝑁𝑙𝜅(𝐫) and 𝐮+𝐿𝑙𝜅(𝐫), 𝐮+𝑀𝑙𝜅(𝐫), 𝐮+𝑁𝑙𝜅(𝐫) is used for the two basic cases with 𝑓𝑙 as the Bessel function 𝐽𝑙 and 𝑓𝑙 as the Hankel function 𝐻𝑙(1), respectively.

For scatterers at

𝐑=(𝑥,𝑦,𝑧)=(𝑚𝑑,0,0),(8) where 𝑑 is the lattice period and 𝑚 runs over the integers, and for an incident plane wave as in (1), the total scattered field 𝐮sc can be written as (cf. [3]) 𝐮sc(𝐫)=𝑃𝑙𝑏𝑙+𝑃𝐑exp𝑖𝐤𝐑𝐮+𝑃𝑙𝜅,(𝐫𝐑)𝑃=𝐿,𝑀,𝑁,(9) where 𝜅 is the 𝑦 component of 𝐤.

The vector 𝐛+={𝑏𝑙+𝑃} is determined by solving the equation system

(𝐈𝑇Ω)𝐛+=𝑇𝐚0,(10) where 𝐈 is the appropriate identity matrix, 𝐚0={𝑎𝑙0𝑃} gives the coefficients for expansion of the incident plane wave in regular cylindrical waves 𝐮0𝑃𝑙𝜅(𝐫), Ω=Ω(𝑘𝑑,𝑝𝑑,𝑞𝑑) is the lattice translation matrix {Ω𝑃𝑃𝑙;𝑙}, and 𝑇={𝑇𝑃𝑃𝑙;𝑙} is the transition matrix for an individual scatterer. Specifically, 𝐛=Ω𝐛+ and 𝐛+=𝑇(𝐚0+𝐛) where 𝐛={𝑏𝑙𝑃} gives the coefficients for expansion in regular cylindrical waves 𝐮0𝑃𝑙𝜅(𝐫) of the scattered field from all scatterers except the one at the origin.

The R/T matrices are obtained, finally, by transforming the expansion (9) to plane waves of the type (1). Specifically, (9) can be rewritten as 𝐮sc(𝐫)=𝐠𝑗=1,2,3Δ𝑑𝐠,𝑗;𝑘𝑑,𝜅𝑑,𝑝𝑑,𝑞𝑑,𝐛+×exp𝑖𝐊±𝐠𝑗𝐫𝐞𝑗.(11) The sign in 𝐊±𝐠𝑗 is given by the sign of 𝑧.

Explicit expressions for the Δ coefficients in (11) are readily obtained from (9) and (4)–(6) by invoking, for 𝑘=𝑝 and 𝑘=𝑞, the relation

𝑚exp𝑖𝑚𝑑𝑘(𝑧+𝑖(𝑥𝑚𝑑))(𝑥𝑚𝑑)2+𝑧21/2𝑙×𝐻𝑙(1)𝑘(𝑥𝑚𝑑)2+𝑧21/2=2(𝑘)𝑙×𝑚𝛾𝑚𝑑1𝑘+2𝜋𝑚𝑑+𝑖𝛾𝑚sgn(𝑧)𝑙𝑖𝑘×exp+2𝜋𝑚𝑑𝑥+𝑖|𝑧|𝛾𝑚,(12) where 𝛾𝑚=[𝑘2(𝑘+2𝜋𝑚/𝑑)2]1/2 with Im 𝛾𝑚0. The relation (12) is valid for 𝑧0, and it can be verified using the Poisson summation formula.

2.2. Computation of the Expansion Coefficients 𝐚𝟎, the Matrices Ω, and the Matrices T

An incident plane compressional wave 𝐮inc(𝐫)=exp(𝑖𝑘𝑝𝐞inc𝐫)𝐞inc, where 𝐞inc=𝑘𝑝1(𝑝sin𝜂inc,𝜅,𝑝cos𝜂inc), can be expanded as

𝐮inc(𝐫)=𝑙i𝑙expi𝑙𝜂inc𝐮0𝐿𝑙𝜅(𝐫).(13) Noting that 𝐮inc(𝐫)=[exp(𝑖𝑘𝑝𝐞inc𝐫)]/𝑖𝑘𝑝, this follows readily from the well-known Bessel function relation exp(𝑖𝛾sin𝜂)=𝑙𝐽𝑙(𝛾)exp(𝑖𝑙𝜂).

An incident plane shear wave of SV or SH type can be expanded by a superposition of two cases. The first case, 𝐮inc(𝐫)=exp(𝑖𝑘𝑠𝐞inc𝐫)𝐞inc=(𝑘𝑠𝑞)1×{×[exp(𝑖𝑘𝑠𝐞inc𝐫)𝐞𝑦]} with 𝐞inc=𝑘𝑠1(𝑞sin𝜂inc,𝜅,𝑞cos𝜂inc) and 𝐞inc=𝑘𝑠1(𝜅sin𝜂inc,𝑞,𝜅cos𝜂inc), can be expanded as 𝐮inc𝑘(𝐫)=𝑠𝑞𝑙𝑖𝑙expi𝑙𝜂inc𝐮0𝑀𝑙𝜅(𝐫).(14)

The second case, 𝐮inc(𝐫)=exp(𝑖𝑘𝑠𝐞inc𝐫)𝐞inc=𝑖/𝑞×[exp(𝑖𝑘𝑠𝐞inc𝐫)𝐞𝑦] with 𝐞inc as before and 𝐞inc=(cos𝜂inc,0,sin𝜂inc), has the expansion 𝐮inc𝑘(𝐫)=𝑠𝑞𝑙𝑖𝑙+1exp𝑖𝑙𝜂inc𝐮0𝑁𝑙𝜅(𝐫).(15)

The lattice translation matrix Ω(𝑘𝑑,𝑝𝑑,𝑞𝑑) can be determined by applying, for 𝑘=𝑝 and 𝑘=𝑞, the relation 𝑚0exp𝑖𝑚𝑑𝑘(𝑧+𝑖(𝑥𝑚𝑑))(𝑥𝑚𝑑)2+𝑧21/2𝑙×H𝑙(1)𝑘(𝑥𝑚𝑑)2+𝑧21/2=𝑛(𝑧+𝑖𝑥)𝑥2+𝑧21/2𝑛Θ𝑙𝑛𝑘𝐽𝑑,𝑘𝑑𝑛(𝑘𝑟),(16) where the last sum is taken over all integers 𝑛 and Θ𝑙𝑘𝑑,𝑘𝑑=𝑖𝑙𝑚>0exp𝑖𝑚𝑑𝑘H𝑙(1)(𝑚𝑘𝑑)+𝑖𝑙𝑚>0exp𝑖𝑚𝑑𝑘H𝑙(1)(𝑚𝑘𝑑).(17) This relation follows from the Graf addition theorem [20] for Bessel functions. The elements {Ω𝑃𝑃𝑙;𝑙} of the lattice translation matrix Ω(𝑘𝑑,𝑝𝑑,𝑞𝑑) vanish unless 𝑃=𝑃. The remaining elements depend on 𝑙 and 𝑙 through |𝑙𝑙| only. Explicitly, Ω𝐿𝐿𝑙;𝑙=Θ𝑙𝑙𝑘,Ω𝑑,𝑝𝑑𝑀𝑀𝑙;𝑙=Ω𝑁𝑁𝑙;𝑙=Θ𝑙𝑙𝑘.𝑑,𝑞𝑑(18) Moroz [21] has published a representation of lattice sums in terms of exponentially convergent series, which has been used in the present work for numerical evaluation of the Θ quantities defined in (17).

For a homogeneous cylindrical scatterer, the interior field and the exterior field can be expanded in cylindrical waves 𝐮0𝑃𝑙𝜅(𝐫) and 𝐮0𝑃𝑙𝜅(𝐫),𝐮+𝑃𝑙𝜅(𝐫), respectively. An equation system for the 𝑇-matrix elements 𝑇={𝑇𝑃𝑃𝑙;𝑙}, which depend on 𝜅, is then readily obtained from the standard boundary conditions concerning continuity of displacement and traction at the cylinder surface. The scatterer as well as the host medium can be either fluid or solid. Because of the circular symmetry with a cylindrical scatterer, scattering only appears to the same 𝑙 component (𝑙=𝑙). Details concerning the case 𝜅=0, for which the 𝑃-SV (𝑃,𝑃𝑀) and SH (𝑃=𝑃=𝑀) solutions decouple, are given by Mei et al. [9].

3. Different Types of Scatterers at the Same Interface

The LMS method is commonly applied for lattices with identical cylindrical scatterers within the same layer. As shown below, however, different types of scatterers within the same layer can also be accommodated. A similar extension for the restriction to in-plane wave propagation is made in Ivansson [22].

An illustration is given in Figure 2. Centered at the same 𝑧 level, at each of the three scatterer interfaces, there are two types of cylindrical scatterers. Each scatterer interface is treated separately. Choosing coordinates appropriately, scatterers of the first type, with transition matrix 𝑇 and scattered-field expansion coefficients denoted 𝐛+, appear at 𝐑=𝑚(𝑑,0,0), for integers 𝑚. Scatterers of the second type, with transition matrix 𝑈 and scattered-field expansion coefficients denoted 𝐜+, appear at points 𝐒 in between, that is, 𝐒=(𝑚+1/2)(𝑑,0,0). The reciprocal lattice vectors become 𝐠=(2𝜋𝑚/𝑑,0,0), where m runs over the integers.

317890.fig.002
Figure 2: The configuration from Figure 1 is extended here, by allowing two types of cylindrical scatterers to appear in an alternating fashion at the same interface depth. The medium is periodic with period d in the x direction, and the cylinder axes are still parallel to the y-axis.

The generalization of the expression (9) for the scattered field becomes

𝐮sc(𝐫)=𝑃𝑙𝑏𝑙+𝑃𝐑exp𝑖𝐤𝐑𝐮+𝑃𝑙𝜅+(𝐫𝐑)𝑃𝑙𝑐𝑙+𝑃𝐒exp𝑖𝐤𝐒𝐮+𝑃𝑙𝜅.(𝐫𝐒)(19) It follows that 𝐛+𝐚=𝑇0+𝐛+𝐛,𝐜+𝐚=𝑈0+𝐜+𝐜,(20) where, for a scatterer of the first type at 𝐑, exp(𝑖𝐤𝐑)𝐛 and exp(𝑖𝐤𝐑)𝐛give the coefficients for expansion in regular cylindrical waves 𝐮0𝑃𝑙𝜅(𝐫𝐑) of the scattered field from all other scatterers of the first and the second types, respectively. The vectors 𝐜 and 𝐜 are defined analogously. For a scatterer of the second type at 𝐒, exp(𝑖𝐤𝐒)𝐜 and exp(i𝐤𝐒)𝐜 thus give the coefficients for expansion in regular cylindrical waves 𝐮𝑙𝜅0𝑃(𝐫𝐒) of the scattered field from all other scatterers of the second and the first types, respectively.

With Ω0=Ω(𝑘𝑑,𝑝𝑑,𝑞𝑑), it follows that 𝐛=Ω0𝐛+ and 𝐜=Ω0𝐜+. For a certain matrix Ωdif, to be determined, 𝐛=Ωdif𝐜+ and 𝐜=Ωdif𝐛+. The equation system for determination of 𝐛+ and 𝐜+ becomes 𝐈𝑇Ω0𝐛+𝑇Ωdif𝐜+=𝑇𝐚0,𝑈Ωdif𝐛++𝐈𝑈Ω0𝐜+=𝑈𝐚0.(21)

In order to form the R/T matrices, incident plane waves with different horizontal wavenumber vectors 𝐤+𝐠inc=(𝑘+𝑔inc,𝜅,0) have to be considered, where 𝑔inc belongs to the set {2𝜋𝑚/𝑑} providing the reciprocal lattice. Noting that the union of the scatterer positions is a small square lattice with period 𝑑/2, the following expression for Ωdif as a difference of Ω matrices is directly obtained Ωdif𝑘=Ω+𝑔inc𝑑2,𝑝𝑑2,𝑞𝑑2Ω0.(22) Only those 𝐠inc in {(2𝜋𝑚/𝑑,0,0)} for which 𝑚 is even are reciprocal vectors for the small lattice with period 𝑑/2. Since a lattice translation matrix Ω is periodic in its first argument with period 2𝜋, there will be two groups of 𝐠inc with different Ωdif matrices according to (22). Specifically, Ωdif,even𝑘=Ω𝑑2,𝑝𝑑2,𝑞𝑑2Ω0(23) pertains to 𝐠inc=(2𝜋𝑚/𝑑,0,0) with even 𝑚, and Ωdif,odd𝑘=Ω𝑑2+𝜋,𝑝𝑑2,𝑞𝑑2Ω0(24) pertains to 𝐠inc = (2π m/d, 0, 0) with odd m.

The transformation of the expansion (19) to plane waves of the type (1) can be done separately for each of the 𝐑 and 𝐒 sums. In the latter case, the translation from the origin causes a sign change for some combinations of incident (𝐠inc) and scattered (𝐠sc) reciprocal lattice vectors. Specifically, with 𝐠inc=(2𝜋𝑚inc/𝑑,0,0) and 𝐠sc=(2𝜋𝑚sc/𝑑,0,0), the double sum corresponding to the one in (11) appears as 𝐮sc(𝐫)=𝐠sc𝑗=1,2,3(1)𝑚inc𝑚scΔ𝑑𝐠sc,𝑗;𝑘𝑑,𝜅𝑑,𝑝𝑑,𝑞𝑑,𝐜+×exp𝑖𝐊±𝐠sc𝑗𝐫𝐞𝑗.(25)

4. Designing 2D PC Slabs with a Transmission Gap

It is well known that band gaps can appear when scatterers with a large density are arranged periodically in a host with a small density. A particular 2D PC slab example with lead cylinders in epoxy was considered in Mei et al. [12]. The epoxy parameters were 2540 and 1159.817 m/s for the compressional- and shear-wave velocities, respectively, and 1.18 kg/dm3 for the density. Corresponding lead material parameters were 2160 and 860.568 m/s, and 11.4 kg/dm3. The slab was formed by sixteen layers of lead cylinders, each with a radius of 3.584 mm, with a spacing between cylinder centers of 11.598 mm in the 𝑥 as well as 𝑧 directions, cf. Figure 1 where there are three layers. A band gap centered at about 60 kHz was found. (Dimensionless frequencies and distances were actually used in the paper, but a specialization is made here.)

Figure 3 is a contour plot of the transmittance, that is, time- (and space-) averaged transmitted energy flux relative to the incident one, up to 120 kHz. The vertical axis is for the angle 𝜃, where the direction vector of the incident compressional plane wave is (sin𝜃cos𝜑, sin𝜃sin𝜑, cos𝜃). In-plane and out-of-plane incidence angles are considered together in Figure 3, where the upper half with 𝜑=0 concerns incidence in the 𝑥𝑧 plane (in-plane propagation) and the lower half with 𝜑=90 concerns incidence in the 𝑦𝑧 plane. Only the first Brillouin zone is involved, since 𝑘𝑥=𝜔/𝛼sin𝜃cos𝜑<2𝜋/𝑑when 𝛼 equals the epoxy compressional-wave velocity 2540 m/s, 𝑑=11.598 mm, and the frequency is less than 𝛼/𝑑=219 kHz. Of course, the band gap from Mei et al. [12], at about 60 kHz, shows up clearly in this kind of plot. One might wonder how the geometrical and/or material parameters of the slab should be modified to achieve a prescribed desired change of the appearance of the gap.

317890.fig.003
Figure 3: Contour plot of transmittance for a PC slab with lead cylinders in an epoxy host. Parameters of the slab are given in the text. The direction vector of the incident compressional wave is (sin 𝜃 cos 𝜑, sin 𝜃 sin φ, cos 𝜃). The black and gray contours are at −150 dB and −10 dB, respectively.

Global optimization methods can be used to design PC slabs with desirable properties. Simulated annealing, genetic algorithms and differential evolution (DE) are three kinds of such methods, that have become popular during the last fifteen years. DE, to be applied here, is related to genetic algorithms, but the parameters are not encoded in bit strings, and genetic operators such as crossover and mutation are replaced by algebraic operators [16].

As a very simple example, to try to position and widen the gap in Figure 3, an objective function for DE minimization is specified as the maximum transmittance in the frequency band 45–65 kHz when the incidence angle 𝜃 is varied between 0 and 20 for the two azimuthal angles 𝜑=0 (in-plane propagation in the 𝑥𝑧 plane) and 𝜑=90 (propagation in the 𝑦𝑧 plane). A configuration with cylinders of two alternating sizes is allowed, as depicted in Figure 2. Still, the slab is formed by sixteen layers of cylinders. The following seven parameters are varied within the indicated search space: scatterer compressional-wave velocity 𝛼 [1500 m/s 𝛼 3000 m/s], scatterer shear-wave velocity 𝛽 [650 m/s 𝛽 1200 m/s], scatterer density 𝜌 [7 kg/dm3 𝜌 14 kg/dm3], layer thickness [8 16 mm], the largest scatterer radius 𝑟max [0.225 𝑟max 0.4 ], the smallest scatterer radius 𝑟min [0.45 𝑟max𝑟min𝑟max], the lattice period 𝑑[8(𝑟max+𝑟min)/3 𝑑 4 (𝑟max+𝑟min)]. The layer thickness is the 𝑧 distance between subsequent scatterer interfaces (cf. Figure 2).

Table 1 shows the optimum obtained with DE, along with corresponding parameter values. A reduction of the transmitted field with more than 160 dB is achieved throughout the band 45–65 kHz and throughout the solid angle intervals 0<𝜃<20 for 𝜑=0, 𝜑=90 for the direction of incidence. High-velocity high-density cylinders seem preferable, since the lead velocities and density are all increased to values close to the upper ends of the corresponding search intervals. In this case, as large as possible values of 𝛼, 𝛽, and 𝜌 could in fact have been fixed from the start. Moreover, the cylinders are almost as densely packed horizontally as allowed by the search interval for 𝑑. Hence, the optimization could in fact have been simplified considerably. With only a few free parameters, a complete search, for example, can be a feasible alternative.

tab1
Table 1: Specification of a PC slab that has been optimized by DE to produce small transmittance in the band 45–65 kHz. Corresponding transmittance results are shown in Figures 4 and 5.

Compared to Figure 3, the contour plot in Figure 4 of the transmittance for the optimized PC slab shows a band gap around 60 kHz that has become significantly wider and also deeper. The allowed increases of 𝛼, 𝛽, 𝜌, and (𝑟max+𝑟min)/𝑑 are certainly essential for producing this effect. Variation of , 𝑟max, and 𝑟min is needed to position the band gap at the desired frequency interval. Although the optimization was performed with the restriction 0<𝜃<20, small transmittance is apparently achieved for all angles 𝜃 (0<𝜃<90).

317890.fig.004
Figure 4: Contour plot of transmittance as in Figure 3 but for the optimized PC slab specified in Table 1.

It turns out that the transmittance for the optimized PC slab remains small within the 45–65 kHz band for all out-of-plane incidence angles (an “absolute” band gap). Figure 5 shows the transmittance in the same way as Figure 4, but for 𝜑=0 changed to 𝜑=30 and 𝜑=90 changed to 𝜑=60. The upper and lower halves exhibit increased symmetry, since the 𝜑 angles involved are closer to one another.

317890.fig.005
Figure 5: Contour plot of transmittance for the optimized PC slab specified in Table 1. The difference to Figure 4 is that two other 𝜑 angles are considered for the direction of the incident compressional wave.

Apparently, the optimized PC slab has large transmittance below about 30 kHz and there are regions with large transmittance above 80 kHz as well. The band gap with small transmittance extends to higher frequencies, than those in the band 45–65 kHz, for plane-wave incidence directions with either large 𝜃 or small 𝜑.

The efficacy of the DE technique can be illustrated by showing the decrease of the maximum transmittance within the specified frequency/angle region as the number of tested parameter settings for the PC slab evolves. Figure 6 shows this decrease for the example from Figure 4. A comparison to the much less efficient MonteCarlo method, with random selection of parameters from the search space, is included. The low efficiency of the MonteCarlo approach shows that the desired PC slabs with small transmittance only appear in a small portion of the search space. For example, only about 0.15% of the random slab selections had a maximum transmittance below −20 dB within the specified frequency/angle region. About 0.01% had a maximum transmittance below −100 dB.

317890.fig.006
Figure 6: Evolution of the maximum transmittance with the number of tested parameter combinations for the DE optimization leading to the PC slab specified in Table 1. A corresponding curve for a brute force Monte Carlo method is also included (gray).

The difference between 𝑟max and 𝑟min in Table 1 is rather small. Optimization was also tried with all cylinders of exactly the same radius, as in Figure 1. The obtained optimum, with cylinders of radius 4.2503 mm and slightly modified values for the other parameters, was almost as good as the one presented in Table 1. The advantages with cylinders of different sizes are expected to be more significant in more complicated filtering cases, involving, for example, more than one frequency band.

5. Conclusions

The layer multiple-scattering (LMS) method is a fast semianalytical technique for computing scattering from layers including periodic scatterer lattices. For the 2D case with cylindrical scatterers and any solid angle direction of an incident plane wave, an extension has been made to scatterer lattices with cylindrical scatterers of two different sizes in the same horizontal plane.

Global optimization methods from inverse theory are useful for designing PC slabs with desirable properties. A differential evolution algorithm has been applied here to position and widen an “absolute” band gap for a certain 2D PC slab. Although only limited angle intervals for the direction of incidence were included for the optimization, small transmittance was achieved for all solid angles specifying the direction of the incident wave.

The possibility to include cylindrical scatterers of two different sizes in the same horizontal plane provides an additional degree of freedom that can be useful for PC design purposes. In the presented example, with its specification of objective function and search intervals for the parameters, however, the difference between the optimal cylinder radii was rather small and the improvement only marginal. The additional flexibility is expected to be more important in more complicated filtering applications. For example, specified regions in frequency-angle space with large transmittance could be desired in addition to specified regions with small transmittance.

Acknowledgment

Alexander Moroz kindly provided his Fortran routine for calculation of lattice sums.

References

  1. T. Miyashita, “Sonic crystals and sonic wave-guides,” Measurement Science and Technology, vol. 16, no. 5, pp. R47–R63, 2005. View at Publisher · View at Google Scholar
  2. N. Papanikolaou, R. Zeller, and P. H. Dederichs, “Conceptual improvements of the KKR method,” Journal of Physics: Condensed Matter, vol. 14, no. 11, pp. 2799–2823, 2002. View at Publisher · View at Google Scholar
  3. I. E. Psarobas, N. Stefanou, and A. Modinos, “Scattering of elastic waves by periodic arrays of spherical bodies,” Physical Review B, vol. 62, no. 1, pp. 278–291, 2000.
  4. Z. Liu, C. T. Chan, P. Sheng, A. L. Goertzen, and J. H. Page, “Elastic wave scattering by periodic structures of spherical objects: theory and experiment,” Physical Review B, vol. 62, no. 4, pp. 2446–2457, 2000.
  5. M. Sigalas, M. S. Kushwaha, E. N. Economou, M. Kafesaki, I. E. Psarobas, and W. Steurer, “Classical vibrational modes in phononic lattices: theory and experiment,” Zeitschrift für Kristallographie, vol. 220, no. 9-10, pp. 765–809, 2005. View at Publisher · View at Google Scholar
  6. R. Sainidou, N. Stefanou, I. E. Psarobas, and A. Modinos, “The layer multiple-scattering method applied to phononic crystals,” Zeitschrift für Kristallographie, vol. 220, no. 9-10, pp. 848–858, 2005. View at Publisher · View at Google Scholar
  7. S. B. Platts, N. V. Movchan, R. C. McPhedran, and A. B. Movchan, “Two-dimensional phononic crystals and scattering of elastic waves by an array of voids,” Proceedings of the Royal Society of London A, vol. 458, no. 2026, pp. 2327–2347, 2002. View at Publisher · View at Google Scholar
  8. S. B. Platts, N. V. Movchan, R. C. McPhedran, and A. B. Movchan, “Transmission and polarization of elastic waves in irregular structures,” Transactions of the ASME, vol. 125, no. 1, pp. 2–6, 2003. View at Publisher · View at Google Scholar
  9. J. Mei, Z. Liu, J. Shi, and D. Tian, “Theory for elastic wave scattering by a two-dimensional periodic array of cylinders: an ideal approach for band-structure calculations,” Physical Review B, vol. 67, Article ID 245107, 7 pages, 2003.
  10. C. Qiu, Z. Liu, J. Mei, and M. Ke, “The layer multiple-scattering method for calculating transmission coefficients of 2D phononic crystals,” Solid State Communications, vol. 134, no. 11, pp. 765–770, 2005. View at Publisher · View at Google Scholar
  11. S. Robert, J.-M. Conoir, and H. Franklin, “Propagation of elastic waves through two-dimensional lattices of cylindrical empty or water-filled inclusions in an aluminum matrix,” Ultrasonics, vol. 45, no. 1–4, pp. 178–187, 2006. View at Publisher · View at Google Scholar · View at PubMed
  12. J. Mei, Z. Liu, and C. Qiu, “Multiple-scattering theory for out-of-plane propagation of elastic waves in two-dimensional phononic crystals,” Journal of Physics: Condensed Matter, vol. 17, no. 25, pp. 3735–3757, 2005. View at Publisher · View at Google Scholar
  13. M. Wilm, A. Khelif, S. Ballandras, V. Laude, and B. Djafari-Rouhani, “Out-of-plane propagation of elastic waves in two-dimensional phononic band-gap materials,” Physical Review E, vol. 67, no. 6, Article ID 065602, 4 pages, 2003.
  14. S. Guenneau, C. G. Poulton, and A. B. Movchan, “Oblique propagation of electromagnetic and elastic waves for an array of cylindrical fibres,” Proceedings of the Royal Society of London A, vol. 459, no. 2037, pp. 2215–2263, 2003. View at Publisher · View at Google Scholar
  15. S. M. Ivansson, “Numerical design of Alberich anechoic coatings with superellipsoidal cavities of mixed sizes,” Journal of the Acoustical Society of America, vol. 124, no. 4, pp. 1974–1984, 2008. View at Publisher · View at Google Scholar · View at PubMed
  16. K. Price, R. Storn, and J. Lampinen, Differential Evolution—A Practical Approach to Global Optimization, Springer, New York, NY, USA, 2005.
  17. S. M. Ivansson, “Sound absorption by viscoelastic coatings with periodically distributed cavities,” Journal of the Acoustical Society of America, vol. 119, no. 6, pp. 3558–3567, 2006. View at Publisher · View at Google Scholar
  18. J. B. Pendry, Low Energy Electron Diffraction, Academic Press, New York, NY, USA, 1974.
  19. B. L. N. Kennett, Seismic Wave Propagation in Stratified Media, Cambridge University Press, Cambridge, UK, 1983.
  20. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, National Bureau of Standards, 1965.
  21. A. Moroz, “Exponentially convergent lattice sums,” Optics Letters, vol. 26, no. 15, pp. 1119–1121, 2001.
  22. S. M. Ivansson, “Extended layer multiple-scattering and global optimization techniques for 2D phononic crystal insulator design,” in Photonics Research Developments, V. P. Nilsson, Ed., chapter 7, pp. 225–245, Nova, 2008.