Математическое моделирование разделения тонкодисперсной фазы в гидроциклоне

ПРОЦЕССЫ И АППАРАТЫ ХИМИЧЕСКОЙ ТЕХНОЛОГИИ

УДК 66.063.62

А. Г. Лаптев, М. М. Башаров, А. И. Фарахова МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАЗДЕЛЕНИЯ ТОНКОДИСПЕРСНОЙ ФАЗЫ

В ГИДРОЦИКЛОНЕ

Ключевые слова: разделение дисперсий, турбулентная миграция, коэффициент переноса, эффективность турбулентной

сепарации.

Рассмотрено турбулентное движение жидкости с дисперсной фазой в гидроциклонах. Используется подход, когда турбулентное осаждение мелкодисперсной фазы представляется как разновидность диффузионного процесса. Получены выражения для определения скорости переноса частиц. Представлены результаты расчетов эффективности турбулентной сепарации гидроциклонов.

Keywords: separation of the dispersions, turbulent migration, transfer coefficients, turbulent separation efficiency.

The turbulent movement of the liquid with dispersed phase in hydrocyclone has been considered. Applied approach assumes presentation of turbulent sedimentation of the fine phase as a variety of diffusion process. Formulae for particle fluence rate determination were received. The results of calculations of the turbulent separation efficiency in hydrocyclone were presented.

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

П! = 1“П

I

где П| - эффективность сепарации за счет I - го механизма.

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

В данной работе предложено определять эффективность турбулентной сепарации по модели идеального вытеснения:

П = 1-exp(- N),

PdF

где N =----- - число единиц переноса; - коэф-

О

фициент переноса дисперсной фазы, м/с; Р - поверхность стенок гидроциклона, м2; О -

производительность гидроциклона, м3/с.

Выражение (2) хорошо согласуется с энергетической теорией пылеулавливания и вероятност-

но-стохастической моделью турбулентной миграции частиц [1, 2].

Для определения коэффициента переноса вс рассмотрим теорию турбулентной миграции

частиц [1 - 4].

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

масштабом несущих их пульсационных вихрей; полидисперсность частиц учитывается пофракционно; при концентрации С <0,2 кг/м3 частицы не соударяются и не коагулируются друг с другом.

В работе [3] сделана оценка размеров частиц, увлекаемых турбулентными пульсациями среды. Весь спектр осаждающихся частиц разделяется на три основные группы:

- I группа - частицы, полностью увлекаемые турбулентными пульсациями среды. Их диаметр должен удовлетворять условию:

d4 < 0,134^ж/рчи* ;

ч ’ а/ мЖ/кч

- II группа - частицы, обладающие некоторой инерционностью по отношению к увлечению турбулентными пульсациями:

0,134^ж/рчи* < с1ч < 13,4^ж/рчы* ;

- III группа - частицы, не увлекаемые турбулентными пульсациями среды, должна удовлетворять условию:

d4 > 13,4^RM*/p4u. ;

^ч " "'’~гм' '^ж/^ч“*>

где Р - радиус контактного устройства (гидроциклона), м; ^ж - вязкость жидкости, Па-с; Рч -

плотность частицы, кг/м3; и* - динамическая скорость трения на стенке, м/с.

Динамическую скорость можно вычислить, используя среднюю скорость диссипации энергии £ по следующим формулам [5]:

^,25

u* =

11,6X£v

- ApSu є =

рж

вх

V.

ж

где X = 0,4 - константа Прандтля; V- кинематическая вязкость жидкости, м2/с; рж - плотность жидкости, кг/м3; АР - гидравлическое сопротивление гидроциклона, Н/м2; Б - площадь входного патрубка, м2; ывх - скорость потока во входном патрубке, м/с; Уж - объем жидкости в гидроциклоне, м3.

Известно, что перенос вещества может происходить как молекулярным путем, так и вследствие турбулентного обмена. Переносом с помощью молекулярного механизма в гидроциклоне можно пренебречь. Для доказательства последнего утверждения сделана оценка коэффициента броуновской диффузии с использованием данных для гидроциклона с диаметром 0,036 м [6] и размером частиц Ьч =2,98-10-6 м.

Коэффициент броуновской диффузии по формуле Стокса-Эйнштейна:

J6p

кбТ

1,38 10 23 • 300

6п^г 6• 3,14 15,910-3 • 2,9810-6 -14,

= 0,463•10_| м2/с

где к б - константа Больцмана, Дж/К; Т - абсолютная температура, К; r = d^2 = 2,98 -10

-6,

м -

радиус частицы, м.

Число Шмидта:

1-6

SC бр == 18,1-10 = 39.1-108

бР D6p 0,463-10-14

Толщина подслоя, в котором преобладает молекулярный перенос:

5бр =

11,6v

11,6 • 18,1-10

,-6

Sc 2/3 u Sc 2/3 s 8

SC6p u*SC6p 0,569 • |39,1-108

= 1,49 •10-10 м,

2/3

где 61 - толщина вязкого подслоя, м.

Диаметр минимальных улавливаемых частиц в гидроциклоне диаметром 0,036м составляет 5,95-10-6 м. Таком образом, диаметр минимальной частицы многократно превышает подслой, в котором происходит

молекулярный перенос. Следовательно, молекулярной диффузией частиц можно пренебречь.

Мерой интенсивности осаждения частиц из турбулентного потока на стенках является скорость турбулентного осаждения частиц, под которой подразумевается количество частиц, осаждающихся из дисперсного потока на единицу поверхности стенок за единицу времени, отнесенной к концентрации частиц: = J/c , где j - удельный поток частиц к

стенке, кг/(м2/с); с - средняя концентрация частиц за пределами пристенного слоя, кг/м3. Так, например, в теории турбулентной миграции аэрозолей Pd = ut - скорость турбулентного осаждения, м/с.

Сопротивление переносу частиц запишем, используя двухслойную модель Прандтля [5]:

1 51 dy y” dy

— = J — + J — Pd 0 v б1 Dd

где у - поперечная координата к стенке гидроциклона, м; - коэффициент турбулентной диффу-

D

зии

2 — Т

частиц, м /с. Известна связь йс =----------, где

1+ Ш^Тр

йТ - коэффициент турбулентной диффузии среды без дисперсной фазы, м2/с; ы^ - угловая частота

и*

энергоёмких пульсаций: ы^ = 0 1р ; Тр - время

РчСи ^

релаксации частицы: Тр = —. Произведение

18Мж

ы^Тр для частиц первой группы очень мало и не учитывается, так как в этом случае ы е т р << 1 [3].

Коэффициент турбулентной диффузии для гидроциклона получен в виде [6] йу = 0,0112У^Г,

где у - окружная скорость жидкости, м/с; г - радиус.

После интегрирования (11) получим коэффициент переноса частиц:

Pd =

R-| (l + wETp )n(R^R1)'l

u*

JE'p 0,0112V

0,0112 V

tu*/('

t

1 + WET

E'pJ

; ,

Р10,112У^ (| + ы^Тр)+ и*!п(да/Р1) где = и*^V - безразмерная координата, так

как при расчете учитывается турбулентный перенос по всему радиусу гидроциклона Р . Согласно двуслойной модели Прандтля: р = 11,6.

Выражение (12) справедливо для частиц I и II группы.

Частицы III группы не увлекаются турбулентными пульсациями среды, и эффективность их удаления вычисляется по другим выражениям.

5

В частном случае с учетом толщины слоя осевших частиц на стенке гидроциклона уравнение (12) получит вид (при R1 > RQC ):

В = 0,0112У^и^

d (R1 - Roc )0,112Vt/(l + ыЕт p )+ u*ln(Rœ/R1 ) ’ где Roc = u*ho^v - безразмерная толщина слоя осевшей дисперсной фазы; hoc - толщина осевшего слоя частиц, м.

Из выражения (13) следует, что шероховатость поверхности гидроциклона, вызванная осевшей дисперсной фазой, повышает значение коэффициента (13) и эффективность разделения (2).

Выполнен расчет эффективности гидроциклона для частиц d4 < 6,0 мкм. Для примера, приведенные в работе [7], получено [8]: перепад давления Др =130783 Н/м2; динамическая скорость и* =0,475 м/с; диссипация энергии 533130 Вт/м3;

тангенциальная скорость Vt =7,1 м/с; безразмерная координата R œ =472; коэффициент переноса дисперсной фазы рd =0,0148 м/с; число единиц переноса N =0,544; эффективность разделения П =0,419. По экспериментальным данным [7] частицы диаметром 5,26 мкм и плотностью рч =2800кг/м3 удаляются из

нефти с эффективностью 42%. Таким образом, подтверждена адекватность рассмотренной математической модели и полученного уравнения (12).

Разработанная математическая модель используется в расчетах промышленных гидроциклонов нефтегазохимического комплекса [9].

Литература

1. Медников, Е.П. Турбулентный перенос и осаждение

аэрозолей. / Е.П. Медников. - М.: Наука, 1980. -176 с.

2. Сугак, Е.В. Очистка газовых выбросов в аппаратах с интенсивными гидродинамическими режимами. . -2-е изд. / Е.В. Сугак, Н.А. Войнов, Н.А. Николаев. - Казань: «Отечество», 2009. -224 с.

3.Лаптев А.Г., БашаровМ.М., Фарахова А.И. Эффективность турбулентной сепарации мелкодисперсной фазы в тонкослойных отстойниках // Энергосбережение и водоподготовка, №5 (73), 2011.- С. 43-46.

4. Лаптев А.Г., БашаровМ.М., Исхаков А.Р. Модель переноса мелкодисперсной фазы в турбулентных газовых средах // Вестник Казан. гос. технол. ун-та. - 2012. -№14. - С.96-99.

5. Лаптев, А.Г. Модели пограничного слоя и расчет тепломассообменных процессов. / А.Г. Лаптев. - Казань: Изд-во Казанск. ун-та, 2007. -500 с.

6. Адельшин, А.Б. Энергия потока в процессах интенсификации очистки нефтесодержащих сточных вод. Ч.1: Гидроциклоны. / А.Б. Адельшин. - Казань: КГАСА, 1996. -200с.

7. Мустафаев, А.М., Гидроциклоны в нефтедобывающей промышленности. / А.М. Мустафаев, Б.М. Гутман. - М.: Недра, 1981. -260с.

8. Башаров, М.М. Устройство и расчет гидроциклонов: учеб. пособие. / М.М. Башаров, О.А. Сергеева; под общ. ред. А.Г. Лаптева. - Казань: Вестфалина, 2012. - 92с.

9. Лаптев А.Г., Фарахов М.И. Энергоэффективное оборудование разделения и очистки веществ в химической технологии // Вестник Казан. гос. технол. ун-та. - 2011. №9. С. 152-158

© А. Г. Лаптев - д-р техн. наук, проф., зав. каф. технологии воды и топлива КГЭУ, tvt_kgeu@mail.ru; М. М. Башаров -канд. техн. наук, дир. по техническому обслуживанию и инжинирингу ОАО «ТАНЕКО»; А. И. Фарахова - асп. каф. технологии воды и топлива КГЭУ.