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

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Attipoe, Wisdom Kwame, Kleefeld, Andreas, Asante-Asamani, Emmanuel
Format: Preprint
Veröffentlicht: 2025
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866918079141773312
author Attipoe, Wisdom Kwame
Kleefeld, Andreas
Asante-Asamani, Emmanuel
author_facet Attipoe, Wisdom Kwame
Kleefeld, Andreas
Asante-Asamani, Emmanuel
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2507_01245
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A fourth-order exponential time differencing scheme with real and distinct poles rational approximation for solving non-linear reaction-diffusion systems
Attipoe, Wisdom Kwame
Kleefeld, Andreas
Asante-Asamani, Emmanuel
Numerical Analysis
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.
title A fourth-order exponential time differencing scheme with real and distinct poles rational approximation for solving non-linear reaction-diffusion systems
topic Numerical Analysis
url https://arxiv.org/abs/2507.01245