A fourth-order exponential time differencing scheme with real and distinct poles rational approximation for solving non-linear reaction-diffusion systems

A fourth-order, L-stable, exponential time differencing Runge-Kutta type scheme is developed to solve nonlinear systems of reaction diffusion equations with nonsmooth data. The new scheme, ETDRK4RDP, is constructed by approximating the matrix exponentials in the ETDRK4 scheme with a fourth order, L-acceptable, non-Padé rational function having real and distinct poles (RDP). Using RDP rational functions to construct the scheme ensures efficient damping of spurious oscillations arising from non-smooth initial and boundary conditions and a straightforward parallelization. We verify empirically that the new ETDRK4RDP scheme is fourth-order accurate for several reaction diffusion systems with Dirichlet and Neumann boundary conditions and show it to be more efficient than competing exponential time differencing schemes, especially when implemented in parallel, with up to six times speedup in CPU time.