Diffusive limit of a two-dimensional well-balanced approximation to a kinetic model of chemotaxis

A (2+2)-dimensional kinetic equation, directly inspired by the run-and-tumble modelingof chemotaxis dynamics is studied so as to derive a both ''2D well-balanced'' and''asymptotic-preserving'' numerical approximation. To this end, exact stationary regimesare expressed by means of Laplace transforms of Fourier-Bessel solutions of associatedelliptic equations. This yields a scattering S-matrix which permits to formulate a timemarchingscheme in the form of a convex combination in kinetic scaling. Then, in thediffusive scaling, an IMEX-type discretization follows, for which the ''2D well-balancedproperty'' still holds, while the consistency with the asymptotic drift-diffusion equation ischecked. Numerical benchmarks, involving ''nonlocal gradients'' (or finite samplingradius), carried out in both scalings, assess theoretical findings. Nonlocal gradients appearto inhibit blowup phenomena.