Paper detail

Dirichlet-Neumann Waveform Relaxation Method for the 1D and 2D Heat and Wave Equations in Multiple subdomains

We present a Waveform Relaxation (WR) version of the Dirichlet-Neumann algorithm, formulated specially for multiple subdomains splitting for general parabolic and hyperbolic problems. This method is based on a non-overlapping spatial domain decomposition, and the iteration involves subdomain solves in space-time with corresponding interface condition, and finally organize an exchange of information between neighboring subdomains. Using a Fourier-Laplace transform argument, for a particular relaxation parameter, we present convergence analysis of the algorithm for the heat and wave equations. We prove superlinear convergence for finite time window in case of the heat equation, and finite step convergence for the wave equation. The convergence behavior however depends on the size of the subdomains and the time window length on which the algorithm is employed. We illustrate the performance of the algorithm with numerical results, and show a comparison with classical and optimized Schwarz WR methods.