This paper deals with the numerical simulation of a turbulent flow around two-dimensional bodies by the finite volume method with non-orthogonal body-fitted grid. The governing equations are expressed in Cartesian velocity components and solution is carried out using the SIMPLE algorithm for collocated arrangement of scalar and vector variables. Turbulence is modeled by the k- ε turbulence model and wall functions are used to bridge the solution variables at the near wall cells and the corresponding quantities on the wall. A simplified pressure correction equation is derived and proper under-relaxation factors are used so that computational cost is reduced without adversely affecting the convergence rate. The numerical procedure is validated by comparing the computed pressure distribution on the surface of NACA 0012 and NACA 4412 hydrofoils for different angles of attack with experimental data. The grid dependency of the solution is studied by varying the number of cells of the C-type structured mesh. The computed lift coefficients of NACA 4412 hydrofoil at different angles of attack are also compared with experimental results to further substantiate the validity of the proposed methodology.