An asymptotically stable discretization for the Euler–Poisson system in the quasi-neutral limit

We are interested in the modeling of a plasma in the quasi-neutral limit using the Euler–Poisson system. When this system is discretized with a standard numerical scheme, it is subject to a severe numerical constraint related to the quasi-neutrality of the plasma. We propose an asymptotically stable discretization of this system in the quasi-neutral limit. We present numerical simulations for two different one-dimensional test cases that confirm the expected stability of the scheme in the quasi-neutral limit.