A numerical investigation of the constrained melting of phase change materials within spherical-like capsule is presented. A single-domain enthalpy formulation is used for simulation of the phase change phenomenon. The solution methodology is verified with the melting process inside an isothermal spherical capsule. Especially, the effect of capsule shape on the heat storage is emphasized. Two shape parameters are considered from the real capsule shape showing good characteristics of heat storage and the effect of these parameters is examined. Early during the melting process, the conduction mode of heat transfer is dominant. Thus the capsule shape with large surface area is desirable. However, the capsule shape with large surface area plays negative role on the strength of buoyancy-driven convection that becomes more important as melting continues.