Выпуск:
2017. Том 3. №3Об авторах:
Аксенов Борис Гаврилович, доктор физико-математических наук, профессор кафедры промышленной теплоэнергетики, Тюменский индустриальный университет; aksenovbg@tyuiu.ruАннотация:
Теплообмен с фазовым переходом традиционно описывается задачей Стефана, которая представляет из себя систему дифференциальных уравнений параболического типа с обычными краевыми условиями и дополнительным условием на границе фазового перехода. Можно формально перейти к одному уравнению типа теплопроводности, но тогда в одном из коэффициентов появляется дельта-функция, отражающая выделение Джоулева тепла при температуре фазового перехода. Широко распространенный метод «сквозного счета», заключающийся в замене дельта-функции на дельтообразную функцию, сводит задачу Стефана к краевой задаче для нелинейного уравнения теплопроводности. Но при таком подходе результатом расчетов является поле температур, по которому трудно идентифицировать положение границы фазового перехода. Между тем расчет часто производится главным образом для установления динамики фазовых превращений. Поэтому разработано большое количество методов решения задачи Стефана, в которых искомой величиной является именно координата фронта. Общим недостатком этих методов является то, что они мало пригодны для ситуаций, когда фронтов несколько, когда эти фронты появляются и исчезают, меняют направление движения, сливаются друг с другом. В этом случае следить за движением каждого фронта становится нелегким делом. В данной статье на примере задачи о промерзании — оттаивании влажного грунта под воздействием сезонных колебаний температуры поверхности изложен метод решения задачи Стефана, позволяющий получать координату фронта как нулевую изотерму. В этом методе исключена необходимость специального контроля за эволюцией каждого фронта. Задача Стефана рассматривается как предельный случай более общей задачи о фазовом переходе в некотором диапазоне температур. Стандартные преобразования и применение функции Грина позволяют записать эту задачу в виде интегрального уравнения. Приближенное решение получается в виде рекуррентной формулы. В статье приводится пример расчета. Полученные результаты достаточно хорошо согласуются с расчетами по методу «сквозного счета». Но здесь мы получаем не поле температур, а зависимость положения границы фронта от времени. Численные эксперименты показывают, что изложенный в настоящей работе метод является удобным способом моделирования многофронтовых задач Стефана. Отметим, что оценки в виде системы функций, поочередно мажорирующих искомое решение сверху и снизу (если таковые нужны), можно получить только для монотонных задач Стефана. Для немонотонных задач это сделать не удается. Здесь требуется дополнительное исследование.
Ключевые слова:
Список литературы:
