The study aims to investigate numerically a two dimensional, steady, heat transfer over a cylinder in porous medium with suspending nanoparticles. Buongiorno model is adopted for nanofluid transport on free convection flow taking the slip mechanism of Brownian motion and thermophoresis into account. Boussinesq approximation is considered to account for buoyancy. The boundary layer conservation equations are transformed into dimensionless, and then elucidated using robust Keller-box implicit code numerically. The numerical results are displayed graphically and deliberated quantitatively for various values of thermo-physical parameters. Our results shows that, increasing Forchheimer parameter, Λ, clearly swamps the nanofluid momentum development, decreasing the flow for some distance near the cylinder viscous region, later its reverse the trend and asymptotically reaches the far field flow velocity. Furthermore, as increases thermophoresis, heat transfer and nanoparticle volume concentration increased in the boundary layer. The present results are validated with the available results of similar study and is found to be in good coincident. The study finds applications in heat exchangers technology, materials processing and geothermal energy storage etc.