Athanasios G. Polimeridis

Samuel P. Groth, Athanasios G. Polimeridis, Alexandra Tambova et al.

Recently, the volume integral equation (VIE) approach has been proposed as an efficient simulation tool for silicon photonics applications [J. Lightw. Technol. 36, 3765 (2018)]. However, for the high-frequency and strong contrast problems arising in photonics, the convergence of iterative solvers for the solution of the linear system can be extremely slow. The uniform discretization of the volume integral operator leads to a three-level Toeplitz matrix, which is well suited to preconditioning via its circulant approximation. In this paper, we describe an effective circulant preconditioning strategy based on the multi-level circulant preconditioner of Chan and Olkin [Numer. Algorithms 6, 89 (1994)]. We show that this approach proves ideal in the canonical photonics problem of propagation within a uniform waveguide, in which the flow is unidirectional. For more complex photonics structures, such as Bragg gratings, directional couplers, and disk resonators, we generalize our preconditioning strategy via geometrical partitioning (leading to a block-diagonal circulant preconditioner) and homogenization (for inhomogeneous structures). Finally, we introduce a novel memory reduction technique enabling the preconditioner's memory footprint to remain manageable, even for extremely long structures. The range of numerical results we present demonstrates that the preconditioned VIE is fast and has great utility for the numerical exploration of prototype photonics devices.

Ilias I. Giannakopoulos, Mikhail S. Litsarev, Athanasios G. Polimeridis

We present a method of memory footprint reduction for FFT-based, electromagnetic (EM) volume integral equation (VIE) formulations. The arising Green's function tensors have low multilinear rank, which allows Tucker decomposition to be employed for their compression, thereby greatly reducing the required memory storage for numerical simulations. Consequently, the compressed components are able to fit inside a graphical processing unit (GPU) on which highly parallelized computations can vastly accelerate the iterative solution of the arising linear system. In addition, the element-wise products throughout the iterative solver's process require additional flops, thus, we provide a variety of novel and efficient methods that maintain the linear complexity of the classic element-wise product with an additional multiplicative small constant. We demonstrate the utility of our approach via its application to VIE simulations for the Magnetic Resonance Imaging (MRI) of a human head. For these simulations we report an order of magnitude acceleration over standard techniques.

Ioannis P. Georgakis, Ilias I. Giannakopoulos, Mikhail S. Litsarev et al.

A stable volume integral equation (VIE) solver based on polarization/magnetization currents is presented, for the accurate and efficient computation of the electromagnetic scattering from highly inhomogeneous and high contrast objects.We employ the Galerkin Method of Moments to discretize the formulation with discontinuous piecewise linear basis functions on uniform voxelized grids, allowing for the acceleration of the associated matrix-vector products in an iterative solver, with the help of FFT. Numerical results illustrate the superior accuracy and more stable convergence properties of the proposed framework, when compared against standard low order (piecewise constant) discretization schemes and a more conventional VIE formulation based on electric flux densities. Finally, the developed solver is applied to analyze complex geometries, including realistic human body models, typically used in modeling the interactions between electromagnetic waves and biological tissue.