A simple tool for the optimization of 1D phononic and photonic bandgap filters

TL;DR

This paper introduces a computational tool for efficiently designing 1D phononic and photonic bandgap filters by simulating wave scattering in composite layered media, enabling rapid optimization of filter configurations.

Contribution

The authors develop a low-dimensional algebraic transfer matrix method that allows fast exploration of filter design parameters, significantly speeding up the optimization process.

Findings

Achieves 40x computational speedup in filter optimization.

GA-optimized configuration reduces filter transmissibility by 40%.

Method accurately matches numerical simulations.

Abstract

We develop an effective computational tool for simulating the scattering of 1D waves by a composite layer architected in an otherwise homogeneous medium. The layer is designed as the union of segments cut from various mother periodic media, which allows us to describe the wavefield in each segment in terms of the ``left'' and ``right'' Bloch waves. For a given periodic medium and frequency of oscillations, the latter are computed by solving the quadratic eigenvalue problem which seeks the wavenumber -- and affiliated eigenstate -- of a Bloch wave. In this way the scattering problem is reduced to a low-dimensional algebraic problem, solved via the transfer matrix approach, that seeks the amplitudes of the featured Bloch waves, amplitude of the reflected wave, and that of the transmitted wave. Such an approach inherently caters for an optimal filter design as it enables rapid exploration of the design space with respect to segment (i) permutations, (ii) cut lengths, and (iii) cut offsets relative to the mother periodic media. Specifically, under (i)--(iii) the Bloch eigenstates remain invariant, so that only the transfer matrices need to be recomputed. The reduced order model is found to be in excellent agreement with numerical simulations. Example simulations demonstrate 40x computational speedup when optimizing a 1D filter for minimum transmission via a genetic algorithm approach that entails $O(10^6)$ trial configurations. Relative to the classical rainbow trap design where the unit cells of the mother periodic media are arranged in a ``linear'' fashion according to their dispersive characteristics, the GA-optimized (rearranged) configuration yields a $40\%$ reduction in filter transmissibility over the target frequency range, for the same filter thickness.