Fast high-order method for multi-dimensional space-fractional reaction–diffusion equations with general boundary conditions

Abstract

To achieve the efficient and accurate long-time integration, we propose a fast and stable high-order numerical method for solving fractional-in-space reaction–diffusion equations. The proposed method is explicit in nature and utilizes the fourth-order compact finite difference scheme and matrix transfer technique (MTT) in space with FFT-based implementation. Time integration is done through the modified fourth-order exponential time differencing Runge–Kutta scheme. The linear stability analysis and various numerical experiments including two-dimensional (2D) Fitzhugh–Nagumo, Allen–Cahn, Gierer–Meinhardt, Gray–Scott and three-dimensional (3D) Schnakenberg models are presented to demonstrate the accuracy, efficiency, and stability of the proposed method.