The dynamics of a dense layer of granular medium in the circulating pneumatic apparatus.pdf Пневматические методы переработки порошковых и гранулированных сред нашли широкое применение при измельчении, дозировании, смешении, сушке, хранении и транспортировании зернистых материалов в современной химической технологии, порошковой металлургии, ядерном топливном цикле, фармацевтической и пищевой промышленности. Разработанная авторами математическая модель применяется для расчета гидродинамики и процесса смешения высококонцентрированной гранулированной среды в пневматическом циркуляционном смесителе (рис. 1, а). Работа аппарата происходит следующим образом. В начальный момент времени ключевой и основной компоненты смеси находятся в несмешанном состоянии. Под действием струи газа из сопла 1 частицы из области I поднимаются вверх по вертикальной трубе 2, затем ударяются об отбойник 3 и равновероятно рассеиваются по поверхности плотного слоя (область II). В результате непрерывной циркуляции компонентов гранулированного материала в корпусе бункера 4 происходит их смешивание по ип б О Рис. 1. Пневматический циркуляционный аппарат: принципиальная схема (а); геометрия цилиндрической области аппарата: б, в - вид сбоку и сверху Fig. 1. Circulating pneumatic apparatus: (a) scheme of work; (б) cylindrical area of the apparatus; (в) top and side view всему объему аппарата. Скорость движения материала регулируется расходом газа из сопла 1. Запыленный поток газа фильтруется через фильтр 5 и возвращается в аппарат. Такой циклический метод смешения удобен тем, что процесс смешивания можно проводить до требуемой однородности смеси. Для повышения интенсивности процесса смешения компонент смеси на цилиндрической стенке аппарата могут дополнительно устанавливаться полки в виде сегментов (рис. 1, б), которые, наряду с осевой компонентой скорости, вызывают появление радиальной и окружной составляющих вектора скорости. В работе приводится исследование по влиянию этих сегментов на процесс смешения зернистой среды. При теоретическом описании плотного слоя зернистой среды в зависимости от скорости сдвига наблюдают два различных режима течения. Известно, что при небольших скоростях сдвига инерцией частиц можно пренебречь, и возникает так называемый квазистатический режим течения, который определяется кулонов-ским трением и описывается в рамках теории предельного равновесия. В другом предельном случае при больших скоростях сдвига имеет место так называемый инерционный режим течения. В этом режиме, с одной стороны, гранулы находятся в постоянном контакте, а с другой стороны, между ними имеются зазоры и взаимодействие между ними обусловлено неупругими соударениями. Описание этого инерционного режима основывается на законах сохранения массы, импульса и энергии хаотического движения частиц [1]. Этот инерционный режим движения высококонцентрированной гранулированной среды относят к научному направлению, получившему название «теория быстрых движений гранулированных сред» [2]. Как показывают экспериментальные исследования [3], движение плотного слоя гранулированной среды в пневматическом циркуляционном смесителе осуществляется в инерционном режиме течения. Анализ результатов экспериментальных и теоретических исследований «быстрых движений гранулированных сред» позволяет при постановке задачи сделать ряд упрощающих допущений. Будем считать, что хорошо сыпучая гранулированная среда состоит из твердых одинаковых шарообразных частиц, находящихся во взаимном контакте, причем объем отдельной гранулы мал по сравнению с объемом сыпучего тела и порозность плотного слоя по объему практически постоянная. Поэтому такую среду можно считать однородной и несжимаемой. Известно, также, что плотный слой гранулированной среды имеет кривую течения, отличную от ньютоновской, и относится к реостабильной дилатантной жидкости [2]. В исследовании [4] показано, что при движении плотного слоя несвязанной, хорошо сыпучей зернистой среды тензор напряжений незначительно отличается от тензора напряжений ньютоновской среды, прежде всего вследствие ее скольжения на стенках в вертикальных трубах. Таким образом, вязкость среды представим в виде суперпозиции постоянного ее значения и нелинейного слагаемого В общем случае значение консистенции материала ^ зависит от второго инварианта тензора скоростей деформаций A, который в цилиндрической системе можно представить в виде В прикладных задачах по реологии, например [5], используется величина интенсивности тензора скоростей деформаций равная J=VA . На основании сказанного выше, и используя степенную модель для нелинейного слагаемого, для тензора внутренних напряжений можно записать ту. = 2(д 0 + , (1) где n - индекс течения, e^ - тензор скоростей деформаций. При моделировании гидродинамики и процессов смешения в пневматическом циркуляционном аппарате будем учитывать только кольцевую цилиндрическую область с плотным слоем, так как она значительно больше по объему конической части аппарата. Безразмерную форму уравнений движения и уравнение неразрывности в цилиндрической системе координат с учетом реологического уравнения (1) можно представить в виде dur + 5u- +1 duru(f+duruz дГ1+bdur 1 дГ 1+b dur 1 дГ1+bdur dt dr r дф dz dr ^ Re dr ) дф^r16 Re дф) dz^ Re dz ) = _dp _и2 + «,2 + 1 db dur + 1 db duz + Гdb du^ _db_u±+b dur_ _2b db dUф bur dr r r Re dr dr Re dz dr Re r ^дф dr дф r dr r дф дф r ) + duru, +1du2 + дуг1+ь dГ1+ь^0= (2) dt dr r дф dz dr^ Re dr ) дф^r2Re дф ) dz^ Re dz ) 1 dp uru 1 ifdbdur db 1 dbdu( 2db dbduz 2bdur , d b 1 =----_2-^+--1--r---uф+-----ur +--z+--r-+b-u_, I ; r дф r Rer ^dr дф dr r дф дф r дф dz дф r дф dr r ) duz duru_ 1 dUфUz du? d Г1+b duz 1 дГ 1+b duz 1 dfi+b du, - +--+-----dt dr r дф dz dr ^ Re dr ) дф^r2 Re дф) dz ^ Re dz dp uu 1 Г db dur 1 db duф db duz b du. - r z-+-I---+---- +--- +- z dz r Re^dr dz rдф dz dz dz r dr dUr 1 диф dUz Ur -L+--- +--+-=0 . dr r дф dz r X Г H Г rn_1. rH.D _pHU0 Здесь b - I J1n_1; J = J-; Re = p- K 0 lU 0) U0 В качестве масштаба скорости и длины выбираются соответственно: среднее значение аксиальной скорости на входе в исследуемую область - U0; расстояние между внутренней транспортной трубой и корпусом аппарата H = r2-r0. Существенной особенностью предложенной модели движения гранулированной среды является постановка граничных условий на твердых поверхностях. На стенке обычно используют условие прилипания среды, которое, по-видимому, является причиной неудовлетворительного соответствия теоретических и экспериментальных данных. Опытные исследования движения гранулированной среды в вертикальных каналах показывают, что скорость перемещения дисперсного материала на стенке не равна нулю, т.е. имеет место некоторое скольжение среды по стенке. Это явление, по-видимому, можно объяснить преобладанием внутренних напряжений в среде по сравнению с силами трения частиц о твердую поверхность. В данной модели для учета эффекта скольжения среды на стенке принимается, что касательные напряжения на стенке пропорциональны скорости движения среды на ней с точностью до некоторой постоянной р, которую назовем коэффициентом скольжения. Граничные условия на твердой поверхности будут иметь вид = -Put| dut dn Un = 0 (3) w w Здесь ut, un - касательная и нормальная составляющие вектора скорости на стенке соответственно, n - нормаль к поверхности стенки. Очевидно, что при р = 0 имеем условие полного скольжения, а в случае р^-да - условие прилипания. Следует отметить, что коэффициент скольжения среды р, вероятно, зависит от многих факторов и определяется из сопоставления результатов расчета модели с опытными данными. При решении задачи использовались следующие граничные условия. На входе задавалось постоянное безразмерное значение аксиальной скорости uz = -1, на стенках аппарата использовались условия скольжения (3), причем на верхних поверхностях горизонтальных полок ставились условия прилипания для моделирования эффекта застойных зон. В выходном сечении ставились условия Неймана для всех составляющих вектора скорости. Решение системы уравнений переноса импульса проводилось методом расщепления в переменных «скорость - давление». Для построения пространственного разностного шаблона использовалась разнесенная шахматная сетка с применением контрольного объема. Разностный аналог уравнений переноса импульса и концентрации решался при помощи метода факторизации. Конвективные и диффузионные слагаемые в этих уравнениях записывались с помощью известной экспоненциальной схемы. Работоспособность и адекватность предложенного подхода проверялась путем сравнения результатов расчета с экспериментальными данными [6] для установившегося течения полистироловых частиц в вертикальной цилиндрической трубе. Результат сравнения численных и опытных данных для вертикальной составляющей вектора скорости представлен на рис. 2. U -0.4- -0.6- -0.8- -1 -0.5 0.5 0 Рис. 2. Распределение вертикальной скорости в круглой трубе при установившемся режиме течения в сравнении с опытными данными [6]. Расчет проведен при параметрах Re = 10; в = 10; n = 1.2 Fig. 2. Comparison of the axial velocity distribution of the steady flow in a cylindrical pipe and experimental data [6]. Flow parameters are Re = 10; в = 10; and n = 1.2 Из графика сравнения видно хорошее совпадение расчетного профиля скорости и обобщенных экспериментальных значений. На рис. 3 показаны схематично сечения контроля составляющих скорости при симметричной установке трех препятствий (а) и одного препятствия (б) на правой стенке в аппарате. Рис. 3. Схема сечений по угловой координате ф при z = z2: сечение 1 - r = ri/2; сечение 2 - r = (ri+r2)/2: а - три симметричных тарелки; б - одна тарелка Fig. 3. Sketch of the sections through the angular coordinate ф at z = z2: section 1, r = r1/2; section 2, r = (r1+r2)/2: (a) three symmetrical plates and (б) one plate На рис. 4 - 6 показано распределение проекций вектора скорости на уровне передней по потоку стенки в зависимости от угловой координаты для аппарата с тремя полками (а) и одной полкой (б). Профили скорости в сечениях 1, 2, приведенных на рис. 3, показывают, что происходит существенное увеличение переноса импульса в радиальном направлении к внутренней стенке аппарата. Полученные симметричные распределения профилей скорости для геометрии с тремя полками позволяют также сделать вывод о достоверности численного расчета. На рис. 7 показаны изолинии проекций скорости в плоскости r - ф при z = z1+z2 (средняя поперечная движению плоскость тарелок), полученные при расчете гидродинамики сыпучего материала с аналогичными параметрами потока. Рис. 4. Профиль радиальной скорости при Re = 20; в = 10; n = 1.2 в сечениях 1 - 2, показанных на рис. 3: а - три симметричных тарелки; б - одна тарелка Fig. 4. Radial velocity profile at Re = 20, в = 10, and n = 1.2 in sections 1-2 shown in Fig. 3: (a) three symmetrical plates and (б) one plate и и ф 0.4 ■ 0- -0.4 ■ -0.8- ~I 2 п -1 2п 0 ф 0.4 0 -0.4 - -0.8 0 Рис. 5. Профиль окружной скорости при Re = 20; в = 10; n = 1.2 в сечениях 1 - 2, показанных на рис. 3: а - три симметричных тарелки; б - одна тарелка Fig. 5. Peripheral velocity profile at Re = 20, в = 10, and n = 1.2 in sections 1 - 2 shown in Fig. 3: (a) three symmetrical plates and (б) one plate Uz u z а б 2 2 -0.4-0.8 -1.2-1.6-2.0- -0.4 -0.8 -1.2 -1.6- -1 2п ~I 2п 0 0 Рис. 6. Профиль аксиальной скорости при Re = 20; в = 10; n = 1.2 в сечениях 1 - 2, показанных на рис. 3: а - три симметричных тарелки; б - одна тарелка Fig. 6. Axial velocity profile at Re = 20, в = 10, and n = 1.2 in sections 1 - 2 shown in Fig. 3: (a) three symmetrical plates and (б) one plate r 0.8 0.4 0.0 -0.4 -0.8 -1.2 -1 -0.6 -0.2 0.2 0.6 1 -1 -0.6 -0.2 0.2 0.6 1 r -1 -0.6 -0.2 0.2 0.6 1 Г Г Рис. 7. Изолинии проекций вектора скорости в аппарате с тремя симметричными полками при Re = 20; в = 10; n = 1.2: а, б, в - проекции скорости ur, иф, uz соответственно Fig. 7. Contours of the velocity vector projections in the apparatus with three symmetrical plates at Re = 20, в = 10, and n = 1.2: (a-в) velocity projections ur, иф, uz, respectively Из представленных графиков и изолиний проекций вектора скорости можно видеть, что при выбранной геометрии аппарата наибольшие градиенты скорости возникают вблизи боковых стенок препятствий. При установке одной тарелки возмущения скорости близки к возмущениям в потоке при установке трех препятствий. Представленные результаты показывают адекватность предложенного метода расчета гидродинамики зернистых материалов при инерционном режиме течения. Разработанная трехмерная модель движения может использоваться для получения полей скорости в пневматических аппаратах с плотным слоем при моделировании процесса смешения и усреднения сыпучих материалов.

Скачать электронную версию публикации