The kinetic energy distribution function satisfying the Boltzmann equation is studied analytically and numerically for a system of inelastic hard spheres in the case of binary collisions. Analytically, this function is shown to have a similarity form in the simple cases of uniform or steady-state flows. This determines the region of validity of hydrodynamic description. The latter is used to construct the phase diagram of granular systems, and discriminate between clustering instability and inelastic collapse. The molecular dynamics results support analytical results, but also exhibit a novel fluctuational breakdown of mean-field descriptions.

15 pages, 4 figures