GEORGE MASON UNIVERSITY
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