Open access peer-reviewed chapter

Phononic Crystals and Thermal Effects

By Arafa H. Aly and Ahmed Mehaney

Submitted: August 23rd 2018Reviewed: October 16th 2018Published: April 15th 2019

DOI: 10.5772/intechopen.82068

Downloaded: 172

Abstract

In this work, we demonstrate a comprehensive theoretical study of one-dimensional perfect and defect phononic crystals. In our study, we investigate the elastic and shear waves with the influences of thermal effects. The numerical calculations based on the transfer matrix method (TMM) and Bloch theory are presented, where the TMM is obtained by applying the continuity conditions between two consecutive sub-cells. Also, we show that by introducing a defect layer in the perfect periodic structures (defect phononic crystals), we obtain localization modes within the band structure. These localized modes can be implemented in many applications such as impedance matching, collimation, and focusing in acoustic imaging applications. Then, we investigate the influences of the incident angle and material types on the number and intensity of the localized modes in both cases of perfect/defect crystals. In addition, we have observed that the temperature has a great effect on the wave localization phenomena in phononic band gap structures. Such effects can change the thermal properties of the PnCs structure such as thermal conductivity, and it can also control the thermal emission, which is contributed by phonons in many engineering structures.

Keywords

  • thermal emission
  • dispersion relation
  • phononic band gap
  • localized modes

1. Introduction

Phononic crystals (PnCs) are new composite materials which can interact, manipulate, trap, prohibit, and transmit the propagation of mechanical waves. Recently, great efforts have been dedicated to study these novel materials to be used in many potential engineering applications. By reference to the meaning of idiomatic, the term “phononic” was derived in analogy to the term “phonon,” which is considered as a quantization of the lattice vibrations. In the previous conscious of science, these vibrations are impossible to be controlled because the atoms in the solid cannot move independently of each other because they are connected by chemical bonds. By appearing PnCs, this assumption was changed and the mechanical wave can be filtered, transmitted, stopped, and localized within specific frequencies called phononic band gaps [1, 2].

Within the phononic band gaps, the mechanical waves of all types are greatly blocked. Actually, the formation of the phononic band gaps is back to the variation in the mechanical properties of the materials that build the PnC structure. Therefore, Bragg interference at the interface between each two materials can be obtained. As a result, the phononic band gaps and transmission bands are formed. Moreover, novel properties such as negative refraction and acoustic metamaterials are presented in PnC structures [3, 4, 5, 6, 7]. These novel properties of PnCs can be utilized in many industrial and engineering applications such as MEMS applications, filters, waveguides, clocking, multiplexers, and sensor applications [8, 9, 10, 11, 12].

As it is well-known, the physical origin of all crystal structures has the same idea of the design. Photonic crystals and semiconductors devices could not develop such types of the previous applications. From the scientific point of view, we should mention that the physical nature of PnCs is different from one of photonic crystals, as well as semiconductors. Generally, the various forms of waves are referred to as: electron waves as scalar waves and optical waves as vector waves, while elastic waves as tensor waves [13, 14, 15, 16, 17, 18, 19].

As a result, to design PnCs having a complete phononic band gap, the mechanical properties must be changed not only in one direction (i.e. x- direction) but also in all three directions of space (i.e., x, y, and z directions). Therefore, PnCs can be classified according to the periodicity into three types, i.e., the one-dimensional (1D), the two-dimensional (2D), and the three-dimensional (3D) PnCs [20, 21, 22, 23, 24]. PnCs can control the entire spectrum of phonons frequency from 1 Hz to THz range. As it is well-known, the waves that propagate through solid media are called elastic waves, while those that propagate through fluids are called acoustic waves. Also, PnCs can control the different types of mechanical waves such as elastic waves, acoustic waves, and surface waves. Therefore, unlike photonic crystals, not all 3D PnCs possess complete phononic band gaps. PnCs must possess band gaps for both elastic and acoustic waves at the same frequency region. Consequently, we have to fabricate PnC structures for both solid and fluid media [25, 26, 27, 28]. Hence after, these different polarizations introduce more challenges to PnCs higher than other crystal, which in turn makes the theoretical manipulation of PnCs more attractive and perspective.

From the previous observations of PnCs structures, several methods were developed to calculate the phononic band gaps such as Plane-Wave-expansion Method (PWM), Bloch-Floquet Method (BFM), and finite different time-domain method (FDTD) [29, 30, 31, 32, 33]. In this chapter, we will depend on the Transfer Matrix Method (TMM) [34] for calculating the reflection coefficient and the transmission coefficient of the one-dimensional PnC structure. Such method considered very suitably for the 1D structure due to its recursive nature, since it allows the continuity of the waveform at the interface between each two layers [35, 36]. Also, by using the TMM, we can obtain the dispersion relations of the mechanical waves through the periodic PnC structures. PnCs and temperature are in mutual influences in several ways [37]. For example, at any temperature, mainly the huge contribution in the thermal conductivity of many materials is dominant by phonon contribution, which is a function of phonon mean free path and the Boltzmann distribution of phonons of any material. Also, the thermal conductivity is depending on the thickness of a material, where optical branches contribute with about 30% of the thermal conductivity of any material. By inhibiting the acoustic phonon population, the optical phonon relaxation is indirectly inhibited by up to 30% and hence limits their contribution to thermal conductivity. Furthermore, in silicon PnC, the thermal contribution of phonons has been reduced to less than 4% of the value for bulk silicon at room temperature [38, 39, 40].

Moreover, the defected PnCs structures have wonderful application in wave guiding and multiplexing. But the produced structures are different from the proposed ideal structures due to the errors and defaults in manufacturing. By removing some layers or materials from the ordered periodic structure in a PnC, we can create a point or a waveguide defects that are able to localize and bend signals. Such mechanism has been implemented similarly in photonic crystals to slow light. Consequently, the produced waveguides in PnCs structures can be used in focusing and collimation of acoustic waves in medical ultrasound applications, sensors, and MEMS applications [2, 41, 42, 43, 44].

In the present work, we introduce the formation of phononic band gaps under the influences of many physical parameters. First, the general case of SH-wave propagation in arbitrary direction will be investigated by using TMM with calculation of the dispersion curves as well. Also, the influences of the incident angle on the band gaps are analyzed and discussed. In addition, this work focuses on correlating and comparing the results of SH-waves with in-plane waves propagating normally to the structure. For the in-plane waves, the reflection coefficients for S- and P-waves are plotted and compared with the dispersion relations curves. Furthermore, we are demonstrating the wave localization phenomenon in PnCs and the effects of the temperature on the band structure of PnCs and the localized modes for both in-plane and SH-waves. Also, the numerical results are presented and discussed to investigate the effect of the defect layer on the wave localization modes inside the structures. Finally, the effects of the thickness and type of the defect layer material on the band gap structure had been discussed.

2. Theoretical analysis and numerical models

2.1 Equation of wave propagation

Figure 1 shows the schematic diagram of the 1D PnC crystal structure. The proposed crystal structure has an infinite number of the periodically arranged unit cells. The unit cell may include two or more layers; here, we propose that it is made by only two materials A and B, respectively. The two materials are labeled by the subscript j=1, 2. Also, the thickness of the unit cell (the lattice constant) is a=a1+a2.The thickness, Lame’ constant, shearing modulus, Poisson’s ratios, mass density, and Young’s modulus of the two layers are denoted by equation aj,λj,μj,νj,ρj,Ej[Ej=λj1+νj12νj/νj, respectively.

Figure 1.

A schematic diagram of a perfect 1D binary PnC structure.

The governing equation of anti-plane shear waves (SH-waves) polarized in the z-direction propagated in the xy-plane can be written in the following form [45, 46, 47].

μj2φjxjyjt+σjtx2φj2xjyjtx2=ρjφ¨jxjyjtj=12,E1

where Tis the temperature variation, βjis the thermal expansion coefficients, φjxjyjtis the displacement components along the z-direction, tis the time, σjtx=EjβjT/12υjis the thermal stress, and 2=/xj2+/yj2is the Laplacian operator. The solution φjxjyjtin the jth layer with time harmonic dependence can be expressed as [48],

φjxjyjt=ϕjxjexpikyjsinθ0t,E2

Where, cis the velocity component of the incident wave, i2=1, k=ω/cis the wave number, θ0is the angle of the incident wave, ωis the angular frequency, and ϕjis the amplitude.

In many cases, it is more convenient to represent the thickness value by the following dimensionless coordinates:

ξj=xja1¯,ηj=yja1¯j=12,E3

where a¯1is the average thickness value of the material A. a¯1is equal to a1for the periodic structures. By inserting Eq. (3) in Eqs. (1) and (2), we can obtain the following wave equations:

φjξjηjt=ϕjξjexpiαηjsinθ0iωt,E4
1+χjd2ϕjdξj2α2sin2θ0αj2α2ϕj=0,E5

where cjrepresents the wave velocity in each material, χj=σjtx/μjis the stress and shearing modulus ratio, α=ka¯1represents the SH-waves dimensionless wave number, αj=kja¯1is the dimensionless wave number and kj=ω/cjis the wave vector of materials Aand B, respectively. Eq. (5) has a general dimensionless solution, and it is given in the following form,

ϕjξj=Ajexpqjξj+Bjexp+qjξj,j=12,E6

where Ajand Bjare unknown coefficient to be determined, and qjis a parameter and has the value qj=11+χjc2/cj2sin2θ0. Therefore, the dimensionless solution φjξjηjtwith time harmonic dependence is:

φjξjηjt=Ajexpqjξj+Bjexp+qjξjexpiαηjsinθ0iωt.E7

2.2 Transfer matrix method

The dimensionless shear stress component is given as follows:

τxzj=μjφja¯1ξjj=12.E8

Assuming that the PnC consists of nunit cells, the boundary conditions at the left and right sides of the two layers in the ith unit cell can be written in the following form [45],

ϕjLi=ϕji0,ϕjRi=ϕjiζj,τxzjLi=μjiϕjia¯1ξj0,τxzjRi=μjiϕjia¯1ξjζji=12.n,E9

where the subscripts L and R denote the left and right sides of the two layers, and 0ξjζj=aj/a¯1j=12are the dimensionless thicknesses of materials Aand B. Substituting Eqs. (7) and (8) into Eq. (9), the following matrix equation can be obtained as follows:

νjRi=Tj'νjLij=12i=12n,E10

where νjRi=ϕjRia¯1τxzjRiTandνjLi=ϕjLia¯1τxzjLiTare the dimensionless state wave vectors at right and left sides of each unit cell and Tj'is the 2 × 2 transfer matrix of each unit cell. The four elements of Tj'can be written in the following forms,

Tj'11=expqjζj+expqjζj2,Tj'12=expqjζjexpqjζj2iαqμjj,Tj'21=iαqμjjexpqjζjexpqjζj2,Tj'22=Tj'11.E11

2.3 Characteristic of the dispersion relation

At the interface between the layers, the following condition is satisfied:

ν1Ri=ν2Li.E12

Thus, the relationship between the right and left sides of the ith unit cell can be obtained from Eq. (10) as follows:

ν2Ri=Tiν1Lii=12n,E13

where Tiis the accumulative transfer matrix of the ith unit cell and can be written in the following form:

Ti=T2'T1'.E14

At the interface between the right side of the unit cell and the left side of the ith unit cell, the following condition is satisfied:

ν1Li=ν2Ri1i=2..n.E15

By equating Eqs. (13) and (15), we can obtain the relationship between the state vectors of the i1th unit cells and the ith unit cell in the following form:

ν2Ri=Tiν2Ri1i=2..n,E16

Therefore, Tirepresents the transfer matrix between each two successive unit cells.

By using Floquet and Bloch theories, we can obtain the displacement and stress fields between each two neighboring unit cells (at the interface) by the following relations [46, 47, 48, 49]:

ϕ2Ri=ϕ2Ri1expika,τxz2Ri=τxz2Ri1expika,E17

where kis the wave number representing the direction of the traveled wave across the structure. Combining the above two equations leads to the following matrix equation:

ν2Ri=ν2Ri1expikai=2..n.E18

By equating Eqs. (16) and (18) which leads to the following eigenvalue problem:

TieikaI=0.E19

Eq. (19) can be rewritten as follows:

Tiν2Ri1=λν2Ri1,E20

where λ=eikais a complex eigenvalue and V2Ri1is a complex eigenvector.

2.4 Band structure formation

The wave number in Eqs. (17) and (18) is a complex number, so it can be a positive or negative number; in general, the wave number can be written in the following form [50]:

k=krealikimaginary,E21

where krealandkimaginaryare the real and imaginary wave numbers, respectively. Therefore, we can deduce from this relation that the complex wave number has two forms, so it can inhibit the incident waves within the phononic band gaps. Consequently, we have two frequency ranges and they are organized as follows:

  1. If k=krealandkreal>0.

From Eq. (18),

v2Ri=v2Ri1eika,v2Ri=v2Ri1eikreala.E22

From Eq. (22), we can deduce that the displacement and stress at the successive unit cells ithandi1thare differ only by a phase factor eikreala. At this condition, the elastic waves are allowed to propagate freely through the PnC structure with the corresponding frequencies and wave number, forming the so-called transmission bands.

  1. If k=ikimaginaryandkimaginary<0

From Eq. (18),

v2Ri=v2Ri1eika,v2Ri=v2Ri1ekimaginarya.E23

In contrast to the above case, we can deduce from Eq. (23) that the displacement and stress at successive unit cells ithandi1thare the same and do not have a phase difference. Moreover, the incident wave has an exponential spatial attenuation of strength magnitude proportional to the value kimaginary. Hence after, at this case, the waves are prohibited from propagation in the PnC structure, which in turn, results in formation of the so-called phononic band gaps or stop bands.

2.5 Temperature influences on PnCs

When a mechanical wave propagates through a PnC structure, resultant vibrations occur which may increase or decrease phonons motions. Therefore, it can heat or cool the PnC structure. Also, the thermal expansion of the constituents materials may change, which will affect the mechanical constants of the material as well. Since the vibrations occur very rapidly, there is no big chance to thermal energies to flow and the elastic constants measured by elastic waves propagation are changed adiabatically. The elastic constants are connected to the isothermal constants by the following relation [51, 52, 53, 54, 55],

λσ=λθ+9β2B2θρCν,μσ=μθ,E24

where the superscripts σand θindicate adiabatic and isothermal constants, βis the thermal expansion coefficient, B is the bulk modulusB=λ+2/3μ, Cνis the specific heat at constant volume, θis the absolute temperature, and ρis the mass density. From Eq. 24, we can deduce that λσandλθare not the same and the difference between them should be considered.

In addition to the above relations, if the temperature is increased, not only the wave velocities will increase but also the thickness of each layer will change by the following relation [53]:

Δa=βaiΔt,E25

where ais the thickness difference of any layer, aiis the original layer thickness, and Δtis the temperature difference. Hence after, these two variables will affect the longitudinal waves speed and the stress component in SH-waves equations. Since the P-wave velocity is cP=λ+2μρ, which, in turn leads to the variation of the band structure and band gaps properties.

3. Numerical examples and discussions

3.1 SH-waves results

First, from a practical point of view, the number of layers in the 1D PnC structure should take a finite number. Therefore, we consider the unit cell of the PnC structure is made from two layers. The two layers are lead and epoxy materials and denoted by the symbols A and B, respectively. Second, dispersion relations are plotted for an infinite number of unit cells because it depends on Bloch theory that manipulates the propagation of waves through infinite periodic structures. Therefore, the dispersive behavior of the periodic materials and structures with an arbitrary chosen unit cell configuration is considered as an example. The constants of the materials used in the calculations can be found in Refs. [46, 49] and in Table 1. Here, we consider the velocity cof the incident wave is 800 m/s and the angle of the incidence is θ0=20°. Also, the ambient temperature is proposed to be T=20°C. The dispersive properties of the inhomogeneous structures are determined by the properties ratios of the constituent materials, here we consider the two materials thicknesses ratio as 1:1.

MaterialsMass density ρ×103 (kg/ m3)Lame’ constant λ×1010 (N/m2)Shearing modulus μ×1010 (N/m2)Young’s modulus E × 1010 (N/m2)Poisson’s ratios νThermal expansion coefficient β × 106 (1/°C)Specific heat Cν×103 (J/kg. °C)
Lead11.43.30.541.5360.4329.50.128
Epoxy1.1800.4430.1590.4350.36822.51.182
Aluminum2.6996.12.56.7520.35523.90.9
Gold19.3215.02.858.1140.4214.20.13
Nylon1.110.5110.1220.3570.4501.70

Table 1.

Material constants.

Figure 2 presents the dispersion relation of SH-waves in the first Brillouin zone, and it is plotted between the nondimensional frequency q=ωa/cBand the nondimensional wave number ξ=k×a, where cBis the wave velocity in the second material (epoxy). We considered the range of the nondimensional frequency is 0q9. From this figure, we can conclude that, first, there are some frequency regions that are plotted with the white color and represent the real valued wave number and pass bands. Second, the other frequency regions corresponding to the imaginary wave number represent the phononic band gaps.

Figure 2.

The dispersion curve of SH-waves incident on the face of the 1D binary perfect PnC structure. Each unit cell consists of lead and epoxy materials (stop-bands shaded by the gray color).

3.1.1 Influence of the angle of incidence on the phononic band gap

Figure 3(a) and (b) show the effects of the incident angle θ0on the band structures of the previous perfect structure. It can be seen that the wave propagation behavior of the perfect PnC changes obviously for different incident angles. For example, the stop bands for θ0= 15° become a pass band for θ0= 45°. With increasing value of the incident angle, the same stop bands became wider in the considered frequency regions. However, in Figure 3(c), we can note a wonderful phenomenon appeared when the angle θ0reached the value 85°. The stop band increased to be the entire band structure of the PnC, and no propagation bands were observed. The explanation of such phenomenon is deduced from the parameter qj=11+χjc2/cj2sin2θ0; if ccjbecome smaller thansinθ0,the displacement field will decay. Therefore, a total reflection to the elastic waves will be occurred, where all the incident wave energy is reflected back to the structure and no waves are allowed to propagate through the PnC structure.

Figure 3.

The calculated dispersion curves of SH-waves in a 1D binary perfect PnC at the incident angles (a) θ0°= 30°, (b) θ0°= 50°, and (c) θ0°= 70°.

3.1.2 Defective mode effect on the band structure of the PnC

As shown in Figure 4, the dispersion relation of the perfect PnC structure will differ than the defected ones. We consider a defect layer from aluminum (the defect layer thickness ad= aAi.e., 1/2 of any unit cell thickness) was immersed after the second unit cell, and the material properties are mentioned in Table 1. From Figure 5, it was indicated that elastic waves can be trapped within the phononic band gap and the band structure changed significantly than the perfect ones. The width of the band gaps increased due to the increment of mismatch between the constituent materials, this back to the insertion of a second interface inside the PnC structure. Moreover, the defect layer acts as a trap inside the PnC structure, so a special wave frequency corresponding to that waveguide will propagate through the structure.

Figure 4.

A schematic diagram of a defect 1D binary PnC.

Figure 5.

The dispersion curve of SH-waves propagating in the defected PnC. A defect layer from aluminum with thickness a d = a A was immersed between the two periodic unit cells (stop-bands shaded with the gray color).

Additionally, in Figure 6(a) (ad= 2aA) and Figure 6(b) (ad= 4aA), we studied the effects of the defect layer thickness on the position and number of the localized modes inside the band structure of the PnC. It was shown that the number and width of the localized modes were increased by increasing the defect layer thickness. Therefore, with increasing the thickness of the defect layer inside the periodic PnC structure, the localized modes within the band gap are strongly confined within the defect layer.

Figure 6.

The calculated dispersion curve of SH-waves propagating in a 1D binary defected PnC. (a) a d = 2 a A (Al), (b) a d = 4 a A (Al), (c) a d = a A (Au), and (d) a d = a A (Nylon) (stop-bands shaded with the gray color).

Not only the thickness of the defect layer has an obvious effect on wave localization but also the type of the defect layer has a significant effect on the localized modes inside the band gaps. Although we introduced two different materials in Figure 6(c) and (d) with the same thickness, the number and width of the localized peaks had greatly changed. In Figure 6(c), we used Au which has mechanical constants higher than the host materials, while in Figure 5(d), we used nylon which has mechanical constants lower than the host materials. As a result, the width and number of the transmission bands are larger in Figure 6(d) than in Figure 6(c). Actually, the nylon is a very soft material like a spring and can introduce more resonant modes inside the PnC structure.

3.1.3 Temperature effects on the band structure of PnC

Now we will investigate the effects of temperature elevation on the PnC structure and phononic band gaps. Two temperature degrees (T = 35and 180°C) were considered. Figure 7 shows the response of the dispersion relations with different temperatures at SH-waves propagation through PnCs. We can note in Figure 7(a) and (b) that, with increasing the temperature from T=35to T=180°C, respectively, the pass/stop band width remains constant due to the opposite increment in the thermal stress, which has a negative value and maintains materials dimensions constant. Therefore, there is not any variation in the thickness of any material, which in turn keeps the width of the band gaps constant. Finally, we can note only band edges changed slightly with increasing the temperatures. Such small shift in the band gaps is related to the different thermal expansion coefficients of the two materials. Based on this result, the effects of temperature on the PnC structure are greatly related to the type of waves that propagate through crystal structure. Hence, we will verify this result by studying the temperature effects on PnCs at plane wave propagation as well.

Figure 7.

The dispersion curve of SH-waves propagates in a PnC from lead and epoxy materials. Different temperatures are considered. (a) T = 35 ° C and (b) T = 180 ° C (stop-bands shaded in gray).

3.2 Plane waves results

As depicted in Figure 1, we can calculate the reflection coefficient of S- and P-waves in the x-direction through the PnC structure. The PnC structure is proposed to be bonded between two semi-infinite materials (nylon material) at the two ends. The subscripts “0” and “e” denote the left and the right of the PnC structure, respectively.

The reflection coefficient of the displacement field through a PnC structure is given by the form [39],

U1U0=T12+E0T11E0EeT21EeT22E0T11EeT21T12EeT22,E26

where U1is the reflected amplitude and Tij=Tijare the elements of the total transfer matrix T=TnTn1TmT1.

Figure 8 shows the relation between the reflectance R versus ωa2PicT(ωis considered the angular frequency in this relation and cT=cSB) for P-wave (red dashed lines) and S-wave (black solid lines) [48]. From Figure 8, we can determine the range of frequencies for which the phononic band gaps can be occurred (reflectance of P- and S-waves is high R ≈ 1). Such high reflectance within the different frequency ranges represents the frequency band gaps of P- and S-waves inside the PnC structure.

Figure 8.

The reflectance R versus ω a / 2 π c T of P-wave (red lines) and S-waves (blue lines) propagated normally through a 1D PnC structure consist of four unit cells. Each unit cell consists of lead and epoxy materials.

From Figure 8, it can be seen that the phononic band gaps described by the reflection coefficient are agreed with those described by the dispersion relations. The frequency regions at R ≈ 1 of P- and S-waves represent the complete band gaps of the plane waves through the PnC structure. There is only a fine difference in the band gaps edges between Figure 8 and those obtained by the dispersion relation. In Figure 8, the relation is plotted between the reflectance and ωa2PicT, and cT=cSBis chosen for both P- and S-waves reflectance. Therefore, we will use the reflection coefficient to describe the effects of the defect layer and temperature instead of dispersion relations in order to plot both P- and S-waves in the same graph and for the practical and industrial purpose as well.

3.2.1 Influences of the defect layer/temperature on the phononic band gaps “plane wave”

In this section, we will study the effects of the defect layer and temperatures on the band structure of PnCs at the propagation of plane waves and compare with those investigated for SH-waves.

3.2.1.1 Defect layer influences on band structure

First, we will use the same defected structure used in Section 3.1.2 with the same materials and conditions, only the angle of incidence will be maintained at θ0=0°. Figure 9 confirms the last results of the effects of the defect layer on the localization modes through the PnC at the plane wave propagation. In Figure 9(b), a number of the localized waves was generated inside the phononic band gaps and was increased by increasing the defect layer thickness as well.

Figure 9.

(a) Aluminum defect layer with a d = a A , (b) Aluminum defect layer with a d = 4 a A .

3.2.1.2 Temperature influences on band structure

In this section, the two temperatures T = 50°Cand 190°Cwere considered in order to illustrate the effects of temperature on the phonic band gaps. We noticed that the phononic band gaps were affected slightly by temperatures at plane wave propagation higher than SH-waves. Therefore, the propagation of elastic waves and localized modes can be affected by temperature elevation. From Figure 10(a), we can note that the reflectance of the P-wave (Red lines) is moved toward the higher frequencies (i.e., band gap at ωa/2πcT= 3.5). Such displacement in the band gap edges is quite noticeable at T= 190°Cin Figure 10(b).

Figure 10.

The reflectance R versus ω a / 2 π c T for P-wave (red lines) and S-waves (blue lines) at temperatures effects (a) T = 35 ° C and (b) T = 180 ° C .

These temperature effects on the band gaps can be explained by two reasons. First, the P-wave velocity is increased according to Eq. (24) because the temperature has a direct effect on the elastic constants. Consequently, temperature makes a notable change in the edges of the phononic band gaps. Second, here in the plane wave case, the thermal stress is absent, so the lattice constants will increase according to Eq. (25) and can make change in the width of the phononic band gaps as well.

4. Conclusions

In this chapter, we used the transfer matrix method and Bloch theory to compute and study the propagation of mechanical waves through different 1D PnC structures with thermal effects. We performed numerical simulations to calculate the reflection coefficient and dispersion curves for the 1D PnC structure. We have discussed and studied the effects of many parameters on the phononic band gaps. The effects of temperature and defect mode on the properties of the PnC and on the localized modes were studied for both SH-waves and in-plane waves. The output results can be stated as follows:

  1. The defect layer has a different and unique effect on the number of localized modes inside the phononic band gaps for both plane and SH-waves. The effects of the incident angle, thickness ratio on the propagation, and localization of elastic waves were studied. The number of localized modes was increased with increasing the defect layer thickness. Also, the material type of the defect layer has an obvious effect on the width and number of resonant defect modes, where the localized modes were increased by using a material with low elastic constant. Inserting a defect layer inside the perfect structure facilitates the process of localized modes generation inside the PnC structure. Also, it can be utilized in many potential applications such as acoustic filter and waveguides and can be applied to produce new types of laser called phonon laser with low thermal effects especially in the hypersonic PnCs regime.

  2. Results revealed that temperature increment has a significant effect on the phononic band gaps and on the number of the localized modes inside the PnC structure, especially for plane wave, where the P-waves velocity and lattice constants are directly related to temperature increment. The contribution of phonons in the thermal conductivity of a solid can be limited by the existence of phononic band gaps, which could be very useful for the thermoelectric devices that convert thermal energy into electricity. Hence after, the performance of many devices such as Peltier thermoelectric coolers, thermocouples, sensors, and thermoelectric energy generators can be enhanced. Moreover, the thermal conductivity of many structures could experience great changes based on the temperature gradients of PnCs. When temperature is increase in PnCs, the phononic band gap experiences a significant change.

Conflict of interest

The authors declare that there is no conflict of interest.

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Arafa H. Aly and Ahmed Mehaney (April 15th 2019). Phononic Crystals and Thermal Effects, Photonic Crystals - A Glimpse of the Current Research Trends, Pankaj Kumar Choudhury, IntechOpen, DOI: 10.5772/intechopen.82068. Available from:

chapter statistics

172total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Metal-Matrix Embedded Phononic Crystals

By Suobin Li, Yihua Dou and Linkai Niu

Related Book

First chapter

Introductory Chapter: A Revisit to Optical Amplifiers

By Pankaj Kumar Choudhury

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us