The slip velocity of the flow at the porous-fluid interface is studied using a 2D porous region, modeled as an array of circular cylinders in square and hexagonal arrangements to resemble the structure of real porous media. The Navier-Stokes and continuity equations are solved using a hybrid FEM-FVM solver by applying appropriate macroscopic boundary conditions. The existence of a slip velocity at the porous-fluid interface is demonstrated and it is observed that the slip velocity increases with increase in porosity and flow Reynolds number.