A fully resolved direct numerical simulation of flow and heat transfer is presented for slender randomly packed bed reactors. The flow and temperature field are solved over a non-body fitted, non-conformalCartesian computational domain. The coupling between fluid and solid (both spherical particles and cylindrical wall) is enforced by a second order accurate, sharp interface immersed boundary method (IBM). The present numerical technique neither requires any challenging volumetric mesh generation process nor demands manipulation of the geometry near the particle-particle and particle-wall contact points. Conjugate heat transfer has been considered where the temperature field is calculated both inside the solid particles and in the fluid. A discrete element method (DEM) is used to generate the random packings of spherical particles in the cylindrical column, and a methodology is proposed to calculate the radial porosity profile of the bed. The column-to-particle diameter ratio (N) is varied from 4 to 8, and two separate cases have been considered where N -> infinity. The particle Reynolds number (Red) is varied from 1 to 500. The numerically obtained pressure drop and overall wall-to-bed heat transfer coefficient for different simulation cases are critically compared with empirical correlations and a good agreement is reported. Moreover, based on the current numerical results, correlations are proposed for the pressure drop and the wall-to-bed heat transfer coefficient. The effect of the column-to-particle diameter ratio (N) on both the flow and heat transfer, as-well-as the effect of the solid to fluid thermal conductivity ratio on the conjugate heat transfer are discussed. Furthermore, the fully resolved accurate numerical simulations have helped to elucidate the detailed pore-scale flow and heat transfer feature in the packed beds.