In this study, three randomly packed beds with varying tube-to-particle diameter ratios (D/d) are constructed using the discrete element method (DEM) and simulated via CFD under low pore Reynolds numbers (Rep < 100). An innovation of this research lies in the application of hydrogen in randomly packed beds, coupled with the consideration of its temperature-dependent thermal properties. The axial analysis of the heat transfer characteristics shows that PB−5 and PB−6 achieve thermal equilibrium 44% and 58% faster than PB−4, respectively, demonstrating enhanced heat transfer efficiency. However, at higher flow rates (0.8 m/s), the large-sized fluid channels in PB−6 severely impact the heat transfer efficiency, slightly reducing it compared to PB−5. Additionally, this study introduces a localized segmentation method for calculating the axial local Nusselt number, showing that the axial local Nusselt number (Nu) not only exhibits an inverse relationship with the axial porosity distribution, but also matches its amplitude fluctuations. The wall effect significantly impacts the flow and temperature distribution in the packed bed, causing notable velocity and temperature oscillations in the near-wall region. In the near-wall region, the average temperature is lower than in the core region, and the axial temperature profile exhibits more intense oscillations. These findings may provide insights into the use of hydrogen in randomly packed beds, which are vital for enhancing industrial applications such as hydrogen storage and utilization.