Analysis of a fully discrete approximation for the classical Keller–Segel model: lower and a priori bounds
share
This paper is devoted to constructing approximate solutions for the classical Keller–Segel model governing chemotaxis. It consists of a system of nonlinear parabolic equations, where the unknowns are the average density of cells (or organisms), which is a conserved variable, and the average density of chemioattranct. The numerical proposal is made up of a crude finite element method together with a mass lumping technique and a semi-implicit Euler time integration. The resulting scheme turns out to be linear and decouples the computation of variables. The approximate solutions keep lower bounds – positivity for the cell density and nonnegativity for the chemioattranct density –, are mass conservative, satisfy a discrete energy law, and have a priori energy estimates. The latter is achieved by means of a discrete Moser–Trudinger inequality. As far as we know, our numerical method is the first one that can be encountered in the literature dealing with all of the previously mentioned properties at the same time.
READ FULL TEXT