Abstract

A standard reaction–diffusion equation consists of two additive terms, a diffusion term and a reaction rate term. The latter term is obtained directly from a reaction rate equation which is itself derived from known reaction kinetics, together with modelling assumptions such as the law of mass action for well-mixed systems. In formulating a reaction–subdiffusion equation, it is not sufficient to know the reaction rate equation. It is also necessary to know details of the reaction kinetics, even in well-mixed systems where reactions are not diffusion limited. This is because, at a fundamental level, birth and death processes need to be dealt with differently in subdiffusive environments. While there has been some discussion of this in the published literature, few examples have been provided, and there are still very many papers being published with Caputo fractional time derivatives simply replacing first order time derivatives in reaction–diffusion equations. In this paper, we formulate clear examples of reaction–subdiffusion systems, based on; equal birth and death rate dynamics, Fisher–Kolmogorov, Petrovsky and Piskunov (Fisher–KPP) equation dynamics, and Fitzhugh–Nagumo equation dynamics. These examples illustrate how to incorporate considerations of reaction kinetics into fractional reaction–diffusion equations. We also show how the dynamics of a system with birth rates and death rates cancelling, in an otherwise subdiffusive environment, are governed by a mass-conserving tempered time fractional diffusion equation that is subdiffusive for short times but standard diffusion for long times.

Highlights

  • Reaction–diffusion partial differential equations are among the most widely used equations in applied mathematics modelling

  • We show how the dynamics of a system with birth rates and death rates cancelling, in an otherwise subdiffusive environment, are governed by a mass-conserving tempered time fractional diffusion equation that is subdiffusive for short times but standard diffusion for long times

  • Over the past two decades there have been large numbers of papers published on numerical methods for nonlinear fractional reaction–diffusion equations

Read more

Summary

Introduction

Reaction–diffusion partial differential equations are among the most widely used equations in applied mathematics modelling These equations govern the time evolution of concentrations, or population densities, of species, at different spatial locations, that are diffusing and reacting. Which can be derived [9,10] from a continuous time random walk (CTRW) [11] with a power law waiting time density Equation (10) can be derived from a CTRW where particles are being removed or added instantaneously at the start of the waiting times between jumps, but only under the contrived. We hope that the examples below will stimulate further activity in this area, where the physical motivation for the modelling equation is stronger

Birth and Death Balance
Fractional Fisher–KPP Equation
Fractional Fitzhugh–Nagumo Equation
Discussion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call