A simulation method to analyze mesoscopic rheology of blood by multi-scale modeling of blood flow was proposed. Blood flow at a microscopic scale was modeled as flowing red blood cells, while blood flow at a macroscopic scale was treated as a continuum where a local variation of viscosity which is a function of a local concentration of red blood cells was taken into account. The simulation results were in good agreement with in vivo results. The paper also discusses the influence of hematocrit on the velocity profile of a macroscale flow.