Архив статей журнала
Цель – установить кинетические закономерности расплавления тепловыделяющего цилиндрического элемента в заведомо надкритических условиях с помощью численного моделирования. Объектом исследования является процесс плавления в однородном образце, выделяющем теплоту за счет протекания реакции или электромагнитного нагрева. Теплофизические свойства образца принимаются постоянными в пределах твердой и жидкой фаз. Основным инструментом исследования является численная модель, построенная на основе нестационарной задачи Стефана в тепловыделяющем теле и включающая описание процессов теплопроводности и плавления. Фазовый переход описывается в энтальпийном представлении. Для выбора параметров численной модели (шагов сетки) проводится исследование точности разностной схемы. В результате проведенных исследований получены расчетные зависимости основных характеристик плавления (время расплавления и максимальная температура образца в момент расплавления) от управляющих параметров (интенсивность тепловыделения, величина теплового эффекта плавления, отношение коэффициентов теплопроводности фаз). С помощью некоторых приближений (усреднение температуры, квазистационарное распределение) получены формулы для оценки времени расплавления исследуемого образца. Расчеты показали, что изменение теплофизических свойств образца (коэффициентов теплопроводности, теплового эффекта) оказывает существенное влияние на скорость его плавления. Установлено, что зависимость времени расплавления от интенсивности тепловыделения и теплового эффекта фазового перехода качественно совпадает с приближенными моделями, но существенно отличается от них количественно, особенно в области малых отклонений от критической интенсивности тепловыделения. Проведенные расчеты могут быть использованы при оценке термомеханической устойчивости материалов с внутренним тепловыделением. Разработанная численная модель дает возможность исследовать процессы плавления в широком диапазоне условий, в том числе при изменении граничных условий.