We will discuss the numerical simulation of acoustic wave propagation with localized inhomogeneities. To do this we will apply a standard finite element method (FEM) in space and explicit time-stepping in time on a finite spatial domain containing the inhomogeneities. The equations in the exterior computational domain will be dealt with a time-domain boundary integral formulation discretized by the Galerkin boundary element method (BEM) in space and convolution quadrature in time.
\\
\\
We will give the analysis of the proposed method, starting with the proof of a positivity preservation property of convolution quadrature as a consequence of a variant of the Herglotz theorem. Combining this result with standard energy analysis of leap-frog discretization of the interior equations will give us both stability and convergence of the method. Numerical results will also be given.