TY - JOUR
T1 - A simple tool for the optimization of 1D phononic and photonic bandgap filters
AU - Salasiya, Prasanna
AU - Guzina, Bojan B.
N1 - Publisher Copyright:
© 2025
PY - 2025/3/1
Y1 - 2025/3/1
N2 - 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 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” (propagating or evanescent) Bloch waves. For a given periodic medium and frequency of oscillations, the latter are 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 matrix 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 for an optimal filter (e.g. rainbow trap) design as it enables rapid exploration of the design space with respect to segment (i) permutations (with or without repetition), (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 (GA) approach that entails O(106) 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.
AB - 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 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” (propagating or evanescent) Bloch waves. For a given periodic medium and frequency of oscillations, the latter are 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 matrix 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 for an optimal filter (e.g. rainbow trap) design as it enables rapid exploration of the design space with respect to segment (i) permutations (with or without repetition), (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 (GA) approach that entails O(106) 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.
KW - Bloch waves
KW - Filter design
KW - Genetic algorithm
KW - Phononics
KW - Photonics
KW - Quadratic eigenvalue problem
UR - http://www.scopus.com/inward/record.url?scp=85214530969&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85214530969&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2025.113723
DO - 10.1016/j.jcp.2025.113723
M3 - Article
AN - SCOPUS:85214530969
SN - 0021-9991
VL - 524
JO - Journal of Computational Physics
JF - Journal of Computational Physics
M1 - 113723
ER -