The two-phase Stefan problems with phase formation and depletion are special cases of moving boundary problems with interest in science and industry. In this work, we study a solidification problem, introducing a front-fixing transformation. The resulting nonlinear partial differential system involves singularities, both at the beginning of the freezing process and when the depletion is complete, that are treated with special attention in the numerical modelling. The problem is decomposed in three stages, in which implicit and explicit finite difference schemes are used. Numerical analysis reveals qualitative properties of the numerical solution spatial monotonicity of both solid and liquid temperatures and the evolution of the solidification front. Numerical experiments illustrate the behaviour of the temperatures profiles with time, as well as the dynamics of the solidification front.