A Petrov-Galerkin finite element interface method for interface problems with Bloch-periodic boundary conditions and its application in phononic crystals

Abstract

In this paper, we propose a Petrov-Galerkin finite element interface method (PGFEIM) to solve the elliptic and elastic interface problems with Bloch-periodic boundary conditions. The main idea of this method is to choose the standard finite element basis independent of the interface to be the test function basis, and choose a piecewise linear function satisfying the jump conditions across the interface to be the solution basis. The grid we use is a non-body-fitted grid. The PGFEIM is capable of dealing with sharp-edged interface problems with matrix coefficients and nonhomogeneous jump conditions. Further, we extend this method to compute the band structure of anti-plane transverse waves and in-plane elastic waves in phononic crystals. Different acoustic impedance ratios, arbitrary complex scatterer geometries and various material properties are considered and discussed. Numerical experiments demonstrate that the PGFEIM is nearly second order accurate in the L∞ norm for interface problems, and it is accurate and efficient for computing the band structure of phononic crystals.