The research on thermal energy storage (TES) systems have received a lot of attention in recent decades for sustainable use of thermal energy in various industrial and residential applications. The existing challenge in designing the TES is the response time of charging and discharging cycles that keeps these systems away from wide utilization in industries. Literature data show that beside the low thermal conductivity of most phase change materials (PCMs) as active media in TES systems, the poor flow distribution may be another factor affecting the response rate. This study aims to considerably reduce the response time by packing the PCMs in a bed of spheres made of high thermal conductivity material. The response rate during the charging cycle is studied numerically by passing hot water at 70 °C over the packed bed of spheres. The numerical analysis is performed using ANSYS Fluent 19. The PCM used in this study is a paraffin and has a melting point of 48 °C. The response rate of the system is studied and it is compared to other similar systems mentioned in literature. The amount of energy storage is also studied by changing the flow rate of water.