A comprehensive model for effective viscosity of nanofluids has been proposed. Particle clustering and interfacial layer formation have been taken into account. PSD analysis has been utilized to characterize the clusters. Effects of clusters, surfactants and particle diameter have been showcased. Model is found to be in good agreement with the experimental data.