For the love of physics walter lewin may 16, 2011 duration. Thomas, simultaneous complete elastic and electromagnetic band gaps in periodic structures, applied physics b. We will show how to use pwe band solver to analyze cases of both normal and oblique incidence, while pointing out other features of the software, namely output normalization and inversion symmetry. Phononic crystal structures for wireless communications and. This paper introduces the concept of the photonic crystal, the research methods and the application of photonic crystals. Atom to transistor epilogue appendixadvanced formalism selected bibliography matlab codes for text figures. As in forbidden electronic energy bands in semiconductors there are certain frequencies for which electromagnetic waves can not travel through those crystals. May 24, 2017 matlab matrix laboratory is a multiparadigm numerical computing environment and fourthgeneration programming language which is frequently being used by engineering and science students. Ultradirectional source of longitudinal acoustic waves based on a twodimensional solidsolid phononic crystal b.
In this case, the blochfloquet theorem for periodic eigenproblems states that the solutions to eq. Phononic and sonic crystals have generated rising scientific interest for very diverse technological applications. Photonic crystal simulation toolkit including band solver. Similar to the frequency bandgap to electromagnetic wave propagation in photonic crystals, it is also possible for phononic crystals to exhibit. Numerical simulation of phonon dispersion relations for. Pdf surface acoustic waves in two dimensional phononic. Yetisen4,5, seokhyun yun4,5 and haider butt1 abstract photonic crystals and band gap materials act as manipulators of light and have a plethora of applications. Bravais lattice for two dimensions square lattice photonic crystal. Comsol multiphysics livelink for matlab is employed with periodic boundary conditions on unit cells. By using the plane strain and mindlin plate physics in comsol, we carried out simulations to determine both the inplane and outofplane phonon dispersion relations. Deymier3 1laboratoire ondes et milieux complexes, umr cnrs 6294, universite du havre, 75 rue bellot, 76058 le havre, france.
The dispersion of the system is spatially dependent and allows the rainbow trapping inside the structure. Controlling the thermal conductivity of semiconductors is of practical interest in optimizing the performance of thermoelectric and phononic devices. Photonic and phononic crystal materials and devices ix. An acoustic metamaterial is a material designed to control, direct, and manipulate sound waves as these might occur in gases, liquids, and solids. Control analysis of the tunable phononic crystal with. Introduction photonic crystals are composed of periodic, dielectric structures that a ect the propagation of electromagnetic waves. Fdtd simulations of acoustic waves in twodimensional. Simultaneous photonic and phononic band gap, defect states. Granular crystals are one type of nonlinear periodic phononic structure and are the focus of this chapter.
Surface and bulk acoustic waves in twodimensional phononic crystal consisting of materials with general anisotr opy t sungt song w u, zigui huang, and s. Phononic crystals are composite materials which are periodically composed in 1. Computing the bandgap of a 2d photonic crystal by comsol matlab scripting this project replicates figure 6. Introduction to phononic crystals and acoustic metamaterials. A twodimensional phononic crystal simulator is developed. Parametric simulations of slanted 1d photonic crystal sensors. Phononic and magnonic dispersions of surface waves on a. Based on your location, we recommend that you select.
Complex band structures of two dimensional phononic crystals. The symmetry points of the two dimensions square lattice2. This approach offers an effective means of contactless tunability of the properties of phononic crystals. Threedimensional adaptive soft phononic crystals sahab babaee,1 pai wang,1 and katia bertoldi1,2 1school of engineering and applied sciences, harvard university, cambridge, massachusetts 028, usa 2kavli institute, harvard university, cambridge, massachusetts 028, usa received 29 april 2015. This thesis presents the design of a twodimensional phononic band gap crystal simulator, and phononic crystal analysis. This prog calculates and plots the wavevector diagram i. Analysis of rayleigh surface acoustic waves propagation on. Achieving minimal heat conductivity by ballistic confinement. Application of the plane wave expansion method to a two. In the present letter, we demonstrate that the band structure of a twodimensional 2d phononic crystal constituted of a magnetoelastic medium can be controlled by application of an external magnetic. Band structures of the 2d phononic crystal with cylindrical holes in a square lattice.
From a historical point of view, we have tried to refer to some of the seminal contributions that have made the field. Mar 26, 2015 i am trying to calculate the phononic band structure of a thin crosshole perforated silicon phononic crystal using the partial differential equation toolbox. Axel scherer for the purchase of this volume in printed format, please visit. The atoms in crystals are arranged as periodical threedimensional arrays. Design of optomechanical cavities and waveguides on a. Three dimensional phononic crystal sonic vacuum without prestress was assembled from 7 vertical cavities arranged in hexagonal pattern in silicone matrix filled with stainless steel spheres. Piezoelectric aluminum nitride aln film is deployed as the interdigital transducers idt to transmit and detect acoustic waves, thus making the whole microfabrication process cmos. In this work, the phonon dispersion relation of 2d phononic crystals was investigated using the finite element method. Numerical simulation of phonon dispersion relations for 2d phononic crystals gaohua zhu, eric dede. Inherent negative refraction on acoustic branch of two. Inset shows a square lattice and its brillouin zone.
A schematic diagram of a one dimensional phononic crystal with 3 periods. Surface and bulk acoustic waves in twodimensional phononic. Phononic crystals and their application to microwave acoustic. A phononic crystal is essentially made up of homogeneous spherical or cylindrical inclusions arranged periodically in a homogeneous host medium. Twodimensional 2d silicon phononic crystal pnc slab of a square array of cylindrical air holes in a 10.
Phononic bandgap guidance of acoustic modes in photonic. The youngs modulus, mass density and poissons ratio of the epoxy are 4. One of the main properties of the phononic crystals is the possibility of having a phononic bandgap. Photonic bands for a 2d photonic crystal file exchange. Under certain conditions, acoustic band gaps can form. The section 4 deals with the proposed algorithm and the numerical results.
Phononic band gap micronanomechanical structures for. Designing phononic crystals with convex optimization. In order to study the propagation of acoustic waves in an apc, we must understand the nature of the eigen modes of the wave in the apc. Seven chains were supported by the single plate placed on the top of piezogauge. The scattering matrix method smm 12 allows the evaluation of. Simulation of 2d photonic crystal with comsol multiphysics. For example, the pncs can prohibit the propagation of acoustic elastic. Phononic crystal, planewave expansion method, phononic band gaps. The band structure reproduce using my matlab code in appendix b, for the same parameters as in ref. These are spectral bands where propagation of waves is forbidden. Phononic crystal, as an elastic and acoustic analogue of the photonic crystal, is a periodic microstructure created from the arrangement of composite elastic materials.
Abstractin this paper, two surface acoustic wave saw. In other words, the wave in a pc is a floquetbloch wave, which is the superposition of many plane waves from all individuals. Design parameters for phononic crystal band gap engineering are outlined. A simple tutorial illustrates how the basic features are used to solve for the modes of. Frequency is normalized to facl, where cl is the longitudinal sound speed of. It has oscillatory front and waves reflected from the cover plate back to the gauge are evident. Watson, sr, laboratory of applied physics, california institute of technology, pasadena, california 91125, usa. Starting from maxwells equations,the functional form of the te mode and tm mode of the 2d crystal is derived. Phononic crystal phononic crystals are synthetic materials that are formed by periodic variation of the acoustic properties of the material i. I am trying to calculate the phononic band structure of a thin crosshole perforated silicon phononic crystal using the partial differential equation toolbox. A phononic crystal pc is a metamaterial engineered material with periodic variations in mass or stiffness.
This book provides an indepth analysis as well as an overview of phononic crystals. Reconfigurable phononiccrystal circuits formed by coupled. Nonlinearity is used to activate the selfdemodulation effect, which is enhanced due to the particular dispersion characteristics of. Bringuier a thesis submitted to the faculty of the department of materials science and engineering in partial fulfillment of the requirements for the degree of masters of science in the graduate college the university of arizona 2011. The observed phononic gaps are considerably larger than those of laterally patterned multicomponent crystals previously reported, mainly a consequence of the high elastic and density. Basic crystal systems were designed following the guide shown in figure 1.
The first acoustic metamaterial, also called as locally resonant sonic materials was demonstrated with negative effective dynamic density. Below is the matlab program that calculates and plots the re. Pdf on sep 30, 2019, yunfeng gu and others published simulation of elastic. Dispersion curves for onedimensional lattice 29 appendix c. A twodimensional squarelattice phononic crystal is chosen for the following examples. Choose a web site to get translated content where available and see local events and offers. The transmitted waves are then refocused at the right side of slab. The insertion of inclusions of nanometer size in a semiconductor is an effective means of achieving such control. At frequencies and wave vectors where the refraction on the acousticbranch passbands is negative, the e ective massdensity and the e ective sti ness tensors of the. When engineering these crystals, it is possible to isolate vibration within a certain frequency range.
The phononic crystal concept originally comes from crystallography. The hereditary line into acoustic metamaterials follows from theory and research in negative index material. Considering a phononic structure consisting of epoxy circular cylinders embedded in a background material of er materials to form the twodimensional pcs with lattice spacing a1 cm. This dependence is a common feature of phononic and photonic crystals 21. Determination of acoustic scattering from a 2d finite. However, i dont know how to setup periodic blochfloquet boundary conditions so that i can loop over a set of kx and ky wavevectors. In this tutorial we design simple bragg grating that consists of alternating layers with permittivity contrast 1 as discussed and analyzed in 1. Pdf simulation of elastic wave transmission in phononic crystal. Phononic crystals phononic crystal are materials that have periodic variations in their. Acoustic waves are emitted from the line source and propagate into the phononic crystal slab, which has negative refractive index. Explicit formulations of the plane harmonic bulk wave and the surface wave dispersion relations in such a general phononic structure are derived based on the plane wave expansion method. The phononic model is based on the acoustic wave equations, and the photonic one in the electromagnetic equations.
A photonic crystal corresponds to a periodic dielectric function. This chapter describes the dynamic behavior of nonlinear periodic phononic structures, along with how such structures can be utilized to affect the propagation of mechanical waves. Twodimensional phonon dispersion curve 31 appendix d. As an example of the possible applications of these pnbg structures, i have.
Sep 14, 2017 the latter is used to derive the endtoend transfer function of a finite phononic crystal as a function of any given parameters. For this value, a complete band gap is known to appear for 1535 ms tutorial. Design of optomechanical cavities and waveguides on a simultaneous bandgap phononicphotonic crystal slab amir h. Modeling phononic band gap materials and structures. Numerical simulation of phonon dispersion relations for 2d. A phononic crystal is an artificially manufactured structure, or material, with periodic constitutive or geometric properties that are designed to influence the characteristics of mechanical wave propagation. A novel micromechanical resonator using twodimensional. The objective of this chapter is to introduce the broad subject of phononic crystals and acoustic metamaterials.
Pdf an analysis is given to the band structure of the two dimensional solid phononic crystal considered as a semi infinite medium. Twodimensional phononic crystal simulation and analysis. Beside the physical properties of the pc crystal component materials, the phononic gap width is found to depend strongly on the lattice symmetry, as well as on the scatterer shape 19, 20. Based on the existing information on the er material preyield rheology, only the electric. Jan 25, 2009 this prog calculates and plots the wavevector diagram i. This book discusses numerous techniques for the analysis of phononic crystals and covers, among other material, sonic and ultrasonic structures, hypersonic planar structures and their characterization, and novel applications of phononic crystals. Dynamic control of phonon propagation in phononic crystal. Furthermore, with acoustic metamaterials, controlling sonic waves can now be extended to the negative refraction domain. Ultradirectional source of longitudinal acoustic waves. On chip complex signal processing devices using coupled.
The normalized dispersion curves are shown for three different poissons ratios. Numerical simulation of phonon dispersion relations for 2d phononic crystals gaohua zhu, eric dede toyota research institute of north america 10032012 excerpt from the proceedings of the 2012 comsol conference in boston. This program calculates and plots the photonic bands for a 2d photonic crystal consisting of cylinders with circular crosssection and infinite height, arranged in a triangular lattice. Intheexampleshownbelow,thephotoniccrystalismadeofn 10 alternating homogeneous layers with. Study of phononic and photonic crystal problems by.
We created an fdtd simulator in the c language that uses matlab as a frontend for providing data to and. Asymmetric propagation using enhanced selfdemodulation in. Simulator operation is validated through comparison with published data. Phononic band gap width control through structural and. Asymmetric propagation of acoustic waves is theoretically reported in a chirped phononic crystal made of the combination of two different nonlinear solids. Digital signal processing and wave mechanics are utilized to analyze fractal and circular inclusion based phononic crystals. Nano express open access parametric simulations of slanted 1d photonic crystal sensors aaron breuerweil1, naif nasser almasoud2, badaruddin abbasi3, ali k. These crystals are made of periodic distributions of scatterers embedded in a matrix. Phononic crystal waveguide transducers for nonlinear elastic. To the left side of the crystal is an acoustic pressure source producing plane waves. Analysis of rayleigh surface acoustic waves propagation on piezoelectric phononic crystal with 3d finite element model honglang li, yahui tian, yabing ke, shitang he institute of acoustics chinese academy of sciences beijing, china email.
Black represent low dielectric and white represent high dielectric. The existence of acoustic bandgaps for outofplane propagation in a twodimensional solidsolid phononic crystal has, however, been recently demonstrated using a planewave expansion approach 10. These are spectral bands where propagation of waves. The analysis reveals intriguing features that explain the evolution of bragg band gaps in the frequency response. Fourier coefficients for the expansion of dielectric. Plane wave expansion method for photonic band gap calculation. Periodic blochfloquet boundary conditions for calculating. Pdf the studies of propagation elastic waves into composite materials have become considerably more frequent in the last years, especially. Surface waves can be easily excited at the surface of a finite size phononic crystal by line source or gaussian beam placed in or launched from the background medium, and they propagate along the. Crystalwave includes a layout editor dedicated to the design of photonic crystal structures and it is packed with a number of highperformance simulation engines, making it the most comprehensive tool for modelling of photonic crystals.
Matlab, and the results are normalized to the transmission of the slab. Phononic crystals pcs with periodic structures are the mechanical analogous of wellknown. Oct 21, 2008 this program calculates and plots the photonic bands for a 2d photonic crystal consisting of cylinders with circular crosssection and infinite height, arranged in a triangular lattice. In this chapter, we employed phononic crystal strip in mems resonators is explained to. Analyzing 1d photonic crystals bragg gratings introduction. Epol and hpol efield and hfield are parallel to the cylinders, respectively. Nonlinear periodic phononic structures and granular crystals. Photonic and phononic crystal materials and devices ix editors. As bloch 48 pointed out, a wave propagating in a periodic structure is different from a plane wave only by a periodic modulation. First band structure of phononic crystals consisting of metal cylinders placed in di. Phononic crystals and acoustic metamaterials sciencedirect. So, the complexvalued magnetic field h can be expressed as.
1388 963 1599 921 578 185 1003 583 1009 598 293 1355 517 561 224 241 1049 1168 332 489 1246 1251 259 111 547 1372 1409 337 1498 1075 353 255 203 1132 1225 602 1422 518