This study aims to numerically investigate a two dimensional and steady heat transfer over a cylinder in a porous medium with suspending nanoparticles. Buongiorno model is adopted for nanofluid transport on a free convection flow taking the slip mechanism of Brownian motion and thermophoresis into account. The Boussinesq approximation is considered to account for buoyancy. The boundary layer conservation equations are transformed into dimensionless and then elucidated using a 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 the Forchheimer parameter, Λ, clearly swamps the nanofluid momentum development, decreases the flow for some distance near the cylinder viscous region, later it reverses the trend, and asymptotically reaches the far field flow velocity. Furthermore, as thermophoresis parameter increases, the heat transfer and nanoparticle volume concentration increase within the boundary layer. The present results are validated with the available results of a similar study and is found to be in good coincidence. The study finds applications in heat exchangers technology, materials processing, and geothermal energy storage etc.