In this paper, a mathematical model of shell-and-tube latent heat thermal energy storage (LUTES) unit of two-dimension of three phase change materials (PCMs) named PCM1, PCM2 and PCM3 with different high melting temperatures (983 K, 823 K and 670 K, respectively) and heat transfer fluid (HTF: air) with flowing resistance and viscous dissipation based on the enthalpy method has been developed. Instantaneous solid-liquid interface positions and liquid fractions of PCMs as well as the effects of inlet temperatures of the air and lengths of the shell-and-tube LHTES unit on melting times of PCMs were numerically analyzed. The results show that melting rates of PCM3 are the fastest and that of PCM1 are the slowest both x, r directions. It is also found that the melting times of PCM1, PCM2 and PCM3 decrease with increase in inlet temperatures of the air. Moreover, with increase in inlet temperatures of the air, decreasing degree of their melting times are different, decreasing degree of the melting time of PCM1 is the biggest and that of PCM3 is the smallest. Considering actual application of solar thermal power, we suggest that the optimum lengths are L-1 = 250 mm, L-2 = 400 mm, L-3 = 550 mm (L = 1200 mm) which corresponds to the same melting times of PCM1, PCM2 and PCM3 are about 3230 s and inlet temperature of the air is about 1200 K. The present analysis provides theoretical guidance for designing optimization of the shell-and-tube LHTES unit with three PCMs for solar thermal power. (C) 2013 Elsevier Ltd. All rights reserved.