Efficient computation of electromagnetic wave fields on unbounded domains using stability-corrected wave functions and Krylov subspace projection methods

We present a Krylov subspace projection framework to model electromagnetic wave propagation in unbounded domains. The extension to infinity is modeled via an optimal complex scaling method and we show that stable time-domain reduced-order models can be efficiently computed via a stability-correction procedure in conjunction with a Lanczos-type reduction algorithm. We show that dominant open resonant modes can be determined via the Lanczos algorithm and illustrate the performance of our technique through a number of numerical examples for two- and three-dimensional configurations.