“…The Euler-Poisson systems for plasmas have been intensively studied in the literature both analytically and numerically (see, e.g., [2,3,6,11,14]). In the zero-relaxation-time limit, the model equations can be approximated by the driftdiffusion equations, which are much simpler to treat, e.g.…”