An adaptation of lattice-gas cellular automata to the simulation of solid-fluid suspensions is described. The method incorporates both dissipative hydrodynamic forces and thermal fluctuations. At low solid densities, theoretical results for the drag force on a single disk and the viscosity of a suspension of disks are reproduced. The zero–shear-rate viscosity has been obtained over a range of packing fractions and results indicate that simulations of three-dimensional suspensions are feasible.

• Received 14 December 1987

DOI:https://doi.org/10.1103/PhysRevLett.60.975