Ekaterinburg, Russian Federation
Ekaterinburg, Russian Federation
Ekaterinburg, Russian Federation
Ekaterinburg, Russian Federation
GRNTI 55.01 Общие вопросы машиностроения
BBK 345 Общая технология машиностроения. Обработка металлов
The scheme of the experimental device for the determining of the velocity of particles in pneumatic transport conditions was given. The method of magnetic marking was used. The pneumatic transport system was implemented on the basis of a pneumatic chamber pump. The particle velocity was measured on a straight section with steady motion. The mode parameters of the two-phase flow were varied: air flow, concentration, pressure. The experiments were carried out with various bulk materials. Comparison of the experimental and calculated values of the relative particle velocity was given in respect to various known dependences.
pneumatic transport, two-phase flow, air flow, relative particle velocity, indicator particle
Как известно, аэродинамическое сопротивление пневмотранспортного тракта складываются из потерь давления на разгон частиц, потерь давления по длине, потерь давления на подъем материала и местных потерь давления. При этом потери давления по длине составляют основную долю. Однако чтобы определить потери давления на длинных участках с установившимся движением необходимо знать относительную скорость движения частиц твердой фазы [1, 2, 3]. Поэтому, вопросу определения относительной скорости частиц посвящено большое количество работ и предложено достаточно большое число различных зависимостей. Так в работах [4, 5] приводится анализ шести зависимостей, предложенных разными авторами. В тоже время работ, посвященных экспериментальному определению скорости частиц в условия пневмотранспорта очень мало. Цели настоящей работы – экспериментальное определение относительной скорости и выбор зависимости для ее расчета.
Схема установки для определения скорости движения частиц в условиях пневмотранспорта приведена на рис.1. и представляет собой систему высоконапорного пневмотранспорта на базе лабораторного пневмокамерного насоса (ПКН) емкостью 30 литров. Транспортный трубопровод имеет внутренний диаметр 16 мм и длину 25 м. Установка включает в себя: 1 – подвод сжатого воздуха из коллектора; 2 – ПКН; 3 – клапан загрузки материала; 4 – пневмотранспортный тракт; 5 – циклон с накопительным бункером, для отделения материала от воздуха; 6 – фильтр санитарной отчистки; 7 – регулировочные вентили расхода воздуха; 8 – стандартная диафрагма, с импульсными трубками; 9,10,11,12 – импульсные трубки для измерения давления; 13 – измерительный комплекс включающий в себя: аналогово-цифровой преобразователь, пять датчиков давления серии «Метран 100», преобразователь интерфейса; 14 – устройство ввода индикаторной частицы; 15 – индукционные датчики; 16 – аналогово-цифровой преобразователь электромагнитных импульсов с частотой дискредитации 44 кГц; 17 – компьютер для обработки данных; 18 – эжектор.
Длина прямого, горизонтального участка перед индукционными датчиками составляет 437 калибров, это позволяет материалу разогнаться и выйти на установившийся режим движения. Во время установившегося режима движения двухфазного потока вводится индикаторная частица, представляющая собой пенопластовый шарик диаметром 10–12 мм, внутри которого помещена малая магнитная частица. Как показали визуальные наблюдения на стеклянной трубке, такой шарик в условиях пневмотранспорта движется с осредненной скоростью твердой фазы и не может обогнать поток, или отстать от него. Прохождение частицы регистрируется индукционными датчиками, расположенными друг от друга на расстоянии l=160 мм. Сигнал от индукционных датчиков фиксируется быстродействующим аналогово-цифровым преобразователем и записывается в файл на компьютере. На рис. 2. представлен временной график сигналов от индукционных датчиков при пролете индикаторной частицы.
Рис. 1. Схема экспериментальной установки для определения относительной скорости частиц
∆τ
Датчик №2
Датчик №1
Рис. 2. Треки прохождения магнитной частицы-индикатора через сечение катушек
Измерив, промежуток времени ∆τ между импульсами, можно вычислить скорость движения индикаторной частицы U, м/с по очевидной зависимости:
(1)
Опыты проводились на двух типах сыпучего материала – монодисперсном и полидисперсном. В качестве монодисперсного материала использовался пропант представляющий собой шарообразные частицы узкой фракции со средневзвешанным диаметром 1243,7 мкм и насыпной плотностью 1940 кг/м3. Кажущаяся плотность – 2750 кг/м3. Полидисперсный полевой шпат имеет размеры частиц в диапазоне 0-400 мкм, со средневзвешенным диаметром 492,5 мкм. Насыпная плотность равна 1440 кг/м3, кажущаяся плотность 2398 кг/м3 (см. табл. 1).
Таблица 1
Гранулометрический состав полевого шпата
|
Сита, мкм |
R, гр |
r, гр |
|
1600 |
0,3 |
0,3 |
|
1000 |
3,1 |
2,8 |
|
630 |
21,1 |
18 |
|
400 |
57,8 |
36,7 |
|
200 |
92,7 |
34,9 |
|
160 |
96 |
3,3 |
|
100 |
98,8 |
2,8 |
|
63 |
99,5 |
0,7 |
|
0 |
100 |
0,5 |
Измерив, массовый расход воздуха qм, кг/с, по стандартной диафрагме, можно вычислить скорость воздушного потока перед индукционными датчиками Wид, м/с, по известной зависимости:
, (2)
где, FТ – площадь сечения трубы, м2; плотность воздуха ρид, кг/м3, определяем по уравнению состояния:
, (3)
где, Pбар – барометрическое давление, Па; Pид – давление в районе установки индикаторных датчиков, Па; R – удельная газовая постоянная для воздуха, Дж/(К·кг), Т – абсолютная температура воздуха, К.
Экспериментальное значение относительной скорости jэкс определяется по зависимости:
. (4)
На рис.3. изображен график изменения скорости в сечениях на трех участках транспотного трубопровода. Вертикальная линия на графике обозначает момент пролета индикаторной частицы через сечение индукционных датчиков.
Рис. 3. Изменение скоростей воздушного потока в транспортной трубе на различных участках:
Wид – скорость перед индукционными датчиками; Wн – скорость в начале тракта; Wк – скорость в конце тракта
При вычислении расчетных значений относительной скорости частиц использовались следующие формулы.
Зависимость (5) Овена, Разумова [2, 3], где vos – скорость осаждения частиц, м/с; w – средняя по сечению скорость потока, м/с.
. (5)
Зависимость (6) Басова [2], где d – диаметр частицы, м.
. (6)
Зависимость (7) Калугина [6], где ζ∂ – коэффициент сопротивления потока; Сx – коэффициент сопротивления частицы; μ – расходная концентрация твердой фазы, кг/кг; ρt, ρ – соответственно плотность частиц и плотность среды, кг/м3; D – диаметр трубы, м.
.(7)
Зависимость (8) Шухарта [2], где С – константа, С=0,014 – 0,05; Fr – критерий Фруда; Fr0 – критерий Фруда, рассчитанный по скорости потока соответствующей осаждению частиц.
. (8)
Зависимость (9) Гастерштадта [7].
. (9)
Формула (10) зональной модели [4], где S1, S2 – отношение площадей зон к площади сечения трубы.
. (10)
В табл. 2 приведены экспериментальные и расчетные данные.
Таблица 2
Соотношение экспериментальных и расчетных данных
|
Режимные параметры |
Монодисперсный (пропант) |
Полидисперсный (шпат) |
|||||||
|
Давление Pc, кПа |
0,452 |
0,342 |
0,353 |
0,370 |
0,344 |
0,437 |
0,287 |
0,447 |
0,456 |
|
Концентрация μ, кг/кг |
27,25 |
22,49 |
23,27 |
21,74 |
22,32 |
24,51 |
16,58 |
24,06 |
30,15 |
|
Скорость Wи, м/с |
20,41 |
16,34 |
15,93 |
15.93 |
16,57 |
11,25 |
19,34 |
18,75 |
21,34 |
|
Экспериментальное значение скорости частиц |
|||||||||
|
Скорость U, м/с |
13,14 |
10,6 |
9,94 |
9,94 |
11,36 |
10,6 |
16,74 |
12,23 |
12,23 |
|
Экспериментальное значение относительной скорости |
|||||||||
|
Скорость jэкс |
0,644 |
0,649 |
0,624 |
0,624 |
0,686 |
0,942 |
0,865 |
0,652 |
0,573 |
|
Расчетные значения относительной скорости по различным зависимостям |
|||||||||
|
Зональная модель |
0,933 |
0,944 |
0,942 |
0,943 |
0,944 |
0,983 |
0,984 |
0,984 |
0,984 |
|
Калугин |
0,988 |
0,988 |
0,988 |
0,988 |
0,988 |
0,992 |
0,992 |
0,992 |
0,992 |
|
Басов |
0,715 |
0,740 |
0,729 |
0,732 |
0,744 |
0,985 |
0,986 |
0,989 |
0,991 |
|
Овен, Разумов |
0,687 |
0,715 |
0,702 |
0,706 |
0,719 |
0,984 |
0,985 |
0,987 |
0,990 |
|
Шухарт |
0,591 |
0,687 |
0,686 |
0,690 |
0,691 |
0,698 |
0,677 |
0,625 |
0,590 |
|
Гастерштадт |
0,871 |
0,895 |
0,892 |
0,894 |
0,896 |
0,945 |
0,941 |
0,931 |
0,924 |
Исходя из представленных в табл.1 данных следует, что наиболее близкие значения к экспериментальным данным дают зависимости для расчета относительной скорости (8) Шухарта и (5) Овена и Разумова. Эти зависимости можно использовать как для полидисперсных, так и для монодисперсных составов. Для монодисперсных материалов, также можно использовать зависимость (6) Басова [2].
1. Shishkin S.F., Shishkin A.S. Dvizhenie dvuhfaznogo potoka v trube postoyannogo secheniya / Fiziko-himiya i tehnologiya oksidno-silikatnyh materialov: materialy mezhdunar. nauch.-tehn. konf. Ekaterinburg, 17-19 fevr. 2000 g. // Ural. gos. tehn. un-t, Ural. otd-nie Assoc. stroit. vuzov RF. Ekaterinburg: GOU VPO UGTU-UPI, 2000. C. 251.
2. Razumov I. M. Pnevmo- i gidrotransport v himicheskoy promyshlennosti. M.: Himiya, 1979. 248s.
3. Sukomel A.S., Cvetkov F.F., Kerimov R.R. Teploobmen i gidravlicheskie soprotivleniya pri dvizhenii gazovzvesi v trubah. M.: Energiya, 1977. 192 s.
4. Gavrilyuk D.N., Shishkin S.F. Otnosi-tel'naya skorost' dvizheniya pnevmotransporta na stacionarnyh uchastkah // Vestnik BGTU im. V. G Shuhova. 2009. №3. S. 122-126.
5. Busroyd R. Techenie gaza so vzveshennymi chasticami. M.: Mir, 1975. 378 s.
6. Kalugin B.F. Poteri napora ot udarov chastic o stenki pri pnevmaticheskom transporte po gorizontal'nym trubam // Inzhenerno-fizicheskiy zhurnal. 1960. T. IV. №7. S. 40-46.
7. Mednikov E. P. Turbulentnyy perenos i osazhdenie aerozoley. M.: Nauka, 1980. 176 s.
