Eric Michielssen

Han Guo, Yang Liu, Jun Hu et al.

A butterfly-based direct combined-field integral equation (CFIE) solver for analyzing scattering from electrically large, perfect electrically conducting objects is presented. The proposed solver leverages the butterfly scheme to compress blocks of the hierarchical LU-factorized discretized CFIE operator and uses randomized butterfly reconstruction schemes to expedite the factorization. The memory requirements and computational cost of the direct butterfly-CFIE solver scale as $O(N\mathrm{log}^2N)$ and $O(N^{1.5}\mathrm{log}N)$, respectively. These scaling estimates permit significant memory and CPU savings when compared to those realized by low-rank (LR) decomposition-based solvers. The efficacy and accuracy of the proposed solver are demonstrated through its application to the analysis of scattering from canonical and realistic objects involving up to 14 million unknowns.

Alexander Litvinenko, Abdulkadir C. Yucel, Hakan Bagci et al.

Computational tools for characterizing electromagnetic scattering from objects with uncertain shapes are needed in various applications ranging from remote sensing at microwave frequencies to Raman spectroscopy at optical frequencies. Often, such computational tools use the Monte Carlo (MC) method to sample a parametric space describing geometric uncertainties. For each sample, which corresponds to a realization of the geometry, a deterministic electromagnetic solver computes the scattered fields. However, for an accurate statistical characterization the number of MC samples has to be large. In this work, to address this challenge, the continuation multilevel Monte Carlo (CMLMC) method is used together with a surface integral equation solver. The CMLMC method optimally balances statistical errors due to sampling of the parametric space, and numerical errors due to the discretization of the geometry using a hierarchy of discretizations, from coarse to fine. The number of realizations of finer discretizations can be kept low, with most samples computed on coarser discretizations to minimize computational cost. Consequently, the total execution time is significantly reduced, in comparison to the standard MC scheme.

Yang Liu, Han Guo, Eric Michielssen

A butterfly-based fast direct integral equation solver for analyzing high-frequency scattering from two-dimensional objects is presented. The solver leverages a randomized butterfly scheme to compress blocks corresponding to near- and far-field interactions in the discretized forward and inverse electric field integral operators. The observed memory requirements and computational cost of the proposed solver scale as O(Nlog^2N) and O(N^1.5 logN), respectively. The solver is applied to the analysis of scattering from electrically large objects spanning over ten thousand of wavelengths and modeled in terms of five million unknowns.