Thermal energy storage units that utilize phase change materials have been widely employed to balance temporary temperature alternations and store energy in many engineering systems. In the present paper, an operational approach is proposed to the Tau method with standard polynomial bases to simulate the phase change problems in latent heat thermal storage systems, that is, the two-dimensional solidification process in rectangular finned storage with a constant end-wall temperature. In order to illustrate the efficiency and accuracy of the present method, the solid-liquid interface location and the temperature distribution of the fin for three test cases with different geometries are obtained and compared to simplified analytical results in the published literature. The results indicate that using a two-dimensional numerical approach can predict the solid-liquid interface location more accurately than the simplified analytical model in all cases, especially at the corners.