Fast solution algorithm for a three-dimensional inverse multifrequency problem of scalar acoustics with data in a cylindrical domain

A new algorithm for stable solution of a three-dimensional scalar inverse problem of acoustic sensing of an inhomogeneous medium in a cylindrical domain is proposed. Data for its solution is the complex amplitude of the wave field measured outside the acoustic inhomogeneities in the cylindrical layer. With the help of the Fourier transform and Fourier series, the inverse problem is reduced to a set of one-dimensional Fredholm integral equations of the first kind. Next, the complex amplitude of the wave field is computed in the inhomogeneity region and the desired sonic velocity field is found in this region. When run on a moderate-performance personal computer, the algorithm takes tens of seconds to solve the inverse problem on rather fine three-dimensional grids. The accuracy of the algorithm is analyzed numerically as applied to test inverse problems at different frequencies, and the stability of the algorithm with respect to data perturbations is investigated.