“…There, it has been possible to prove that BGK equations retain most of the significant features of the Boltzmann equations: for instance, they preserve the correct collision invariants, conservation equations, and local equilibria (6) with mass action law (7). In addition, such a relaxation-time approximation, in spite of the heavy machinery of computing the auxiliary parameters in terms of the actual fields, turns out to be reliable in comparison to actual deterministic kinetic calculations in space homogeneous conditions [32], and to be effective in dealing with classical space dependent problems like steady shock structure [25] and Riemann problem [33]. However, there are some drawbacks which affect the previous BGK approach, and which confine its application to a partial range of the various physical parameters, on one side, and prevent a deeper mathematical investigation, on the other.…”