Analysis of a fully discrete approximation for the classical Keller--Segel model: lower and a priori bounds
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2020
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866911774293360640 |
|---|---|
| author | Gutiérrez-Santacreu, Juan Vicente Rodríguez-Galván, José Rafael |
| author_facet | Gutiérrez-Santacreu, Juan Vicente Rodríguez-Galván, José Rafael |
| contents | This paper is devoted to constructing approximate solutions for the classical Keller--Segel model governing \emph{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 chemoattractant.
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 chemoattractant density --, are bounded in the $L^1(Ω)$-norm, satisfy a discrete energy law, and have \emph{ 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. Furthermore, some numerical examples are carried out to support and complement the theoretical results. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2003_08867 |
| institution | arXiv |
| publishDate | 2020 |
| record_format | arxiv |
| spellingShingle | Analysis of a fully discrete approximation for the classical Keller--Segel model: lower and a priori bounds Gutiérrez-Santacreu, Juan Vicente Rodríguez-Galván, José Rafael Numerical Analysis 35K20, 35K55, 65K60 This paper is devoted to constructing approximate solutions for the classical Keller--Segel model governing \emph{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 chemoattractant. 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 chemoattractant density --, are bounded in the $L^1(Ω)$-norm, satisfy a discrete energy law, and have \emph{ 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. Furthermore, some numerical examples are carried out to support and complement the theoretical results. |
| title | Analysis of a fully discrete approximation for the classical Keller--Segel model: lower and a priori bounds |
| topic | Numerical Analysis 35K20, 35K55, 65K60 |
| url | https://arxiv.org/abs/2003.08867 |