In this work we quantitatively assess, via instabilities, a Navier–Stokes-order (small-Knudsen-number) continuum model based on the kinetic theory analogy and applied to inelastic spheres in a homogeneous cooling system. Dissipative collisions are known to give rise to instabilities, namely velocity vortices and particle clusters, for sufficiently large domains. We compare predictions for the critical length scales required for particle clustering obtained from transient simulations using the continuum model with molecular dynamics (MD) simulations. The agreement between continuum simulations and MD simulations is excellent, particularly given the presence of well-developed velocity vortices at the onset of clustering. More specifically, spatial mapping of the local velocity-field Knudsen numbers ( $K{n}_{u} $ ) at the time of cluster detection reveals $K{n}_{u} \gg 1$ due to the presence of large velocity gradients associated with vortices. Although kinetic-theory-based continuum models are based on a small- $Kn$ (i.e. small-gradient) assumption, our findings suggest that, similar to molecular gases, Navier–Stokes-order (small- $Kn$ ) theories are surprisingly accurate outside their expected range of validity.

A moment method with closures based on Gaussian quadrature formulas is proposed to solve the Boltzmann kinetic equation with a hard-sphere collision kernel for mono-dispersed particles. Different orders of accuracy in terms of the moments of the velocity distribution function are considered, accounting for moments up to seventh order. Quadrature-based closures for four different models for inelastic collision-the Bhatnagar-Gross-Krook, ES-BGK, the Maxwell model for hard-sphere collisions, and the full Boltzmann hard-sphere collision integral-are derived and compared. The approach is validated studying a dilute non-isothermal granular flow of inelastic particles between two stationary Maxwellian walls. Results obtained from the kinetic models are compared with the predictions of molecular dynamics (MD) simulations of a nearly equivalent system with finite-size particles. The influence of the number of quadrature nodes used to approximate the velocity distribution function on the accuracy of the predictions is assessed. Results for constitutive quantities such as the stress tensor and the heat flux are provided, and show the capability of the quadrature-based approach to predict them in agreement with the MD simulations under dilute conditions.