A Krylov Stability-Corrected Coordinate-Stretching Method to Simulate Wave Propagation in Unbounded Domains

TL;DR

This paper introduces a Krylov subspace method with stability correction and spectral accuracy for simulating wave propagation in unbounded domains, outperforming traditional FDTD methods on large time scales.

Contribution

It develops a novel Krylov-based approach with stability correction and spectral PMLs for accurate, stable reduced-order modeling of wave problems in unbounded domains.

Findings

Unconditionally stable reduced-order models for wave scattering.

Spectrally accurate perfectly matched layers using Zolotarev approximation.

Outperforms finite-difference time domain methods on large time intervals.

Abstract

The Krylov subspace projection approach is a well-established tool for the reduced order modeling of dynamical systems in the time domain. In this paper, we address the main issues obstructing the application of this powerful approach to the time-domain solution of exterior wave problems. We use frequency independent perfectly matched layers to simulate the extension to infinity. Pure imaginary stretching functions based on Zolotarev's optimal rational approximation of the square root are implemented leading to perfectly matched layers with a controlled accuracy over a complete spectral interval of interest. A new Krylov-based solution method via stability-corrected operator exponents is presented which allows us to construct reduced-order models (ROMs) that respect the delicate spectral properties of the original scattering problem. The ROMs are unconditionally stable and are based on a renormalized bi-Lanczos algorithm. We give a theoretical foundation of our method and illustrate its performance through a number of numerical examples in which we simulate 2D electromagnetic wave propagation in unbounded domains, including a photonic waveguide example. The new algorithm outperforms the conventional finite-difference time domain method for problems on large time intervals.