MECHANICS OF THE SUSPENDED LAYER WITH POLYDISPERSE PARTICLES IN CONTINUOUS AIR MIXERS
Abstract and keywords
Abstract (English):
The quality of mixing of dry building mixes largely depends on the specialized equipment that is used in this process. For the best mixing of dry building mixes, continuous pneumatic mixers are used. Continuous air mixers are the most effective equipment for producing dry building mixes with a higher degree of uniformity. This indicator allows to achieve a better final product. Currently, there are a large number of continuous air mixers that are actively being modernized both in our country and abroad. The article presents a technique that describes the mechanics of a suspended layer with polydisperse particles in the working volume of the displacement palette of a continuous air mixer. The method is based on research conducted by the authors of the article, taking into account experimental data obtained by other scientists. This paper summarizes an approach to numerical modeling of a suspended layer with particles of different dispersion (for example, for dry building mixes) and density. The dependences of the layer boundary, its porosity, and speed parameters on the particle diameters, geometric parameters of the mixing chamber, and other technological values that affect the suspended layer inside the chamber are established. The article presents an algorithm for calculating the mechanics of the weighted layer, which proved itself in numerical modeling of the layer with high accuracy of calculations, which is later confirmed by the necessary amount of experimental data.

Keywords:
mixing chamber, flow, particle, braking parameter, circulation, tangential velocity, layer porosity
Text
Publication text (PDF): Read Download

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

Коллективом авторов разработаны конструкции пневмосмесителей непрерывного действия [1, 2, 3], которые отвечают современным требованиям, предъявляемым к технологическому оборудованию для гомогенизации полидисперсных компонентов: низкая металлоемкость конструкции (за счет изготовления корпуса смесительного агрегата из мягкого полиуретана); отсутствие в готовом продукте металлических включений; высокая производительность и коэффициент однородности смеси.

За последние годы как у нас в стране [4, 5] так и за рубежом [6, 7] эффективно используются агрегаты с пневматическим принципом перемешивания порошкообразных масс: вертикального исполнения, горизонтального типа, с псевдоожижением слоя компонентов смеси и др. В этом особое место занимают пневмосмесители горизонтального типа непрерывного действия с камерой смешения переменного поперечного сечения. Для использования пневмосмесителей при производстве сухих строительных смесей различной дисперсности, в конструкции которых присутствует камера в виде параболического корпуса, очень важно уделить внимание динамике полидисперсного потока [8]. В виду сложности всех динамических моделей на различных этапах движения компонентов в пневмосмесителе необходимо поддерживать полидисперсный поток во взвешенном состоянии. Именно поэтому очень важно определить характер поведения частиц различного диаметра в рабочем объеме пневмосмесителя. Для этого необходимо описать механику взвешенного слоя в камере смешения с целью управления эффективными режимами эксплуатации пневмосмесителя [9].

Основная часть. Условия механики взвешенного слоя сформулируем так: если критический диаметр частицы dкр больше максимального dmax  для имеющихся в слое частиц, то все частицы будут взвешены [10]. В этом случае аргумент функции распределения частиц τкр≫1  и D(τкр)≈1  поэтому  из (13) параметр S=1 , т.е. полностью взвешенный слой не будет тормозить поток о криволинейную поверхность камеры. С уменьшением такого торможения возрастает значение торможения потока с частицами о торцевые днища камеры. Этот вид торможения играет более существенную роль для мелких частиц, взвешенных в газообразном объеме, в камере с параболическими коническими стенками [11]. В соответствии с [12] для полностью взвешенного слоя будем учитывать торможение потока о торцевые поверхности. Рассматривается однородный взвешенный слой цилиндрической формы длиной L , наружным радиусом Rн  и внутренним Ri [13]. Если M  - масса частиц в слое, то доля объема, занятая средой, т.е. порозность слоя ε , будет:

ε=VcV=1-VV=1-Mπ(Rн2-Ri2)Lρp ,   (1)

где V  – объем слоя; Vc  – объем занимаемый средой; V  – объем всех частиц. Тогда среднюю плотность слоя можно записать так:

ρс=ρVc+ρpVpV=ρp1-ε+ρε             (2)

Рассмотрим взаимодействие кольцевого взвешенного слоя радиусом  r  и шириной δr  с торцевыми  стенками камеры смешения. Сила воздействия двухфазного потока на единицу площади торцевой поверхности определяется скоростным напором Cfρcvc2  , где Cf  – коэффициент трения потока о стенку [14]. Тогда момент сил взаимодействия кольцевого слоя на две торцевые стенки камеры равен:

=2Cfρcvc222πrrdr=2πCfρcГ2 dr ,   (3)

где Г= vcr  – циркуляция потока, проходящего через слой; vc  – тангенциальная скорость среды и частиц слоя [15].

Этот момент сил приводит к уменьшению потока момента количества движения средыΩ=GГ , проходящей через слой dr :

=GdГ                          (4)

Исключая  из (3) и (4), получаем уравнение

GdГ=2πCfГ2ρсdr ,

после интегрирования которого при граничном условииГRн= Гск на наружном радиусе Rн  слоя имеем [16]:

Г=Гск1+2πCfρcRн-rГскG.                 (5)

Так как наружный радиус слоя Rн  может быть меньше Rк , то в областиRкrRн  циркуляция будет постоянна [17]. Поэтому можно записать

Гск=vскRн ,                          (6)

Отсюда

Гск=4πCfρcГкLG+1-12πCfρcLG                   (7)

где Гк=vкRк;vк  – тангенциальная скорость на периферии ненагруженной камеры смещения.

Выведенные соотношения зависят от коэффициента трения Cf . В [11] результате сопоставления расчетов с экспериментом полученоCf=0,003 , а [8] для псевдоожиженного слоя по аналогии с дисперсно-кольцевым потоком в трубах предлагается Cf=0,005 . Такого порядка значения для коэффициента сопротивления следуют из формулы Прандтля [8]

Cf=0,077Re0,2                                 (8)

при 5∙105<Re<107 , которая обобщает эксперименты по сопротивлению гладкой пластины. Для камеры смешения Re=vкRкv .

Методика исследований. Последовательность действий для расчёта взвешенного слоя в камере смешения пневмосмесителя:

1. Исходными данными является геометрия камеры смешения, массовый расход полидисперсной воздушной среды, её свойства, свойства частиц и масса слоя М.

2. Определяется максимальный и минимальный диаметр частиц. При известном нормально-логарифмическом распределении эти диаметры рассчитываются с 5%-ной обеспеченностью. Тогда согласно [18]:

τ=lgd,d50         lgσ                      (9)

где τ  – аргумент функции распределения ; σ-дисперсия; d-диаметр частицы; d50- медианный диаметр частицы; и с учётом того, что D1.65=1-D-1.65=0.95 , получаем

dmax=d50σ1.65              (10)
dmin=d50σ-1.65             (11)

3. По экспериментальной зависимости в [19]

1SM1-p1-Dτкрff+1            (12)

где s – параметр торможения потока зондом или слоем частиц; М – средняя масса слоя в смесительной камере; Р – порозность слоя; D – функция распределения от τкр(τ) ; ff-  коэффициент трения частиц о поверхность камеры;

определяем параметр торможения S. Если S 1, это свидетельствует, что частицы не соприкасаются с поверхностью камеры.

 

Таблица 1

Входные данные

Наименование параметра

Обозначение

Величина

Массовый расход газа

G

.24000E-01(кг/с)

Радиус камеры

Rk

.80000Е-01(м)

Радиус выходного отверстия

Ri

.27500Е-01(м)

Длина камеры

L

.20000Е+00(м)

Площадь щелей завихрителя

Fvx

.36000Е-03(м**2/с)

Угол наклона

Psivx

.60000Е+02(град)

Вид завихрителя

VID

1

Кинематическая вязкость

ν

.15000Е-04(м**2/с)

Плотность среды

ρc

.12600Е+01(кг/м**3)

Плотность частиц

ρp

.23200Е+04(кг/м**3)

Медианный диаметр частиц

d50

.10000Е-04(м)

Мин. диаметр частиц

dmin

.96272Е-05(м)

Макс. диаметр частиц

dmax

.10387Е-04(м)

Логарифм дисперсии распр.

σ

.10000Е-01

Коэф. Отставания частицы

β

.10000Е+01

Масса слоя

M

.14600Е+00(кг)

Структура слоя для частиц с dmin  
(1-одна ветвь, 2-две ветви)

KOD3

1

Относительная точность

EPS

.10000Е-01

 

Таблица 2

Результаты расчетов

Наименование параметра

Обозначение

Величина

Скорость у стенки камеры

vск

.67976Е+01(м/с)

Средняяпорозность слоя

ε

.97354Е+00

Возможные параметры слоя при  А=.45926Е+01

Диаметр частицы на R1

dR1

.95721Е-05(м)

Тангенциальная скорость

v1

.49266Е+01(м/с)

Минимальная танг. скорость

vmin

.39925Е+01(м/с)

Радиус vmin

Rmin

.48710Е-01(м)

Диаметр частицы наRmin

dRmin

.11812Е-04(м)

Диаметр частицы наRк

dRк

.69374Е-05(м)

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

Диаметр частицы

dCmin

.96272Е-05(м)

Внутренний радиус слоя

R1

.27762Е-01(м)

Тангенциальная скорость

vCR1

.48984Е+01(м/с)

Параметры для средней частицы D50

Диаметр D50

d50

.10000Е-04(м)

На внутренней ветви

Радиус траектории

R50

.29633Е-01(м)

Тангенциальная скорость

v50

.47158Е+01(м/с)

На наружной ветви

Радиус траектории

R50

.67786Е-01(м)

Тангенциальная скорость

v50

.47158Е+01(м/с)

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

Диаметр

dcRн

.10387Е-04(м)

Наружный радиус слоя

Rн

.65626Е-01(м)

Тангенциальная скорость

vcRn

.45401Е+01(м/с)

 

4. Задаются границы слоя Rн=R1  и Ri=R1  (где Rн -начальный радиус камеры, R1 -радиус выходного отверстия)

5. По (1) и (2) рассчитывается плотность слоя pc

6. По (10) определяется коэффициент трения потока о поверхность Cf .

7. Необходимая тангенциальная скорость частиц в слое на определённом радиусе Ri  находится из условия Ks=Ri  (где Ks  – коэффициент сепарации), из которого из стоксовских частиц получаем:

vc=9vQρπL(1-ρ)d2                    (13)

здесь v – климатическая вязкость; L – длина камеры; D – диаметр частицы [20].

Как видно, скорость зависит от диаметра частицы и явно не зависит от радиуса Ri

8. Согласно (9) определяется циркуляция Гск  на периферии слоя

9. Из (5) определяется радиус Ri=RiRк  на котором частицы имеют скоростьvc :

Ri=(1+ARн(1+АRн)2-4Аvc2A              (14)


где

A=2πCfpcRнτск              G                (15)

vc=vcvск                              (16)
 

Здесь vск  – тангенсальная скорость слоя на поверхности камеры; vc- среднее значение тангенсальной скорости частиц в камере; А – геометрический коэффициент потока со взвешенными частицами; τск  – циркуляция слоя в камере. Rн  – среднее значение начального радиуса камеры.

10. Дляdmin  и  dmax  по (13) и (14) определяются радиусы орбит этих частиц rmin  и  rmax

11.Задаются новые границы слоя Rн=rmax  и  Rк=rmin , и расчёт повторяется с п.4. Расчёты проводятся до совпадения результатов с необходимой точностью. В итоге данная методика расчёта взвешенного слоя в камере пневмосмесителей позволяет определить границы слоя, его порозность и скоростные параметры его вращения. По данной методике получены расчётные данные о механике взвешенного слоя, которые представлены в таблице 2. Исходные данные для расчёта показаны в таблице 1.

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

References

1. Uvarov V.A, Orekhova T.N., Gordienko S.I., Kachaev A.E. Continuous pneumatic mixer for the production of dry building mixes. Patent RF, no. 20101140830, 2010.

2. Uvarov V.A., Orekhova T.N., Klyuyev S.V., Kachaev A.E. Countercurrent pneumatic mixer for the production of dispersed-reinforced mixtures. Patent RF, no. 2013159013, 2013.

3. Uvarov V.A., Orekhova T.N., Gordeev S.I., Kachaev A.E. Pneumatic mixer of multi-component dry building mixes. Patent RF, no. 2011151913, 2012.

4. Korneev V.I. Dry building mixes (composition, properties). [Sukhie stroitelnye smesi: sostav, svoystva]. Moscow: RIF "Building Materials", 2010. 320 p. (rus)

5. Banit F.G. Operation, repair and installation of equipment for the construction materials industry. [Ekspluatatsiya, remont i montazh oborudovaniya promyshlennosti stroitelnykh materialov]. Moscow: Stroizdat, 1991. 368 p. (rus)

6. Arratia P.E., Duong Nhat-hang, Muzzio F.J., Godbole P., Reynolds S. A study of the mixing and segregation mechanisms in the Bohle Tote blender via DEM simulations. Powder Technology. 2006. Vol. 164. Pp. 50–57.

7. Berthiaux H., Mizonov V. Applications of Markov Chains in Particulate Process Engineering: A Review. The Canadian Journal of Chemical Engineering. 2004. Vol. 85. No. 6. Pp. 1143–1168.

8. Polyakov E.P., Kardash A.P. Hydrodynamics of a vortex hyperbolic chamber in the presence of a solid phase [Gidrodinamika vihrevoj giperbolicheskoj kamery pri nalichii tvyordoj fazy]. Izv.SO RAN SSSR. Ser.tekhn. nauk. 1984. Vol. 2. No. 10. Pp. 90–98. (rus)

9. Anderson L.A., Hasinger S.H., Turman B.N. Two-component vortex flow studiess of the colloid core nuclear rocket. J. Spacecrafit and rock. 1972. Vol. 9. No. 5. Pp. 311–317.

10. Smulsky J.J. The Aerodynamics and processes in vortex chambers. [Aerodinamika i protsessy v vikhrevyh kamerakh]. Novosibirsk: VO "Nauka". Siberian publishing company, 1992. 301 p. (rus)

11. Uvarov V.A. The counterflow mixer for receiving the disperse reinforced composines. Research Jornal of Applied Sciences. 2014. Vol. 9. No. 12. Pp. 1211–1215.

12. Romanovich A.A., Orekhova T.N., Meshcheryakov S.A., Prokopenko V.S. Technology of producing mineral admixtures. [Tehnologia polychenia mineralnih dobavok]. Bulletin of BSTU named after V.G. Shukhov. 2015. No. 5. Pp. 188–192. (rus)

13. Orekhova T.N., Uvarov V.A. Determination of the particle velocity of the material of the pneumatic mixer of dry building mixes [Opredelenie skorosti chastits materiala pnevmosmesitelya sykhikh stroitelnykh smesey]. Fundamental research. 2013. No.4. Pp. 592–596. (rus)

14. Uvarov V.A., Orekhova T.N. The structural analysis of pnevmonita for the production of dry construction mixtures. [Analiz konstryktsii pnevmosmesiteley dlya proizvodstva sykhikh stroitelnykh smesey]. Interstroymeh. 2010. Pp. 91–96. (rus)

15. Orekhova T.N., Agarkov A.M., Golubyatnikov A.A. Directions of constructive and technological improvement of pneumatic mixers for the production of building materials. [Napravleniya konstryktivno-tekhnologicheskogo sovershenstvovaniya pnevmosmesiteley dlya proizvodstva stroitelnykh materialov]. Scientific almanah. 2015. No. 3. Pp. 124–127. (rus)

16. Orekhova T.N., Krasnov V.V., Demushkin N. P. Pneumatic Mixer for the production of dispersed-reinforced mixtures. News of Science and Education. 2018. Vol. 1. No. 4. Pp. 24–27.

17. Orekhova T.N., Krasnov V.V., Demushkin N.P. Pneumatic vertical mixer. News of Science and Education. 2018. Vol. 1. No. 4. Pp. 21–23.

18. Orekhova T.N., Krasnov V.V., Demushkin N.P. Trends in the development of modern technology and technologies for mixing bulk materials. News of Science and Education. 2018. Vol. 1. No. 4. Pp. 17–20.

19. Orekhova T.N., Krasnov V.V., Demushkin N.P. Classification of mixing equipment for the construction industry. News of Science and Education. 2018. Vol. 1. No. 4. Pp. 13–16.

20. Orekhova T.N., Krasnov V.V., Demushkin N.P. Improvement of devices for mixing dry building mixes. News of Science and Education. 2018. Vol. 1. No. 4. Pp. 9–12.


Login or Create
* Forgot password?