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

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

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

Итак,

 

 

 

w d t ~

= Jöl — - i - grad p ô / + coVô/ — 2 [йл. x ô»] wdt

или, имея в виду, что:

 

\. w dt -^- = wöw wxôwx

-\- Wybwy ~\- wzôwz ô ^ - ^ - j ;

2. 75/ =

*6x +

+ 0• 5z =

ô ( - v 2 " ^ 2 ) = ô ( - Ç ) ,

где вектор г имеет проекции „ѵ и у, так как перпендикулярен оси вращения z, а вектор о/ имеет проекции ôx, ôy и ôz;

3.

[ЮІ - Х М Л и; = 0,"••так

как вектор

xxw]

j _ w;

4.

/о/ = ôU — потенциал

массовой

силы;

 

5. grad pôf = у - Ôp,

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

Интегрирование этого уравнения вдоль линии тока от точки / до точки 2 (рис. 43) дает

 

9 О

О О

О О

 

 

2

 

 

 

 

 

 

1

Для

несжимаемой

жидкости

уравнение

энергии в форме Бер-

нулли при относительном движении имеет вид

 

ш-

и" ,

 

і

Р

.

 

+

+

— = const,

где и =

cor — переносная скорость

(окружная скорость точки ко­

леса, с которой в данный момент

совпадает частица жидкости);

U = —gy + с — потенциал

массовых сил.

 

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

тельного движения

примет вид

 

т

- т + т е г • - ! - « » » '•

<7 6 >

Применение уравнения Бернулли к анализу течения газа в лопаточных машинах

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

90

Рассмотрим течение газа в ступени осевой турбины, состоящей из соплового аппарата и рабочего колеса (рис. 44). Для средней линии тока можно приближенно записать, что

"1 = «2, т. е. частицы газа выходят из рабочего колеса на том же радиусе, на котором

они в

него входят.

 

 

Следовательно, параметры газа на входе / и на выходе 2 из рабочего колеса

будут

связаны уравнением

 

 

 

2 + k 1 ' Pi

2

k 1 ' p2 '

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

течения газа в ступени осевой турбины

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

В самом деле, из уравнения Бернулли следует, что при рг < pi относитель­ ная скорость W2 на выходе из рабочего колеса должна быть больше относительной

скорости wi на входе в рабочее колесо,

т. е.

W2 >

Oil.

Так как в обычных ступенях плотность р« газа на выходе немного ме ньше чем плотность рі на входе, а площадь проходного сечения немного больше, чем Fi, но так, что FiPi я« F2рг, то из уравнения расхода G = piFiWia = p-iF2W2a следует, что осевые скорости газа на входе и на выходе из рабочего колеса при­ мерно равны одна другой, т. е. шіаг=»Ш2а-

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

В сопловом аппарате также уменьшается давление и увеличивается скорость газа; в данном случае имеет место соотношение ѵоа як ѵіа (индексом 0 отмечено сечение на входе в сопловой аппарат).

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

сора. Только в этом случае надо учесть, что в компрессоре давление повышается,

91

т. е. рг >

pi, и следовательно, Дог <СДОі-ИЗ уравнения расхода также следует, что

w-iaœ

w \ a .

Таким образом, приходим к выводу, что в рабочем колесе компрес­

сора

происходит поворот потока с замедлением. Такая решетка называется диф-

фузорноіі.

 

 

 

На рис. 46 изображены решетки профилей компрессорной ступени и треуголь­

ники скоростей для рабочего колеса.

 

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

только от разности квадратов относи­

тельных скоростей юг и ші.

 

Для

радиальной

турбомашины

(центростремительной

турбины или

 

 

 

Рис.

46. Треуголь­

 

 

 

ники [скоростей н

Рис. 45. Векторы скоростей решетки

 

форма решеток осе­

соплового аппарата осевой ступени

 

вой

компрессорной

турбины

 

 

ступени

 

центробежного

компрессора) в уравнении

Бернулли сохраняется

член

с окружной

скоростью и:

 

 

 

 

 

 

2

k 1 p2

 

 

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

 

ttlj — U>2

V Р а

Р і /

 

[ ( o u f — и>1) + (и\ — и\)].

Поэтому в одной ступени радиальной турбомашины можно получить (срабо­ тать) больший теплоперепад, чем в одной ступени осевой турбомашины.

92

Уравнение Бернулли для течения с трением и теплообменом

Выведенное ранее уравнение Бернулли было получено для идеальной жидкости при интегрировании уравнения Эйлера вдоль линии тока.

Рассмотрим течение газа в канале при наличии трения и тепло­ обмена через стенки канала (рис. 47). Будем считать, что течение одномерное и установившееся.

Пусть за время At масса газа из положения /—2 переместится в положение )'—2'. Применяя закон сохранения энергии к этой массе газа, запишем

AQ = Де + ALßl! + AK,

(77)

где A Q количество

теплоты,

подве­

 

денное

к газу

за время

At;

Де = е 2

">

г1 — изменение

внутренней

энергии

Л

массы газа за время At;

A L B I I

работа,

р

совершенная

газом за время

Дг!; АК

=

 

= /Сг— Кі — изменение

кинетической

 

энергии

массы

газа

за

время

At.

 

 

Так

как

движение

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

Рис.

то параметры газа в объеме

Г2'

не

 

изменились

и

 

 

 

 

 

 

 

 

2

2'

47.

Движение

газа

в

канале

 

е 2

Еі

G At

вТа

с-/Га);

 

K2—Ki

= GM

(vi

—vf);

 

 

ALm

=

A L T

p +

p 2 / 2 A / 2

pxfx Al,

±

ALe;

 

 

 

 

 

 

 

AQ=(±AQBH4rAQTp),

 

 

 

 

 

 

где G — расход газа

в кг/с;

+ Д ^ е

— полезная работа,

совершен­

ная газом за время А^.

 

 

 

 

 

 

 

 

GAt,

Подставляя эти равенства в выражение

(77)

и деля на

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

±

QB„ + QTP = №

-

сѵТг) +

L T P

+

 

 

 

 

 

•P2J /2 A4

•Pi

GM

-L(vl-vV)±Le

 

 

 

 

 

 

 

GM

 

 

 

 

 

 

 

 

Так как

Q^J == у

>а

{^ѴТ + ~ ^ " )

представляет

собой

тепло­

содержание

і, то

 

 

 

 

 

 

 

 

 

 

 

±QBH + <2тр LTp = (k — h) + -g" (Vl — ѴЬ ± Le-

Учитывая, что работа

трения

переходит в теплоту, т. е. A Q T P =

= A L T P , имеем

следующее

выражение

для

уравнения

энергии:

±

<?вн = (І2

I I ) +

"n- (V2

Vi)

± L e

(78)

93

Если же течение газа происходит без энергообмена с внешней

средой,

то L e = О, QB1I = 0 н уравнение (78)

принимает

вид

 

і + -д- =

const

 

 

или

 

 

 

 

 

 

Т

+ 2с„

const.

 

(79)

Это

уравнение по форме

не отличается от

уравнения

(72), ио

по содержанию они совершенно различны. В уравнение (79) вхо­

 

 

дит

действительная

температура

 

 

газа

Т,

в то время

как

в

уравне­

 

 

ние Бернулли для идеальной жид­

 

 

кости входит температура пзо-

 

 

энтропического

течения

ТІІ3.

На

 

 

рис. 48

изображены

два

процесса

 

 

расширения

от

давления

р1

до

 

 

давления

р2.

Точке

Г 2

соответст­

 

 

вует конец расширения

в действи­

 

Пиз

тельном

 

процессе, а точке

Тгм

 

конец

воображаемого

процесса

 

 

Рис. 48. Полнтропический и изо

расширения,

идущего

пзоэнтро-

пически.

Площадка

под

линией

энтропический

процессы расшире

ния в

Т—s-диаграмме

Т^Т,

равна

работе

сил

трения.

 

 

Уравнение

(79)

в

дальнейшем

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

под

Т

действительную температуру

газа в

потоке.

 

 

 

 

 

 

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

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

dQ: d(±QB!l^Qrp) = ds + pd(jr)

или, используя понятие теплосодержания,

à ( ± QBH - f QT p ) = dl — ~ dp.

Интегрируя это выражение от состояния / до состояния 2Г найдем

о

± QBII + QT p = с.„ 2 - Тх) - J -J- dp,

94

вычисляя интеграл

при политропическом

процессе между

pup

и

 

 

kR

 

 

 

 

 

 

 

 

 

 

 

заменяя

ср =

k _ {

,

получим

 

 

 

 

 

 

 

 

 

*

 

-

Т і ) -

ІГЬТ

R (Т,

-

Т,)

= L T P ±

QB „,

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Il

 

k

^ т р

±

QDH

 

 

(80)

 

 

 

п — 1 — k — 1

Я (7Y— Гі) "

 

 

 

 

 

 

 

Работу сил трения подсчитывают по формуле

 

 

 

где £ — коэффициент

потерь,

определяемый

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

V — скорость

газа,

к

которой

отнесен

коэффициент

потерь

t..

Для

процесса

расширения

Тг < Тг

 

и п < /г, для

процесса

сжатия

Г 3 >

7\

и /г >

k.

написанного

для

течения

с трением,

Из уравнения

Бернулли,

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

 

71! =

Т2.

При этом имеется в виду, что вся кинетическая энергия струйки,

т. е. энергия

поступательного и

вихревого движения, переходит

в теплоту при изоэнтропическом торможении потока.

Рассмотрим как изменяется давление торможения при течении

с трением без

энергообмена.

 

Из термодинамики известно, что рассеяние энергии, имеющее место при необратимых процессах (трение), оценивается измене­ нием энтропии.

Изменение энтропии газа, подчиняющегося уравнению состоя­ ния ~- = RT, выражается формулой

As = s, — s, == с. In — .

_Pj_

Так как параметры торможения связаны со статическими параметрами изоэнтропы уравнением, то

Р/Р* = Р*/Р**,

но

v

RT* '

95

тогда

P'

откуда

Далее, так как

7 i = г ; ,

а для реальных (необратимых) процессов As > О,

то для течения с трением давление торможения уменьшается, т. е.

РІ<р\-

Следовательно, течение с трением характеризуется постоян­ ством температуры торможения Т* и падением давления торможе­ ния р*.

Потенциальное течение газа всегда будет изоэитропическим.

Уравнение движения в форме Крокко

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

Для вывода уравнения Крокко рассмотрим уравнение (63) Громеко—Лэмба, уравнение состояния совершенного газа

-^- = RT

(81)

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

s =

Су In

h const.

(82)

Вычисляя оператор у

выражений

(81) и (82), получим

соот­

ветственно

 

 

 

 

f

- f

ѴР =

«ѴГ

 

Исключая член с ур, имеем

после несложных

преобразований

с учетом уравнения (79)

 

 

 

f + TVs = T ^ T R V T

= T ^ T RWT* -

у 4 - .

(83)

96

Подставляя из полученного выражения - ^ - в уравнение Гро-

меко—Лэмба и приводя

подобные члены,

учитывая,

что ур =

= grad р, получим

 

 

 

 

 

g - + 2 (ш X Ъ) ="/ + 7 > -

j ^

j RyT*.

(84)

Полученное

выражение

называется

уравнением

движения

в энергетической

форме

или

уравнением

движения совершенного

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

вым

(особый случай винтового движения, когда

ш 11 w, в данном

случае не рассматривается).

 

§ 12.

УРАВНЕНИЕ СОХРАНЕНИЯ ЭНЕРГИИ

 

Закон сохранения энергии позволяет связать газодинамиче­

ские

процессы с одновременно протекающими

энергетическими

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

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

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

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

самого движения).

 

Энергию системы можно изменить с помощью

теплообмена

с окружающей средой вследствие выделения или

поглощения

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

Рассматривая изменение полной энергии системы, закон сохра­ нения энергии можно представить так: приращение полной энер-

7 B . C . Бекнев

97

гии Э энергетически неизолированного объема жидкости за еди­ ницу времени t равно работе L, которую внешние силы совершили над жидкостью, и притоку теплоты Q за ту же единицу времени:

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

Виды энергии, которые составляют полную энергию Э и ко­ торые следует принимать во внимание при рассмотрении процес­ сов в ГТ и КУ, таковы:

1)внутренняя тепловая энергия единицы массы (для перехода

квнутренней энергии единицы объема е надо умножить на

плотность p) de

cv dT;

 

2) кинетическая энергия -^- единицы массы.

В

условиях,

характерных

для рассматриваемых установок,

можно

пренебречь изменением

потенциальной энергии и даже

энергии излучения.

 

В величину полной энергии Э жидкости не следует включать энергию электромагнитного поля, имеющую объемную плотность

так как эта энергия присуща не веществу, а

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

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

Связь энергии электромагнитного поля с мощностью Ej дает формула Умова—Пойнтинга.

Для вывода этой формулы используем соотношение Максвелла rot Я

и получим, что Ej = Е rot Я

Первый член в правой части преобразуем при помощи вектор­ ного тождества: div (а X Ь) = b rot а — a rot b.

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

firot// = —div(£ X Я) + Я rot Е. По соотношению Максвелла

98

поэтому

 

 

 

£ rot 77 =

—div(£ X 77)

~

 

 

f'ma Ol

 

тогда окончательно получаем, что мощность

 

F / = - d i v ( £ x t f ) - ^ ( ^ + ^ ) .

Работа Е\, которую

электромагнитное

поле

совершает в еди­

ницу времени над электрическими зарядами, заключенными в еди­

нице

объема, слагается из скорости

уменьшения

объемной

плот-

ноети энергии электромагнитного поля

/ еа £2

,

 

В°-

\

и

из

у -

 

1

 

j

входящего в объем потока вектора Умова—Пойнтинга

Е

X

Н.

Положительным

направлением

вектора

ЕхН

считается

направ­

ление

из рассматриваемого объема.

 

 

 

 

 

 

 

 

 

В

условиях

рассматриваемых

КУ

энергия

электрического

 

8 Е2

 

энергии

-~

Вг

магнитного

поля.

поля —%— существенно меньше

 

 

стационарен

dB

п

то

 

 

 

 

 

 

 

Если процесс

—^- = 0,

 

 

 

 

 

 

 

Ё] = div(I X 77).

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

Работа вектора плотности поверхностной силы р„, действую­ щей на единицу площади за время dt по перемещению объема,

имеющего единичную поверхность, равна рп

(v dt),

где v dt —

путь, пройденный

этим объемом.

 

 

 

 

Работа вектора

плотности

массовой

силы

fm

по перемещению

единицы массы за

время dt

равна fm

(vdt)

=

-р—/V

{оdt).

При расчете ГТ необходимо рассматривать развиваемую или затраченную лопаточной машиной механическую (внешнюю) ра­ боту.

Механическую работу единицы массы в единицу времени обо­

значим через L m , тогда

механическая работа

L v единицы объема

в единицу времени

L v

= pLm

электромагнитного

Выделим работу

массовых (объемных) сил

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

7*

99

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