Compressed absorbing boundary conditions for the Helmholtz equation

Absorbing layers are sometimes required to be impractically thick in order to offer an accurate approximation of an absorbing boundary condition for the Helmholtz equation in a heterogeneous medium. It is always possible to reduce an absorbing layer to an operator at the boundary by layer-stripping elimination of the exterior unknowns, but the linear algebra involved is costly. We propose to bypass the elimination procedure, and directly fit the surface-to-surface operator in compressed form from a few exterior Helmholtz solves with random Dirichlet data. We obtain a concise description of the absorbing boundary condition, with a complexity that grows slowly (often, logarithmically) in the frequency parameter. We then obtain a fast (nearly linear in the dimension of the matrix) algorithm for the application of the absorbing boundary condition using partitioned low rank matrices. The result, modulo a precomputation, is a fast and memory-efficient compression scheme of an absorbing boundary condition for the Helmholtz equation.

Department dc:contributor.department

Massachusetts Institute of Technology. Department of Mathematics.

Grantor dc:publisher

Massachusetts Institute of Technology

Year dc:date.issued

2014

Author dc:creator

• Bélanger-Rioux, Rosalie

Advisor dc:contributor.advisor

• Laurent Demanet.

• Mathematics.

Statement dc:rights

• M.I.T. theses are protected by copyright. They may be viewed from this source for any purpose, but reproduction or distribution in any format is prohibited without written permission. See provided URL for inquiries about permission.

Licence dc:rights.uri

• http://dspace.mit.edu/handle/1721.1/7582

Language dc:language.iso

eng

Handle dc:identifier.uri

http://hdl.handle.net/1721.1/90179

OAI identifier oai:identifier

oai:dspace.mit.edu:1721.1/90179