Dirichlet-Neumann and Neumann-Neumann Waveform Relaxation Methods for PDEs with Time Delay
Deeksha Tomer, Bankim Chandra Mandal
TL;DR
This work extends Dirichlet-Neumann Waveform Relaxation (DNWR) and Neumann-Neumann Waveform Relaxation (NNWR) to PDEs with time delay, analyzing parabolic, hyperbolic, and neutral models. By applying Laplace-transform techniques to the error equations, it derives explicit convergence factors: for symmetric two-subdomain partitions DNWR converges in two iterations at $\theta=\tfrac{1}{2}$ and NNWR at $\theta=\tfrac{1}{4}$, with linear convergence for other $\theta$ values; analogous results hold for hyperbolic and neutral PDEs, with corresponding expressions $\hat{h}^k(s)=(1-2\theta)^k\hat{h}^0(s)$ or $\hat{g}^k(s)=(1-4\theta)^k\hat{g}^0(s)$. Numerical illustrations for reaction-diffusion, wave, and neutral PDEs corroborate the theory and demonstrate effective parallelization, including multi-subdomain scenarios, and show superiority over classical/optimized Schwarz methods for optimal $\theta$ choices. The findings establish DNWR and NNWR as efficient, scalable tools for time-delayed PDEs, offering rapid convergence and parallel performance across a range of models. This contributes a practical framework for solving Delay PDEs in large-scale computations where interface-driven iteration is advantageous.
Abstract
We introduce and compare two domain decomposition based numerical methods, namely the Dirichlet-Neumann and Neumann-Neumann Waveform Relaxation methods (DNWR and NNWR respectively), tailored for solving partial differential equations (PDEs) incorporating time delay. Time delay phenomena frequently arise in various real-world systems, making their accurate modeling and simulation crucial for understanding and prediction. We consider a series of model problems, ranging from Parabolic, Hyperbolic to Neutral PDEs with time delay and apply the iterative techniques DNWR and NNWR for solving in parallel. We present the theoretical foundations, numerical implementation, and comparative performance analysis of these two methods. Through numerical experiments and simulations, we explore their convergence properties, computational efficiency, and applicability to various types of PDEs with time delay.
