Численные методы решения задачи Стефана

Abstract

All of the most well-known numerical methods for solving the Stefan’s problem, as well as a new method developed by the author are considered for the purpose of choosing the most efficient of them from a perspective of accuracy and computational speed. A comparison is carried out on the results of solving the problem for the boundary motion of “ice-water” phase transition around the vertical well passing through the thickness of permafrost. The conclusions, which are distributed to other multidimensional and multi-front statements of the Stefan’s problem, are made. The mathematical model, the brief description of the considered numerical methods and the boundaries of their applicability are presented. The comparison shows the advantages and disadvantages of different methods. It is demonstrated that the use of the explicit scheme leads to a marked increase in computation time, the six-point symmetric scheme may have oscillating solution; therefore, the implicit scheme is the most preferred. It is concluded that the most efficient method for one-dimensional and one front Stefan’s problems is the method of catching the front in the grid node using the implicit scheme, and the most efficient method for multi-dimensional and multi-front Stefan’s problems is the enthalpy method using the implicit scheme, which has been developed by the author.

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