Численное моделирование одномерных многофронтовых задач Стефана

Abstract

Heat transfer with the phase transition is traditionally described with the Stefan problem which consists of a system of parabolic differential equations with usual boundary conditions and an extra condition at the phase transition front. It is formally possible to use one equation of heat transfer type but then a delta-function appears in one of the coefficients which corresponds to the latent heat of phase transition emission at the temperature of phase transition. The widely used method of “continuous calculation” involves the substitution of delta-function by a delta-shaped function. It transforms the Stefan problem into a boundary problem for non-linear equation of heat transfer. Yet, in this approach the result of calculations is the temperature field, and identifying the position of the phase transition front is difficult. Meanwhile, the point of calculations is often the estimation of the dynamic of phase transition. That is why there have been developed many methods of the Stefan problem solution to find the front coordinates. A general disadvantage of these methods is the fact that they cannot be used when there are several fronts appearing, disappearing, changing directions, and joining. Following the dynamic of each front is difficult in this case. In this paper we examine the Stefan problem as the problem of moist soil freezing and melting. A method of solving the Stefan problem which identifies the front as a zero isotherm is developed in this article. This method excludes the necessity to follow the evolution of each front. The Stefan problem is considered as the limiting case of the more general problem of phase transition in some temperature range. A number of standard transformations and Green’s function help to present the Stefan problem as an integral equation. Approximate solution is given by a recurrence formula with an example of numerical simulation. The results correspond with the results of the “continuous calculation” method. However, here we have not the temperature field, but the evolution of the front in time. The numeric simulation shows that the method presented here is convenient for multi-front Stefan problems solution. It should be noted that estimators in form of a system of functions, which majorize the required solution alternately above and below (if we need such estimators), may be obtained only for monotonous Stefan problems. For non-monotonous problems, it is not possible. This issue requires additional investigating.

Теплообмен с фазовым переходом традиционно описывается задачей Стефана, которая представляет из себя систему дифференциальных уравнений параболического типа с обычными краевыми условиями и дополнительным условием на границе фазового перехода. Можно формально перейти к одному уравнению типа теплопроводности, но тогда в одном из коэффициентов появляется дельта-функция, отражающая выделение Джоулева тепла при температуре фазового перехода. Широко распространенный метод «сквозного счета», заключающийся в замене дельта-функции на дельтообразную функцию, сводит задачу Стефана к краевой задаче для нелинейного уравнения теплопроводности. Но при таком подходе результатом расчетов является поле температур, по которому трудно идентифицировать положение границы фазового перехода. Между тем расчет часто производится главным образом для установления динамики фазовых превращений. Поэтому разработано большое количество методов решения задачи Стефана, в которых искомой величиной является именно координата фронта. Общим недостатком этих методов является то, что они мало пригодны для ситуаций, когда фронтов несколько, когда эти фронты появляются и исчезают, меняют направление движения, сливаются друг с другом. В этом случае следить за движением каждого фронта становится нелегким делом. В данной статье на примере задачи о промерзании – оттаивании влажного грунта под воздействием сезонных колебаний температуры поверхности изложен метод решения задачи Стефана, позволяющий получать координату фронта как нулевую изотерму. В этом методе исключена необходимость специального контроля за эволюцией каждого фронта. Задача Стефана рассматривается как предельный случай более общей задачи о фазовом переходе в некотором диапазоне температур. Стандартные преобразования и применение функции Грина позволяют записать эту задачу в виде интегрального уравнения. Приближенное решение получается в виде рекуррентной формулы. В статье приводится пример расчета. Полученные результаты достаточно хорошо согласуются с расчетами по методу «сквозного счета». Но здесь мы получаем не поле температур, а зависимость положения границы фронта от времени. Численные эксперименты показывают, что изложенный в настоящей работе метод является удобным способом моделирования многофронтовых задач Стефана. Отметим, что оценки в виде системы функций, поочередно мажорирующих искомое решение сверху и снизу (если таковые нужны), можно получить только для монотонных задач Стефана. Для немонотонных задач это сделать не удается. Здесь требуется дополнительное исследование.