A non-iterative domain decomposition time integrator for linear wave equations
This work provides an incremental improvement for computational scientists and engineers dealing with wave propagation simulations, adapting existing methods from parabolic problems to wave equations.
The authors tackled the problem of efficiently solving linear acoustic wave equations by proposing a non-iterative domain decomposition integrator that enables parallelization across space without iterations per time step, achieving second-order accuracy in time and global convergence of O(h + τ²) with larger time steps than explicit schemes.
We propose and analyze a non-iterative domain decomposition integrator for the linear acoustic wave equation. The core idea is to combine an implicit Crank-Nicolson step on spatial subdomains with a local prediction step at the subdomain interfaces. This enables parallelization across space while advancing sequentially in time, without requiring iterations at each time step. The method is similar to the methods from Blum, Lisky and Rannacher (1992) or Dawson and Dupont (1992), which have been designed for parabolic problems. Our approach adapts them to the case of the wave equation in a fully discrete setting, using linear finite elements with mass lumping. Compared to explicit schemes, our method permits significantly larger time steps and retains high accuracy. We prove that the resulting method achieves second-order accuracy in time and global convergence of order $\mathcal{O}(h + Ï^2)$ under a CFL-type condition, which depends on the overlap width between subdomains. We conclude with numerical experiments which confirm the theoretical results.