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

We develop an effective computational tool for simulating the scattering of one-dimensional (1D) waves by a composite layer architected in an otherwise homogeneous medium. The is designed as the union of segments cut from various mother periodic media, which us to describe the wavefield in each segment in terms of the "left" and "right" (propagating evanescent) Bloch waves. For a given periodic medium and frequency of oscillations, the latter computed by solving the quadratic eigenvalue problem (QEP) which seeks the (real- or complex valued) 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 approach, that seeks the amplitudes of the featured Bloch waves (two per segment), amplitude of the reflected wave, and that of the transmitted wave. Such an approach inherently caters an optimal filter (e.g. rainbow trap) design as it enables rapid exploration of the design with respect to segment (i) permutations (with or without repetition), (ii) cut lengths, 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 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 (GA) approach that entails O(106) trial configurations. Relative classical rainbow trap design where the unit cells of the mother periodic media are arranged 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, the same filter thickness.