CC BY
45
УДК 532.527
Р. Р. Усманова, Г. Е. Заиков, Р. Я. Дебердеев
МОДЕЛИРОВАНИЕ ЗАКРУЧЕННЫХ ТУРБУЛЕНТНЫХ ПОТОКОВ ПРИ ОПТИМИЗАЦИИ ПРОЦЕССА РАЗДЕЛЕНИЯ ДИСПЕРСНЫХ СИСТЕМ
Ключевые слова: численное моделирование; динамический газопромыватель; гидродинамика; Ansys CFX; компоненты
скорости; уравнение движения.
В данной работе предложен алгоритм моделирования процесса сепарации дисперсной фазы в газовом потоке с орошением жидкостью. Проведенные расчеты позволяют определять потенциальные возможности динамического газопромывателя при использовании его в качестве аппарата для очистки газовых выбросов. Верификация полученных расчетом данных, а также оценка параметров, определяющих возможность сепарации дисперсной фазы на капли орошающей жидкости, проводится путем моделирования процесса течения газожидкостного потока в пакете вычислительной гидродинамики Ansys CFX.
Keywords: Numerical Modelling; A Dynamic Spray Scrubber; Hydrodynamics; Ansys CFX;Blending Agents Of Speed; The Motion
Equation.
In the given work the algorithm of modelling of process of separation of a dispersoid in a gas stream with irrigation by a liquid is offered. The carried out calculations allow to define potential possibilities of a dynamic spray scrubber at its use in the capacity of the apparatus for clearing of gas emissions. Verification of the data had by calculation, and also an estimation of the parametres defining possibility of separation of a dispersoid on drops of an irrigation water, is spent by means of modelling of a gas-liquid stream in the program of computing hydrodynamics Ansys CFX.
1. Состояние вопроса, актуальность
С самого начала появления газопромывателей (а также всех прочих технических устройств) перед инженерами, которые их проектировали, стояла задача - спрогнозировать параметры работы создаваемой ими машины до того, как чертежи будут отданы в производство. Применительно к аппаратам ротационного действия задача осложнялась тем, что основные параметры их работы (сопротивление, эффективность, потребляемая мощность и др.) определяются картиной течения потоков внутри проточной части, а течение это по самой природе многофазной среды настолько сложно, что до сих пор нередко единственным надежным методом исследования в гидродинамике остается эксперимент. Лишь в последние годы был достигнут существенный прогресс в создании средств моделирования и расчета течений многофазных сред, позволяющих выполнять расчет со столь высокой достоверностью получаемых результатов, что необходимый объем эксперимента во многих случаях сводится до минимума. В отличие от эксперимента, численный подход дает возможность варьировать ряд важных параметров задачи, таких как вязкость, угловая скорость вращения модели, существенно влияющих на формирование и поведение вихревых течений.
Цели данной работы:
1. Разработать алгоритм моделирования процесса сепарации твердых частиц в динамическом газопромывателе, позволяющий определять потенциальные возможности аппаратов газоочистки.
2. Провести анализ численной схемы, построить численную модель в пакете вычислительной гидродинамики Ansys CFX. Выявить наличие вихревых движений и получить распределение скорости газодисперсного потока.
3. Для расчета движения частиц в закрученном потоке составить математическую модель процесса сепарации дисперсных частиц на капли орошающей жидкости. Сравнить результаты расчетов для различных режимов течения. Произвести верификацию полученных расчетом данных с данными, полученными в программе Ansys CFX.
2. Имитационное моделирование движения газодисперсных потоков
Для описания гидродинамики и процессов переноса используется физико-математическая модель [1], которая включает:
- уравнения Навье - Стокса, осредненные по Рейнольдсу, для описания гидродинамики и уравнения неразрывности;
—(rpvrvr)+—(rpvrv,)
dr dr
dvr Л д
dr { ^ dr j dz [ r^ dz
dvr VI dP
- -d(rPvrvT)+-d(rPv,vT) r dr dz ,
dv.
—I rr- I +--1 УУ-T —
dr [ dr I dz 1 T dz
dv,
PUU. . Pv ,
-л r 2 +
dr r r
PU, PV,Vr ---
r r
d~r (rPvrU )+dz (грул)
d ( dv Л d ( dv
dr [ Г dr ) dz dz divpv = 0
dP_
'dz.
где - скорость потока вдоль оси; иг- скорость потока в радиальном направлении; -
тангенциальная скорость потока; р - плотность смеси; коэффициент турбулентной вязкости; Р -давление; и - вектор скорости.
Для замыкания системы уравнений используется двухпараметрическая модель
д
д
турбулентности к - е [8], для неё решается два дополнительных уравнения переноса с целью определения к - турбулентной кинетической энергии и е - турбулентной энергии диссипации.
Модель турбулентности:
^Pckc_ а + а
dt
Эх.
Эис
= а-а Рс'£c +
. дх
Ърсес_ +РА ас =
Эх.
dt
Эх,
„ £с Эис £с Э = C —■- т—с— C —— +--
£ k 1 Эх, £ k Эх,
ц'с 1 Экс
а +t J К
а *+i I i
к2
* £
где кс- турбулентная кинетическая энергия газовой фазы; ок- турбулентное число Прандтля для уравнения кинетической энергии; цси -
молекулярная и турбулентная вязкость газовой фазы; £с - скорость диссипации турбулентной кинетической энергии; оЕ - турбулентное число Прандтля для уравнения диссипации кинетической энергии; ту -декартовы компоненты тензора напряжений: С^= 0.09, Cs 1 = 1.44, Се2 = 1.92, ак= 1.0, а= 1.3.
Расчеты показывают, что вблизи твердых стенок происходит весьма резкое изменение параметров к и s. Для надлежащего разрешения этих изменений приходится использовать весьма густую расчетную сетку.
Первый этап подготовки исходных данных для расчета течения - создание твердотельной геометрической модели, имитирующей объем, внутри которого происходит исследуемое течение [2]. Поскольку элементы проточной части динамического газопромывателя [3] (закручивающее устройство, подводящие патрубки, каплеотбойники) часто имеют весьма сложную форму, создание их твердотельной модели является неординарной задачей. 3D модель аппарата была выполнена в пакете твердотельного моделирования Solid Works (см. рис.1).
Рис. 1 - Геометрическая модель газопромывателя
В качестве начальных условий для расчета задается величина внешнего давления и температура окружающей среды. Величина внешнего (окружающего) давления составляет 1 атм., температура окружающего воздуха 250 С. Полагается, что на боковых стенках выполняется
условие прилипания итац=0, на верхней границе и в области стока задаются значения скорости:
U inlet = ^'aaiai
иах1а1 —щ
U radial = 0
Uswirl = 0
U inlet Щах1а
и , = ахга1 —и2
^radial = 0
и , = _ swirl 0
i + и„
J + Uswirlk
Исследуемые задачи решаются в осесимметричной постановке (не учитывается зависимость от азимутальной координаты ф), течение жидкости предполагается турбулентным и описывается системой управляющих уравнений в размерной формулировке. Решение системы уравнений совместно с граничными условиями было реализовано численным методом.
Рис. 2 - Проекции окружных скоростей в потоке
Из рисунка (см. рис.2) видно, что тангенциальная скорость возрастает от центра к периферии. Величина тангенциальной скорости, полученная численным решением, качественно согласуется с экспериментом. Профиль
тангенциальной скорости имеет форму параболы с максимумом, расположенным ближе к цилиндрической стенке, что можно объяснить действием центробежных сил. Экспериментальные результаты измеренных скоростей потока в осевом направлении показывают, что вблизи выхлопного патрубка в пристенном слое имеет место достаточно интенсивное течение. В центральной части потока осевая скорость имеет фактически постоянное значение. Численное моделирование по расчету осевых компонент скоростей потока показывает качественное согласие с опытными данными. Проведенное исследование работы газопромывателя позволят анализировать его работу с целью уменьшения энергозатрат при сохранении качества газоочистки. Разработанная модель помогает быстро и наглядно смоделировать движение запыленного газового потока с учетом внесенных в геометрию аппарата изменений. Таким образом, модель может
Э
применяться для оптимизации динамического газопромывателя.
конструкции
3. Численное моделирование движения газодисперсных потоков
Для расчета движения частиц в закрученном потоке была составлена математическая модель процесса сепарации дисперсных частиц на капли орошающей жидкости. Учитывались следующие силы, действующие на частицу: сила тяжести, сила Кориолиса, сила аэродинамического сопротивления среды и центробежная сила [5, 6].
Рис. 3 - Силы, действующие на частицу в потоке газа
Центробежные силы, возникающие при вращении ротора, обеспечивают дробление жидкости на мелкие капли, что обуславливает интенсивный контакт газов и улавливаемых частиц с жидкостью.
При разработке модели были приняты следующие упрощающие допущения о потоке и движении частиц [4]: течение воздуха
рассматривается в рамках модели вязкой несжимаемой жидкости; продольная (осевая) компонента усредненной скорости движения воздуха постоянна во времени и неизменна во всех точках потока; тангенциальная компонента усредненной скорости движения воздуха иф не зависит от продольной координаты; течение воздуха рассматривается как установившееся; усредненная скорость движения воздуха постоянна по времени и по сечению, эффект проскальзывания частиц относительно газа отсутствует; затуханием и переформированием тангенциальной скорости ифвдоль потока пренебрегаем; осевая скорость частицы равна осевой скорости воздуха в точке, в которой она находится; концентрация частиц пыли равномерна по времени и сечению.
Система уравнений движения частицы может быть записана в виде:
= 1 л.к и
(г
4 Р„ й„
■(Ух - их ) + и = ^ х (г
1Р.. к и
(г 4 Рч (ч а
■(У,-и,) +
со ■ г
(у
и, = т г = г-
^ = ^£ и отн .(Уг - иг)+ т2 г;
(г 4 рч (ч г
иг = (г
г (г
Аналогично записывается система уравнений движения капель жидкости:
= — . Ь.. к з .(У- W) + 8; Ъ = —
(г 4 рж (К отн ( х х) 8; 3 (г
(г
= —.Ь- ~ 3 . (У - ъ ) + 4 Рж (к 3отн У 3) +
т. г
3т =т. г = г—
у йг
(2)
(г
Л. Р.. к 3
4 Рж (к °
ТТ7 ч т. г о .(Уг - Ъг) —— + т г;
3 =
(г
Из системы уравнений движения частицы следует, что её траектория зависит от следующих факторов: (ч- диаметра частиц (капель (к); рч-плотности частицы (жидкости рж); лг- динамической вязкости газа (жидкости л ж); г- радиуса (характерного размера) аппарата; иго- начальной тангенциальной скорости газа, а также от геометрии завихрителя а, 6, 2.
Предложен подход, основывающийся на теории геометрического подобия пылеуловителей, и позволяющий учесть принудительную закрутку потока и геометрические параметры завихрителя.
В качестве критериев подобия предлагаются следующие безразмерные параметры:
относительный диаметр частиц (./Во, относительная плотность частиц ррг; число Рейнольдса Ке=уч(г/уг; число Струхаля БН= учитывающее угловую скорость ротора, характерный размер лопастного завихрителя В число определяющих параметров не включены, как не имеющие обычно существенного значения, ускорение свободного падения и параметры, характеризующие сжимаемость газа и тепловые процессы в нем. Запыленность газа в предположении независимого движения частиц пыли также не входит в число определяющих параметров.
Если ввести коэффициент (характерный размер лопастного завихрителя), то систему дифференциальных уравнений (1-2) можно перезаписать так:
(Ш, = (И,и _и |-■
Рч (ч I Уч
ёг
(и
(г (и
_ 18^г ГОп " Р.. (I I я1
.У'
- + 8
г В-,.1 -у,1+и±
г г
18«г ( О
= 3' г Г А ;£_ ^ ;В ;... 1 - Уг 1 - ^ + т* г;
(г Рч(ч2 [«О г I Я7 Я7 Я7 Я7 ) г I г
Введем безразмерные цилиндрические координаты: □=г/Яп, ф'=ф.
Реальные скорости заменим безразмерными:
и-=£ и-и-1^'=£ и. безразмерное время:
г
(
п,
1
г=%.,
я3
Для безразмерных величин сохраняются определения скорости:
Сделав соответствующие подстановки
их = ^, и■ К, и.= г, х йг' 9 йг' г йг'
обозначим: с = С = У ™ в
р г
я
После преобразования получим систему дифференциальных уравнений для безразмерных величин:
йи' ( С
= 18Ср \ с- - и\ йг' \ Г ч
йи',
= 18Ср (и' г (а;г' ; г' ; О';...)-V\) + ^
гч
-Л8Ср (и'г (а;г' ; Г ;Б';...)-W\)
(3)
йг' йи '
йг'
При начальных условиях:
Гч |г=0 = Г0' 2ч| г=0 = г0' 9 ч | г=0 = 0
-для цилиндрических координат,
и\[=0и' г0 , и;| г=0 = и' г0 , и'г| г=0 = ^0--для скорости частицы и двух безразмерных параметрах Ср и Сг решение этой системы определяет траекторию частицы в циклоне.
Система (3) является системой обыкновенных дифференциальных уравнений второго порядка. Зная радиальную Ог и аксиальную скорости их газового потока, её можно проинтегрировать и получить траекторию частицы.
Решение системы (3) при одинаковом виде функций и'ф(а; г'; г'; О';...) и и'г(а; г'; г'; Б';...), т.е. при одинаковой форме пылеуловителей, определяется безразмерными параметрами Ср и Сг, которые зависят от режимных и геометрических параметров. При равенстве этих параметров траектория частицы будет одинаковой. Таким образом, Ср и Сг являются критериями подобия для движения частиц пыли в газопромывателе. Для полного подобия движения частиц необходима также тождественность функций и'ф(а; г'; г'; Б';...) и и'г(а; г'; г'; Б';...) (геометрическое подобие).
Видно, что критерий Ср определяется физическими параметрами пылегазовой смеси и удельной производительностью газопромывателя. Критерий Сг определяется размерами конического завихрителя и углом наклона его лопастей, от которых зависит соотношение между тангенциальной и радиальной скоростями, а следовательно, соотношение сил инерции и силы сопротивления.
Уравнения движения, также как и уравнение Лапласа, решаются численно. Отметим, что система уравнений движения частицы является жесткой, так как содержит резко отличные значения производных
[9]. Матрица системы уравнений первого приближения для уравнений движения имеет вид:
(4)
V 18К, - — и 0 -18^ т 0
г г г2 г2
2и г -18К, 0 18К, V - ' йг и2 " г2 18К, V йг
0 0 -18К, 18К, ^ ' йг 18К, тг йг
0 1 0 0 0
0 0 1 0 0
Степень
её
обусловленности
приблизительно 10 , что говорит о жесткости системы уравнений движения. Использование явных методов решения системы (4) в результате возможной неустойчивости вычислительного процесса накладывает сильные ограничения на величину шага интегрирования. Решить эту задачу можно путем применения неявных методов [9]. При этом на каждой итерации по времени необходимо решать нелинейные уравнения, что существенно усложняет расчет. Поэтому было решено применять полунеявный метод решения системы уравнений движения. В используемом методе каждое расчетное разностное уравнение неявно по той переменной, по которой производится дифференцирование, при этом по остальным переменным уравнения остаются явными.
Для интегрирования уравнений движения частицы (3) были использованы два метода Гира: второго и первого порядка (который совпадает с неявным методом Эйлера). Если обозначить значение переменной на п-м шаге интегрирования как хп, а на (п+1 )-м как хп+1, то разностный аналог (4) при интегрировании полунеявным методом Эйлера примет вид (с учётом (4)):
и"+' - ип = 18К (К - ип+1 ^ - и+Г дг ' I гп I гп
V"+1 -V"
дг
W"+1 - W" дг
= 18К, V (г", г" )-V"+1 К^Тт" = 18К, (¿¥г (г", г" )-W"+1)
(5)
дг
- = V"
- = W"-
дг
"+1 -9" = и"+
Следует заметить, что применение полунеявного метода обеспечивает значительно более широкие границы устойчивости вычислительного процесса, чем явные методы (при сохранении простоты вычислений). Так верхняя грань устойчивого шага интегрирования по времени повышается на два порядка в случае использования метода Эйлера.
При расчете траектории частиц их скорость очень неравномерна, при этом шаг интегрирования выбирается автоматически в зависимости от скорости частицы. В соответствии с (5) можно записать следующую итерационную формулу для
и \ V' ч
г
г"+1- г"
численного решения уравнении движения и построения траекторий движения частиц пыли:
'un+1 = (un + 18KtAtn/rn )/(l + (l8Kt + Vn/rn )Atn) Vn =Vn + (l 8KtVz (rn, zn )+(un )2/ rn )Atn )/(l + 18Kt Atn) (6) Wn = (wn +18KtWz (rn, zn )Atn V(l + 18KtAtn)
,„n + 1
r-- = r» + Vn+1 Atn zn+1 = zn + Wn+1 Atn (pn+1 = pn + un+1 A tn I rn+1
Atn+1 = C VV^T^W")2
Картина движения газодисперсного потока, полученная путем решения уравнений (6) достаточно сложная. По сечению аппарата можно выделить несколько характерных областей (см. рис.4).
Рис. 4 - Характерные газодисперсного потока
зоны движения
В приосевой зоне возникает замкнутая циркуляционная зона 2, характеризующаяся интенсивным возвратным течением вдоль оси. Обнаружена крупномасштабная неоднородность течения по всему объему газопромывателя, которая хорошо регулируется путем изменения режимно -конструктивных параметров процесса. В конической части корпуса возникает тороидальный вихрь 4, смещенный к периферии. Газодисперсный поток проходит между циркуляционными контурами 2 и 4 и движется к выхлопному патрубку. В выхлопной зоне может быть реализован интенсивный сток очищенного газа. В зонах 3 и 4 линии тока замыкаются на стенке корпуса, что способствует интенсивной сепарации дисперсных частиц в зоне 1.
4. Верификация результатов расчета
Величина тангенциальной скорости, полученная численным решением, качественно согласуется с с результатами расчета в программном комплексе Ату8 С¥Х (см. рис.5). Профиль тангенциальной скорости имеет форму параболы с максимумом, расположенным ближе к цилиндрической стенке, что можно объяснить действием центробежных сил. Полученные на
модели результаты скоростей потока в осевом направлении показывают, что вблизи выхлопного патрубка в пристенном слое имеет место достаточно интенсивное течение. В центральной части потока осевая скорость имеет фактически постоянное значение.
Рис. 5 - Структура движения газодисперсного потока в плане
Численное моделирование осевых компонент скоростей потока показывает качественное согласие с расчетными данными. При этом расчетные данные хорошо согласуются с результатами в Ansys CFX при значениях удельного орошения до величины L/G =0,9 и существенно различаются при увеличении подачи жидкости. Это объясняется тем, что с ростом величины орошения реальное течение, исследуемое в численном эксперименте, теряет осесимметричность, а в расчетах в программном комплексе мы ставим условие осесимметричности потока.
Выводы
1. Разработан алгоритм моделирования процесса сепарации твердых частиц в динамическом газопромывателе. Выполнен анализ численной схемы и построена численная модель в пакете вычислительной гидродинамики Ansys CFX.
2. Выявлено наличие вихревых движений и получено распределение скоростей газодисперсного потока по сечению аппарата.
3. Результаты расчетов для различных режимов течения верифицированы с данными, полученными в программном комплексе Ansys CFX.
4. Разработанная модель помогает быстро и наглядно смоделировать движение запыленного газового потока с учетом внесенных в геометрию аппарата изменений. Таким образом, модель может применяться для оптимизации конструкции динамического газопромывателя.
Литература
1. D. C. Wilcox Multiscale Model for Turbulent Flows // In AIAA 24th Aerospace Meeting / American Institute of Aeronautics and Astronautics. (1986).
2. F. R. Menter, T. Esch Advanced Turbulence Modelling in CFX// CFX Update; Spring 2001. 20. 4-5. (2001).
3. Пат. РФ 2339435 (2008).
4. А. Ю. Вараксин Турбулентные течения газа с твердыми частицами. ФИЗМАТЛИТ, Москва, 66-170. (2003).
5. М. Г. Лагуткин, Д. А. Баранов. Оценка действия силы Кориолиса в аппаратах с закрученным потоком. ТОХТ.1. 9-13. (2004).
6. А. В. Старченко, А.М. Бубенчиков, Е.С. Бурлуцкий. Математическая модель неизотермического турбулентного течения газовзвеси в трубе. Теплофизика и аэромеханика, 1. 59-70. (1999).
7. M. Boivin, O. Simonin, K. D. Squires. Direct numerical simulation of turbulence modulation by particles in turbulence. Fluid Mech. Vol.375. 235-263. (1998).
8. M. Boivin, O. Simonin, K. D. Squires. On the prediction of gas-solid flows with two-way coupling using large eddy simulation. Phys. Fluids. Vol. 12. 8. 2080-2090. (2000).
9. C. T. Crowe, T. R. Troutt, J. N. Chung. Numerical models for two-phase turbulent flows. Annual Review Fluid Mech. Vol.28. 11-43. (1996).
10. O. Simonin, P. Fevrier, J. Lavieville. On the spatial distribution of heavy particle velocities in turbulent flow: From continuous field to particulate chaos. Turbulence. Vol. 3. 40-57. (2002).
© Р. Р. Усманова - канд. техн. наук, доц. Уфимского госуд. авиационного технич. ун-та; [email protected]; Г. Е. Заиков - д-р хим. наук, проф. каф. ТПМ КНИТУ, [email protected]; Р. Я. Дебердеев - д-р техн. наук, проф., зав. каф. ТППКМ КНИТУ.
