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

Journal of Computational and Applied Mathematics, vol. 485, pp. 117587

Abstract

A novel fourth-order, L-stable, exponential time differencing Runge–Kutta type scheme is developed to solve non-linear systems of reaction–diffusion equations, particularly with non-smooth 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. The L-acceptable rational approximation ensures efficient damping of spurious oscillations arising from non-smooth initial and/or boundary conditions. The real and distinct poles of the rational function eliminate the presence of complex arithmetic in the solution of linear systems and thus make the scheme more attractive for parallelization in low-level languages like Fortran and C++. 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 when implemented in parallel, with up to five times speed-up in CPU time.

Authors 3

1. W.K. Attipoe

   Clarkson University

   Affiliation as printed

   Clarkson University, 8 Clarkson Ave, Potsdam, 13676, NY, USA

2. Andreas Kleefeld Aachen Faculty of Medical Engineering and Technomathematics

   FH Aachen · Forschungszentrum Jülich · Jülich Supercomputing Centre

   Affiliation as printed

   Faculty of Medical Engineering and Technomathematics, University of Applied Sciences, Aachen, 52428, Jülich, Germany

   Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH, 52425, Jülich, Germany

3. Emmanuel Owusu Asante-Asamani corresponding

   Clarkson University

   Affiliation as printed

   Clarkson University, 8 Clarkson Ave, Potsdam, 13676, NY, USA