Научная статья на тему 'ОСОБЕННОСТИ КОНВЕКТИВНОГО ТЕПЛОПЕРЕНОСА В РАСПЛАВЕ МНОГОЭЛЕКТРОДНОЙ ДУГОВОЙ ПЕЧИ'

ОСОБЕННОСТИ КОНВЕКТИВНОГО ТЕПЛОПЕРЕНОСА В РАСПЛАВЕ МНОГОЭЛЕКТРОДНОЙ ДУГОВОЙ ПЕЧИ Текст научной статьи по специальности «Физика»

CC BY
91
11
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДУГОВАЯ ПЕЧЬ / МАГНИТНАЯ ГИДРОДИНАМИКА / ЭЛЕКТРОВИХРЕВЫЕ ТЕЧЕНИЯ / ТЕРМОГРАВИТАЦИОННАЯ КОНВЕКЦИЯ / ARC FURNACE / MAGNETOHYDRODYNAMICS / ELECTROVORTEX FLOWS / THERMOGRAVITATIONAL CONVECTION

Аннотация научной статьи по физике, автор научной работы — Кухарев Алексей Леонидович

Современным направлением совершенствования технологии производства стали в высокомощных дуговых печах является интенсификация магнитогидродинамических эффектов для перемешивания расплава. В связи с этим актуальным является исследование особенностей теплопереноса в расплаве печи с учетом взаимодействия электровихревой и термогравитационной конвекции. Результаты получены с помощью трехмерной математической модели магнитогидродинамических и тепловых процессов, построенной при использовании безындукционного приближения с учетом k-e модели турбулентности. В качестве теплогенерирующих источников модель учитывает тепловые потоки от электрических дуг и джоулев нагрев. Обработка результатов произведена с использованием методов визуализации вихревых структур. Предложена конструкция печи, содержащая три сводовых и три подовых электрода, обеспечивающая формирование дополнительных электровихревых течений в расплаве. Показано, что при заданных условиях моделирования и токах в электродах 80 кА в результате взаимодействия электровихревой и термогравитационной конвекции в расплаве печи формируется многовихревая структура течений. Установлено, что вблизи оси ванны доминирует электровихревая конвекция, термогравитационная конвекция, обусловленная неравномерным нагревом расплава, приводит к уменьшению размеров основного электровихревого течения и формированию дополнительного течения вблизи боковых стенок печи. Зафиксированы максимальные значения скорости 2 м/с в областях расплава, находящихся под электрическими дугами. При этом скорость нисходящего течения под электрическими дугами уменьшается, а скорость восходящего течения в области подовых электродов увеличивается. Выявлено влияние термогравитационной конвекции на азимутальные течения расплава в основном в области подовых электродов, что приводит к увеличению азимутальной скорости и смещению центров вихрей к центру ванны. Верификация предложенной модели выполнена путем сопоставления результатов расчета с данными экспериментов, полученных на лабораторных установках с различными схемами расположения электродов. Результаты предполагается использовать для дальнейшего улучшения энергетических и конструктивных параметров дуговых печей.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

CHARACTERISTICS OF CONVECTIVE HEAT TRANSFER IN THE MELT OF MULTI-ELECTRODE ARC FURNACE

The modern direction of improving the technology of steel production in high-power arc furnaces is the intensification of magnetohydrodynamic effects for melt mixing. In this regard, it is relevant to study the characteristics of heat transfer in the melt of this furnace, taking into account the interaction of electrovortex and thermogravitational convection. The results were obtained using a three-dimensional mathematical model of magnetohydrodynamic and thermal processes, constructed using a non-inductive approximation, taking into account the k-e turbulence model. As heat-generating sources, the model takes into account the heat flows from electric arcs and Joule heating. Processing of the results was carried out using visualization methods of vortex structures. A furnace design containing three arched and three bottom electrodes and providing the formation of additional electrovortex flows in the melt is proposed. It is shown that under the given simulation conditions and currents in 80 kA electrodes a multivortex flow is formed in the furnace melt as a result of the interaction of electrovortex and thermogravitational convection. Electrovortex convection dominates near the bath axis. Thermogravitational convection, due to uneven heating of the melt, leads to a reduction in the size of the main electrovortex flow and the formation of an additional flow near the side walls of the furnace. Maximum speeds of 2 m/s are fixed in the melt areas under electric arcs. In this case, the speed of the downward flow under the electric arcs decreases, and the speed of the upward flow in the region of the bottom electrodes increases. The effect of thermogravitational convection on the azimuthal melt flow is manifested mainly in the region of the bottom electrodes, leading to an increase in the azimuthal velocity and displacement of the vortices to the center of the bath. The verification of the proposed model was carried out by comparing the calculation results with the experimental data obtained in laboratory installations with various electrode arrangements. The results will be used to further improve the energy and design parameters of arc furnaces.

Текст научной работы на тему «ОСОБЕННОСТИ КОНВЕКТИВНОГО ТЕПЛОПЕРЕНОСА В РАСПЛАВЕ МНОГОЭЛЕКТРОДНОЙ ДУГОВОЙ ПЕЧИ»

15. Chuchueva, I.A., Inkina, N.E. Nauka i obrazovanie: nauchnoe izdanie MGTU im. N.E. Baumana, 2015, no. 8, рр. 195-238.

16. Heinhold, I. Ingeniur statistic. München; Wien: Springier Verlag, 1964. 352 p.

Бубнов Кирилл Николаевич,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», магистрант, e-mail: kirill.bubnov.96@mail.ru Bubnov Kirill Nikolayevich,

Ivanovo State Power Engineering University, Master student, e-mail: kirill.bubnov.96@mail.ru Барочкин Алексей Евгеньевич,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры тепловых электрических станций, телефон (4932) 26-99-13, e-mail: acorp27@yandex.ru Barochkin Aleksei Yevgenyevich,

Ivanovo State Power Engineering University, Candidate of Engineering Sciences (PhD), Associate Professor, Thermal Power Plants Department, telephone (4932) 26-99-13, e-mail: acorp27@yandex.ru

Жуков Владимир Павлович,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор, заведующий кафедрой прикладной математики, телефон (4932) 26-97-45, e-mail: zhukov-home@yandex.ru Zhukov Vladimir Pavlovich,

Ivanovo State Power Engineering University, Doctor of Engineering Sciences, Professor, Head of Applied Mathematics Department, telephone (4932) 26-97-45, e-mail: zhukov-home@yandex.ru

Ледуховский Григорий Васильевич,

ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, доцент, заведующий кафедрой тепловых электрических станций, телефон (4932) 26-99-31, e-mail: lgv83@yandex.ru Ledukhovsky Grigory Vasilievich,

Ivanovo State Power Engineering University, Doctor of Engineering Sciences, Head of Thermal Power Plants Department, telephone (4932) 26-99-31, e-mail: lgv83@yandex.ru

УДК 669.187.2; 537.84

ОСОБЕННОСТИ КОНВЕКТИВНОГО ТЕПЛОПЕРЕНОСА В РАСПЛАВЕ МНОГОЭЛЕКТРОДНОЙ ДУГОВОЙ ПЕЧИ

А.Л. КУХАРЕВ

ГОУВПО «Луганский национальный университет имени Владимира Даля», г. Луганск, Украина

E-mail: alex.kuharev@mail.ru

Авторское резюме

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

Материалы и методы. Результаты получены с помощью трехмерной математической модели магнитогидродинамических и тепловых процессов, построенной при использовании безындукционного приближения с учетом k-e модели турбулентности. В качестве теплогенерирующих источников модель учитывает тепловые потоки от электрических дуг и джоулев нагрев. Обработка результатов произведена с использованием методов визуализации вихревых структур.

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

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

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

Ключевые слова: дуговая печь, магнитная гидродинамика, электровихревые течения, термогравитационная конвекция

CHARACTERISTICS OF CONVECTIVE HEAT TRANSFER IN THE MELT OF MULTI-ELECTRODE ARC FURNACE

A.L. KUKHAREV Lugansk Vladimir Dalh National University, Lugansk, Ukraine E-mail: alex.kuharev@mail.ru

Abstract

Background. The modern direction of improving the technology of steel production in high-power arc furnaces is the intensification of magnetohydrodynamic effects for melt mixing. In this regard, it is relevant to study the characteristics of heat transfer in the melt of this furnace, taking into account the interaction of electrovortex and thermogravitational convection.

Materials and methods. The results were obtained using a three-dimensional mathematical model of magnetohydrodynamic and thermal processes, constructed using a non-inductive approximation, taking into account the k-e turbulence model. As heat-generating sources, the model takes into account the heat flows from electric arcs and Joule heating. Processing of the results was carried out using visualization methods of vortex structures.

Results. A furnace design containing three arched and three bottom electrodes and providing the formation of additional electrovortex flows in the melt is proposed. It is shown that under the given simulation conditions and currents in 80 kA electrodes a multivortex flow is formed in the furnace melt as a result of the interaction of electrovortex and thermogravitational convection. Electrovortex convection dominates near the bath axis. Thermogravitational convection, due to uneven heating of the melt, leads to a reduction in the size of the main electrovortex flow and the formation of an additional flow near the side walls of the furnace. Maximum speeds of 2 m/s are fixed in the melt areas under electric arcs. In this case, the speed of the downward flow under the electric arcs decreases, and the speed of the upward flow in the region of the bottom electrodes increases. The effect of thermogravitational convection on the azimuthal melt flow is manifested mainly in the region of the bottom electrodes, leading to an increase in the azimuthal velocity and displacement of the vortices to the center of the bath.

Conclusions. The verification of the proposed model was carried out by comparing the calculation results with the experimental data obtained in laboratory installations with various electrode arrangements. The results will be used to further improve the energy and design parameters of arc furnaces.

Key words: arc furnace, magnetohydrodynamics, electrovortex flows, thermogravitational convection

DOI: 10.17588/2072-2672.2020.2.013-022

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

превышает 100 МВА при токах в электродах 50-150 кА [1-2]. Для снижения нагрузок на электроды и распределения вводимой электрической мощности такие печи, как правило, выполняются многоэлектродными, причем для усиления магнитогидроди-

намических эффектов при перемешивании расплава предложены различные схемы расположения сводовых и подовых электродов [1, 3].

В патенте [4] нами предлагается конструкция дуговой печи, которая содержит три подовых и три сводовых электрода, причем угол между смежными осями сводовых и подовых электродов составляет 60°. Такая компоновка электродов в ванне, как показано в [5], при питании постоянным током или током низкой частоты (0,01-1 Гц) обеспечивает подавление скин-эффекта и улучшение равномерности распределения электромагнитных сил в расплаве.

В [6] показано, что основными механизмами теплопереноса в расплаве таких печей является электровихревая и термогравитационная конвекция. Однако вопросы взаимодействия электровихревой и термогравитационной конвекции для многоэлектродных печей все еще недостаточно проработаны.

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

Методы исследования. Предлагаемая шестиэлектродная дуговая печь в данном исследовании имела ванну цилиндро-сфероконической формы емкостью 180 т. Геометрическая модель и профиль ванны представлены на рис. 1 с указанием основных геометрических размеров ванны.

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

Шлаковый слой в модели не учитывался, электроды и дуги учитывались только при расчете электромагнитного поля [5]. В гидродинамических и тепловых расчетах влияние электрических дуг учитывалось посредством введения соответствующих граничных условий в области менисков.

т

6350

б)

Рис. 1. Геометрическая модель ванны печи: а - общий вид; б - профиль ванны в направлении оси У; 1, 2, 3 - сводовые электроды; 4, 5, 6 - подовые электроды

При разработке математической модели в качестве основного допущения принято магнитогидродинамическое приближение, при котором рассматриваемые процессы описываются системой уравнений магнитной гидродинамики [1, 7], при этом не учитываются химические реакции, а среды в пределах выделенных зон считаются однородными, изотропными и немагнитными. Так как магнитное число Рей-нольдса Рет <1, то перенос магнитного поля потоком электропроводной жидкости не учитывался [1, 7].

Ввиду различных особенностей ведения технологического процесса в печах на данном этапе не учитывался теплообмен между шлаком и расплавом, а также влияние топливно-кислородных горелок. Так как число Рейнольдса для течений расплава составляет Ре = 104-105, то для описания турбулентности используем к-е модель, которая при ее относительной простоте обес-

печивает наиболее близкие к экспериментальным данным распределения скоростей [1, 8]. При принятых допущениях базовые уравнения модели имеют следующий вид: • уравнение энергии

'дГ

РСР — ] = "РСр 0? • V7") + X3-VZT + Q;

(1)

• уравнение Навье - Стокса

Р

dv

~dt

= -Vp + rigV v + pg + F; (2)

• уравнение неразрывности потока У-\? = 0; (3)

• уравнения к-г модели турбулентно-

сти

д(рк)

+ V (pvk) = V

д t

+Gk+Gb- рв; <3(ps)

у

д t

- + V (pvs) = V

Г)

ц-

Цт

к

\

Цт

Vk

Vs

(4)

(5)

+C1s s (Gk + C3sGb ) - C2sP у;

уравнение Колмогорова-Прандтля

Цт = PCL

к2

(6)

• уравнение для турбулентной теплопроводности

_ СрПт.

Хт =

PrT

(7)

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

Р = Р0[1-Р(Г-Г0)], (8)

где переменные модели: V- скорость движения расплава, м/с; р - давление, Па; Т -температура, К; + -эффективный

коэффициент динамической вязкости, Пас; = А + А у - эффективный коэффициент

теплопроводности расплава, Вт/(мК); к -турбулентная кинетическая энергия, м2с2; г - диссипация кинетической энергии турбулентности, м2/с3; Ок - специфическая генерация турбулентности от средних градиентов скорости; Оь - специфическая генерация турбулентности от силы Архимеда; С3г - коэффициент, зависящий от направлений потока и гравитационной силы.

Силы Лоренца в расплаве F = JxB и

плотность джоулева нагрева Q =

определялись на основании расчета электромагнитного поля, базовые уравнения и граничные условия для электромагнитного поля подробно описаны в [5]. Моделирование печи проводилось при значениях токов в электродах 80 кА.

Характеристики электрических дуг определялись нами на основании численного решения уравнения Эленбааса - Геллера [9].

Постоянные параметры модели, определяемые свойствами расплава жидкой стали и постоянными модели турбулентности, приведены в табл. 1.

Таблица 1. Постоянные параметры модели

Параметр Обо- Единица Значе-

зна- измере- ние

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

чение ния

Плотность Ро кг/м3 6900

Характерная температура То К 1900

Коэффициент

динамическои Ц ПаС 0,007

вязкости

Коэффициент

молекулярной теплопровод- X Вт/(мК) 35

ности

Теплоемкость Ср Дж/(кгК) 792

Коэффициент объемного р К-1 0,00014

расширения

Турбулентное

число Pry 0,85

Прандтля

Cle 1,44

Постоянные C2e 1,92

модели турбулентности Cl 0,09

CTs 1,0

CTfc 1,3

Используемые граничные условия (ГУ) приведены в табл. 2. Отметим, что в мощных дуговых печах в ходе плавки формируют шлаковый слой, покрывающий верхнюю поверхность расплава. Таким образом, на этой поверхности применялось условие теплоизолированности.

2

s

Таблица 2. Граничные условия

Наименова- Накладываемые ГУ

ние поверх- гидроди- тепловые

ности намические

Нижняя и Условие Условие

боковая прилипа- первого рода:

стенки ния: V = 0 Тг = 1900 К

ванны

Области Условие Условие

менисков прилипа- второго рода: Я = 6,43107 Вт/м2

электриче- ния: V = 0

ских дуг [8]

Верхняя по- Условие Условие тепло-

верхность прилипа- изолированности:

расплава ния: V = 0 Яг = 0

Результаты исследования. На первом этапе были исследованы течения расплава без учета термогравитационной конвекции, при этом уравнение (8) не использовалось, а в уравнениях (4) и (5) полагалось Оь = 0.

Результаты моделирования течений в вертикальном сечении ванны, проходящем через ось У, без учета термогравитационной конвекции приведены на рис. 2, где сплошными линиями показаны линии тока (линии скорости), а стрелками - направления потоков.

Рис. 2. Траектории течения расплава в вертикальном сечении ванны вдоль оси У без учета термогравитационной конвекции

Анализ полученных результатов (рис. 2) показывает, что структура течений расплава - многовихревая. В зоне под электрическими дугами формируются мощные струи расплава, направленные под углом приблизительно 20° (относительно осей электрических дуг) к подине печи. Ширина этих струй в средней по глубине части ванны составляет около 800 мм. Максимальная скорость в струе расплава достигает 2,2 м/с. В нижней части ванны расплав растекается к периферии со скоростями 0,2-0,5 м/с. Далее расплав медленно (V < 0,2 м/с) поднимается на поверхность на участке, удаленном приблизительно на 2500 мм от оси ванны, и далее

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

В областях между электрическими дугами формируются восходящие вихревые потоки с максимальной скоростью вблизи подовых электродов (0,5-0,7 м/с). Ширина струи восходящего течения в средней части ванны составляет около 1200 мм. Отчетливо наблюдается второй контур циркуляции, в котором расплав вблизи оси ванны движется вверх, а на поверхности медленно движется к периферии.

Кроме того, в районе пятна дуги между первым и вторым контуром имеется малый вихрь, движущийся в противоположную сторону относительно первого контура циркуляции. Под электрической дугой течение в этом контуре направлено вниз, но уже примерно в середине глубины ванны оно разворачивается вверх под действием течения второго вихря. На поверхности расплава его движение направлено от центра ванны к области электрической дуги с максимальной скоростью около 0,5-1 м/с непосредственно под мениском дуги.

Учет в модели термогравитационной конвекции, формируемой вследствие неравномерного нагрева расплава (рис. 3), несколько изменяет структуру течений.

3 400 3 200 3 000 2 800 2 600 2 400 2 200 2 000 1 800

г -----У

-2

3

У, м

а)

т, к

б)

Рис. 3. Графики распределения температуры вдоль оси У: 1 - при Ж = -70 мм; 2 -Ж = -310 мм; 3 - Ж = -715 мм; 4 -Ж = -1320 мм

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

Рис. 4. Упрощенная схема взаимодействия электровихревой и термогравитационной конвекции

Также для оценки влияния термогравитационной конвекции нами использовался критерий Ричардсона:

% (9)

га = ■

Ре

График распределения локального числа Ричардсона вдоль диаметра ванны приведен на рис. 5. Его анализ показывает, что вблизи оси ванны Р « 1, это говорит о том, что влияние электровихревой конвекции велико, а влияние сил плавучести незначительно. Вблизи боковых стенок Р > 3 (перепад температур около 50 К), что на структуре течений в вертикальной плоскости проявляется в уменьшении размеров первого контура циркуляции и формировании четвертого контура (рис. б). При этом интенсивность течений во втором и третьем контурах увеличивается.

Рис. 6. Траектории течений в вертикальном сечении ванны вдоль оси У с учетом термогравитационной конвекции

Численное изменение скоростей при учете термогравитационной конвекции поясняется графиками распределения скорости (рис. 7), где сплошными линиями показаны профили скорости с учетом термогравитационной конвекции, а точечным пунктиром - без учета термогравитационной конвекции.

Рис. 5. Распределение локального числа Ричардсона вдоль оси У при 2 = -75 мм

Рис. 7. Графики распределения скорости вдоль оси У: а - при 2 = -70 мм; б -2 = -715 мм; в - 2 = -1320 мм

Анализ графиков (рис. 7) показывает, что в областях под электрическими дугами значения скоростей нисходящих потоков при учете термогравитационной конвекции

уменьшаются, а скорости восходящих потоков в приосевой области ванны и в области подовых электродов увеличиваются. В среднем изменение скоростей составляет 15-25 %.

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

а)

б)

Рис. 8. Траектории течений расплава в горизонтальных сечениях ванны с учетом термогравитационной конвекции: а - при Z = -70 мм; б -Z = -1320 мм

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

Нами были проанализированы и сопоставлены графики азимутальных скоростей и электромагнитных сил вблизи поверхности расплава и вблизи подовых электродов. Графики построены по угловой координате Ф вдоль диаметра распада электродов. За начало отсчета 0° принято направление оси X (рис. 8). Вблизи поверхности расплава максимальные значения V,,;, наблюдаются в областях боковых поверхностей пятен электрических дуг (рис. 9,а), где V, составляет около 1,1 м/с (при 83°, 95°, 203°, 215°, 323°, 325°). Изменение знака азимутальной скорости вблизи координат 90°, 210°, 330° соответствует областям расхождения потоков, а вблизи координат 30°, 150°, 270° - областям замыкания потоков.

В нижней части ванны (рис. 9, в) V, достигает 0,3 м/с вблизи боковых поверхностей подовых электродов (вблизи координат 0°, 60°, 120°, 180°, 240°, 300°, 360°), Области расхождения потоков соответствуют координатам 90°, 210°, 330°, а области замыкания потоков - координатам 30°, 150°, 270°.

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

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

ченных на разработанной нами шестиэлек-тродной установке [10], а также на двух-электродных экспериментальных установках Объединенного института высоких температур РАН [7] и Массачусетского технологического университета [8].

1.5

0.5

Г

[

ф, фзд-

а)

б)

в)

10000 7500 5000 2500 0

-2500 -5000 -7500 -Ю000

1 1

\ \ \

\ \ \

\ __J \ \ J

\ < " \ г \ Г

\ \ \

\ \ \

1 __„

ф, град

г)

Рис. 9. Графики азимутальных компонент скорости и электромагнитных сил: а, б - при Z = -70 мм; в, г - Z = -1320 мм;

Отмечено удовлетворительное качественное и количественное совпадение результатов расчета скорости течения и температуры, при чем среднее отклонение экспериментальных данных от расчетных не превышает 15 %. В качестве примера на рис. 10 показано сопоставление расчетных и экспериментальных данных на графике распределения скорости по глубине ванны вдоль ее оси (результаты получены на установке Объединенного института высоких температур, сплав - индий-галлий-олово, температура - 283 К, ток - 250 А [7]).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Рис. 10. График распределения скорости по

глубине ванны вдоль ее оси: - - расчет;

• - эксперимент

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

Выводы. На основании численного исследования конвективного теплоперено-са в предложенной шестиэлектродной печи при питании ее постоянным током показано, что структура течений в такой печи -многовихревая. На основании полученного распределения локального числа Ричардсона показано, что вблизи оси ванны доминирует электровихревая конвекция, а вблизи боковых стенок - термогравитационная конвекция, обусловленная неравномерным нагревом расплава.

Влияние термогравитационной конвекции обусловливает уменьшение шири-

ны основного электровихревого течения и формирование дополнительного течения вблизи боковых стенок печи, при этом скорость нисходящего течения под электрическими дугами уменьшается, а скорость восходящих потоков в области подовых электродов увеличивается. В среднем указанное изменение скоростей составляет 15-25 %.

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

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

1. Недопекин Ф.В., Казак О.В., Семко А.Н.

Моделирование гидродинамических и тепло-физических процессов в сталеплавильных печах постоянного тока. - Донецк: ДонНУ, 2014. -315 с.

2. Адати Т., Селлан Р. Сверхмощная 420-тонная электродуговая печь компании Tokyo Steel, Япония // Металлургическое производство и технология. - 2012. - № 2. - С. 8-17.

3. Multi-Electrode Arc Furnace Technology with Improved Metal Processing Capability Using Current Driven Mixing / H. Molder, J. Jarvik, T. Vaimann, R. Gordon // Proc. 11th International Conference on Environment and Electrical Engineering (EEEIC). - Venice, 2012. - Р. 313-316. DOI: 10.1109/EEEIC. 2012.6221394

4. Пат. 119279 Украина МПК Н05В 7/144, F27D 11/10. Система электропитания многоэлектродной дуговой печи / А.Л. Кухарев, А.А. Ковальчук; опубл. 25.09.2017, Бюл. № 18.

5. Кухарев А.Л. Электромагнитное поле в ванне многоэлектродной печи // Вестник Луганского национального университета имени Владимира Даля. - 2018. - № 5(11). - С. 247-251.

6. Кухарев А.Л. Течения расплава в многоэлектродной дуговой печи в условиях электровихревой и естественной конвекции // Вестник Луганского национального университета имени Владимира Даля. - 2019. - № 7(25). -С.217-223.

7. Ивочкин Ю.П. Исследование механизмов термогидродинамических и МГД процессов с жидкометаллическими рабочими телами: дис. ... д-ра техн. наук: 01.04.14; Нац. исслед. ун-т МЭИ. - М., 2015. - 407 с.

8. Murthy A., Szekely J., El-Kaddah N. Experimental measurement and numerical compu-

tation of velocity and turbulence parameters in a heated liquid metal system // Metallurgical Transactions B. - 1988. - Vol. 19, No. 0. - Р. 765-775.

9. Кухарев А.Л., Корсунов К.А., Серги-енко С.Н. Расчет параметров электрической дуги в электропечных установках для производства стали и ферросплавов // Сборник научных трудов Донбасского государственного технического университета. - 2017. - № 5(48). -С. 52-57.

10. Кухарев А.Л., Петренко В.В. Экспериментальная многоэлектродная установка для исследования энергетических процессов в токо-проводящих расплавах // Ресурсосберегающие технологии производства и обработки давлением материалов в машиностроении. - 2018. -№ 4(25). - С. 64-71.

11. Дементьев С.Б., Чайковский А.И., Чудновский А.Ю. Формирование электровихревых течений в ваннах с многоэлектродным токоподводом // Магнитная гидродинамика. -1988. - № 1. - С. 85-89.

References

1. Nedopekin, F.V., Kazak, O.V., Semko, A.N. Modelirovanie gidrodinamicheskikh i teplofizi-cheskikh protsessov v staleplavil'nykh pechakh postoyannogo toka [Modeling of hydrodynamic and thermophysical processes in direct current steelmaking furnaces]. Donetsk: DonNU, 2014. 315 p.

2. Adati, T., Sellan, R. Sverkhmoshchnaya 420-tonnaya elektrodugovaya pech' kompanii Tokyo Steel, Yaponiya [Ultra high power 420-ton electric arc furnace at Tokyo Steel, Japan]. Metal-lurgicheskoe proizvodstvo i tekhnologiya, 2012, no. 2, pp. 8-17.

3. Molder, H., Jarvik, J., Vaimann, T., Gordon, R. Multi-Electrode Arc Furnace Technology with Improved Metal Processing Capability Using Current Driven Mixing. Proc. 11th International Conference on Environment and Electrical Engineering (EEEIC). Venice, 2012, pp. 313-316. doi: 10.1109/EEEIC. 2012.6221394

4. Kukharev, A.L., Koval'chuk, A.A. Sistema elektropitaniya mnogoelektrodnoy dugovoy pechi [Power supply system of multi-electrode arc furnace]. Patent Ukraina, no. 119279, 2017.

5. Kukharev, A.L. Elektromagnitnoe pole v vanne mnogoelektrodnoy pechi [Electromagnetic field in the bath of multielectrode furnace]. Vestnik Luganskogo natsional'nogo universiteta imeni Vla-dimira Dalya, 2018, no. 5(11), pp. 247-251.

6. Kukharev, A.L. Techeniya rasplava v mnogoelektrodnoy dugovoy pechi v usloviyakh el-ektrovikhrevoy i estestvennoy konvektsii [Melt flow in the multielectrode arc furnace under conditions of electro-vortex and natural convection]. Vestnik Luganskogo natsional'nogo universiteta imeni Vladimi-ra Dalya, 2019, no. 7 (25), pp. 217-223.

7. Ivochkin, Yu.P. Issledovanie mekhanizmov termogidrodinamicheskikh i MGD protsessov s zhidkometallicheskimi rabochimi telami. Diss. dokt. tekhn. nauk [Study of the mechanisms of thermohy-drodynamic and MHD processes with liquid metal working bodies. Dr. tech. sci. diss.]. Moscow, 2015. 407 p.

8. Murthy, A., Szekely, J., El-Kaddah, N. Experimental measurement and numerical computation of velocity and turbulence parameters in a heated liquid metal system. Metallurgical Transactions B, 1988, vol. 19, no. 10, pp. 765-775.

9. Kukharev, A.L., Korsunov, K.A., Sergien-ko, S.N. Raschet parametrov elektricheskoy dugi v elektropechnykh ustanovkakh dlya proizvodstva stali i ferrosplavov [Calculation of electric arc parameters in electric furnace for the production of steel and ferroalloys]. Sbornik nauchnykh trudov

Donbasskogo gosudarstvennogo tekhnicheskogo universiteta, 2017, no. 5(48), pp. 52-57.

10. Kukharev, A.L., Petrenko, V.V. Eksperi-mental'naya mnogoelektrodnaya ustanovka dlya issledovaniya energeticheskikh protsessov v to-koprovodyashchikh rasplavakh [Experimental mul-tielectrode installation for the study of energy processes in conducting melts]. Resursosberegay-ushchie tekhnologii proizvodstva i obrabotki davleniem materialov v mashinostroenii, 2018, no. 4(25), pp. 64-71.

11. Dement'ev, S.B., Chaykovskiy, A.I., Chudnovskiy, A.Yu. Formirovanie el-ektrovikhrevykh techeniy v vannakh s mnogoel-ektrodnym tokopodvodom [Generation of electro-vortex flows in liquid-metal baths with a multielec-trode current input]. Magnitnaya gidrodinamika, 1988, no. 1, pp. 85-89.

Кухарев Алексей Леонидович,

ГОУВПО «Луганский национальный университет имени Владимира Даля», кандидат технических наук, доцент кафедры горной электромеханики и транспортных систем, e-mail: alex.kuharev@mail.ru Kukharev Aleksei Leonidovich,

Lugansk Vladimir Dahl National University, Candidate of Engineering Sciences (PhD), Associate Professor of Mining Electromechanics and Transport Systems Department, e-mail: alex.kuharev@mail.ru

УДК 66.047

УСЛОВИЯ И КОНСТРУКТИВНЫЕ МЕРЫ МАСШТАБНОГО ПЕРЕХОДА ОТ ЛАБОРАТОРНЫХ УСТАНОВОК К ПРОМЫШЛЕННЫМ ОБРАЗЦАМ ПРИ ПРОЕКТИРОВАНИИ СУШИЛЬНЫХ АГРЕГАТОВ ВИХРЕВОГО ТИПА

А.И. СОКОЛЬСКИЙ, Е.В. ГУСЕВ, О.Б. КОЛИБАБА ФГБОУВО «Ивановский государственный энергетический университет имени В.И. Ленина»,

г. Иваново, Российская Федерация E-mail: tevp@tvp.ispu.ru; gusev_pcm@mail.ru; koli-baba@mail.ru

Авторское резюме

Состояние вопроса. Аппараты вихревого и циклонного типов находят все большее распространение для интенсивного протекания совмещенных технологических процессов обработки (тонкого измельчения, классификации, сушки, смешения, гидрофобизации и др.) дисперсных материалов. Одним из видов данного класса оборудования с активным гидродинамическим режимом являются аппараты с восходящим (или с нисходящим) спирально-закрученным газодисперсным потоком. Исследованию подобных устройств посвящены работы А.Ф. Фролова, В.И. Муштаева, В.М., Ульянова, Б.С. Сажина, которые подтверждают высокую эффективность сушки высоковлажных дисперсных материалов. Одной из задач масштабного перехода от лабораторных установок к промышленным образцам является разработка конструктивных особенностей, определяющих подобие гидродинамических режимов протекающих тепломассообменных процессов.

Материалы и методы. Для определения гидродинамических характеристик газодисперсного потока использованы методы, основанные на измерении расходов, скоростей, гидравлических сопротивлений. Для снятия кривых сушки и температурных кривых полидисперсного материала (золы гидроудаления ТЭС) различной фракции при различных условиях конвективного подвода теплоты применен

i Надоели баннеры? Вы всегда можете отключить рекламу.