Bandgap calculation of two-dimensional mixed solid-fluid phononic crystals by Dirichlet-to-Neumann maps

A numerical method based on the Dirichlet-to-Neumann (DtN) map is presented to compute the bandgaps of two-dimensional phononic crystals, which are composed of square or triangular lattices of circular solid cylinders in a fluid matrix. The DtN map is constructed using the cylindrical wave expansion in a unit cell. A linear eigenvalue problem, which depends on the Bloch wave vector and involves relatively small matrices, is formulated. Numerical calculations are performed for typical systems with various acoustic impedance ratios of the solid inclusions and the fluid matrix. The results indicate that the DtN-map based method can provide accurate results for various systems efficiently. In particular it takes into account the fluid-solid interface conditions and the transverse wave mode in the solid component, which has been proven to be significant when the acoustic impedance of the solid inclusions is close to or smaller than that of the fluid matrix. For systems with an acoustic impedance of the inclusion much less than that of the matrix, physical flat bands appear in the band structures, which will be missed if the transverse wave mode in the solid inclusions is neglected.