Asymptotic-preserving and well-balanced schemes for the 1D Cattaneo model of chemotaxis movement in both hyperbolic and diffusive regimes (Q663648): Difference between revisions

In this paper the original Well-Balanced (WB) framework of \textit{J. M. Greenberg} and \textit{A.-Y. Le Roux} [SIAM J. Numer. Anal. 33, No. 1, 1--16 (1996; Zbl 0876.65064)] and \textit{L. Gosse} [Math. Comput. 71, No. 238, 553--582 (2002; Zbl 0997.65108)] relying on Non-Conservative (NC) productsis set up in order to efficiently treat the so-called Cattaneo model of chemotaxis in 1D. It proceeds by concentrating the source terms onto Dirac masses: this allows to handle them by NC jump relations based on steady-state equations which can be integrated explicitly. A Riemann solver is deduced and the corresponding WB Godunov scheme completed with the standard Hoff-Smoller theory for the diffusion-reaction equation ruling the evolution of the chemoattractant concentration is studied in detail. Later, following former results of \textit{L. Gosse} and \textit{G. Toscani} [C. R., Math., Acad. Sci. Paris 334, No. 4, 337--342 (2002; Zbl 0996.65093)] and, a simple rewriting of the NC jump relations allows to generate another version of the same Godunov scheme which is well adapted to the parabolic scaling involving a small parameter \(\epsilon \). The standard BV framework is used to study the uniform stability of this Asymptotic-Preserving (AP) scheme with respect to \(\epsilon \) allows to pass to the limit and derive a simple centered discretization of the Keller-Segel model. Finally, results by \textit{F. Filbet} [Numer. Math. 104, No. 4, 457--488 (2006; Zbl 1098.92006)] permit to pass to the complementary limit when the space-step \(h\) is sent to zero. Numerical results are included to illustrate the feasibility and the efficiency of the method.