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