Authors: Luke Hatcher, Shannon Horrigan and Zak Lawrence
Supervisors: Marie-Therese Wolfram, Charles Elliott
Our research study group is interested in microorganisms that respond to extracellular changes in their environment by biased random motion towards attractants. Such movement is commonly referred to as chemo-taxis. One example would be the bacterium E.Coli. When the receptors detect a chemo-attractant gradient, the bio-motors are activated in such a way that aggregation occurs .
The Keller-Segel equations can be used to model this behaviour and are given as follows
with Neumann boundary conditions, and specified initial conditions. Here represents the population density and represents the concentration of the chemoattractant.In one dimension it is known that global in time solutions exist, but in two dimensions it is possible for blow-up to occur . For the simulation below a Forward-Euler finite difference method was used where , and . The initial conditions were given by
This leads to blow-up.
We restricted ourselves to considering only the Parabolic-Elliptic equations (i.e. ), and carried out some numerical simulations to investigate stability of solutions as well as proving our numerical scheme is convergent.
We acknowledge and thank the help of our supervisors Dr Marie-Therese Wolfram and Prof Charles Elliott.
We acknowledge the funding body EPSRC and the support from MASDOC CDT.