Белгород, Белгородская область, Россия
В горноперерабатывающей промышленности значительная часть технологических процессов сопровождаются загрузками сыпучего материала в различные емкости и бункера. При загрузке бункеров силосного типа саморазгружающимися тележками проблема выбивания пыли из загрузочных проемов стоит особенно остро. Работа посвящена совершенствованию методов расчета аспирации при загрузках бункеров силосного типа полидисперсным материалом. Предложен новый статистический подход к учету условий стесненности на коэффициент лобового сопротивления частиц в условиях свободной струи падающего материала. При падении свободной струи материала частицы более крупных фракций затеняют мелкие и пылевые частицы, аэродинамическое сопротивление которых в таких условиях не значительно. Был предложен аналитический вывод вероятности активного воздействия на воздух частиц, находящихся вне аэродинамической тени, что позволяет полу-чить методику оценки эжектирующей способности потока частиц при обеспыливании бункеров силосного типа. При дальнейшем развитии данного подхода возможен вывод вероятности аэродинамического взаимодействия частиц полидисперсного материала. Учет дисперсности перегружаемого материала позволит рассчитать оптимальное количество аспирационного воздуха и снизить энергопотребление систем аспирации.
аспирация бункеров, загрузка бункеров, полидисперсный материал, эжектирование воздуха, динамика частиц, аэродинамика, обеспыливающая вентиляция
Введение. Значительная часть технологических процессов переработки горных пород сопровождаются загрузками сыпучего материала в различные емкости и бункера.
Расход аспирационного воздуха, удаляемого от местных отсосов является основным параметром, определяющим энергоемкость и стоимость эксплуатации вытяжной вентиляции.
Существующие методики определения расхода эжектируемого воздуха, предполагают замену реального полидисперсного материала неким монофракционным материалом, имеющим частицы некторого среднего диаметра. Большая часть обрабатываемых материалов являются полидисперсными, что приводит к необходимости введения опытных коэффициентов.
В случае значительной площади загрузочных проемов, что характерно для бункеров силосного типа как правило руководствуются защитной скоростью 0,5–1 м/с, обеспечивающей невыбивание запыленного воздуха из бункера. Такой подход приводит к значительным расходам аспирационного воздуха, не учитывающих как особенностей самого перегружаемого материала, так и аэродинамических процессов протекающих внутри бункера.
В связи с этим уточнение существующих методик расчета и проектирования систем аспирации на предприятиях горноперерабатывающей отрасли является весьма актуальной задачей.
Фундаментальные основы теории эжекции воздуха равноускоренным потоком падающих частиц были заложены С.Е. Бутаковым [1] и развиты О.Д. Нейковым [2]. В последние годы проблемой снижения объемов аспирируемого воздуха в России занимались В.А. Минко, И.Н. Логачев [3], К.И. Логачев [4], изучившие закономерности движения сыпучих материалов по желобам. Так же проблемами аспирации бункеров занимались Семинеко А.С., Гольцов А.Б. [5, 6]. За рубежом известны работы Олифера В.Д. [7], ZeQinLiu [8]. Однако в трудах этих ученых рассматриваются главным образом монофракционные потоки в то время как большинство перерабатываемых материалов полифракционые.
Основная часть. Аэродинамическое взаимодействие коллектива частиц в струе материала и эжектируемого воздуха, как известно [4], определяется суммой аэродинамических сил всех частиц этого коллектива в единице объема :
где – аэродинамическая сила одной частицы, Н; – плотность воздуха, кг/м3; – счетная концентрация частиц, 1/м3; – элементарный объем двухкомпонентной среды «твердые частицы – воздух», м3; – площадь миделевого сечения частицы, м2; – скорость эжектируемого воздуха, м/с; – коэффициент лобового сопротивления частицы, зависящий в общем случае, от геометрической формы, числа Рейнольдса относительного движения и объемной концентрации этих частиц.
В отличие от достаточно полно изученного механизма межкомпонентного взаимодействия в потоке усредненных по крупности частиц мы будем рассматривать аэродинамическое взаимодействие потока частиц с разными эквивалентными диаметрами , т.е. поток полифракционного сыпучего материала, состоящего из N фракций.
Будем рассматривать ускоренный поток свободно падающих частиц, объемная концентрация которых уменьшается по мере перемещения материала от "насыпной" (во время транспортировки конвейером или в момент сбрасывания)
до "разреженной" (при падении)
где – массовый расход падающих частиц, кг/с; – площадь поперечного сечения струи частиц, м2; – скорость падения частицы, м/с; , – соответственно, плотности материала частиц и их насыпная плотность, кг/м3.
Чтобы не вводить эмпирические поправки на коэффициент *, связанные со стесненностью частиц в процессе их падения используем вероятностный подход. Так как обычно эти поправки связаны с объемной концентрацией потока частиц условной монофракции, а не реального потока частиц разной крупности. Передоложим, что коэффициент лобового сопротивления частиц i-й фракции пропорционален вероятности активного аэродинамического взаимодействия этой частицы :
где – коэффициент пропорциональности (в дальнейшем коэффициент пропорциональности принимаем равным единице); – вероятность нахождения частицы в аэродинамической тени; – коэффициент лобового сопротивления одиночной (свободной, не затененной) частицы i-й фракции.
Таким образом под вероятностью будем понимать отношение количества частиц, находящихся вне аэродинамической тени к общему количеству частиц в рассматриваемом объеме двухкомпонентной смеси. Ее величина в процессе падения частиц будет возрастать практически от нуля (в момент сбрасывания, когда все частицы находятся в аэродинамической тени соседних частиц) до величины, близкой к единице (при весьма разреженном коллективе частиц, падающих с большой скоростью).
Рассмотрим случай падения коллектива одинаковых по массе частиц.
Полагаем, что место положения частицы в элементарном объеме равновероятно и объемная концентрация этих частиц равна:
где – массовый расход частиц, кг/с; – поперечное сечение струи, м2; – плотность частиц, кг/м3; – скорость частиц, м/с.
В силу этого вероятность нахождения частицы в аэродинамической тени других частиц в этом элементарном объеме составляет:
где – счетная концентрация частиц в элементарном объеме, шт/м3, очевидно равна
– объем одной частицы, м3, – объем аэродинамической тени одной частицы, м3
Тогда с учетом (6) имеем:
Вероятность активного воздействия частиц на воздух, которую можно рассматривать как противоположное событие нахождению частиц в аэродинамической тени, равна:
где
Объем аэродинамической тени зависит от формы частиц и числа Рейнольдса. При Re > 105 по данным академика Скочинского А.А., приведенных в справочнике И.Е. Идельчика [9]: а) для случая поперечного обтекания бруса сечением a*b (вектор скорости воздуха направлен вдоль длинной стороны a): при a/b =2,2 ([9], стр. 407); б) для случая поперечного обтекания цилиндра диаметром d:; в) для случая поперечного обтекания диска диаметром d и толщиной (вектор скорости воздуха направлен перпендикулярно плоскости диска):, ([9], стр. 4067) (при δ = 0,1 d), ([9], стр. 397).
В данном случае в качестве длины аэродинамической тени принималось наименьшее расстояние между двумя одинаковыми телами, когда суммарный коэффициент лобового сопротивления становится равным удвоенному коэффициенту лобового сопротивления одиночного тела. Если принять более жесткие условия для аэродинамической тени, например, если принять, что при нахождении частицы в "плотной" тени суммарный коэффициент равен коэффициенту сопротивления одиночной частицы, то отношение – для случая "а"; – для случая "б"; – для случая "в".
Примем средневзвешенное значение для рассмотренных тел при двукратном значении суммы коэффициентов местного сопротивления в качестве реперного расчета .
Очевидно, вероятность нахождения в аэродинамической тени двух и более частиц (n > 2) будет равна произведению вероятности
или т.к. β1
, (11)вероятность нахождения частицы в "плотной" тени второй частицы, скажем для случая "б" (при), составляет. Можем этот случай считать маловероятным, поскольку вероятность его по сравнению со случаем нахождения в тени первой частицы по крайней мере на порядок больше.
Однако в случае потока полифракционных частиц этой вероятностью пренебрегать нельзя. Дальнейшая работа должна быть направлена на определение вероятностей активного аэродинамического взаимодействия для полидисперсного материала.
Заключение. Поток падающего материала условно разделяется на аэродинамически активные частицы и частицы, находящиеся в аэродинамической тени и не участвующие в силах аэродинамического взаимодействия при эжектировании воздуха в свободных струях падающего материала.
Выводы. Произведен аналитический вывод вероятности активного воздействия частиц на воздух в свободной струе падающего материала, позволяющий не прибегая к эмпирическим соотношениям определить эжектирующую способность этой струи. Усовершенствованная методика расчета расхода аспирационного воздуха при загрузках бункеров силосного типа, позволит снизить энергоемкость систем аспирации.
*Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект # 14-41-08005 р офи м), а так же в рамках НИР «Разработка методик расчета систем обеспыливания и исследование условий загрузки бункеров с учетом дисперсности материалов» НИР: Б8/13.
1. Бутаков С. Е. Аэродинамика систем промышленной вентиляции. М.: Профиздат, 1949. 268 с.
2. Нейков О. Д., Логачев И. Н. Аспирация и обеспыливание воздуха при производстве порошков. М.: Металлургия, 1981. 192 с.
3. Обеспыливающая вентиляция: монография / В.А. Минко, И.Н. Логачев, К.И. Логачев и др.; под общ. ред. В.А. Минко. Белгород: изд-во БГТУ, 2010. 565 с.
4. Logachev K.I., Averkova O.A., Tolma-cheva E.I., Logachev A.K., and Dmitrienko V.G., Modeling of Air and Dust Flows in the Range of Action of a Round Suction Funnel Above an Impermeable Plane. Refractories and Industrial Ceramics. 2016. 57. Р. 103-107.
5. Попов Е.Н., Семиненко А.С. Определе-ние объемов аспирационного воздуха при обеспыливании бункеров // Наука и молодежь в начале нового столетия: Материалы III Международной науч.-практ. конф. Губкин: ИП Уваров В.М. 2010. С. 71-75.
6. Гольцов А.Б., Киреев В.М., Попов Е.Н. Проблемы комплексного обеспыливания при переработке рудных материалов // Сборник трудов № 4. Воронеж: изд-во ВГАСУ, 2007. С. 123-129.
7. Гервасьев А. М., Олифер В. Д. Некоторые результаты исследования процессов, происходящих при перегрузке сыпучих материалов по вертикальным желобам. - В об. Обеспыливающая вентиляция / ВНИОТ г. Свердловск. 1973. 3-9 с.
8. Liu Ze Qin, 2003. Air entrainment in free falling bulk materials, Doctor of Philosophy thesis, Faculty of Engineering, University of Wollongong.
9. Идельчик И. Е.. Справочник по гидравлическим сопротивлениям. М.: Книга по Требованию, 2012. 466 с.
10. Пулачев В.С. Теория вероятности к математической статистике. Москва: Наука, 1979. С. 496
