ВВЕДЕНИЕ 4
ГЛАВА 1. АНАЛИЗ НАУЧНО-ТЕХНИЧЕСКОЙ ИНФОРМАЦИИ И ПОСТАНОВКА ЗАДАЧ ИССЛЕДОВАНИЯ 10
1.1. Влияние углового отрыва на производительность ОК и согласование ступеней 10
1.2. Теоретическое описание трехмерного отрыва потока от гладкой стенки 12
1.3. Обзор мирового опыта в области исследований трехмерного углового отрыва
в осевых турбомашинах 18
1.4. Обзор исследований угловых отрывов в осевых турбомашинах
с помощью CFD методов 19
1.5. Способы визуализаций течений в турбомашинах 20
1.6. Влияние радиального зазора, распределения входных параметров и относительного
удлинения венца на производительность, загромождение межлопаточного канала и образование углового отрыва 23
1.7. Способы аэродинамической доводки лопаток осевых компрессоров 29
1.8. Выводы по первой главе 43
1.9. Постановка задач исследования 44
ГЛАВА 2. ПОСТАНОВКА ЗАДАЧИ ЧИСЛЕННОГО ИССЛЕДОВАНИЯ ТЕЧЕНИЯ В ПРОТОЧНОЙ ЧАСТИ ОСЕВОГО КОМПРЕССОРА 46
2.2. Анализ научно-исследовательской литературы по верификации
течений в осевых турбомашинах 55
2.3. Выводы, основанные на мировом опыте верификации CFD моделей осевых
компрессоров и турбин 65
2.4. Верификация течения потока в двухступенчатом осевом
низкоскоростном компрессоре 66
2.5. Выводы по второй главе 84
ГЛАВА 3. РАЗРАБОТКА КРИТЕРИЕВ ДЛЯ ИДЕНТИФИКАЦИИ УГЛОВОГО ОТРЫВА ПОТОКА В МЕЖЛОПАТОЧНОМ КАНАЛЕ ОК И ЕГО УСТРАНЕНИЕ С ПОМОЩЬЮ ТАНГЕНЦИАЛЬНОГО НАВАЛА 85
3.2. Проверка критериев для идентификации углового отрыва потока на
направляющем аппарате второй ступени исследуемого ОК 95
3.3. Исследование влияния тангенциального навала на снижение размеров
углового отрыва потока в направляющих аппаратах ОК 97
3.4. Определение эффективности использования лопаток с тангенциальным навалом
для расширения диапазона устойчивой работы ОК 104
3.5. Выводы по третьей главе 105
ГЛАВА 4. РАСЧЕТНОЕ И ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ТЕЧЕНИЯ В ОСЕВОМ КОМПРЕССОРЕ НИЗКОГО ДАВЛЕНИЯ ГАЗОТУРБИННОГО ДВИГАТЕЛЯ ТИПА ДН80Л1 107
4.2. CFD анализ потока в венцах многоступенчатого осевого КНД установки
ДН80Л1 с использованием параметров WSa и Т1а 110
4.3. Оптимизация направляющего аппарата второй ступени КНД с целью снижения
размеров углового отрыва потока 113
4.4. Экспериментальное подтверждение результатов оптимизации
направляющего аппарата второй ступени 116
4.5. Выводы по четвертой главе 120
ЗАКЛЮЧЕНИЕ 122
СПИСОК СОКРАЩЕНИЙ И УСЛОВНЫХ ОБОЗНАЧЕНИЙ 123
СПИСОК ЛИТЕРАТУРЫ 125
ПРИЛОЖЕНИЕ
Актуальность темы. Газотурбинные двигатели (ГТД) широко применяются в авиации, энергетике, корабельных силовых установках, в качестве привода центробежных компрессоров газоперекачивающих агрегатов, а также в технологических циклах некоторых химических производств.
Основными требованиями ко всем ГТД, независимо от их назначения, являются высокая эффективность, эксплуатационная надежность и длительный ресурс. Одним из ключевых факторов, определяющих эффективность работы ГТД, является качество проектирования проточной части осевого компрессора (ОК) и турбины. Поскольку течение потока в осевых турбомашинах в общем случае является пространственным, большое значение имеет рассмотрение трехмерной модели течения в межлопаточных каналах их проточных частей.
Важным узлом в ГТД является компрессор и проектирование его лопаточного аппарата напрямую влияет на эффективность работы установки, а в частности на КПД и на топливную эффективность. В связи с чем повышения эффективности лопаточного аппарата (ЛА) осевого компрессора возможно добиться с помощью создания аэродинамически оптимальных форм перьевой части рабочих и направляющих лопаток компрессорных ступеней при учете пространственного характера течения в них.
Распространенным явлением, которое возникает в ОК с увеличением степени повышения давления в ЛА, является угловой отрыв потока. Он возникает в области перехода корневой полки в перо лопатки, приводя к загромождению межлопаточного канала и снижению запаса газодинамической устойчивости (ГДУ) ОК. Для того чтобы снизить размер углового отрыва возможно выполнить аэродинамическую оптимизацию пера лопатки с помощью тангенциального навала. Данный способ позволяет повысить запас ГДУ ЛА за счет снижения градиентов давления и увеличения значения касательных напряжений на поверхности лопаток у корневой полки.
В первых ступенях осевых компрессоров современных газотурбинных двигателей вероятность возникновения углового отрыва потока при согласованной работе ступеней крайне мала, при этом толщина пограничного слоя на торцевых поверхностях вследствие торможения потока об ограничивающие поверхности составляет 1 -2% от высоты межлопаточного канала, что не вносит существенных изменений в структуру потока, а отклонение фактических треугольников скоростей от расчетных является незначительным. По мере продвижения потока к промежуточным и последним ступеням происходит существенное нарастание торцевого пограничного слоя, что связано с увеличением углов атаки на лопатки у меридиональных обводов. При этом величина данного пограничного слоя к последней ступени многоступенчатого ОК может достигать 20% от высоты канала. Эффект углового отрыва вносит нестационарность в поток перед венцом и снижает скорость потока в торцевых областях лопатки, увеличивая размер торцевого погранслоя.
Сложность замера полей параметров потока в межлопаточном канале для определения градиентов давления и скоростей, особенно у поверхности лопатки, позволяющих установить пространственный характер течения газа, заставлет все больше уделять внимание методам вычислительной газовой динамики. Использование численного метода для решения задач газовой динамики позволяет количественно и качественно установить вихревые потоки, что в экспериментальном исследовании, в данном случае, является невозможным.
Актуальность работы, определяющая цели и задачи исследования, заключается в необходимости совершенствования стационарных и транспортных ГТУ, в том числе путем разработки и реализации мероприятий в области проектирования, доводки и модернизации проточной части ОК с применением современных программных комплексов.
Степень разработанности темы исследования. Исследованиями механизма образования трехмерного углового отрыва в межлопаточном канале и влиянием его на производительность лопаточного венца занимаются порядка 50 лет. Проведенные исследования в области контроля отделения потока у торцов лопатки позволили установить, что в потоке присутствуют отрицательные градиенты давления, которые отвечает за отрыв пограничного слоя, а сам угловой отрыв является эффектом совместного взаимодействия обратного градиента давления и вторичных вихрей. Полученные результаты предыдущих исследований, дают возможность сконцентрировать свое внимание на структуре потока и способах идентификации углового отрыва в межлопаточном канале.
Современные способы исследования движения потоков газа, которые основаны на методах вычислительной газовой динамики, позволяют более детально изучить поток в области углового отрыва, а в частности, области критических точек и предельных линий тока, что дает возможность спрогнозировать работу лопаточного венца в зависимости от течения потока в межлопаточном канале.
Цель работы: Разработка критериев для идентификации углового отрыва потока в лопаточных венцах при решении задач численного моделирования течения в проточных частях осевых компрессоров ГТД и способа его устранения при проектировании новых и модернизации существующих конструкций. Для достижения указанной цели поставлены и решены следующие задачи:
• анализ научно-технической литературы, связанной с исследованиями трехмерных
течений в компрессорных решетках, а также анализ и систематизация основных подходов к повышению эффективности компрессорных решеток с помощью изменения пространственной формы пера лопатки;
• моделирование работы ОК на основе численного решения системы уравнений Навье-Стокса в трехмерной постановке; верификация расчетной модели по данным экспериментальных исследований для выбора оптимальных настроек CFD кода в данном исследовании;
• разработка критериев для идентификации углового отрыва потока в межлопаточном канале и определения режима оптимизации лопаточного венца с помощью тангенциального навала;
• практическое применение разработанных критериев определения углового отрыва потока при аэродинамическом совершенствовании проточной части ОК натурной ГТУ.
Научная новизна работы определяется тем, что:
1. Установлено значение коэффициента для определения вихревой вязкости потока при CFD расчете, обеспечивающее наилучшее совпадение результатов CFD расчета с экспериментальными данными на различных режимах работы ОК;
2. Разработаны критерии для идентификации углового отрыва в межлопаточном канале ОК при решении задачи численного моделирования;
3. Разработан критерий для определения размеров углового отрыва;
4. Предложен способ определения условий работы венца на основе разработанных критериев идентификации углового отрыва.;
5. Предложен и физически обоснован подход к параметрическому проектированию формы линии совмещения плоских профилей лопаточного венца ОК в тангенциальном направлении с помощью кривой Безье третьего порядка;
6. Разработаны критерии для выбора режима работы лопаточного венца, на котором должна проводиться оптимизация формы ЛА с помощью тангенциального навала для достижения максимального эффекта;
Методология и методы диссертационного исследования. В работе использованы современные методы вычислительной газовой динамики, которые позволяют изучить структуру вихрей в угловом отрыве потока и идентифицировать образование данного явления в межлопаточном канале лопаточного венца осевого компрессора. Также в работе представлены экспериментальные методы исследования потока, с помощь которых идентифицировано наличие углового отрыва потока в проточной части осевого компрессора натурного газотурбинного двигателя.
Теоретическая и практическая значимость работы для отрасли турбостроения состоит в следующем:
• Сформулированы рекомендации для стационарных RANS методов, которые позволяют приблизить результаты CFD расчета низкоскоростного ОК к эксперименту на режимах работы с устойчивой вихревой структурой в потоке;
• Сформулированы рекомендации по идентификации углового отрыва и определения его размеров с использованием разработанных критериев при решении задачи численного моделирования для различных режимов работы ОК;
• Предложены рекомендации по выбору режима для оптимизации лопаточного венца с целью повышения его аэродинамической эффективности во всем диапазоне режимов работы ОК;
• Разработана и впервые апробирована на натурной ГТУ методика определения степени распространения углового отрыва, как устойчивой структуры потока.
• Выработаны рекомендации по оптимизации направляющего аппарата второй ступени с помощью введения тангенциального навала и изменения угла установки лопаток для снижения размеров углового отрыва потока в компрессоре низкого давления (КНД) ГТД типа ДН80Л1 и предложены рекомендации для линии совмещения плоских профилей, а также рекомендации по выбору наиболее эффективного способа оптимизации венцов осевого компрессора.
Положения, выносимые на защиту:
1. Разработанные рекомендации по повышению точности CFD расчета в стационарном RANS методе;
2. Критерии для идентификации углового отрыва в межлопаточном канале осевого компрессора;
3. Разработанный способ построения линии совмещения плоских профилей по высоте лопатки;
4. Критерий выбора режима для оптимизации лопаточного аппарата при анализе результатов численного моделирования работы ОК с целью снижения размеров углового отрыва потока;
5. Результаты проведенных численных и оптимизационных исследований;
6. Результаты апробации данных исследования при аэродинамической оптимизации ОК натурного ГТД типа ДН80Л1.
Степень достоверности результатов определяется:
• использованием в процессе выполнения работы наиболее современных апробированных и научно обоснованных программ, методик численного трехмерного расчета течений потока газа в лопаточных аппаратах турбомашин;
• высокой точностью совпадения результатов численного моделирования течения с экспериментальными данными и физическим обоснованием полученных расхождений в каждом конкретном случае;
• промышленной апробацией разработанных принципов при оптимизации направляющего аппарата осевого компрессора низкого давления натурной ГТД типа ДН80Л1;
• метрологическим обеспечением и точностью измерительной аппаратуры для обработки результатов экспериментальных исследований.
1. Предложены и обоснованы внесенные коэффициенты для определения вихревой вязкости, которые могут быть применены с целью улучшения совпадения результатов CFD расчета с экспериментом для установки качественной и количественной картины течения в осевом компрессоре с низкими числами Маха. Установлено, что наилучшее совпадение результатов CFD расчета и эксперимента достигается при влиянии на вихревую вязкость с помощью данного коэффициента, который может варьироваться в диапазоне от 0.27 до 0.4;
2. Разработаны два критерия (WSa и TIa) для идентификации углового отрыва в межлопаточном канале ОК на этапе CFD моделирования: Первый критерий идентифицирует отрыв потока от поверхности лопатки и показывает, степень удаленности режима работы венца от своей границы ГДУ; Второй критерий позволяет определить размеры области всего вихревого движения потока в межлопаточном канале. В результате исследования выделены критические значения для критерия №1 и №2, которые находятся в диапазоне от 1.6 до 1.8 и от 4 до 7% соответственно. При этом получена диаграмма для определения работоспособности венца, которая отражает режим его работы на основе загромождения канала и степени устойчивости венца к отрыву потока;
3. Предложена и подтверждена эффективность использования линии совмещения плоских профилей по высоте лопатки с помощью кривой Безье 3 порядка. Для проверки эффективности применения тангенциального навала использованы установленные критерии №1 и №2. Определен максимальный прогиб пера лопатки, при котором достигается максимальная эффективность с точки зрения снижения размеров углового отрыва. В результате применения тангенциального навала для направляющего аппарата исследуемого низкоскоростного двухступенчатого ОК увеличился запас газодинамической устойчивости на 5%;
4. Реализован и апробирован способ идентификации углового отрыва потока в межлопаточном канале при анализе результатов CFD расчета натурного ОК ГТД типа ДН80Л1. Для снижения размеров углового отрыва выполнена двухэтапная оптимизация направляющего лопаточного аппарата второй ступени в CFD программе с помощью тангенциального навала и с помощью изменения угла установки всего венца. На натурном ОК реализовано изменение угла установки всего венца НА2 в качестве первого этапа оптимизации данной установки.
Перспективы дальнейшей разработки темы исследования связаны с расширением диапазона исследуемых граничных условий на входе и выходите расчетной модели лопаточного венца компрессорной ступени; с разработкой методики оптимизации профиля пера лопатки с помощью осевого навала.
1. Абдельвахид, М. Б. Расчёт влияния радиального зазора на характеристику кнд двигателя РД-33/Абдельвахид, М. Б., Фёдоров, Р. М., Федечкин К. С.// Вестник Самарского государственного аэрокосмического университета им. академика С.П. Королёва (национального исследовательского университета) - 2012. - 23-27 cc.
2. Ануров Ю.М. Численный анализ возможностей 2-D и 3-D методов проектирования осевых турбомашин / Ю. М. Ануров, С. О. Беляева, В. А. Коваль // Восточно-Европейский журнал передовых технологий. - 2009. вып. 4.
3. Баранник В. С Пространственная аэродинамическая оптимизация направляющей решетки осевой турбины: диссертация на соискание ученой степени кандидата технических наук: 05.05.16 / Баранник Валентин Сергеевич. - Харьков, 2016.
4. Батурин О.В. Расчетное изучение рабочего процесса в ступени компрессора, спроектированной при курсовом проектировании по курсу «Теория и расчет лопаточных машин»: учеб. пособие/ О.В. Батурин и др. - Самара: Изд-во Самар. гос. аэрокосм. ун-та,
2011. - 94 с.
5. Белов, И.А. Моделирование турбулентных течений: учебное пособие / И.А. Белов, С.А. Исаев. - СПб.: Балт. Гос. Техн. Ун-т., 2001. - 108 с.
6. Блинов, В. Л. Разработка принципов параметрического профилирования плоских решеток осевых компрессоров ГТУ на основании результатов многокритериальной оптимизации: диссертация на соискание ученой степени кандидата технических наук: 05.04.12 / Блинов Виталий Леонидович. - Екатеринбург, 2015.
7. Бродов Ю. М. Выяснение причин недостаточной прочности лопаток компрессора низкого давления ГТУ типа ДН80, ДУ80 и выработка мероприятий по их устранению: отчет о НИОКР / ООО УНН11 «Теплоэнергетик» ; рук. Ю. М. Бродов ; исполн.: О. В. Комаров, А. В. Скороходов, В. А. Седунин, В. Л. Блинов. Екатеринбург 2016. 110 с.
8. Галаев, С.А. Численное моделирование течения вязкого газа в решетках осевых турбомашин: методика и результаты применения современных программных средств: автореферат диссертации на соискание ученой степени кандидата технических наук: 05.04.12 / Галаев Сергей Александрович. - Санкт-Петербург, 2006.
9. Галеркин, Ю.Б. Турбокомпрессоры: учеб. пособие / Ю.Б. Галеркин, Л.И. Козаченко. - СПб.: Изд-во Политехн. ун-та, 2008. - 374 с.Про Fd
10. .Галеркин, Ю.Б. Анализ эффективности пространственных лопаточных решеток осевых компрессоров по данным продувок плоских решеток / Ю.Б. Галеркин, Ю.А. Попов // Компрессорная техника и пневматика. - 2005. - №3. 158
11. Галеркин, Ю.Б. Оптимизация проточной части осевых компрессоров на стадии вариантного расчета. Часть 1 / Ю.Б. Галеркин, Ю.А. Попов // Компрессорная техника и пневматика. - 2009. - №5.
12. Галеркин, Ю.Б. Расчетный анализ характеристик осевых компрессорных ступеней / Ю.Б. Галеркин, Ю.А. Попов // Компрессорная техника и пневматика. - 2005. - №5.
13. Гарбарук А.В. Моделирование турбулентности в расчетах сложных течений: учебное пособие / А.В. Гарбарук, М.Х. Стрелец, М.Л. Шур - СПб: Изд-во Политехн. ун-та, 2012. - 88 с.
14. Гостелоу, Дж. Аэродинамика решеток турбомашин: пер. с англ. Н.М. Савина, А.П. Кадетова. Под ред. В.Т, Митрохина / Дж. Гостелоу. - М.: Мир, 1987. - 391 с.Про Fd
15. Дейч, М.Е. Техническая газодинамика / М.Е. Дейч. - М.: Госэнергоатомиздат, 1961. - 669 с.
16. Иноземцев, А.А. Газотурбинные двигатели / А.А. Иноземцев, В.Л. Сандрацкий. - Пермь: Изд-во ОАО «Авиадвигатель», 2006 - 1204 с.
17. Ковалева, Е.А. Влияние радиального зазора на границу срыва осевых компрессоров/ Ковалева Е.А., Коваль В.А., Романов В.В.// НАСОСЫ. ТУРБИНЫ. СИСТЕМЫ - 2013. - №3(8) - 27-34 cc.
18. Комаров, А.П. Аэродинамические характеристики 133 компрессорных решеток, составленных из профилей одного семейства (при малых скоростях набегающего потока) / А.П. Комаров. - М.: ЦИАМ, 1955. - 79 с.
19. Комаров, А.П. Исследование плоских компрессорных решеток / А.П. Комаров // Лопаточные машины и струйные аппараты. Сб. статей. - М.: Машиностроение, 1967. - Вып. 2. - С. 67-110.
20. Лапин, Ю.В. Внутренние течения газовых смесей / Ю.В. Лапин, М.Х. Стрелец. - М.: Наука, 1989. - 356с.
21. Самохвалов, Н.Ю. Верификация расчетных трехмерных моделей образования вторичных вихревых течений в межлопаточных каналах турбин газотурбинных двигателей по данным испытаний тестового профиля Т106 / Н.Ю. Самохвалов, В.Т. Хайрулин, А.С. Тихонов // Вестник ПНИПУ. Аэрокосмическая техника. - 2017. - вып. 49.
22. Плотников, А. И. Определение оптимальных и критических зависимостей параметров ступеней осевого компрессора для повышения запасов устойчивой работы / А. А. Иноземцев, Е. Т. Гузачев, В. А. Катаев, А. И. Плотников // Вестник Пермского национального исследовательского политехнического университета. Аэрокосмическая техника. - 2015. - №42 - 5-19 cc.
23. Прокопец, А.О. Повышение эффективности работы приводных стационарных газотурбинных установок в условиях эксплуатации ООО «Газпром трансгаз Югорск»: диссертация на соискание ученой степени кандидата технических наук: 05.04.12 / Прокопец Алексей Олегович. - Екатеринбург, 2012.
24. Разработка, исследование и доводка ГТУ, компрессоров, нагнетателей и их элементов: сб. науч. тр. / НПО по исследованию и проектированию энер. оборудования им. И.И. Ползунова // Труды ЦКТИ. - 1990. - Вып. 270.
25. Рассохин, В.А. Расчет трехмерного течения в ступенях малорасходных турбин / А.А. Епифанов, А.И. Кирилов, В.А. Рассохин // Научно-технические ведомости СПбГПУ. -
2012. - №142. - С. 65-70.
26. Рассохин, В.А. Расчет характеристик лопаточных решеток малорасходных турбин / А.И. Кириллов, В.А. Рассохин, А.А. Епифанов // Научно-технические ведомости СПбГПУ. - 2012. - №1(142).
27. Рассохин, В.А. Устойчивость потока в рабочем колесе при различных условиях на входе в ступень / А.В. Захаров, Г.Г. Шпензер, В.А. Рассохин // Научно-технические ведомости СПбГПУ. - 2013. - №3(178). - С. 28-31.
28. Ревзин, Б.С. Осевые компрессоры газотурбинных газоперекачивающих агрегатов: учебное пособие / Б.С. Ревзин. - 2-е изд., стер. - Екатеринбург: Изд-во УГТУ-УПИ, 2000. - 90 с
29. Роджерс Д. Математические основы машинной графики / Д. Роджерс, Дж. Адамс. - М.: Мир, 2001. - 604 с.
30. Роуч, П. «Вычислительная гидродинамика/ Роуч, П.// М.: Мир, 1980. -618 с
31. Седунин, В.А. Исследование и разработка методов повышения эффективности работы первой ступени осевого компрессора ГТУ с регулируемым входным направляющим аппаратом: диссертация на соискание ученой степени кандидата технических наук: 05.04.12 / Седунин Вячеслав Алексеевич. - Екатеринбург, 2011
32. Серков, С. А. Особенности верификации СЕО-модели осевого компрессора. Часть 1: моделирование локальных эффектов и подбор параметров для совпадения с экспериментальными данными/ Ю. М. Бродов, О. В. Комаров, В. А. Седунин, В. Л. Блинов, С. А. Серков//Компрессорная техника и пневматика. - 2018. - № 4. - С. 4-28
33. Серков, С. А. Особенности верификации СЕО-модели осевого компрессора. Часть 2: анализ интегральных параметров / Ю. М. Бродов, О. В. Комаров, В. А. Седунин, В. Л. Блинов, С. А. Серков//Компрессорная техника и пневматика. - 2018. - № 4. - С. 29-28
34. Серков, С. А. Верификация задачи численного моделирования течения воздуха в осевой компрессорной ступени / О. В. Комаров, В. А. Седунин, В. Л. Блинов, С. А.
Серков//Вестник Московского государственного технического университета им. Н.Э. Баумана. Серия «Машиностроение. - 2016. - № 4. - С. 54-68
35. Серков, С. А. Optimization approach and some results for 2d compressor airfoil/ O. V. Komarov, V.A. Sedunin, V. L. Blinov, S. A. Serkov//International Journal of Gas Turbine, Propulsion and Power Systems December 2016, Volume 8, Number 3 С. 39-46
36. Серков, С. А. Расчётное исследование влияния закона закрутки и тангенциального навала при оптимизации лопаточного аппарата осевого компрессора / В. А. Седунин, В. Л. Блинов, С. А. Серков//Молодежный Вестник УГАТУ, №3(12), 2014, С. 86-93
37. Серков, С. А. Влияние тангенциального навала лопаток высоконагруженного осевого компрессора на его характеристики и запас газодинамической устойчивости / Ю. М. Бродов, О. В. Комаров, В. А. Седунин, В. Л. Блинов, С. А. Серков//Всероссийской научно¬технической конференция молодых ученых и специалистов «Новые решения и технологии в газотурбостроении» -2015- С. 34-38
38. Серков, С. А. Параметрическое исследование тангенциального навала лопаток осевого компрессора / В. А. Седунин, В. Л. Блинов, С. А. Серков //Энерго- и ресурсосбережение. Энергообеспечение. Нетрадиционные и возобновляемые источники энергии. -2015 - С. 17¬22