Theoretical study of two-dimensional phononic crystals with viscoelasticity based on fractional derivative models

Wave propagation in two-dimensional phononic crystals (PCs) with viscoelasticity is investigated using a finite-difference-time-domain (FDTD) method. The viscoelasticity is evaluated using the Kelvin–Voigt model with fractional derivatives (FDs) so that both the dispersion and dissipation are considered. Numerical approximation of FDs is integrated into the FDTD scheme to simulate wave propagation in such PCs. All the constituent materials are treated as isotropic and homogeneous. The gaps are substantially displaced and widened and the attenuation is noticeably enhanced due to the dispersion and dissipation of host material and the complicated multiple scattering between scatterers. These results indicate that the viscoelasticity of the damping host has significant influence on wave propagation in PCs and should be considered.