Abstract
Equations of motion for model colloidal particles are derived, which include a mean-field local density representation for the many-body-hydrodynamics (MBH) interactions. These are implemented as a modification of Ermak's 1975 Brownian dynamics free-draining algorithm, with a local density prescription for the friction coefficient. The computed long-time self-diffusion coefficients derived from this model agree well with experiment at all volume fractions. The hydrodynamics model relative and dynamic viscosities are larger than those of the free-draining algorithm, reflecting slower structural relaxation processes taking place in the new model.