Научный журнал
Успехи современного естествознания
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

ИССЛЕДОВАНИЕ АЭРОДИНАМИКИ ПОЛИФРАКЦИОННОГО ПЛОТНОГО СЛОЯ РУДНОГО СЫРЬЯ

Бобков В.И. 1 Орехов В.А. 1
1 ФГБОУ ВО «Национальный исследовательский университет «МЭИ»
В статье рассмотрены вопросы аэродинамики плотного моно- и полифракционного слоя при изотермических и неизотермических условиях. Приведены результаты экспериментального исследования аэродинамического сопротивления шаров различного диаметра, а также их смесей двухфракционных диаметров при различных весовых содержаниях мелкой фракции. Описана схема экспериментальной установки для исследования аэродинамики плотного слоя при высоких температурах, в диапазонах действующих обжиговых конвейерных и агломерационных машин. На примере полифракционной засыпки кускового фосфоросодержащего рудного сырья рассмотрено влияние температурного распределения по высоте на сопротивление слоя в диапазоне температур действующих обжиговых агрегатов. При высокотемпературных процессах, протекающих в зернистых засыпках, и сложностях, обуславливаемых структурой слоя, добавляются факторы, связанные с изменением свойств потока и материала, а в большинстве случаев с изменением структуры слоя под воздействием химических реакций, протекающих при высоких температурах полифракционных засыпок при неизотермических условиях. Представлена методика расчёта сопротивления слоя при неизотермических условиях на основе позонного разделения плотного слоя. Выявлено влияние на сопротивление слоя при высоких температурах свойств потока газа теплоносителя, особенно вязкостью, так как она зависит от состава газового потока и сильно от температуры. Зависимость кинематической вязкости от температуры представлена в виде многочлена второй степени по температуре. Исследована аэродинамика полифракционных засыпок в изотермических и неизотермических условиях на примере смесей шаров и кускового фосфоросодержащего рудного сырья.
рудное сырьё
плотный слой
аэродинамическое сопротивление
температура
структурная модель
полифракционная засыпка
1. Пучков А.Ю., Лобанева Е.И., Култыгин О.П. Алгоритм прогнозирования параметров системы переработки отходов апатит-нефелиновых руд // Прикладная информатика. 2022. Т. 17. № 1.
2. Леонтьев Л.И., Григорович К.В., Костина М.В. Фундаментальные исследования как основа создания новых материалов и технологий в области металлургии. Часть 1 // Известия высших учебных заведений. Черная металлургия. 2018. Т. 61. № 1. С. 11–22.
3. Новичихин А.В., Шорохова А.В. Процедуры управления поэтапной переработкой железорудных отходов горнопромышленных районов // Известия высших учебных заведений. Черная металлургия. 2018. Т. 60. № 7. С. 565–572.
4. Курилин С.П., Соколов А.М., Прокимнов Н.Н. Компьютерная программа для моделирования показателей технического состояния электромеханических систем // Прикладная информатика. 2022. Т. 17. № 2. С. 105–119. DOI: 10.37791/2687-0649-2022-17-2-105-119.
5. Бобков В.И., Мищенко М.Н. Исследование теплофизических характеристик окомкованного фосфатного материала // Современные наукоемкие технологии. 2016. № 7–1. С. 26–29.
6. Бобков В.И. Энергосбережение в технологии сушки материала в плотном слое на основе интенсификации тепломассообмена // Современные наукоемкие технологии. 2015. № 12–4. С. 585–589.
7. Ильин И.В., Лёвина А.И., Калязина С.Е. Function-oriented approach to mining enterprise automation // Прикладная информатика. 2022. Т. 17. № 2. С. 5–19. DOI: 10.37791/2687-0649-2022-17-2-5-19.
8. Tian Y., Qin G., Zhang Y., Zhao L., Yang T. Experimental research on pellet production with boron-containing concentrate. Characterization of Minerals, Metals, and Materials. 2020. P. 91–102.
9. Тимофеева А.С., Никитченко Т.В., Федина В.В. Определение комкуемости железорудной шихты с целью прогнозирования прочностных свойств окатышей // Современные наукоемкие технологии. 2015. № 8. С. 53–57.
10. Бобков В.И., Дли М.И., Панченко С.В. Обобщённая структурно-функциональная модель инжиниринга и управления экологически безопасной переработкой отвалов горно-обогатительных комбинатов апатит-нефелиновых руд // Успехи современного естествознания. 2019. № 9. С. 48–52.
11. Kossoy A. Effect of thermal inertia-induced distortions of DSC data on the correctness of the kinetics evaluated. Journal of Thermal Analysis and Calorimetry. 2021. Т. 143. № 1. С. 599–608.
12. Dli M.I., Vlasova E.A., Sokolov A.M., Morgunova E.V. Creation of a chemical-technological system digital twin using the Python language. Journal of Applied Informatics. 2021. Vol. 16. No. 1 (91). P. 22–31.
13. Бобков В.И. Оптимизация химико-технологического процесса сушки в стационарном режиме многослойной массы фосфоритовых окатышей по критерию энергоресурсоэффективности // Современные наукоёмкие технологии. 2018. № 5. С. 25–29.
14. Kavchenkov V.P., Kavchenkova E.V., Chernenkov I.D. Modeling of the relationship between the earth population growth and the electric energy production processes. Journal of Applied Informatics. 2021. Vol. 16. No. 4 (94). P. 110–121.
15. Matkarimov S.T., Berdiyarov B.T., Yusupkhodaev A.A. Technological parameters of the process of producing metallized iron concentrates from poor raw material. International Journal of Innovative Technology and Exploring Engineering. 2019. No. 8 (11) P. 600–603.

В современной промышленности часто встречаются устройства, работающие с загрузочными зернистыми слоями, как стационарные, так и подвижные [1, 2]. К таким устройствам, в частности, относятся обжиговые конвейерные и агломерационные машины для обогащения рудного сырья в металлургической, химической и горно-обогатительной промышленности. Расчёт таких устройств требует знания аэродинамических характеристик подвижного плотного слоя кускового и окомкованного рудного сырья с перекрёстной подачей газа-теплоносителя [3, 4]. По гранулометрическому составу, как правило, чаще встречаются полифракционные системы, включающие широкие классы по размерам [5, 6]. Геометрическая структура таких систем будет сильно зависеть от степени неоднородности или полидисперсности. Аэродинамические параметры динамического плотного слоя определяются главным образом геометрическими характеристиками полифракционной засыпки рудного сырья и режимом функционирования обжиговых конвейерных машин [7, 8].

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

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

Материалы и методы исследования

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

Bobkov001.wmf,

или среднеповерхностный

Bobkov002.wmf,

где di – средний диаметр i-й фракции, gi – весовая доля i-й фракции. Порозность определялась экспериментально

ε = 1 – ρн / ρм,

где ρн – насыпная плотность, кг/м3, ρм – плотность материала, кг/м3. Удельная поверхность засыпки определялась соотношением:

Bobkov003.wmf,

где s0 – удельная поверхность элементов слоя, м2, D – диаметр рабочего участка, м. Для монофракционных засыпок шаров удельная поверхность определяется как s0 = 6 / d. Для полифракционных засыпок из шаров удельная поверхность слоя определяется соотношением Bobkov004.wmf. Для полифракционных засыпок из частиц неправильной формы определение удельной поверхности основывается на методике использования коэффициента формы зерна:

Bobkov005.wmf, Bobkov006.wmf, Bobkov007.wmf,

где dc – среднеситовой, ds – эквивалентный по поверхности, dv – эквивалентный по объёму диаметр зерна. Эквивалентные диаметры определяются соответственно:

Bobkov1.pdf и Bobkov1.pdf,

где sз, vз – соответственно поверхность и объём зерна. Тогда удельная поверхность зерна полифракционного слоя будет определяться как Bobkov010.wmf.

Особенностью полифракционных засыпок является способность мелких фракций заполнять промежутки между крупными элементами, тем самым резко снижая среднюю порозность слоя. Это явление характеризуется способом укладки и степенью неоднородности засыпки. В различных источниках под степенью неоднородности понимается или отношение диаметра крупной фракции к диаметру мелкой фракции или наоборот. Нами принят коэффициент неоднородности Bobkov011.wmf.

Движение газа в слое подчиняется уравнению Навье – Стокса. Учитывая, что во многих промышленных установках с зернистым слоем скорость фильтрации газовой среды невелика, то можно считать газ несжимаемым. Однако сложность геометрической структуры затрудняет выполнить аналитически граничное условие: отсутствие скольжения на поверхности [9, 10]. В большинстве случаев можно пренебречь гравитационной составляющей уравнения Навье – Стокса.

С принятыми допущениями для усреднённых значений давления и скоростей газа-теплоносителя выражение для перепада давления вдоль потока будет иметь вид

Bobkov012.wmf. (1)

Для анализа уравнения (1) был предложен ряд моделей.

В области преобладания сил вязкости уравнение движения на основе капиллярной модели преобразуется к уравнению Козени – Кармана:

Bobkov013.wmf, (2)

где K – константа Козени – Кармана, определяемая экспериментально для шаров 4÷4.5, s0 – удельная поверхность элементов слоя, Н – высота слоя.

Для слоя шаров:

Bobkov3.pdf.

Тогда коэффициент гидравлического сопротивления слоя:

Bobkov1.pdf,

где v – кинематическая вязкость потока, Reэ – эквивалентное число Рейнольдса.

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

Для модели из ансамбля шаров Хаппелем получено уравнение

Bobkov016.wmf,

которое отличается от (2) поправкой ε3. Введение такой поправки позволяет получить плавный переход к закону Стокса для омывания одиночного шара, что получается при ε → 1.

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

Bobkov1.pdf. (3)

Обычно ограничиваются двумя первыми членами суммы (3), причём показатель степени для первого слагаемого, характеризующего вязкостный режим, получается равным нулю n1 = 0, а n2 = 2, характеризующий второй предельный случай – инерциальный режим. На практике чаще встречаются интерполяционные уравнения вида

Bobkov018.wmf, (4)

а для коэффициента аэродинамического сопротивления:

Bobkov1.pdf,

тогда уравнение (4) преобразуется к виду

Bobkov1.pdf, (5)

где Ки – инерционная составляющая сопротивления.

Эксперименты проводились на лабораторной установке, адаптированной авторами для проведения настоящих исследований, состоящей из рабочего участка высотой 0,5 м с внутренним диаметром для шаров 0,225 м, для кусковых засыпок – 0,3 м. Для выравнивания скорости газа на выходе в участок установлен ресивер с соплом на выходе с коэффициентом неравномерности 0,95. Для подогрева воздуха служит нагреватель, для получения более высоких температур использовалась камера сгорания, работающая на сжиженном газе – смеси пропана и бутана. Расходы воздуха и газа измерялись диафрагмами «четверть круга», температуры и температурные напоры регистрировались потенциометрами. Перепады давления в слое измерялись наклонными микроманометрами. Для усреднения давления по сечению установлены по три импульсные трубки, всего замеры проводились по четырём сечениям. Воздух нагнетался вентилятором высокого давления.

Перед исследованием полифракционных засыпок были проведены опыты с монофракционными слоями шаров с диаметром 16,67; 14,15; 9,51; 7,97 и 6,75 мм. Коэффициенты аэродинамического сопротивления определяли как

Bobkov2.pdf,

где w – скорость потока газа-теплоносителя на свободное сечение, м/с.

После обработки результатов по коэффициенту аэродинамического сопротивления слоя шаров по методу наименьших квадратов усреднённая зависимость получилась:

Bobkov2.pdf,

что удовлетворительно совпадает с известными зависимостями, при этом константа Козени – Кармана получилась равной K = 4, а инерционная составляющая Ки = 0,54 при средней величине порозности ε = 0,36.

При исследовании полифракционных засыпок использовались смеси шаров вышеуказанных диаметров.

Изменение структуры слоя хорошо заметно по зависимости порозности слоя от содержания мелкой фракции и степени неоднородности. Для смеси шаров диаметром 16,67 и 6,35 при содержании мелких шариков 30 % по весу наблюдается минимальное значение порозности 0,295. Такое резкое снижение порозности объясняется соотношением диаметров. При кубической укладке шаров большего диаметра промежуток между шарами составляет Bobkov024.wmf, при диаметре 16,67 мм это будет 6,7 мм, то есть мелкие шарики свободно могут размещаться между крупными, а при степени неоднородности близкой к 0,5 зависимость порозности от весовой доли мелкой фракции получается более плавной, без резкого минимума.

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

При степени неоднородности f менее 0.5:

Bobkov2.pdf.

При степени неоднородности f сравнимой с 0.5:

Bobkov2.pdf.

Чем больше соотношение между диаметрами крупной и самой мелкой фракций, тем сильнее влияние на геометрическую структуру.

Результаты исследования и их обсуждение

Разработана методика расчёта сопротивления слоя при неизотермических условиях на основе позонного разделения слоя. Аэродинамическое сопротивление отдельных зон рассчитывается по уравнению

Bobkov027.wmf, (6)

где hi – высота зоны, м; ρ0i – плотность газа в зоне, т/м3; w0i – скорость газа в зоне, м/с; vti – кинематическая вязкость газа при средней температуре в зоне, м2/с; Ti – средняя температура в зоне, К; k1i, k2i – коэффициенты аэродинамического сопротивления отдельных зон,

Bobkov2.pdf, Bobkov2.pdf,

где c1, c2 – константы в интерполяционной формуле для коэффициента аэродинамического сопротивления Bobkov030.wmf, dэ – эквивалентный диаметр зерна; ψ – коэффициент формы зерна.

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

Сопротивление слоя при высоких температурах во многом определяется свойствами потока, особенно вязкостью. Общеизвестно, что вязкость зависит от состава газового потока и сильно от температуры [13, 14]. Часто зависимость кинематической вязкости от температуры представляют в виде многочлена второй степени по температуре: Bobkov032.wmf.

Предлагается методика расчёта сопротивления неизотермического слоя как условно изотермического по температуре эквивалентной, которая определяется как

Tэкв Bobkov033.wmf. (7)

Температура, рассчитанная по этой формуле, на 5–15 % выше среднеинтегральной, следовательно, не рекомендуется рассчитывать сопротивление по среднеинтегральной температуре.

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

Таблица 1

Гранулометрический состав фосфоросодержащего рудного сырья

di, мм

60

45

35

25

12,5

2,5

ξi

0,0595

0,201

0,293

0,289

0,152

0,005

Прогрев слоя осуществляется продуктами сгорания пропан-бутановой смеси следующего состава:

Bobkov034.wmf, Bobkov035.wmf,

Bobkov036.wmf, Bobkov037.wmf.

Динамическая вязкость смеси рассчитывалась по формуле Хеннинга и Цепперера:

Bobkov2.pdf,

где хi – молярная доля, μi – вязкость, Мi – молекулярная масса, Ткрi – критическая температура i-го компонента.

Изменение динамической вязкости от температуры рассчитывалось по формуле Улыбина:

Bobkov2.pdf,

где

Bobkov040.wmf.

Коэффициенты многочлена Аi, Вi, Сi для различных газов представлены в табл. 2.

Таблица 2

Постоянные аппроксимационного многочлена в формуле динамической вязкости

Газ

A

B

C

Bobkov041.wmf

CO2

0,942

2,1

0,56

115,5

O2

0,76

1,4

0,64

70,2

N2

0,753

1,3

0,4

59,5

H2O

1,082

0,25

0

108

Расчёт аэродинамического сопротивления полифракционного слоя при неизотермических условиях производился по уравнению (5), при этом принимались постоянными константа Козени – Кармана, полученная при температуре T0 = 273 K, и инерционная составляющая, равные соответственно К = 7 и К = 1,2.

Сопротивление рассчитывалось численным интегрированием по высоте слоя:

Bobkov042.wmf. (8)

Использование условно изотермического слоя по эквивалентной температуре, рассчитанной по (7) даже для слабонеизотермических распределений, даёт значительное отклонение от уравнения (8). Эквивалентные температуры значительно превышают среднеинтегральные и в некоторых случаях превышают максимальную температуру в слое, что лишено физического смысла [15].

Следует рассчитывать сопротивление в неизотермическом полифракционном слое по уравнению (8).

Заключение

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

1. Структура слоя из смеси шаров зависит от степени неоднородности. Минимальная порозность слоя соответствует весовой доле мелкой фракции 30 %.

2. Получены значения константы Козени – Кармана и инерционной составляющей для двух значений коэффициента неоднородности.

3. При неизотермических условиях расчёт сопротивления слоя необходимо проводить с учётом изменения свойств потока газа-теплоносителя. Расчёт по среднеинтегральным температурам условно изотермического слоя приводит к значительным погрешностям.

Исследование выполнено за счет гранта Российского научного фонда № 22-11-00335, https://rscf.ru/project/22-11-00335/.


Библиографическая ссылка

Бобков В.И., Орехов В.А. ИССЛЕДОВАНИЕ АЭРОДИНАМИКИ ПОЛИФРАКЦИОННОГО ПЛОТНОГО СЛОЯ РУДНОГО СЫРЬЯ // Успехи современного естествознания. – 2022. – № 9. – С. 67-72;
URL: https://natural-sciences.ru/ru/article/view?id=37894 (дата обращения: 18.06.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674