Полный текст

В современной металлургии титана широкое распространение получил способ магниетермического восстановления губчатого титана из его тетрахлорида. Способ основан на химической реакции взаимодействия тетрахлорида титана (ТХТ) с расплавленным магнием, в результате которого ТХТ восстанавливается до чистого титана с образованием побочного продукта – хлорида магния. В нашей стране имеются два предприятия – производителя губчатого титана – это АВИСМА – филиал корпорации ВСМПО-АВИСМА (г. Березники) и Соликамский магниевый завод (СМЗ). Подробно технология получения титана изложена в специальной литературе [1, 2]. Здесь ограничимся упоминанием о том, что в современной отечественной практике используются аппараты с цикловым съемом 4,8 и 7,0 т губчатого титана, которые имеют следующие геометрические размеры: внутренние диаметры 1,55 и 1,70 м и длины 3,82 и 4,30 м соответственно. Одним из основных параметров, характеризующих протекание процесса восстановления, является температура в реакционной зоне и по всей высоте аппарата. Вместе с тем значительные размеры основного оборудования, агрессивная среда и высокие температуры во внутреннем объеме аппарата не позволяют организовать достаточно полный контроль над данным параметром. Общая навеска магния, загружаемая в реактор для проведения процессов восстановления, составляет для аппарата АВИСМА 8 т, для аппарата СМЗ – 12 т. При таких количествах жидкого металла влияние конвекции на технологию ведения процесса восстановления весьма существенно. На практике оценить и спрогнозировать этот фактор возможно лишь приблизительно, по косвенным признакам. В предыдущей работе [3], посвященной вопросам разработки математической модели, описана теплофизика процесса восстановления в технологии получения губчатого титана, приводятся постановка задачи, методика решения, а также результаты вычислительного эксперимента. Сравниваются кривые показаний термопар, снятые с реального промышленного процесса, и значения температур, полученные расчетным путем при выполнении вычислительного эксперимента. Показана сходимость реальных и расчетных показаний. Также некоторые результаты тестирования и проверки адекватности разработанной математической модели приведены в работе [4]. Эти результаты были получены в ходе начального этапа работы, когда задача решалась в безконвективном приближении, без учета конвективной составляющей. В данном приближении разработанная модель дает представление о распределении температур в неподвижном расплаве магния в отсутствие естественного конвективного перемешивания. Однако значительный объем жидкой фазы в аппарате и наличие температурных градиентов не оставляют сомнений в том, что влияние конвекции на скорость разогрева или ход самого процесса восстановления при подаче ТХТ в аппарат будет существенным. Возникает задача разработки математической модели, адекватно описывающей неравновесные теплофизические процессы в реакторе восстановления с учетом химизма протекающих реакций и гидродинамики расплавленного металла. Постановка задачи тепловой конвекции. Аппарат восстановления, установленный в шахтную электропечь, рассматривается в цилиндрической осесимметричной системе координат. Общий вид аппарата приведен на рис. 1, а, на рис. 1, б представлена расчетная схема реального объекта. Расчетная область ограничена слева осью симметрии, справа – стенкой реторты, в нижней части – сварным швом крепления сферического днища к обечайке реторты, сверху – заглублением крышки в реторту. Формулировка задачи тепловой конвекции сводится к системе трех дифференциальных уравнений: переноса тепловой энергии, завихренности и Пуассона, которые в обозначениях работ [5, 6] имеют вид (1) где а – коэффициент температуропроводности, а = λ0/(ρс), м2/с, ρ – плотность, кг/м3, с – теплоемкость, Дж/(кг∙К); λэф – эффективный коэффициент теплопроводности гетерогенной системы, Вт/(м∙К); λ0 – масштабное значение коэффициента теплопроводности, Вт/(м∙К); t – температура, К; r, z – соответственно радиальная и осевая координаты, м; q1 – сток тепла парообразования, Вт/м3; q2 – источник тепла от химических реакций, Вт/м3; u, ν – компоненты скорости в проекциях на оси r и z соответственно, м/с; β – коэффициент объемного расширения, К–1, g – ускорение свободного падения, g = 9,8 м/с2. Для решения системы уравнений (1) заданы следующие начальные и граничные условия: t (τ = 0) = t0, ω(τ = 0) = ψ(τ = 0) = 0, где t0 – температура в реторте в начальный момент времени. Рис. 1. Аппарат восстановления: а – общий вид аппарата и печи восстановления: 1 – крышка; 2 – реторта; 3 – шахтная электропечь; 4 – верхний нагреватель; 5 – теплообменник (область охлаждения воздухом реторты в зоне реакции); 6 – нижний нагреватель; 7 – каналы печи для подачи и отвода охлаждающего воздуха; 8 – уровень расплава магния в реторте; б – расчетная схема: Н – расстояние от 0 до крышки аппарата; Нж – расстояние от 0 до поверхности расплава магния; Н–Нж – высота заполнения аппарата газообразным аргоном; ВНВ, ВНН, Вт – высоты верхнего, нижнего нагревателей и теплообменника печи восстановления соответственно; ННВ – расстояние до верхнего нагревателя; Нт – расстояние до теплообменника; R – внутренний радиус реторты Условие теплообмена на правой границе: – в зоне верхнего нагревателя, r = R, ННВ < z < ННВ + ВНВ: ; (2) – в зоне обдува реторты воздухом (теплообменник), r = R, Нт < z < Нт + Вт: ; (3) – в зоне нижнего нагревателя, r = R, 0 < z < ВНН: (4) где , , – соответственно коэффициенты теплоотдачи на наружной поверхности стенки реторты в зонах нагрева верхнего нагревателя, действия теплообменника и нижнего нагревателя печи, Вт/(м2∙К); t, tп – температуры расплава в аппарате и на внутренней поверхности стенки реторты; – плотности тепловых потоков верхнего и нижнего нагревателей соответственно, Вт/м2. Условие теплообмена на верхней границе, z = Н: (5) где tв – температура на наружной поверхности крышки; αкр – коэффициент теплоотдачи от наружной поверхности крышки аппарата в окружающую среду, Вт/(м2∙К); qр.кр – плотность теплового потока от расплава к крышке, Вт/м2. Парогазовую смесь в аппарате условно считаем прозрачной. Условие теплообмена на нижней границе, z = 0: (6) где tдн – температура на наружной поверхности днища; αдн – коэффициент теплоотдачи наружной поверхности днища аппарата, Вт/(м2∙К); qдн – плотность теплового потока от днища аппарата к подине печи, Вт/м2. На оси симметрии – адиабатная граница. Для функции скорости на всех границах в зоне расположения расплава магния принято условие полного прилипания вязкой среды. Методика определения коэффициентов теплоотдачи, входящих в граничные условия (2)–(6), составлена на основании источников [5–7]. Методика позволяет учесть источники и стоки тепла при нагреве аппарата электронагревателями печи, охлаждении зоны реакции аппарата, тепловой эффект от химического взаимодействия ТХТ с магнием в зависимости от расхода тетрахлорида, тепловые потери через крышку аппарата в окружающую среду. Реализация метода решения и результаты вычислительного эксперимента. Разработана математическая модель и компьютерная программа, позволяющие анализировать неравновесные теплофизические процессы в реакторе аппарата для получения губчатого титана. На первом этапе работы при помощи численных методов решено дифференциальное уравнение переноса энергии в бесконвективном приближении, получено физически правдоподобное решение при моделировании различных периодов ведения процесса восстановления. Уравнения переноса тепловой энергии и завихренности решались методом продольно-поперечной прогонки, уравнение Пуассона – методом последовательной линейной верхней релаксации. Основы и реализация указанных методов достаточно детально рассмотрены в специальной литературе [5, 6]. Разработанная математическая модель, учитывающая гидродинамические явления в расплаве магния, тестировалась в соответствии с методикой, указанной в диссертации Ю.А. Селянинова [8], где приводится случай свободной тепловой конвекции дистиллированной воды без изменения ее агрегатного состояния. Для выполнения численного моделирования был выбран период процесса, когда в аппарате находится только магний. Это начальная подготовительная стадия процесса, в ходе которой в аппарат заливаются несколько порций расплавленного магния. После заливки магния в количестве, достаточном для получения заданного циклового съема по титану, выполняется технологический отстой для прогрева и усреднения навески залитого металла. В отечественной промышленной практике в зависимости от цикловой производительности в аппарат может быть залито от 1 до 4 порций магния. Наличие промежуточных технологических операций по транспортировке и заливке магния приводит к тому, что изначально все порции заливаемого металла имеют различную температуру (в пределах 30–40 °С). В результате после выполнения заливки в реторте аппарата располагаются несколько слоев расплавленного металла с различными свойствами. В данной работе проведено моделирование гидродинамики навески магния, состоящей из двух слоев, имеющих различные начальные температуры. Расчетная прямоугольная область высотой 2,937 м и шириной 0,775 м заполнена расплавом магния и условно разделена пополам. Температуры магния принимаются равными 670 и 710 °С, более разогретый слой магния расположен в верхней половине расчетной области. Вертикальные границы изотермические, имеют температуры t1 = 700 °C и t2 = 870 °C, горизонтальные – адиабатные. На рис. 2 представлены расчетная схема, изотермы и линии тока, полученные на модели в процессе установления. В результате нагрева образуется двухвихревая структура течения расплавленного металла против часовой стрелки. В разогретом до 710 °С верхнем слое магния образуется более мощный контур течения жидкой фазы, чем в слое, имеющем начальную температуру 670 °С. Разогретый у правой стенки металл выносится в верхнюю часть канала, где охлаждается и опускается вдоль левой границы. В ходе свободной тепловой конвекции наблюдалось постепенное перемешивание холодного и разогретого слоев металла с усреднением температуры по всей высоте расчетной области. Два вихря сливаются, образуя один замкнутый контур. а б в Рис. 2. Расчетная схема (а), изотермы (б) и линии тока (в), полученные на модели в процессе установления в условиях свободной конвекции, расчет на сетке N×M = = 30×115, t1 = 700 °С и t2 = 870 ºС, начальные температуры магния tMg1 = 670 °С и t Mg2 = 710 ºС На рис. 3 представлены изотермы и линии тока после выхода решения на стационарный режим. На следующем этапе работы математическую модель следует развивать в направлении описания тепломассопереноса в блоке губчатого титана, осаждаемого в нижней части реторты и являющегося пористым телом. а б Рис. 3. Изотермы (а) и линии тока (б) расплава магния в условиях свободной конвекции, расчет на сетке N×M = 30×115, t1 = 700 °С и t2 = 870 ºС Таким образом, разработана математическая модель и компьютерная программа, позволяющие анализировать неравновесные теплофизические процессы в реакторе аппарата для получения губчатого титана. На первом этапе работы при помощи численных методов решено дифференциальное уравнение переноса энергии в бесконвективном приближении, получено решение при моделировании различных периодов ведения процесса восстановления. Учет гидродинамических явлений в расплаве магния позволил уточнить температурную обстановку, выявив концентрацию температурных градиентов у стенки реторты. Модель следует развивать с учетом фильтрации жидкой фазы через пористую среду осаждаемого губчатого титана. Разрабатываемая математическая модель может оказаться полезной для прогнозирования динамики загрузки компонентов, условий теплообмена, обеспечивающих экономию энергоресурсов при максимальной производительности процесса получения губчатого титана. Список литературы

Об авторах

Владимир Николаевич Нечаев

ОАО «Российский научно-исследовательский и проектный институт титана и магния», г. Березники

Email: nechvladimir@yandex.ru
618421, г. Березники, Пермский край, пр. Ленина, 101 аспирант ПНИПУ, заведующий лабораторией термии титана и магния, ОАО «Российский научно-исследовательский и проектный институт титана и магния»

Алексей Иванович Цаплин

Пермский национальный исследовательский политехнический университет

Email: tai@pstu.ru
614990, г. Пермь, Комсомольский пр., 29 доктор технических наук, профессор, Пермский национальный исследовательский политехнический университет

Список литературы

Статистика

Просмотры

Аннотация - 56

PDF (Russian) - 26

Ссылки

• Ссылки не определены.