This study explores the impact of tube rotation on the melting performance of a multi-tube latent heat energy storage system through numerical analysis. The system involves a phase change material (PCM) contained within a tube, which absorbs heat from hot water flowing through a tube surrounding the PCM, as well as two tubes enveloped by the PCM. The investigation considers the angle between the line connecting the centers of the two inner tubes and the horizon (α = 0◦, 45◦, 90◦, and 135◦), along with the direction of rotation of inner tubes. Results are compared with data from scenarios involving only outer tube rotation and stationary systems. Temporal variations in temperature and liquid fraction contours, along with the temporal evolution of average temperature and average liquid fraction of PCM are presented. It was determined that for each α, the melting performance varied from best to worst across the system configurations as follows: the system with all three rotating tubes, the system with only the outer rotating tube, and the stationary system. The minimum time required to complete PCM melting for cases at α = 0◦, 45◦, 90◦, and 135◦ is respectively 44.92 %, 42.11 %, 41.32 %, and 40.62 % less than that of the stationary system.