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

The Krylov subspace projection approach is a well-established tool for the reducedorder 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 stabilitycorrected 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 two-dimensional electromagnetic wave propagation in unbounded domains, including a photonic waveguide example. The new algorithm outperforms the conventional finitedifference time-domain method for problems on large time intervals.

[1]  N. Moiseyev,et al.  Quantum theory of resonances: calculating energies, widths and cross-sections by complex scaling , 1998 .

[2]  David J. Goodman,et al.  Personal Communications , 1994, Mobile Communications.

[3]  Gunilla Kreiss,et al.  Perfectly Matched Layers for Hyperbolic Systems: General Formulation, Well-posedness, and Stability , 2006, SIAM J. Appl. Math..

[4]  Patrick Joly,et al.  Stability of perfectly matched layers, group velocities and anisotropic waves , 2003 .

[5]  Athanasios C. Antoulas,et al.  Approximation of Large-Scale Dynamical Systems , 2005, Advances in Design and Control.

[6]  James Demmel,et al.  Applied Numerical Linear Algebra , 1997 .

[7]  Vladimir Druskin,et al.  On convergence of Krylov subspace approximations of time-invariant self-adjoint dynamical systems , 2012 .

[8]  J. Combes,et al.  Spectral properties of many-body Schrödinger operators with dilatation-analytic interactions , 1971 .

[9]  R. Freund,et al.  Software for simplified Lanczos and QMR algorithms , 1995 .

[10]  Vladimir Druskin,et al.  Application of the Difference Gaussian Rules to Solution of Hyperbolic Problems , 2000 .

[11]  Michael E. Taylor,et al.  Partial Differential Equations II: Qualitative Studies of Linear Equations , 1996 .

[12]  Vladimir Druskin,et al.  Gaussian Spectral Rules for the Three-Point Second Differences: I. A Two-Point Positive Definite Problem in a Semi-Infinite Domain , 1999, SIAM J. Numer. Anal..

[13]  Serkan Gugercin,et al.  Interpolatory projection methods for structure-preserving model reduction , 2009, Syst. Control. Lett..

[14]  V. Edwards.  Scattering Theory , 1973, Nature.

[15]  W. Koch,et al.  On resonances in open systems , 2004, Journal of Fluid Mechanics.

[16]  Frank Olyslager.  Discretization of Continuous Spectra Based on Perfectly Matched Layers , 2004, SIAM J. Appl. Math..

[17]  Israel Michael Sigal,et al.  Introduction to Spectral Theory , 1996 .

[18]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[19]  Joseph E. Pasciak,et al.  The computation of resonances in open systems using a perfectly matched layer , 2009, Math. Comput..

[20]  Murthy N. Guddati,et al.  On Optimal Finite-Difference Approximation of PML , 2003, SIAM J. Numer. Anal..

[21]  D. Wulbert.  The Rational Approximation of Real Functions , 1978 .

[22]  J. Joannopoulos,et al.  High Transmission through Sharp Bends in Photonic Crystal Waveguides. , 1996, Physical review letters.

[23]  Barry Simon,et al.  Resonances in n-Body Quantum Systems With Dilatation Analytic Potentials and the Foundations of Time-Dependent Perturbation Theory , 1973 .

[24]  L. Knizhnerman,et al.  Two polynomial methods of calculating functions of symmetric matrices , 1991 .

[25]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[26]  Vladimir Druskin,et al.  Optimal finite difference grids and rational approximations of the square root I. Elliptic problems , 2000 .

[27]  J. Combes,et al.  A class of analytic perturbations for one-body Schrödinger Hamiltonians , 1971 .

[28]  Anne Greenbaum,et al.  Using Nonorthogonal Lanczos Vectors in the Computation of Matrix Functions , 1998, SIAM J. Sci. Comput..

[29]  Murthy N. Guddati,et al.  Absorbing boundary conditions for scalar waves in anisotropic media. Part 2: Time-dependent modeling , 2010, J. Comput. Phys..

[30]  L. Knizhnerman,et al.  Extended Krylov Subspaces: Approximation of the Matrix Square Root and Related Functions , 1998, SIAM J. Matrix Anal. Appl..

[31]  Patrick Joly,et al.  Mathematical Modelling and Numerical Analysis on the Analysis of B ´ Erenger's Perfectly Matched Layers for Maxwell's Equations , 2022 .

[32]  Maciej Zworski,et al.  Resonance expansions of scattered waves , 2000 .

[33]  Thomas Hagstrom,et al.  On generalized discrete PML optimized for propagative and evanescent waves , 2012, 1210.7862.

[34]  Lothar Reichel,et al.  Error Estimates and Evaluation of Matrix Functions via the Faber Transform , 2009, SIAM J. Numer. Anal..

[35]  Brian Davies,et al.  Partial Differential Equations II , 2002 .

[36]  Jan S. Hesthaven,et al.  Spectral Methods for Time-Dependent Problems: Contents , 2007 .

[37]  James Demmel,et al.  Model Reduction for RF MEMS Simulation , 2004, PARA.

[38]  Stefan Hein,et al.  Fano resonances in acoustics , 2010, Journal of Fluid Mechanics.

[39]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[40]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .