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

摘要

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.