CC BY
16
УДК 621.926
ОПИСАНИЕ ПРОЦЕССОВ ИЗМЕЛЬЧЕНИЯ И КЛАССИФИКАЦИИ СЫПУЧИХ МАТЕРИАЛОВ НА ОСНОВЕ УРАВНЕНИЯ БОЛЬЦМАНА
ЖУКОВ В.П., OTWINOWSKI H., доктора техн. наук, БЕЛЯКОВ А.Н., URBANIAK D. кандидаты техн. наук
Предложено обобщенное описание совмещенных процессов движения, измельчения и классификации сыпучих материалов на основе уравнения Больцмана. Приведено сравнение численных и аналитических результатов расчета процесса гравитационной классификации порошков.
Ключевые слова: уравнение Больцмана, измельчение, классификация, транспорт, сыпучие материалы.
DESCRIPTION OF GRINDING OPERATION AND GRANULAR MATERIALS CLASSIFICATION BASED ON BOLTZMANN EQUATION
V.P. ZHUKOV, H. OTWINOWSKI, Doctors of Engineering, AN. BELYAKOV, D. URBANIAK, Candidates of Engineering
The authors consider the general description of combined movements, grinding and classification processes of granular materials based on the Boltzmann equation. The article shows the comparison of numerical and analytical calculation results of the gravitational classification process of powders.
Keywords: Boltzmann equation, grinding, classification, transportation, granular materials.
Для описания на единой методической основе совмещенных процессов измельчения, классификации и движения сыпучих материалов предлагается использовать уравнение Больцмана [1].
Кинетическое уравнение Больцмана определяет изменение во времени распределения плотности вероятности в выбранном фазовом пространстве. Для порошкообразных материалов важной характеристикой является размер частиц, который предлагается использовать в качестве дополнительной координаты фазового пространства. В качестве искомой функции рассматривается плотность распределения частиц по координатам, скоростям и размерам f (г,у, 5, t). Произведение функции
Цг,\/, 5,0 и фазового объема
= dx1dx2dx3dv1dv2dv3d5 показывает вероятность в момент времени + dt) частицы размера (5,5 + d 5), находящейся в точке с координатами (х1,х1 + dx1), (х2,х2 + dx2), (х3,х3 + dx3), двигаться со скоростью (v1,v1 + dv1), (V2^2 + dv2), ^3^3 + dv3). В мельнице функция
распределения изменяется за счет разрушения и движения частиц, в классификаторе - только за счет движения.
В общем случае изменение функции распределения обусловлено, во-первых, физическим перемещением частиц divr V), во-вторых, изменением скорости частиц divV(af) и, в-третьих, переходом частиц в другой класс крупности за счет разрушения 4. В дифференциальной форме уравнение Больцмана для совмещенных процессов движения, классификации и измельчения принимает вид
df
dt '' или
df dt +
dv1f dv2f dx1 dx2
dv3f da1f
___3___i____!__
dx3
da2f dv1 dv2
da3f dv 3
= f.
dvkf
дХь
dai = fc dvu c
(k = 1, 2, 3),
(1)
где f(r,V, 5,0 - искомая плотность распределения частиц по координатам г(х1,х2,х3), скоростям у^^2^3) и размерам; а(а1,а2, а3) - ускорение частицы; повторение индекса к вверху и внизу дроби показывает суммирование по этому индексу.
Величина 4 в правой части уравнения описывает изменение функции распределения за счет ухода разрушенных частиц наблюдаемого размера и прихода осколков разрушения более крупных частиц. Используя для описания измельчения селективную модель [2], уравнение (1) представляется в виде
дf дvkf да^ 5 Т*, ,
— + —^ + Г fSbdе,
дt я\/ "
Лк к 5
где Б, Ь - селективная и распределительная функции разрушения; 5, е - наблюдаемый и текущий размеры частиц.
Для уравнения (1) предлагается численный метод решения [3], построенный с использованием математического аппарата теории цепей Маркова. Обозначим через N число координат или размерность задачи. Вдоль каждой координаты процесса число выбранных ячеек (или шагов) обозначим п,, где индекс i = 1, 2,..., N показывает порядковый номер координаты. Общее число анализируемых ячеек п будет определяться произведением N чисел п = п1*п2* ...*%. Вектор состояния системы Б = {Б,} размера 1*п составляется из вероятностей попадания в ,-ю ячейку: Б, = f^. Сумма этих вероятностей из условия
нормировки, которое имеет смысл закона сохранения массы, равна единице: £ Б, = 1.
Для описания эволюции системы во времени согласно (1) для каждой ячейки введем характеристику возможных переходов - случайную величину У,, которая определяется множеством значений, которые она может принимать, и множеством соответствующих им вероятностей:
У = {Уь У2,...Ут,},
О, = {Я1,Я2,-Ят/},
где у1,у2,...ут, - номера ячеек, куда возможен переход из ячейки /' с вероятностями д1,д2,...дт, ; т1 - число ячеек, в которые возможен переход из ячейки ¡. В соответствии с условием нормировки сумма вероятностей возможных переходов равна единице: £ д = 1.
В общем случае для описания перехода может быть использована матрица размера п*п, которую, однако, в силу большого числа нулевых элементов целесообразно использовать для одномерных или двумерных задач. Предлагаемый метод расчета содержит адреса ненулевых переходов и соответствующие им вероятности и требует минимального объема машинной памяти, что особенно актуально при решении многомерных задач.
Если переход в многомерном фазовом пространстве разрешен только в соседние ячейки, адреса возможных переходов записываются следующим образом:
У, = {,,, ± 1,, ± п1,, ± п1п2,...,, ± п1п2...п т.-1}. Адрес
означает переход в свою ячейку, адрес У ± 1 - переход в правую или левую ячейки, адрес 1 ± п1 -переход в ячейки, расположенные на строку выше или ниже ячейки с номером ^ и т. д. Вероятности переходов д определяются из физических условий протекания процессов.
Для описания перехода системы из одного состояния в другое используется специальный матричный оператор
Бк+1 = р(У, Бк), (2)
где верхний индекс к показывает номер шага по времени. Матричный оператор р в выражении (2) может быть определен (заменен) следующим тождественным суммированием возможных переходов в рассматриваемую ,-ю ячейку потоков вероятности из всех (п) ячеек системы:
б/+1 = £ ,
]=1, п
где р] - вероятность перехода из ]-й ячейки в ,-ю, которая определяется согласно
Гду 1 при У] 1 =,,
Pii 1 0 при У]1 ф ,, где ]1 = 1,т;.
В качестве тестового примера рассматривается одномерное равноускоренное движение материала под действием силы тяжести. Предлагаются два описания этого процесса: стохастический с помощью уравнения Больцмана и детерминированный в рамках механики Ньютона [4].
Для движения вдоль координаты х уравнение Больцмана (1) переписывается в виде дf д^ дaf
— +----+-----= 0. (3)
дt дх дv
Считается, что материал может переходить только в соседние ячейки. Вероятности переходов рассматриваются в двухмерном фазовом пространстве, которое допускает использование более наглядной двойной индексации ячеек. Для равноускоренного одномерного движения вероятность перехода из ,-й в соседнюю по координате (М^ю ячейку за время At находится из выражения рх = VА^Ах , где V - текущая скорость; Ах - размер ячейки вдоль х. Вероятность перехода в соседнюю по скорости (/+1)-ю ячейку за время At определяется выражением pv = аАЦAv , где Av - размер ячейки вдоль V; а = д - ускорение свободного падения. Уравнение баланса массы для текущей ячейки с номером ] позволяет получить выражение для искомой функции в следующий (к+1) момент времени:
]=Г, - ^ - а*. 1++] а*. .
,, ] ,, ] I Ах Av ) , 1, ] Ах ,, ] 1 Av
(4)
Эволюция плотности вероятности во времени определяется согласно (2). На рис. 1 представлены результаты численных экспериментов по определению закона движения и закона изменения скорости материала при детерминированном (х = х0 + v0t + gt2/2, V = уэ + gt, рис. 1, линии) и стохастическом (согласно (2), рис. 1, точки) описании движения.
В качестве второго примера рассматривается движение частиц материала в среде с линейным сопротивлением (коэффициент сопротивления к). Ускорение материала в уравнении (1) определяется как разность ускорения свободного падания и силы сопротивления, отнесенной к массе: а = д - к^ /т. На рис. 2 представлены результаты расчетного анализа, выполненного для детерминированного и стохастического подходов. На рис. 2 решение дифференциального уравнения движения (х = х0 + (т2/к2)(д - ^0/т)(в~м/т - 1) + тд^к, V = (т/к)(д - ку0т)е_к‘/т) представлено линиями, а результаты численного решения уравнения Больцмана - точками. Анализ зависимостей, приведенных на рис. 1, 2, показывает, что уравнение Больцмана и метод его решения дают результаты, совпадающие с детерминированным описанием движения.
V, м/с :
t. с
°0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 ¡.¡'4
Рис. 1. Сопоставление зависимостей координаты х и скорости V от времени при равноускоренном движении для детерминированного (линии) и стохастического (точки) описания
Рис. 2. Сопоставление зависимостей координаты х и скорости V от времени при движении в среде с линейным сопротивлением для детерминированного (линии) и стохастического (точки) описания
Следует отметить, что предложенный подход на единой методологической основе позволяет описывать и более сложные совмещенные процессы, анализ которых более подробно дадим на примере процессов движения и классификации частиц (по плотности и размеру) в гравитационном аппарате.
Классификация материала по плотности в неподвижной жидкости осуществляется следующим образом. Более плотные частицы двигаются вниз быстрее и, соответственно, быстрее оседают. Менее плотные частицы могут вообще не оседать (например, дерево в воде). На частицу при движении в вертикальном канале действуют сила тяжести и сила Архимеда. Ускорение частицы определяется выражением а = д(Рм - Рж) / Рм, где рм, рж - плотности материала и жидкости соответственно. На контрольной глубине устанавливается поглощающий экран (адсорбер) для улавливания частиц. Для выбранного времени процесса согласно (2) определяется распределение плотности вероятности по фазовому пространству. Вероятность поглощения экраном (адсорбером) определяется через отношение поглощенной экраном массы к массе подаваемой порции. Зави-
симости вероятности поглощения от плотности материала для различных характерных времен процесса представлены на рис. 3.
Рис.3. Зависимости вероятности осаждения от плотности частицы при различном времени процесса: 1 - t = 0,2; 2 -0,25; 3 - 0,3; 4 - 0,35; 5 - 0,45 с
0,005 0,015 0,025 0,035
Рис. 4. Зависимости вероятности осаждения от крупности частицы при различном времени процесса: 1 - t = 0,2; 2 -0,25; 3 - 0,3; 4 - 0,35; 5 - 0,45 с
Процесс классификации материала по крупности частиц описывается согласно уравнению (2) с учетом известной зависимости коэффициента сопротивления от поперечного сечения сферической частицы [4]. На рис. 4 представлены результаты расчета разделения частиц по крупности в соответствии с законом сопротивления Ньютона [4].
Найденные из решения уравнения Больцмана зависимости (рис. 3, 4) являются, по-существу, кривыми Тромпа или кривыми разделения [4], которые позволяют полностью рассчитывать процесс классификации по крупности или плотности частиц в анализируемом аппарате.
Список литературы
1. Вулис Л.А. и др. Теория и расчет магнитогазодинамических течений. - М.: Атомиздат, 1971.
2. Mizonov V., Zhukov V., Bernotat. Simulation Of Grinding: New Approaches. - Ivanovo, 1997.
3.Свидетельство о государственной регистрации программы для ЭВМ «Расчет многомерных совмещенных процессов измельчения, классификации в сыпучих средах
№2010612671 от 19 апреля 2010 года» / А.Н. Беляков, 4. Мизонов В.Е., Ушаков С.Г. Аэродинамическая
В.П. Жуков, А.А. Власюк, А.Е. Барочкин. классификация порошков. - М.: Химия, 1989.
Жуков Владимир Павлович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор кафедры прикладной математики, телефон (4932) 26-97-45, е-mail: zhukov@ispu.ru
Беляков Антон Николаевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, докторант кафедры прикладной математики телефон (4932) 26-97-45.
Otwinowski Henryk,
Ченстоховский политехнический институт (Польша), доктор технических наук, профессор, телефон +48-34-250-579.
Urbaniak Darek
Ченстоховский политехнический институт (Польша), кандидат технических наук телефон +48-34-250-579.
