Utilizing phase change materials in thermal energy storage systems is commonly considered as an alternative solution for the effective use of energy. This study presents numerical simulations of the charging process for a multitube latent heat thermal energy storage system. A thermal energy storage model, consisting of five tubes of heat transfer fluids, was investigated using Rubitherm phase change material (RT35) as the. The locations of the tubes were optimized by applying the Taguchi method. The thermal behavior of the unit was evaluated by considering the liquid fraction graphs, streamlines, and isotherm contours. The numerical model was first verified compared with existed experimental data from the literature. The outcomes revealed that based on the Taguchi method, the first row of the heat transfer fluid tubes should be located at the lowest possible area while the other tubes should be spread consistently in the enclosure. The charging rate changed by 76% when varying the locations of the tubes in the enclosure to the optimum point. The development of streamlines and free-convection flow circulation was found to impact the system design significantly. The Taguchi method could efficiently assign the optimum design of the system with few simulations. Accordingly, this approach gives the impression of the future design of energy storage systems.