A latent heat storage system for concentrated solar plants (CSP) is numerically examined by means of CFD simulations. This study aims at identifying the convective flows produced within the melted phase by temperature gradients and gravity. Simulations were carried out on experimental devices for applications to high temperature concentrated solar power plants. A shell-and-tube geometry composed by a vertical cylindrical tank, filled by a Phase Change Material (PCM) and an inner steel tube, in which the heat transfer fluid (HTF) flows, from the top to the bottom, is considered. The conjugate heat transfer process is examined by solving the unsteady Navier Stokes equations for HTF and PCM and conduction for the tube. In order to take into account the buoyancy effects in the PCM tank the Boussinesq approximation is adopted. The results show that the enhanced heat flux, due to natural convective flow, reduce of about 30% the time needed to charge the heat storage. A detailed description of the convective motion in the melted phase and the heat flux distribution between the HTF and PCM are reported. The effect of the mushy zone constant is also investigated. (C) 2015 Elsevier Ltd. All rights reserved.