La viscosità delle sospensioni non Browniane di sfere rigide immerse in una matrice Newtoniana sottoposte a flusso di shear semplice nel range delle concentrazioni dilute e semidilute è stata prevista numericamente modellando la sospensione con la Dinamica Stokesiana. In particolare, la viscosità della sospensione può essere sviluppata in serie di Taylor in funzione della concentrazione delle sfere e il coefficiente di ordine due è stato determinato con l’analisi numerica. Le previsioni numeriche correttamente degenerano nel limite lineare predetto teoricamente da Einstein, ed al secondo ordine sono in accordo qualitativo con i dati sperimentali di Pasquino et al..