Файл: Зысина-Моложен, Л. М. Теплообмен в турбомашинах.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 132

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

5 =

Ч

= Рг, + 5Рг, У Ц - {JL (1 - Рг,) [ £ + -§- <1 - РГ,)] +

 

2 St

 

 

 

+ (

w

-

1) +

, n [ 1 +

i (

^

-

1) ]

По данным работы [237] S

Р х —0,57- по экспериментальным дан-

 

ным S

0,9 ч-1,1.

В ра­

 

 

 

 

 

 

боте [74] для этого пара­

 

 

 

 

 

 

метра

получены значения,

 

 

 

 

 

 

представленные на рис. 42,

 

 

 

 

 

 

из

которого видно,

что

S

 

 

 

 

 

 

зависит от М иф =

TwITaw

 

 

 

 

 

 

и значительно

отличается

 

 

 

 

 

 

от единицы при обтекании

 

1,0

?,о

з,

 

сжимаемым потоком.

 

 

 

 

В работе

[239]

на ос­

 

Рис. 42. Зависимость S от ф

 

 

нове обобщения

результа­

 

 

 

тов

прецизионных

опыт­

 

 

 

 

 

 

ных данных предлагаются следующие формулы для расчета моди­

 

фицированного фактора аналогии Рейнольдса а, представляющего

 

собой

отношение

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

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

обмена

 

тепла

 

к коэффициенту турбулентной вязкости а = К/А:

 

 

 

 

1)

при у+ <

„+

 

 

 

 

 

 

 

 

 

г/?

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а = Р г—0,5

 

 

 

 

 

2)

при y t

<

у+ < y t

 

п(у+

 

 

 

 

 

 

а-

 

Рг°’5+

 

 

' У^т)

 

 

 

 

 

 

1 -h (Рг0-5 — 1) sin

] ‘

 

 

 

 

2 Рг0,5

 

 

 

 

 

 

 

3)

при у,+ > Ш

 

 

 

 

 

 

 

 

а1.

При этом

У? = 4; УТ — Ю; y L yt + yt ytm yt — yt

В работе [228] предлагается уточнение формулы Прандтля для длины пути перемешивания I путем введения демпфирующего множителя Ван-Дриста:

1) при у!б < 0,1

■f = 0,41 [ l - e x p ( - - g - ) ] - g - ;

9

131


Рис. 43. Зависимость 6л/6 от 1|) и М<ю

Рис. 44. Зависимость

т)0 от г|э при

следующих

значениях Мга:

при R = 105

О — 6,8; •

— 5,0; V — 3,0;

X — 2,0; J

— 1.5

 

2)

 

при

0,2 <

у/8 < 0,6

 

-

4 -= 0,41 1 — еХр ( — |^ ) ]

— 1,53506

 

 

 

+ 2,75625 ( - |---- 0,1 ) 3

1,88425 —

3)

при

ylb >

0,6

0,089.

 

 

 

 

//8 =

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

денных на рис. 43 кривых [74] видно,

что при ф =

T J T aW — 2 и

Мсо *=» 5 толщина подслоя

6Л может

составлять

50%

толщины

всего

пограничного слоя

8. Безразмерная

толщина

подслоя

г|о =

yvjv, являющаяся при обтекании пластины

несжимаемым

потоком величиной постоянной и равной ri0 =

10,65 (штриховая

линия на рис. 44), с изменением температурного фактора в слу­ чае сжимаемого потока, как видно из рис. 44 [74], заметно изме­ няется. Линия на рисунке соответствует Мга г» 0.

Все эти данные говорят, что основные допущения, принимае­ мые при построении приближенных методов расчета турбулент­ ного пограничного слоя (Pr, = 1; т)0 = const; S = 2St/cf = 1; l = ху и т. д.), становятся практически неприемлемыми для рас­ чета теплообмена при значении ф, значительно отличающемся от единицы, и при больших значениях числа Мт .

21. Расчет турбулентного пограничного слоя при обтекании с продольным градиентом давления

Задача о градиентном обтекании поверхности потоком сжимае­ мого газа (обтекание крылового профиля и криволинейных поверх­ ностей, течение в межпрофильных каналах решеток турбинных и компрессорных лопаток и т. п.) значительно осложняется тем, что при сверхзвуковых течениях сильные продольные градиенты давления или изменение кривизны обтекаемой поверхности могут сопровождаться возникновением скачков уплотнения и вызывае­ мых ими отрывов пограничного слоя. Неоднородность основного потока, непостоянство давления поперек пограничного слоя и другие факторы, обусловленные взаимодействием скачка уплот­ нения с пограничным слоем, приводят к невозможности примене­ ния уравнений пограничного слоя в обычном виде. Описанные в литературе методы расчета как ламинарного, так и турбулент­ ного пограничного слоя у поверхностей, обтекаемых сжимаемым газом с продольным градиентом давления, основаны на допущении

133


об идеализированном (бесскачковом и безотрывном) обтекании. Это нужно иметь в виду, используя эти методы для решения прак­ тических задач.

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

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

Напишем интегральное соотношение импульсов в виде

т

Дб** + 6** (2 + Я )U^' +

d p '

РV*

d x

d x _

Пусть профиль скорости в пограничном слое соответствует закону одной седьмой. Тогда 8*/8** 1,28. Соответственно допустим, что для т можно использовать выражение [118]

—0,25

т = 0,0128pt/2 ( - ^ —)

Введя эти значения в интегральное соотношение импульсов после преобразований получим

d R * *

d x

у

_

0,0128_______ у_

2,29UR**

 

2,29 (R**)1,25

(IV.98)

d x '

Если принять, что обратный переход происходит при тех же значениях R**, что и обычный переход от ламинарного к турбу­

лентному течению, т.

е. при R**

360,

то

 

d R**

у

 

 

v d U

(IV.99)

d x

825U 3,55-10~6—

Тр Чх '

Поскольку при обратном переходе

 

величина

dR**

 

d должна

уменьшаться, т. е. dR** < 0, то условием обратного перехода

будет

-jjs- > 3,5 • 10~®.

Существующие в настоящее время полуэмпирические методы

расчета теплового и динамического пограничных слоев при гра­

134


диентном обтекании можно разбить на две группы. К первой группе относятся методы, базирующиеся на решении дифференциальных уравнений в частных производных [уравнения (VII)—(X) в гл. I] с привлечением эмпирических соотношений для зависимостей так называемых эффективных коэффициентов переноса (вязкость, теплопроводность, число Прандтля и т. д.). Ко второй группе относятся методы, часто называемые эмпирическими, основанные на решении обыкновенных дифференциальных уравнений (обычно интегральные соотношения импульсов и энергии).

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

Из работ первой группы в этом плане особого внимания заслу­ живает метод Патанкора и Сполдинга [226]. Эта работа является дальнейшим развитием идей метода Сполдинга и Чи [241 ], описан­ ного в предыдущем параграфе. Широкое применение ЭВМ позво­ лило авторам произвести подробный численный эксперимент и сопоставить результаты расчетов с экспериментальными данными с целью отбора гипотез, определяющих коэффициенты турбулент­ ного обмена, использование которых дает наилучшее совпадение расчетных результатов с непосредственными измерениями для конкретных задач. На основе этих сопоставлений проведено дальнейшее развитие метода и усовершенствование основных физических гипотез, а также усовершенствование математического аппарата с точки зрения повышения его точности и экономичности. Метод численного решения основной системы дифференциальных уравнений пограничного слоя является конечно-разностным, он подразумевает использование ломаных профилей из прямолиней­ ных отрезков. Независимой переменной в поперечном сечении пограничного слоя служит безразмерная функция тока. Непре­ рывное изменение эффективной вязкости поперек погранич­ ного слоя характеризуется модифицированной гипотезой ВанДриста [252]:

\Ьи = И + Р*У 1 — ехр - У

ди I

(IV. 100)

д у г

 

 

где величина

ди

P*V ду

называется турбулентным аналогом вязкости.

Метод решения задачи запрограммирован на ЭВМ, программа записана на алгоритмическом языке «Фортран-IV», подробное ее

изложение приведено в книге

[226], недавно вышедшей на рус­

ском языке.

. . . . . . .

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


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

Так, известный метод Солодкина для расчета обтекания сжи­ маемым потоком представляет собой обобщение известного одно­ параметрического метода расчета турбулентного пограничного слоя в несжимаемом газе Лойцянского [115]. Распространение идей этого метода расчета на тепловой пограничный слой шло двумя путями [60, 150].

В работе [60] разработан приближенный полуэмпирический метод, позволяющий рассчитывать единообразно все три области теплового пограничного слоя (ламинарного, переходного и тур­ булентного режимов течения), возникающие при обтекании по­ верхности сжимаемым газом. Рассматривается обтекание поверх­ ности плоским потоком вязкой сжимаемой жидкости при наличии теплообмена в предположении постоянства температуры поверх­

ности. Интегральное соотношение энергии

 

 

■+ <£

и

То Nil*

(IV.101)

d\

_L

 

и

Тц) Р г R*

 

решается в переменных Дородницына. Здесь переменные | и т)

определяются по (1.53), 6Т — по (1.77), t и to — по (1.39); Nu* = ах/Х\ Pr = v*la*.

Использовав соотношения (1.54), (1.56), (1.61), можно уравне­

ние (IV. 101) привести к

виду

ft

 

 

I и ' я**

£I ш 0

Nu,

 

2\ f t - I

(IV.101a)

+ 7Г6т

 

■оЗ)

Pr Rjc

Вводятся параметры:

сV

~Ц~ От

.

(IV. 102)

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

136