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

2014 Математика и механика № 2(28)

УДК 532.517.3

А.В. Шваб, М.Ю. Попп

МОДЕЛИРОВАНИЕ ЛАМИНАРНОГО ЗАКРУЧЕННОГО ТЕЧЕНИЯ В ВИХРЕВОЙ КАМЕРЕ1

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

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

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

Физическая и математическая постановка

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

1 Работа поддержана грантом РФФИ 13-08-00367 А.

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

E

f

( В2

В

I/ //////////////У/У ///////.

3

B

Рис. 1. Расчетная область

ZJ\

2

1

0

Вз

A

A

Таким образом, аэродинамику закрученного течения в зоне сепарации (рис. 1, область ББ1Б2Бз), показанной на рис. 1, будем проводить с учетом возможного вращения всех стенок вихревой камеры. Движение закрученного потока несущей среды описывается уравнениями Навье - Стокса, которые с учетом осевой симметрии в цилиндрической системе координат в безразмерном виде можно представить как

U. + d(u2) , d(uU) =_dP + J_

dz dr Re

дт

dr

д \ dr2

д 2u„

1 du„

du_L + d(UrUz) + d(uz2) =_dp + j_

дт dr dz dz Re

dz

d 2Uz dr2

2 r dr

d 2u

„ +1 д dz2 r dr

(1)

(2)

duL ldiuruy) ld{uzuy) = dp + J_ dr dr dz dr Re

d 2u„

dr 2

d 2 u,„

dz 2

+1 d^L-u r dr r

2uruL

(3)

dur duz ur —- +—- + — = 0 . dr dz r

(4)

Здесь безразмерная форма уравнений получена путём введения масштабов длины R (радиус вихревой камеры) и скорости U0 (среднерасходное значение скорости на входе в вихревую камеру), т = tU0/R - безразмерное время, Re = UR/v - число Рейнольдса, h = H/R, где H - высота вихревой камеры. В силу небольших скоростей, плотность газа считается постоянной и безразмерная форма давления определяется с помощью постоянной плотности и квадрата скорости U0.

Для получения единственного решения ставятся следующие граничные условия. На стенках камеры радиальные и аксиальные составляющие вектора скорости равны нулю. На оси симметрии задаются условия Неймана для осевой составляющей скорости и равенство нулю для радиальной компоненты скорости. Определение граничных условий для окружной компоненты скорости на входе в вих-

ревую камеру и на вращающихся стенках даёт два дополнительных критерия: = Я/и0 и ЯЛ = Я/и0. Здесь - средняя угловая скорость вращения газа на входе в вихревую камеру и - угловая скорость стенок зоны сепарации. Таким образом, закрученное течение в вихревой камере определяется тремя критериями: Яе, и ЯЛ. Последние два критерия, по сути, есть обратные числа Россби.

Решение в переменных функция тока, вихрь, окружная скорость

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

1 дТ 1 дТ

иг =--; иг =---; (5)

г дг г дг

диг диг

ю = —г----. (6)

дг дг

После подстановки составляющих вектора скорости (5) в зависимость определения вихря (6) получим уравнение Пуассона для определения функции тока Т:

д2 Т д2 Т 1 дТ

2 ■ = юг +--. (7)

дг2 дг2 г дг

Эллиптическое уравнение (7) можно перевести в параболическое относительно времени, так как рассматривается стационарное решение задачи. В результате получим

дТ

5ij

Г д2 т д2 Тл

v дг2 дг2 ,

1 дТ

= -ю г--— . (8)

г дг

Здесь XI - безразмерное время, которое фактически является итерационным параметром при численном решении уравнения (8).

Перекрестным дифференцированием уравнений (1) и (2) с последующим вычитанием второго из первого получим уравнение переноса вихря ю:

дю д(игю) д(игю) = 1 (д2ю 1 д2ю 1 1 дю ю ^ 1 2иф диф

-+ 4 ' ' +-дт дг да Re

- + —^ + -

дг2 да2 г дг

2

г

+ (9)

г да

Таким образом, решение задачи сводится к решению дифференциальных уравнений (3), (8), (9). Для получения единственного решения задаются следующее граничные условия. На входе в вихревую камеру: функция тока определяется интегрированием осевой скорости по зависимости (5); вихрь равняется нулю в силу постоянства осевой составляющей скорости; окружная скорость определяется формулой мф = г Rg. На оси симметрии: у = 0, ю = 0, иф = 0. На стенках вихревой камеры имеем: у = const, вихрь определяется формулой Тома [4], иф = г-Rd. На выходе из вихревой камеры для всех искомых функций используется условие Неймана д/дг = 0.

Решение в физических переменных скорость - давление

Для проверки достоверности полученных результатов поставленная задача решается также в переменных «скорость - давление». Решение безразмерной системы уравнений (1) - (4) проводится методом физического расщепления полей

ип+1 - и

скорости и давления с использованием разнесенной разностной сетки [5]. Представим кратко методику расчета на основе этого метода. Уравнение переноса импульса символически в векторном и безразмерном виде можно записать как

— + (и-У)и —- У2и = -Ур. (10)

дт Яе

При применении метода физического расщепления полей скорости и давления задача разбивается на два этапа:

+ (н-У)и, = — V2и* - Урп; (11)

5т v 7 Яе

Й ■ = -У(5р). (12)

5т

Здесь сумма уравнений (11) и (12) аппроксимирует уравнение (10), символом 5 обозначается приращение функции, а звездочкой обозначается промежуточный вектор скорости. Умножая на градиент уравнение (12) и учитывая выполнение уравнения неразрывности для вектора скорости на п+1 временном слое, получим уравнение Пуассона для определения поправки к давлению 5р:

V2 (5р. (13)

5т

После расчета уравнения (11) определяется значение промежуточного вектора скорости, затем из уравнения Пуассона (13) определяется поправка к давлению. Зная величину вектора скорости на промежуточном временном слое и поправку к давлению, можно перейти к расчету скоростей и давления на п+1 временном слое с помощью зависимости (12):

ип+1 = и* -У(5р)5т; рп+1 = рп +5р.

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

Численный метод решения

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

дФ

-= Лг (Ф) + Лг (Ф) + Ш ;

дт

д дФ д , д дФ д , ч

Л г (Ф) = — (А —) --(иг Ф); Л г (Ф) = - (А —) --(иг Ф).

дг дг дг дг дг дг

В результате расщепления задача сводится к решению системы уравнений:

^ - 1 К (5ф*) = /п ;

5х 2 г

5Ф** 1 5Ф*

---Лг (5Ф ) =-;

5т 2 5т

фп+1 =фп +5Ф**,

где /п = Лг (Фп) +Лг(Фп) + Ш .

Результаты численных расчетов

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

0

0,2 0,4 0,6 0,8

0

0,1

0,2

Рис. 2. Сравнение результатов численного решения двумя подходами: вихрь -функция тока и скорость - давление при параметрах закрученного течения Яе = 10, = 1, Яа = 1, к = 1,5

Численные исследования поля окружной скорости показали, что вблизи стенок вихревой камеры образуются существенные градиенты окружной скорости при отсутствии вращения стенок, что может стать причиной попадания крупных частиц в мелкий продукт процесса фракционного разделения. В случае же вращения стенок камеры поле окружной скорости существенно выравнивается не только вблизи стенок аппарата, но также во всей области (примерно при г > 0,4), где собственно и происходит сепарация крупных частиц. Для иллюстрации сказанного на рис. 3 показано распределение изолиний окружной скорости при вращающихся стенках камеры, а на рис. 4 - без вращения стенок при тех же параметрах потока. На рис. 5 показано распределение линий тока в вихревой камере.

г

г

г 1.8-

0.2 0.4 0.6 0.8 г

г 1.8

0.2 0.4 0.6 0.8 г

г 1.8

0.2 0.4 0.6 0.8 г

Рис. 3. Изолинии окружной Рис. 4. Изолинии окружной Рис. 5. Распределение ли-

составляющей скорости в составляющей скорости в ний тока в вихревой камере

вихревой камере при пара- вихревой камере при пара- при параметрах потока Яе =

метрах Яе = 10, Rg = 1, метрах Яе = 10, Rg = 1, = 10, Rg = 1, Яа = 1, к = 1,5

Яа = 1, к = 1,5 Яа = 0, к = 1,5

При уменьшении в два раза высоты вихревой камеры показаны аналогичные распределения изолиний окружной скорости при вращении стенок камеры (рис. 6) и без вращения стенок (рис. 7), а также распределение линий тока (рис. 8) при тех же параметрах закрученного течения.

0.2 0.4 0.6 0.8 г

0.2 0.4 0.6 0.8

0.8 г

Рис. 6. Изолинии окружной Рис. 7. Изолинии окружной Рис. 8. Распределение ли-

составляющей скорости в составляющей скорости в ний тока в вихревой каме-

вихревой камере. Яе = 10, вихревой камере. Яе = 10, ре. Яе = 10, = 0, Яа = 1,

= 1, Яа = 1, к = 0,75 = 1, Яа = 0, к = 0,75 к = 0,75

1.6

1.6

1.6

1.4

1.4

1.4

1.2

1.2

1.2

1

0.8

0.8

0.8

0.6

0.6

0.6

0.4

0.4

0.4

0.2

0.2

0.2

г

0.4

0.2

0.2 0

Г

На рис. 9 показано изменение осевой составляющей вектора скорости в сечениях 1, 2, 3 (рис. 9, а) вихревой камеры и в сечениях Е - Е, 4, 5 (рис. 9, б), местоположение сечений показано на рис. 1. Из рис. 9, а хорошо видно, что в области сепарации, т.е. примерно при г > 0,4 распределение осевой скорости существенно не меняется, а значение радиальной скорости существенно меньше осевой. Таким образом, добиваясь стабильности окружной составляющей скорости за счет вращения стенок камеры, получаем в области сепарации постоянное поле скорости, что несомненно, должно привести к повышению эффективности процесса фракционного разделения частиц.

Рис. 9. Изменение вертикальной составляющей вектора скорости: а - в средних сечениях расчетной области и б - в выходной части вихревой камеры. Re = 10, Rw = 1, Rd = 1, h = 1,5

Заключение

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

ЛИТЕРАТУРА

1. Ушаков С.Г., ЗверевН.И. Инерционная сепарация пыли. М.: Энергия. 1974. 166 с.

2. Патент РФ № 2407601, ПК В07В7/083. Способ воздушно-центробежной классификации порошков и устройство для его осуществления / П.Н. Зятиков, А.Т. Росляк, И.М. Васе-нин, А.В. Шваб, А.А. Демиденко, Ш.Р. Садретдинов. Опубл. Б.И. № 36, 27.12.2010.

3. Шваб А.В., Хайруллина В.Ю. Исследование закрученного турбулентного течения между вращающимися профилированными дисками // Теоретические основы химической технологии. 2011. Т. 45. № 5. С. 557-565.

4. Роуч П. Вычислительная гидромеханика. М.: Мир, 1977. 618 с.

5. Патанкар С., ПетуховБ.С. Численные методы решения задач теплообмена. М.: Энерго-атомиздат, 1984. 124 с.

6. Андерсон Д., Таннехилл Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. М.: Мир, 1990. Т. 2. 337 с.

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

Shvab A.V., Popp M.Yu. MODELING OF THE LAMINAR SWIRLING FLOW IN A VORTEX CHAMBER. In this paper, numerical modeling of the swirling laminar viscous flow in the vortex chamber used for the separation of fine powders into small and large fractions by a given particle size is considered. To enhance the efficiency of particle separation, it is proposed to additionally tighten the walls of the vortex chamber to equalize the circumferential velocity field and thereby to stabilize the balance of centrifugal and aerodynamic forces acting on the separated particles. In this formulation of the problem, the choice of geometrical parameters of the vortex chamber is of great importance. Reliability of the results and conclusions presented in the paper is based on test calculations, as well as on a comparison of the numerical solutions obtained in the vortex - stream function - circumferential velocity variables and in the speed-pressure variables.

Keywords: swirling flows, speed, pressure gradient, vorticity, stream function, laminar flow, vortex chamber.

SHVAB Alexander Veniaminovich (Doctor of Physics and Mathematics, Prof., Tomsk State University, Tomsk, Russian Federation ) E-mail: avshvab@inbox.ru

POPP Marina Yurievna (Tomsk State University, Tomsk, Russian Federation) E-mail: popp.marina@yandex.ru

REFERENCES

1. Ushakov S.G., Zverev N.I. Inertsionnaya separatsiya pyli. Moskow, Energiya Publ., 1974. 166 p.

2. Zyatikov P.N., Roslyak A.T., Vasenin I.M., Shvab A.V., Demidenko A.A., Sadretdinov Sh.R. Sposob vozdushno-tsentrobezhnoy klassifikatsii poroshkov i ustroystvo dlya ego osushchestvleniya. Patent RF № 2407601, PK V07V7/083.

3. Shvab A.V., Khayrullina V.Yu. Issledovanie zakruchennogo turbulentnogo techeniya mezhdu vrashchayushchimisya profilirovannymi diskami (2011) Teoreticheskie osnovy khimicheskoy tekhnologii. V.45. No. 5, pp. 557-565.

4. Roache Patrick J. Computational fluid dynamics. Albuquerque, Hermosa Publishers, 1976. 616 p.

5. Patankar S., Petukhov B.S. Chislennye metody resheniya zadach teploobmena. Moskow, En-ergoatomizdat Publ.,1984. 124 p.

6. Tannehill John C., Anderson Dale A., Pletcher Richard H. Computational fluid mechanics and heat transfer . Taylor and Francis, 1984. 792 p.