Abstract
We study the inverse boundary value problem for the Helmholtz equation using the Dirichlet-to-Neumann map at selected frequencies as the data. A conditional Lipschitz stability estimate for the inverse problem holds in the case of wavespeeds that are a linear combination of piecewise constant functions (following a domain partition) and gives a framework in which the scheme converges. The stability constant grows exponentially as the number of subdomains in the domain partition increases. We establish an order optimal upper bound for the stability constant. We eventually realize computational experiments to demonstrate the stability constant evolution for three-dimensional wavespeed reconstruction.
Original language | English (US) |
---|---|
Pages (from-to) | 3962-3983 |
Number of pages | 22 |
Journal | SIAM Journal on Mathematical Analysis |
Volume | 48 |
Issue number | 6 |
DOIs | |
State | Published - 2016 |
Keywords
- Helmholtz equation
- Inverse problems
- Stability and convergence of numerical methods
ASJC Scopus subject areas
- Analysis
- Computational Mathematics
- Applied Mathematics