Приведены математическая модель и численное решение задачи нестационарной теплопроводности многослойных сред с учетом фазовых превращений на границе раздела фаз. При этом учитывается перемещение фронта фазовых превращений по сечению многослойных сред. Разработанная математическая модель применяется для расчета полей температур изо-тропных и анизотропных сред при различных граничных условиях.