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

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

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

Система

(256) содержит п уравнений с (п -)- 1) неизвестными: п

величин ут

в серединах участков

и значение •фдг, причем точка N

взята

на контуре БЛ И \\>N =

const,

так как контур лопатки яв­

ляется

линией тока.

 

 

 

 

 

Для замыкания системы

(256) используют постулат Чаплы­

гина—Жуковского, согласно которому на задней

кромке А ло­

патки

вихревая плотность уА =

0.

Тогда

(п. -|-

уравнение

имеет

вид

 

 

 

 

 

 

 

 

= Ѣ\о - 1 7

S УтаАм

As,„.

(257)

m m=l

При практических расчетах функцию тока невозмущенного потока т|)0 удобно представить через функцию тока обычно задан­ ного потока "Фх перед решеткой и искомую циркуляцию Г вокруг лопатки. В самом деле, для равномерного потока комплексный потенциал

W0 = w0e~iaz — w0 (cos a — i sin a) (x - j - iy) = cp0 -(- i\\>0,

откуда функция тока равномерного потока я|)0 = wXoy — wUox или, пользуясь треугольником скоростей (рис. 139), согласно которому

- Г л

получим

причем знак плюс соответствует компрессорной решетке, а минус— турбинной.

Циркуляция скорости вокруг лопатки

 

Г я =

J V (s) ds =

2 Ут As„

 

 

л =

 

 

 

 

ш = 1

 

выражается

через те же значения ут

и As,,,, которые входят в си­

стему (256).

 

 

 

 

Рассмотрим поведение коэффициентов системы aNM

в отдель­

ных особых случаях, когда точки N и M или совпадают, или

сближаются. Совпадение точек имеет место п раз, когда

«бегущая»

по контуру

точка M попадает на фиксированную точку N. Сбли­

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

числяются коэффициенты aQA

или аРА.

На участках

AsP

и AsQ

в этом

случае

распределение у (s)

можно

считать

линейным:

 

 

Y (s) = -fïE.s .

 

 

 

На

остальных

участках у,п

можно

считать

постоянной

и рав­

ной ее значению

в середине

участка

Asm .

 

 

 

260

Рассмотрим

коэффициент

амм

в выражении

 

 

VmßMMA s,„ =

lim

ymaNM ds =

 

 

M->N

J

 

 

2

 

 

 

= 2 у т 1 і т

Г ln [ch *N

. е м 2Я c o s ^ . % ! 2n

ds. (258)

M-+N oj

Преобразуем подынтегральное выражение, пользуясь разло­

жением в ряды cos x

и ch x для малых

углов, получим

 

 

 

 

 

\2

ch » t

™ 2я = 1 +

w

, а « 2я

cos

2 я = 1 - 4 - (

^ 7

^ 2 я ) \

Следовательно, можно

записать,

что

 

ch ' Ѵ 7 І Л 1

2n - co s Х " ~ П м 2л =

= i(xN—i/Vj)2_f~(№v — 'Пм)2]

^

— ~ ( 2 ~ s2>

где s — расстояние по прямой между точками M я N. Вычислим правую часть выражения (258).

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

т J i n - ^ - s 2 d s = 2 y m A s m ( l n - ^ - + 1 п я — 1

g-In 2 ) ,

о

 

 

 

 

откуда

 

 

 

 

 

 

aMM = 2(\n^fL-

0,202).

(259)

Для

участка, прилегающего к задней кромке,

запишем

yPapAAsP=

J • ^ - 5 І п 2 ^ - з 2 й ( з = у р А 5 р ( 2 1 п - ^ - - | - 1 п 2 я 2 — l j ,

 

о

р

 

 

откуда

 

 

 

 

 

 

а Р А = 2 1 п - ^ - +

1,98.

(260)

Обычно на практике рассматривается обтекание решетки тон­ ких профилей произвольной формы. В этом случае при произволь-

261

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

Такой случай был впервые рассмотрен проф. В. В. Уваровым в виде дополнения к разобранному выше методу П . В. Мелентьева. Рассматривая участок лопатки, включающий переднюю

кромку (точка

О, рис. 1 4 0 ) , можно заметить,

несмотря на то, что

у0 стремится к

бесконечности, интеграл | у a

ds по участку Asx

As,

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

Следовательно, указанное на рис. 140 распределение у (s) можно заме­ нить прямолинейной зависимостью вида

 

 

 

 

Y(s) =

Yo (YO Y I ) - Ц - S .

 

 

 

 

Вычисляя теперь для

участка Asx

 

 

 

интеграл

вида

 

 

 

 

 

 

 

 

J Уаым ds,

 

 

 

 

 

 

 

As,

 

 

 

 

 

получим

 

 

 

 

Рис.

140.

Распределение

вихре­

As,

 

 

 

 

вой

плотности в области

перед­

 

 

(Yo —Yi)

 

X

ней

кромки при произвольном

 

 

às.

 

угле

натекания потока

 

 

 

 

 

 

 

X In 2 -p-s2 ds = yi_ Asx

2 In

+

1,98 —Yo

ASJL-

Следовательно, в случае «ударного» входа на тонкую лопатку для участка Asx получим

y.La01 àSi = Yi Asx ( 2 l n - ^ - + 1,98) — y0 Asb

Причем величина y0 входит только в первое уравнение си­ стемы (256). Для замыкания системы (256) составляют дополни­ тельное уравнение для точки О (передняя кромка).

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

Циркуляция вокруг лопатки

262

Угол выхода ß 2

потока определяют по одной из рассмотрен­

ных выше формул,

например,

 

 

 

t g ß 2 = tgßi

Гл

 

 

wr

t '

 

 

где знак плюс для турбины, а минус для компрессора. Распределение скоростей по контуру лопатки в решетке,

а следовательно, и давлений, находят по известным соотношениям

 

 

 

до..

ду

 

11

тУ = -

дх

 

где ір =

wXly

— [wai

± -^)

 

Л

— ізг

J V l n

ch 2л •

— cos 2л. y

1 1

 

 

 

 

 

 

 

 

 

Выполняя

дифференцирование,

получим

 

 

 

WxN

= W x ,

'

 

 

 

/Ѵ/ѴІ

 

 

 

wyN

 

 

 

 

 

 

 

 

где

 

 

 

 

.

,

 

 

 

 

 

 

 

 

 

 

 

°NM — •

- r

-

sin

(4f,v т ) Л І )

 

 

 

XN

 

ем

• cos 2л уы — Ли

 

 

 

ch 2л

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

,

,

 

 

 

 

 

 

s h —

(х.,ѵ-

IM)

 

 

 

слгм

ch 2я

XN — i .

• cos

!/Л' Члі

 

 

 

 

t

 

Л1

 

 

 

 

 

 

 

 

 

t

Вихревая

плотность y

(s)

была

найдена

выше.

§ 30. МЕТОД ГОДОГРАФА СКОРОСТИ ДЛЯ РЕШЕНИЯ ЗАДАЧ ПОТЕНЦИАЛЬНОГО

ТЕЧЕНИЯ ЖИДКОСТИ С ОБРАЗОВАНИЕМ СТРУЙ

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

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

Во всех этих задачах форма струй заранее неизвестна, ее нужно определить в ходе решения. Зная форму струи, можно найти коэффициент сужения и, следовательно, расход жидкости через щель.

263

На поверхности струи, вытекающей в пространство с постоян­ ным давлением, давление также постоянно.

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

"Г + Т = c o n s t '

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

Для плоского течения несжимаемой жидкости существует ком-

IV/ / ч dw плексныи потенциал W (z), производная которого —^- дает так

называемую сопряженную

комплексную

скорость

[см.

выраже­

ние (203)].

число х—іѵу)

 

 

 

Обозначим комплексное

через

£,

тогда

ѵх— іѵу

= l = I + ni,

 

(261)

а плоскость комплексного переменного £ назовем плоскостью

годографа

сопряженной скорости.

В самом деле, каждому значению скорости в произвольной

точке M

в плоскости потока соответствует точка М' плоскости

годографа (рис. 141), причем в силу сопряженности скорости

получается зеркальное отображение относительно оси

так как

I = ѵх, т| = — ѵу.

 

Скоростям точек поверхности струи от А (В) до С, очевидно, будет соответствовать дуга окружности А'С В' радиуса ѵх.

Ранее было установлено, что проекции скорости ѵх и ѵу внутри жидкости не могут иметь экстремального значения. Максимум достигается только на границах. Следовательно, все линии тока

264

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

границы полукруга

А'С'В',

подобно линии

О'М'С.

Следовательно, задача отыскания линий тока в

плоскости тече­

ния z решается для плоскости годографа сопряженной

скорости £.

Определим каким-либо способом линии

тока

в

плоскости £

и найдем комплексный потенциал этого течения W (£).

Тогда, пользуясь выражениями (203) и

(261)

и,

помня, что

W (z) = W (Q -f- С,

получим

dW (z)ldz =

£.

Следовательно,

можно найти форму линий

тока

в плоскости потока z, т. е.

 

г=\™&-.

 

 

 

(262)

Итак, задача определения формы струи сведена к решению интеграла, содержащего комплексный потенциал в плоскости £.

Переіідем к решению задач.

Задача 1. Определить форму струи и коэффициент сужения струи при исте­ чении из щели в плоской стенке (рис. 141).

Перенесем линии тока из плоскости потока г в плоскость £ годографа. Если сосуд с жидкостью имеет неограниченную емкость, то можно считать, что доста­ точно далеко от щели жидкость находится в покое. Тогда это состояние жидкости в плоскости годографа изобразится точкой О'. Линии тока вблизи щели перейдут в кривые, выходящие из точки 0' и входящие в точку С , лежащую на окружности радиуса vi. В плоскости потока точке С соответствует то сечение струи, в котором скорости выровнялись и равны Vi.

Следовательно, в плоскости £, надо осуществить поток, у которого линиями тока служат указанные кривые, диаметр А'В' и дуга А'С'В'.

Очевидно, для-этого необходимо поместить в точку О' источник обильности 2Q (пока неизвестной), а в точку С и ей симметричную С" — стоки с обильностью 2Q каждый. В результате получим внутреннее обтекание круга радиусом vi с обте­ каемым диаметром А'В'.

Комплексный потенциал такого течения

^ ( 0 = - § - ш С — ü - l n ( S - ° i ) — і - 1 п ( 5 + °і)

или

W = [In £ - ln-(S - ѵЛ - In (£ + 0 l ) ] •

Определим дифференциал комплексного потенциала, входящего под знак интеграла в формуле (262):

dW(Ç) = ^л

-±-}dl.

Я .

 

Картину линий тока в плоскости потока найдем при решении интеграла

или

z=*

I'

f

с dt,

 

 

 

 

Второй и третий интегралы находят разложением на рациональные дроби: 1 _ А В

— ~Г'Т,

265

или Л£ — A vi - j - Bt, = 1, откуда Л - f S = О и Л = — 1/і>і, тогда ß = 1/оі, а интеграл

С ( S - о і ) = - - M n £ + - M n ( £ - t ; 1 ) .

Аналогично получим для третьего интеграла

S(S + »i)

»,

Следовательно,

 

- { + A l n c _ J _ l n ( s _ , 1 ) _

J - I n g + - A in (5 + ^ ) 1 + С ,

 

 

_ i

+

_ L i n l ± 3

 

(263)

Для определения формы струи выделим в полученном выражении действи­

тельную и мнимую части.

 

 

 

 

 

 

Рассмотрим

плоскость

годографа.

 

 

 

 

Пусть £ определяет скорость в какой-то точке на границе струи.

Найдем логарифмы в выражении (263):

 

 

 

In

-J- Di) = In ri -f-

tor,

In (Ç oi)

= In гг +

і«2,

где

 

 

 

 

 

 

 

' i = У (S + * i ) 2 + л2 ; r2 = y (g - ^ ) 2

+ л2 ;

 

a ^ a r c t g ^ j L - ;

a ^ a r c t g - ^ . .

 

Тогда

 

 

 

 

 

 

 

Х + *У = 1Г — Iа , " I + І П

+ — — « о )

2-

 

I 2 - f i f

1 t>!

 

/"a

u i

-

/ j i i i

 

 

 

 

 

 

 

(264)

Выделяя коэффициенты при мнимой единице і,

получим

 

 

У = "I [рф|г +

К

- «*)] + С2 .

 

Определим

коэффициент

сужения

струи

 

 

 

Поскольку величина Сг в выражении для у

неизвестна, то исключим ее, составив

разность

 

 

 

 

 

 

 

Q Г

il

^

i

1

i

чі

УА - Ус = т

[ ¥

Т

+ —

( * г -

<н)\А. -

0_ і ф р + ^ К - ^ ) ] , , .

266

или, подставляя координаты точек А' и С и значения углов a i и аг, получим

Ул~Ус = 1 Г

4}

1

4

« ) - - £ - ( - * ) .

 

 

 

 

ИЛИ УА—УС

=

Q/яоі .

 

 

2Q/nvi,

 

 

i'id, имеем hd = 2d/n.

Откуда находим Л—d =

или, заменяя Q =

Следовательно,

ц = d/Л =

я / ( я -(- 2) = 0,611.

 

 

Это значение коэффициента сужения довольно хорошо совпадает с определен­

ным из опыта.

 

 

 

 

 

 

Расход

жидкости через

 

щель

подсчитывается

по

формуле

Q = v1d = v 1 h 1 ^ ,

где і>і находят по уравнению Бернулли. Форму струи определяют по выра­ жению (263), где g H г) берут с окружности радиуса t>i.

Рис. 142. Истечение через плоский насадок Борда

Задача 2. Определить коэффициент сужения струи при истечении из плоского насадка Борда с отрывом от стенок. Поступая по указанной методике, перенесем линии тока из плоскости z потока во вспомогательную плоскость £ годографа.

Очевидно, картина линий тока будет соответствовать указанной на рис. 142. Такое течение можно осуществить, помещая источник обильности Q в начало координат О' и двойной сток в точку С на окружности радиуса vi в плоскости

годографа £.

Комплексный потенциал этого течения имеет вид

W(0

[In £ - 2 1 ^ - ^ ) ] .

Найдем дифференциал

этого

комплексного потенциала:

 

 

І _ 2 — L

dt,.

Картина линий тока определится при решении интеграла:

Используя решение предыдущей задачи, получим:

Q

1 I

2

,п

+

ІП £ - » 1

267

Рассмотрим плоскость годографа. Для точки на поверхности свободной струп можно написать

 

 

Іп £ = ln/i - f - 'di;

In (£ — i^) =

In r 2

+

('a2 .

 

 

Выделим действительную и мнимую части в выражении

для г:

 

 

 

 

Q Г

 

£ — и і

,

2

 

г,

,

2 .

 

i - a 2 )

С.

2 = ^ + 1 « / = - ^ [ - ^ - ^

 

+ - 1 - 1 , ^ + — /

(

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(265)

Откуда

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У =

U

aт т а + 7 г ( а і а

г ) 1

 

- Н е ­

 

 

 

 

 

 

+

if

 

 

 

 

 

 

 

 

 

 

 

 

определим

коэффициент сужения

струи,

для

чего

найдем

разность

 

 

 

 

 

 

 

Q

Г 2

,

,

 

.

2

 

/

 

 

л

\

 

 

^

-

 

Ус =

-ШГ

[ 1 7

( + я

-

*> -

1Г,

( -

Т

)

 

или

 

 

 

 

 

 

_

 

Q

л

_

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

і/л —

Ус — "2л" —

-

 

 

 

 

 

 

Обозначим

2уд

=

Л; 2^/с =

d; тогда

Q =

vid.

Следовательно,

Л — d =

d.

Откуда находим коэффициент сужения струи

ц =

 

djh

=

0,5.

 

Форму

струи

находим

по уравнению

(265).

 

 

 

 

 

 

 

 

§ 31. ПЛОСКОЕ СВЕРХЗВУКОВОЕ ТЕЧЕНИЕ ГАЗА

Для анализа плоского течения газа воспользуемся уравнением движения идеальной жидкости в форме. Эйлера

dv

1

1

j

d7- =

 

f - - g r a ä >

Спроектировав

уравнение

Эйлера

на

оси

ординат,

пренебре­

гая силой тяжести

и считая

движение

установившимся,

получим

 

 

дѵх

 

î

dp .

 

 

 

дх '

 

p

dx

'

 

 

 

дѵу .

 

î

dp

 

 

 

1

дх

 

p

dy

 

 

 

 

 

 

 

 

 

Уравнение неразрывности для рассматриваемого случая за­ пишем в виде

Выше было показано, что только для безвихревого, т. е. по­ тенциального течения газа справедливо уравнение адиабаты

-4- = const. р к

268

После этих предварительных рассуждений можно построить уравнение, которому должен удовлетворять потенциал скорости

ср {х, у).

В самом деле, в уравнения движения и расхода входят частные производные от давления и плотности по независимым перемен­

ным X и

у.

 

 

 

 

 

 

адиабаты р =

р (р), получим

Имея в виду,

что

по уравнению

 

 

 

dp

_

dp

dp

_

q 2

gp .

 

 

 

 

 

дх

 

dp

дх

 

 

дх

'

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dp

_

dp

др

 

2

др

 

 

 

 

 

 

~ду

~

~dp~dy ~

 

ду

'

 

 

где а

скорость

звука в данной

точке потока.

 

 

Подставим эти значения

др/дх

и дріду

в уравнения движения:

 

 

 

дѴѵ

,

 

дѵх

 

 

о?

dp

 

,cicc\

 

 

v.. -~

+

v., -rr±- =

p

-J-

 

(266)

 

 

A

dx

1

y

dy

 

 

dx

v

;

 

 

 

dvy

 

 

dvy

_

_

ü l

І £ .

 

(267)

 

 

v*^7-rVy

'

"u

dy

 

 

 

 

dx

 

 

p

dy

 

 

Умножим уравнение (266) на vx, a (267) на vy и сложим их:

Заменим сумму справа из уравнения неразрывности, тогда, перенося все в левую часть и используя условие потенциальности, получим

(vi

-

а2 )

+ хѵу *Ъ

+

(4 - а2 )

=

0.

(268)

Но для

потенциального

течения

 

 

 

 

 

 

-r

ô x

^

<3</

 

 

 

Следовательно, окончательный вид уравнения для потенциала

скорости примет вид

 

 

 

 

 

 

 

(<4 -

а*) g

+

&

 

+ (4 - «'> £

-

».

0»)

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

Это уравнение нелинейно, так как его коэффициенты сами зависят от искомой функции ср (точнее, от ее частных производных

по X я у).

269

Соседние файлы в папке книги из ГПНТБ