Abstract In this study, two phase simulation of nanofluid flow and heat transfer between parallel plates is investigated. The important effects of Brownian motion and thermophoresis have been included in the model of nanofluid. The governing equations are solved via homotopy perturbation method. According to comparison with previous works, this method has good accuracy to solve this problem. The semi analytical investigation is carried out for different governing parameters namely; the squeeze number, Hartmann number, Schmidt number and Eckert number. The results indicate that absolute skin friction coefficient decreases with increase of Hartmann number and squeeze number. Also it can be found that that Nusselt number is an increasing function of Hartmann number, Eckert number and Schmidt number but it is a decreasing function of squeeze number.