To investigate the microstructural evolution dependency on the apparent viscosity in shear-thickening fluids (STFs), a hybrid mesoscale model combined with stochastic rotation dynamics (SRD) and molecular dynamics (MD) is used. Muller-Plathe reverse perturbation method is adopted to analyze the viscosities of STFs in a two-dimensional model. The characteristic of microstructural evolution of the colloidal suspensions under different shear rate is studied. The effect of diameter of colloidal particles and the phase volume fraction on the shear thickening behavior is investigated. Under low shear rate, the two-atom structure is formed, because of the strong particle attractions in adjacent layers. At higher shear rate, the synergetic pair structure extends to layered structure along flow direction because of the increasing hydrodynamics action. As the shear rate rises continuously, the layered structure rotates and collides with other particles, then turned to be individual particles under extension or curve string structure under compression. Finally, at the highest shear rate, the strings curve more severely and get into two-dimensional cluster. The apparent viscosity of the system changes from shear-thinning behavior to the shear-thickening behavior. This work presents valuable information for further understanding the shear thickening mechanism. (C) 2017 Elsevier Inc. All rights reserved.