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

УДК 621.45.037

РАСЧЁТ ХАРАКТЕРИСТКИ РАБОЧЕГО КОЛЕСА В ДВУХМЕРНОЙ ПОСТАНОВКЕ С ИСПОЛЬЗОВАНИЕМ ИМИТАЦИОННОГО МОДЕЛИРОВАНИЯ

© 2011 А. Б. Михайлова, Д. А. Ахмедзянов, Ю. М. Ахметов, А. Е. Михайлов

Уфимский государственный технический университет

Разработана и реализована в системе имитационного моделирования COMPRESSOR_2D методика решения прямой задачи в двухмерной постановке для осевого компрессора. Проведено исследование адекватности предложенного комплексного критерия возникновения срывных зон.

Осевой компрессор, имитационное моделирование, характеристика компрессора, граница устойчивой работы, углы атаки и отставания.

Рабочий процесс компрессора - сложное пространственное явление, протекающее в вязком нестационарном потоке.

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

Таким образом, рабочий процесс в межлопаточном канале определяется углом поворота и потерями полного давления на различном уровне по высоте лопаточного венца. Потери полного давления выступают функцией параметров геометрии лопаточного венца, числа Маха на входе в межлопа-точный канал, коэффициента диффузорно-сти, радиального положения в лопаточном венце и осевого положения лопаточного венца. Угол поворота потока в межлопаточ-ном канале является функцией параметров геометрии лопаточного венца, числа Маха на входе в межлопаточный канал, радиального положения в лопаточном венце [1].

Наиболее важным этапом проектирования компрессора в газодинамическом аспекте является определение допустимой аэродинамической нагрузки профилей. Это было отмечено на ранней стадии развития методов проектирования компрессоров, в

работах A.R. Howell, которые до настоящего времени используются при проектировании компрессоров в промышленности.

Наиболее важным критерием оценки аэродинамической нагрузки решеток является коэффициент или степень диффузорности (diffusion factor), предложенный впервые в работах S. Lieblein. Согласно Lieblein аэродинамическая нагруженность межлопаточ-ной решетки определяется распределением давления по профилю лопатки и его влиянием на отрыв потока. В дозвуковой области режимов работы основным фактором, влияющим на отрыв пограничного слоя, является градиент давления на профиле лопатки. Более удобным для использования по сравнению с распределением давления является распределение скоростей по профилю лопатки. Распределение скоростей по спинке лопатки является основанием для определения потерь в межлопаточном канале и ограничения аэродинамической нагруженности лопаточного венца.

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

В наиболее простой форме степень диффузорности определяется для ротора в следующей форме:

w

D =_т« ,

^2

где мтах - максимальная местная скорость на спинке лопатки; ^2 - среднее значение скорости на выходе из межлопаточного канала.

Альтернативные выражения для вычисления степени диффузорности были предложены Koch&Smith:

^ г ■ с. — г ■ с.

= 1 —^ +

w,

(ri + r2 )•

) • W1

'R

Koch&Smith описали способ учета в выражении коэффициента диффузорности влияния сужения струйки тока (которое соответствует отношению осевых плотностей тока), а также изменения радиуса струйки тока между входом и выходом.

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

Скорости, приведённые в формулах, представлены на треугольниках скоростей (рис. 1).

Рис. 1. Треугольник скоростей и решётка профилей рабочего колеса (ротора)

Существует и альтернативный подход к определению коэффициента диффузорности - так называемый модифицированный коэффициент диффузорности (mofified diffusion factor), предложенный Schoibeiri [2]. Модифицированный коэффициент диффузорности позволяет учесть изменение осевой плотности тока в лопаточном венце (за счет изменения осевой скорости потока), изменение окружной скорости в лопаточном венце (за счет изменения радиуса струйки тока), влияние сжимаемости потока в межлопаточ-ном канале (изменение чисел Маха и образо-

вание местных сверхзвуковых областей потока на спинке профиля).

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

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

к

СТ* _ ^к-1

Awi

Xr =

w2

t*

V wl 0

w2

Ръ1 — Р1

В упрощённой трактов ке полные потери складываются из профильных и волновых потерь. Зависимость профильных потерь от коэффициента диффузорности и относитель-

(

ной высоты лопатки

Hr =

ГК ГВТ 0

приве-

дена на рис. 2. В литературе [2] аналогичная зависимость указывается для определения полного коэффициента потерь, однако данное утверждение может быть справедливо лишь для дозвуковых режимов работы лопаточного венца.

■о-Нг=0,0 -0-НЫ).1 и 1.0 -А-Нг=0,2 и 0,9 -О-Иг-0.3 и 0.8 -ш-Нг=0.4-0 7

Рис. 2. Распределение коэффициента профильных потерь по высоте пера лопатки в зависимости от коэффициента диффузорности (ротор)

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

модеиствием потоков, сходящих со спинки и корытца лопатки, образованием вихревой структуры течения за профилем.

Резкий рост коэффициента потерь полного давления при значительном повышении коэффициента диффузорности (как правило, более 0,5-0,6) сопровождается отрывом пограничного слоя со спинки лопатки в связи со значительным повышением толщины потери импульса. Отличное от нуля значение коэффициента потерь давления при нулевом значении коэффициента диффузорности определяется вязкостным трением (наличием касательных напряжений в потоке) и в меньшей степени ненулевой толщиной выходной кромки лопатки.

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

^рЯ^)™ Ь2

(Re). Разработка критерия

2s

2s

= a + eb' Dr+c+d'Dr + f • Dr .

Кроме профильных потерь в лопаточном венце имеются волновые потери. Волновые потери возникают при наличии околозвуковых и сверхзвуковых областей течения в межлопаточном канале. На рис. 3 приведена зависимость коэффициента волновых потерь от числа Маха на входе в решётку.

В качестве граничного значения для определения волновых потерь в решетке принято значениеМк1 > 0.8 (Мс2 > 0.8).

Рассмотрение волновых потерь при дозвуковой скорости потока на входе в межлопаточный канал связано с образова-

нием сверхзвуковой области течения на спинке профиля, что приводит к образованию скачков уплотнения, вызывающих потери давления.

«,(!)■ «ІЯЯЛ

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

Для ротора и статора предложена следующая аппроксимирующая зависимость вида:

ХрЩЕ)5™ Ь2

0,8 0.85 0,9 0.95 1 1,05 1.1 1.15 1.2 1.25 1,3 1.35 1.4 1,45 1,5 1,55

о эксперимент —аппроксимирующая кривая (м 2 у

Рис. 3. Зависимость коэффициента волновых потерь от числа Маха на входе в решётку

Взаимодействие скачков уплотнения с пограничным слоем является существенным фактором, приводящим к отрыву

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

(39.174832 + 28.219423- Ыы

v+ 0.710458 • Miw1

(51.603992 + 44.186608 • Mw1

v- 3.356873 • M21

- 0.732542.

Для статора разработана зависимость, аналогичная с заменойMw1 на Mc2 .

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

Х = Xs + Xp •

Для решения прямой задачи для осевого компрессора в двухмерной постановке авторами разработан алгоритм, реализованный в СИМ COMPRESSOR 2D.

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

Для верификации разработанной методики была выбрана аэродинамически высо-конагруженная ступень компрессора, типичная для первой ступени многоступенчатого вентилятора двухконтурного турбореактивного двигателя с низкой степенью двухкон-турности [3]. Основные параметры ступени в проектной точке приведены в табл. 1.

Таблица 1. Основные параметры ступени

Параметр Значение

GВ, кг/с 84,867

п, об/мин 11110

* ПРК 2,000

Р Лрк 0,887

Р ПСТ 1,936

Р п ст- 0,842

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

обеспечивающего осевой выход потока. На рис. 4 приведена характеристика рабочего колеса с положением расчетных точек. На рис. 5 приведена структурная схема рабоче-

го колеса при работе в СИМ СОМРЯЕЗЗСЖ 2Б [4].

Рис. 4. Характеристика рабочего колеса

Из особенностей системы моделирования С0МРКЕ880К_2Б можно отметить возможность задания различных типов аэродинамического профиля, возможность задания различных моделей профильных и волновых потерь в лопаточных венцах, в том числе индивидуальных альтернативных моделей по результатам обработки результатов испытаний экспериментальных ступеней-прототипов, возможность задания произвольного числа струек тока.

Некоторые параметры задаются в виде файла распределения по радиусу. В алгоритме системы предусмотрено осреднение параметров по двум различным методикам. После проведения расчетов создается набор файлов, содержащих в себе распределение основных параметров по высоте пера лопатки.

Расчеты экспериментальной ступени проводились в трех точках ветви характеристики (ппр = 100% ) согласно рис. 4.

Точка 1 (зона устойчивой работы) на характеристике обладает следующими параметрами:

GВ = 83,06 кг/с; п = 11088 об/мин;

* *

7ТРК = 1,8926; ^РК = 0,8823.

На рис. 6,7 представлены радиальные расчетные распределения степени повышения давления и КПД в сравнении с результатами испытаний ступени. На рис. 8 представлено радиальное распределение углов атаки, полученных расчетным путем и экспериментально с нанесением расчетных границ зарождения срывных зон и развитого срыва со спинки профиля.

Нг

9 кгаерим вит *о* расч т

Рис. 6. Радиальное распределение степени повышения полного давления в РК

»./Г. П.Н І'.НГі 0,4 П.чг. Лрк

«•"экспери мент -о-р асчет

Рис. 7. Радиальное распределение полного адиабатического КПД в РК

ных зон и развитого срыва предложены соотношения [5].

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

Точка 2 (граница устойчивой работы) на характеристике обладает следующими параметрами:

GВ = 78,79 кг/с; п = 11095 об/мин;

ярК = 2,0231; ^РК = 0,8573.

На рис. 9-10 представлены радиальные расчетные распределения степени повышения давления и КПД в сравнении с экспериментом. На рис. 11 представлено радиальное распределение углов атаки, полученных расчетным путем и экспериментально с нанесением расчетных границ зарождения срывных зон и развитого срыва.

2.ІН ■:,! ¿ля ТГр,

-«-:н.чіг-ртміт -оря.'чгт

Рис. 9. Радиальное распределение степени повышения давления в РК

Рис. 8. Радиальное распределение углов атаки в РК

Авторами предложена методика определения границы устойчивой работы лопаточного венца ступени компрессора. Для оценки устойчивости рабочего процесса в компрессоре предлагаются два безразмерных комплекса. Для границы зарождения срыв-

0,65 0,7 0,75 0,8 0,85 0,9 0,95 Лрк

-♦-эксперимент -о-расчет

Рис. 10. Радиальное распределение КПД в РК

-•-экспери мент

Рис. 12. Радиальное распределение степени повышения давления в РК

0,6 0,65 0,7 0,75 0,8 0,05 0,9 0,95 Лрк

эксперимент -о-расчет

Рис. 11. Радиальное распределение углов атаки в РК

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

Точка 3 (зона запирания) на характеристике обладает следующими параметрами:

GВ = 83,61 кг/с; п = 11095 об/мин;

ярК = 1,7771; ^РК = 0,8534.

На рис. 12-13 представлены радиальные расчетные распределения степени повышения давления и КПД в сравнении с экспериментом.

На рис. 14 представлено радиальное распределение углов атаки, полученных при расчете и в эксперименте с нанесением расчетных границ начала срыва и срыва.

Рис. 13. Радиальное распределение КПД в РК

Рис. 14. Радиальное распределение углов атаки в РК Выводы

1. Разработана методика решения прямой задачи в осевом компрессоре в двухмерной постановке.

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

3. Проверена адекватность разработанной системы и методики.

4. Разработана методика оценки запасов газодинамической устойчивости компрессора в двухмерной постановке.

Библиографический список

1. Johnsen, I.A. Aerodynamic design of axial-flow compressors. NASA SP-36 Report.3 / I.A. Johnsen, R.O. Bullock. - Cleveland : NASA-Lewis Research Center, 1965. - 526 p.

2. Schobeiri, M. Turbomachinery Flow Physics and Dynamic Performance / M. Schobeiri. - Berlin : Springer-Verlag Berlin Heidelberg, 2005. - 522 p.

3. Sulam, D.H. Single stage evalution of highly-loaded high-Mach-number compressor stages. NASA technical report CR-72694 / D.H. Sulam, M.J. Keenan, J.T. Flynn. - Maryland : Pratt&Whitney Aircraft Division (PWA-3772), 1974. - 108 p.

4. Козловская, А.Б. Система моделирования компрессоров авиационных ГТД (COMPRESSOR_2D) № 2010611711 / А.Б. Козловская, ДА. Ахмедзянов, А.Е. Михайлов. - М.: Роспатент, 2010.

5. Расчет характеристик осевых

компрессоров совместно с границей устойчивой работы [Текст] / А.Б.

Козловская, Д. А. Ахмедзянов [и др.]: Вестник УГАТУ - Уфа : УГАТУ, 2010. Т. 14

- №5 (40). - С. 15-26.

TWO-DIMENTIONAL AXIAL ROTOR COMPRESSOR OFF-DESIGN PERFORMANCE PREDICTION MODEL WITH THE USE OF SIMULATION METHODS AND TOOLS

© 2011 A. B. Mikhailova, D. A. Akhmedzyanov, Yu. M. Akhmetov, A. Ye. Mikhailov

Ufa State Aviation Technical University

The method of solving the two-dimensional direct problem for axial compressor is developed and realized in the simulation system COMPRESSOR_2D. The study on adequacy of proposed all-inclusive criterion for surge zones emergence prediction is executed.

Axial compressor, simulation, compressor performance, surge line, stall, incidence and deviation angles

Информация об авторах

Михайлова Александра Борисовна, аспирант Уфимского государственного авиационного технического университета. E-mail: [email protected]. Область научных интересов: осевые и центробежные компрессоры, расчет характеристик осевых компрессоров, имитационное моделирование, устойчивость.

Ахмедзянов Дмитрий Альбертович, доктор технических наук, профессор Уфимского государственного авиационного технического университета. E-mail: [email protected]. Область научных интересов: авиационные двигатели, рабочие процессы, динамика ГТД, автоматика и регулирование, имитационное моделирование.

Ахметов Юрий Мавлютович, кандидат технических наук, доцент Уфимского государственного авиационного технического университета. Область научных интересов: устойчивость, динамика ГТД, автоматика и регулирование, имитационное моделирование.

Михайлов Алексей Евгеньевич, аспирант Уфимского государственного авиационного технического университета. E-mail: [email protected]. Область научных интересов: рабочие процессы, устойчивость, динамика ГТД, автоматика и регулирование, имитационное моделирование.

Mikhailova Alexandra Borisovna, post-graduate student of Ufa State Aviation Technical University. E-mail: [email protected]. Area of research: axial and centrifugal compressors, off-design performance prediction, simulation, gas-dynamic stability.

Akhmedzyanov Dmitry Albertovich, doctor of technical science, professor, of Ufa State Aviation Technical University. E-mail: [email protected]. Area of research: aircraft engines, working processes, dynamics of GTE, controlling and automatics, simulation.

Akhmetov Yury Mavlyutovich, candidate of technical science, associate professor of Ufa State Aviation Technical University. Area of research: gas-dynamic stability, dynamics of GTE, controlling and automatics, simulation.

Mikhailov Alexey Yevgenevich, post-graduate student of Ufa State Aviation Technical University. E-mail: [email protected]. Area of research: working processes, gas-dynamic stability, dynamics of GTE, controlling and automatics, simulation.