In this paper we propose a model and numerical scheme for quantifying the motion of the phase transition front and the temperature fields in layered materials using boundary conditions of the third kind at the surface.