Математический анализ физико-химических процессов в электродуговой сталеплавильной печи с подачей материалов через донные фурмы

Размещено на http://www.allbest.ru/

УДК 658.512.011; 519.711.3; 669.18

Математический анализ физико-химических процессов в электродуговой сталеплавильной печи с подачей материалов через донные фурмы

А.А. Арсеньева

Описаны результаты математического анализа процесса получения жидкого железа прямого восстановления, выполненного на базе уравнения термодинамического состояния веществ. Определены параметры оптимизации печи и установлена наиболее рациональная конструкция.

Ключевые слова: физико-математическое моделирование, железо прямого восстановления, тепло- массоперенос, оптимизация печи.

The results of mathematical analysis of the process of producing liquid DRI based on the thermodynamic state equation substances. The parameters of optimization of furnaces and the most rational design.

Keywords: physical-mathematical modeling, direct reduced iron, heat mass transfer, optimization of the furnace.

Проблемы отечественной металлургии требуют решения задач, направленных на разработку эффективных технических решений, поиск оптимальных режимов управления существующими процессами, а также создание принципиально новых технологий и конструкций агрегатов. Одним из новых является процесс одновременного получения железа прямого восстановления и газификации угля в электродуговой сталеплавильной печи с загрузкой материала через донные фурмы, рис.1

ДСП имеет корпус, в котором наводится металлическая и шлаковая ванны, три графитовых электрода, три донные фурмы для подачи оксида железа, три фурмы для подачи угольной пыли сквозь толщу металлической ванны, и фурма для подачи кислорода.

Так как экспериментальное исследование в металлургии исключительно дорого, то до начала выполнения работ по созданию предложенного агрегата необходимо выполнить теоретический анализ, целью которого является определение оптимального расположения его элементов, в первую очередь электродов и фурм.

Критериями оценки являются термодинамическое состояние материалов в печи и химический состав расплава. Для обеспечения максимальной энергетической эффективности необходимо минимизировать рабочую температуру агрегата. Так как эта температура ограничена температурой плавления, то задача сводится к определению расположения элементов, при котором распределение температуры будет равномерным. Эта задача решена методом компьютерного моделирования физико-химических явлений. железо термодинамический печь

Так как процесс протекает во времени и пространстве, то термодинамическое состояние непрерывно изменяется вследствие выделения теплоты электрической дугой и химическими реакциями.

Химический состав также непостоянен вследствие протекающих реакций и подачи элементов во время плавки. Физико-математическое моделирование процесса заключается в решении системы дифференциальных уравнений энергии и массопереноса.

Рис.1. Устройство дуговой сталеплавильной печи энергометаллургического комплекса: 1 - корпус печи; 2 -электроды; 3 - фурмы для подачи железной руды; 4 - фурмы для подачи угольного порошка; 5 - летка слива железа прямого восстановления; 6 - летка слива шлака; 7 - кислородная фурма

При описании физических явлений использована ортогональная декартова система координат.

Пространство моделирование условно разделено на области: E - область угольных электродов; D - область электрических дуг; R - область расплавленного шлака; M - область расплава металла; F - область футеровки печи; G - область газовой среды. Каждая из областей в зависимости от свойств веществ и процессов, происходящих в ней, описывается специальной системой уравнений.

Движение расплава. Подача двуокиси железа и угольного порошка, выделения угарного газа при их химическом взаимодействии ведет к возникновению давления. Это давление вызывает перемещение расплава. Течение жидкости описывается уравнением Навье-Стокса, которое в декартовой системе координат имеет вид [1]

,

, (1)

,

где vx, vy, vz - составляющие скорости течения в направлении соответствующих координат, - плотность расплава, p - давление в данной точке пространства, - динамическая вязкость.

Распределение давления в расплаве определяется при решении уравнения неразрывности

, (2)

где E - модуль упругости.

Плотность расплавам , содержащего двуокись железа, угольный порошок и угарный газ, рассчитывается по концентрации С этих компонент в расплаве

, (3)

где Fe, FeO, C, СО - плотности расплава, двуокиси железа, углерода и угарного газа.

Начальные условия для решения уравнения Навье-Стокса: t=0, vx=0, vy=0, vz=0, p=pj.

Граничные условия уравнения Навье-Стокса:

- на поверхностях FRFM соприкосновения расплава с футеровкой и с металлошихтой ХRХM принято условие прилипания - vx=0, vy=0, vz=0.

На поверхности соприкосновения расплава с газовой средой GRGM и областью дуги DRDM принята свободная граница для движения расплава -

Подача материалов в области действия донных фурм учтена заданием на диаметрах dF, dC отверстий фурм скоростей движения расплава, равной значению -

,

где G - массовый поток подаваемого материала, кг/с, - плотность материала.

Распределение концентрации элементов в шлаковой и металлической ваннах. Концентрация элементов в расплаве изменяется при плавлении компонентов шихты и перехода ее элементов в расплав, химического взаимодействия элементов в расплаве и перемещения продуктов взаимодействия потоком жидкости. Изменение концентрации C всех элементов расплава описывается уравнением переноса. Это изменение определяется минимальной концентрацией компонента в соответствии с химической реакцией С(t)=min(C1(t), C2(t). Для основной реакции восстановления железа Fe2O3+3C=2Fe+3CO.

а) изменение концентрации оксида железа

;

б) изменение концентрации углерода

; (5)

в) изменение концентрации угарного газа

где D - коэффициент диффузии данного элемента в жидком железе, vx, vy, vz - скорости движения расплава, определяемые из решения уравнения Навье-Стокса, С(t) - изменение концентрации вещества вследствие течения химических реакций, A - скорость химической реакции (1 с-1), R(t) - скорость изменения концентрации.

Граничные условия уравнения переноса:

на поверхностях соприкосновения расплава со стенками печи и газовой средой используется условие непроницаемости этих поверхностей для жидких компонент расплава

для газообразных компонент расплава на поверхности соприкосновения расплава с газовой средой используется условие полного удаления газа из расплава

.

на выходе донных фурм заданы потоки двуокиси железа и углерода

.

Начальные условия: в момент начала цикла (после слива части металла и шлака, когда подача материалов в печь прекращается) принято, что в печи нет непрореагировавших двуокиси железа и углерода

Удельные мощности выделения и поглощения теплоты. Основным источником теплоты является электрическая дуга. Подача холодных материалов (двуокиси железа и угольного порошка) рассматривается как сток теплоты. Большое значение имеет теплота химических реакций, главной из которых является эндотермическая реакция восстановления железа углеродом из двуокиси железа.

Электрическая дуга. Допустимо считать, что вся мощность дуги равномерно выделяется в круге, радиус которого больше радиуса электродов Re на длину дуги:

, (7)

где Ue, Ie - ток и напряжение дуги, gradU - градиент потенциала в столбе дуги, равный 3 В/мм для дуги, горящей в парах металла.

Допустимо считать, что удельная теплота взаимодействия кислорода, подающегося через сопло, и железа, находящегося в расплаве, равномерно выделяется в круге.

(8)

Расход теплоты на нагревание поступающих материалов учитывается автоматически при решении уравнения теплопроводности скоростными слагаемыми теплопереноса и граничными условиями уравнения Навье-Стокса, которые задают скорость движения материала через донные фурмы.

Реакция восстановления железа углеродом с выделением угарного газа идёт в объёме металлической M ванны. Интенсивность реакции определяется концентрациями двуокиси железа и углерода в рассматриваемой точке металлической ванны.

Длительность химической реакции мала по сравнению с общим временем плавки. Поэтому интенсивность поглощения теплоты определяется концентрацией элементов в расплаве. Удельная мощность поглощения теплоты реакции

(9)

где Fe2O3, C - распределение концентрации реагента, кг/м3, Q - энергия химической реакции, Дж/кг.

Удельная мощность выделения теплоты реакции окисления

(10)

Теплота эндотермической реакции восстановления составляет 3,07 МДж в расчете на 1 кг Fe2O3.

Общая интенсивность источника теплоты.

Модель термодинамического состояния и теплопереноса. Во всех указанных областях печи протекает нестационарный термодинамический процесс, который описывается изменением энтальпии H(t) множества точек пространства во времени t. Нестационарное линейное уравнение теплопроводности в декартовой системе координат x, y, z имеет вид [1]

, (11)

где T - температура точек пространства, - коэффициент теплопроводности среды, зависящий от координат расположения точки в пространстве, типа вещества и температуры в этой точке, vx, vy, vz - скорости движения вещества в направлении соответствующих координат, qi - удельные значения мощности выделения и поглощения теплоты в данной точке пространства.

(12)

Энтальпия и температура в этом уравнении связаны нелинейными функциями T(H), которые учитывают теплоёмкость и теплоты фазовых и агрегатных превращений вещества в каждой из выделенных зон пространства. Коэффициент теплопроводности зависит от температуры и нужно учитывать его различие в разных зонах печи [2-5].

Начальными условиями принято, что все точки пространства в начальный момент времени имеют одинаковую температуру T0

- .

Граничные условия учитывают теплообмен печи с внешней средой. На внешней поверхности футеровки F0 имеется теплоотдача, создающая в футеровке градиент температуры

-

где b - коэффициент теплоотдачи, F - коэффициент теплопроводности футеровки.

Оптимизация конструкции печи

Для нормального протекания процесса плавления материалов в металлургической печи требуется обеспечить равномерное распределение температуры по объему печи. Равномерное распределение температуры в расплаве железа, необходимое для нормального протекания процесса, обеспечивается при расположении зон выделения теплоты на поверхности расплава над зонами нисходящих вихревых потоков, возникающих в расплаве при подаче руды и угля и струй кислорода. Этого можно достичь путем нахождения оптимального расположения электрических дуг и фурм для подачи реагирующих веществ.

Таким образом, параметрами оптимизации являются: а) оптимальное расположение диаметров центров донных фурм для подачи материала; б) оптимальное расположение диаметров центров электрических дуг.

Определив среднеквадратичную дисперсию температуры, можно выявить диапазон значений оптимальных диаметров центров.

, (13)

где ST - среднеквадратичная дисперсия температуры, V - объем ванны с расплавом, T0 - среднее значение температуры.

. (14)

Оптимальное расположение электрических дуг и донных фурм для подачи материала определено решением вариационной задачи

.

Распределение зависимости дисперсии температуры от конструктивных характеристик показано на рис. 2 и рис. 3.

Рис. 2. Зависимость дисперсии температуры от диаметров центров фурм

Рис. 3. Зависимость дисперсии температуры от диаметров центров электрических дуг

Конечные результаты компьютерного моделирования представлены текущими распределениями по объёму расплава: интенсивности объёмного тепловыделения; энтальпии и температуры; плотности расплава; скоростей движения расплава в направлении каждой из координат; концентраций каждого из веществ.

Распределение давления на поверхности расплава определяется расположением фурм для подачи кислорода. В донной зоне основным источником давления являются потоки порошка углерода и оксида железа. Различие давлений создаёт потоки в расплаве.

Сумма составляющих течения в направлении ортогональных координат в плоскости поперечного сечения является вращательным движением расплава, обусловленным подачей струй кислорода. Подача порошка углерода и оксида железа вызывает движение расплава в виде трёх вертикальных вихрей [6]. Движение расплава вызывает перемещение углерода и оксида железа в объёме расплава и определяет распределение их концентрации CcCFeO.

Концентрация углерода и оксида железа быстро убывает при удалении от подающих донных фурм. Взаимодействие этих веществ вызывает выделение монооксида углерода, которое наиболее интенсивно в областях, равноудалённых от донных фурм.

Электрические дуги и окисление углерода нагревают приповерхностную область расплава, соответственно температура поверхностных слоёв выше. В придонные слои поступают холодные углерод и оксид железа, реакция взаимодействия между которыми поглощает тепло, что снижает температуру расплава.

Теплоперенос осуществляется исключительно потоками расплава, которые выравнивают температуру по объёму. Под действием подачи железной руды, углерода через донные фурмы (3, 4, рис. 4в) и струй кислорода (2, рис. 4а, б), возникает давление, которое на поверхности имеет ярко выраженные зоны растяжения (А, рис.4б) и сжатия (Б, рис.4б). Это давление создаёт сложное вихревое движение расплава, направленное вверх от донных фурм и вниз в зоне воздействия струй кислорода. При сдвиге кислородных фурм на 30о, касательно направленных на поверхность расплава, зона тепловыделения реакции окисления углерода кислородом оказывается в зоне нисходящего потока расплава. Этот поток уносит расплав, нагретый этой реакцией и электрическими дугами, в зону эндотермической реакции восстановления железа.

Рис. 4. Распределение температур (а), давления (б) на поверхности расплава и концентрации монооксида углерода (в) на глубине 1/3 толщины расплава при диаметре центров электродов (1) 200 см, диаметре центров кислородных фурм (2) и фурм для подачи угольного порошка (3) и руды (4) 420 см.

Разработанная компьютерная модель гибридного процесса получения жидкого железа прямого восстановления в дуговой сталеплавильной электропечи энергометаллургического комплекса позволила решить задачу оптимизации конструкции элементов ДСП и технологии ведения данного процесса. Полученные результаты численного эксперимента показали, что наиболее рациональной является конструкция дуговой плавильной печи, при которой центры нижних фурмы для подачи измельчённых руды и угля расположены от оси печи на расстоянии 0,71 её внутреннего диаметра. Центры зон воздействия кислородных струй и электрических дуг сдвинуты относительно центров донных фурм на угол 30о. Центры электрических дуг желательно расположить ближе к оси печи, на расстоянии, равном 0,35 - 0,5 внутреннего радиуса печи.

Общие выводы

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

2. Течение расплава описано уравнением Навье-Стокса, в котором движущей силой является гравитационное давление, возникающее при выделении угарного газа при химическом взаимодействии двуокиси железа и угольного порошка, подача которых учтена в граничных условиях.

3. Изменение концентрации всех элементов расплава описывается уравнением переноса, учитывающим диффузию и движение веществ, а также их химическое взаимодействие. На поверхностях соприкосновения расплава со стенками печи и газовой средой использовано условие непроницаемости этих поверхностей для жидких компонент расплава. Для газообразных компонент расплава на поверхности соприкосновения расплава с газовой средой использовано условие полного удаления газа из расплава.

4. Теплоперенос описан нелинейным уравнением теплопроводности, учитывающим кондуктивный и конвективный теплоперенос, в котором энтальпия и температура связаны нелинейными функциями. Уравнение учитывает объёмные источники теплоты, создаваемые реакция восстановления железа углеродом и мощность которых зависит от их концентрации. Граничные условия учитывают воздействие электрических дуг и теплоту реакции окисления железа на поверхности расплава.

5. При выполнении численного эксперимента исследовали влияние расположения фурм для подачи руды, угля, кислорода и электродов на распределения по объёму расплава скоростей его движения, концентраций окиси железа, углерода и моноокиси углерода, интенсивности выделения и поглощения теплоты и температуры расплава.

6. Показано, что равномерное распределение температуры в расплаве железа, необходимое для нормального протекания процесса, обеспечивается при расположении зон выделения теплоты на поверхности расплава над зонами нисходящих вихревых потоков, возникающих в расплаве при подаче руды и угля и струй кислорода.

7. Установлено, что наиболее рациональной является конструкция дуговой плавильной печи энергометаллургического комплекса, при которой центры нижних фурмы для подачи измельчённых руды и угля расположены от оси печи на расстоянии 0,71 её внутреннего диаметра. Центры зон воздействия кислородных струй и электрических дуг сдвинуты относительно центров донных фурм на угол 30о. Центры электрических дуг желательно расположить ближе к оси печи, на расстоянии, равном 0,35 - 0,5 внутреннего радиуса печи.

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

1. Тихонов А. Н., Самарский А. А. Уравнения математической физики. - М.: Наука, 1972.- 735 c.

2. Тихонов, А. Н. Математическое моделирование технологических процессов и метод обратных задач в машиностроении /А. Н. Тихонов, В. Д. Калько, В. Б. Гласко. - М.: Машиностроение, 1990. - 264 с.

3. Белащенко, Д. К. Компьютерное моделирование жидких и аморфных веществ /Д. К. Белащенко. - М.: МИСИС, 2005. - 408 с.

4. Рябов А. В. Современные способы выплавки стали в дуговых печах / А. В. Рябов, И. В. Чуманов, М. В. Шишимиров - М.: Теплотехник, 2007. - 192 с.

5. Modelling, Optimization and Control of an Electric Arc Furnace / Richard MacRosty. - Hamilton: McMaster University, 2005. - 160 p.

6. Арсеньева А.А. Оптимизация конструкции электродуговой печи энергометаллургического комплекса методом компьютерного инженерного анализа // Известия ТулГУ. Технические науки. Вып. 11: в 2 ч. Ч. 1. - Изд-во ТулГУ. - 2014. - С. 142-150.

Размещено на Allbest.ru