Energy Dissipative Local Discontinuous Galerkin Methods for Keller–Segel Chemotaxis Model

Document Type


Publication Date



© 2018, Springer Science+Business Media, LLC, part of Springer Nature. 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