Numerical modeling was performed to simulate the melting process of a fixed volume/mass phase-change material (PCM) in different shell-and-tube type latent thermal energy storage units with identical heat transfer area. The effect of liquid PCM natural convection (NC) on the latent heat storage performance of the pipe and cylinder models was investigated using a 3D numerical model with FLUENT software. Result shows that NC can cause a non-uniform distribution of the solid-liquid interface, which accelerates PCM melting rate. The PCM melting rate and heat storage rate in the horizontal cylinder model are higher than those in the horizontal pipe model because of the combined effects of heat conduction and NC. A comparative study was conducted to determine the effects of horizontal and vertical shell-and-tube models with different heat transfer fluid (HTF) inlets including the effects of NC. The results indicate that the vertical model with an HTF inlet at the bottom exhibits the highest PCM melting rate and heat storage rate for the pipe models. For the cylinder models, the horizontal model and the vertical model with an HTF inlet at the bottom can achieve nearly the same completed melting time. In addition, NC has minimal effect on any model with an HTF inlet at the top.