CC BY
91
УДК 621.438
МЕТОД РАСЧЕТА ИНТЕНСИВНОСТИ ТУРБУЛЕНТНОСТИ ПОТОКА ГАЗА НА ВЫХОДЕ ИЗ КАМЕРЫ СГОРАНИЯ ГАЗОТУРБИННОГО ДВИГАТЕЛЯ
© 2007 А. В. Григорьев, В. А. Митрофанов, О. А. Рудаков, А. В. Соловьева
ОАО «Климов», г. Санкт-Петербург
Обоснована формула для оценки величины интенсивности турбулентности потока газа на выходе из камеры сгорания газотурбинного двигателя и показана связь этой величины со значениями ее конструктивных параметров.
В настоящее время при решении задачи гидроаэродинамики вычислительными методами для газотурбинного двигателя возникает проблема оценки параметров турбулентности потока на выходе из такого элемента конструкции, как камера сгорания.
В работах [1, 2, 3] авторами предложен интегральный подход к оценке характеристик турбулентности потока газа, связывающий параметры турбулентности с размерами канала переменной или постоянной геометрии, в котором происходит течение потока сплошной среды.
Запишем уравнение Навье-Стокса для установившегося турбулентного течения, представив истинную скорость потока сплошной среды как сумму средней скорости потока (с осреднением по Рейнольдсу) и ее пульсационной составляющей в форме
Э W Э W Э W
P и-----+ P v---+ P w----+
Э x Э у Э z
,Э W' ,Э W' ,Э W'
+ P и---------+ P v--------+ P w--------=
Э у
(э 2w э 2w э 2w л
Э x2 Э y 2 Э z2
/эр + др +эрл
Э x Э у Э z
правую части на величину ё¥=ёхёуё1, а также вводя понятие элементарного массового расхода ёО жидкости или газа через элемент поверхности, нормальный к линии тока, являющейся диагональю параллелепипеда со сторонами йхйуйг, ограничивающими этот элемент:
ёЕ = dydzcos (Ь,х) = dxdzcos (Ь,у) = ёхЛуео8 (Ь^);
= W; dG = p W dF
cos (L,x) cos (L,y)
cos L,z
(L,z )
где dF - площадь элемента поверхности, нормального к линии тока; L обозначает направление вдоль линии тока.
Преобразованное уравнение с использованием формулы Грина, а также принимая
jlj = const в рассматриваемом объеме, интегрируем в следующих пределах:
G W2
Г dG Г
Gl W1
dG d W = m,- l^dS -
где W - вектор средней скорости потока сплошной среды (жидкости, газа); Р, р - давление и плотность потока сплошной среды; п,\, Н - проекции вектора средней скорости потока на декартовы оси координат, цг - коэффициент динамической вязкости жидкости или газа; знак «'» обозначает пульсационные составляющие скорости потока (знак, обозначающий осреднение по времени, опускаем).
Произведем преобразование написанного выше уравнения, умножив его левую и
d W С ГЭU
X
А
Э n
„SP ггг/ар +ар +эр V -
X JJ Дэ x Э у Э z f-
W'
- Fm Г P W'dW';
-W'
где Оь - массовый расход жидкости (газа)
через трубку тока; G - массовый расход жидкости (газа) через рассматриваемый объем
Ур; 8р - площадь поверхности, охватываю-
x
z
щей рассматриваемый объем; - площадь
элемента поверхности, охватывающей рассматриваемый объем; Г - площадь проходного сечения канала рассматриваемого объема; 5 - изменение величины параметра от выхода из рассматриваемого объема до вхо-
— Ш - Р да в этот объем; и = —— ;Р = —- - коэффи-
о Ш о Р
циенты, учитывающие неравномерность распределения параметров в пространстве рассматриваемого объема (параметры в числителе представляют собой местные значения);
п =
-1 +
cos (п,х) cos (п,у) cos (п,2 )
_ х _ у г
х = —; у = —; г = —; п - нормаль к поверхнос-XXX
ти, охватывающей рассматриваемый объем; С - характерный размер рассматриваемого объема; 1, _], к - единичные орты; т - средняя величина параметра в диапазоне его изменения от выхода из рассматриваемого объема ко входу в этот объем; индексы «1» и «2» обозначают соответственно вход в рассматриваемый объем и выход из него.
Изменение величины давления на протяжении области турбулентного движения определим в соответствии с [2]:
5 * Ї + в І-'
Яв
*
здесь К, В - коэффициенты пропорциональности; I - относительная длина рассматри-
р
ваемого объема.
Произведя интегрирование и преобразования, получим
—5 *
Ур
1 - САи^—
Ур
X
+ Р (* У Fm =
Кр т*„
/ 5 * 4 2
їїт
I
+ в-^
Яв
Ар;
где Аи =-
т* ш
Л_____.. А--Ур-
ЭР ЭР ЭР
—И +—3 +—3
Э х Э у Э 2
СхсуССг
функции, учитывающие влияние неравномерности распределения параметров в пространстве рассматриваемого объема; С - коэффи-
циент пропорциональности; комплекс
X
—
ур
представляет собой соотношение скоростей молекулярного и конвективного переноса ко -личества движения в потоке [(кг/с м3)/
(кг/с м3)], с другой стороны
X
—
ур
, _Р_
Яв .
Правая часть написанного выражения характеризует работу сил давления, обусловленных переменной геометрией канала и сил трения, затрачиваемую на изменение количества движения потока в пространстве рассматриваемого объема и «порождение» турбулентности; слагаемые в левой части характеризуют перенос потоком импульса с учетом вязкости среды и кинетическую энергию турбулентных пульсаций скорости потока в пространстве рассматриваемого объема.
Окончательно после преобразований получаем
є2 = КАр
5 *
\2 -------------
+ в
Я.е
Рт X
— Ь - А*Сі ^
*т * 1Ке
(1)
Левая часть полученного соотношения
(шО2
есть критерий Кармана е
*
фициент пропорциональности.
В общем виде для потока с неравномерным распределением параметров в пространстве рассматриваемого объема Ар, А^ являются функцией гидродинамических, гео-
5 * Т
метрических критериев
і
Ур
И
*т Яв’ ^т X
начальной неравномерности распределения
і
2
2
2
параметров на входе в рассматриваемым объем.
Соотношение (1) подтверждено опубликованными экспериментальными данными по измерению параметров турбулентности в конфузорном канале, смешивающихся спут-ных турбулентных струях и в пограничном слое, набегающего на пластину потока [2].
В работе [4] приведены данные по измерению интенсивности турбулентных пульсаций потока воздуха при втекании в него поперечных струй на модели смесителя камеры сгорания с варьированием следующих
параметров: —L - отношения массовых рас-^2
ходов воздуха в поперечных струях и в ос-
новном потоке;
SmFa
дей поперечного сечения канала для основного потока и суммарной эффективной площади отверстий для подвода воздуха в основной поток (т - коэффициент расхода); t -шага между отверстиями подвода воздуха; l -длины от отверстий до сечения измерения интенсивности турбулентных пульсаций. Измерения производились при условии изо-термичности взаимодействующих потоков и постоянства геометрии канала для основного потока с H = const, где H - высота исследуемого канала.
Пределы изменения варьируемых параметров в [З]:
GF —^ » 0,2...0,45; —Рж—
G2 ЇИРОЖ
З...1З;
— » 0,44... 1,02; — » 0,5... 1,8. HH
Экспериментальные данные в [3] описываются выражением (1), которое после их обработки приобретает вид
г2 =(1 +1)0,5
0 ,029
Ґд WV
W
v m
+ 0,0748
S W
Wm
+ 0,002; (2)
здесь
—_ Рж -1
sw —2 їиРож
W
Vp t / l
p = const; l = —.
H
На рис. 1 представлено сопоставление расчетных (обозначение на рисунке - «расч.») и экспериментальных данных (обозначение на рисунке - «эксп.»).
Перейдем к интегральной оценке интенсивности турбулентных пульсаций в камере сгорание, используя (2). Для всего объема жаровой трубы камеры сгорания, где происходят процессы турбулентного смешения,
с учетом того, что —- = 1, выражение (2) приобретает вид
- отношения площа- г2 = (1 + lr )-°':
0,029
-1
- +1
+ 0,0748
-1
-+1
+ 0,002;
(З)
где /г - относительная длина от последнего
ряда отверстий жаровой трубы до сечения выхода из камеры сгорания (относительная длина газосборника жаровой трубы камеры сгорания).
Рис. 2 иллюстрирует влияние величи-
РЖ
ны параметра у ^ на интенсивность тур-
И^ОЖ
булентных пульсаций в камере сгорания. Расчет произведен по (3) при /г = 1,1. По оси
ординат графика на рис. 2 отложены значения интенсивности турбулентных пульсаций в процентах.
График на рис. 3 иллюстрирует влияние на «затухание» интенсивности турбулентных пульсаций длины газосборника жаро-
вой трубы при
Ж
9.
Расчетные соотношения (2) и (3) могут быть использованы для оценки интенсивности турбулентных пульсаций на входе в сопловой аппарат турбины при решении дифференциальных уравнений газовой динамики численными методами применительно к газотурбинному двигателю.
F
F
Ж
Ж
F
F
Ж
Ж
є2 (зксп.)
0,05
*
ж /
і *
У □
0,04
0,03 0,02 0,01 0
О 0,01 0,02 0,03 0,04 0,05
е2(расч.)
Рис.1. Сопоставление расчетных и экспериментальных данных по измерению интенсивности турбулентных пульсаций
I Р О -
А — = 1,02; = 13,2;= уаг; / = 1,1;
• - Н 02
р.
О =,
А — = 0,44; Ж = 3,25; ^ = уаг; / = 1,1; А - н т,р о
№
02
ґ р о -
п _ = 0,44;——— = 6,5;^ = уаг;/ = 1,1; и - Н 02
ґ
01
А — = 0,44;—= 4;—1- = 0,12;/ = 1,1; * - Н ХтРож -
02
ґ р о -
Ж - Т7 = 0,44; „ Ж = 4;^- = 0,22;/ = уаг_
Н
О2
- экспериментальные точки,
Ож и2
соответствующие различным геометрическим и гидродинамическим параметрам
испытываемых моделей [4]
интенсивность турулентных пульсаций %
р
Ж/
'ЦРс
ОЖ
Рис. 2. Влияние на интенсивность турбулентных пульсаций относительной суммарной площади отверстий жаровой трубы камеры сгорания
интенсивность турулентных пульсаций %
Гг
Рис. 3. Влияние на интенсивность турбулентных пульсаций длины газосборника жаровой трубы камеры сгорания
Список литературы
1. Митрофанов В. А., Рудаков О. А. Определение параметров турбулентности потока газа в камере сгорания на основе интегрирования и преобразований уравнений Рейнольдса и Навье-Стокса// Вестник СГАУ
- 2003. - №2. - С.92-96.
2. Митрофанов В. А., Рудаков О. А Критериальная модель турбулентных явлений в потоке сплошной среды при ее движении по каналу переменной геометрии // Пром. теплотехника.- 2004. - №3.
3. Григорьев А. В., Митрофанов В. А., Рудаков О. А., Сигалов Ю. В. Теоретическое обоснование критического числа Рейнольдса как функции гидродинамических и геометрических критериев потока // Уткинские чтения. Материалы конференции. Т.2. - СПб.: БГТУ, 2005.
4. Безменов В. Я., Онищик И. И. Экспериментальное исследование процесса смешения газовых потоков в моделях смесителей кольцевых камер сгорания ГТД / Технический отчет. - ЦИАМ, 1965.
COMPUTATIONAL METHOD OF THE GAS FLOW TURBULENCE INTENSITY AT THE GAS TURBINE ENGINE COMBUSTION CHAMBER OUTLET
© 2007 A. V. Grigoriev, V. A. Mitrofanov, O. A. Rudakov, A. V. Soloveeva
JSC Klimov
Formula for evaluation of the gas flow turbulence intensity value at the gas turbine combustion chamber outlet was substantiated and relation of this value to its design parameters was shown.
