You are here: MIMS > EPrints
MIMS EPrints

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.

There is a more recent version of this eprint available. Click here to view it.

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
2212 Kb


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. Güttel, 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
MIMS number:2016.53
Deposited By:Stefan Güttel
Deposited On:04 November 2016

Available Versions of this Item

Download Statistics: last 4 weeks
Repository Staff Only: edit this item