Multifrequency inverse obstacle scattering with unknown impedance boundary conditions using recursive linearization

DOI10.1007/S10444-021-09915-1MaRDI QIDQ2070290zbMATH OpenOpenAlexFDO

Publication date 24 January 2022

Published in Advances in Computational Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2104.13489

boundary integral equations Helmholtz equation impedance boundary condition inverse obstacle problem recursive linearization

Mathematics Subject Classification ID

Integral representations, integral operators, integral equations methods in two dimensions (31A10) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Ill-posed problems for PDEs (35R25) Inverse problems for PDEs (35R30) Inverse problems for integral equations (45Q05) Numerical methods for ill-posed problems for boundary value problems involving PDEs (65N20) Numerical methods for inverse problems for boundary value problems involving PDEs (65N21) Numerical methods for inverse problems for integral equations (65R32)

Abstract: We consider the reconstruction of the shape and the impedance function of an obstacle from measurements of the scattered field at receivers outside the object. The data is assumed to be generated by plane waves impinging on the obstacle from multiple directions and at multiple frequencies. This inverse problem is reformulated as the optimization problem of finding band-limited shape and impedance functions which minimize the

L^{2}

distance between the computed value of the scattered field at the receivers and the data. The optimization problem is non-linear, non-convex, and ill-posed. Moreover, the objective function is computationally expensive to evaluate. The recursive linearization approach (RLA) proposed by Chen has been successful in addressing these issues in the context of recovering the sound speed of a domain or the shape of a sound-soft obstacle. We present an extension of the RLA for the recovery of both the shape and impedance functions. The RLA is a continuation method in frequency where a sequence of single frequency inverse problems is solved. At each higher frequency, one attempts to recover incrementally higher resolution features using a step assumed to be small enough to ensure that the initial guess obtained at the preceding frequency lies in the basin of attraction for Newton's method at the new frequency. We demonstrate the effectiveness of the method with several numerical examples. Surprisingly, we find that one can recover the shape with high accuracy even when the measurements are from sound-hard or sound-soft objects. While the method is effective in obtaining high quality reconstructions for complicated geometries and impedance functions, a number of interesting open questions remain. We present numerical experiments that suggest underlying mechanisms of success and failure, showing areas where improvements could help lead to robust and automatic tools.

Recommendations

Cites work

Cited in

(9)

Describes a project that uses

Uses Software

FaIMS

This page was built for publication: Multifrequency inverse obstacle scattering with unknown impedance boundary conditions using recursive linearization

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2070290)