Laminar separation bubble and transitional flow over the NACA0012 are investigated at a moderate range of Reynolds numbers. A Reynolds-Averaged Navier-Stokes code is coupled with an empirical transition model that can predict transition onset points and the length of transition region. Without solving the boundary layer equations, approximated e-N method is directly applied to the RANS code and iteratively solved together. The computational results are compared with the experimental data for the NACA0012 airfoil. Results of transition onset point and the length are compared well with experimental data and Xfoil prediction. The present RANS results show at high angles of attack better agreement with experimental data than Xfoil results using the boundary layer equations.