On a chemotaxis model with saturated chemotactic flux

Abstract

We propose a PDE chemotaxis model, which can be viewed as a regularization of the Patlak-Keller-Segel (PKS)system. Our modification is based on a fundamental physical property of the chemotactic flux function---itsboundedness. This means that the cell velocity is proportional to the magnitude of the chemoattractant gradientonly when the latter is small, while when the chemoattractant gradient tends to infinity the cell velocitysaturates. Unlike the original PKS system, the solutions of the modified model do not blow up in either finiteor infinite time in any number of spatial dimensions, thus making it possible to use bounded spiky steady statesto model cell aggregation. After obtaining local and global existence results, we use the local and globalbifurcation theories to show the existence of one-dimensional spiky steady states; we also study the stabilityof bifurcating steady states. Finally, we numerically verify these analytical results, and then demonstrate thatsolutions of the two-dimensional model with nonlinear saturated chemotactic flux function typically develop verycomplicated spiky structures.