“…In particular, if the spatial gradients present in the system are weak, the Navier-Stokes (NS) constitutive equations for the fluxes of mass, momentum, and energy have been derived (with explicit expressions for the transport coefficients) for the model of inelastic hard spheres characterized by constant coefficients of normal restitution α ij . Most of the early derivations were restricted to the quasielastic limit (α ij ≈ 1), thus assuming an expansion around Maxwellians at the same temperature [3,4,5,6,7,8]. However, the nonequipartition of energy becomes significant beyond the quasi-elastic limit, as confirmed by kinetic theory [9,10,11,12], computer simulations [11,13,14,15,16,17,18,19,20,21], and real experiments [21,22,23].…”