We study suspensions of deformable (viscoelastic) spheres in a Newtonian
solvent in plane Couette geometry, by means of direct numerical simulations. We
find that in the limit of vanishing inertia the effective viscosity μ of
the suspension increases as the volume-fraction occupied by the spheres Φ
increases and decreases as the elastic modulus of the spheres G decreases;
the function μ(Φ,G) collapses to an universal function, μ(Φe),
with a reduced effective volume fraction Φe(Φ,G). Remarkably, the
function μ(Φe) is the well-known Eilers fit that describes the rheology
of suspension of rigid spheres at all Φ. Our results suggest new ways to
interpret macro-rheology of blood