Vladimir Druskin, Rob Remis
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 design
efficient frequency independent perfectly matched layers with controlled
accuracy and introduce a new Krylov-based solution method via
stability-corrected operator exponents. This approach 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, which enables us to efficiently
compute the time-domain solution. 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.
View original:
http://arxiv.org/abs/1202.0424
No comments:
Post a Comment