“…As formal theoretical asymptotic results are obtained by a moment method , assuming that the distribution function of each species i is at a local Maxwellian state with a small velocity for any ( t , x ) ∈ ℝ + × Ω, we apply the same approach in order to derive a numerical scheme which nicely behaves when ϵ tends to 0. More precisely, the results are obtained under the following ansatz, for any 1 ≤ i ≤ p , The same ansatz is also made on the initial condition .…”