We describe a direct simulation method that effectively determines the appropriate hindered settling behavior of polydisperse particle suspensions. Our method allows adequate representation of the hydrodynamic interactions as well as system specific colloidal interactions. Simulation results are in good agreement with experimental data obtained by MRI imaging. Our results demonstrate the importance of hydrodynamic diffusion and Brownian diffusion in the hindered settling suspensions.