Combining Plane Wave Expansion and Variational Techniques for Fast Phononic Computations

Abstract

In this paper, the salient features of the plane wave expansion (PWE) method and the mixed variational technique are combined for the fast eigenvalue computations of arbitrarily complex phononic unit cells. This is done by expanding the material properties in a Fourier expansion, as is the case with PWE. The required matrix elements in the variational scheme are identified as the discrete Fourier transform coefficients of material properties, thus obviating the need for any explicit integration. The process provides succinct, closed-form expressions for all the matrices involved in the mixed variational method. The scheme proposed preserves both the simplicity of expression that is inherent in the PWE method and the superior convergence properties of the mixed variational scheme. This paper presents numerical results and comments upon the convergence and stability of the current method. The current representation renders the results of the method stable over the entire range of the expansion terms as allowed by the spatial discretization. When compared with a zero-order numerical integration scheme, the present method results in greater computational accuracy of all eigenvalues. A higher-order numerical integration scheme comes close to the accuracy of the present method but only with significantly more computational expense.