Fractional reaction-diffusion equation for species growth and dispersal
نویسنده
چکیده
Reaction-diffusion equations are commonly used to model the growth and spreading of biological species. The classical diffusion term implies a Gaussian dispersal kernel in the corresponding integro-difference equation, which is often unrealistic in practice. In this paper, we propose a fractional reaction-diffusion equation where the classical second derivative diffusion term is replaced by a fractional derivative of order less than two. The resulting model captures the faster spreading rates and power law invasion profiles observed in many applications, and is strongly motivated by a generalised central limit theorem for random movements with power-law probability tails. We then develop practical numerical methods to solve the fractional reactiondiffusion equation by time discretization and operator splitting, along with Partially supported by NSF grants DMS-0139927 and DMS-0417869 and the Marsden Fund administered by the Royal Society of New Zealand. Boris Baeumer Department of Mathematics and Statistics, University of Otago, P.O. Box 56, Dunedin 9001, New Zealand Tel.: +643-479-7763 Fax: +643-479-8427 E-mail: [email protected] Mihály Kovács Department of Mathematics and Statistics, University of Otago, P.O. Box 56, Dunedin 9001, New Zealand Tel.: +643-479-7889 Fax: +643-479-8427 E-mail: [email protected] Mark M. Meerschaert Department of Mathematics and Statistics, University of Otago, P.O. Box 56, Dunedin 9001, New Zealand Tel.: +643-479-7889 Fax: +643-479-8427 E-mail: [email protected] 2 Baeumer, Kovács, Meerschaert some existing methods from the literature on anomalous super-diffusion. In the process, we establish the mathematical relationship between the discrete time integro-difference and continuous time reaction-diffusion analogues of the model, along with error bounds. Our general approach also applies to other alternative non-Gaussian dispersal kernels, and it identifies the analogous continuous time evolution equations for those models.
منابع مشابه
Fractional reproduction-dispersal equations and heavy tail dispersal kernels.
Reproduction-Dispersal equations, called reaction-diffusion equations in the physics literature, model the growth and spreading of biological species. Integro-Difference equations were introduced to address the shortcomings of this model, since the dispersal of invasive species is often more widespread than what the classical RD model predicts. In this paper, we extend the RD model, replacing t...
متن کاملNumerical Solution of Caputo-Fabrizio Time Fractional Distributed Order Reaction-diffusion Equation via Quasi Wavelet based Numerical Method
In this paper, we derive a novel numerical method to find out the numerical solution of fractional partial differential equations (PDEs) involving Caputo-Fabrizio (C-F) fractional derivatives. We first find out the approximation formula of C-F derivative of function tk. We approximate the C-F derivative in time with the help of the Legendre spectral method and approximation formula o...
متن کاملSpreading speed, traveling waves, and minimal domain size in impulsive reaction-diffusion models.
How growth, mortality, and dispersal in a species affect the species' spread and persistence constitutes a central problem in spatial ecology. We propose impulsive reaction-diffusion equation models for species with distinct reproductive and dispersal stages. These models can describe a seasonal birth pulse plus nonlinear mortality and dispersal throughout the year. Alternatively, they can desc...
متن کاملA numerical investigation of a reaction-diffusion equation arises from an ecological phenomenon
This paper deals with the numerical solution of a class of reaction diffusion equations arises from ecological phenomena. When two species are introduced into unoccupied habitat, they can spread across the environment as two travelling waves with the wave of the faster reproducer moving ahead of the slower.The mathematical modelling of invasions of species in more complex settings that include ...
متن کاملNumerical solutions for fractional reaction-diffusion equations
Fractional diffusion equations are useful for applications where a cloud of particles spreads faster than the classical equation predicts. In a fractional diffusion equation, the second derivative in the spatial variable is replaced by a fractional derivative of order less than two. The resulting solutions spread faster than the classical solutions and may exhibit asymmetry, depending on the fr...
متن کامل