ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 06.04.2021
Просмотров: 1393
Скачиваний: 4
ЕНЕРГЕТИЧНІ
ТА
ТЕПЛОТЕХНІЧНІ
ПРОЦЕСИ
Й
УСТАТКУВАННЯ
............................... 143
............................ 150
............................................................................................ 154
........................................................ 168
.............................................................................................. 178
8’2012
5
ЕНЕРГЕТИЧНІ
ТА
ТЕПЛОТЕХНІЧНІ
ПРОЦЕСИ
Й
УСТАТКУВАННЯ
УДК
621.165
А
.
В
.
БОЙКО
,
д
-
р
техн
.
наук
;
проф
.
НТУ
«
ХПИ
»;
Ю
.
Н
.
ГОВОРУЩЕНКО
,
канд
.
техн
.
наук
;
с
.
н
.
с
.
НТУ
«
ХПИ
»;
М
.
В
.
БУРЛАКА
,
канд
.
техн
.
наук
;
н
.
с
.
НТУ
«
ХПИ
»;
В
.
С
.
БАРАННИК
,
инженер
2-
й
кат
.
НТУ
«
ХПИ
»
ПРОСТРАНСТВЕННАЯ
ОПТИМИЗАЦИЯ
ПРОФИЛЯ
ЛОПАТКИ
ПРЯМОЙ
ТУРБИННОЙ
РЕШЕТКИ
Приведена
постановка
задачи
исследования
с
описанием
варьированных
геометрических
параметров
профиля
и
методов
,
которые
использовались
.
Также
,
приведены
результаты
оптимизации
пространственного
профиля
направляющей
турбинной
решетки
с
использованием
различных
граничных
условий
на
торцевых
поверхностях
канала
.
Рассмотрено
методику
,
при
которой
достигается
стабилизация
потока
за
выходными
кромками
.
Приведена
постановка
задачі
дослідження
з
описанням
варійованих
геометричних
параметрів
профілю
та
методів
,
що
використовувалися
.
Також
,
приведені
результати
оптимізації
просторового
профілю
напрямної
турбінної
решітки
з
використанням
різних
граничних
умов
на
торцевих
поверхнях
каналу
.
Розглянуто
методику
,
при
якій
досягається
стабілізація
потоку
за
вихідними
кромками
.
The formulation of research task with the specification of the varied geometrical parameters of profile and
methods that is resulted. Also, the results of optimization of spatial profile of directing turbine grate with the
using of different maximum terms on the butt-end surfaces of channel are showed. The methodology at that
stabilizing of stream is arrived at after initial edges is considered.
Как
известно
,
характер
обтекания
турбинных
лопаток
в
ядре
потока
и
на
концах
решеток
существенно
отличается
.
В
то
же
время
постановка
задачи
оптимизации
турбинных
решеток
должна
учитывать
отличия
сложных
физических
процессов
,
имеющих
место
при
обтекании
решеток
профилей
по
высоте
лопаток
.
В
связи
с
выше
изложенным
,
задача
оптимизации
турбинных
решеток
с
функцией
цели
минимума
интегральных
потерь
должна
включать
возможность
применения
сложного
тангенциального
навала
решеток
с
изменением
формы
самого
профиля
по
высоте
лопатки
или
на
отдельных
ее
участках
(
у
корня
и
периферии
лопатки
).
Однако
оптимизация
лопатки
в
такой
постановке
представляет
сложную
задачу
и
требует
варьирования
большого
количества
параметров
и
,
следовательно
,
большого
количества
CFD
расчетов
.
Анализ
физических
процессов
,
имеющих
место
в
решетках
с
одновременным
изменением
многих
параметров
,
очень
сложен
и
не
позволяет
оценить
влияние
отдельных
параметров
на
потери
в
решетках
.
В
связи
с
этим
представляется
целесообразным
поэтапное
решение
указанной
задачи
:
1)
оптимизация
формы
профилей
в
двухмерной
постановке
;
2)
оптимизация
формы
профиля
прямой
решетки
в
трехмерной
постановке
;
3)
оптимизация
пространственной
формы
лопатки
направляющей
решетки
по
нескольким
сечениям
без
применения
сложного
тангенциального
навала
;
4)
оптимизация
сложного
тангенциального
навала
направляющей
решетки
совместно
с
формой
профиля
лопатки
по
нескольким
сечениям
.
По
первому
этапу
были
проведены
некоторые
исследования
,
представленные
в
работе
[1].
Данная
статья
посвящена
решению
второго
этапа
поставленной
оптимизационной
задачи
.
8’2012
6
ЕНЕРГЕТИЧНІ
ТА
ТЕПЛОТЕХНІЧНІ
ПРОЦЕСИ
Й
УСТАТКУВАННЯ
Постановка
задачи
оптимизации
Цель
исследования
–
повышение
аэродинамической
эффективности
профиля
лопатки
прямой
турбинной
решетки
при
заданных
параметрах
на
входе
и
выходе
с
сохранением
расхода
рабочего
тела
через
решетку
и
площади
профиля
.
Объектом
исследования
является
прямая
турбинная
решетка
по
своим
характеристикам
соответствующая
направляющему
аппарату
второй
ступени
цилиндра
высокого
давления
турбины
К
-310-23,5.
Высота
лопатки
равна
25,5
мм
,
площадь
профиля
– 0,000322
м
2
,
хорда
b
= 53,17
мм
,
относительный
шаг
t
/
b
= 0,8406,
длина
входного
участка
– 0,3
b
,
длина
выходного
участка
– 0,35
b
.
В
CFD
расчетах
использовалась
трехмерная
сетка
размерностью
166
×
72
×
62
на
один
межлопаточный
канал
.
Расчетный
объем
с
названиями
граничных
поверхностей
приведен
на
рис
. 1,
где
«pressure side»
и
«suction side»
это
соответственно
стороны
давления
и
разрежения
соседних
лопаток
в
турбинной
решетке
.
Построение
параметрической
модели
профиля
и
варьирование
параметров
производилось
с
помощью
,
созданной
на
кафедре
турбиностроения
,
программы
TOpGrid
.
Описание
параметрической
модели
профиля
и
входящие
в
нее
параметры
описаны
в
работе
[1].
При
оптимизации
варьировались
следующие
параметры
:
β
у
–
угол
установки
профиля
,
α
1
г
–
геометрический
угол
выхода
,
δ
–
угол
отгиба
,
ω
0
–
угол
заострения
входной
кромки
,
ω
1
–
угол
заострения
выходной
кромки
.
Функция
цели
–
суммарные
потери
(
сумма
профильных
и
концевых
).
При
расчетах
накладывались
ограничения
на
площадь
поперечного
сечения
лопатки
и
расход
рабочего
тела
.
В
качестве
минимально
допустимой
площади
поперечного
сечения
была
величина
площади
исходного
профиля
ТС
-1
А
,
а
максимальной
– 0,000386
м
2
(
увеличение
на
20
%).
Технологические
ограничения
в
данной
работе
не
рассматривались
.
Рис
. 1.
Граничные
условия
расчетной
области
Определение
оптимальной
формы
профиля
проводилось
с
помощью
теории
планирования
эксперимента
и
ЛП
τ
-
поиска
по
алгоритму
,
представленному
в
работе
[1].
Во
всех
расчетах
задавались
одни
и
те
же
граничные
условия
:
полное
давление
на
входе
– 16,15
МПа
;
полная
температура
на
входе
– 457,54
°
С
;
статическое
давление
на
выходе
– 13,72
МПа
;
степень
турбулентности
потока
на
входе
– 1 %;
модель
турбулентности
–
k
-
ω
SST
;
расчетная
схема
TVD
второго
порядка
аппроксимации
с
локально
первым
в
местах
больших
градиентов
;
y
+ < 1;
на
границах
расчетной
области
«hub», «shroud»
и
«blade»
задавалось
условие
стенки
с
прилипанием
потока
.
В
работе
[2]
по
ряду
причин
производилась
оптимизация
пространственной
формы
профиля
без
учета
кромочного
следа
,
так
как
выходная
плоскость
измерения
параметров
располагалась
непосредственно
за
выходной
кромкой
.
Однако
,
при
оптимизации
формы
профиля
прямой
турбинной
решетки
,
необходимо
учитывать
потери
,
возникающие
в
кромочном
следе
,
так
как
они
являются
неотъемлемой
составляющей
потерь
.
Для
определения
влияния
кромочных
потерь
на
результат
оптимизации
были
проведены
расчетные
исследования
в
двух
постановках
:
с
учетом
и
без
учета
кромочных
потерь
.
Влияние
кромочных
потерь
на
результаты
оптимизации
.
Результаты
расчетов
исходного
и
оптимальных
вариантов
приведены
в
таблице
1.
Как
видно
из
таблицы
,
кромочные
потери
оказывают
существенное
влияние
8’2012
7
ЕНЕРГЕТИЧНІ
ТА
ТЕПЛОТЕХНІЧНІ
ПРОЦЕСИ
Й
УСТАТКУВАННЯ
на
результат
оптимизации
.
Более
того
,
на
рис
. 2
можно
видеть
,
что
форма
оптимальных
профилей
,
полученных
при
различных
постановках
,
существенно
отличаются
друг
от
друга
.
То
есть
разница
в
значениях
суммарных
потерь
с
учетом
кромочных
потерь
и
без
их
учета
имеет
не
только
количественный
,
но
и
качественный
характер
.
Из
этого
следует
,
что
при
оптимизации
формы
профиля
по
суммарным
потерям
необходимо
учитывать
кромочные
потери
,
то
есть
располагать
плоскость
определения
выходных
параметров
на
некотором
расстоянии
от
выходных
кромок
.
Таблица
1
Результаты
оптимизации
С
учетом
кромочных
потерь
Без
учета
кромочных
потерь
Профиль
Параметры
Исходный
Оптимальный
Исходный
Оптимальный
Потери
, %
4,6033
4,48287
3,13186
3,09235
Расход
,
кг
/
с
2,08142
2,082
2,08142
2,08129
Площадь
,
м
2
0,000322 0,000384 0,000322 0,000386
В
ряде
работ
[3, 4]
по
экспериментальным
исследованиям
прямых
турбинных
решеток
рекомендуется
располагать
плоскость
измерения
выходных
параметров
на
расстоянии
1/7
хорды
от
выходных
кромок
решетки
или
величины
горла
.
Такие
рекомендации
обусловлены
тем
,
что
плоскость
замера
выходных
параметров
должна
располагаться
в
области
,
где
наступает
стабилизация
потока
и
кромочных
потерь
.
Однако
,
при
постановке
задачи
принятой
в
нашем
расчетном
исследовании
,
стабилизация
потока
не
происходит
вплоть
до
границы
расчетной
области
(
рис
. 3),
хотя
она
располагается
на
значительно
большем
расстоянии
,
чем
рекомендуется
в
указанных
работах
.
исходный
;
оптимальный
а
)
с
учетом
кромочных
потерь
б
)
без
учета
кромочных
потерь
Рис
. 2.
Формы
профилей
,
полученные
при
оптимизации
Рис
. 3.
Распределение
интегральных
потерь
за
выходными
кромками
Рис
. 4.
Структура
формирования
пограничного
слоя
на
торцевых
поверхностях
канала
Дать
объяснение
этого
эффекта
можно
,
если
рассмотреть
меридиональное
сечение
расчетного
канала
и
формирование
в
нем
пограничного
слоя
на
торцах
(
рис
. 4).
Из
рисунка
видно
,
что
на
торцевых
стенках
(«hub»
и
«shroud»)
имеет
место
постоянно
накапливающийся
пограничный
слой
,
формирующийся
от
начала
расчетного
канала
и
до
его
конца
.
Поэтому
потери
в
любом
поперечном
сечении
канала
будут
постоянно
увеличиваться
по
ходу
движения
потока
рабочего
тела
и
,
следовательно
,
стабилизация
потерь
,
имеющая
место
в
экспериментах
не
происходит
.
Очевидно
,
что
при
схеме
8’2012
8
ЕНЕРГЕТИЧНІ
ТА
ТЕПЛОТЕХНІЧНІ
ПРОЦЕСИ
Й
УСТАТКУВАННЯ
задания
граничных
условий
,
описанной
в
рассмотренной
постановке
задачи
,
стабилизация
потока
не
произойдет
в
не
зависимости
от
длины
выходного
участка
.
Рис
. 5.
Построения
расчетной
сетки
с
помощью
трех
объемов
Рис
. 6.
Распределение
интегральных
потерь
за
выходными
кромками
В
экспериментальных
исследованиях
прямых
турбинных
решеток
,
как
правило
,
используются
так
называемые
отсечные
пластины
,
отсекающие
накопленный
пограничный
слой
на
подводящих
к
исследуемой
решетке
каналах
.
В
свою
очередь
на
отсечных
пластинах
формируется
новый
пограничный
слой
.
Длина
этих
пластин
значительно
меньше
длины
границ
расчётной
области
«hub»
и
«shroud».
В
связи
с
выше
изложенным
потребовалось
изменить
схему
задания
граничных
условий
так
,
чтобы
пограничный
слой
накапливался
только
на
участках
границ
«hub»
и
«shroud»
примерно
соответствующим
размерам
отсечных
пластин
,
а
граница
«inlet»
и
плоскость
измерения
выходных
параметров
располагались
за
переделами
этих
участков
.
Кроме
этого
,
при
численных
расчетах
течений
в
турбинных
решетках
из
соображений
точности
и
сходимости
обычно
рекомендуется
,
чтобы
длина
входного
участка
была
не
менее
100 %
хорды
,
а
выходного
участка
не
менее
200 %
хорды
.
Таким
образом
,
была
принята
новая
схема
задания
граничных
условий
и
формирования
расчетной
области
(
рис
. 5),
удовлетворяющая
указанным
выше
требованиям
.
Как
видно
из
рис
. 5
новая
расчетная
область
состоит
из
трех
объемов
:
входной
(40
×
40
×
40),
собственно
решетка
(166
×
72
×
62)
и
выходной
(166
×
72
×
62).
На
границах
«back»
и
«front»
входного
и
выходного
расчетных
объемов
задавалось
условие
стенки
со
скольжением
потока
,
а
прилипание
на
этих
границах
только
в
объеме
турбинной
решетки
.
При
такой
схеме
задания
граничных
условий
,
пограничный
слой
формируется
только
в
объеме
турбинной
решетки
,
по
аналогии
с
отсечными
пластинами
в
эксперименте
.
Межлопаточный
канал
,
входной
и
выходной
объемы
были
построены
с
использованием
программы
TOpGrid
.
Характер
изменения
суммарных
потерь
в
зависимости
от
расстояния
за
решеткой
приведен
на
рис
. 6.
Как
видно
из
рисунка
,
в
новой
постановке
имеет
место
стабилизация
потерь
,
более
того
расстояние
стабилизации
в
достаточной
степени
согласуется
с
экспериментальными
данными
.
Таким
образом
,
для
оптимизации
формы
профилей
прямых
турбинных
решеток
необходимо
задавать
граничные
условия
как
представлено
на
рис
. 5.
Входной
участок
и
межлопаточный
канал
связывались
между
собой
через
границы
«inlet»
и
«left»
с
использованием
интерфейса
переноса
параметров
рабочего
тела
,
выходной
участок
с
межлопаточным
каналом
связывались
границами
«outlet»
и
«right 2»
с
использованием
того
же
граничного
условия
.
Границы
межлопаточного
канала
,
используемые
при
этом
,
аналогичны
приведенным
на
рис
. 1.
8’2012
9