The present study is a 2-D numerical study which discusses the thermohydraulic performance of packed bed duct under local thermal non-equilibrium and steady state conditions exposed to forced convection boundaries (BC-6). The numerical model is well validated with the experimental work reported in the literature and found to be accurate enough to perform further parametric investigations. The analysis is done to see the effect of different ball diameter (5 mm, 9 mm and 11 mm), bed porosity, heat transfer fluid (air, water and engine oil - Pr = 0.70–645) and ball material (EPS, steel and bronze) on heat transfer and fluid flow characteristics in turbulent flow region of Reynolds number ranging from 900 to 14,320. The numerical results obtained using commercial CFD software COMSOL shows that, the heat transfer coefficient and pressure drop in packed bed increases with increase in ball diameter, thermal conductivity of ball and bed porosity which is exactly reverse as reported in literature for BC-1 to BC-5. The maximum thermal performance factor with bronze particles is 2.45 and 24.5 times more than stainless steel and EPS particles respectively for larger particle size and bed porosity. Engine oil exhibits significantly higher heat transfer coefficient (9.7 × 105 W/m2K) and pressure drop (3.3 × 106 Pa) compared to water (40,126 W/m2K and 2028 Pa) and air (600 W/m2K and 250 Pa) respectively. Overall, the combination of water as heat transfer fluid along with bronze particles of larger diameter and larger bed porosity emerges as the optimal choice for enhancing heat transfer in the packed duct exposed to forced convection boundary condition BC-6.