We present a simulation technique for the direct simulation of the motion of spherical particles suspended in fluids at low Reynods numbers. We use this method to simulate bidisperse suspensions at varying size ratios of the suspended spheres and compare to analytic results of Batchelor for the low concentration regime. We find that the pair correlation functions indicate that three body hydrod...