Analysis of the Hessian for inverse scattering problems: II. Inverse medium scattering of acoustic waves

We address the inverse problem for scattering of acoustic waves due to an inhomogeneous medium. We derive and analyze the Hessian in both Hölder and Sobolev spaces. Using an integral equation approach based on Newton potential theory and compact embeddings in Hölder and Sobolev spaces, we show that the Hessian can be decomposed into two components, both of which are shown to be compact operators. Numerical examples are presented to validate our theoretical results. The implication of the compactness of the Hessian is that for small data noise and model error, the discrete Hessian can be approximated by a low-rank matrix. This in turn enables fast solution of an appropriately regularized inverse problem, as well as Gaussian-based quantification of uncertainty in the estimated inhomogeneity. (Some figures may appear in colour only in the online journal)

[1]  R. Kress,et al.  Integral equation methods in scattering theory , 1983 .

[2]  Omar Ghattas,et al.  Analysis of the Hessian for inverse scattering problems: I. Inverse shape scattering of acoustic waves , 2012 .

[3]  David Higdon,et al.  Adaptive Hessian-Based Nonstationary Gaussian Process Response Surface Method for Probability Density Approximation with Application to Bayesian Solution of Large-Scale Inverse Problems , 2012, SIAM J. Sci. Comput..

[4]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[5]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[6]  F. B. Hildebrand,et al.  Methods of applied mathematics , 1953 .

[7]  E. Zeidler.  Nonlinear functional analysis and its applications , 1988 .

[8]  Helmut Harbrecht,et al.  Compact gradient tracking in shape optimization , 2008, Comput. Optim. Appl..

[9]  R. Kanwal.  Linear Integral Equations , 1925, Nature.

[10]  Bart G. van Bloemen Waanders,et al.  Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations , 2011, SIAM J. Sci. Comput..

[11]  Laurent Demanet,et al.  Matrix probing: A randomized preconditioner for the wave-equation Hessian , 2011, 1101.3615.

[12]  W. McLean.  Strongly Elliptic Systems and Boundary Integral Equations , 2000 .

[13]  Tobias J. Hagge,et al.  Physics , 1929, Nature.

[14]  A. Kirsch.  An Introduction to the Mathematical Theory of Inverse Problems , 1996, Applied Mathematical Sciences.

[15]  Helmut Harbrecht,et al.  Coupling of FEM and BEM in Shape Optimization , 2006, Numerische Mathematik.

[16]  Todd Arbogast,et al.  Methods of Applied Mathematics , 2008 .

[17]  R. Kress.  Linear Integral Equations , 1989 .

[18]  George Biros,et al.  FaIMS: A fast algorithm for the inverse medium problem with multiple frequencies and multiple sources for the scalar Helmholtz equation , 2012, J. Comput. Phys..

[19]  Ken Kreutz-Delgado,et al.  The Complex Gradient Operator and the CR-Calculus ECE275A - Lecture Supplement - Fall 2005 , 2009, 0906.4835.

[20]  Nonlinear functional analysis and its applications, part I: Fixed-point theorems , 1991 .

[21]  James Martin,et al.  A Stochastic Newton MCMC Method for Large-Scale Statistical Inverse Problems with Application to Seismic Inversion , 2012, SIAM J. Sci. Comput..

[22]  J. Navarro-Pedreño.  Numerical Methods for Least Squares Problems , 1996 .

[23]  Peter Monk,et al.  An analysis of the coupling of finite-element and Nyström methods in acoustic scattering , 1994 .