In this study, the utilisation of composite phase change materials (PCMs) slabs with fins in a thermal energy storage (TES) system is numerically analysed. The storage system consists of three PCM layers at different melting temperatures. First, the predicted results are verified against the analytical solution that is taken from literature. Afterwards, the influence of the design parameters such as fin length and fin spacing on time for complete melting of the multi-layer TES are examined. Custom field functions are integrated into the ANSYS-FLUENT software to assess the local and volume averaged entropy generation within the system. The results show that at constant fin spacing, that increasing fin length from 20 mm to 80 mm enhances the melting rate up to 50%. In addition, when the fin spacing is doubled, the elapsed time for complete melting decreased about 37%. Moreover, the fin arrangement has a considerable effect on entropy generation.