Flow behavior of colloidal suspensions in complex fluids has been the subject of intense research [1][2][3][4][5][6][7][8][9][10]. The widespread use of suspensions in industry has for a long time provided motivation for these investigations. Early theoretical research was based on hydrodynamic models. For example, the low shear rate viscosity of a suspension with low volume fraction, , of solid spherical particles can be estimated by the often-quoted Einstein equation [11] where η m represents the linear viscosity of the matrix. Underlying these theoretical models is the assumption of continuum viscoelasticity; i.e. the colloidal particle is saturated by the dispersing medium and is large enough such that the non-hydrodynamic contributions such as Brownian motion, surface forces, and van der Waals interactions between particles are negligible. As a result, models derived from hydrodynamic theories are inherently independent from the size of filler particles and the nature of interfacial bonding. It is therefore not surprising that these equations appear to explain only the results of micron-sized particles [5,12,13]. They are indeed entirely inadequate to elucidate the results obtained from the reinforcing fillers of colloidal and sub-colloidal size such as silica and graphitized carbon black in non-Newtonian fluids [12,14,15].