Viscometric effects in a dense hard-sphere fluid

A recently proposed Monte Carlo simulation method for the Enskog equation is applied to uniform shear flow. This state is characterized by uniform density n and temperature T and a linear velocity profile: u(r)=a•r, aαβ=aδαxδβy. In each unit time step Δt, the peculiar velocities {Vi} of N particles are updated in two stages. In the free streaming stage, the velocity Vi is changed into Vi→Vi−a•Vi Δt; in the collision stage, with probability equal to , where is a random unit vector, σ is the diameter of the spheres, being the velocity of a random partner j, and χ(n) is the equilibrium pair correlation function at contact. The kinetic and collisional transfer contributions to the pressure tensor are calculated as functions of density. The Navier-Stokes shear viscosity is seen to agree with the theoretical value. Furthermore, the nonlinear Burnett coefficients associated with normal stresses are obtained.