Published In

Computers & Mathematics with Applications

Document Type


Publication Date



Numerical analysis, Chemotaxis, Finite element method


We investigate nonnegativity of exact and numerical solutions to a generalized Keller–Segel model. This model includes the so-called “minimal” Keller–Segel model, but can cover more general chemistry. We use maximum principles and invariant sets to prove that all components of the solution of the generalized model are nonnegative. We then derive numerical methods, using finite element techniques, for the generalized Keller–Segel model. Adapting the ideas in our proof of nonnegativity of exact solutions to the discrete setting, we are able to show nonnegativity of discrete solutions from the numerical methods under certain standard assumptions. One of the numerical methods is then applied to the minimal Keller–Segel model. Recalling known results on the qualitative behavior of this model, we are able to choose parameters that yield convergence to a nonhomogeneous stationary solution. While proceeding to exhibit these stationary patterns, we also demonstrate how naive choices of numerical methods can give physically unrealistic solutions, thereby justifying the need to study positivity preserving methods.


Copyright © 2013 Elsevier Ltd. All rights reserved.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.


NOTICE: this is the author’s version of a work that was accepted for publication in Computers & Mathematics with Applications. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Computers & Mathematics With Applications, 66(3), 356-375 (2013).



Persistent Identifier

Included in

Mathematics Commons