Abstract
A numerical study is made of the melting of ice in a rectangular porous cavity heated from above. The Landau transformation is used to immobilize the ice-water interface, and the Darcy-Boussinesq equations are solved by a finite-difference technique. Results are analyzed in terms of the heating temperature and the aspect ratio of the cavity. A comparison is made with the case of melting from below. It was found that melting from above is more effective than melting from below when the heating temperature is between 0 and 8°C: convection arises earlier, the melting process is faster, and the total melt at steady state is thicker. The critical time for onset of convection is minimum when the upper boundary is heated at 6°C. At this heating temperature, one also obtains a maximum heat transfer rate (Nusselt number).