Dirichlet–Neumann waveform relaxation methods for parabolic and hyperbolic problems in multiple subdomains

Abstract

In this paper, a new waveform relaxation variant of the Dirichlet–Neumann algorithm is introduced for general parabolic problems as well as for the second-order wave equation for decompositions with multiple subdomains. The method is based on a non-overlapping decomposition of the domain in space, and the iteration involves subdomain solves in space-time with transmission conditions of Dirichlet and Neumann type to exchange information between neighboring subdomains. Regarding the convergence of the algorithm, two main results are obtained when the time window is finite: for the heat equation, the method converges superlinearly, whereas for the wave equation, it converges after a finite number of iterations. The analysis is based on Fourier–Laplace transforms and detailed kernel estimates, which reveals the precise dependence of the convergence on the size of the subdomains and the time window length. Numerical experiments are presented to illustrate the performance of the algorithm and to compare its convergence behaviour with classical and optimized Schwarz Waveform Relaxation methods. Experiments involving heterogeneous coefficients and non-matching time grids, which are not covered by the theory, are also presented.