Combining plane wave expansion and variational techniques for fast phononic computations

TL;DR

This paper introduces a combined approach using Plane Wave Expansion and variational techniques to efficiently compute phononic band structures, achieving high accuracy and stability with reduced computational cost.

Contribution

It develops a novel method that integrates Fourier-based material property expansion with a mixed variational scheme, eliminating explicit integration and enhancing convergence.

Findings

Method provides stable, accurate eigenvalues across expansion ranges.

Outperforms zero-order numerical integration in accuracy.

Achieves similar accuracy to higher-order schemes with less computational effort.

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 allows us to provide succinct and closed form expressions for all the matrices involved in the mixed variational method. The scheme proposed here preserves both the simplicity of expression which is inherent in the PWE method and the superior convergence properties of the mixed variational scheme. We present numerical results and comment upon the convergence and stability of the current method. We show that 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.