Published In

Computers & Mathematics with Applications

Document Type

Post-Print

Publication Date

2013

Subjects

Numerical analysis, Chemotaxis, Finite element method

Abstract

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.

Description

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).

doi:10.1016/j.camwa.2013.05.014

DOI

10.1016/j.camwa.2013.05.014

Persistent Identifier

http://archives.pdx.edu/ds/psu/10605

Included in

Mathematics Commons

Share

COinS