VEROYATNOSTNO-STATISTICHESKIY PODHOD K OPISANIYU AERODINAMICHESKOGO VZAIMODEYSTVIYA KOLLEKTIVA PADAYUSCHIH CHASTIC S VOZDUHOM
Abstract and keywords
Abstract (English):
In the mountain industry a considerable part of technological processes are followed by loadings of bulk in various capacities and the bunker. When loading bunkers of silage type self-dumping carts the problem of knocking-out of dust from loading apertures is particularly acute especially. Work is devoted to improvement of methods of calculation of aspiration when loadings bunkers of silage type by polydisperse material. New statistical approach to the accounting of conditions of constraint on coefficient of front resistance of particles in the conditions of a free stream of the falling material is offered. When falling a free stream of material of a particle of larger fractions shade small and dust particles which aerodynamic resistance in such conditions isn´t considerable. Analytical conclusion of probability of active impact on air of the particles which are out of an aerodynamic shadow that allows to receive a technique of assessment of ezhektiruyushchy ability of a stream of particles at dust removal of bunkers of silage type has been offered. At further development of this approach a conclusion of probability of aerodynamic interaction of particles of polydisperse material is possible. The accounting of dispersion of the overloaded material will allow to calculate optimum amount of aspiration air and to reduce energy consumption of systems of aspiration.

Keywords:
aspiration of bunkers, loading of bunkers, unequigranular material, ejection of air, loudspeaker of particles, aerodynamics, dedusting ventilation
Text

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

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

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

В случае значительной площади загрузочных проемов, что характерно для бункеров силосного типа как правило руководствуются защитной скоростью 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.

References

1. Butakov S. E. Aerodinamika sistem promyshlennoy ventilyatsii. M.: Profizdat, 1949. 268 s.

2. Neykov O. D., Logachev I. N. Aspiratsiya i obespylivanie vozdukha pri proizvodstve poroshkov. M.: Metallurgiya, 1981. 192 s.

3. Obespylivayushchaya ventilyatsiya: monografiya / V.A. Minko, I.N. Logachev, K.I. Logachev i dr.; pod obshch. red. V.A. Minko. Belgorod: izd-vo BGTU, 2010. 565 s.

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. R. 103-107.

5. Popov E.N., Seminenko A.S. Opredele-nie ob´´emov aspiratsionnogo vozdukha pri obespylivanii bunkerov. Nauka i molodezh´ v nachale novogo stoletiya: Materialy III Mezhdunarodnoy nauch.-prakt. konf. Gubkin: IP Uvarov V.M. 2010. S. 71-75.

6. Gol´tsov A.B., Kireev V.M., Popov E.N. Problemy kompleksnogo obespylivaniya pri pererabotke rudnykh materialov. Sbornik trudov № 4. Voronezh: izd-vo VGASU, 2007. S. 123-129.

7. Gervas´ev A. M., Olifer V. D. Nekotorye rezul´taty issledovaniya protsessov, proiskhodyashchikh pri peregruzke sypuchikh materialov po vertikal´nym zhelobam. - V ob. Obespylivayushchaya ventilyatsiya / VNIOT g. Sverdlovsk. 1973. 3-9 s.

8. Liu Ze Qin, 2003. Air entrainment in free falling bulk materials, Doctor of Philosophy thesis, Faculty of Engineering, University of Wollongong.

9. Idel´chik I. E.. Spravochnik po gidravlicheskim soprotivleniyam. M.: Kniga po Trebovaniyu, 2012. 466 s.

10. Pulachev V.S. Teoriya veroyatnosti k matematicheskoy statistike. Moskva: Nauka, 1979. S. 496


Login or Create
* Forgot password?