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

книги из ГПНТБ / Жуковский М.И. Расчет обтекания решеток профилей турбомашин

.pdf
Скачиваний:
17
Добавлен:
30.10.2023
Размер:
15.1 Mб
Скачать

здесь рг = рз/Рь и 3 — угол между направлением кромок и осью решетки.

Величины Х2, 8 и рг связаны при этом равенством

При выводе формул (91) и (92) предполагалось, что в выходном сечении параметры газа достигают критических значений, так что

Рь =Р*-

Формулы (89) или (90) неприменимы, когда волна разрежения достигнет выходного среза решетки. В этом случае в косом срезе имеет место полное расширение, т. е.

М2 sin р2 = 1.

(93)

Такое число М2, которое может быть вызвано расчетным или пре­

дельным, найдется из равенства (89) и (93)

 

= tq(M2np) ’

(94)

Учитывая равенства (6) и (85), получим

Значение угла отклонения потока 8яр, соответствующее полному

расширению, получим, используя равенство

(93),

 

Зяр = arcsin — _ .

--2(1_fe)

----- — Ъ*Ф-

(96)

1/

/ М 1 + *___

 

r k — 1

k t j

k—\

 

Если давление за решеткой р2

ниже давления р2пр, соответствую­

щего полному расширению, то дальнейшее расширение,

сопровож­

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

расхода для следа за решеткой

Ч(М2яр) sin ( 2J0 + Зяр) = q (М2) sin р2.

Будем иметь

„.•„о _

Ч (№гпр) sin (Риф Ъпр)

(97)

sin Р2""

7<м5

'

198

Значения величин q (М2пр) и s’n (Р25ф 4 drip) находим, исполь­

зуя формулы (95) и (96). Далее, задавая давление р2 за решеткой,

определяем q (М2) при помощи таблиц

газодинамических функций

и из формулы (97) вычисляем угол р2.

Значения 8

вычисленные

с использованием формулы (97), показаны на фиг. 90.

косом срезе

Приведенные способы расчета отклонения угла в

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

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

28. ПРИМЕНЕНИЕ ГИДРОГАЗОАНАЛОГИИ (ГГА) ДЛЯ ИССЛЕДОВАНИЯ ТЕЧЕНИЯ ГАЗА

Аналогия между потоком газа и течением тонкого слоя тяжелой несжимаемой жидкости основана на идентичности уравнений, опи­ сывающих эти движения. Случай одномерного движения рассмотрен с учетом сил вязкости Н. Е. Жуковским [63]. Аналогия для пло­ ских потоков газа указана Жуге [64 ] и применена для двухмерных

задач Прайсверком [65], Рябушинским [66] и др. Из основных

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

что аналогия имеет место при условии, когда:

1)газ и жидкость не обладают вязкостью;

2)газ имеет показатель изэнтропы k = 21;

3)вертикальная составляющая ускорения в потоке воды равна

нулю;

4)поверхностное натяжение не проявляется;

5)отсутствуют ударные явления.

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

ния по гидрогазоаналогии [63—71 ]2 посвящены главным образом обоснованию аналогии и рассмотрению способов, расширяющих возможности ее использования для практических целей.

Напишем

уравнения

установившегося 3

движения идеальной

жидкости в открытом канале с горизонтальным дном (фиг. 91)

 

 

W

дш

 

 

dwr

 

 

dwr

1

dp'

 

 

~л— +

w'—^.----

г

дг

 

 

х

дх

 

у

dy

 

?'

dx ’

 

 

'

dw'y

,

,

dw'y

,

,

dw

1

dp'.

 

 

wr

дх

+ W

 

. ..

w

-

 

 

х

 

у ду

1

2

<Z?

P'

dy ’

w

 

dw,

 

dw,

 

dw,

 

 

1

—--- p ™

—r

-I-

W,—x— = —

 

 

 

дх

у

dy

 

x

dz

 

 

 

1 Такой несуществующий газ принято называть гипотетическим.

3 Следует отметить также результаты, полученные Г. Т. Таушкановым В. С. Березниковым и В. Н. Сунцовым.

3 Аналогия справедлива и для неустановившегося движения.

199

Здесь и везде ниже штрихом отмечены величины, относящиеся к потоку воды.

Предположим, что движение происходит с малым изменением

глубин

слоя воды

/?, т.

е.

что

уклон

свободной

поверхности

мал.

 

 

 

 

 

 

 

 

 

В этом случае можно пренебречь ускорением по оси z по сравне­

нию с

ускорением

земного тяготения,

т.

е.

 

 

 

,

dw,

 

dw,

dw'

(98)

 

W

—-Т

h W„—4-----h W—= 0.

 

 

*

dz

у

dy

1

2 dz

 

После интегрирования последнего уравнения системы получаем

 

 

 

Р’—Рат =

P'^(^ — Z).

(99)

Таким образом, при выполнении условия (98) давление в слое воды изменяется линейно. Можно легко показать при этом, что wx и w? бу­

 

дут функциямитолькохи у.

 

Уравнения

неразрыв­

 

ности для дифференциаль­

 

ного элемента

жидкости

 

с основанием dx, dy и вы­

 

сотой h запишутся в виде

 

<?(<A)

 

d(w'yh)

 

 

дх'

+

 

ду'

~

 

Сравнивая

с

уравнением

 

неразрывности

для

пло­

 

ского движения газа

 

Фиг. 91. Изменение глубины слоя в потоке.

<3(^хР)

 

d(wtf) _ „

 

дх

“Г

ду

 

находим, что при наличии геометрического подобия между потоками

несжимаемой жидкости и; газа с показателем изэнтропы k - 2 в соответствующих точках плотность газа пропорциональна глубине

слоя воды

Р

h

(ЮО)

Ро

 

 

 

Рассмотрим уравнение энергии

для несжимаемой жидкости

Р' + -у- ^'2 + ?gz = р'о + ?'gz0.

Используя формулу (99), получим следующее выражение:

■^=2g(h0-h). (Ю1)

200

В случае течения газа уравнение энергии запишется так:

х>2 = 2g(i0 — i).

(102>

Сравнения уравнения (101) и (102), находим, что отношение тепло-

содержаний в газе пропорционально отношению глубин

i

h

г0

^0

или

Т_ _

ПОЗ)

То

h0

 

Из уравнения изэнтропы имеем

 

1

?_ _ /

1

Ро

' \ Ти )

(Ю4)

 

Соотношения (100) и (103) совместимы только в случае, когда k = 2; при этом для давления в газе с показателем изэнтропы k = 2

имеет место

равенство

 

 

_Р_ =

(105)

 

Ро

 

X hd J

Скорость перемещения прямого скачка уплотнения в газе опре-

деляется из

следующего равенства

Для скорости перемещения гидравлического прыжка имеем соответственно

• с'-~ У(107)

Здесь индексы 1 и 2 относятся к величинам до и за скачком уплот­ нения или гидравлическим прыжком.

Для слабых возмущений получаем из соотношений (106) и (107}

(Ю8>

а' = Vgh-

(109)

Таким образом, скорости звука а в газе соответствует скорость-

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

Таким образом, точность получаемых на основе гидрогазоаналогии результатов зависит от выбора допустимых значений глубин.

201

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

«* = /2^7-4 Z<”

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

верхностью, определяемой

по

формуле

 

 

«* =

 

На­

 

значение критической глубины /г* получим из равенства

?*

 

 

 

 

 

_ /

2

\

2

_h*_

~ \ k i

1’ У

3

hu

при k = 2.

 

 

 

 

 

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

 

 

 

 

 

 

 

/г*

=

 

(1 Ю)

Если h < h*, то

течение

воды будет стремительным и

JVV = —~ > 1. Этот случай соответствует сверхзвуковому движе-

V gh

нию газа. При h > h* число М' < 1, и течение будет спокойным.

Такое течение воды соответствует дозвуковому потоку гипотетиче­

ского газа.

ДЕ в гидравлическом прыжке

Потеря • механической энергии

равна

 

-л1 z —

сСУ о

--^-=-ДЕ. (111)

Для случая скачка уплотнения имеем

(112)

Несоответствие между уравнениями (111) и (112) объясняется тем, что переходящая в потоке воды в тепло механическая энергия не вхо­ дит, в отличие от случая течения газа (уравнение (112), в уравнение энергии.

Следовательно, аналогия нарушается как в случае процессов

движения воды с большими потерями, вследствие трения, так и при волновых потерях энергии.

Так, для прямого гидравлического прыжка имеем следующее уравнение связи между скоростями до и после прыжка [65].

3-/.f

(ИЗ)

Л2 — 8

202

В случае прямого скачка в газе имеем:

М2 — 1>

(ИЗ')

 

здесь

 

При числах М < 1,6 ч 1,8 указанные отклонения между скач­

ком уплотнения в газе и гидравлическим прыжком в воде не очень

велики, и соотношение (ИЗ) мало отличается от (ИЗ'). В этом случае изменение величин h и w' происходит на достаточно коротком участке, и гидравлический прыжок (прямой и косой) может рассматриваться как поверхность разрыва ве­

личин

h

и w.

скачка

 

 

/

В

случае косого

 

 

/

(фиг. 92) могут быть с неко-

 

 

/

торым приближением исполь----------------------------------„

 

зованы

при М < 1,6 -г- 1,8

w>

 

 

известные формулы газодина-

_______________

_______

ми ки.

 

 

 

Формула (6) связи между 1

/

 

X и М'

при k = 2 имеет вид

 

/

 

 

ЗМ'2

,, ,

Фиг. 92.

Схема косого

скачка.

 

А

=2 + М-'2‘

<114>

 

 

 

Теоретический анализ влияния вязкости, выполненный Бекке­ ром [68], показал, что ее влияние при глубинах /г>10: 15 мм

вцелом невелико.

Вцелях уменьшения действия вязкости дно канала экспери­ ментальной установки ГГА выполняется из полированного стекла.

Рекомендуется также придавать дну подводящей части канала небольшой уклон. Исследование поля скоростей в слое при наличии уклона дна выполнено В. С. Березниковым.

Вопрос о выборе допустимых глубин hdon тесно связан с обеспе­ чением условия о равенстве нулю вертикального ускорения, Для выполнения этого условия необходимо выбирать размеры модели,

аследовательно, и экспериментальной установки достаточно боль­ шими, так как в этом случае изменение глубин h будет происходить более плавно.

Величина hdon зависит от длины исследуемого тела, с увеличе­ нием которого она растет. В случае очень малых h, когда длины

волн i h, скорость распространения возмущений не зависит от длины волны, т. е. выполняется формула а = ]/gh. В случае боль­

ших размеров модели это соотношение может выполняться

и при

больших глубинах.

К возможно большим значениям hdon

следует

стремиться также и

потому, что при этом уменьшаются искажения,

связанные с силами

поверхностного натяжения, которые

следует

уже учитывать при

h 5 мм.

 

203

Подробное исследование этого вопроса выполнено Лейтоном [67 ], Д. Рябушинским [66], Иппеном и Харлеманом [69], В. Н. Сун­ цовым и др.

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

Н-1

/г— 1

(115)

k + 1

Фиг. 93. Зависимость F/F* от М при &= 1,4 и /г = 2.

Если через это же сопло протекает газ с показателем изэнтропы k’,

то

можно

написать очевидное равенство

 

 

 

 

 

fe+i

 

 

 

 

М k k-\- 1 Г k+ 1

)

 

 

 

 

 

k'+i

 

 

 

 

= М-A...+

М'2\2

 

(И5')

 

 

М' \k' 4- A k' + 1

)

 

Из

соотношения (115') находится связь М = f (М')

при

задании

чисел k и

k'

 

 

 

 

Таким образом, если известны значения М' из опыта на уста­

новке ГТА,

то при задании значения

k, отличного от

2,

из (115')

определяем

М.

 

 

 

на

Зависимость (115) при k = 2 и k,-= 1,4 изображена графически

фиг. 93.

 

 

 

 

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

204

Воспользуемся для этой цели соотношением, определяющим пово­

рот сверхзвукового потока при обтекании угла в зависимости от числа М

о -./£4-1 .

1 /(Z? — 1) (М2 — 1) .

1 it

= ]/ F=Tarcs,n

V 2 + (£ —1)М2 + arCsin

"м~ ~ ~2~ ' (llb>

Полагая в уравнении (116) k = 2, получим аналогичное соотно­

шение для гипотетического

газа

 

fl = р'З arcsin

+ arcsin -------J- .

(117)

В соответственных точках контура, обтекаемого гипотетическим газом (k = 2) и потоком газа с реальным значением k, углы пово­ рота потока р одинаковы, т. е.

р(М', 2) = Р(М, k).

Обтекание

криволинейной

1,8

 

1,7

 

стенки может быть рассмотрено,

 

 

как предельное значение обте-

 

 

кания контура,

составленного

 

 

из бесконечно малых отрезков 7,

 

прямых.

Поэтому

применение

ф

 

формулы

(118)

с

некоторым

1,31,41,51,61,71$ 1,92fl2J 22ZJ12,4Z/iM

приближением

позволяет опре-

V

делять при обтекании криво­

 

Фиг. 94. Зависимость М от М'.

линейного

контура

значения

 

числаМ для газа с показателем

 

М' для гипотетического газа

изэнтропы

k, если

известно число

с показателем изэнтропы k = 2. Зависимость М = F (W) при k = 1,4 приведена для случая Мх — 1 на фиг. 94. Численные значения М и М' даны в табл. 25

Таблица 25

II £ го

1

1,096

1,160

• 1,214 1,268 1,316 1,366 1,416

М

k = 1,4

1,081

1,133

1,176

1,218

1,258

1,294

1,330

II 3?

1,465

1 1,515

1,565

1 1,617

1,6691,720

м

М'

k = 1,4

k = 2

1,365

1,772

1,400

1,826

1,435

1,889

1,469

1,942

1,503

2,002

1,537

2,065

 

2,130

*a-

1,571

1,605

1,638

1,672

1,706

1,741

1,775

М'

м

( k = 2

£=1.4

2,190

1,810

2,260

1,845

2,332

1,880

2,405

1,931

2,488

1,950

2,560

1,986

j

Если скорость набегающего потока

> Уgh

(Мм >1), то

угол поворота потока р — Pj при расширении будет равен

Р -pi = (9- %) + (а - а1),

(119)

205

где

Pj = Oj <- aj---- 1---- условный

угол поворота при расшире­

нии

от Moo = 1 до

заданного значения Мю, вычисляемый по фор­

муле (116), и

 

 

 

 

 

 

6 = 1/HZ агСс.п

-|Z(fe—1)(М^—1)

 

I'

k— 1

rcsl

V

2 + (k— 1)М2 ’

 

 

 

я

 

.

1

 

 

 

arcsin

-,-г .

 

 

 

 

 

 

М

Приведенные соотношения для пересчета чисел М, получаемых

на основе

гидрогазоаналогии, на показатель изэнтропы k

1,4

применены

Г. П. Таушкановым.

 

С некоторым приближением в экспериментальной установке ГГА могут производиться измерения углов выхода потока из решетки. Возможность определения углов выхода потока основана на том,

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

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

tw\h-i sin jij tw9h2 sin p o. (120)

Разделив обе части этого равенства на произведение критических значений величин и”) и /г,*, найдем

'^1 _ 3 I

и>1 А*

2

w9 hiit

Следовательно, уравнение (120) может быть записано таким образом:

 

/

1]

\

з

 

</(A')sinp

q' (k2)sin fi2,

(121)

=

'*01

/

*

'

 

 

где обозначено

з

 

■'

; '2 \

 

 

 

<122>

q’ (лэ=4 а'

(1-44-

206

Если пренебречь потерями энергии, то равенство (121) примет вид

 

д' (aJ)

sin 0, = д' (X’) sin 0,.

(123)

В

случае течения газа

имеем

 

 

д (XJ sin

= д (Х2) sin 02.

 

Из

последних двух

равенств найдем

 

 

.

о

 

<?'(Х2) <7 Рч)

(124)

 

 

 

sin rf, =--------------Sin fi2.

q' (>-i) q (>••>)

Фиг. 95. Рабочий

участок экспериментальной установки ПЛ.

Функции 0', д'(/.'),

а также / , и М'

в зависимости от X' даны

в табл. 26.

'0

X.) и 02 являются извест­

Из опыта в установке ГГА величины X',

ными. Воспользовавшись графиками фиг. 93 или фиг. 94 и формулой (122), находим значения Xj и Х2, после чего по формуле (124) вычис­ ляем угол 02 при показателе изэнтропы k = 1,4.

Точность метода ГГА существенно повышается при соблюдении приведенных выше условий, а также при осуществлении пересчета данных, получаемых в водяном канал.е. Несмотря на приближенный

характер аналогии, ее использование является весьма полезным,

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

вом потоке.

Установка для исследования обтекания решеток профилей пока­

зана на фиг. 95.

207