MathSciDoc: An Archive for Mathematician ∫

In a previous work [P. Crispel, P. Degond, and M.-H. Vignal, J. Comput. Phys., 223 (2007), pp. 208–234], a new numerical discretization of the Euler–Poisson system was proposed. This scheme is “asymptotic preserving” in the quasineutral limit (i.e., when the Debye length $\varepsilon$ tends to zero), which means that it becomes consistent with the limit model when $\varepsilon \to 0$. In the present work, we show that the stability domain of the present scheme is independent of $\varepsilon$. This stability analysis is performed on the Fourier transformed (with respect to the space variable) linearized system. We show that the stability property is more robust when a space-decentered scheme is used (which brings in some numerical dissipation) rather than a space-centered scheme. The linearization is first performed about a zero mean velocity and then about a nonzero mean velocity. At the various stages of the analysis, our scheme is compared with more classical schemes and its improved stability property is outlined. The analysis of a fully discrete (in space and time) version of the scheme is also given. Finally, some considerations about a model nonlinear problem, the Burgers–Poisson problem, are also discussed.

stiffness, Debye length, electron plasma period, Burgers–Poisson, sheath problem, Klein–Gordon