D01:10.6060/snt.20216501.0009 УДК 542.223
КИНЕТИКА ИСТИРАНИЯ ЗЕРНИСТЫХ МАТЕРИАЛОВ В ПСЕВДООЖИЖЕННОМ СЛОЕ
С.В. Воробьев, М.Ю. Колобов, И.В. Постникова, В.Н. Блиничев
Сергей Владимирович Воробьев, Михаил Юрьевич Колобов, Ирина Викторовна Постникова, Валерьян Николаевич Блиничев
Кафедра технологических машин и оборудования, Ивановский государственный химико-технологический университет просп. Шереметевский,7, Иваново, Российская Федерация, 153000
E-mail: ch4rlythec4t@gmail.com, mikhailkolobov@rambler.ru, poirvic@mail.ru
Представлены результаты экспериментального исследования процесса истирания зернистых материалов, на основании которых предложено уравнение, позволяющее рассчитать массу продукта в зависимости от конструктивных параметров газораспределительной решетки, технологических параметров и физико-механических свойств материала с возможностью применения его для расчета производительности аппаратов непрерывного действия одного -, двух- и трехступенчатого типов.
Ключевые слова: псевдоожиженный слой, измельчение, кинетика истирания
THE KINETICS OF ABRASION OF GRANULAR MATERIALS IN A FLUIDIZED
BED LAYER
S.V. Vorobiev, M. Yu. Kolobov, I. V. Postnikova, V.N. Blinichev
Vorobyov Sergey Vladimirovich, Kolobov Michael Yurievich, Postnikova Irina Viktorovna, Blinichev Valerian Nikolaevich
Department of technological machines and equipment, Ivanovo state University of chemical technology, prosp. Sheremetyevo, 7, Ivanovo, 153000
The results of an experimental study of the process of abrasion of granular materials are presented, on the basis of which an equation is proposed that allows calculating the mass of the product depending on the design parameters of the gas distribution device, technological parameters and physical and mechanical properties of the material with the possibility of using it to calculate the performance of continuous-action apparatus of one-, two-and three-stage type.
Keywords: fluidized bed layer, grinding, abrasion kinetics.
Введение
Физическую модель процесса измельчения твердых частиц в псевдоожи-женном слое можно представить следующим образом: частицы твердой фазы разрушаются вследствие двух параллельно протекающих процессов: истирания частиц в псевдоожиженном слое и разрушения их в результате соударений.
Таким образом, в псевдоожиженном слое можно выделить две зоны, в которых измельчение частиц происходит двумя разными способами. Первая зона - надреше-точная (до 0,05 м), в которой совместно с истиранием частицы испытывают значительные ударные нагрузки вследствие «соплового эффекта струй» энергоносителя. В связи с чем, процесс разрушения частиц на крупные части протекает, в основном,
лишь в надрешеточной зоне псевдоожи-женного слоя.
Вторая зона охватывает остальную область псевдоожиженного слоя, лежащую выше первой зоны. Уменьшение линейного размера частиц 5 во второй зоне происходит вследствие истирания материала в результате интенсивной циркуляции частиц твердой фазы. Разрушение частиц на крупные части в этой зоне не происходит вследствие сравнительно небольшой кинетической энергии их движения.
Наглядное отображение этих зон показано на рис. 1, на котором представлен результат численного моделирования распределения поля скоростей газовой фазы в программе Solidworks Flow Simulations. В качестве исходных данных принято: диаметр аппарата dan = 0,13 м; диаметр отверстий газораспределительной решетки dome = 0,003 м; значение живого сечения решетки ^ = 5%; рабочая скорость газовой фазы ьг = 3 м/c; энергоноситель - воздух.
J
55.00 50.00 45.00 40.00 35.00 30.00 25.00 20.00 15.00 10.00 5.00 0
Скорость [mis] Скорость [m/s;
Рис. 1. Результаты численного моделирования распределения поля скоростей газовой фазы в объеме
псевдоожиженного слоя: зоны измельчения:1 - зона разрушения и истирания; 2 - зона истирания
Значимость той или иной зоны для процесса измельчения может изменяться в зависимости от параметров псевдоожижаю-щего агента, физико-механических свойств материала и конструкции газораспределительной решетки.
Методика эксперимента
Для исследования влияния конструкции газораспределительной решетки на процесс истирания зернистых материалов в псевдоожиженном слое нами были проведены эксперименты по истиранию материалов на решетках с различными живыми сечениями 1 %, 2 %, 3 %, 5 %, 10 % и 20 %. Для каждого живого сечения варьировались диаметры отверстий йотв: 0,001 м; 0,002 м; 0,003 м; 0,004 м; 0,005 м; 0,006 м и
0,008 м. В качестве модельного материала использовались частицы известняка различных фракций. Замеры скорости воздуха внутри аппарата производились с помощью термоанемометра смарт-зонда Testo 405 i. Виртуальное моделирование гидравлики внутри аппарата проводилось в виртуальной среде Solidworks Flow Simulations.
Результаты и их обсуждение
Из литературных источников известно, что интенсивность и стабильность многих процессов химической технологии, осуществляемых в аппаратах с псевдоожи-женным слоем, определяются гидродинамической обстановкой, которая в значительной степени зависит от работы газорас-
пределительной решетки. Таким образом, характер истирания зернистых материалов и преобладание той или иной зоны измельчения в аппарате комбинированного действия во многом будет зависеть от конструкции газораспределительных решеток.
0,006 П
0,005 -
0,004 0,003 -0,002 -|
0,001 -
На рис. 2 приведены кинетические кривые истирания СаСОз с изменениями живого сечения ф газораспределительных решеток при постоянной скорости псевдо-ожижащего воздуха ьг=3 м/с.
1 2
3
4
~1-1-1-1-1-Г-
600 1200 1800 2400 3000 3600 4200 4800
т, с
Рис. 2. Влияние живого сечения ф газораспределительной решетки на максимальную скорость истирания Утах частиц СаСОз; исходные данные: иг=3м/с; йотв=0,002 м; У=0,002^0,0025 м; #0=0,1 м;
живое сечение ц>\ 1-3%; 2-5%; 3-10%, 4-20%.
Как видно из графика, при уменьшении живого сечения (р решетки от 20 % до 3% максимальная скорость истирания Утах частиц слоя возрастает (под скоростью истирания У будем понимать отношение приращения пылевидного продукта в процентах по весу слоя к единице времени - секунде). Это связано со значительным ростом скорости истечения воздушных струй из отверстий решетки, а значит и с увеличением кинетической энергии Екин, воспринимаемой частицами от среды; возрастанием при этом подвижности частиц, числа и силы соударений между ними в надрешеточной области, в результате чего возникает непрерывное дробление частиц на более мелкие части и интенсивное истирание поверхности как мелких, так и крупных частиц.
Помимо живого сечения ф на интенсивность истирания частиц большое влияние также оказывает диаметр отверстий ^отв решетки. Исследование влияния диаметра отверстий ё.отв решетки на скорость истирания У зернистых материалов проводилось при значении живого сечения ^=3%, рекомендуемом рядом авторов, как оптимальное для цилиндрических аппаратов.
Обработка экспериментальных данных показывает, что с увеличением диаметра отверстий йотв от 0,001 до 0,0045 м максимальная Утах и постоянная УС0П51 скорости истирания частиц слоя непрерывно растут (рис. 3). Это объясняется повышением дальнобойности и мощности струй воздуха, вытекающих из отверстий решетки.
0
0
0,025 -,
0,02 -
0,015 -
0,01 -
0,005
0 0,0005 0,001 0,0015 0,002 0,0025 0,003 0,0035 0,004 0,0045 дотв, м
Рис. 3. Влияние диаметра отверстий йоте газораспределительной решетки на изменение скорости истирания У частиц СаСОз; исходные данные: иг=3 м/с, ^ =3%; йтв=0,002^0,0025 м, #0=0,1 м; скорость истирания У: 1- Утах; 2-УС0ПХ{
При максимальном диаметре отверстий dome =0,007^0,008 м мощность струй воздуха достигает такого значения, что слой из псевдоожиженного состояния переходит в режим фонтанирующего псевдоожижения. При тонких слоях (Н0 =0,1^0,15 м) можно визуально наблюдать фонтанирование слоя, причем, обычно число фонтанов равно числу отверстий решетки.
При уменьшении диаметра dome, увеличиваются скорость истечения и мощность струй энергоносителя, истекающих из от-
верстий решетки, что приводит к более интенсивному энергообмену фаз, увеличению силы и числа соударений между частицами. Соответственно, с увеличением диаметра отверстий dome решетки растет КПД псевдоожиженного слоя - как измельчителя.
Результаты анализа полученных экспериментальных данных можно обобщить и представить в виде функциональной зависимости:
AG
йтв, ртв, (
слоя ®сж> ^отв'
V ,т)
(1)
где: АС - количество пыли, кг; Сслоя- вес слоя, кг; ртв и рс - плотность твердого материала и плотность среды, кг/м3; д - динамическая вязкость псевдоожижающей среды, м2/с; g - ускорение свободного падения,
м/с2
предел прочности материала на
сжатие, МПа.
Метод анализа размерностей и обработка экспериментальных данных позволяет привести функциональную зависимость к критериальному виду:
AG
1,75 • 103 Fr1,3 Ho0,6n71,3 )
1 \Рт„/
1,3 0,21Йотв
10 ^тв
18ф
(2)
В критериях Fr =
V2 Ьг X
-и Но = — за
9 ^тв ^тв
определяющий линейный размер принят средний начальный диаметр частицы dТв. п1 = ^ - комплекс, отражающий соотно-
шение между пределом прочности материалов на сжатие и кинетической энергией потока; 10 "тв - комплекс, отражающий особенности конструкции газораспределительной решетки.
1
2
0
Формула (2) проверена в пределах: Fr = 40 + 3000; Но = 0 + 6000000; ^ =
^тв
0,2 + 5,0; ф = 0,006 + 0,05; п1 = 1,4 • 105 - 1,5 • 105; — = 0,07 + 63%. Раскры-
вая значение критериев и комплексов, уравнение (2) можно привести к более простому виду:
AG =
90ФвслояУг
5,8
Рс2,610~
0,21da
-18ф
(3)
1,9
-0,6
а 1,3 о 1,3
"сж Итв
Для того чтобы уравнение (3) можно было использовать для расчета процесса истирания частиц в аппарате непрерывного действия, необходимо ввести в него поправочный коэффициент К, учитывающий отклонения величины действительного времени пребывания, а также вместо т подставить среднее время пребывания частицы в аппарате.
Поправочный коэффициент К для аппарата непрерывного действия равен: одноступенчатый аппарат - К = 0,894; двухступенчатый аппарат - К = 0,947; трехступенчатый аппарат - К = 0,967.
Выводы
1. Надрешеточная зона является поставщиком не только тонкодисперсных частиц, но также и большого количества крупных осколков неправильной формы разрушившихся частиц, энергия которых расходуется на интенсивный процесс истирания острых углов осколков в объеме вихревого движения частиц псевдоожиженного слоя.
2. Обобщая исследования газораспределительных решеток, необходимо отметить, что если нужно иметь сравнительно однородное псевдоожижение и минимальное истирание частиц, то в цилиндрическом аппарате должна быть решетка с живым сечением 3-5 % и с возможно меньшим диаметром отверстий.
3. Проведена обработка экспериментальных данных по влиянию конструктивных параметров газораспределительной решетки, технологических параметров и физико-механических свойств материала. Получена эмпирическая зависимость, позволяющая рассчитать массу продукта
ЛИТЕРАТУРА
1. Петров В.Н., Шабалин А.С., Сопин В.Ф., Петров С.В., Малышев С.Л. Исследование взаимодействия турбулентной осесимметричной струи с преградой // Вестник Казанского технологического университета, т. 20. №. 2. 2017. С. 93-96.
2. Giersemehl, M., G. Plihal. Fine Grinding System with Impact Classifier Mill and Cyclone Classifier. Company Affiliation: Neuman & Esser GmbH Mahl-und Sichtsysteme. Journa l: powder handling & processing. 1999 V. 11. N: 3 P. 269.
3. Berthiaux H., Dodds J. Modelling fine grinding in a fluidized bed opposedjet mill Part 1: Batch grinding kinetcs. Powder Techn. 1999. N 106, P. 78-87.
4. Heng P.W.S., Chan L.W., Lee C.C. Ultrafine grinding using a fluidized bed opposedjet mill: effect of process parameters on the size distribution of milled particles // S.T. P. PharmaSciences. 2000. N 10, P 445451.
5. Vogel L., Peukert W. Modelling of grinding in an Air Classifier Mill Based on a Fundamental Material Function // KONA. 2003. N 21. P 109-120
6. Агапов Ю.Н., Лукьяненко В.И., Бырдин А.П., Стогней В.Г. Влияние высоты псевдоожижен-ного слоя на параметры газораспределительной решетки // Вестник ВГТУ. Серия: Энергетика. Т. 2. №6. 2006. С. 139-143.
7. Агапов Д.Ю., Надеев А. А., Родионов К.Н., Стогней В.Г., Клейников И.Ю. Экспериментальное сравнение газораспределительных устройств для формирования центробежного слоя. Вестник Воронежского государственного технического университета. V. 8. N. 7-1. 2012. P. 118-121.
8. Калашникова Н.В., Волженцев А.В. Совершенствование технологического процесса сушки зерна пшеницы с обоснованием параметров сушилки с псевдоожиженным слоем // Вестник аграрной науки. V. 16. N. 1. 2009. P. 44-45.
9. Волженцев А.В. Определение качества ожижения зерна в сушилках // Образование, наука и производство. N. 4 (9). 2014. P. 4-7.
10. Гвоздев Д.В., Сальников А.А., Фомичев А.Г. Расчет и конструирование аппаратов с псевдоожи-женным слоем зернистого материала. Газораспределительные решетки // Изв. МВО СССР, Серия Химия и химич. тех. № 4. 1963. с.38.
11. Грек Ф. З. О геометрических характеристи-
ках газораспределительных решеток. // Хим. пром. 1966. с.67-68
12. Афанасьева О. В., Мингалеева Г. Р. Аппарат для получения активированного угля на мини-тэс. Известия Томского политехнического университета // Инжиниринг георесурсов. Т. 330. №. 9, 2019. P. 4250.
13. Тодес О. М. Аппараты с кипящим зернистым слоем: (Гидравл. и тепловые основы работы) Л.: Химия: Ленингр. отд-ние. 1981. 296 с.
14. Яковлев В. А., Сахибгареев А. Ф. Аэродинамические особенности проектирования смесепод-готовительных систем щелевых ГГУ с принудительной подачей воздуха. Academy. N. 7 (34). 2018. P. 2325.
15. Белоглазов И.И., Шариков Ю.В., Белоглазов И.Н. Изучение распределения времени пребывания продуктов обжига в печи кипящего слоя Записки Горного института, V. 169. 2006. P. 61-62.
Работа выполнена при финансовой поддержке гранта РФФИ № 19-03-00787
REFERENCES
1. Petrov V.N., Shabalin A.S., Sopin V.F., Petrov S.V., Malyshev S.L. Investigation of the interaction of a turbulent axisymmetric jet with an obstruction Vestnik of Kazan Technological University, vol. 20. No. 2. 2017. S. 93-96.
2. Giersemehl, M., G. Plihal. Fine Grinding System with Impact Classifier Mill and Cyclone Classifier. Company Affiliation: Neuman & Esser GmbH Mahl-und Sichtsysteme. Journa l: powder handling & processing. 1999 V. 11. N: 3 P. 269.
3. Berthiaux H., Dodds J. Modelling fine grinding in a fluidized bed opposedjet mill Part 1: Batch grinding kinetcs. Powder Techn. 1999. N 106, P. 78-87.
4. Heng P.W.S., Chan L.W., Lee C.C. Ultrafine grinding using a fluidized bed opposedjet mill: effect of process parameters on the size distribution of milled particles. S.T. P. PharmaSciences. 2000. N 10, P 445451.
5. Vogel L., Peukert W. Modelling of grinding in an Air Classifier Mill Based on a Fundamental Material Function. KONA. 2003. N 21. P 109-120
6. Agapov Yu.N., Lukyanenko V.I., Bardin A.P., Stogney V.G. The influence of the height of the fluidized bed on the parameters of the gas distribution grid. Bulletin of VSTU. Series: Energy. vol. 2. No. 6. 2006.S. 139-143.
7. Agapov D.Yu., Nadeev A.A., Rodionov K.N., Stogney V.G., Kleinikov I.Y. Experimental comparison of gas distribution devices for the formation of a centrifugal layer. Bulletin of the Voronezh State Technical University. V. 8. N. 7-1. 2012. P. 118-121.
8. Kalashnikova N.V., Volzhentsev A.V. Improving the technological process of drying wheat grain with the justification of the parameters of the fluidized bed dryer. Bulletin of agricultural science. V. 16. N. 1. 2009. P. 44-45.
9. Volzhentsev A.V. Determination of the quality of grain liquefaction in dryers. Education, science and production. N. 4 (9). 2014.P. 4-7.
10. Gvozdev D.V., Salnikov A.A., Fomichev A.G. Calculation and design of apparatuses with a fluidized bed of granular material. Gas distribution grilles. Izv. MVO USSR, Series Chemistry and Chemistry. No. 4. 1963. p. 38.
11. Greek F.Z. On the geometric characteristics of gas distribution grilles. Chem. prom 1966. p. 67-68
12. Afanasyeva O. V., Mingaleeva G. R. Apparatus for producing activated carbon on a mini-thermal power plant. News of Tomsk Polytechnic University. Engineering of geo-resources. T. 330. No. 9, 2019. P. 42-50.
13. Todes O. M. Devices with a fluidized granular layer: (Hydraulic and thermal fundamentals of work) L.: Chemistry: Leningrad. Separation. 1981. 296 p.
14. Yakovlev V. A., Sakhibgareev A. F. Aerodynamic design features of mixture-preparation systems for slotted gas-turbine units with forced air supply. Academy, N. 7 (34). 2018.P. 23-25.
15. Beloglazov I.I., Sharikov Y.V., Beloglazov I.N. Study of the distribution of the residence time of firing products in a fluidized bed furnace. Zapiski Mining Institute, V. 169. 2006. P. 61-62.