2016.53: Compressing variable-coefficient exterior Helmholtz problems via RKFIT
2016.53: Vladimir Druskin, Stefan Güttel and Leonid Knizhnerman (2016) Compressing variable-coefficient exterior Helmholtz problems via RKFIT.
This is the latest version of this eprint.
Full text available as:
|PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader|
The efficient discretization of Helmholtz problems on unbounded domains is a challenging task, in particular, when the wave medium is nonhomogeneous. We present a new numerical approach for compressing finite difference discretizations of such problems, thereby giving rise to efficient perfectly matched layers (PMLs) for nonhomogeneous media. This approach is based on the solution of a nonlinear rational least squares problem using the RKFIT method proposed in [M. Berljafa and S. Guettel, SIAM J. Matrix Anal. Appl., 36(2):894--916, 2015]. We show how the solution of this least squares problem can be converted into an accurate finite difference grid within a rational Krylov framework. Several numerical experiments are included. They indicate that RKFIT computes PMLs more accurate than previous analytic approaches and even works in regimes where the Dirichlet-to-Neumann functions to be approximated are highly irregular. Spectral adaptation effects allow for accurate finite difference grids with point numbers below the Nyquist limit.
|Item Type:||MIMS Preprint|
|Uncontrolled Keywords:||finite difference grid, Helmholtz equation, Dirichlet-to-Neumann map, perfectly matched layer, rational approximation, continued fraction|
|Subjects:||MSC 2000 > 30 Functions of a complex variable|
MSC 2000 > 35 Partial differential equations
MSC 2000 > 65 Numerical analysis
|Deposited By:||Stefan Güttel|
|Deposited On:||08 November 2016|
Available Versions of this Item
- Compressing variable-coefficient exterior Helmholtz problems via RKFIT (deposited 08 November 2016) [Currently Displayed]