A high-order algorithm for time-domain scattering in three dimensions (Q6168050)

Let \(\Omega\) be a bounded domain in \(\mathbb{R}^n\) and \(\Omega^c := \mathbb{R}^n \setminus \Omega\) be its exterior domain. The authors give the numerical methods to solve a linear wave equation in the unbounded domain \([0,T] \times \Omega^c\) with homogenous initial conditions and inhomogeneous boundary conditions such as the Dirichlet, Neumann or Robin boundary conditions. Their strategy is not using the well-known truncation into the bounded domain, but using the surface integral operator with the fundamental solution for the wave equation. Remark that the space-time domain \([0,T] \times \partial \Omega\) for the surface-integration is bounded. First, they derive the surface integral equation, and they derive the fully discrete scheme by discretizing the surface integral equation by decomposing the space-time domain. They apply the spectrally accurate Galerkin method in the discretization with respect to the space variable and the convolution quadratic Runge-Kutta method in the discretization with respect to the time variable, which is a combination of the Laplace transform and the Runge-Kutta method. Lastly, they give numerical examples to verify the method's superiority by comparing the previous results.