Energy Dissipative Local Discontinuous Galerkin Methods for Keller–Segel Chemotaxis Model
Document Type
Article
Publication Date
3-1-2019
Department
Department of Mathematical Sciences
Abstract
In this paper, we apply the local discontinuous Galerkin (LDG) method to solve the classical Keller–Segel (KS) chemotaxis model. The exact solution of the KS chemotaxis model may exhibit blow-up patterns with certain initial conditions, and is not easy to approximate numerically. Moreover, it has been proved that there exists a definition of free energy of the KS system which dissipates with respect to time. We will construct a consistent numerical energy and prove the energy dissipation with the LDG discretization. Several numerical experiments in one and two space dimensions will be given. Especially, for solutions with blow-up (converge to Dirac delta functions), the densities of KS model are computed to be strictly positive in the numerical experiments and the energies are also numerically observed to be strictly positive and decreasing as are seen in the figures. Therefore, the scheme is stable for the KS model with blow-up solutions.
Publication Title
Journal of Scientific Computing
Recommended Citation
Guo, L.,
Li, X.,
&
Yang, Y.
(2019).
Energy Dissipative Local Discontinuous Galerkin Methods for Keller–Segel Chemotaxis Model.
Journal of Scientific Computing,
78(3), 1387-1404.
http://doi.org/10.1007/s10915-018-0813-8
Retrieved from: https://digitalcommons.mtu.edu/michigantech-p/4957