Особенности движения радиальных потоков в центробежных экстракторах безнапорного типа

УДК 66.061.35

А. А. Салин, Н. С. Гришин, С. И. Поникаров

ОСОБЕННОСТИ ДВИЖЕНИЯ РАДИАЛЬНЫХ ПОТОКОВ В ЦЕНТРОБЕЖНЫХ ЭКСТРАКТОРАХ БЕЗНАПОРНОГО ТИПА

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

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

Keywords: centrifugal extractor, radial two-phase flows, the opening angle of the nozzle.

The relationship between the velocity vectors of individual elementary volumes of fluid in free-flow centrifugal machines was determined. An continuity equation was discussed with regard to the unidirectional radial flow description in a centrifugal extractor. The features of liquid flows in the plano-cylinder rotor with diff erent types of nozzles, working in the field of centrifugal forces were considered, using the analogy with nozzles for the gravitational field.

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

Центробежные силы, необходимые для развития большой поверхности межфазного взаимодействия и механического разделения потоков могут генерироваться не только путем подвода дополнительной внешней механической энергии, но и за счет особенностей движения самих взаимодействующих потоков [1, 2].

Для установления зависимости между векторами скорости отдельных элементарных объемов потока жидкости в центробежных аппаратах безнапорного типа можно использовать цилиндрическую систему координат. Выделим внутри потока жидкости в роторе элементарный объем, ограниченный тремя парами смежных координатных поверхностей (рис. 1).

Рис. 1 - Изменение массы элементарного объема жидкости в цилиндрической системе координат

Примем, что составляющие скорости потока вдоль осей координат - Уг, Уф, Уг, а плотность движущейся жидкости - р. Тогда через грань

[г-йф-йг] выделенного участка за единицу времени т вдоль оси г войдет масса жидкости, равная Мг:

Мг = Р  VI-  Г  Ср С—  с1т .

На противоположной грани выделенного объема скорость и плотность жидкости изменятся, и

дVГ

будут соответственно равны VI +---Сг и

dr

dp

Р +--Сг . Тогда через эту грань [(г+ёг)-ёа-ёг] за

дг

то же время йт выйдет масса жидкости равная Мг+йг: Г д(р- VI) Л

М|+Сг = Р- V +——- Сг

V д )

Пренебрегая членом высшего порядка малости, приращение массы в элементарном объеме вдоль оси г составит:

Л. \/

d| = МГ -V^:

p- Vr d(p- Vr) +—

r

- dr- r - dp- dz- dT.

По аналогии находим приращение массы в выделенном объеме вдоль осей ф и г:

dMp =-

dMz =-

d(p - V

dp d

dz

dp- dr - dz - dT;

- dz - r - dp- dr - dT .

Полное приращение массы в элементарном объеме за время йт будет равно сумме найденных приращений:

VV ао> и) 1 ^ V дР  V—^ —г +-- +---— +-—

dM=-

v r dr r dp dz J

-r - dp- dz- dT - dr.

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

dM= I p + — I  r  d^-dz dT  dr- p  r  d^-dz dT  dr =

d ^ .

Эр

=— r  d^ dz  dT  dr Эт

dp

Для установившегося движения — = 0,

Эт

масса элементарного объема жидкости не изменяется, следовательно:

р  Vr д

+-— +---—+-

r

r dr r d^ dz

Для капельных жидкостей, которые практически несжимаемы (р = const) уравнение неразрывности можно записать в виде:

Vr dVr 1 d

+ -r dr

+ - r

= 0.

Эй

Рассмотрим уравнение сплошности для однонаправленного радиального движения жидкости (вдоль оси г) в центробежном экстракторе (рис. 2).

Рис. 2 - Схема движения жидкости в роторе центробежного экстрактора

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

V. с1Уг

— + —- = 0 . (10) г Сг

Разделив переменные и интегрируя, получим: + 1пг + С = 0 , откуда V. • г = С 2.

о С2

Из рисунка видно, что (V)r=R =

nRH R

следовательно C2 =-

Q

Полученное уравнение

2-п-Н

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

Движение потоков в плоскоцилиндрическом роторе центробежного аппарата (рис. 3) можно рассматривать аналогично движению потоков в диффузоре гравитационного насадка [3-

5], который имеет относительно большие гидравлические потери (в среднем коэффициенты Цс= ц =0,45 и ^ = 4).

Рис. 3 - Движение радиальных потоков в центробежном экстракторе

Рассмотрим несколько случаев движения потоков в плоскоцилиндрическом роторе с различными видами насадок, работающих в условиях наложения на поток поля центробежных сил, интенсивность воздействия которого переменна по радиусу. Различного вида насадки можно получить, изменяя угол раствора а1 от 0° до 180°: от расширяющегося типа при а1 = 0°, когда проходное сечение насадка растёт с увеличением радиуса и до случая истечения жидкостей в осевом направлении через плоскую цилиндрическую щель

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

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

а = а' + а" + а",

где а' - конический сходящийся насадок для гравитационного поля; а" - цилиндрический насадок для центробежного поля, обеспечивающий равенство сечений канала для

плоскоцилиндрического ротора; а'" - угол, обеспечивающий соответствие профиля потока в насадке (сопле).

Для гравитационного поля установлено [6, 7], что наиболее благоприятным для

гидродинамического режима истечения жидкостей является конический сходящийся насадок с углом конусности а' = 13°, коэффициент скорости у которого равен 0,97, а коэффициент расхода 0,95.

Однако оптимальный угол раствора для гравитационного насадка (конфузора) не будет соответствовать радиальному потоку в насадке центробежного экстрактора. Это объясняется тем, что наличие центростремительного ускорения с увеличением радиуса способствует росту скорости потока, тогда, согласно условию неразрывности потока сечение его с изменением радиуса будет также меняться. Следовательно, оптимальный угол будет несколько отличаться от значения величины 13° в зависимости от числа оборотов ротора и скоростей фаз.

Для определения второй составляющей угла раствора а'', будем исходить из условия равенства сечений потока в плоскостях Я0 и Я1 (рис. 3).

Согласно данному условию, принимаем:

2 • п • R0 • hQ = 2 • п • R1 • h1:

откуда

Г R ^

или

h0 = h1

tgal =

R

v 0 у

2 •

Подставляя значение h0, из уравнения (3),

получим:

tgal =

Г Rj^ h' ; 1 - h'

v R0 jy

h'

2-(R- R0) 2 • R0

Согласно имеющимся данным [8-13], в центробежном поле при истечении жидкостей через цилиндрические сопла и распылители замечено, что при длине более 3 мм происходит отрыв струи от стенок сопла, т.е. при дальнейшем движении жидкость не заполняет сечение сопла полностью. Это происходит за счет изменения (возрастания) центробежной силы, действующей на текущую в канале сопла жидкость. Следовательно, к условию цилиндрического насадка необходимо добавить условие, чтобы угол конусности соответствовал профилю потока жидкости в насадке. Для упрощения считаем, что сечение потока уменьшается пропорционально увеличению скорости. В данном случае потери скорости на трение при движении по конусу минимальны и в расчете ими пренебрегаем. Согласно этому запишем из условия неразрывности:

UR0- fR0 = UR2- fR2 = UR1- fR1 = const, (16) где f" - сечение струи; f - сечение канала из условия равенства сечений канала; UR0 и URi - скорость жидкости на радиусах R0 И R:.

Скорость в сечении R0 можно взять как отношение расхода тяжелой жидкости к площади

сечения Я0, тогда скорость в сечении Я1 определим из следующих соображений. Поскольку нас интересует качественная картина истечения, а не количественная его сторона, то рассмотрим два случая:

1) При однофазном истечении жидкости из цилиндрического насадка в центробежном поле, как отмечено выше, происходит отрыв струи от стенок и после 2-3 мм происходит свободное течение с ускорением ш2Я, то есть сечение струи будет изменяться пропорционально значению радиуса:

U

f'= f т i = f0

R0

Ui uuc • Ri

Mhi= h0

R0 ^ UR0

R ' U

Ri

где

f0 = 2п • R0 • h0 !

f¡"= 2п • Rj • h¡

h¡ - высота

потока жидкости на R¡ радиусе.

Когда насадок выполнен в виде плоскопараллельного ротора, то к выражению (17) добавится геометрический фактор постоянство сечения по направлению радиуса:

2-п-R= 2-n-R|-h1 = 2-п-R • =cons.

Отсюда следует, чтобы выполнить это условие «цилиндричности» необходимо изменять высоту проходного сечения по зависимости:

hí =

R0 • h0

R

В итоге, мы имеем две составляющие угла раствора (а1), которые будут определяться из соотношения:

(В V и >

(hi)1 = h' + h¡¡= h0

0

v Ri yv

1+-

U

U

Ri У

откуда

tga1

= 0.5 • h

Г ^ V

1 -

v Ri yv

U

1 +

R0

URi У

0

Ri - R0

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

h

tga=-

0

ГR V'

1-

0

U

Л

R0

v Ri A URi У

+ 0.23,

где 1да' = 1д13 = 0.23 .

2) Двухфазное истечение жидкостей. Рассмотрим плоскоцилиндрический ротор дифференциально-контактного экстрактора как насадок с противоточным движением в нем двух фаз. Учитывая, что скорость движения дисперсной фазы на порядок больше по величине сплошной фазы, то принимаем допущение, что уменьшение скорости и2 незначительно, подтверждаемое как экспериментальными, так и литературными данными [12].

Отсюда можно предположить, что незанятая струей полость насадка будет просто

заполнена другой фазой. Следовательно, выводы по первому случаю можно использовать и для второго, внося, соответственно, поправку на скорость в сечении Я1 в уравнении (22), получим:

К

tga = -

'О

(Ч- *О>

1-

^О V р1)

1+-

и.

Л"

иЬоп - W )

+ 0.23'

Данное выражение (23) и будет гипотетичным предложением изменения высоты проходного сечения внутрироторного пространства для радиальных двухфазных потоков, с точки зрения рациональной формы гидравлического насадка. Однако это условие не является оптимальным с точки зрения массоотдачи и эффективности работы центробежного экстрактора.

Таким образом, становится возможным определить оптимальный угол насадки исходя из гидравлических условий движения потоков в проходном сечении центробежного аппарата безнапорного типа. Дальнейшие исследования в данной области предполагают оценку влияния изменения геометрии насадка на процесс массообмена.

Литература

1. И.Р. Калимуллин, А. В. Дмитриев, А. Н. Николаев, Вестник Казанского технологического университета, 15, 11, 153-155 (2012).

2. Л.Н. Москалев, С. И. Поникаров, И. И. Поникаров, Вестник Казанского технологического университета, 14, 14, 235-241 (2011).

3. Б. Гебрахт, Свободно-конвективные течения, тепло- и массообмен: пер. с англ. Мир, Москва, 1991. 678 с.

4. Е.В. Аметистов, Тепло-массообмен: Справочник. Энергоиздат, Москва, 1982. 512 с.

5. Н.З. Френкель, Гидравлика. Госэнергоатомидат, Москва, 1956. 456 с.

6. И.Е. Идельчик, Аэродинамика технологических аппаратов. Машиностроение, Москва, 1983. 351 с.

7. М.Е. Дейч, А.Е. Зарянкин, Гидрогазодинамика: учебное пособие для вузов. Энергоатомиздат, Москва, 1984. 384 с.

8. Гришин, Н.С. Исследование гидродинамических процессов в центробежном экстракторе и иксообразной насадкой: дис. ... канд. техн. наук / Н.С.Гришин. - 1975.

- 152 с.

9. Бочкарев, В.Г. Экспериментальное исследование течения жидкостей через отверстия в контактных элементах центробежных экстракционных аппаратов: дис. ... канд. техн. наук / В.Г. Бочкарев. - Казань, 1969.

- 146 с.

10. Ластовцев, А.М. Пропускная способность вращающихся распылителей / А.М, Ластовцев // Тр. МИХМ. - Вып. 11. - 1957. - С. 71-82.

11. Ластовцев, А.М. Уравнение дробления жидкости вращающимися распылителями / А.М, Ластовцев // Тр. МИХМ. - Вып. 13. - 1957. - С. 29-42

12. Шарифуллин, Р.Г. Исследование процесса истечения жидкостей из насадков в поле действия центробежных сил: дис. ... канд. техн. наук / Р.Г. Шарифуллин. -Казань, 1981. - 112 с.

13. Ш.А. Шибазов, В.И. Ашихмин, Изв. вузов. Сер. «Нефть и Газ», 9, 61-64 (1972).

© А. А. Салин - асс. каф. машин и аппаратов химических производств КНИТУ, c888aa@mail.ru; Н. С. Гришин - д-р техн. наук, проф. той же кафедры; С. И. Поникаров - д-р техн. наук, проф., зав. каф. машин и аппаратов химических производств КНИТУ, mahp_kstu@mail.ru.