Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

THE MODELING OF THE TW0-DIMENSIONAL DISTRIBUTION FUNCTION FOR THE PROBABILITY FOR THE DESPERTION OF THE PERCUSSION MIXING OF THE FRAIBLE MEDIUM

Kapranova A.B. 1 Verloka I.I. 1 Lebedev A.E. 1
1 FGBOU VPO «Yaroslavl State Technical University»
Proposed a accidental method for the forming of the two-dimensional distribution function for the probability for the rarefied flows of the components particles according to its diameters and angles of the deflection in the torches. The translatory and rotatory motions of every stream’s particles are considered at that, including due to its secondary collisions. The semi normal law used possessing of the modeling differential probability for the reflected particles by dispersion angles. In this range of constructive-regime parametres of the device obtained limits tilting flat bump. The resulting expression for the heterogeneity can be used for the formation of an engineering method for calculating the percussion device for mixing bulk materials.
process of the mixing of the components
equipment of the impact action
breaking surface
rarefied flow
secondary collisions
design and regime set-up’s parameters of the equipment
distribution function for the probability by components particles diameters and dispersion angles
semi normal law for the probability
1. Bakin M.N., Kapranova A.B., Verloka I.I. Fundamentalnye issledovanija. 2014, no. 11, pp. 928–933, available at: www.rae.ru/fs/?section=content&op=show_article&article_id= 10003288
2. Bakin M.N., Kapranova A.B., Verloka I.I. Fundamentalnye issledovanija, 2014, no. 5, pp. 687–691, available at: www.rae.ru/fs/?section=content&op=show_article&article_id= 10003246
3. Bakin M.N., Kapranova A.B., Verloka I.I. Fundamentalnye issledovanija, 2014, no. 5, pp. 923–927, available at: www.rae.ru/fs/?section=content&op=show_article&article_id= 10003287
4. Zajcev A.I., Bytev D.O. Udarnye processy v dispersno-plenochnyh sistemah [Shock processes of dispersion-film systems]. Moscow, Chemistry, 1994. 176 p.
5. Kapranova A.B., Bakin M.N., Lebedev A.E. Izv. VUZov. Himija i himicheskaja tehnologija. Ivanovo, 2013. Vol. 56, no. 6, pp. 83–86.
6. Kapranova A.B., Vasilev A.M., Petrov A.A. Matematicheskie metody v tehnike i tehnologijah MMTT-25 (Mathematical Methods in Technics and Technologies MMTT-25). Saratov, 2012, pp. 85–86.
7. Kapranova A.B., Lebedev A.E., Bytev D.O. Izv. VUZov. Himija i himicheskaja tehnologija. Ivanovo, 2004. Vol. 47, no. 6, pp. 99–101.
8. Kapranova A.B., Lebedev A.E., Nikitina Ju.V. Teoreticheskie osnovy jenergo-resursosberegajushhih processov, oborudovanija i jekologicheski bezopasnyh proizvodstv (Theoretical Foundations of energy and resource-saving processes, equipment and environmental security). Ivanovo, 2010, pp. 91–93.
9. Protodjakonov N.O., Bogdanov S.R. Statisticheskaja teorija javlenij perenosa v processah himicheskih tehnologii [Statistical theory of transport phenomena in the processes of chemical technology]. Moscow, Chemistry, 1983. 400 p.

Эффективность процесса смешивания сыпучих материалов в аппаратах с образованием разреженных потоков компонентов значительно повышается при использовании дополнительных отбойных элементов. Актуальная проблема оценки коэффициента неоднородности получаемой смеси может быть решена с помощью стохастического подхода [4, 7, 9] с учетом эффекта ударного взаимодействия формируемых разреженных факелов [2] с отбойной поверхностью заданной формы. Пусть установка имеет вращающийся барабан с радиально расположенными гибкими элементами при частоте вращения n, который имеет некоторый зазор с подвижной лентой. Данная лента одновременно служит транспортером – для сыпучих компонентов при их послойной подаче в указанный зазор и для готовой смеси после отражения факелов разносортных материалов от установленного под углом γ плоского отбойника (рис. 1). Условно можно выделить два вида потоков – налетающие на отбойник после «срыва » твердых частиц дисперсных компонентов с гибких лопастей вращающегося барабана и отраженные от него.

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


Описание распределения частиц сыпучих компонентов в зависимости от их диаметров

Базируясь на стохастическом подходе [1, 3], предлагается ввести двумерные функции распределения вероятностей

kapranova01.wmf i = 1, ..., r

для каждого разреженного потока частиц сыпучих компонентов по их диаметрам Di и углам отклонения δ в i факелах. Каждая из указанных функциональных зависимостей Fi(Di) и Fβs(δ) моделируется отдельно. Считается, что согласно теории ударных процессов в дисперсных системах [4] дифференциальная функция распределения для частиц i-го компонента по их диаметрам моделируется в форме

kapranova02.wmf (1)

для фазового пространства с элементом dΓi = dWxidDi, соответствующим совокупности одной компоненты скорости Wxi центра масс i-й сферической частицы и ее диаметра Di.

Отыскание функций вида Fi(Di) соответствует экспоненциальному убыванию числа частиц каждого из набегающих потоков для смешиваемых компонентов сыпучих материалов dNi в элементе фазового объема dΓi в зависимости от стохастической энергии частицы Ei [5, 6]

kapranova03.wmf (2)

Входящие в выражение (1) суммы кинетических энергий при поступательном и вращательном движениях соответствуют энергиям частиц каждого i-го потока. Дополнительно учитываются случайные моменты импульсов описанных частиц в виде kapranova04.wmf и kapranova05.wmf, где обозначены: их масса mi, диаметр Di и скорость Wi. Заметим, что моменты инерции рассчитываются при вращениях частиц – после «срыва» с гибкого элемента барабана, вращающегося со скоростью n1, и за счет вторичных столкновений [8]. Константа Ai задается нормировкой kapranova06.wmf. Следовательно, учитывая ρTi – истинные плотности смешиваемых компонентов и обозначения kapranova07.wmf; kapranova08.wmf с2i = 10ai, имеем

kapranova09.wmf (3)

где kapranova10.wmf

kapranova11.wmf.

Входящие в (3) параметры имеют размерности [E0i] = Дж; [κ1i] = м4с–2, [κ2i] = кг2м2 и вычисляются из системы энергетических уравнений.

В частности, E0i и κ1i – из уравнений в виде равенств – энергий потоков до столкновения со стенкой (в момент удара каждого из i разреженных исходных компонентов об отбойную поверхность), а также – при падении отраженного факела на горизонтальную ленту (рис. 1).

pic_31.tif

Рис. 1. Упрощенная схема устройства ударного действия: 1 и 2 – дозаторы; 3 и 4 – слои сыпучих материалов; 5 – лента; 6 – барабан; 7 – гибкие элементы; 8 – отбойник; 9 – дисперсная смесь

 

Параметр κ2i – из уравнения сохранения потоков при равенстве концентрации столкнувшихся частиц из i компонентов kapranova12.wmf и разности между концентрациями – полного числа частиц для каждого ингредиента kapranova13.wmf и для общего числа частиц i сортов, не участвующих в столкновениях kapranova14.wmf, когда выполняется условие

kapranova15.wmf

Значение kapranova16.wmf рассчитывается для единичного объема области, в которой предполагаются возможными их взаимодействия, т.е. Zβs – зоны наиболее вероятного присутствия частиц каждого из потоков – налетающих и отраженного. Например, в качестве такой может быть выбран сектор окружности с центром в точке предполагаемого падения усредненного потока, радиусом h и углом раствора π.

Описание распределения частиц сыпучих компонентов в зависимости от угла отражения от плоского отбойника

При выборе полунормального закона для Fβs(δ) – распределения отраженных частиц по углам с учетом нормировочного множителя F0s

kapranova17.wmf (4)

наиболее вероятной областью отражения смешиваемых частиц считается часть описанного сектора для βs = β1 + β2 – суммы углов падения β1 и отражения β2 усредненных векторов перемещений потоков. Таким образом, суммы указанных концентраций kapranova18.wmf и kapranova19.wmf рассчитываются с учетом коэффициента восстановления kv; концентраций числа не столкнувшихся частиц i-го налетающего потока с элементами j-го налетающего, а также – отраженного.

Формирование критерия качества смеси

Критерий качества смеси выбран в форме коэффициента неоднородности

kapranova20.wmf (5)

когда i = 1, 2 и объемные концентрации равны

kapranova21.wmf (6)

kapranova22.wmf (7)

где kapranova23.wmf – усредненная скорость налетающего потока; Zβs = hL/2; Zβs = πh2/2; L – расстояние между точкой касания барабана с поверхностью сыпучего компонента на ленте и точкой крепления отбойника на горизонтальной поверхности.

Входящие в выражение (7) интегральные выражения с учетом вида Fβs(δ) – распределения частиц по углам в применении соответственно к налетающему и отраженному потокам имеют вид

kapranova24.wmf (8)

kapranova25.wmf (9)

где значения нормировочных параметров F01 и F02 равны

kapranova26.wmf (10)

kapranova27.wmf (11)

Подстановка в формулу (5) для коэффициента неоднородности выражений (8), (9) для указанных распределений частиц с учетом (10), (11) приводит к следующему представлению:

kapranova28.wmf (12)

где kapranova29.wmf

kapranova30.wmf

когда при ϑ1 = ϑ2 = 30, ϑ3 = 20, ϑ4 = 10, ϑ5 = 4 введены обозначения

kapranova31.wmf

Результаты исследования и их обсуждение

Введение безразмерного параметра ε = D/Dmin2 позволяет согласно (3), (4) получить следующее представление для Qi(ε, δ) – общей функции распределения частиц по диаметрам и углам (рис. 2), применяемой в соответствии с (5)–(7) для расчета Vc (рис. 3) в зависимости от режимных параметров аппарата – частоты вращения барабана n1 и угла наклона отбойной поверхности γ. Для сравнения на рис. 4 приведем семейство кривых для qi(ε) – одномерных дифференциальных функций распределения по параметру ε.

pic_32.tif

Рис. 2. Зависимости Qi(ε, δ) и Qm(ε, δ): 1 – каолин ГОСТ 21235-75 (i = 1); 2 – технический углерод П803 ГОСТ 7885-86 (i = 2); 3 – для их смеси

pic_33.tif

Рис. 3. Зависимость для коэффициента неоднородностей Vc(n, γ) при смешении каолина ГОСТ 21235-75 и технического углерода П803 ГОСТ 7885-86 в устройстве ударного действия

pic_34.tif

Рис. 4. Зависимости qi(ε) и qm(ε): 1 – каолин ГОСТ 21235-75 (i = 1); 2 – технический углерод П803 ГОСТ 7885-86 (i = 2); 3 – для их смеси

Пересечение диапазонов диаметров частиц для рабочих сыпучих материалов (каолина ГОСТ 21235-75 и технического углерода П803 ГОСТ 7885-86) приводит к промежуточному положению соответствующей кривой qm(ε) (кривая 3, рис. 4) или поверхности Qm(ε, δ) (поверхность 3, рис. 2) для смеси порошков. При этом наблюдается некоторое снижение значений двумерных распределений вероятностей Qi(ε, δ) и Qm(ε, δ) по сравнению с одномерными – qi(ε) и qm(ε), что объясняется наличием дополнительного перераспределения частиц описанных потоков по углу γ. Полученные распределения qi(ε) и Qi(ε, δ) позволяют различными способами оценить качество приготавливаемой смеси при выборе коэффициента неоднородностей в виде критерия (5).

Конструктивными и режимными показателями установки являются: D1 = 7,5∙10–2 м; l = 4,00∙10–2 м; L = 0,45 м; n = 300 об/мин; Wi = 10–2 м/с; α = 40°; γ = 75°. С увеличением угла наклона отбойника до значения π/2 происходит значительное улучшение качества продукта, что можно объяснить расширением условной зоны перекрещивания потоков компонентов.

Заключение

Примеры вычисления Vc по формуле (5) удовлетворительно описывают процесс смешения порошков, что подтверждается опытными данными: при n = 300 об/мин и изменении угла наклона отбойной поверхности в пределах γ = (0,54–1,72) рад экспериментальные значения имеют диапазон Vc = (19,4 – 23,2). Выражения (5)–(7) могут быть использованы при формировании инженерной методики расчета устройства ударного действия для смешивания сыпучих материалов.

Рецензенты:

Бачурин В.И., д.ф.-м.н., профессор кафедры «Высшая и прикладная математика» Ярославский филиал, ФГБОУ ВПО «Московский государственный университет путей сообщения (МИИТ)», г. Ярославль;

Мурашов А.А., д.т.н., профессор кафедры математики, Ярославское высшее военное училище противовоздушной обороны, г. Ярославль.