One dimensional numerical modeling was carried out for the melting behavior of dry snow and the unsaturated flow when heat was supplied from the bottom surface. Discrepancy between the previous experimental data and the present numerical results is substantially reduced by considering the density change of water permeation layer due to the infiltration of meltwater. In the parametric study for effective thermal conductivity, it was found that the effect of this parameter to the behavior of snow melting is minor. Sensitivity analysis showed that the melting time is most sensitive to changes in supplied heat flux, snow temperature, and bulk density, whereas snow bulk density and residual saturation have a significant effect on the height of water permeation layer in snow.