УДК 532.6
А. Н. Аксенов, А. Б. Шабаров
ДЕМПФИРОВАНИЕ КОЭФФИЦИЕНТА ТУРБУЛЕНТНОЙ ВЯЗКОСТИ ПРИ РАСЧЕТЕ ТЕЧЕНИЙ В ОСЕВЫХ ТРАНСЗВУКОВЫХ КОМПРЕССОРАХ
Приведены результаты моделирования течения в трансзвуковом осевом компрессоре НАСА Rotor 37 с помощью ks-модели турбулентности и сравнение расчетных распределений степени повышения полного давления и полной температуры с экспериментальными данными. Для повышения точности моделирования предлагается демпфирование турбулентной вязкости в зависимости от отношения генерации турбулентной энергии к скорости ее диссипации.
E-mail: muon@inbox.ru; ashabarov@utmn.ru; tu1487@mail.ru
Ключевые слова: трансзвуковые компрессоры, численное моделирование, неравновесная турбулентность, демпфирование турбулентной вязкости.
Результаты численного моделирования течений в проточной части высоконагруженных осевых компрессоров значительно отличаются от экспериментальных данных (в пределах 4 % по интегральным характеристикам). Во многом указанное отклонение обусловлено расчетом турбулентных характеристик, поэтому совершенствование моделей течения — актуальная задача.
Трансзвуковой осевой компрессор НАСА Rotor 37 (рис. 1) характеризуется степенью повышения полного давления п^ = 2,056 и адиабатическим КПД nad = 87,6 % (при массовом расходе воздуха 20,74 кг/с). Лопатки компрессора изготовлены из мартенситно-стареющей стали maraging 200, обладающей высоким сопротивлением малым пластическим деформациям, коррозионной стойкостью и размерной стабильностью. Выбор данного компрессора в качестве объекта исследования обусловлен доступностью экспериментальных данных [1]-[3] и
____i___L
Рис. 1. Проточная часть компрессора Rotor 37
теоретическими работами [4, 5], позволяющими провести проверку методики расчета и оценить адекватность моделирования.
Математическая постановка задачи. Во вращающейся с постоянной угловой скоростью Ш системе координат уравнения сохранения массы, момента импульса и энергии имеют вид [4]
dt + v- (pW ) = 0;
(1)
dpW dt
+ v • (pW • W) + p(2ÜÜ x W + üj x ÜÜ x r) = -Vp + Vfr; (2)
dpEr dt
+ V • (pWHr) = V • ((X + cpptJPrt)VT + fr • W). (3)
В уравнениях (1)-(3) р — плотность; V = IV + Ш х г — вектор абсолютной скорости; IV — вектор скорости в относительной системе координат; г — радиус-вектор; р — статическое давление; Л — коэффициент теплопроводности; Т — температура; тг = (р + рг) • ('УШ + (V Ш)т) — тензор вязких и турбулентных напряжений; р — ламинарная вязкость; рг — турбулентная вязкость; ср — удельная изобарная теплоемкость; Рг — турбулентное число Прандтля. Значения коэффициентов теплопроводности Л и ламинарной вязкости р считаются постоянными во всей расчетной области. В уравнении сохранения энергии (3) обозначено:
Р
Hr = Er +—;
P
2„2
p W2 — ü2r
Ег = срТ---+
р 2
Для замыкания системы (1)-(3) использовалось уравнение состояния для совершенного газа
(cp cv
P
(4)
где сь — удельная изохорная теплоемкость. Расчет турбулентной вязкости рг проводили с помощью к—е-модели турбулентности [6]:
C k2 Pt = p—;
dpk
+ V • (pWk) = V
P + p ) Vk
+ Pk — pe;
(5)
(6)
dpe ~öt
+ V • (pWe) = V •
P + P ) Ve
+ k (CiPk — C£2pe); (7)
Pk — ptS ; S — sj2SijSij; Sij —
1 iöWi dWj
2 V дх^
+
dxt
С = 0,09; Ое1 = 1,46; Се2 = 1,83; ак = 1; = 1,3,
где к — кинетическая энергия турбулентных пульсаций; г — скорость диссипации энергии турбулентности; Рк — генерация турбулентной энергии сдвигом скорости; — компоненты осредненного тензора скоростей деформации.
При определении граничных условий на входе принимались полное давление (1,013 • 105 Па), полная температура (288,15 К), степень турбулентности (Ти = ^2к/3/(# + Л х г)2 = 0,05) при входе и массовый расход (20,51 кг/с) при выходе из рабочего колеса. Условие прилипания Ш = 0 использовалось на поверхностях рабочих лопаток и диске компрессора, на втулке и корпусе V = 0. На всех непроницаемых границах ставилось условие отсутствия теплового потока через стенку. Для сокращения числа расчетных узлов вблизи твердых стенок применялись масштабируемые пристеночные функции [6, 7].
Решение системы (1)-(7) проводилось методом установления по времени (Д£ ~ 2 • 10-5 с) с консервативной аппроксимацией потоков по схемам второго порядка точности (реализация АшуБ СБХ11). За
Рис. 2. Фрагмент расчетной сетки в межлопаточном канале (250 тыс. расчeтных узлов)
критерий сходимости принята среднеквадратичная невязка 5-10-6. Интегрирование проводилось на структурированной многоблочной сетке (рис. 2) с максимальным общим числом расчетных узлов 1,5 • 106.
Результаты численного моделирования. Расчетные распределения п* и Тв*ых/Тв*х при выходе из рабочего колеса (сечение А-А, см. рис. 1) приведены на рис. 3, а и рис. 3, б соответственно. Моделирование проводилось на сетках с различным разрешением межлопаточного канала (0,25, 0,5 и 1,2 млн расчетных узлов). Для сравнения приведены результаты, полученные с помощью SST (Shear Stress Transport) модели Ментера [8].
Результаты вычислений по модели k—e характеризуются завышенными значениями п* и Тв*Ь1х/Тв*х; расчеты с использованием уравнений SST лучше согласованы с экспериментальными данными по Тв*ых/Тв*х и имеют тенденцию к заниженным оценкам п*.
Демпфирование турбулентной вязкости в к—^-модели турбулентности. Для повышения точности расчета сдвиговых течений Роди [9] предложил считать коэффициент турбулентной вязкости C^ функцией отношения генерации турбулентной энергии G = Pk/р (за счет деформации профиля средней скорости) к скорости диссипации e этой энергии в теплоту:
C = C0 • g((G/e)*),
где C0 = 0,09; (G/e)* — характерное интегральное значение неравновесности генерации и скорости диссипации энергии турбулентности в расчетной области.
В разработанной Лаундером и другими [10] модели ke2 использование зависимости g((G/e)*), представленной на рис.4, улучшает результаты расчета осесимметричных течений.
Численное решение стандартной системы (1)-(7) показало, что течение газа в трансзвуковых осевых компрессорах характеризуется высокой степенью неравновесности: G/e ~ 2... 3 в отрывных областях и до 103 вблизи скачка уплотнения (рис. 5). Таким образом, необходимость демпфирования турбулентной вязкости (рис. 6) согласуется с результатами, приведенными на рис. 3 и 4.
Калибровка коэффициента турбулентной вязкости на течение в компрессоре Rotor37 позволила получить C = 0,065 при (G/e)* = = 1,35, что находится в полном соответствии с результатами Лаунде-ра и Роди [10] (см. рис.4). При расчете газодинамических характеристик компрессора применялось демпфирование, основанное на слабом локально-неравновесном приближении [11] (рис.7,а, б):
G/e
C = min<¡ C0, A • exp
(G/s)*
(8)
Рис. 3. Распределения степеней повышения полной температуры (а) и полного давления (б), рассчитанные при различном разрешении межлопаточного канала по стандартным моделям турбулентности
1,5
1,0
0,5
\ Настоя \ L щая работа -aunder ке2 - = 0,06Е •
\
_____________
0,5
1,0
1,5
G/e
Рис. 4. Локальное и интегральное демпфирование коэффициента турбулентной вязкости
Рис. 5. Неравновесность генерации и скорости диссипации энергии турбулентности в межлопаточном канале компрессора (в среднем по высоте сечении)
Рис. 6. Поле турбулентной вязкости в межлопаточном канале компрессора
н
Ж,
к
и
Эксперимент к-е (1200 К)
к-е(Сц = 0,065,1200 К) к—е(Сц(д), 1200 К)
Рис. 7. Распределения степеней повышения полной температуры (а) и полного давления (б), рассчитанные при локальном и интегральном демпфировании коэффициента турбулентной вязкости
где A = C0/exp[-((G/e)T1; ((G/^T = -nC0/C; C = 2 x
x
1
CR
2
C2 — TT" (C2 — 1 + CR) CR
3
; C2 = 0,6; CR = 2,486; n = 0,45.
Заключение. Установлена существенная неравновесность генерации и скорости диссипации энергии турбулентности в отрывных зонах со стороны спинки рабочей лопатки и на скачках уплотнения. Интегральная поправка Лаундера-Роди повышает точность расчета газодинамических характеристик компрессора. Предложенное локально-неравновесное демпфирование осложняет сходимость численного интегрирования, но улучшает согласованность с экспериментальными данными и применимо для случаев с большими значениями G/e.
СПИСОК ЛИТЕРАТУРЫ
1. R e i d L. and M o o r e R. D. Design and overall performance of four highly loaded, high-speed inlet stages for an advanced high-pressure-ratio core compressor // NASA Technical Paper 1337, 1978.
2. M o o r e R. D. and L o n n i e R. Performance of single-stage axial-flow transonic compressor with rotor and stator aspect ratios of 1.19 and 1.26, respectively, and with design pressure ratio of 2.05 // NASA Technical Paper 1659, 1980.
3. S u d e r K. L. Blockage development in a transonic, axial compressor rotor // NASA Technical Memorandum 113115, 1997.
4. Dunham J. (editor). CFD validation for propulsion system components // AGARD Report AR-355, 1998.
5. C h i m a R. V. and L i o u M. Comparison of the AUSM+ and H-CUSP schemes for turbomachinery applications // NASA Technical Memorandum 212457, 2003.
6. L a u n d e r B. E. and S p a l d i n g D. . The numerical computation of turbulent flows // Computer Methods in Applied Mechanics and Engineering. - 1974. - Vol. 3. - P. 269-289.
7. Grotjans H. and M e n t e r F. R. Wall functions for general application CFD codes / In K.D. Papailiou et al., editor // ECCOMAS 98 Proceedings of the Fourth European Computational Fluid Dynamics Conference, 1998. - P. 1112-1117.
8. Menter F. R. Two-equation eddy-viscosity turbulence models for engineering applications // AIAA Journal. - 1994. - Vol. 32. - P. 269-289.
9. R o d i W. The prediction of free turbulence boundary layers by use of a two-equation model of turbulence / Ph. D. dissertation, University of London, 1972.
10. L a u n d e r B. E. et al. Prediction of free shear flows: A comparison of the performance of six turbulence models // NASA Report SP-321. - 1973. - P. 361-422.
11. P o p e S. B. Turbulent flows. - Cambridge University Press, Cambridge, 2000. -750 p.
Статья поступила в редакцию 23.04.2009
Андрей Николаевич Аксенов родился в 1985 г., окончил Тюменский государственный университет в 2006 г. Сотрудник Тюменского государственного университета, ПИИ ОАО "Газ-ТурбоСервис". Автор восьми научных работ в области аэродинамики, процессов горения, турбомашиностроения.
A.N. Aksenov (b. 1985) graduated from the Tyumen' State University in 2006. Worker of the Tyumen' State University and JSC "GasTurboServis". Author of 8 publications in the field of aerodynamics, combustion processes, turbo-machine building.
Александр Борисович Шабаров родился в 1944 г., окончил МВТУ им.Н.Э.Баумана в 1967г. и МГУ им.М.В.Ломоносова в 1969 г. Д-р техн. наук, профессор Тюменского государственного университета. Автор свыше 200 научных работ в области механики многофазных систем, газотурбостроения, теплофизики, нефтегазовых технологий.
A.B. Shabarov (b. 1944) graduated from the Bauman Moscow Higher Technical School in 1967 and Lomonosov Moscow State University in 1969. D. Sc. (Eng.), professor of the Tyumen' State University. Author of more than 200 publications in the field of mechanics of multiphase systems, gas-turbo-building, thermal physics, oil and gas technologies.