The present research proposes a pressure based approach along with Langmuir slip condition for predicting microscale fluid flows. Using this method, gaseous slip flows in 2 -dimensional microchannels are numerically investigated. Compared to the DSMC simulation, statistical errors could be avoided and computing time is much less than that of the aforementioned molecular approach. Maxwell slip boundary condition is also studied in this research. These two slip conditions give similar results except for the pressure nonlinearity at high Knudsen number regime. However, Langmuir slip condition seems to be more promising because this does not need to calculate the streamwise velocity gradient accurately and to calibrate the empirical accommodation coefficient. The simulation results show that the proposed method using Langmuir slip condition is an effective tool for predicting compressibility and rarefaction in microscale slip flows.