Melting in porous media within a rectangular enclosure with the presence of natural convection is simulated using an interfacial tracking method. This method combines the advantages of both the deforming and fixed-grid methods. Convection in the liquid region is modeled using the Navier-Stokes equation with Darcy's term and Forchheimer's extension. The results are obtained by using the interfacial tracking method, which is 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.thermal conductivities, k=k ' K = modified dimensionless thermal conductivity K s' = ratio of thermal conductivities, k s =k ' k ' = thermal conductivity in liquid, W=m K Nu = Nusselt number at heated wall, hH=k ' P = dimensionless pressure, pH 2 = 2 l Pr = Prandtl number of liquid phase-change materials, ' = ' p = pressure, Pa Ra = Rayleigh number, Gr Pr S = dimensionless location of solid-liquid interface, s=H Ste = Stefan number, c ' T h T m =h s' s = location of solid-liquid interface, m S 0 = dimensionless location of solid-liquid interface at last time step T = temperature, K t = time, s U = dimensionless velocity in x direction, uH= ' u = velocity component in x direction, m=s U I = dimensionless solid-liquid interfacial velocity, u I H= ' u I = solid-liquid interfacial velocity, m=s, @s=@t V = dimensionless velocity in y direction, vH= ' , or volume, m 3 v = velocity component in y direction, m=s W = width of enclosure, m X = dimensionless coordinate, x=H x = dimensional coordinate, m Y = dimensionless coordinate, y=H y = dimensional coordinate, m = thermal diffusivity, m 2 =s = liquid fraction in Eq. (2) = liquid fraction in Eq. (3) " = porosity = dimensionless temperature, T T m =T h T m = permeability, m 2 = viscosity of the liquid phase-change materials, kg=ms = kinematic viscosity, m 2 =s = density, kg=m 3 = dimensionless time, ' t=H 2 = heat capacity ratio, c=c ' Subscripts c = cold E = east e = east face of control volume eff = effective f = fluid h = heated I = interface i = initial ' = liquid m = melting point N = north n = north face of control volume p = porous matrix ref = reference S = south s = solid W = west w = west face of control volume