Rayleigh–Bloch waves along elastic diffraction gratings

Rayleigh–Bloch (RB) waves in elasticity, in contrast to those in scalar wave systems, appear to have had little attention. Despite the importance of RB waves in applications, their connections to trapped modes and the ubiquitous nature of diffraction gratings, there has been no investigation of whether such waves occur within elastic diffraction gratings for the in-plane vector elastic system. We identify boundary conditions that support such waves and numerical simulations confirm their presence. An asymptotic technique is also developed to generate effective medium homogenized equations for the grating that allows us to replace the detailed microstructure by a continuum representation. Further numerical simulations confirm that the asymptotic scheme captures the essential features of these waves.


Introduction
For scalar wave systems such as those in linear water waves, acoustics, polarized electromagnetic, flexural waves or shear-horizontal polarized elastic waves, it has long been known that one can find a class of surface waves created by periodicity in the geometry of the surface or equivalently by a periodic diffraction grating embedded within the medium. These waves propagate along a surface, exponentially decay in amplitude perpendicular to the surface and are created by geometric periodic corrugations, or perturbations, to the surface [1][2][3]   would otherwise not exist. Although often studied for surface waves, they also propagate along diffraction gratings and provide a mechanism for energy transport and guiding along gratings [4] and arrays [5]. Owing to their ubiquitous nature in these different settings, they have been discovered and rediscovered several times and appear as edge waves [6] for water waves localized to periodic coastlines, spoof surface plasmon polaritons (SPPs) [7,8] in modern applications of plasmonics, array-guided surface waves [9] in Yagi-Uda antenna theory, Rayleigh-Bloch surface waves [10,11] for scalar wave systems with diffraction gratings, and in elastic plates [12,13] among other areas. The name Rayleigh-Bloch (RB hereinafter) waves seems most descriptive as surface waves are typically called Rayleigh waves and Bloch waves arise due to periodicity. RB waves have been of particular interest in terms of their existence [14,15] as they require either Neumann boundary conditions upon the grating, or penetrable scatterers (see [16] among others), to exist and they have implications, as non-zero eigenfunctions of the system, in terms of uniqueness of solution and connection to trapped modes; the Dirichlet cases cannot support RB waves [10]. Curiously, given this substantial literature, the importance of RB waves in applications, their connections to trapped modes [17] and the universal nature of diffraction gratings [18] there has been no investigation of whether RB waves occur, or even exist, for elastic diffraction gratings for the in-plane vector elastic system; we address the diffraction grating here. For the scalar case of, say, an array of Neumann cylindrical voids [19,20] (figure 1) one finds that there is a dispersion relation connecting the phase shift along the array to the frequency and RB waves exist beneath the so-called light-line (the linear dispersion relation from the surrounding bulk medium; figure 2a). In the low-frequency limit, the RB dispersion curve approaches the bulk dispersion relation and in many situations, for instance, comb-like surfaces [2,6], one can have multiple RB waves beneath the light-line.
In isotropic linear elasticity, the bulk medium now supports two waves, compressional and shear, with differing wave speeds and therefore two light-lines exist (e.g. figure 2b), moreover the natural boundary conditions in elasticity for a realistic grating are either rigid or traction free; neither are simply Neumann or Dirichlet and the fields are coupled at these boundaries, that is, conversion of pressure (resp. shear) to shear (resp. pressure) waves occur at the boundaries. It is therefore not obvious whether the theory underlying scalar RB waves will translate to the elastic system. We identify the boundary conditions that can support such waves, numerical simulations confirm their presence and an asymptotic technique is developed to generate effective medium equations for them. It is worth noting that, recently, RB waves have been studied in the setting of the full Maxwell's equations for electromagnetic waves propagating along a linear array of dielectric spheres [21]. Yet, despite several papers analysing scattering problems for arrays of spheres and cylinders in full vector elasticity (e.g. [22][23][24][25][26]), to the best of our knowledge there has been no study of the existence or behaviour of RB waves for the Lamé system.
We will couch our study in the language of homogenization theory as we can then simultaneously generate RB waves and an asymptotic theory for them. We shall view the  figure 1). Here, Ω is radian frequency and k is the Bloch wavenumber. (a) For the acoustic case (e.g. anti-plane shear waves) of an array of Neumann cylindrical voids; (b) for the elastic case (coupled in-plane shear and pressure waves) where the boundary conditions are traction free. The voids have radius 0.8 and = 1 in both cases. The (red) dotted lines show the light-lines in each case. The reader's attention is drawn to the different scales on the ordinate axes. We also emphasize that the lowest dotted line in (b) corresponds to in-plane shear waves, whereas the upper one which is for in-plane pressure waves is markedly different from that of acoustics waves in (a). (a) Acoustic case and (b) elastic case.
diffraction grating, as shown in figure 1, as consisting of two length scales; a short-scale set by the periodicity and a long scale. The long-scale can be thought of as the overall size of the grating, 1 or the typical macroscopic wavelength dominating the dynamic response (cf. §4a), or some other physical length that is large compared with the microscopic cell size. The precise physical nature of the long scale is unimportant for the analysis which follows; it is sufficient that there exists some length L such that /L 1. This separation of scales is then a natural setting for multiple scales with a fast scale and slow scale in space and thus for homogenization. Classically, homogenization theory is developed for long waves that are of low frequency, and for bulk media (not surfaces) and is detailed in many monographs, [27][28][29][30], and essentially relies upon the wavelength being much larger than the microstructure which is usually assumed to be perfectly periodic: this limits the procedure to low frequencies/long waves and a quasi-static situation. For many applications, particularly in optics such as photonics, the limitation to low frequencies is a serious deficiency, nonetheless, the attraction of having an effective equation for a microstructured medium where one need no longer model the detail of each individual scatterer is highly attractive [31]. The classical theory has been recently extended in three directions: to higher frequencies where the high-frequency homogenization of [32] provides the necessary methodology (although limited to scalar waves) with short-scale structure within the solution captured by standing Bloch waves and this then modulated by a long-scale envelope function that satisfies a partial differential equation. The upshot of the theory is that an effective equation is deduced entirely on the long scale with the microstructure encapsulated by spatially dispersive coefficients (e.g. [33]). Recently, the theory has also been extended to in-plane elasticity for bulk media [34,35] and to scalar RB waves [20]. Therefore, the theory is ideally suited to the exploration of elastic RB waves.
We begin with a brief discussion of RB waves in §2, which is followed by a formulation of the problem in the asymptotic setting of high-frequency homogenization ( §3). The asymptotic theory then leads naturally to an effective condition for the diffraction grating, as an ordinary differential equation on the long scale, that supports waves that travel along the array. These RB waves are then found numerically in §4 for an array of cylindrical voids that have a traction-free boundary condition upon them. Finally, some concluding remarks are drawn together in §5.

Rayleigh-Bloch waves
As discussed earlier, RB waves have a long history for scalar wave systems that are equivalent to the Helmholtz equation. In the bulk medium, which is dispersionless, there is a linear relation between wavenumber and frequency; this linear bulk dispersion relation is the light-line shown in figure 2a as the dotted red line. For a periodic linear array, it is conventional to relate the phase shift along the array, characterized by the Bloch wavenumber (defined here by the Bloch-Floquet quasi-periodicity condition such that u| ∂S − = u| ∂S + e i2k where ∂S + denote the boundaries of the elementary cell as illustrated in figure 1) to frequency (Ω). In this paper, it is convenient to work with a non-dimensional Bloch wavenumber and to this end, we take the semi-width of the cell as a natural unit and set = 1. To have eigensolutions that decay, as one tends to infinity in directions perpendicular to the array, one typically needs to be below the light-line to prevent coupling into radiation. There are some counterexamples to this, in particular, there are embedded RB waves [17] but we will not discuss them or their analogues in elasticity here. For a linear array of Neumann cylindrical inclusions, with radius 0.8, for the Helmholtz equation the RB dispersion curves are shown in figure 2a. Notably, as the radius increases further beyond 0.81, another RB mode emerges from the light-line and this feature is also seen in the elastic case; the standing wave frequencies (Ω when k = π/2 and, hence, u| ∂S − = −u| ∂S + ) for the lowest four modes are shown as a function of radius in figure 3b with one emerging near 0.8. A good intuitive guide to the possible existence of RB waves for an array or grating is to consider a doubly periodic rectangular array of inclusions, say, of fixed width, 2, in x 1 and of height d in x 2 . The limit of d = 2 is that of a doubly periodic array of cylinders in a square array and the limit of d → ∞ is that of the linear array we consider here. In the scalar case, the square array is an oft-studied system [36] and the Dirichlet case has a zero frequency stop-band, that is, the dispersion curves lie above the light-line for d = 2 and never cross beneath it as d increases. Unsurprisingly, therefore one observes no RB waves for a linear array of Dirichlet cylinders. Conversely, for a square array of Neumann cylinders there is a low-frequency linear dispersion relation through the origin, beneath the light-line, and as the radius increases an additional mode passes beneath the light-line-both remain as d increases. We now parallel this intuitive process for elastic inclusions: the two cases of primary interest are rigid inclusions or traction-free inclusions.
On the microstructure, shown in figure 1, S 0 represents the inclusion and we could set either traction-free or rigid conditions on its boundary (∂S 0 ) depending on the physical setting under consideration; for the sake of clarity, we do not consider the case of penetrable scatters in the  present paper. Figure 4a shows the first four dispersion curves for a doubly periodic array of rigid voids for two-dimensional in-plane elasticity (e.g. [35]). In this case, there is a zero frequency stop-band and, hence, we would not expect to see RB-type modes; this is analogous to the scalar Dirichlet case mentioned earlier.
On the other hand, parts (b)-(d) of the same figure show the corresponding curves for the same geometry, but for voids with traction-free boundaries, and for cells of increasing height d = 2, 5 and 10, respectively. In these cases, we observe two acoustic modes emerging from the origin, both of which, eventually lie beneath the shear light-line in the vicinity of k = π/2 and that resemble the lowest two dispersion curves in figure 2b. These also have a similarity with the RB dispersion curve shown in the acoustic case (cf. figure 2a).
As the height of the elementary cell increases, that is, as the doubly periodic geometry approaches that of the diffraction grating, the higher frequency modes collapse down onto the shear light-line (which is the lower one), with a key difference between the elastic and acoustic cases, being the persistence of a RB mode near the compressional light-line as illustrated in figure 2b. The fifth mode can also be observed in figure 4d where the doubly periodic array of voids tends to the grating geometry.
It is notable that the numerical results were produced using the commercial finite-element software Comsol Multiphysics and that the infinite strip is created computationally using (elastic) perfectly matched layers; these are standard numerical approaches. There is, however, an issue with a large number of spurious modes detected that lie in-between the light-lines. We have carefully tracked the decaying RB modes and have left gaps in the third and fifth RB modes in figure 2b where the accuracy was diminished by the proximity of spurious modes. However, as we shall see, the asymptotic method developed herein will allow us to extract all the salient information regarding the behaviour of RB waves despite the numerical difficulties associated with the full computation.

Diffraction gratings: an asymptotic analysis
The time-harmonic deformation of a linearly elastic solid is governed by the elastodynamic wave equation [37] where repeated indices are summed over hereafter unless otherwise stated and subscript commas followed by lowercase letters indicate differentiation with respect to the spatial variables. The stresses are denoted by σ ij , displacements by u i , and the density and angular frequency are ρ and ω, respectively. We choose to use a Cartesian coordinate frame, x i , with the diffraction grating orientated in the x 1 direction. As we assume homogeneous isotropic linearly elastic media, the constitutive equation relating stress and strain is with λ, μ as the Lamé coefficients and δ ij as the Kronecker delta.
We are concerned with in-plane elastodynamic waves propagating along a linear microstructured grating such as that shown in figure 1. In particular, it is assumed that there are two natural length scales: (i) a short scale > 0 characterizing the dimension of the microstructure and (ii) a long scale L > 0 characterizing the overall size of the macrostructure (or macrocell). This being the case, it is natural to introduce a slow variable X = x 1 /L along the array and a fast variable . An important feature of the asymptotic approach is that we introduce these two scales in space and that they are now treated as independent and therefore this builds the mismatch of scales into the equations.
In these new variables, the elastodynamic wave equation is then expressed as where subscript commas followed by lower case letters indicate differentiation with respect to the fast variables ξ j and subscript commas followed by X indicate differentiation with respect to the slow variable X. Subscript commas followed by multiple characters indicate successive differentiation. In the rescaled variables, the constitutive relation (3.2) is Combining the rescaled equations (3.3) and (3.4) yields the governing equation of time-harmonic deformation in fast and slow variables where the non-dimensional frequency-squared Ω 2 = 2 ρω 2 /μ and the ratio Λ = λ/μ have been introduced.
(a) General, high-frequency, homogenization The form of (3.5) suggests the following anstatz for the displacement field and frequency: whence a hierarchy of equations in ascending orders of is obtained (there is no need to go beyond second order for the present theory): and k,k δ ij and ξ ∈ S, where S is the elementary cell which contains some microstructural element S 0 (see figure 1, for example).
The RB waves travel along an infinite, periodic or quasi-periodic, array and decay exponentially in the directions perpendicular to the array. The approach here is to create a perturbation about standing wave solutions which are either in-phase or out-of-phase across the elementary cell ( figure 1). In this case, the elementary cell is the infinite strip S = {ξ : |ξ 1 | < 1, |ξ 2 | < ∞} \ S 0 which contains some microstructural element S 0 with boundary ∂S 0 . The exterior boundaries of the cell are denoted by ∂S ± = {ξ : ξ 1 = ±1} (see figure 1, for example). The in-phase and out-of-phase conditions on the edge of the cell are then with + (resp. −) for the in-phase (resp. out-of-phase) solutions. Here, n j are the components of the outward pointing normal. As discussed in §2, we require that the microstructure is traction free for RB waves, and then σ ij n j | ∂S 0 = 0.
1,X n i (3.13) and 1,X n i . (3.14) Returning to the O(1) problem and assuming that Ω 2 0 is a simple eigenvalue (eigenvalues of multiplicity greater than one will be discussed later in the context of low-frequency homogenization), the leading-order solution admits the decomposition It initially appears counterintuitive that the long-scale behaviour reduces to a single scalar function f (0) (X): as we are working with the vector system of in-plane elasticity we might naturally expect two scalar functions analogous to the shear and compressional potentials. Nevertheless, the leading-order problem formally permits the aforementioned decomposition and the elastic high-frequency homogenization theory of [35] also generates a similar situation. As we shall see in §3b, two scalar functions emerge naturally from the asymptotic analysis near the repeated root at zero frequency, which yield the expected shear and compressional potentials and, hence, the two light-lines emerging from the origin in figure 2b.
In order to determine Ω 2 1 , the following equality, formed by multiplying the O(1) equation by u (1) i and subtracting the O( ) equation multiplied by U subject to the boundary condition where we introduce the short-scale stress Σ (n) (ξ ) as Σ k,k δ ij . The next-toleading order vector U (1) is not required for the determination of f (0) and, hence, is irrelevant for our purposes here.
Moving to the second-order problem and forming a similar solvability condition as for the O( ) case the following homogenized differential equation for the envelope function f (0) (X) is obtained On the long-scale, the effective medium is then entirely characterized by (3.18) which is equivalent to that of a one-dimensional string of stiffness T = A/B. It is emphasized that the problem for the scalar envelope function (3.18) is cast entirely on the long-scale, with the short-scale properties of the microstructure encoded in the T-values. The overall leading-order behaviour is the rapid short-scale oscillations described by the leading-order eigenvector U (0) , modulated by the longscale oscillations of f (0) . Indeed, for waves of frequency Ω propagating along the diffraction grating, the amplitude of the short-scale oscillations are modulated by a long-scale oscillation of wavelength λ = 2π T/(Ω 2 − Ω 2 0 ). The homogenized differential equation (3.18) can be used to obtain the asymptotic dispersion relation in the vicinity of the standing wave frequencies at k = π/2. The Bloch-Floquet quasiperiodicity condition u| ∂S − = u| ∂S + e i2k requires that the solution be anti-periodic with respect to ξ 1 at k = π/2. As U (0) (ξ ) is anti-periodic by construction, we choose f (0) (X) = exp(i[π/2 − k]X/ ) which, after substitution into (3.18), yields whence the asymptotic dispersion equation is

(b) Classical, low-frequency, homogenization
In the low-frequency regime, that is when Ω 2 = O( 2 ), the above procedure is still valid but there are some important technical differences. The boundary value problem for the leading order is ij,j = 0, for ξ ∈ S and σ (0) ij n j = 0, for ξ ∈ ∂S 0 , (3.23) and periodic conditions are imposed on the lateral boundaries of the elementary strip. The leading-order boundary value problem has an eigenvalue Ω 0 = 0 of multiplicity two and, hence, admits a solution of the form where repeated indices are summed over and U which admits solutions of the form Finally, the spectral parameter appears in the second-order problem and σ (2) ij n j = − The solvability condition yields the following homogenized differential equation for f (0, ) (X): where and S = S dξ . For localized modes, that is for U j,k → 0 as ξ 2 → ∞, the final term vanishes yielding (3.32) Hence, returning to dimensional form, the homogenized continuum has two characteristic wave speeds (c ij = μT ij /ρ) (3.33) corresponding to compressional and shear waves, which yields two light-lines (e.g. figure 2b). Moreover, the diagonal entries of the tensor T can be identified with the P-wave modulus (T 11 ) and the shear modulus (T 22 ); the off-diagonal terms of T vanish. Thus, the homogenized behaviour of the elastodynamic system is governed by two scalar functions (alternatively a single vector function with two components), as expected [27].

An infinite array of cylindrical voids
As an example, we apply the high-frequency asymptotic procedure to an infinite line of tractionfree cylindrical voids, as shown in figure 1. The material and geometrical parameters used are μ = 1, λ = 1.5 (hence ν = 1 3 ), ρ = 1, = 1, and the radius of the voids S 0 is 0.8. This choice of radius is made because it is very close to the radius at which the fourth mode appears (cf. figure 3b). There exist five RB-like localized modes below the compressional light-line; the dispersion curves for each mode are shown in figure 3a. Figure 3b shows the variation in standing wave frequencies with respect to increasing radius of the voids. The dashed line indicates the intersection of the shear light-line with the line k = π/2 where we expect to see standing waves. For the radius chosen (r = 0.8), figure 3b indicates that there are four RB-type modes underneath the shear light-line. For smaller radii, the number of modes decreases. It is interesting to note that the fourth RB mode emerges as we reach the chosen radius of r = 0.8; that is, r ≈ 0.8 is the critical radius at which a fourth RB mode emerges and descends from the shear light-line (cf. figure 3b). The dispersion curve is constrained to remain below the shear light-line and hence appears to be very short in figure 3a. Moreover, as the fourth mode lies very close to the light-line, it is relatively weakly localized. Indeed, figure 10 illustrates that the RB wave is shear dominated and the leading order fields resemble a bulk elastic wave. Although it may appear from figure 10b,c that the field does not decay in the x 2 direction it, in fact, does but very slowly. A similar effect is observed in the scalar case [20].
The two rays emanating from the origin have slopes as defined in (3.33) and represent the shear and compressional light-lines, which are obtained from the low frequency asymptotics as described in §3b. Moving to the end of the Brillouin zone (k = π/2), where the displacement fields in neighbouring cells are out of phase with each other, we find five localized modes at frequencies away from Ω = 0. As can be seen from figure 3a, the asymptotic procedure captures the behaviour of the dispersion curves very well in the vicinity of the standing wave frequencies.
The standing wave frequencies, together with the associated T-values, computed using the high-frequency asymptotic procedure as described in §3 are detailed in table 1. The T-value for the fourth mode is omitted as it emerges very close to the shear light-line. The T-values represent the stiffness of the homogenized material, which replaces the array of voids on the long-scale.
The implementation uses a bespoke algorithm developed for the commercial finite-element package Comsol Multiphysics in order to solve the required eigenvalue problems for U (0) (ξ ) and U (1) (ξ ) and, hence, obtain the T-values shown in table 1. It is worth mentioning that other appropriate methods do exist, which could be used to investigate the problem at hand. In particular, an alternative, semi-analytical method, would be to use multipole expansions about the centre of the circular voids in a similar approach to that used in scattering problems for elasticity (e.g. [25,26]). The use of multipole expansions may also lead to further simplification of the expressions for the T-values. Figures 5-11 show the leading-order displacement fields U (0) (ξ ), together, with the leadingorder shear ψ (0) (ξ ) and compressional ϕ (0) (ξ ) potentials. It is clear from figures 5-11 that, although the problem is fully coupled, in some cases the shear potential is larger than the compressional   shear potential. However, it should be noted that even in cases where χ differs significantly from unity, the problem remains fully coupled through the boundary conditions on the inclusion ∂S 0 and the edges of the elementary cell ∂S ± (via conversion from pressure to shear and vice versa). Moreover, the reader should bear in mind that the displacement field is obtained from linear combinations of the partial derivatives of the potentials, rather than the potentials themselves.
(a) The macroscopic behaviour where Ω 0 is the standing wave frequency, P.V. denotes the principal value, and the negative sign is chosen to prevent growth of the macroscopic field at infinity. Clearly, the long-scale behaviour depends on the sign of the argument of the square root: if (Ω 2  negative (resp. positive). It is emphasized that the beat frequency of the macroscale envelope function depends both, on the perturbation away from the standing wave frequency, and the microscale properties encoded on the dynamically homogenized stiffness T. Figures 7-13 illustrate the overall behaviour of the RB modes. The illustrative computations were produced using the commercial finite-element software Comsol Multiphysics and simulated RB waves, excited by a remote source, and propagating along an infinite array of traction-free cylindrical voids. Figures 7, 12 and 13 correspond to modes one, two and three, respectively; parts (a), (c) and (e) are the results of the numerical simulations, while parts (b), (c) and (f ) are obtained from the high-frequency homogenization scheme. It is clear from these figures, firstly that inplane elastic RB modes exist, and secondly that the high-frequency homogenization scheme captures their behaviour very well. Both, the rapid short-scale oscillations, and the macroscale envelope functions are accurately reproduced by the homogenization scheme.

Concluding remarks
The recent proposal of geometrically induced (or spoof) SPP with structured metal interfaces [7] has renewed the interest of physicists and mathematicians alike in RB waves. The spoof SPP concept was first realized experimentally on planar geometries [38,39], but it has been also extended to wedge diffraction gratings [40], where strong confinement of electromagnetic waves has been demonstrated in the terahertz regime. Potential applications of RB waves in plasmonics lie in photovoltaics via photo-thermal effects, but recent experiments on light and sound interplay through Brillouin scattering in arrays of fibres with defects [41,42] suggest that the existence of RB waves in elastodynamics might open new vistas in fast opto-elastic switches.
The most significant outcome of this paper is the demonstration that RB waves do exist for diffraction gratings in the setting of fully coupled (in-plane) vector elasticity. Despite significant interest for scalar problems (acoustics, flexural waves, polarized electromagnetic waves, etc.) and their numerous applications, this paper is the first study of RB waves for elastic diffraction gratings in the scientific literature. By couching this study in the language of homogenization, we were also able to develop an efficient asymptotic homogenization theory for RB waves. Indeed, owing to the large number of spurious modes generated by the numerical implementation of perfectly matched layers and infinite-element domains, there are significant computational challenges in extracting RB waves numerically. However, the behaviour of the RB waves can be conveniently obtained using the asymptotic procedure developed herein. Moreover, using the homogenization theory, long diffraction gratings which are computationally expensive to simulate (such as those illustrated in figures 7-13) can be replaced by an effective continuum as detailed in §3, whence the two-dimensional diffraction grating problem reduces to a onedimensional problem on a line-segment. Figures 7-13 are a testament to the efficacy of this approach.
Beyond efficient numerical modelling, one can envisage many practical applications of elastic RB waves. For example, the long-scale wavelength can be altered by the introduction of a perturbation or defect in the diffraction grating and this effect could be used to filter or guide waves of specific frequency: see figure 8, wherein the removal of two voids within the grating leads to a confinement of the RB wave in the defect so created. It should be noted that the length of the defect can be adjusted to support a larger number N of beats of the RB wave, i.e. one needs to remove one inclusion every 8N to 9N inclusions. Such a localization result, in conjunction with recent advances in light and heat confinement in phononic crystals [42], suggests a whole new range of applications with elastic gratings. From the point of view of non-destructive evaluation, one could infer properties of a diffraction grating (e.g. weldment or riveted joint) from the scattered waves; waves could even be guided along structured joints during ultrasonic testing. Similar problems have been considered, by Thompson and co-workers [43,44], for acoustical diffraction gratings. In [43], Linton et al. studied the excitation of RB waves along semi-infinite diffraction gratings consisting of sound-hard scatterers. In the later paper, Thompson & Linton [44] developed the earlier analysis to study scattering by clusters of defects in infinite arrays of Neumann scatterers; in this case, the defects were formed by removing a finite number of scatters from the array. There the authors found that the defects excited RB waves and, moreover, that the scattered RB waves could be trapped between two widely spaced scatterers.
More fundamentally, the existence of RB waves for elastic diffraction gratings raises interesting questions with regard to other wave phenomena, which exist in scalar problems but have yet to be studied in elasticity. Given its wide applicability (see [20,31,32,34,35,45] among others), the highfrequency homogenization procedure would seem to be an excellent tool with which to further examine possible correspondences between phenomena across a wide range of physical systems.