On a numerical method for solving one-dimensional Stefan-type problems

A finite-difference method for solving one-dimensional one- and two-phase Stefan-type problems with multiple phase boundaries is proposed. The essence of the method is that the problem under study is solved step by step when moving from one time layer to another. The time step value, the spatial grid of the time layer, the corresponding nodal values of temperature, and the position of interphase boundaries are simultaneously determined at each step.