УДК622.794.7
Оптимизация геометрических параметров гидровихревого инерционного
стратификатора Вентури
В.Н.МАКАРОВ, АВ.УГОЛЬНИКОВИ, Н.В.МАКАРОВ
Уральский государственный горный университет, Екатеринбург, Россия
Применение наноразмерных частиц в качестве модифицирующих добавок открывает новые возможности в создании материалов с уникальными свойствами. В основе эффективного механизма качественного совершенствования структуры ВВП России лежит утилизация техногенных минеральных образований (ТМО) и производство высокотехнологичной продукции. Исследования показывают, что эффективность утилизации ограничена высокими требованиями к фракционному составу, медианным размерам и дисперсности частиц ТМО, несовершенством техники и технологии и их классификации. При разработке способов классификации по дисперсии медианных размеров микрочастиц ТМО в условиях вероятностного распределения физико-механических параметров исходного сырья необходимо учитывать жесткие требования к классификации. Приведенные в статье исследования основаны на положении о существенно большем влиянии инерционных сил на траекторию гидродинамически не установившегося движения дисперсной системы микрочастица -капля жидкости в процессе гидровихревой стратификации по отношению к аэродинамическим силам их перемещения в кипящем слое. Показано, что в диапазоне кинетической энергии поступательного движения капель жидкости, обеспечивающем преодоление аэродинамического барьера коагуляции гидрофобных частиц ТМО, минимальный диаметр поглощаемых микрочастиц в процессе гидровихревой коагуляции зависит только от величины угловой скорости вращения капель жидкости.
Получены уравнения для критериев Эйлера и Рейнольдса, их средних значений, времени релаксации капель жидкости с интегрированными в них микро- и наночастицами ТМО в зависимости от их медианного размера в процессе гидровихревой стратификации. Разработанная математическая модель гидровихревой инерционной стратификации в кипящем слое тонкодисперсных ТМО позволяет определять оптимальные геометрические параметры и энергетические характеристики стратификатора Вентури, его аэратора, положения приемных бункеров.
Проведенные испытания подтвердили возможность классификации тонкодисперсных отходов горнометаллургического производства в диапазоне медианных размеров (0,5-5)-10-6 м по фракциям с дисперсией не более 20 %.
Ключевые слова: утилизация; стратификация; циркуляционная гетерокоагуляция; критерии Рейнольда и Эйлера; гидрофобность; медианный размер; дисперсия; труба Вентури; силы инерции
Благодарность. Авторы признательны А.Е.Пелевину, доктору технических наук, профессору кафедры обогащения полезных ископаемых Уральского государственного горного университета за содействие в проведении аналитических исследований, С.Г.Хронусову, директору ООО «ЭМПИКО», за практические исследования и внедрение материалов статьи.
Как цитировать эту статью: Макаров В.Н. Оптимизация геометрических параметров гидровихревого инерционного стратификатора Вентури / В.Н.Макаров, А.В.Угольников, Н.В.Макаров // Записки Горного института. 2019. Т. 240. С. 638-648. DOI: 10.31897/РМ1.2019.6.638
Введение. Обеспечение конкурентоспособности горно-металургического комплекса (ГМК) РФ в глобальном экономическом пространстве невозможно без ускорения темпов реконструкции действующих и создания новых производств с учетом передовых достижений науки, внедрения современного высокотехнологического оборудования для производства изделий и материалов с уникальными функциональными свойствами. В значительной степени это касается технологических переделов подготовки, классификации исходного сырья, существенно влияющих на качество конечного продукта.
Наиболее эффективной стратегией улучшения структуры ВВП России является внедрение технологии утилизации, т.е. использования техногенных минеральных образований (ТМО) в производстве высокотехнологичной продукции.
Применение наноразмерных частиц в качестве модифицирующих добавок, а также в качестве индивидуального материала открывает новые возможности применения известных веществ. Так, использование нанопорошков для реализации нового комплекса функциональных свойств не имеет альтернативы при создании тугоплавких дисперсно-упрочненных композиционных материалов [2].
Возвратная наносодержащая глиноземная пыль в подобном случае представляет собой оборотный балласт, массовая доля которого достигает 7-14 % от общего количества получаемого глинозема. Учитывая, что годовой объем производства глинозема в РФ оценивается в 11,5 млн т, масса оборотной глиноземной пыли составляет значимую величину. Глиноземная пыль может быть использована в качестве легирующих элементов при производстве тугоплавких материалов и в электронной промышленности [3, 11, 16].
Многочисленные исследования показывают, что эффективность утилизации ограничена высокими требованиями к фракционному составу, медианным размерам и дисперсности частиц ТМО. В большинстве случаев требуемый размер микро- и наночастиц и их дисперсия находятся в диапазонах: dч = (0,1 -6)10 6 м; 3о = 0,2dm. Одним из сдерживающих факторов повышения эффективности утилизации мелкодисперсных ТМО является высокая энергоемкость и недостаточное совершенство техники и технологии их классификации [5, 10, 12].
Существенным недостатком используемых средств классификации, является недостаточная эффективность формирования узкого диапазона улавливаемых фракций частиц микро- и нано-размера [13-15].
Жесткие требования классификации по дисперсии медианных размеров тонкодисперсных ТМО приводят к необходимости поиска способов и технических средств, которые в условиях вероятностного распределения физико-механических, геометрических, кинематических параметров наночастиц могут их удовлетворять.
Постановка задачи. Для обеспечения качественным сырьем производства материалов с уникальными свойствами необходима технология, в которой управляющее внешнее воздействие на процесс классификации по дисперсии медианного размера будет автомодельно, т.е. независимо от вероятностных характеристик физико-механических свойств ТМО.
Статья посвящена разработке математической модели гидровихревой стратификации в кипящем слое микро- и наночастиц ТМО и устройства для ее реализации в виде гидровихревого стратификатора Вентури.
Задача заключается в разработке способа сепарации частиц ТМО размером (0,5-6)10-6 м с дисперсией медианных размеров ±5 % и эффективностью содержания заданной фракции в общей массе отбора по кривой Тромпа 95-97 %, а также создании устройств гидровихревой сепарации, при которой стратификация частиц в заданном диапазоне размеров гарантирована их коагуляцией вращающимися каплями жидкости. Селективность стратификации должна обеспечиваться полученными с помощью предложенной математической модели параметрами стратификатора Вентури, аэратора, гидровихревых форсунок, а также энергетическими характеристиками вращающихся капель жидкости и кипящего слоя мелкодисперсного ТМО.
Методология и объект исследования. Гидровихревая инерционная стратификация - процесс выделения частиц заданного медианного размера и дисперсии из ТМО, находящегося в состоянии кипящего слоя, путем формирования траекторий их движения к приемному бункеру под действием инерционных сил неустановившегося процесса гидровихревого взаимодействия с вращающимися каплями жидкости за период не более времени релаксации.
Для установления зависимости между энергией дисперсной системы капля жидкости - микрочастица ТМО и потоком энергии вертикального ее движения в кипящем слое и геометрическими параметрами стратификатора Вентури в процессе инерционной гидровихревой стратификации в фазе релаксации необходимо получить уравнения движения частиц в зависимости от действующих на них сил.
Используемый метод аддитивного аэрогидродинамического расчета гидровихревой стратификации отличается раздельным учетом влияния скорости движения частиц в кипящем слое и скорости инерционного неустановившегося гидровихревого движения дисперсной системы капля жидкости - микрочастица ТМО. Такой подход позволяет исследовать влияние кинематических и энергетических параметров вращающихся капель жидкости на траекторию движения частиц, определить зависимость эффективности сепарации по фракциям от геометрических и энергетических параметров стратификатора.
Методически решение задачи разделено на два этапа. На первом этапе изучено гидродинамическое взаимодействие капель жидкости и частиц ТМО при больших числах Рейнольдса
в режиме установившегося движения с использованием основных положений теорий подобия для получения зависимости коэффициента сопротивления ki от критериев Рейнольдса и Эйлера. Для построения адекватной математической модели процесса гидровихревой стратификации использованы уравнения установившегося стоксовского движения в процессе классической коагуляции в форме уравнения Буссинеска [1, 7].
На втором этапе установлены основные закономерности изменения полученных на первом этапе индикаторов подобия их корреляции с временем релаксации т и коэффициентом сопротивления ki в условиях неустановившегося движения. Рассмотрение неустановившегося движения капли жидкости с последующим осреднением аэрогидродинамических параметров ее движения по времени релаксации с учетом закономерности неустановившегося надстоксовского движения позволяет решить задачу с высокой точностью осреднения значений коэффициента сопротивления ki [7, 8]. Таким образом, динамически активный участок гидровихревой коагуляции в стра-тификаторе Вентури, характеризующийся неустановившимся надстоксовским режимом течения, в котором сила лобового сопротивления существенно зависит от критериев Рейнольдса и Эйлера, может быть представлен математической моделью с осредненными значениями критериев, описывающих данный процесс.
На базе известной модели кинетической гетерокоагуляции частицы пыли каплей жидкости при угловой скорости жидкости юж = 0 в статье [6] предложена графоаналитическая модель гидровихревой ортокинетической гетерокоагуляции, построена система уровней, описывающая физический процесс поглощения твердых частиц вращающимися каплями жидкости с участием дополнительной энергии присоединенного вихря.
С использованием предложенной графоаналитической модели гидровихревой инерционной ортокинетической гетерокоагуляции в статье [9] получено уравнение для минимального диаметра поглощаемой твердой компоненты в системе жидкое - твердое при вращении капли жидкости с угловой скоростью юж и эффективного краевого угла смачивания в зоне контакта жидкой и твердой фазы с учетом дополнительной энергии обусловленной вращением капли жидкости.
Для абсолютно гидрофобных частиц, каковыми являются микрочастицы ТМО с диаметром менее 6-10-6 м, математически и экспериментально подтверждена гипотеза о прямой корреляции минимального диаметра полностью поглощаемых частиц с угловой скоростью вращения капель жидкости при гидровихревой инерционной кинематической гетерокоагуляции. При этом за порогом аэродинамического барьера физико-механические и химические свойства частиц не влияют на процесс смачивания в диапазоне диаметров поглощаемых частиц, соответствующих вихревой и классической коагуляции: dчш min < dч < dHmin [9].
Полученные уравнения использованы в данной статье для создания математической модели гидровихревой стратификации частиц компонентов ТМО по фракциям с заданной дисперсией. Уравнение для минимального диаметра полностью поглощаемой гидрофобной наночастицы при гидровихревой коагуляции (при юж = 0 м) можно представить в виде
d чш min = 24 5- (cos eV 1 - K Ш ®ж - Sin 0K¿2 ®ж \ (1)
(P ч -P гГж
где рч, рг - плотность частицы и газа соответственно, кг/м3; Уж - скорость капли жидкости, м/с; 5ж-г - коэффициент поверхностного натяжения на границе раздела двух сред жидкость-газ, Дж/м2; 9 - краевой угол смачивания на границе раздела двух сред жидкость-газ, рад; Кш - коэффициент влияния угловой скорости вращения капли жидкости на минимальный диаметр поглощаемых частиц; рж - плотность капли жидкости, кг/м3,
K _ ЛРжd4min Sin¿4 e
ш 85ж-г cose •
По результатам экспериментальных исследований установлено, что эффективное влияние закрутки на коагуляцию наблюдается при относительной угловой скорости вращения капли жидкости Кш < ш2 < 0,3. Коэффициент вариаций медианного диаметра поглощаемых наночастиц компонентов ТМО от угловой скорости вращения капель жидкости можно представить в следующем виде
К* = 48-
5„
(Р ч -Р гГж
К^008 е(шж - КШ®ж)-Кш 81П ешж.
(2)
Из уравнения (2) следует, что угловая скорость вращения капель жидкости может быть эффективным управляющим параметром в процессе гидровихревой стратификации гидрофобных частиц ТМО. Отличительной особенностью гидровихревой стратификации является ее высокая чувствительность к дисперсии медианного размера микрочастиц, так как их разделение основано на гидровихревой коагуляции, притом что размер капель жидкости существенно больше поглощаемых ими микрочастиц.
Для практического решения задачи утилизации наночастиц ТМО в статье предложен гидровихревой стратификатор Вентури, обеспечивающий их разделение по фракциям с заданным медианным размером и дисперсией (рис.1).
Конструктивно гидровихревой стратификатор Вентури включает в себя пневмотранспортный трубопровод для вертикального перемещения частиц сыпучих материалов ТМО и устройство для гидровихревого разделения частиц по фракциям за счет инерционной гетерокоагуляции их вращающимися каплями жидкости, состоящее из трубы Вентури, по оси которой в критическом сечении установлен аэратор с вихревыми форсунками, а по периметру - приемный кольцевой бункер.
Гидровихревой стратификатор Вентури содержит загрузочный питатель 9, установленный над коллектором 1 стратификатора. В смесительной камере 11 установлены пористая газораспределительная перегородка 10 и патрубок 12 для подачи сжатого воздуха и формирования на входе в коллектор 1 кипящего слоя сыпучего материала ТМО. На входе в трубу Вентури 7 установлен хонейкомб 2 для выравнивания скорости движения частиц. По оси трубы Вентури установлен вращающийся аэратор с гидровихревыми форсунками 5, а по его периметру расположен коллектор стратификации 6 с приемным бункером 4 сбора частиц компонентов ТМО по фракциям с дренажным клапаном 3, на выходе сепаратора установлен бункер с дренажным клапаном 8 для отходов.
ТМО из загрузочного питателя 9 непрерывно направляют в смесительную камеру 11, ограниченную газораспределительной перегородкой 10. По патрубку 12 под слой ТМО подают сжатый газ. Энергией сжатого газа аэрируют сыпучий материал до псевдосжиженного состояния и подают через коллектор 1, выравнивающий хонейкомб 2 на вход в трубу Вентури 7. Закрученные капли жидкости из форсунок аэратора 5 смачивают частицы кипящего слоя. Кинетическая энергия и скорость вращения капель жидкости обеспечивают гарантированную коагуляцию частиц ТМО с заданным минимальным диаметром.
Диаметр dа аэратора и диаметр трубы Вентури dВ в критическом сечении, а также энергия капель жидкости выбирают таким образом, чтобы коагуляция капель жидкости и частиц сыпучего материала ТМО с заданным минимальным медианным диаметром происходила в зоне до 80 % от диаметра трубы Вентури. Таким образом, в кольцевой зоне 20 % диаметра трубы Вентури по ее периметру частицы компонента ТМО в кипящем слое поднимаются в бункер отходов, не смачиваясь в силу недостаточности энергии капель жидкости для преодоления аэродинамического энергетического барьера. Это необходимо для обеспечения гарантированной эффективности стратификации частиц по фракциям с заданной дисперсией за счет исключения из процесса стратификации зоны неустойчивой коагуляции вблизи входа в бункер стратификации. Положение коллектора стратификации 6 и бункера стратификации 4 по отношению к плоскости расположения форсунок гидровихревого аэратора 5 определяется в соответствии с траекториями движения капель жидкости с интегрированными в них частицами сыпучего материала ТМО
12
Рис. 1. Схема гидровихревого стратификатора Вентури
Fм
Fд
Fri <■■
в процессе инерционной гидровихревой гетерокоагуляции, рассчитанными по предложенной математической модели. Через бункер с дренажным клапаном 8 происходит удаление влаги и последующий сбор готового продукта.
В целях повышения эффективности стратификации за счет управления влиянием силы Магнуса на скорость вертикального перемещения частиц ТМО, поглощаемых вращающимися каплями жидкости в стратификаторе Вентури, использован вращающийся с угловой скоростью юа аэратор.
Траектория движения наночастиц ТМО определяется инерционным взаимодействием наночастиц, вращающихся капель жидкости и потока энергии кипящего слоя.
Результаты исследований. Для построения математической модели неустановившегося гидродинамического взаимодействия капли жидкости с частицами компонентов в условиях больших чисел Рейнольдса примем допущение о том, что в процессе движения капли жидкости она сохраняет сферическую форму dж той же плотности, что и жидкость капли, у которой аэродинамические характеристики движения в газовой среде соответствуют фактическим характеристикам движения капель при тех же числах Рейнольдса. Диаметр dж будем считать аэродинамическим диаметром капли [7].
Для обеспечения однозначности и определенности решения с учетом предложенной научной гипотезы примем условие, при котором энергии поступательного движения капель жидкости достаточно для преодоления аэродинамического энергетического барьера в области кипящего слоя не более 80 % от диаметра трубы Вентури, т.е. диаметра расположения входа в приемный бункер.
Для построения системы уравнений движения капли жидкости с интегрированной в нее на-ночастицей компонента ТМО (рис.2) введем понятия приведенных аэродинамического диаметра и плотности наночастиц и представим их в виде
Fn
FА
Рис.2. Схема движения капли жидкости с интегрированной в нее частицей компонента ТМО и действующих на нее сил в условиях гидровихревой стратификации: - ось стратификатора Вентури; - координаты нижней и верхней границ коллектора стратификации
1 - капля жидкости с интегрированной в нее частицей компонента ТМО;
2 - координата расположения гидровихревой форсунки аэратора;
3 - траектория движения капли жидкости с частицей компонента ТМО;
4 - плоскость кипящего слоя
d Е =
Рх =
й ж р ж + й чр ч
й 3 + й:
(3)
где йж - диаметр капли жидкости, м; йч - диаметр микрочастицы ТМО, м.
Уравнение движения /-частицы при ее поглощении каплей жидкости в проекции на ось 0г в плоскости расположения гидровихревых форсунок аэратора (в соответствии с классическим уравнением Ньютона), учитывая рис.2, запишем в виде [4]:
дV ■
дК" ~ (4)
ш,
д?
= К,.
где Кп - сила сопротивления движению частицы компоненты ТМО в гидровихревом стратифи-каторе,
рп = к, 4 й Х/Р Л2,;
к, - коэффициент сопротивления ,-й частицы; йх , - диаметр ,-й частицы, м; VZi - скорость ,-й частицы, м/с.
п
Уравнение движения i-й частицы в проекции на ось 0z имеет вид
дК2
= -РА, - Ра + Рдг - РМ, , (5)
где РА - сила Архимеда, направленная вниз и действующая на ,-ю частицу, являющуюся аналогом силы тяжести,
ра, =1 лd I, (Ре, -Р г к; 6
Рс, - сила сопротивления Стокса, обусловленная вязкостью воздуха и физическими свойствами компонентов исходного сырья,
рс, = ф ц глd ъУ*; Рд - сила давления сжатого газа, создающего кипящий слой,
Рд, = СсР К2 % dе2, ;
РМ, - сила Магнуса, направление которой определяется угловой скоростью вращения аэратора,
% з
РМ, = -РгdЕ,®Е,Уф, ;
о
Сс - коэффициент силы давления сжатого газа, создающего кипящий слой, действующий на ,-ю частицу; Кг, У^ - скорость сжатого газа, создающего кипящий слой, и вертикальная составляющая скорости ,-й частицы, м/с; g - ускорение свободного падения, м/с2; ф, - коэффициент формы частицы в законе Стокса; р - коэффициент динамической вязкости воздуха, кг/мс; - окружная скорость движения капли жидкости с интегрированной в нее /'-частицей в плоскости кипящего слоя.
Уравнение (4) идентично уравнению Буссинеска, описывающему гидродинамически неустановившийся режим инерционного движения капли жидкости до и после коагуляции наноча-стицы компонента ТМО [7]
дК„ К
Г1 _ Г1
дt х,,
(6)
где тч; - время релаксации ,-й частицы,
х . = XV = к ч, Е,
2 + 3ц dе2, (РЕ, -Рг). , 3 + 3Ц 18цг ;
ц = Цж.
Ц г
Уравнение Буссинеска позволяет определить время релаксации как капли жидкости, так и наночастицы и связывает его с коэффициентом сопротивления. При этом коэффициент сопротивления в уравнении Буссинеска соответствует к, в уравнении (4).
Экспериментально установлено, что сила сопротивления возрастает существенно нелинейно с ростом числа Рейнольдса на участке надстоксовского движения в отличии от линейного роста ее при числах Re < 1 и при одновременном уменьшении времени релаксации т, что существенно усложняет нахождение его фактической величины, препятствуя тем самым использованию классических уравнений стоксовского движения при коагуляции [8].
Поскольку в уравнении аэрогидродинамики капли жидкости в горизонтальной плоскости стратификатора Вентури (4) существенной переменной величиной является коэффициент сопротивления газовой среды движению капли жидкости к,, установим его зависимость от физических величин, характеризующих гидродинамический процесс инерционного движения капли жидко-
сти под действием ускорения с начальной скоростью ¥0 в форме безразмерного симплекса в условиях установившегося движения. Уравнение зависимости коэффициента к, от независимых переменных в виде безразмерной степенной зависимости запишем в виде
ki _ СцгрГд/V2 + 0,25dj sin 02Уd
fdV- р -р
ri г ж г
ЛФ
-t
р г
V7.
(7)
Для определения индикаторов подобия воспользуемся матрицей независимых размерностей, ранг которой равен трем [1]. С учетом однородной системы линейных уравнений, составленной из показателей степеней уравнения (7), матрица индикаторов подобия примет вид
г
Ц г р г ■V + 0,25юЖd2 Sin 0ю d ж -Vri р ж -р г V)
-t рг
п1 - 2 2 2 2 0 0
п 2 0 0 0 1 1 -2
П3 - 2 2 0 2 0 + 2
(8)
Раскрывая определители матрицы (8), получим три индикатора подобия физического процесса установившегося аэрогидродинамического движения капли жидкости в газовой среде при больших числах Рейнольдса
п _р2 Vr2 + 0,25шЖdjSin0co)dz2i _ 2 . (9)
2 - Re jrnl. (9)
Ц 2
п P -Vtd2i (р Ж -р г) . (10)
п 2 _rV2 _-TV2-= Eu i' (10)
р V0 р V0
Р Х ^ 2/ па2 /114
% 3 =-2- = Ке^0/ . (11)
^ 2
Таким образом, уравнение (7) в критериальной форме примет вид
к, = с(ке2югЕй,6(яе2о, )с. (12)
Численные значения коэффициента пропорциональности С и показатели степени а, Ъ, с в уравнении (12) определяем из уравнения установившегося движения капли жидкости в газовой среде
2
= Ей ,Р Л2 (13)
о? 4
С учетом уравнений (4, 13) выражение для коэффициента аэродинамического сопротивления получим в виде
к■ =2 2Еи' 2 . (14)
После соответствующих преобразований уравнений (12-14) относительно числа Рейнольдса получим уравнение для определения времени релаксации дисперсной системы капля жидкости -микрочастица ТМО в виде
2 d j (ра -рг)Re2 9 Eui-цг Reji
_„2i\j~ 2l ГТ/ ~0i СЧ
_-—--;п-. (15)
Учитывая, что критерии Рейнольдса и Эйлера функционально связаны между собой и непрерывно и существенно изменяются по длине инерционного пробега капли жидкости в условиях надстоксовского режима, рассмотрим возможность решения задачи посредством осреднения
кинематических параметров. Классическая теория гидродинамического движения в условиях установившегося режима позволяет в квадратурах получать выражения для времени релаксации капель жидкости и частиц компонентов в зависимости от кинематических параметров течения. В работе [8] получено выражение для изменения времени релаксации при больших числах Рей-нольдса путем осреднения их значений. С учетом указанных данных выражение для критерия Эйлера в условиях неустановившегося гидродинамического движения на этапе релаксации дисперсной системы капля жидкости - микрочастица ТМО можно представить в виде
Ей = (1 + 0Д5КеГ7). (16)
ке0,
Таким образом, с учетом данных, приведенных в работах [7, 8], среднее значение критерия Эйлера на участке релаксации дисперсной системы капля жидкости - микрочастица ТМО жидкости представим в виде
6 1„ ---- 0,687 |
Еиср, = — (1 + 0,07Яе0Г7). (17)
кеог-
Учитывая соотношение между критерием Рейнольдса и коэффициентом сопротивления движению дисперсной системы капля жидкости - микрочастица ТМО в газовой среде, выражение для среднего значения коэффициента сопротивления в уравнении (7) запишем в виде
кер, = ^ (1 + 0,07Ке0;687). (18)
Ке0 ,
После соответствующих преобразований получим выражение для среднего значения времени релаксации дисперсной системы капля жидкости - микрочастица ТМО
= 4ё2 (рв-рг)(1 + 0^01Ке°);687 )(3 + зр) = 2. (19)
Еср' ъ Яе0, (2 + 3Ц)ц г '
Из уравнения видно, что среднее время релаксации капель жидкости с интегрированными в них компонентами ТМО является функцией квадрата диаметра дисперсной системы капля жидкости - микрочастица ТМО, что позволяет использовать данный факт для разработки технологии эффективной классификации мелкодисперсных сыпучих ТМО.
Учитывая однородность уравнений поступательного и вращательного движения тел в широком диапазоне чисел Рейнольдса, полученное соотношение для средних значений коэффициента сопротивления и времени релаксации может быть применено к угловой скорости вращения капель жидкости с интегрированными в них частицами ТМО, обусловленными вращением аэратора. При этом скорость вертикального перемещения дисперсной системы капля жидкости - микрочастица ТМО в кипящем слое может быть представлено в виде
1 Сс Р гУг2 ё ъ - 1 ё 2 (Ръ -Р г ^ - 4-10-4 ®0 ®а Р г Ръ 1 +
¿а
у* = --3-^^, (20)
Ф^ г
где Ръ . ъ р
гж
Из анализа уравнения (20) видно, что при выборе направления вращения аэратора можно добиться существенного замедления скорости вертикального перемещения частиц ТМО, изменяющейся пропорционально квадрату их диаметра, что повышает эффективность стратификации.
Таким образом, предложенный вариант поэтапного осреднения значений коэффициента сопротивления движению дисперсной системы капля жидкости - микрочастица ТМО в газовой среде и времени ее релаксации позволяет с достаточной степенью точности использовать уравнения классической аэрогидродинамики установившегося движения капли жидкости в диапазоне
чисел Рейнольдса до 104 на длине их свободного инерционного пробега для построения математической модели инерционной гидровихревой стратификации.
В техническом задании на проектирование гидровихревого стратификатора Вентури, как правило, указана производительность утилизации по массе сыпучих материалов ТМО Q, т/ч.
Диаметр стратификатора Вентури йВ при заданной его производительности определяем по формуле
Л В = • (21)
где Кр - коэффициент плотности кипящего слоя частиц ТМО, Кр = рч / ркс; ркс - плотность частиц ТМО в состоянии кипящего слоя.
Скорость капель жидкости в процессе стратификации на выходе из гидровихревых форсунок аэратора, характеризующего энергетические параметры, определяем по формуле
у - 3,7/15^xose+ db-d^. (22)
М 2 PzdEmm XZ
Угловая скорость вращения капель жидкости и коэффициент вариации медианного диаметра поглощаемых микрочастиц ТМО, определяющие требуемую фракцию в процессе стратификации, находим по формулам (1, 2).
Координаты входного коллектора стратификатора z\, z2 приемного бункера при заданном медианном диаметре йт и его дисперсии от определяем по формулам:
-2А 2 / „ \„,2 с 1Г1-4.
CcP TVrdE max - 2 (Pz max -P г kdZmax - 5 ^O^ ®0 ®а P acPsLx d4max ^ax f! +
z -_3_-_^lL X • (23)
z! _ у Zcp max ' V /
10фЦг ln , 0 !
V65ж-г C0s ePzmaxdZmax
1 CcP rVr2 d Z min - 1 (PZ min -P г )gd Zmin - 5 '!0 _ ^ ®0 ®aP ж Pz!nin d 4min d Zmin f! + ^ |
z2 - --3---XZ ■ (24)
2 у Zcp min' V /
ФЦ г ln / 0 ! !
V65ж-г C0s ePZmindZmin
гДе dZmin - dZm - 3Zm • dZmax - dZm + 3Zm ■
На рис.3-6 представлены результаты экспериментальных исследований и их сравнения с расчетами по предложенной математической модели.
Данные экспериментальных исследований (рис.3) в сравнении с расчетом по предложенной математической модели подтверждают эффективность управления минимальным диаметром поглощаемых гидрофобных частиц ТМО с помощью угловой скорости вращения капель жидкости в процессе гидровихревой стратификации.
Результаты экспериментов и расчетов (рис.4, 5) показывают широкие возможности управления геометрическими параметрами и расположением входного коллектора приемного бункера с помощью угловой скорости вращения капель жидкости для обеспечения требуемой эффективности стратификации микрочастиц ТМО по фракциям.
Существенное изменение времени релаксации как по результатам эксперимента, так и по расчетным данным, приведенным на рис.6, позволяет добиваться высокой эффективности классификации при заданной дисперсии медианных размеров микрочастиц ТМО вследствие существенного влияния изменения номинального диаметра их на инерционные силы, определяющие траекторию движения дисперсной системы капля жидкости - микрочастица ТМО.
Результаты экспериментальных исследований гидровихревого стратификатора Вентури, проведенного в «СМК-ТЕСТ», подтвердили достаточную для инженерных расчетов сходимость с предложенной математической моделью.
z2■10-1, м 2 _
5 6 ^Юж, с 1
Рис.3. Зависимость медианного диаметра смачивания частиц компонента ТМО от угловой скорости вращения капель жидкости при гидровихревой стратификации 1 - уголь; 2 - оксид кремния; 3 - оксид алюминия
Рис.4. Зависимость координаты входного коллектора стратификации от медианного диаметра частиц стратификации 1 - Юа = 50 с-1; 2 - Юа = 100 с-1; 3 - Юа = 150 с-1
5
1
2
0
1
2
3
4
йт10-6, м
4
1
2
3
5
Дz■10-2, м 5
4
3
а
т, с 1
0,9 0,8 0,7 0,6 0,5
Рис.5. Зависимость ширины входного коллектора стратификатора от дисперсии медианного размера частиц стратификации 1 - Юа = 50 с-1; 2 - Юа = 100 с-1; 3 - Юа = 150 с-1
б
т, с 1
0,9
0,95
0,8
~т-Г
20
2
7
1 2 3 4 5 6 7 8 ат-10 , м
10-6, м
0
1
2
3
4
5
Рис.6. Зависимость среднего времени релаксации капли жидкости с интегрированной в нее: а - частицей оксида алюминия от критерия Рейнольдса и б - частицей компонента ТМО от медианного диаметра 1 - й5;3,5-10~4, м; 2 - й^^-Ю-4, м; 3 - ¿^■Ю"4, м; 4 - уголь; 5 - оксид кремния; 6 - оксид алюминия; Яе0 = 40
Заключение. Результаты выполненной работы позволяют сделать научные и практические выводы:
1. Диаметр полностью поглощаемых твердых частиц зависит от угловой скорости вращения капель жидкости в процессе гидровихревой гетерокоагуляции, что позволяет использовать данный физический процесс для классификации микро- и наночастиц ТМО.
2. Время релаксации капель жидкости с интегрированными в них микро- и наночастицами ТМО в процессе гидровихревой стратификации, а следовательно, и инерционные силы неустановившегося взаимодействия с вращающимися каплями жидкости зависят от медианного диаметра, что является основным определяющим фактором достижения высокой эффективности классификации.
3. Геометрические параметры гидровихревого стратификатора Вентури определяются потребной производительностью и энергетическими характеристиками его аэратора.
4. Гидровихревая стратификация позволяет сепарировать микрочастицы ТМО в диапазоне 0,5-5-10-6 м с дисперсией медианных размеров не более 20 %.
ЛИТЕРАТУРА
1. Веников В.А. Теория подобия и моделирования. Применительно к задачам электроэнергетики. М.: Либроком, 2014. 439 с.
2. Влияние добавок легирующих керамических наночастиц на структурные параметры и свойства твердых сплавов / Ю.И.Гордеев, А.К.Абкарян, Г.М.Зеер, А.А.Лепешев // Вестник Сибирского государственного аэрокосмического университета им. академика М. Ф. Решетнева. 2013. № 3. С. 174-181.
3. Давыдов С.Я. Улавливание наноразмерной фракции частиц глиноземного производства / С.Я.Давыдов, Р.А.Апакашев, В.Н.Корюков // Новые огнеупоры. 2016. № 2. С. 12-15. DOI: 10.17073/1683-4518-2016-2-12-15
4. Лойцянский Л.Г. Механика жидкости и газа. М.: Дрофа, 2003. 840 с.
5. Макаров В.Н. Теоретические основы повышения эффективности вентиляции в технологических процессах на промышленных предприятиях / В.Н.Макаров, С.Я.Давыдов // Новые огнеупоры. 2015. № 2. C. 59-63. DOI: 10.17073/1683-45182015-2-59-63
6. Перспективный способ повышения эффективности высоконапорного гидрообеспыливания / В.Н.Макаров, Н.В.Макаров, В.В.Потапов, Э.М.Горшкова // Вестник Забайкальского государственного университета. 2018. Т. 24. № 5. С. 13-20. DOI: 10.21209/2227-9245-2018-24-5-13-20
7. ФроловА.В. Основы гидрообеспыливания / А.В.Фролов, В.А.Телегин, Ю.А.Сечкерев // Безопасность жизнедеятельности. 2007. № 10. С. 1-24.
8. Фукс Н.А. Механика аэрозолей. М.: Изд-во АН СССР, 1955. 352 с.
9. Эффективная локализация взрывов угольной пыли с использованием гидровихревой коагуляции / В.Н.Макаров, Н.П.Косарев, Н.В.Макаров, А.В.Угольников, А.В.Лифанов // Вестник Пермского национального исследовательского политехнического университета. Геология. Нефтегазовое и горное дело. 2018. Т. 18. № 2. С. 178-189. DOI: 10.15593/22249923/2018.4.7
10. Davydov S.Ya. Capturing Nanoparticles in Alumina Production / S.Ya.Davydov, R.A.Apakashev, V.N.Korukov // Refractories and Industrial Ceramics. 2016. Vol. 57. Iss. 1. Р.9-12. DOI: 10.1007/S11148-016-99-17-6
11. Davydov S.Ya. Utilization of Alumina Calcining Furnace Dust Containing Nanoparticles / S.Ya.Davydov, R.A.Apakashev, V.N.Korukov // Refractories and Industrial Ceramics. 2014. Vol. 55. Iss. 4. Р. 291-294. DOI: 10.1007/511148-014-97-11-2
12. Justification of environmental technologies and means for dust control of tailing dumps surfaces of hydrometallurgical production and concentrating plants / V.I.Lyashenko, A.A.Gurin, F.F.Topolniy, N.A.Taran // Metallurgical and mining industry. 2017. Iss. 4. Р. 8-17.
13. Makarov V.N. Theoretical basis for increasing ventilation efficiency in technological processes at industrial enterprises / V.N.Makarov, S.Ya.Davydov // Refractories and Industrial Ceramics. 2015. Vol. 56. Iss. 1. P. 103-106. DOI: 10.1007/s11148-015-9791-7
14. Novakovskiy N.S. Numerical simulation of shock-free strong compression of 1D gas layer's problem subject to conditions on characteristic / N.S.Novakovskiy, S.P.Bautin // Journal of Physics: Conference Series. 2017. Vol. 894. № 1. P. 1-8. DOI: 10.1088/1742-6596/894/1/012067
15. Numerical modeling of heat and mass transfer during hot and cool air mixing in a supply shaft in underground mine / N.I.Alymenko, А.А.Kamenskikh, A.V.Nikolaev, A.I.Petrov // Eurasian mining. 2016. № 2. P. 45-47. DOI: 10.17580/em.2016.02.11
16. Reverse circulation drilling method based on a supersonic nozzle for dust control / D.Wu, K.Yin, Q.Yin, X.Zhang, J.Cheng, D.Ge, P.Zhang // Applied Sciences. 2017. Vol. 7. Iss. 1. P. 5-20. DOI: 10.3390/APP7010005
Авторы: В.Н.Макаров, д-р техн. наук, профессор, uk.intelnedra@gmail.com (Уральский государственный горный университет, Екатеринбург, Россия), А.В.Угольников, канд. техн. наук, заведующий кафедрой, ugolnikov@yandex.ru (Уральский государственный горный университет, Екатеринбург, Россия), Н.В.Макаров, канд. техн. наук, заведующий кафедрой, gmf.gm@ursmu.ru (Уральский государственный горный университет, Екатеринбург, Россия).
Статья поступила в редакцию 03.06.2019.
Статья принята к публикации 26.07.2019.