An efficient implementation of a numerical method for a chemotaxis system

  • Authors:
  • Michael W. Smiley

  • Affiliations:
  • Department of Mathematics, Iowa State University, Ames, IA, USA

  • Venue:
  • International Journal of Computer Mathematics - RECENT ADVANCES IN COMPUTATIONAL AND APPLIED MATHEMATICS IN SCIENCE AND ENGINEERING
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

A system of partial differential equations modelling the attraction of a population of cells to a biochemical concentration in its environment is considered. The system incorporates convective and diffusive effects, either of which may dominate. A numerical method is presented that allows for both possible features, which conserves positivity of solutions and their mass in the absence of sources or sinks. An efficient implementation of the method in a two-dimensional setting is described that uses fast Poisson solvers for either 5-or 9-pt stencils approximating the Laplacian. It is also shown that the computational work for either stencil is essentially the same.