This paper extends the volume filling chemotaxis model [18, 26] by taking into account the cell population interactions. Theextended chemotaxis models have nonlinear diffusion and chemotactic sensitivity dependingon cell population density, which is a modification of the classical Keller-Segel model inwhich the diffusion and chemotactic sensitivity are constants (linear). The existence andboundedness of global solutions of these models are discussed and the numerical patternformations are shown. The further improvement is proposed in the end.
