ТЕХНИЧЕСКИЕ НАУКИ
УДК 533.601
С.В. Карпов, А.А. Загоскин
К МЕТОДИКЕ АЭРОДИНАМИЧЕСКОГО РАСЧЕТА ЦИКЛОННЫХ УСТРОЙСТВ С ВНЕШНЕЙ РЕЦИРКУЛЯЦИЕЙ ГАЗОВ
Авторами [9] было предложено использовать новый принцип организации внешней рециркуляции газов в циклонных устройствах (топках, печах и др.), при котором транспортировка теплоносителя с боковой поверхности циклонной камеры, где температура его наиболее высокая, в приосевую область осуществляется за счет перепада давления между пристенной и приосевой зонами циклонного потока. Статья посвящена сравнению различных подходов к аналитическому описанию распределения тангенциальной скорости по радиусу в подобных устройствах.
Аэродинамика, циклонная камера, внешняя рециркуляция.
The authors suggest the new principle of the external gas recirculation in cyclone devices. According to this principle, the heat carrier can be moved from the walls of the device to the axis due to the radial pressure gradient. The paper considers analytical investigation of the aerodynamics of cyclone devices with external gas recirculation.
Aerodynamics, cyclone chamber, external gas recirculation.
Обозначения:
wф, wr, wz - тангенциальная, радиальная и осевая компоненты полной скорости;
г, ф, z - цилиндрические координаты;
Яю г я, Гфт, гвых - радиусы: циклонного устройства, ядра потока, положения максимума тангенциальной скорости, выходного отверстия;
^ф.вых, ^фя - тангенциальная скорость на радиусах
Гвых и Гя;
/вх, /вых - площади входа и выхода потока из циклонной камеры;
V, vt - молекулярная и турбулентная вязкости;
Qвх, Qpец - объемные расходы газа через циклонную камеру и рециркуляционный канал;
кс = Qpец / Qвх - коэффициент саморециркуляции.
Введение. Постановка задачи
Закрученный высокотурбулентный поток газов широко используется в технических устройствах энергетической, машиностроительной и металлургической промышленности благодаря увеличению доли конвективного теплообмена, выравниванию температуры поверхности заготовок, стабилизации процессов горения, увеличению полноты выгорания топлива и др. Применение циклонных нагревательных устройств позволяет увеличить скорость нагрева металла и уменьшить требуемый для этого температурный напор.
Простейшее циклонное устройство (рис. 1) состоит из цилиндрической камеры 1 и одного или двух тангенциально расположенных вводов 5, отвод газов обычно осуществляется через отверстие 6 в одном из его торцов. Благодаря взаимодействию входящей струи газа с боковой стенкой образуется закрученный поток.
Тангенциальная компонента скорости течения газов в пределах ядра закрученного потока является наибольшей, и поэтому ее следует считать основной. Из уравнения радиального равновесия потока можно сделать вывод [8], [10], что необходимым условием существования кругового течения в циклоне является равенство центробежной силы и радиального градиента давления, имеющего направление от оси к боковой стенке камеры. В циклонных устройствах статическое давление на стенке камеры может быть по величине соизмеримо со статическим давлением во входных каналах, в то время как на оси камеры наблюдается отрицательное избыточное статическое давление.
Если в конструкции циклонного устройства предусмотреть канал 3, соединяющий боковую стенку и приосевую зону, пе- ^
репускным вентилем 4 можно изменять перепад давления в циклоне, создавать внешний переток, газов из периферийной в приосевую зону ядра течения и, тем самым, регулировать его аэродинамику без изменения основных конструктивных параметров рабочего объема.
Экспериментальное исследование аэродинамики циклонных устройств с
Qb
ft
Qb
4
Qрец
Рис. 1. Схема экспериментальной установки
внешней саморециркуляциеи газов выполнено на металлической модели диаметром Ок = 201 мм, основные геометрические характеристики которой приведены в работе [9]. Для определения коэффициента саморециркуляции кс модель была дополнена камерой смешения 2, в мерном сечении которой определялся расход 0рец. В экспериментах кс изменялся в диапазоне 0.. .0,6.
Аналитическое решение
Основой аэродинамического расчета циклонных устройств является нахождение распределения тангенциальной скорости по радиусу, так как именно оно определяет важнейшие аэродинамические характеристики устройства [8]: величину крутки потока, максимальную тангенциальную скорость wtfm, радиус гфт распределения статического и полного давления по радиусу, коэффициенты аэродинамического сопротивления, аэродинамической эффективности и др.
Аналитическое решение задачи о распределении тангенциальной скорости по радиусу закрученного ядра потока в циклонном устройстве выполним при следующих допущениях:
1) ядро потока осесимметричное, стационарное, величины давления, тангенциальной и радиальной скоростей не меняются по высоте ядра. Таким образом: джф / дг = 0 , джф / дф = 0, джф / д2 = 0 ,
джг / дф = 0, др / дф = 0 . Справедливость этих допущений подтверждается многочисленными экспериментальными данными [8] для устройств с диаметрами выходного отверстия dвЬIX = ^вых/^к = 0,2...0,4 и /вх = 4/вх/ О > 0,04;
2) течение изотермическое, газ несжимаемый;
3) гипотеза Буссинеска адекватно описывает турбулентные свойства потока, что позволяет заменить Рейнольдсовы напряжения турбулентной вязкостью, значительно превосходящей молекулярную;
4) турбулентная вязкость не меняется по радиусу ядра.
Использование такого допущения, например, в [4], [5], значительно упрощает поставленную задачу и не приводит к большим отклонениям расчетного профиля скорости wtí> от экспериментального. В дальнейшем при расчетах будет использоваться ос-редненное по радиусу значение турбулентной вязкости.
Запишем уравнения движения для тангенциальной компоненты полной скорости и сплошности потока в цилиндрических координатах с учетом принятых допущений:
dw w w
^ yv ф yv ф r
wr-- + —--
dr r
( d 2 w„
1 dw„
d(rwr) dr
dr2
dwz
dz
r dr
0.
(1)
(2)
Для решения уравнения (1) необходимо предварительно задать распределение wr по радиусу. Наиболее целесообразны следующие пути решения:
1. Предположить, что wz во внутренней зоне ядра (0 < r < гвых) возрастает линейно по высоте от глухого торца к диафрагмированному, а во внешней зоне (гвых < r < гя) она равна нулю. Следуя данному предположению, из уравнения сплошности (2) можно найти распределение wr по радиусу и решить уравнение (1) раздельно для внешней и внутренней зон ядра потока.
2. Принять, что течение в ядре носит спиралевидный характер с постоянным углом закручивания а; в таком случае можно записать: wr = a w9, где a = ctg a = wr / w9 = const.
Решение системы уравнений (1, 2) первым способом основано на подходах, изложенных в работах [4], [5]. Перепишем уравнение (1) через wvr, приведем к безразмерному виду, используя выражения Wф = Wф/Wф.Bых, У = r/rBbK, и введем в рассмотрение безразмерный комплекс - радиальное турбулентное число Рейнольдса Rer = wrs/^ /vt. В результате получим следующие выражения для внутренней и внешней зон ядра потока:
d
(w9 W):
1
dy
Л
Re r • w
при r < r,
(w d2 (w) _ d(w^w)
v dw2 dW у
(3)
l{wg>w) . dw
Re Г
(w d 2 (w<Pv) d ( w) dw2
Л
dw
при Гвых< r < Гя.
(4)
Дифференциальные уравнения (3) и (4) решаются простой заменой переменных.
Опуская промежуточные преобразования, запишем:
- 1
(
wa =-
w
Re Г
• exp
- 1
w ф = —
w
Сз • V
i-D 2 ^
Re r • w 2
Rer+2
Л
+ C2
при r < Гвых; (5)
Re r + 2
- + C 4
при Гвых < Г < Г я. (6)
Постоянные интегрирования С1...С4 можно найти из граничных условий: при г = 0 wф = 0, при г = гвых = 1 и из условия непрерывности функции ^ (ж? у | / d у в области г = гвых получим:
C =
Rer
exp (Rer /2)_ 1
С2 = _
exp (Rer /2)_ 1
C3
C4 = 1 _
Rer • exp (Rer /2) exp (Rer /2)_1 , Rer • exp (Rer /2) (exp (Rer /2) _ 1) (Rer+ 2).
1
w
ф
+
v
2
1
r
При использовании фиксированного значения циркуляции тангенциальной скорости Г0 = wфяtyя на границе ядра потока постоянные интегрирования будут иметь следующий вид:
С, = Г о
( exp(Re, /2) -1 + exp(Re, / 2)(у^+2 -1) V1
Re,
Re, + 2
(7)
С,
С2 = -—Ч С3 = С, • exp(Re, /2): Re,
С = Г0 - Cl exp(Re, /2) • Vя
Re, + 2
(8)
При использовании констант С1...С4 в виде (7) -(8) в качестве скоростного масштаба можно выбрать среднюю скорость газа во входных каналах ует, что более удобно для практических расчетов.
На рис. 2а приведено радиальное распределение тангенциальной скорости, полученное по выражениям (5), (6) при различных значениях Яег и фиксированном значении Г0 на границе ядра. Из анализа профилей видно, что в зависимости от Яег положение максимума может находиться как во внутренней (при |Яег| > 2,513), так и во внешней (при |Яег| < 2,513) зонах ядра потока. Для связи величины гфт с Яег можно использовать следующую формулу, являющуюся аппроксимацией экстремума функций (5), (6):
фт
( 1.513 Л |Re ,1 -1
(9)
Решение уравнения (1) возможно при предположении о постоянстве угла закрутки [6]. Используя условие на границе ядра течения /w^ = w,/w9 = = const, уравнение (1) преобразуем к виду:
d ()
dr w,n ,я wp/ •, / ,я
d 2 ( wp, ) d ( wp, )
d,2
d,
. (10)
Выделив в правой части Яег = wrяrя/vt и, используя переменные = ^ф/^фя, Е = г/гя, представим уравнение (10) в безразмерном виде:
d
(ф%)
1
( d2 (ф%) d (ф%)
d% Re,Wф^
Л Л
d%2
d%
. (11)
Решение дифференциального уравнения (11) выглядит следующим образом:
— 1 1
w^ =-(Re,2 + 4)2 х
Re, I 1 ' )
1 ( 1 1 : th -1 (Re,2 + 4)2 С2 -(2C1Re,2 + 4)2 ln(%)
//
- 2, (12)
где С и С2 - постоянные интегрирования.
Выполним замену (2С^ег2 + 4)12 = у и упростим уравнение (12), используя (13):
(1 л ^С2 )_ 1
л [ 2 (( _ >""« »)-'-^Щ^1 =
вуС %-Y +1 - 2 = 1 - 2
еУС2 %-Y
г+1 (С2/ % ) +1'
С2 ) у + 2 + Re, wpE
Е J у-2 - Re,
(13)
(14)
Из граничных условий задачи найдем постоянные интегрирования С1 и С2: 1) при Е = 0 и = 0 у = 2; 2) при Е = 1 и = 1 С22 = - (4 + Яег) / Яег. После подстановки С1 и С2 в уравнение (14) получим:
4%
Re, (1 -Е2)+ 4'
(15)
На рис. 2б показано радиальное распределение тангенциальной скорости, рассчитанное по формуле (15) при различных значениях Яег.
Найдем положение и величину максимума тангенциальной скорости, определив экстремум функции (15), т.е. (й ^ф)/ йЕ = 0:
4 w 2 J-4/Re, -1 —— 1; = ^-:-. (16)
,я V Re, wp„„ Re, + 4
Используя выражения (16), перейдем к безразмерным координатам w = w9/ wфm, n = ,/,фт:
- 2 J-4/Re, -1 = 4^-4/Re, -1 =
Re,
4 Re, (1 - n2 (-4/Re, -1)) + 4
4nV-4/Re, -1 (n2 +1)(Re,+ 4).
(17)
Формула (17) после сокращений преобразуется к виду известной аппроксимации V Вулиса-Устименко [3] (такая же зависимость получена экспериментально при изучении вращательного движения в вихревых воронках [1]) и становится независимой от Rer:
- 2n w —
n2 +1
(18)
Выражения для безразмерного радиуса пя и числа Rer (16) можно представить в виде:
Пя ^ —=
Re,
4 + Re,
;Re,
-4n2 n2 +1
(19)
0.38
,
вых
Зависимость (18) представим в обобщенном виде
[11]
2П
П2 +1
(20)
где п - показатель степени, определяемый условиями генерации закрутки; при использовании условия достижения максимума циркуляции Г = ю-ц на границе П2. +1 2
ядра
п2 -1
Яе.. + 2
Как видно из рис. 2, в обоих случаях с уменьшением Яег максимум юф смещается к периферии, а общий уровень снижается. В предельном случае при Яег=0 отсутствует квазипотенциальное течение, и весь объем газа вращается как квазитвердое тело. Существование данного режима подтверждается результатами численного моделирования циклонного течения [7] при кс = 1, а также экспериментальными данными [2]. Из рисунка видно, что при одинаковом
значении Яег из-за различия в допущениях, сделанных при анализе системы уравнений Навье-Стокса, метод, основанный на постоянном угле закрутке, дает меньший уровень вращательных скоростей, что необходимо учитывать при дальнейшем построении методики аэродинамического расчета рассматриваемых устройств.
Сравнение с экспериментальными данными
Использование радиального турбулентного числа Рейнольдса удобно для анализа аэродинамики циклонных рециркуляционных устройств с отбором газа с их боковой поверхности, так как: 1) исходя из уравнений (9), (19) и эмпирических зависимостей для гя и гфт, приведенных в работе [8], оно должно определяться соотношением площадей входа /вх и выхода /вых потока; 2) величина радиального расхода 'Гягя изменяется обратно пропорционально кс; 3) изменение величины V/; также можно представить в виде функции кс; 4) при кс = 1 значение Яег = 0.
а) . \
X
\
-----
0.2
0,4
0.6
0.8
г1г, "0
б) ч
/ ч
/
0.2
0.4
0.6
0.8
г /г,
Рис. 2. Распределение безразмерной тангенциальной скорости по радиусу циклонной камеры при различных |Яег| и методах расчета: а) по (5, 6) (^вых = 0,2); б) по (15)
1^1 3.5
2,5
1.5 I
3.5 3
2.5
0.3
а) а)
^------ о
) 1 -о .2-0
1.5
0-4 0.5 0,6 0.7 /„¡/т,(1~К)
б)
\
1 -гг 2 -О ? ~ *
0.5
0.2
в) /о
р ^ \ 4 \в
\
□ / 3
1? )
0-3 0.4 0.5 0.6 /„//„(1-0 0 0. I 0,2 0.3 0.4 0.5 0.6 г/К,
Рис. 3. Изменение |Яег| в зависимости от кс: а) расчет по формулам (5), (6): 1 - ^вых = 0,2; 2 - ^вых = 0,25; б) расчет по формуле (19): 1 - 5вых = 0,2; 2 - ^вых = 0,3; 3 - ^вых = 0,4; в) сопоставление расчетных и экспериментальных профилей = '</увх; эксперимент: 1 - ^вых = 0,3, кс = 0,47; 2 - 5вых = 0,2, кс = 0,26; расчет: 3 - по (20);
4 - по (5), (6)
Расчет числа Яег для метода, основанного на постоянном угле закрутки, выполнялся по формуле (19) с привлечением экспериментальных значений гя и гфт, а для метода, основанного на гипотезе о линейном возрастании осевой скорости по высоте камеры, - исходя из наименьшего расхождения расчетного и экспериментального профилей wч,. Значения Яег для каждого случая показаны на рис. 3а, 3б и обобщены выражениями вида
(
|Rer| = A А-(1 -.
/в
(21)
В формуле (21) А, а, в - опытные константы: при использовании метода с постоянным углом закрутки и/вх//вых < 0,52 А=4,5, а=1,6, в=0,31; при> 0,52 А=3,7, а=1,6, в=0,04; при использовании второго метода А=3,7, а=0,7, в=1. Расчетные кривые показаны на рис. 3а, б линиями, экспериментальные данные -точками.
В таблице приведены среднеквадратичные отклонения расчетных по (формулам 5, 6, 20) и экспериментальных профилей Как видно из представленных данных, расчет по (5), (6) дает меньшее отклонение при йвых < 0,3, а расчет по (20) - при йвых > 0,3. Средневзвешенные значения отклонений составили по первой методике 0,053, по второй - 0,057.
Таблица
Среднеквадратичные отклонения расчетных и экспериментальных значений wv
Среднеквадратичное отклонение
dвых kc Расчет по формулам (5), (6) Расчет по формуле (20)
0,2 0,26 0,021 0,034
0,2 0,30 0,020 0,041
0,2 0,57 0,028 0,072
0,3 0,22 0,097 0,083
0,3 0,35 0,075 0,063
0,3 0,51 0,040 0,058
0,4 0,20 0,070 0,033
0,4 0,35 0,068 0,062
0,4 0,41 0,059 0,064
На рис. 3в в качестве примера показано сопоставление расчетных и опытных радиальных распределений в объуме циклонной рециркуляционной камеры при йвых=0,2, Лс=0,26 и йвых=0,3, &с=0,47. Как видно, результаты эксперимента в целом хорошо согласуются с расчетными кривыми.
Выводы:
1. Выполнено сравнение различных подходов к аналитическому расчету профиля тангенциальной скорости. Получены теоретические зависимости вида w9=/(r, Rer/) для двух методов расчета и определены целесообразные для их применения диапазоны варьирования основных геометрических характеристик циклонной камеры.
2. Предложена методика расчета радиального распределения тангенциальной скорости в циклонном устройстве с боковым отбором газа на рециркуляцию, основанная на использовании радиального турбулентного числа Рейнольдса. Представлены уравнения для расчета Rer при различных значениях соотношения /,х/,ых и коэффициента kc.
Литература
1. Альтшуль, А.Г. Обобщенная формула распределения окружных скоростей в вихревых воронках / А.Г. Альтшуль, М.Ш. Марголин // Инж.-физ. журн. - 1970. - Т. 18.-№ 4. - С. 731 - 733.
2. Брук, Ю.Г. Исследование движения газов и конвективного теплообмена в вертикальных рециркуляционных печах / Ю.Г. Брук, А.У. Пуговкин // Оптимизация металлургических процессов. - М., 1971. - Вып.5. - С. 351 - 355.
3. Вулис, Л.А. Об аэродинамике циклонной топочной камеры / Л.А. Вулис, Б.П. Устименко // Теплоэнергетика. -1954. - № 9. - С. 3 - 10.
4. Гольдштик, М.А. Вихревые потоки / М.А. Гольд-штик. - Новосибирск, 1981.
5. Гольдштик, М.А. К теории эффекта Ранка (закрученный поток газа в вихревой камере) / М. А. Гольдштик // Известия АН СССР. Отд. техн. наук. Механика и машиностроение. - 1963. -№ 1. - С. 132 - 137.
6. Гришанин, К.В. Динамика донного течения в стационарном водовороте // Сборник трудов Ленинградского института водного транспорта / К.В. Гришанин. - 1964. -Вып. 61. - С. 36 - 43.
7. Загоскин, А.А. О физическом и математическом моделировании аэродинамики циклонных устройств с внешней саморециркуляцией газов / А.А. Загоскин, С.В. Карпов, Э.Н. Сабуров // Проблемы газодинамики и тепломассообмена в энергетических установках: Труды XIX Школы-семинара молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева (20-24 мая 2013 г., г. Орехово-Зуево). - М., 2013. - С. 102 - 105.
8. Карпов, С.В. Высокоэффективные циклонные устройства для очистки и теплового использования газовых выбросов / С.В. Карпов, Э.Н. Сабуров. - Архангельск, 2002.
9. Карпов, С.В. К вопросу о физическом моделировании внешней саморециркуляции газов в циклонных устройствах / С.В. Карпов, Э.Н. Сабуров, А.А. Загоскин // Сб. трудов VII Международной науч.-техн. конференции «Автоматизация и энергосбережение». - Вологда, 13-15 марта 2012 г. - С. 174 - 177.
10. Сабуров, Э.Н. Теория и практика циклонных сепараторов, топок и печей / под ред. Э.Н. Сабурова / Э.Н. Сабуров, С.В. Карпов. - Архангельск, 2000.
11. Штым, А.Н. К аэродинамике закрученного потока в циклонно-вихревых камерах / А.Н. Штым, П.М. Михайлов // Известия вузов. Энергетика. - 1965. - № 11. - С. 50 - 53.