Численное исследование процесса дегазации в гидроциклонах Текст научной статьи по специальности «Физика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2012 Математика и механика № 4(20)

УДК 621.928.37

О.В. Матвиенко, М.В. Агафонцева

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ПРОЦЕССА ДЕГАЗАЦИИ В ГИДРОЦИКЛОНАХ

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

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

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

Одной из главных причин многочисленных сбоев и проблем в работе систем теплоснабжения является проблема качества воды, с которой в систему поступают газы. Присутствие их в теплоносителе может вызывать большое количество проблем в системах отопления - коррозию, грязь, шум, проблемы циркуляции, ухудшение теплопередачи и т.д. Примерно 22,1 мл/л воздуха входит в систему в растворенном состоянии [1]. Присутствие азота вызывает образование пузырей и пробок, в то время как присутствие кислорода и углекислого газа - коррозию. Свободный газ в виде пузырьков может существенно нарушить циркуляцию. В результате этого происходит снижение производительности насосов или их поломка.

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

- снижение скорости потока,

- увеличение скорости подъема пузырьков,

- центробежный эффект,

Теоретические основы гидромеханики гидроциклонов и связанных с ней процессов разделения развиваются уже несколько десятилетий. С основными этапами этого развития можно ознакомиться в [2, 3]. На начальном этапе становления теории ограничивались подбором на основе экспериментов эмпирических формул для расчета сепарационных характеристик гидроциклона [3] или при более или менее достоверных данных о характере течений в аппарате стремились к выводу достаточно простых по форме аналитических выражений. Как правило, при этом рассматривалось движение одиночных частиц в детерминированной постановке [4, 5]. Работы в этом направлении выявили ряд особенностей движения частиц в

гидроциклоне и позволили оценить в целом влияние некоторых параметров аппарата и смеси на сепарационные характеристики. В исследованиях [6-9] предложен механизм сепарации твердой фазы в аппарате, в основе которого лежит представление о турбулентной диффузии частиц, вызывающей поток твердых частиц, противоположный седиментационному, обусловленному действием центробежных сил. Математические модели, ориентированные на широкое применение современной компьютерной техники, даны в [10, 11]. Здесь обсуждаются многие важные детали гидромеханики жидкости и движения одиночных частиц в аппарате. Отметим, что подавляющее большинство публикаций относится к изучению процессов разделения твердых частиц. Исследованию движения воздушной фазы в гидроциклонах посвящено заметно меньшее количество работ [11, 12]. В основном они относятся исследованию движения жидкости в гидроциклонах с низким содержанием воздушной фазы.

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

Математическая модель процесса гидроциклонирования

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

Расчетные параметры аппарата имели значения, соответствующие экспериментам (рис. 2.) [10]:

ё„ = 75 мм, = 25 мм,

= 25 мм,

= 12,5 мм,

Рис. 1. Схема течения в гидроциклоне

Ь1 = 75 мм, Ь2 = 200 мм, Ьъ = 25 мм, 11 = 100 мм,

12 = 50 мм.

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

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

Для описания гидродинамики и процессов переноса в гидроциклоне использовалась физико-математическая модель, которая включает:

- двумерные уравнения Навье - Стокса, осредненные по Рейнольдсу для описания гидродинамики суспензии (уравнения На-вье - Стокса записаны с использованием цилиндрических координат, которые наилучшим образом подходят для описания осесимметричного режима течения);

- модификацию к - е-модели турбулентности, учитывающую неизотропность турбулентности и влияние центробежных сил на процессы генерации/диссипации турбулентности;

- модель дрейфа пузырьков с учетом их турбулентной диффузии.

¿ої

[Г

¿иї

Рис. 2. Внешний вид и схема гидроциклона

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

(1)

дри 1 друг -*— +—-— = 0;

дх г дг

дри 1 дриуг

дх

дг

др д

дх дх

1 д

г дг

Неїї I 2

( ди

ди

дх

ду

Неїїг I ^ + ТТ \дг дх

ди 1 дуг

+

дх г дг

(2)

дриу + 1 дру г

дх

дг

др д

дг дх

+-—

г дг

Неїї г

,ду 2 ( ди 1 дуг

2--------1 — +----------

дг 3 V дх г дг

ду ди

+

дх дг

Неїї ^ + —

(3)

дрuw + 1 друwг = д дх г дг дх

дрик + 1 друкг д дх г дг дх

дw 1 д 1 ї д ( w ^ руw

Не^ _ дх _ г 2 дг _^Гф дг V г ) г

Неї дк + 1 д Неї г дк

ак дх г дг 1 к 0 1

дх

+

г

друєг д Неїї дє 1 д Неїї дє

дг дх _ СТє дх. г дг _ СТє дг _

+ О _ рє ;

+ (О_ Є2рє));

(4)

(5)

(6)

а

с

I

Ь

а

2

Ь

2

Ь

3

(dv Л2 I v + idr +

fdu dvЛ (dwЛ f dw/r' \dr dx) v dx) f dr ,

В к -£ -модели турбулентности коэффициент турбулентной вязкости | определяется как I = С^рк 2е-1.

Значения констант выбираются в соответствии с рекомендациями [13]:

С = 1,44, С2 = 1,92 (1 - C3Ri), С3 = 0,001, ak = 1, стЕ = 1,3,

r<^>

= 2,5.

Ri =

k2 w d(wr)

2 2 Є r

dr

C = 0,09.

I-1

Для определения скорости дрейфа пузырьков относительно несущей фазы предполагалось равновесие между действующими на пузырек массовыми силами и силой сопротивления [14]. В соответствии с принципом динамического баланса сил скорость движения дисперсной фазы относительно несущей жидкости может быть определена как

{рь -рИд) ^ ^ =[„К)!

Kd = 4d2b- с r

3 №liqCd Rerel

vbwb

ab =

(7)

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

CD =

16

Re

rel

1+

16

Re

rel

1 (1 9 V

+ 3,315Re

rel

-1/2'I

Rerel = pliq Kell db / H-liq .

Баланс массы дисперсной фазы определялся с помощью уравнения турбулентной дисперсии, которое описывает конвективный перенос пузырьков осреднен-ным потоком и стохастическое движение пузырьков вследствие турбулентных пульсаций (турбулентную диффузию).

др (u + urel )Mj 1 dp(v + vrel) rMj

dx

+

r

dr

д

dx

dM,

PDp' —

+11 r dr

dM, rpDp'—

(8)

Коэффициент турбулентной диффузии дисперсной фазы определяется следующим образом:

Dr, =±

Pb

где |at - коэффициент турбулентной вязкости несущей фазы, у - параметр, характеризующий инерционность дисперсной фазы. Для низких значений турбулентного числа Рейнольдса Ret = [p4kdb) /) < 1 параметр инерционности рассчитывается следующим образом:

у = 1 -

1 db 2а L

3 -exp I -а

db

1-exp I -а

db

где а — 18—Яе, 1 - безразмерный параметр, обратный времени релаксации пу-

Рь

зырьков, Ь — -

к

3/2

- масштаб турбулентности.

В переходной области (1 < Яе1 < 103) можно использовать зависимость

=1+^76

р Ь

1 1

1

где

ї - (бЯе-2/3 +1)ехрІ —р —\-1, ї2 = бКе-2/3

ї3 — б

1

— +1 2Яе -1 Яе, ^3 +— І ехр б І 1 І 1 б) *\ Яе, рьаь І б

Яеі Рь аь 12 р Ь

-1,5

-2/3

12 Р — . ,

ехр І------- ------\-1,

Яеі Рь аь

ї4 - б

1 +1 2Яе, - ІЯе;2/3 +1 1 ехр (— б ^ ‘ І ‘ б\ ЯЯЄі Рь ¿ь І б

-2/3

-1.

В случае высоких значений турбулентного числа Рейнольдса (е, > 1000) параметр инерционности вычисляется как

V-1 -1 [ ТІ1п

Рь { ¿1

1 + б-Р|—1 +18 (Р—

Рь ¿ь

При расчете движения совокупности пузырьков предполагалось, что их начальное распределение по размеру подчиняется функции плотности распределения Розина - Раммлера - Шперлинга - Боннета (КЯ8Б) [1б]. Перейдем от непрерывного распределения пузырьков по размерам к дискретному. Для этого выберем число рассматриваемых фракций Ы, характерный диаметр пузырька, характеризующий эту фракцию 5г0 и ширину каждой фракции: 5-0 <5г0 < 5-0.

Массовую концентрацию пузырьков каждой фракции можно определить следующим образом:

и

| 50 ехр

М (5г0)-■

Я \т-1

1 а 50

5*

0ехр

Я \т-1

1 а-0

5*

(9)

-1/П ■

где т - параметр распределения, 5* - характерный размер пузырька, связанный

с медианным диаметром 5м следующим соотношением: 5* =(1п2)

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

є

йь =80|^

1/3у

(10)

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

( Д Л NМ,

■ -ТМ,

Р =

і=1

+ ^

і=1 Ръ

Зависимость вязкости среды от концентрации воздушной фазы может быть учтена, например, с помощью формулы Томаса [17].

р N ( р Л2 N ( _ N Л

: л* , г\ г/

к

N Г п V N

= 1 + 2,5Р-£М} + 0,0275[РМ £М2

Н-г Ръ і=1 \Ръ / і=1

• ехр

1Рг_ ръ 1=1

16,6-^- £ М}

Эффективная вязкость (цей-) определяется как сумма молекулярной вязкости среды (цт) и турбулентной вязкости (^ ).

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

Таким образом, граничные условия на входе формулируются в виде

Уб1Г

а

к = Ти • м-

к

3/2

8 = -

М1 = І

4пКСКт Р 4пКСКтР Р' П^С

Здесь у = 0,15, п = 0,005, Ти = 0,03 - константы модели, 0>т - массовый расход, ЯС = ВС /2 - радиус гидроциклона, Я1п - радиус подводящего патрубка, 3п

- площадь сечения подводящего патрубка.

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

г = 0, V = 0

ди

дг

дк_

дг

=0

д8

дг

=0

ддМі

дг

= 0 .

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

V = 0, и = 0, м = 0

кпм =

N С

8_, =

к 3/2с 3/4 кпм

М

дп

= 0.

V =

Т

Здесь к = 0,4 - константа фон Кармана, % - напряжения на стенке, гпм - расстояние от стенки гидроциклона до ближайшего пристеночного узла, отмеченного индексом - пм.

На выходе из гидроциклона (как в верхнем, так и в нижнем сливе) осевые составляющие градиента тангенциальной скорости, а также турбулентных характеристик к и є предполагаются равными нулю. Значения радиальной скорости V в выходных сечениях берутся равными нулю. Давление р в верхнем сливе определяется из предположения о радиальном равновесии потока, в то время как давление в нижнем сливе предполагается равным атмосферному. Для осевой составляющей скорости в выходных сечениях граничные условия не формулируются, а задается распределение давления, которое и определяет неявным образом распределения осевой скорости. Таким образом, в выходных сечениях граничные условия можно записать в виде

V = 0 , — = РМ (верхний слив), р = рет! (нижний слив),

дг г

dw _ о дк _ о де _ ^ дм' дх дх дх дк

Анализ результатов

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

На основе этой математической модели было проведено численное исследование структуры течения в гидроциклоне. Расчеты проводились для пузырьков, распределенных в соответствии с функцией распределения Розина-Рамлера. Значения параметров распределения полагались следующими: dm = 0,3 мм и m = 2 .

Процедура решения описанной выше математической модели основывается на решении системы уравнений Рейнольдса в динамических переменных. Конечноразностный аналог системы дифференциальных уравнений получен интегрированием их внутри контрольного объема конечно-разностной сетки. В расчетах использовались смещенные сетки со 100 узлами в радиальном и 300 узлами в осевом направлении. При аппроксимации конвективных членов использовались разности против потока. При моделировании диффузионных членов использовалась степенная аппроксимация. Система конечноразностных уравнений является нелинейной, и для ее решения применялся итерационный метод. На каждой итерации применялась продольно-поперечная прогонка. Давление рассчитывалось с помощью итерационной процедуры SIMPLE.

Картина линий тока (рис. 3, а) упрощает анализ течения жидкости в азимутальной плоскости. Здесь можно обнаружить тороидальный вихрь. Этот вихрь находится в центре цилиндрической части гидроциклона и обеспечивает возврат по-

тока по стенке к зоне ввода, а также отвечает за транспорт жидкости в нижней части аппарата вниз по конической стенке в направлении нижнего слива и подъем жидкости в приосевой зоне по направлению к верхнему патрубку.

Расчеты показали, что в приосевой области в верхней части аппарата давление становится ниже давления окружающей среды, что находится в хорошем соответствии с экспериментальными данными [10]. Когда одно (или оба) сливных отверстия открыты, в эту зону происходит всасывание воздуха и формирование воздушного столба. В периферийной области давление резко увеличивается и достигает максимального значения в окрестности стенки. Изменение давления в осевом направлении мало по сравнению с изменением в радиальном направлении, поэтому изобары практически параллельны оси гидроциклона.

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

Рис. 3. Распределение основных параметров в гидроциклоне: а - линии тока, б - изолинии турбулентной кинетической энергии

На рис. 4 приведены экспериментальные и теоретические распределения осевой скорости в двух сечениях гидроциклона.

&

§

о

3

я

О

Радиальная координата г, м

Рис. 4. Радиальное распределение осевой скорости в гидроциклоне: 1, 3 - эксперимент; 2, 4 - расчет; 1, 2 - 60 мм от верхней крышки; 3, 4 - 175 мм

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

о

а

§

О

те

(Й

К

л

К

к

2

к

(Й

н

Радиальная координата г, м

Рис. 5. Радиальное распределение тангенциальной скорости в гидроциклоне: 1, 3 - эксперимент; 2, 4 - расчет; 1, 2 -40 мм от верхней крышки; 3, 4 - 60 мм

Анализ радиального распределения тангенциальной скорости (рис. 5) показал,что в ядре потока она плавно возрастает от нуля на оси течения до некоторого максимального значения в периферийной части потока, далее, вследствие прилипания к стенке, происходит падение V . Отметим, что в ядре потока радиальное распределение тангенциальной скорости можно моделировать зависимостью вида

V = агп, где а и п - некоторые константы.

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

Рис. 6. Изолинии концентрации пузырьков воздуха в гидроциклоне: а - ¿ь = 0,1 мм, шаг изолиний - 0,01 мг/кг; б - ¿ь = 0,5 мм, шаг изолиний - 0,1 мг/кг; в - ¿ь = 1 мм, шаг изолиний - 1 мг/кг

В результате концентрация воздушной фазы с ростом диаметра пузырьков в верхнем сливе увеличивается (рис. 6, б, в). В то же время доля пузырьков, покидающих гидроциклон через нижний слив, резко падает и для пузырьков с диаметром более ёъ = 0,1 мм становится практически равной нулю.

Рассчитывая потоки каждой фракции пузырьков воздуха через верхнее и нижнее сливные отверстия, можно получить кривую разделения (зависимость доли определённой фракции, попадающей в верхний слив от размера пузырьков этой фракции: q(dъ0) = Qoí(dъ0)/Qln(dъ0), где q(dъ0) - доля пузырьков заданного размера, покидающих гидроциклон через верхний слив, Qof(dъ0), Qln(dъ0) - массовые потоки пузырьков воздуха в верхнем сливе и подводящем патрубке. Результаты таких расчётов приведены на рис. 7. q

0,98 0,96 0,94 0,92 0,90

0,88

10-5 dъo 10-4

Рис. 7. Кривые разделения (выход в верхний слив) пузырьков воздуха различного начального диаметра

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

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

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

ЛИТЕРАТУРА

1. Gase in kleinen und mittleren WasserheiznetzenTechnische Universität Dresden, Institut für Energietechnik, koordinierter Schlussbericht, AiF Forschungsthema Nr. 11103 B, November 1998.

2. Svarovsky L. Hydrocyclones. London: Technomic Publishing Co., 1984.

3. ПоваровА.И. Гидроциклоны на обогатительных фабриках. М.: Недра, 1978.

4. Баранов Д.А., Кутепов А.М., Лагуткин М.Г. Расчет сепарационных процессов в гидроциклонах // ТОХТ. 1996. Т. 30. № 2. C. 117-122.

5. Кутепов А.М., Лагуткин М.Г., Баранов Д.А. Метод расчета показателей разделения суспензий в гидроциклонах // ТОХТ. 1994. Т. 28. № 3. С. 207-211.

6. Neeße Th. and Schubert H. Modellierung und verfahrenstechnische Dimensionierung der turbulenten Querstromklassierung. Th. 3. 1976. Chem. Techn. V. 28. No. 5. P. 273-277.

7. Schubert H. and Neesse Th. A hydrocyclone separation model in consideration of the turbulent multiphase flow // Proc. Int. Conf. on Hydrocyclones. Cambridge. 1980. P. 23-27.

8. Matvienko O., Dück J., and Neeße Th. Numerishe Simulation der Strömungen in einem Hydrozyklon // Gessellschaft für Angewandte Mathematik und Mechanik, Annual Meeting, Bremen, April 6-9, 1998, Book of Abstracts. S. 85.

9. Matvienko O., Dück J., and Neeße Th. Hydrodynamics and particle separation in the hydrocyclone // Proc. 2nd Int. Symp. on Two-Phase Flow Predictions and Experimentation. May 23-26, 1999. Pisa, Italy. V. 2. P. 923-929.

10. Monredon T.C., Hsien K.T., andRajamani R.K. Fluid flow model of the hydrocyclone: an investigation of device dimensions // Int. J. Mineral Process. 1992. V. 35. No. 1. P. 65-83.

11. Pericleous K.A. and Rhodes N. The hydrocyclone classifier — a numerical approach // Int. J. Mineral Process. 1986. V. 17. No. 1. P. 23-43.

12. Pericleous K.A.Mathematical simulation of hydrocyclones // Applied Mathematical Modelling. 1987. V. 11. Issue 4, P. 242-255.

13. Гупта А., Лилли Д., Сайред Н. Закрученные потоки. М.: Мир, 1987.

14. Boysan, F., Ayers W.H., Swithrnbank J. A fundamental mathematical modelling approach to cyclone design // Trans. Inst. of Chemical Engineers. 1982. V. 60. P. 222-230.

15. Tryggvason G., Bunner B., Esmaelli, et. al. A front-tracking method for the computations of multiphase flow // J. Comput. Phys. 2001. P. 169, 708-759.

16. Островский Г.М. Прикладная механика неоднородных сред. СПб.: Наука, 2000. 359 с.

17. Dueck J., Matvienko O., and Neesse Th. Numerical modelling of hydrocyclone dynamics for process control // Advances in Filtration and Separation Technology, Science & Technology of Filtration and Separations for the 21st Century / Eds. Shiao-Hung Chiang and Samuel E. Lee, American Filtration&Separation Society. V.15. Pittsburgh, 2001.

Статья поступила 19.12.2011 г.

Matvienko O. V., AGAFONTSEVA M. V. NUMERICAL SIMULATION OF THE DEGASSING PROCESS IN HYDROCYCLONES. In this paper, the degassing process in a hydrocyclone was investigated. The possibility to apply the hydrocycloning method for air phase removing from a liquid is shown.

Keywords: hydrocyclone, degassing, turbulence, mathematical modeling,

MATVIENKO Oleg Viktorovich (Tomsk State University)

E-mail: [email protected]

AGAFONTSEVA Margarita Vladimirovna (Tomsk State University) E-mail:[email protected]