Melting in porous media within a rectangular enclosure with presence of natural convection is simulated using an interfacial tracking method. This method combines the advantages of the both deforming and fixed grids methods. Convection in the liquid region is modeled using the Navier-Stokes equation with Darcy's term and Forchheimer's extension. The results obtained by using the interfacial tracking method validated by comparing with the existing experimental and numerical results. The results show that the interfacial tracing method is capable of solving natural convection controlled melting problems in porous media at both high and low Prandtl numbers.