Mathematical analysis and numerical simulation of a fractional reaction-diffusion system with Holling-type III functional response

In recent years, many investigators have questioned the use of convectional diffusion equation to model many physical or real life situations. As a result, fractional space derivatives have been proposed to model anomalous diffusion or related processes, where a particle plume spreads at inconsistent rate with the classical Brownian motion model. By replacing the second derivative in the classical diffusion model with fractional derivative, results to enhance a process known as superdiffusion. A high-dimensional predator-prey reaction-diffusion system with Holling-type III functional response, where the usual second-order derivatives give place to a fractional derivative of order α with 1 < α ≤ 2. Analysis of the main equation guides in the correct choice of parameter values. We established the condition for local and global stabilities. We also show that the system undergoes a Hopf bifurcation subject to a small perturbation of the steady-state solution. The complexity of fractional derivative at some instances of order α for the superdiffusive scenario is demonstrated with some numerical experiments in one, two and three dimensions. The effectiveness of the numerical method is demonstrated through numerical simulations to confirm the theoretical results.