“…Hyperbolic models eliminate this undesirable feature, see [18,22,43]; discrete velocity kinetic models can be handled by means of similar techniques [27,28,40]. Advances in numerical analysis allow to simulate efficiently much more detailed kinetic models [2,16,23,41,42,48]: consult for instance [11,50]. Finally, we stress again the similarity between this class of distinguished numerical schemes and the ones derived for the nonlinear parabolic equations by extending the ideas of the classic Scharfetter-Gummel scheme [44] for which numerical fluxes are computed by solving a Dirichlet problem for the stationary equations, see [20].…”