An almost fourth order uniformly convergent domain decomposition method for a coupled system of singularly perturbed reaction–diffusion equations

We consider a system of M ( ≥ 2 ) singularly perturbed equations of reaction–diffusion type coupled through the reaction term. A high order Schwarz domain decomposition method is developed to solve the system numerically. The method splits the original domain into three overlapping subdomains. On two boundary layer subdomains we use a compact fourth order difference scheme on a uniform mesh while on the interior subdomain we use a hybrid scheme on a uniform mesh. We prove that the method is almost fourth order ε -uniformly convergent. Furthermore, we prove that when ε is small, one iteration is sufficient to get almost fourth order ε -uniform convergence. Numerical experiments are performed to support the theoretical results.