CC BY
13
УДК 621.928
Моделирование механических процессов в струйной мельнице кипящего слоя на основе уравнения Больцмана
Жуков В.П., Otwinowski H., доктора техн. наук, Беляков А.Н., Urbaniak D., кандидаты техн. наук
На основе уравнения Больцмана предложено обобщенное описание совмещенных процессов движения, измельчения и классификации сыпучих материалов. Приведено сравнение результатов расчетных и экспериментальных исследований совмещенных процессов измельчения и классификации материала в струйной мельнице кипящего слоя.
Ключевые слова: уравнение Больцмана, механические процессы, измельчение, классификация, движение, сыпучие материалы.
Modelling of Mechanical Processes in Jet Mill of Boiling Bed based on Boltzmann equation
V.P. Zhukov, H. Otwinowski, Doctors of Engineering, A.N. Belyakov, D. Urbaniak, Candidates of Engineering
The authors suggest a generalized description of the combined processes of motion, grinding and classification of granular materials based on the Boltzmann equation. The authors compare the calculated and experimental results of the combined processes of grinding and material classification in the jet mill of boiling bed.
Key words: the Boltzmann equation, mechanical processes, grinding, classification, transportation, granular materials.
Спецификой работы струйной мельницы кипящего слоя [1, 2] является совместное протекание в рабочем объеме процессов движения, измельчения и классификации сыпучих материалов. Традиционно [3] описание измельчения и классификации выполняется раздельно, а описание всей установки синтезируется из моделей отдельных процессов. Такой подход оправдан для линейных независимых процессов, для которых возможна такая суперпозиция моделей. Для более общего случая нелинейных и взаимосвязанных процессов предлагается использовать уравнение Больцмана [4], в котором искомая функция распределения массы по выбранным определяющим координатам автоматически учитывает взаимное влияние процессов. Ниже предлагается описание взаимосвязанных процессов в струйной мельнице на основе уравнения Больцмана.
Объектом исследования является струйная мельница кипящего слоя, схема потоков в которой показана на рис.1,а. Исходный материал загружается в размольную камеру, туда же компрессором через сопла (рис. 1,а, стрелки) нагнетается воздух. Потоком воздуха частицы разгоняются и приобретают необходимую для измельчения энергию. При столкновении частиц с другими частицами и со стенками аппарата происходит их измельчение. Мелкие осколки выносятся из слоя и направляются в бункер готового продукта. Крупные частицы остаются в мельнице и продолжают измельчаться заданное время или до достижения заданных технологических параметров.
Рис. 1. Схема потоков (а) и структура фазового пространства (б) в струйной мельнице кипящего слоя
Для построения модели мельницы выбирается набор определяющих координат для искомой функции распределения. Считается, что частицы материала в поперечном сечении мельницы интенсивно перемешиваются, поэтому распределение материала в поперечном сечении не рассматривается. В качестве определяющей геометрической координаты выбирается вертикальная ось х. В качестве второй координаты фазового пространства выбирается скорость движения частицы в вертикальном направлении V. Преобразование размера частиц при измельчении описывается изменением функции распределения вдоль третьей координаты размера зерен 5. Структура трехмерного фазового пространства для искомой плотности распределения в струйной мельнице кипящего слоя приведена на рис.1,б.
Искомая плотность распределения зависит от трех выбранных фазовых координат и времени - і(ху,8,Т) , а уравнение Больцмана [5] для рассматриваемого случая принимает вид
ді
дt'
дукІ
дхи
дв^І = .
дУи
8 тах
ІБ + | ІБЬ6є,
(1)
где Б, Ь - селективная и распределительная функции разрушения [3]; а - ускорение; 8, є -наблюдаемый и текущий размеры частиц.
Решением кинетического уравнения Больцмана является функция распределения массы материала в выбранном фазовом пространстве. Произведение функции І(х,у,8,Т) на элементарный объем = дх ду d8 показывает вероятность в момент времени (Т, t + д) частицы размера (8,8 + д 8), находящейся в точке с координатами (х,х + дх), двигаться со скоростью (у,у + ду). Известный вид функции распределения позволяет определить гранулометрические составы продуктов и их массовые потоки в заданном сечении или объеме, то есть найти технологические характеристики процесса.
Уравнение (1) записывается для получения численного решения в конечно-разностной форме
t+1 t ( V, Л aЛt
+1 = 11 --1-------------— - в Лt
Лх Лу
/, 1,к '/,1 ,к
V: Лt
+п
і-1,1, к'
:
Лх
(2)
хТ ВЛТ ^ _т ,
Іі ’ 1 -1к "ЛУ + ^ 1'і ’ 1 к+Лвк+ЛЬк+л, к
ЛТ ,
Л=1, п-к
строкой, значение каждого элемента которой находится произведением искомой плотности распределения на объем ячейки.
Обозначим через N число координат или размерность задачи. Вдоль каждой координаты процесса число выбранных ячеек (или шагов) обозначим п/, где индекс / = 1,2,.., N показывает порядковый номер координаты. Общее число анализируемых ячеек п будет определяться произведением N чисел: п = п1*п2* ...*п№ Вектор состояния системы S = {Б/} размера 1*п составляется из вероятностей попадания материала в /-ю ячейку: Б; = fД V. Сумма этих вероятностей из условия нормировки, которое имеет смысл сохранения массы, равна единице: ^ Б, = 1.
/'
Для описания эволюции системы во времени согласно (4) для каждой ячейки вводится случайная величина Y¡, которая определяется множеством значений, которые она может принимать, и множеством соответствующих им вероятностей:
V; = У 2,...Уті
Q/ = ^ Я2,..-Ят,},
(5)
где Дх - размер ячейки вдоль оси х; Дv - размер ячейки вдоль оси V; Д5 - размер ячейки вдоль оси 5; нижние индексы /, }, к - номер ячейки вдоль фазовых координат х, V, 5 соответственно; верхний индекс t - номер шага по времени.
Вводя новые обозначения
рх = VДt|Дх, рv = aДt|Дv, р5= БДt (3)
и подставляя их в (2), получаем
1 = fi,1 к (1 - РХ,1 к - Р11 к - Р/51 к) +
+f'-1,],кр/-1,у,к + fi, 1 -1.кр/',у-1,к + fi, 1 к+Др/51,к+ДЬк+Д,к ,
Д=1, п-к
где рх, р¥ ,р5 - вероятности перехода материала из ячейки вдоль координат х, V, 5.
Отрицательные слагаемые в правой части (4) определяют уход материала из ячейки, а положительные - приход в нее частиц из других ячеек.
Для решения уравнения Больцмана (1) с учетом разностной схемы (4) разработан универсальный метод [6]. Многомерное фазовое пространство делится на ячейки аналогично рис. 1,б. Полученные ячейки перестраиваются вдоль прямой таким образом, что адрес ячейки определяется одним, а не тремя индексами. Состояние моделируемой системы задается матрицей-
где у1,у2,.. ут. - номера ячеек, куда возможен переход из ячейки / с вероятностями цьц2,...цт.; - число ячеек, в которые возможен переход из
ячейки /. В соответствии с условием нормировки сумма вероятностей возможных переходов равна единице ^ д = 1.
С другой стороны, для описания эволюции системы во времени может быть использована цепь Маркова и переходная матрица размера п*п, которую, однако, в силу большого числа нулевых элементов целесообразно использовать только для одномерных или двумерных задач. Предлагаемый метод расчета оперирует с адресами только ненулевых переходов и соответствующими им вероятностями и требует минимального объема машинной памяти, что особенно актуально при решении многомерных задач.
Если переход в многомерном фазовом пространстве разрешен только в соседние ячейки, адреса возможных переходов записываются следующим образом: У/ = {/,/ ± 1,/ ± п1,/ ±п1п2,...,/ ±п1п2...п т.-1}.
Адрес / означает остаться в наблюдаемой ячейке, адрес / ± 1 - переход в правую или левую ячейки, адрес / ± п1 - переход в ячейки, расположенные в фазовом пространстве на строку выше или ниже ячейки с номером / (рис. 1,б), и т. д. Вероятности переходов д определяются из физических условий протекания процессов.
Для описания перехода системы из одного состояния в другое за время Д, используется специальный матричный оператор
іТ+1
= Р(У, эт ),
(6)
где верхний индекс t показывает номер шага по времени.
Матричный оператор р в выражении (6) может быть определен (заменен) тождественным суммированием возможных переходов в рассматриваемую /-ю ячейку из всех п ячеек системы:
V+1 =1 ,
]=1, п
где рц - вероятность перехода из ]-й ячейки в -ю, которая определяется согласно
(7)
Рц
Я] 1 при уц 1 = /, 0 при у]1 ф /,
(8)
8;
5/.
7+1
5„
§п+1-0
Рис. 2. Схема материальных потоков (вероятностей) при истирании частиц ]-ой фракции (п - номер самой мелкой фракции)
Р/, ] =
1 - ^ 5] 5к]+1
]+1 + 0 5] 8*п
/ = ],
/ = ] +1,
/ = п,
(9)
]+1
0,
где ]1 = 1,ту.
Для проведения расчетов для струйной мельницы по предложенной методике для каждой ячейки определяются вероятности переходов согласно (3). На частицу в вертикальном канале действуют две силы: тяжести и сопротивления. При известном расходе газа и скорости воздуха в каждом сечении аппарата w(x) сила сопротивления находится из выражения [7]
Гв = к{ ^(х) - и )п ,
где коэффициент сопротивления ^ и параметр п определяются законом сопротивления [7]. Выражение для ускорения частицы из второго закона Ньютона записывается в виде
а = д - ^(х) - и)п.
т
При измельчении материала истиранием считается, что от исходной частицы откалываются только мелкие осколки, а масса частиц фракции после разрушения переходит только в два класса крупности: мелкий и соседний (рис. 2). Вероятности перехода частицы при разрушении из фракции ] во фракцию I за время At или распределительная функция разрушения определяются соотношением [8]
Л+1
/ Ф ], ] + 1 п ,
где к - параметр, соответствующий форме разрушаемых частиц: к = 1 для частиц линейной структуры, к = 2 - пластинчатой, к = 3 - кубической; индексы показывают номера фракций.
Селективная функция измельчения для закона Риттенгера, согласно [3], находится из выражения Б = а8..
Экспериментальные исследования совмещенных процессов проводились для периодического измельчения кварцевого песка в струйной мельнице кипящего слоя. В мельницу загружалась порция материала с известным гранулометрическим составом и измельчалась там в течение 10 мин. После этого определяли массу и гранулометрический состав готового продукта.
При проведении расчетов считалось, что материал вдоль осей х и V может переходить только в соседние ячейки. Число шагов (ячеек) вдоль оси х - 10, число классов скоростей - 10, число классов крупности - 33. Общее число анализируемых ячеек составило 10х10х33 = 3300.
Вероятность перехода частиц по оси х определяется согласно (4) с учетом знака скорости. Если скорость отрицательная, то частицы движутся вниз и полученная вероятность перехода определяет долю частиц, перешедших в нижнюю ячейку, а для положительной скорости - в верхнюю ячейку. Вероятность переходов вдоль оси V определяется согласно (3) с учетом (7). Для определения вероятности перехода по оси 5 используется выражение (8).
Исходный материал подается в ячейку с х = 0 и V = 0. Расчетное время составляет 600 с.
Готовый продукт определяется по содержанию верхних ячеек при положительной скорости движения частиц.
В качестве параметра идентификации модели используется коэффициент селективной функции а. Расчетная доля выноса материала в циклон подбирается равной экспериментальному выносу (18,3 %) за счет выбора значения параметра идентификации.
Сопоставление расчетных и экспериментальных результатов приводится на рис. 3. Результаты численных экспериментов согласно (6) показаны линиями, результаты эксперимента -точками. Штриховая линия соответствует гранулометрическому составу исходного порошка. Результаты численных экспериментов удовлетворительно описывают гранулометрический состав готового продукта.
5
Рис. 3. Сопоставление расчетных (сплошная линия) и экспериментальных (точки) гранулометрических составов (Р(5) -проход через контрольное сито) измельченного материала. Штриховой линией показан гранулометрический состав материала до измельчения
Сопоставление объема оперативной компьютерной памяти, необходимой для решения рассмотренной задачи, с помощью предлагаемого метода и цепей Маркова проводилось следующим образом. Матрица переходных вероятностей при использовании цепей Маркова имеет размер 3300х3300 = 10890000 элементов. Если для хранения каждого элемента требуется 6 байт памяти, то объем оперативной памяти для хранения всей матрицы составляет почти 66 Мбайт. При использовании предлагаемого метода решения объем памяти составляет 3300х6х6х2 = 0,2 Мбайт. То есть
машинные ресурсы при использовании предлагаемого метода экономятся более чем в 100 раз.
При увеличении размерности фазового пространства преимущество предлагаемого метода еще более увеличивается, что делает принципиально возможным решение сложных многомерных задач.
Список литературы
1. Palaniandy S., Azizli K., Hussin H., Hashim S. Effect of operational parameters on the breakage mechanism of silica in a jet mill // Minerals Engineering. - 2008. - 21.
2. Fukunak T., Golmanb B., Shinohara K. Batch grinding kinetics of Ethenzamide particles by fluidized-bed jet-milling // International Journal of Pharmaceutics. - 2006. - 311.
3. Mizonov V., Zhukov V., Bernotat. Simulation Of Grinding: New Approaches. - Ivanovo: Ispeu Press, 1997.
4. Вулис Л.А. Теория и расчет магнитогазодинамических течений. - М.: Атомиздат, 1971.
5. Жуков В.П., Беляков А.Н. Моделирование и расчет совмещенных процессов на основе уравнения Больцмана // Изв. вузов. Химия и химическая технология. - 2010. - Т. 53. -Вып. 11. - С. 114-117
6. Свидетельство о государственной регистрации программы для ЭВМ «Расчет многомерных совмещенных процессов измельчения, классификации в сыпучих средах №2010612671 от 19 апреля 2010 года» / А.Н. Беляков, В.П. Жуков, А.А. Власюк, А.Е. Барочкин.
7. Мизонов В.Е., Ушаков С.Г. Аэродинамическая классификация порошков. - М.: Химия, 1989. - 160 с.
8. Огурцов А.В. Метод расчета струйной мельницы кипящего слоя // Изв. вузов. Химия и химическая технология. -2004. - Т. 47. - Вып. 10. - С. 122-124.
Беляков Антон Николаевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, докторант кафедры прикладной математики, телефон (4932) 26-97-45.
Жуков Владимир Павлович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор кафедры прикладной математики, телефон (4932) 26-97-45, e-mail: zhukov@ispu.ru
Otwinowski Henryk.,
Ченстоховский политехнический институт (Польша), доктор технических наук, профессор, e-mail: otwinowski@kkt.pcz.czest.pl
Urbaniak Darek,
Ченстоховский политехнический институт (Польша), кандидат технических наук, e-mail: otwinowski@kkt.pcz.czest.pl
