An asymptotic preserving scheme for front propagation in a kinetic reaction-transport equation
This work addresses the computational challenge of simulating front propagation in kinetic equations for researchers in applied mathematics and physics, offering an efficient numerical method for stiff regimes.
The authors propose an asymptotic preserving scheme for a kinetic reaction-transport equation with KPP-type nonlinearity, enabling efficient simulation of front propagation in the sharp interface regime. Numerical tests demonstrate the scheme's stability, discrete maximum principle, and ability to handle singular behaviors.
In this work, we propose an asymptotic preserving scheme for a non-linear kinetic reaction-transport equation, in the regime of sharp interface. With a non-linear reaction term of KPP-type, a phenomenon of front propagation has been proved in [9]. This behavior can be highlighted by considering a suitable hyperbolic limit of the kinetic equation, using a Hopf-Cole transform. It has been proved in [6, 8, 11] that the logarithm of the distribution function then converges to the viscosity solution of a constrained Hamilton-Jacobi equation. The hyperbolic scaling and the Hopf-Cole transform make the kinetic equation stiff. Thus, the numerical resolution of the problem is challenging, since the standard numerical methods usually lead to high computational costs in these regimes. The Asymptotic Preserving (AP) schemes have been typically introduced to deal with this difficulty, since they are designed to be stable along the transition to the macroscopic regime. The scheme we propose is adapted to the non-linearity of the problem, enjoys a discrete maximum principle and solves the limit equation in the sense of viscosity. It is based on a dedicated micro-macro decomposition, attached to the Hopf-Cole transform. As it is well adapted to the singular limit, our scheme is able to cope with singular behaviors in space (sharp interface), and possibly in velocity (concentration in the velocity distribution). Various numerical tests are proposed, to illustrate the properties and the efficiency of our scheme.