Fourier Spectral Methods with Exponential Time Differencing for Space-Fractional Partial Differential Equations in Population Dynamics

Abstract

Physical laws governing population dynamics are generally expressed as differential equations. Research in recent decades has incorporated fractional-order (non-integer) derivatives into differential models of natural phenomena, such as reaction-diffusion systems. In this paper, we develop a method to numerically solve a multi-component and multi-dimensional space-fractional system. For space discretization, we apply a Fourier spectral method that is suited for multidimensional PDE systems. Efficient approximation of time-stepping is accomplished with a locally one dimensional exponential time differencing approach. We show the effect of different fractional parameters on growth models and consider the convergence, stability, and uniqueness of solutions, as well as the biological interpretation of parameters and boundary conditions.