Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Жаров Г.Г. Судовые высокотемпературные газотурбинные установки

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
17.85 Mб
Скачать

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

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

Обтекание кольцевых и плоских решеток практически одинаково, поэтому при исследованиях часто испытывают плоские решетки, а полученные результаты распространяют и на кольцевые. Так как работа рабочих лопаток отличается от работы сопловых лопаток их движением в потоке газа, то, естественно, различается и их теплообмен. По данным работы [20], согласно опытам коэффициент теплоотдачи для рабочих лопаток в реальных условиях турбины на 20—30% выше, чем для лопаток неподвижной решетки. Однако эти опытные данные не единственные. Имеются и другие результаты, которые весьма противоречивы, что можно объяснить особенностью постановки опытов и проведением их в различных условиях.

Теплообмен от газа к межлопаточным участкам ротора и статора представляет собой сложный процесс. Особенно это относится к начальному участку криволинейного канала, где скорости газа ме­ няются по его длине. Правда, поскольку площадь межлопаточных каналов в сравнении с общей площадью проточной части составляет незначительную величину (5—8%), то принятые значения среднего коэффициента теплоотдачи для лопаток мало скажутся на вычисле­ нии теплоотвода по всей проточной части. Для точного вычисления коэффициента теплоотдачи в межлопаточных каналах можно ввести поправки относительно среднего коэффициента теплоотдачи, которые учитывали бы криволинейность канала и влияние начального участка.

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

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

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

Таким образом, для различных узлов характерны тот или другой вид теплообмена и его интенсивность. В зависимости от задач иссле­ дования определяют локальные (местные) или средние значения коэф­ фициентов теплоотдачи.

§ 30. Локальные значения

коэффициентов

теплоотдачи от газа

к охлаждаемым

лопаткам

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

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

Наиболее эффективным, на наш взгляд, теоретическим методом расчета теплового пограничного слоя по профилю турбинной лопатки является метод, изложенный в работе [22]. Характерная особенность этого метода состоит в том, что расчет теплового и динамического по­ граничных слоев при всех трех режимах течения (ламинарном,

9*

131

переходном и турбулентном) единообразен. Это позволяет сократить расчеты при определении локальных значений коэффициентов тепло­ отдачи.

В работе [22] получена формула для определения интенсивно­ сти теплообмена на профиле в зависимости от параметров погранич­ ного слоя:

 

Nu = PrRe4--

б "

Щ

Г а "

'

 

 

 

 

 

 

rcosp\,

w2co

7*

 

 

где і

— шаг решетки;

 

 

 

 

 

 

 

Р

— периметр

профиля

лопатки;

 

 

 

 

б** —• толщина

потери энтальпии;

 

 

 

 

 

Р 2

угол между направлением

скорости

и нормалью

к оси

w3

решетки

в сечении,

 

где сливаются пограничные слои;

— скорость

на задней

 

кромке

профиля;

 

 

 

ш 2 т

— скорость

потока в сечении,

находящемся

далеко

за ре­

 

шеткой;

 

 

 

 

 

 

 

 

Т2 со температура потока

 

в сечении,

находящемся далеко за

 

решеткой;

 

 

 

 

 

 

 

Го — разность

температур

 

потока

(за

границей

пограничного

 

слоя) и

стенки.

 

 

 

 

 

 

 

Толщину потери

энтальпии определяют

полуэмпирическим

мето­

дом [21 ], одинаковым для всех режимов течения. Этот метод можно использовать, только зная точки перехода, которые находят прибли­ женно. Так как для определения толщины потери энтальпии и точек перехода требуется большая трудоемкая работа, то разработанная программа для ЭЦВМ М-220 позволяет сократить время расчета. Расчет одного распределения коэффициентов теплоотдачи по профилю занимает 15—20 мин машинного времени. Полуэмпирический метод расчета обеспечивает уже при проектировании машин наиболее рацио­ нальный выбор аэродинамического профиля охлаждаемых лопаток и системы их охлаждения. К тому же данный метод можно применить при расчете вращающихся решеток. Однако препятствием для его использования в этом случае служит отсутствие методов расчета распределения скоростей по профилю во вращающихся решетках.

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

Указанным методом можно пользоваться для приближенных и быстрых расчетов.

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

Достоверность данного метода подтверждается хорошей сходи­ мостью его с результатами эксперимента. Максимальное расхожде­ ние находится в пределах 10—15%. Следует отметить, что такой ре­ зультат можно получить, если точки перехода при расчете брать из

эксперимента.

Сравнение экспериментальных данных, полученных

Ж. Вильсоном

и И. Поупом [108], с расчетом по методике [22] для

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

Таким образом,

разработанная в [22] методика позволяет теоре­

тически определять

теплоотдачу по профилю охлаждаемой

лопатки

с удовлетворительной для практики точностью.

 

Помимо теоретического определения локальных значений

коэффи­

циентов теплоотдачи по контуру охлаждаемой лопатки, были выпол­ нены экспериментальные исследования. Наиболее фундаменталь­ ные из них были проведены М. Н. Бодуновым [3] на специальной экспериментальной установке. Локальный коэффициент теплоот­

дачи определялся

по формуле

 

 

 

 

 

 

а,. -

0 , 8 6 4 Г К

 

 

 

 

 

 

 

 

 

где Ї

— сила электрического тока, проходящего через

пластинку;

 

V

— напряжение

на концах

измерительного участка;

 

а —• расстояние

между местами замера напряжения;

 

В — ширина

пластинки;

 

 

 

°\,

62

— ширина

промежутков с обеих сторон пластинки;

 

Тот — средняя

температура

пластинки;

 

Для

Тг

•— температура

торможения вне пограничного слоя лопатки.

канедой полоски находили локальное значение числа

Нуссельта

 

 

 

 

Nu = axb

 

 

где

b — хорда профиля;

 

 

 

 

А, — коэффициент

теплопроводности

воздуха.

 

На

каждом угле атаки подсчитывали

величины

 

где

w, р, yi определяли

по параметрам

воздуха

за решеткой;

 

 

" _ Рпр — р а

і

 

 

 

 

и —

*

 

 

где

р п р — статическое

давление

на

профиле;

р 2

— статическое

давление за решеткой;

р*—давление

торможения

за решеткой.

Результаты опытов были представлены в виде зависимостей

Nu = f(x)

при изменении Re в пределах 2,16 -105—6,34 • 10Б и различных углах атаки.

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

Nu = cRe".

Для входной кромки за определяющую температуру принимали тем­ пературу торможения воздуха перед решеткой, за определяющий размер — удвоенный радиус входной кромки. Число Re подсчиты­ вали по теоретической скорости за решеткой. Для сравнения тепло­ отдачи на кромках лопатки приведены результаты исследования теплоотдачи по профилю цилиндра.

Для входной кромки профиля лопаток опытные точки хорошо укладываются на кривую

 

 

Nu = 0,265Re0 '58

(88)

при і «=* 0, Re =

5-Ю3 —5-Ю4 .

 

Для

выходной

кромки

 

 

 

Nu = 0,003Re0-93

(89)

при і =

10---10,5о; Re = 6 • 103 н-3 • 10*.

 

Таким образом,

разработанные методики расчета

и эксперимен­

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

§ З І . Средние значения коэффициентов теплоотдачи от газа

кохлаждаемым лопаткам

Вбольшинстве задач, связанных с передачей

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

делении количества

тепла, отводимого в процессе охлаждения.

Если при расчете

температуры стенки лопатки с использованием

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

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

Рассмотрим одно из первых экспериментальных исследований по определению среднего коэффициента теплоотдачи охлаждаемой ло­ патки (1940 г.). Рабочим телом при проведении опытов был слабо перегретый пар. Лопатки были выбраны активного типа и охлажда­ лись водой. Пакет состоял из трех лопаток симметричного типа. При исследованиях угол натекания не менялся. Все измерения проводили на среднем значении высоты лопатки. По наружному контуру ло­ патки заделывали 14 термопар, а по охлаждаемому — четыре тер­ мопары. Такое количество термопар позволяло определить темпера­ турное поле лопатки и удельные потоки тепла и даже локальные зна­ чения коэффициентов теплоотдачи. С целью определения числа Рейнольдса перед пакетом лопаток и после него замеряли расход пара и его параметры. Опытные данные были обработаны в виде зависимости среднего по сечению числа Нуссельта от среднего числа Рейнольдса. Число Рейнольдса измеряли в пределах (1,2—6,0) * 105. Коэффициент теплопроводности пара определялся по температуре стенки, а коэф­ фициент вязкости — по средней по объему температуре потока. За определяющий размер был принят гидравлический диаметр.

При этих условиях было рекомендовано соотношение Nu = 1000 + 0,006Ще.

Однако полученные результаты оказались завышенными. По-види­ мому, это можно объяснить либо влиянием стенок канала, либо конденсацией пара на поверхности лопаток.

Наиболее значительные эксперименты, проведенные по опреде­ лению среднего коэффициента теплоотдачи, изложены в работе [20]. Опыты выполняли на плоском пакете, состоящем из трех лопаток и обтекаемом воздухом. Расход воздуха определялся при помощи спе­ циальной мерной диафрагмы.

Тепловой поток имел направление, обратное его протеканию при работе турбины. Количество тепла, отведенного от лопатки, подсчи­ тывали по затрате мощности на ее нагревание. По профилю лопатки заделывала 10 термопар. Число Рейнольдса изменялось в пределах (0,7—3,0)-105. При обработке опытных данных определяющим раз­ мером считали эквивалентный диаметр. В критериальное уравнение вводили отношение средних температур стенки и воздуха, а также относительный шаг. Во время опытов находили среднее по периметру профиля значение коэффициента теплоотдачи, которое относилось к разности температур торможения газа и лопатки. При этом темпе­ ратуру лопатки определяли как среднеарифметическое значение из измеренных значений температур в десяти точках профиля. Все физические константы вычисляли по средней температуре в канале, а скорость потока — по среднеарифметическому значению между скоростью потока на выходе и входе в канал. Опыты проводились при использовании лопаток с различными профилями: с большой реактивностью и различной толщиной входной кромки.

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

Nu = , R e « ( £ ) « ( l ) ° > " \

, 9 0 )

где с — коэффициент, учитывающий влияние формы профиля ло­

Тг

патки (с = 0,135н-0,160);

 

— средняя

температура

газа;

 

Тст

— средняя

температура

поверхности лопатки;

 

 

+ w9

 

 

 

wcp =

1 '—-

средняя относительная скорость газа.

 

Из работ зарубежных ученых в этой области следует привести

работы К. Баммерта и X. Элленброка.

 

В

работе

[2] критериальное

уравнение для реактивной

лопатки

получено в

виде

 

 

 

 

 

 

Nu = 0,72Re0 '5 4 5 .

(91)

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

В работе [111] для подсчета среднего коэффициента теплоотдачи -предлагается зависимость

Nu = cRe-'Pr1'3,

Таблица 19

Критериальные зависимости, полученные при определении среднего коэффициента теплоотдачи от газа к лопатке

А в т ор

Ф о р м у л а

R e

[ и с т о ч н и к ]

 

 

Г. С. Жирицкий,

 

Nu

=

 

wdep

=

0,130Re°'6 6 e

0 '5

В. И. Локай

[20]

 

 

Решетка 1

 

 

 

 

 

 

 

Г. С. Жирицкий,

 

Nu

=

 

wdep

=

0,154

Re°-6 S 0

0 '5

В. И. Локай

[20]

Iі

 

Решетка 2

 

 

 

 

 

М. Н. Бодунов

Nu = 0,271Re0 , 6 3 1

wxbp

[3]

 

 

Решетка 1

 

V-

И н т е р в а л

О п р е д е л я ю ­

 

 

і lb

 

 

Pi/Ps-

щ и е

 

M

lib

b

и з м е н е н и я R e

 

 

т е м п е р а т у р ы

 

 

 

 

 

 

7-Ю4 —3-Ю6

V- (<c),

1,21

0,5

0,715

— -

100

56/34

 

P Vc, P)

 

 

 

 

 

 

7-Ю1 —3-Ю5

1-І Cc).

1,21

0,5

0,715

54/34

P Vc, P),

 

X ( * c )

 

 

 

 

 

 

1,5-105 —4-Ю5

v(0-

1,17

0,5

0,721

1,26

95,4

60/28

P{t, P),

М. Н. Бодунов

Nu =

0,0646Re0 '7 3

Wjbp

1,5- 10s—5- 10Б

p(t, p),

1,17

0,5

0,576

1,5

79,5

40/30

[3]

Решетка 2

V-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А. Смит [106]

 

 

Wibp

4-Ю'1 —2-Ю6

1-1(4).

~ 1

0,75

35/54

Nu =

0,0205Re0 '8 8

Pit. p)>

 

 

 

¥•

 

 

 

 

 

 

 

 

3

со

А п т ор [ и с т о ч н и к ]

С. А. Эндрюс,

Р.С. Бредли

Е.Польман

Д. Эйнли

Ж- Вильсон, И. Поуп

К- Баммерт [2]

ХТГЗ

Ф о р м у л а

Nu = 0,169Re0 '6 6

Nu = 0,045Re0 '7 7

Nu = 0,100Re0-7

Nu = 0,150Re0 '6 5

Nu = 0,72Re0 -5 4 5

Nu = 0,394Re0 '5 1

Re

w2bp

wdep

ал, ftp fi

w2bp

P-

P-

wdep

P-

Продолжение табл. 19

Интервал

О п р е д е л я ю ­

rB /r*

 

 

 

 

и з м е н е н и я R e

щ и е

M

t/b

lib

b

 

т е м п е р а т у р ы

 

 

 

 

 

4- 10і—4-106

8-10"—2-Ю5

8 • 10*—2 • 105

1,5-106 —7-Ю6

106—2.8-106

2-Ю4 —6-Ю4

P- ((a).

0,455

0,2—

0,62

2,0

25,4

45/45

Pit,

P),

\(tB)

 

1,0

 

 

 

 

 

 

 

 

 

 

V

(tc),

1,21

0,1

0,715

31/31

P Uc, P),\

K(tc)

 

 

 

 

 

 

Hit*),

0,82

0,6

0,687

2,4

26,6

45/33

Pit,

P),

Ч'в)

 

 

 

 

 

 

Pit,

P),

1,04

0,25

0,625

1,67

152

60/30

Ч ' в )

 

 

 

 

 

 

H

ik),

1,2

0,56

0,741

3,0

65,2

53/30

P ik, Pi),

M ' I )

 

 

 

 

 

 

M ' I )

0,93

0,25

0,805

1,87

36,4

75/28

P ihP),

М'в)

где для активных лопаток при Re = 3-Ю3 —4-Ю5

 

 

п

= 0,65ч-0,74 ^

0,7

и с = 0,103;

для

реактивных

лопаток при Re =

7-Ю4 —2-Ю5

 

 

 

я = 0,57 и с = 0,398.

Критериальные зависимости для определения среднего коэффи­

циента

теплоотдачи

получены и рядом других исследователей. Наи­

более

характерные

и применяемые

из

них при расчетах представ­

лены

в

табл. 19. Широкое распространение из всех приведенных

зависимостей получили формулы (90) и (91).

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

нениям при их использовании. Это

 

 

 

 

 

объясняется

тем, что условия прове- о,п

 

 

 

 

ден и я опытов и способы их обработки

1

 

 

 

 

бывают различные. Следует отметить,

д д

 

 

 

 

что условия

теплообмена, как прави-

'

 

ft'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л О , О п р е д е л я е т Д Л И Н а ПО К О Н Т у р у O X - Of

 

 

 

 

лаждаемой лопатки

ламинарного

и

0,4

 

1

 

 

турбулентного

пограничных

слоев.

 

 

 

Опыты обрабатываются

в зависимо­

0.2

 

\

 

 

сти от чисел

Нуссельта

и от

чисел

 

 

 

Рейнольдса, при этом

не учитывается

 

 

 

 

 

расстояние

до

точки

перехода

из

0,2

0,k

0,6 0,8

1,0

1,2 1,* г У

ламинарного

слоя в

турбулентный.

Рис.

85.

Зависимость

коэффициен­

В наилучшей степени данный пе­

та

с и показателя

степени п от

реход учитывает параметр у =

^

||2

,

параметра у.

 

 

который представляет собой, по существу, отношение ширин

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

Nu = (0'0^805 0,0022) PrvsRe o,74Y M3, до)

Эта формула обобщает 17 критериальных зависимостей вида Nu = cRe"Pr"!.

На рис. 85 представлены зависимости коэффициента с и показа­ теля степени п от параметра у. Максимальные отклонения, вычислен­ ные для крайних значений чисел Re, составляют в девяти случаях менее 10%, в шести случаях 10 и 15% и в двух случаях 15—20%.