DEPARTMENT OF MATHEMATICAL SCIENCES

APPLIED AND COMPUTATIONAL MATHEMATICS SEMINAR

**Speaker:** Yekaterina Epshteyn, Mathematics, Carnegie Mellon University.

**Title: ***
New Discontinuous Galerkin Methods for the Chemotaxis Model and Closely Related Biomedical Problems
*

**Abstract:**
In this work, first, we propose a family of new interior penalty discontinuous Galerkin methods for the Keller-Segel chemotaxis model. This model is described by a system of two
nonlinear PDEs: a convection-diffusion equation for the cell density coupled with a reaction-diffusion equation for the chemoattractant concentration. It
has been recently shown that the convective part of this system is of a mixed hyperbolic-elliptic type, which may cause severe instabilities when the
studied system is solved by straightforward numerical methods. Therefore, the first step in the derivation of our new methods is made by introducing the new
variable for the gradient of the chemoattractant concentration and by reformulating the original Keller-Segel model in the form of a
convection-diffusion-reaction system with a hyperbolic convective part. We then design interior penalty discontinuous Galerkin methods for the rewritten
Keller-Segel system. Our methods employ the central-upwind numerical fluxes, originally developed in the context of finite-volume methods for hyperbolic
systems of conservation laws.

We consider Cartesian grids and prove error estimates for the proposed high-order discontinuous Galerkin methods. Our proof is valid for
pre-blow-up times since we assume boundedness of the exact solution. We also show that the blow-up time of the exact solution is bounded from above by the
blow-up time of our numerical solution.

In the numerical tests that will be presented,
we first compare three different discontinuous Galerkin schemes applied to the Keller-Segel model:

1) primal discontinuous Galerkin methods applied to the original formulation of the Keller-Segel model,

2) primal discontinuous Galerkin methods with the standard upwind numerical fluxes for the reformulated Keller-Segel model and,

3) the new discontinuous Galerkin methods.

We show, that compare to the new discontinuous Galerkin methods, the first two schemes fail to give the accurate, oscillation free solutions for the classical Keller-Segel chemotaxis model.\\
Second, we consider application of the new discontinuous Galerkin schemes, to some biomedical problems.

This work is based upon collaboration with Alexander Kurganov, Tulane University.

**Time:** Friday, February 29, 2008, 1:30-2:30 p.m.

**Place:** Research I, Room 301

Department of Mathematical Sciences

George Mason University

4400 University Drive, MS 3F2

Fairfax, VA 22030-4444

http://math.gmu.edu/

Tel. 703-993-1460, Fax. 703-993-1491