This paper is concerned with the numerical investigation of a macroscopic model for complex fluids in “1+2” dimension case. We consider the planar pressure driven flow where the direction of the molecules is constrained in the shear plane. The modified Crank-Nicolson finite difference scheme satisfying a discrete energy law will be developed. By using this scheme, it is observed numerically that the direction of the molecules will tumble from the boundary layer and later on the inner layer with a much longer time period. This is consistent with the theoretical prediction. Moreover, we find some complex phenomena, where the tumbling rises from boundary layer and is then embedded into the interior area more clearly when the viscosity coefficient µ of the macro flow has a larger value. The norm of the molecular director d will endure greater change as well. This implies that the viscosity of flow plays the role of an accelerator in the whole complex fluids. Comparing these results with the theoretical analysis, we can find that the gradient of the velocity has direct impact on the tumbling phenomena. These results show that the proposed scheme is capable of exploring some physical phenomena embedded in the macro-micro model.