Plane Wave-Perturbative Method for Evaluating the Effective Speed of Sound in 1D Phononic Crystals

Flores Méndez, J.; Salazar Villanueva, M.; Ambrosio Lázaro, R. C.; Calixto Sirene, B.; Mota González, M. L.; Candia García, F.

doi:https://doi.org/10.1155/2016/3017835

Advances in Materials Science and Engineering

On this page

Abstract Introduction Discussion Conclusions Acknowledgments References Copyright Related Articles

Special Issue

Metamaterials and Smart Structures in a Big Data Era

View this Special Issue

Research Article | Open Access

Volume 2016 | Article ID 3017835 | https://doi.org/10.1155/2016/3017835

Plane Wave-Perturbative Method for Evaluating the Effective Speed of Sound in 1D Phononic Crystals

J. Flores Méndez,^1,2M. Salazar Villanueva,¹R. C. Ambrosio Lázaro,¹B. Calixto Sirene,¹M. L. Mota González,³and F. Candia García¹

Academic Editor: Tiejun Liu

Received06 Jul 2016

Revised01 Oct 2016

Accepted03 Nov 2016

Published29 Nov 2016

Abstract

A method for calculating the effective sound velocities for a 1D phononic crystal is presented; it is valid when the lattice constant is much smaller than the acoustic wave length; therefore, the periodic medium could be regarded as a homogeneous one. The method is based on the expansion of the displacements field into plane waves, satisfying the Bloch theorem. The expansion allows us to obtain a wave equation for the amplitude of the macroscopic displacements field. From the form of this equation we identify the effective parameters, namely, the effective sound velocities for the transverse and longitudinal macroscopic displacements in the homogenized 1D phononic crystal. As a result, the explicit expressions for the effective sound velocities in terms of the parameters of isotropic inclusions in the unit cell are obtained: mass density and elastic moduli. These expressions are used for studying the dependence of the effective, transverse and longitudinal, sound velocities for a binary 1D phononic crystal upon the inclusion filling fraction. A particular case is presented for 1D phononic crystals composed of W-Al and Polyethylene-Si, extending for a case solid-fluid.

1. Introduction

At present there is a great interest in fabricating artificial materials, having extraordinary properties, which considerably extend those of natural materials. Such a new class of materials is known as metamaterials. From the beginning, the photonic metamaterials, possessing negative index of refraction, have intensively been investigated. It was established that such an unusual optical property is found in periodic structures (photonic crystals) whose dielectric function is spatially modulated. A peculiarity of the photonic metamaterials is the high dielectric contrast between the components. The negative refraction has been observed in double negative metamaterials with simultaneously negative effective permittivity and permeability [1–6]. However, this phenomenon is also manifested in anisotropic materials, as simple as a one-dimensional (1D) photonic crystal or superlattice, having elements in the effective permittivity tensor of different sign [7–10]. Analogously, metamaterials with uncommon acoustic properties have also searched among phononic crystals (PCs), that is, materials with periodic modulation of their elastic properties. As in the photonic case, the calculation of the effective parameters of a phononic crystal is an important task, for designing of resonant elastic metamaterials and acoustic lenses; there is a great interest in investigating the effective mass density and compliance tensors, as well as the effective sound velocity.

Several homogenization theories, which are valid when the acoustic wavelength is much longer than the lattice constant of the phononic crystal, have been proposed [11–15]. Among homogenization theories, we can identify two commonly used approaches. One of them provides effective acoustic parameters within the framework of multiple scattering [11–13] in which the effective parameters are not obtained according to the direction of propagation; for example, for a sonic crystal the solid material is simply considered rigid; the study is done on the propagation of sound in the air of the structure’s matrix, neglecting the modes that can propagate in the solid. The second, which is based on the Fourier formalism, makes use of the expansion of the microscopic acoustic field into plane waves and allows calculating the effective parameters as a function of the direction of propagation [14, 15]. The latter approach is of particular interest because it can be applied to different geometries of the inclusions inside the unit cell. Here, we will precisely use the Fourier formalism to calculate the effective velocities of elastic wave propagation in 1D solid-solid and fluid-solid phononic crystals. The very first known experimental observation of 1D phononic crystals was when Narayanamurti et al. investigated the propagation of high-frequency phonons through a GaAs/AlGaAs superlattice; other domains in which 1D phononic crystals have potential applications are crystal sensors, waveguide, and resonant transmission [16–19]. As in the case of 1D photonic crystals, such an inherently anisotropic elastic system is a potential metamaterial. However, the main goal of the present work is to establish the basis for the development of a general homogenization theory of (1D, 2D, or 3D) elastic phononic crystals based on the plane wave expansion.

The work is divided as follows: in Section 2, the Fourier formalism applied to a binary elastic superlattice is presented. The proposed method for calculating the effective acoustic parameters for a phononic superlattice, which is based upon the derivation of the wave equation for the macroscopic displacements field, is described in Section 3. Finally, we shall apply the derived explicit formulas for the effective sound velocities to binary superlattices composed of W-Al and Polyethylene-Si and thus study the behavior of this parameter effective in crystals with high and low contrast between the mechanical properties of its components; also a case solid-fluid is presented.

2. One-Dimensional Phononic Crystals

Let us consider a one-dimensional phononic crystal (or elastic superlattice) composed of alternating layers of isotropic elastic materials, A and B (Figure 1). Their thicknesses are, respectively, and , where is the lattice constant. Assuming that the -axis is parallel to the superlattice growth direction, the mass density , as well as the longitudinal and transverse sound velocities, and , turn out to be functions of the -coordinate only.

In this case, the second Newton’s law for the displacement vector , in the absence of external forces, acquires the form [20, 21]:Here, the Cauchy stress tensor is expressed in terms of the strain tensoraccording to the relationwhere is Kronecker’s delta. Substituting (2) and (3) into (1), the wave equation for an elastic inhomogeneous medium is obtained [22, 23]:

Due to the periodicity of the functions , , and , we can expand them into Fourier series aswhere and are the longitudinal and transversal elastic constants for an isotropic elastic medium, and the components for the vectors of the reciprocal lattice are expressed asThe component for an isotropic medium is given by the relation:

According to our geometry (see Figures 1 and 2), the Fourier coefficients , appearing in (5), can be calculated with the formulawhere is the form factor of the inclusion A,and the ratio is its filling fraction. Similarly, the Fourier coefficients and of the expansions (6) and (7) can, respectively, be written as

For a periodic system, the solution of (4) should satisfy the Bloch theorem [24]. Therefore, we can write the displacement vector in the formwhere is a periodic function with the same period as for the superlattice and is the frequency. We should note that the component is here the Bloch wave number and is an independent vector parallel to the plane. After expanding the function into Fourier series, an expression can be obtained It is worth emphasizing that the Fourier formalism for solving the wave equation (4) allows us to avoid the application of boundary conditions for the displacement vector , since the structural information of the system is explicitly contained in the Fourier coefficients . This is owing to the fact that the boundary conditions for are consistent with the wave equation (4).

Substituting (5)–(7) and (14) into (4), we obtainWe can rewrite (15) in dyadic form aswhereand is the unit dyadic .

The inverse υ of matrix in Eq. (17) satisfies the relationand its elements are directly calculated by using the expressionMultiplying (17) by (19), we getwhere we have introduced the dyadic :The system (20) represents a homogeneous system of algebraic equations, which has a nontrivial solution if the determinant of the associated matrix is identical to zero; that is,

The condition (22) provides the dispersion relation for vibrational eigenmodes in the PC. It should be noted that such a phononic dispersion parametrically depends upon the wave vector component . Thus, the phononic dispersion can be straightforwardly calculated by numerically diagonalizing the matrix . In doing that, we should use a finite matrix of sufficiently large size to guarantee accurate results.

3. Effective Sound Velocity

In the case when the wave length of sound is much larger than the lattice constant of the phononic crystal, this artificial periodic heterostructure can be modeled as a homogeneous medium with effective acoustic parameters. In this section, we shall calculate the effective sound velocities for both transverse and longitudinal vibrational modes, propagating along the growth direction of a superlattice like that considered in the previous sections.

From the Bloch theorem, given by (14), we can write the displacements field asAssuming and , the first term on the right-hand side of (23) describes the smooth variations of the displacement vector as a function of the coordinates, whereas the second term therein is associated with the rapid oscillations (fluctuations) over the unit cell. So, we can define the average (macroscopic) displacements field asAccording to (20), the Fourier coefficients satisfy the system of equationswhere we have introduced the dyadicOne can express the coefficients in terms of by using (25) for . Thus,Here, is a submatrix, obtained from (26) after eliminating its block rows (columns) with . As it was demonstrated in [25], the inverse of the submatrix is related with the inverse of the matrix according to the formula:where is a block, which is obtained from the original matrix (26) of infinite size, and symbolizes the inverse of the matrix block.

Let us rewrite (25) for asSubstituting (27) and (28) into (29), we obtain the macroscopic wave equation:where the effective matrix is given byFrom the macroscopic wave equation (30), we can define an effective dyadic representing the Christoffel tensor, divided by the mass density, asThus,In the situation, considered here, of acoustic waves propagating along the growth direction, the effective tensor (33) turns out to be diagonal; that is, with . The principal values and are, in fact, the effective sound velocities for transverse and longitudinal modes, respectively. It should be emphasized that formula (33) is valid in the long-wavelength limit, that is, when both and .

4. Numerical Results and Validation

Let us apply the derived formula (33) for calculating the effective sound velocities of modes propagating along the growth direction of specific binary 1D phononic crystals. First, we shall consider a Polyethylene-silicon superlattice having a period cm. The parameters used in the calculations are mass densities gr/cm³ and gr/cm³, transverse sound velocities cm/s and cm/s, and longitudinal sound velocities cm/s and cm/s for silicon and Polyethylene, respectively (the material parameters were taken from [24]). In the numerical calculations, 300 plane waves were needed to achieve good convergence of the results.

Figures 3 and 4 (continuous line) show the dependencies of the transverse and longitudinal effective sound velocities upon the silicon filling fraction . As it is seen, both transverse and longitudinal effective sound velocities slowly vary with increasing the silicon filling fraction up to . In contrast, at , the slopes of the curves and are relatively large, because of the large contrast between the acoustic impedance and () of the system.

Figures 5 and 6 (continuous line) exhibit numerically calculated effective parameters, and , for a W-Al 1D phononic crystal. The parameters used here are [24] gr/cm³ and gr/cm³, transverse sound velocities cm/s and cm/s, and longitudinal sound velocities cm/s and cm/s for aluminium and tungsten, respectively. In this case, the effective sound velocity for transverse [longitudinal] modes decreases with the Al filling fraction up to its minimum value at . For larger values of (i.e., in the interval ) the effective parameter [] increases with .

In order to verify our numerical results in Figures 3–6 (continuous line), we use the effective nonlocal-response matrix in the quasistatic limit (). The block structure of the matrix for such systems has the form [26]where and are the effective mass density and compliance tensors within the long-wavelength limit.

Figure 7 exhibits graphs of the nonzero elements in the matrices (only principal value ) and (only the effective stiffness constants and ) in (34) for a W/Al and Polyethylene/Si 1D phononic crystals versus the inclusion filling fraction . The principal -axis of the 1D phononic crystal has been oriented parallel to periodicity direction. Hence, the homogenized phononic crystal acquires tetragonal symmetry in the interval , for which the effective sound velocity for transverse and longitudinal vibrational modes in the low frequency limit propagating along the growth direction is given according to the formulas (see [27], where Christoffel’s equation for the tetragonal system can be factored to obtain the solutions for the transverse and longitudinal acoustic velocity):Applying formulas (35) and (36), the effective sound velocity for transverse and longitudinal vibrational modes is shown in Figures 3–6 (squares). We have found that our numerical results given by formula (33) coincide with those predicted by equations (35) and (36). The effective nonlocal-response matrix method is based on matrix inversion of size (: number of plane waves); the number of plane waves increases as the contrast between the mechanical properties of the components of phononic crystal is large and this represents a great disadvantage for the numerical convergence of the method. Formula (33) shows that the sums over and converge more rapidly than the matrix inversion with finite number of plane waves; therefore numerical calculation of the sum is not a time consuming procedure at all.

(a)

(b)

(c)

Otherwise, formula (33) can be applied for determining the effective sound velocity of 1D-phononic crystals with a type fluid component, which has zero shear modulus; this is possible by using a very small value of () in order to avoid numerical errors. With this mathematical artifice, we have calculated the effective sound velocity for the longitudinal vibrational mode for a water-aluminium superlattice having a period cm (see Figure 8). The parameters used for water are gr/cm³ and longitudinal sound velocity cm/s. The graph shows the dependence of the effective longitudinal sound velocity upon the aluminium filling fraction ; note that until filling fractions near 0.8 the values of the effective velocity remain close to the water, for later ascending quickly until reaching the values of the solid, thanks to the large mismatches at the interfaces of the water-aluminium building-block (due to high contrast of the acoustic impedance of the materials in the unit cell).

5. Discussion

In the developed method in this work, the homogenization is achieved in the limit , and the effective Christoffel tensor is calculated from the wave equation for an unidimensional elastic inhomogeneous medium; this tensor relates the elastic constants with the propagation velocities of the elastic waves according to different directions of the principal axes that are entirely dependent on the normal wave on the direction of but not its magnitude, so that the phase velocity is defined by the square roots of the eigenvalues of the effective matrix . We solve the determinant of the effective matrix taking into account the solution method proposed by [27]. It should be noted that [27] does not present a methodology of homogenization for different layers of materials but only develops the general theory of wave propagation. In this way, we obtain the relationship between the propagation velocity in the direction and the elastic and mechanical homogenized properties of the PC. Therefore, our formula simultaneously calculates the effective velocities to which the transverse and longitudinal elastic waves propagate in the PC. On the other hand, when the 1D-PC has fluid type components, close to zero values are ingressed for (transversal elastic constant) to calculate the effective sound velocity for the longitudinal vibrational mode. This allows establishing a more general homogenization methodology for calculating the effective velocity of propagation of sound waves in phononic crystals with elastic constituents of solid and fluid type, which is useful to define a methodology towards the development of a theory for two- and three-dimensional systems. In contrast to the methodology presented in [14, 15] which presents a separate analysis to study the propagation of acoustic waves in fluid-fluid and solid-solid phononic crystals, in such studies, the first case is based on solving the equation sound wave propagating in a gas for the longitudinal modes; in the second case the equation of motion for an elastic medium is solved considering only transverse modes. The effective sound velocity in the homogenized system is obtained by taking the limit ω, ; in this limit the group velocity is equal to the phase velocity, which is calculated as . With the above considerations, the authors obtain for each case the formula for the effective sound velocity; also they do not show results for a solid-liquid PC.

6. Conclusions

We have derived explicit formulas for the calculation of the effective sound velocities in a 1D phononic crystal in the long-wavelength limit. The formulas were applied for analyzing the dependence of the effective, transverse and longitudinal, sound velocities upon the inclusion filling fraction for binary superlattices composed of Polyethylene-Si and W-Al. In the latter case, the contrast of material parameters is relatively larger and, as a result, at Al filling fractions , the effective sound velocity for both transverse () and longitudinal () modes, propagating along the superlattice growth direction, takes values smaller than the sound velocity for each component (W or Al). For this reason, we can say that the homogenized 1D W-Al phononic crystal behaves as a metamaterial. We have verified that our results in the quasistatic limit for Polyethylene-Si and W-Al 1D phononic crystals coincide with the results predicted by the effective mass density and compliance tensors method. Finally, formula (33) can also be applied for determining the effective sound velocity of 1D-phononic crystals with a type fluid component having zero shear modulus.

Although the homogenization theory developed here is valid only for phononic crystals with one-dimensional periodicity and isotropic inclusions, it shows the usefulness of the plane wave expansion method to obtain explicit expressions for theoretical results of the effective sound velocity. The generalization of this approach to 2D and 3D periodic elastic structures with anisotropic inclusions and arbitrary contrast of the materials parameters will be presented elsewhere.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was partially supported by VIEP-BUAP (Grant FOMJ-ING16-I). The authors also thank the support of PRODEP program (Grant BUAP-PTC-384, Agreement no. DSA/103.5/15/7449).

References

J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Transactions on Microwave Theory and Techniques, vol. 47, no. 11, pp. 2075–2084, 1999.
View at: Publisher Site | Google Scholar
D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Physical Review Letters, vol. 84, no. 18, pp. 4184–4187, 2000.
View at: Publisher Site | Google Scholar
L. Zhang, G. Tuttle, and C. M. Soukoulis, “GHz magnetic response of split ring resonators,” Photonics and Nanostructures—Fundamentals and Applications, vol. 2, no. 2, pp. 155–159, 2004.
View at: Publisher Site | Google Scholar
C. Enkrich, M. Wegener, S. Linden et al., “Magnetic metamaterials at telecommunication and visible frequencies,” Physical Review Letters, vol. 95, no. 20, Article ID 203901, 2005.
View at: Publisher Site | Google Scholar
S. Zhang, W. Fan, N. C. Panoiu, K. J. Malloy, R. M. Osgood, and S. R. J. Brueck, “Experimental demonstration of near-infrared negative-index metamaterials,” Physical Review Letters, vol. 95, no. 13, Article ID 137404, pp. 1–4, 2005.
View at: Publisher Site | Google Scholar
J. A. Reyes-Avendaño, U. Algredo-Badillo, P. Halevi, and F. Pérez-Rodríguez, “From photonic crystals to metamaterials: the bianisotropic response,” New Journal of Physics, vol. 13, Article ID 073041, pp. 1–33, 2011.
View at: Publisher Site | Google Scholar
V. A. Podolskiy and E. E. Narimanov, “Strongly anisotropic waveguide as a nonmagnetic left-handed system,” Physical Review B—Condensed Matter and Materials Physics, vol. 71, no. 20, Article ID 201101, 2005.
View at: Publisher Site | Google Scholar
B. Wood, J. B. Pendry, and D. P. Tsai, “Directed subwavelength imaging using a layered metal-dielectric system,” Physical Review B, vol. 74, no. 11, Article ID 115116, pp. 1–8, 2006.
View at: Publisher Site | Google Scholar
A. L. Rakhmanov, V. A. Yampol'Skii, J. A. Fan, F. Capasso, and F. Nori, “Layered superconductors as negative-refractive-index metamaterials,” Physical Review B, vol. 81, no. 7, Article ID 075101, 2010.
View at: Publisher Site | Google Scholar
B. Zenteno-Mateo, V. Cerdáan-Ramìrez, B. Flores-Desirena, M. P. Sampedro, E. Juáarez-Ruiz, and F. Páerez-Rodrìguez, “Effective permittivity tensor for a metal-dielectric superlattice,” Progress in Electromagnetics Research Letters, vol. 22, pp. 165–174, 2011.
View at: Publisher Site | Google Scholar
J. Mei, Zh. Liu, W. Wen, and P. Sheng, “Effective mass density of fluid-solid composites,” Physical Review Letters, vol. 96, no. 2, Article ID 024301, pp. 1–4, 2006.
View at: Publisher Site | Google Scholar
D. Torrent, A. Håkansson, F. Cervera, and J. Sánchez-Dehesa, “Homogenization of two-dimensional clusters of rigid rods in air,” Physical Review Letters, vol. 96, no. 20, Article ID 204302, 2006.
View at: Publisher Site | Google Scholar
D. Torrent and J. Sánchez-Dehesa, “Effective parameters of clusters of cylinders embedded in a nonviscous fluid or gas,” Physical Review B—Condensed Matter and Materials Physics, vol. 74, no. 22, Article ID 224305, 2006.
View at: Publisher Site | Google Scholar
A. A. Krokhin, J. Arriaga, and L. N. Gumen, “Speed of sound in periodic elastic composites,” Physical Review Letters, vol. 91, no. 26, Article ID 264302, pp. 1–4, 2003.
View at: Google Scholar
L. N. Gumen, J. Arriaga, and A. A. Krokhin, “Metafluid with anisotropic dynamic mass,” Low Temperature Physics, vol. 37, no. 11, pp. 1221–1224, 2011.
View at: Google Scholar
V. Narayanamurti, H. L. Störmer, M. A. Chin, A. C. Gossard, and W. Wiegmann, “Selective transmission of high-frequency phonons by a superlattice: the “dielectric” phonon filter,” Physical Review Letters, vol. 43, no. 27, pp. 2012–2016, 1979.
View at: Publisher Site | Google Scholar
R. Lucklum, J. Li, and M. Zubtsov, “1D and 2D phononic crystals sensors,” Procedia Engineering, no. 5, pp. 436–439, 2010.
View at: Google Scholar
M. Xiao, G. Ma, Z. Yang, P. Sheng, Z. Q. Zhang, and C. T. Chan, “Geometric phase and band inversion in periodic acoustic systems,” Nature Physics, vol. 11, no. 3, pp. 240–244, 2015.
View at: Publisher Site | Google Scholar
H. Sanchis-Alepuz, Y. A. Kosevich, and J. Sánchez-Dehesa, “Acoustic analogue of electronic bloch oscillations and resonant zener tunneling in ultrasonic superlattices,” Physical Review Letters, vol. 98, no. 13, Article ID 134301, 2007.
View at: Publisher Site | Google Scholar
L. D. Landau and E. M. Lifshitz, Theory of Elasticity, Butterworth Heinemann, New York, NY, USA, 1995.
L. M. Brekhovskikh and O. A. Godin, Acoustics of Layered Media I: Plane and Quasi-plane Waves, Springer, Berlin, Germany, 1998.
M. S. Kushwaha, P. Halevi, G. Martínez, L. Dobrzynski, and B. Djafari-Rouhani, “Theory of acoustic band structure of periodic elastic composites,” Physical Review B, vol. 49, no. 4, pp. 2313–2322, 1994.
View at: Publisher Site | Google Scholar
M. S. Kushwaha, “Classical band structure of periodic elastic composites,” International Journal of Modern Physics B, vol. 10, no. 9, pp. 977–1094, 1996.
View at: Publisher Site | Google Scholar
Ch. Kittel, Introduction to Solid State Physics, John Wiley & Sons, New York, NY, USA, 1995.
V. Cerdán-Ramírez, B. Zenteno-Mateo, M. P. Sampedro, M. A. Palomino-Ovando, B. Flores-Desirena, and F. Pérez-Rodríguez, “Anisotropy effects in homogenized magnetodielectric photonic crystals,” Journal of Applied Physics, vol. 106, no. 10, Article ID 103520, 2009.
View at: Publisher Site | Google Scholar
J. Flores-Méndez and F. Pérez-Rodríguez, “Metasolid with anisotropic mass density,” Europhysics Letters, vol. 103, no. 5, Article ID 54001, 2013.
View at: Publisher Site | Google Scholar
B. M. Lempriere, Ultrasound and Elastic Waves, Academic Press, New York, NY, USA, 2002.

Copyright

Copyright © 2016 J. Flores Méndez et al. 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.

PDF Download Citation

Download other formats

Order printed copies

Views

1376

Downloads

994

Citations