Abstract

Reaction–diffusion equations (RDEs) are often derived as continuum limits of lattice-based discrete models. Recently, a discrete model which allows the rates of movement, proliferation and death to depend upon whether the agents are isolated has been proposed, and this approach gives various RDEs where the diffusion term is convex and can become negative (Johnston et al., 2017), i.e. forward–backward–forward diffusion. Numerical simulations suggest these RDEs support shock-fronted travelling waves when the reaction term includes an Allee effect. In this work we formalise these preliminary numerical observations by analysing the shock-fronted travelling waves through embedding the RDE into a larger class of higher order partial differential equations (PDEs). Subsequently, we use geometric singular perturbation theory to study this larger class of equations and prove the existence of these shock-fronted travelling waves. Most notable, we show that different embeddings yield shock-fronted travelling waves with different properties.

Highlights

  • Reaction–diffusion equations (RDEs) are widely used to study population dynamics in cell biology and ecology [1]

  • For RDEs established from the continuum limit of stochastic models, a solution of the RDE shows the macroscopic evolution of U(x, t), but it reflects how microscopic behaviour of individuals influences the macroscopic outcomes [2,3,4,5,6,7]

  • During each time step of duration τ, isolated agents attempt to move to nearest neighbour lattice sites with a probability Pmi, attempt to proliferate to form new agents in neighbour sites with a probability Ppi and to die with a probability Pdi

Read more

Summary

Introduction

Reaction–diffusion equations (RDEs) are widely used to study population dynamics in cell biology and ecology [1]. Considering these different behaviours of isolated and grouped agents, including motility, proliferation and death events, an RDE with a nonlinear diffusivity function and a logistic or Allee type reaction term was derived as the continuum limit. Note that with the nonlinear diffusivity function D(U) centred around 2/3 given by (2) we only observe shock-fronted travelling wave solutions with the strong Allee effect. Note that in this case, GSPT has to be extended since the critical manifold loses normal hyperbolicity near a fold point.

Non-local regularisation
Preliminary observations
Layer problem
The construction of the heteroclinic orbit in the singular limit
Viscous relaxation
Reduced problem
Persistence analysis
Regularisations and the lattice-based stochastic model
Generalisations
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