Asymptotic high order mass-preserving schemes for a hyperbolic model of chemotaxis

We introduce a new class of finite difference schemes for approximating the solutions to an initial-boundary value problem on a bounded interval for a one-dimensional dissipative hyperbolic system with an external source term, which arises as a simple model of chemotaxis. Since the solutions to this problem may converge to nonconstant asymptotic states for large times, standard schemes usually fail to yield a good approximation. Therefore, we propose a new class of schemes, which use an asymptotic higher order correction, second and third order in our examples, to balance the effects of the source term and the influence of the asymptotic solutions. Special care is needed to deal with boundary conditions to avoid harmful loss of mass. Convergence results are proved for these new schemes, and several numerical tests are presented and discussed to verify the effectiveness of their behavior.