Numerical simulation of chemotaxis models on stationary surfaces

In this paper we present an implicit finite element method for a class of chemotaxis models, where a new linearized flux-corrected transport (FCT) algorithm is modified in such a way as to keep the density of on-surface living cells nonnegative. Level set techniques are adopted for an implicit description of the surface and for the numerical treatment of the corresponding system of partial differential equations. The presented scheme is able to deliver a robust and accurate solution for a large class of chemotaxis-driven models. The numerical behavior of the proposed scheme is tested on the blow-up model on a sphere and an ellipsoid and on the pattern-forming dynamics model of Escherichia coli on a sphere.