A domain integral method employing a specific Green's function (i.e., incorporating some features of the global problem of wave propagation in an inhomogeneous medium) is developed for solving direct and inverse scattering problems relative to slab-like macroscopically inhomogeneous porous obstacles. It is shown how to numerically solve such problems, involving both spatially-varying density and compressibility, by means of an iterative scheme initialized with a Born approximation. A numerical solution is obtained for a canonical problem involving a two-layer slab.
IntroductionThis work was initially motivated by two problems: i) the design problem connected with the determination of the optimal profile of a continuous and/or discontinuous spatial distribution of the material/geometric properties of porous materials for the absorption of sound [7] and ii) the retrieval of the spatially-varying mechanical and geometrical parameters of bone for the diagnosis of diseases such as osteoporosis [10].Such inverse problems [20] can be decomposed into two sub-problems: i) the determination of the constitutive and conservation relations linking the various spatially-variable mechanical parameters of the porous medium to its response to an acoustic solicitation, and ii) the resolution of the * Correspondence to: J.-P. Groby, Akoestieke en Thermische Fysica, KULeuven, Celestijnenlaan 200D , 3001 Heverlee, Belgium † E-mail: jeanphilippe.groby@fys.kuleuven.be 1 wave equation in an inhomogeneous porous medium (for instance, within the Biot, or rigid frame approximations). Here we focus on the second point.In [17], it is shown that the wave equation describing the propagation in a macroscopicallyinhomogeneous porous medium in the rigid frame approximation can formally take the form of the usual acoustic wave equation in a macroscopically-inhomogeneous fluid (in which the microscopic features of the porous medium are homogenized) with spatial (and frequency) dependent compressibility κ e (x, ω) and density ρ e (x, ω).The present work deals with a method of resolution of direct problems involving acoustic wave propagation in a macroscopically-inhomogeneous fluid medium, whose density and compressibility are both space dependent, this being a prerequisite to the resolution of related inverse problems.This topic is also of great interest in quantum physics (inverse potential scattering [2,37,38,39]), ocean acoustics [9,8,32,11,36,12] (detection of inhomogeneities, sediment exploration, influence of seawater and seafloor composition and heterogenity on the long-range propagation of acoustic waves in the sea, ...), seismology [1,40,42] (determination of the internal structure and composition of the Earth via seismic waves,...), geophysics [42,25,48,44] (characterization of soil, detection of geological features such as hydrocarbon reservoirs, ...), optics and electromagnetism [47,41] (design and characterization of materials having specified response to waves, detection of flaws,...).The wave equation in an inhomogeneous medi...