A fast volume integral equation method for the direct/inverse problem in elastic wave scattering phenomena

A fast method for solving the volume integral equation is introduced for the solution of forward and inverse multiple scattering problems in an elastic 3-D full space. For both forward and inverse scattering analysis, the volume integral equation in the wavenumber domain is used. By means of the discrete Fourier transform, the volume integral equation in the wavenumber domain can be dealt with as a Fredholm equation of the 2nd kind with respect to a non-Hermitian operator on a finite dimensional vector space. The Bi-CGSTAB method is employed to construct the Krylov subspace in the wavenumber domain. The current procedure establishes a fast and simplified method without requiring the derivation of a coefficient matrix. Several numerical results validate the accuracy and effectiveness of the current method for both forward and inverse scattering analysis. According to the numerical results, the reconstruction of inhomogeneities of the wave field is successful, even for multiple scattering of several cubes.