Methods to simulate colloid suspension structure, dynamics, and rheology are presented and distinguished by whether the solvent is modeled by a continuum or by meso-scale particles, and whether or not a mesh is used. The unmeshed/continuum solvent methods include various forms of Stokesian Dynamics, while the meshed continuum solvent methods include the Arbitrary Lagrangian–Eulerian method and the Immersed Boundary Method. Mesoscopic particles are used on a mesh in the Lattice–Boltzmann (LB) method and are mesh-free in Dissipative Particle Dynamics (DPD). These and other methods are reviewed and guidance provided on which ones are preferred in various situations. For spherical hard particles at high concentrations, without inertia, when accurate treatment of hydrodynamic interactions is needed, versions of Stokesian dynamics set benchmarks for accuracy and efficiency. Other methods become preferable depending on the speed needed for computation, the numbers of particles to be simulated, how well hydrodynamic interactions need to be modeled, whether Brownian motion or particle inertia is important, the shapes and deformability of the particles, the complexity of the flow domain, and whether the solvent is non-Newtonian. Information in this chapter indicates how well the different methods perform in their present state of development.