Numerical analysis of a one-dimensional model of a melting-freezing snowpack

The article is devoted to a numerical study of a one-dimensional non-stationary problem on thermomechanical processes in a snowpack with account of effects of melting and freezing. Snow is modeled as a continuous medium consisting of water, air and porous ice skeleton. The governing equations of snow are based on the fundamental conservation laws of continuum mechanics. A finite-difference algorithm is constructed and a series of numerical experiments is fulfilled. The results of the computations correspond well to laboratory observations.