The dynamics for a system of hard spheres with dissipative collisions
is described at the levels of statistical mechanics, kinetic theory, a
nd simulation. The Liouville operator(s) and associated binary scatter
ing operators are defined as the generators for time evolution in phas
e space. The BBGKY hierarchy for reduced distribution functions is giv
en, and an approximate kinetic equation is obtained that extends the r
evised Enskog theory to dissipative dynamics. A Monte Carlo simulation
method to solve this equation is described, extending the Bird method
to the dense, dissipative hard-sphere system. A practical kinetic mod
el for theoretical analysis of this equation also is proposed. As an i
llustration of these results, the kinetic theory and the Monte Carlo s
imulations are applied to the homogeneous cooling state of rapid granu
lar flow.