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

Белозеров В.И. Учебное пособие по курсу Техническая термодинамика (исправлено)

.pdf
Скачиваний:
97
Добавлен:
04.12.2020
Размер:
4.85 Mб
Скачать

Åñëè w <<w , òî

12

 

 

 

 

 

ª

 

 

 

 

k 1

 

 

k

 

 

§

P2

· k

 

 

 

«

 

w

2

 

 

Pv

1

 

¨

 

¸

 

 

 

 

 

2

 

k 1

1 1

«

 

P1

 

 

 

 

«

 

©

¹

 

 

 

 

 

 

¬

 

 

 

 

 

º

»». (12.4.3)

»¼

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

w тем больше, чем меньше величина отношений давлений.

2

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

V v2G,

где v – удельный объем в выходном сечении сопла.

2

С другой стороны, величина V может быть определена как

V Sw2 ,

где S – площадь выходного сечения сопла, откуда получаем

G Sw2 .

(12.4.4)

v2

Заменяя в этом соотношении v с помощью уравнения адиабаты,

2

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

1

1 §¨ P2 ·¸k 1 , v2 © P1 ¹ v1

получаем

 

 

 

 

 

1

 

 

 

Sw

§ P

·

 

 

 

 

k

 

G

 

2

¨

2

¸ .

(12.4.5)

v1

 

P1

 

 

©

¹

 

 

 

Подставляя сюда значение w из уравнения (12.4.3), имеем

2

 

 

 

 

ª§ P

2

 

 

 

k 1

 

k

 

P

·

 

§ P ·

 

 

 

 

k

k

G S 2

 

 

1

Ǭ

2

¸

 

¨

2

¸

 

 

k 1 v

P

P

 

«

¹

 

©

¹

 

 

 

 

1

«©

1

 

1

 

 

 

 

 

 

¬

 

 

 

 

 

 

 

 

º

»

,

êã

.

(12.4.6)

 

»

 

c

 

»

 

 

 

 

¼

 

 

 

 

161

G

 

 

 

С помощью этого уравнения

 

 

 

 

 

 

 

 

 

 

можно решить и обратную задачу –

 

k

 

 

по заданным расходу, начальным и

G

 

 

 

конечным параметрам газа найти

max

 

 

 

 

 

 

 

площадь выходного сечения сопла.

 

 

 

 

Анализ характера зависимости

 

 

 

 

расхода G, задаваемой уравнением

 

 

 

 

(12.4.6), от величины Р /Р (обозна-

 

 

 

 

Ñ

1

 

 

 

 

 

чим эту величину через ), показы-

 

 

 

 

 

 

 

0

 

 

 

вает, что эта зависимость имеет

 

(P /P )

1

 

 

 

 

2

1 êð

 

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

 

 

 

 

Ðèñ. 12.4.2

 

(кривая 1-k-0). Очевидно, что при

 

 

 

 

= 1 G = 0. При уменьшении ве-

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

значении . Затем, согласно (12.4.6) G уменьшается, обращаясь в

íîëü ïðè = 0.

 

 

 

 

 

Экспериментальные исследования показали, что при 0 δ δ êð

уменьшение давления среды за соплом не влияло на величину рас-

хода газа через сопло; расход газа G оставался постоянным. Для

объяснения этого расхождения теории с экспериментом в 1839 г. Сен-

Венан выдвинул гипотезу, что при расширении газа в суживающем-

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

давления истечения P , соответствующего максимальному расхо-

 

 

 

êð

 

 

 

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

за соплом, меньших P , давление газа в выходном сечении остает-

 

 

 

êð

 

 

 

ся постоянным и равным P . Дальнейшие исследования доказали

 

 

 

 

êð

 

 

 

 

 

 

правильность этой гипотезы.

 

 

 

 

 

 

Сказанное иллюстрируется графи-

 

 

 

 

ком, представленным на рис. 12.4.2.

P

 

 

 

Как видно из рис.12.4.3, при P > P

2

 

 

 

 

 

 

 

 

2

êð

 

 

 

 

давление газа в выходном сечении со-

 

 

 

 

пла равно давлению среды, в которую

 

 

 

 

истекает газ из сопла. При P < P

äàâ-

 

 

 

 

c

êð

 

 

 

 

 

ление газа в выходном сечении сопла P

 

 

 

 

 

 

2

 

 

 

 

остается постоянным и равным крити-

P

 

 

 

ческому давлению P .

 

 

êð

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

êð

 

 

45

î

 

P

Выясним, что происходит в сопле при

 

 

 

 

 

 

 

 

ñ

 

 

 

 

 

 

 

снижении давления среды до величины

Ðèñ. 12.4.3

 

P . С этой целью исследуем уравнение

 

 

 

 

êð

 

 

162

 

 

 

 

 

 

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

2

 

2

1

 

 

k 1

 

1

 

 

 

\k

 

 

\k

0.

 

 

k

k

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначая величину \ для G = G

через \ , имеем

 

 

 

 

 

 

 

 

 

ìàêñ

êð

 

 

 

 

 

 

§ 2

·

k

 

 

 

 

\êð

k 1

.

 

 

 

¨

 

 

¸

 

 

 

(12.4.7)

 

 

 

 

 

 

 

 

 

 

 

© k 1

¹

 

 

 

 

 

Чтобы получить из (12.4.3) выражение для скорости истечения газа из сопла при максимальном расходе, нужно подставить в уравнение (12.4.3), найденное по уравнению (12.4.7), значение \ :

 

 

 

 

êð

c w

2

k

Pv .

(12.4.8)

 

êð

 

k 1 1 1

 

Соответственно выражение для величины максимального расхода через сопло получим из (12.4.6) с учетом (12.4.7):

 

 

 

 

 

 

 

 

2

 

 

 

 

k

 

P1

§

2

 

·

 

 

 

 

 

 

 

k 1

 

Gìàêñ

S 2

 

 

 

¨

 

 

¸ .

(12.4.9)

 

 

 

 

 

 

 

k 1 v1 © k 1

¹

 

 

 

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

1C

Р /Р с величиной \ , определяемой уравнением (12.4.7). Если

Ñ 1 êð

 

P

!\ ,

 

 

C

то P =P и расчет истечения следует вести по уравнени-

 

P1

 

êð

2 C

 

 

 

ÿì (12.4.3) è (12.4.6).

 

P

\ ,

 

 

 

 

 

Åñëè

C

òî P = P = P

и для расчета следует применить

P1

 

êð

2

C

êð

 

 

 

 

 

 

 

 

уравнения (12.4.8) и (12.4.9).

 

 

 

 

 

P

 

 

 

Hаконец, если

C

\ ,

òî P = P

> P и расчет ведется также

 

 

 

 

êð

 

2 êð

C

 

 

 

P1

 

 

 

по уравнениям (12.4.8) и (12.4.9).

163

âûõ 1 2

Расчет скорости истечения реальных газов, как было отмечено ранее, проводится по уравнениям (12.1.7) или (12.1.8) с помощью hS-

диаграмм или P-, v-, T-данных. Площадь выходного сечения сопла подсчитывается по уравнению (12.4.4), которое справедливо для течения любых сред.

12.5. Адиабатное течение с трением

Рассмотрим истечение газа с учетом трения о стенки канала. Процесс течения считаем адиабатным. Вместе с тем очевидно, что

этот процесс необратим – при течении выделяется тепло трения (q )

òð

и энтропия потока увеличивается:

dq

dS

òð

 

.

(12.5.1)

T

Рассмотрим, как изображается процесс истечения с трением на

hS- и TS-диаграммах (рис. 12.5.1, 12.5.2).

Если бы истечение было обратимым, без трения, то процесс изоб-

ражался бы на hS- и TS-диаграммах отрезком изоэнтропы S = S =

1 2

= const, заключенным между изобарами P и P (между точками 1 и

1 2

2), а скорость на выходе из сопла w определялась бы как (h – h ). Вследствие необратимых потерь при трении энтропия газа в процессе истечения возрастает и действительная адиабата отклоняется от изоэнтропы вправо (рис. 12.5.1). Далее, поскольку расширение газа в потоке с трением и без трения происходит до одного и того же

h

 

 

1

 

P = const

 

1

 

 

h

 

 

1

 

 

 

 

P = const

 

 

2

h

 

 

 

 

h

 

 

2

2

 

 

 

S

S

S

2

Ðèñ. 12.5.1

 

T

 

 

 

P = const

1

 

1

 

 

T

 

 

1

 

 

 

 

P = const

 

 

2

T

 

T

 

 

2

 

 

2

 

 

I

II

S

Ðèñ. 12.5.2

 

164

давления на выходе из сопла p , то очевидно, что точка действитель-

2

ного процесса будет лежать на той же изобаре, но правее точки 2

(точка 2д относится к истечению с трением), поскольку S > S .

2ä 2

Поскольку изобары на hS-диаграмме имеют положительный наклон, то

h > h ,

2ä 2

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

при течении без трения. Поскольку w < w, то можно записать

ä

w Μw,

(12.5.2)

ä

 

ãäå Μ – так называемый скоростной коэффициент, величина которого меньше единицы. Величина Μ лежит в пределах 0,95–0,98.

Потеря энергии потока на преодоление трения (обозначим ее E )

òð

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

 

 

Eòð

 

w2 w2

,

 

 

 

(12.5.3)

 

 

 

 

 

 

ä

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

откуда с учетом (12.5.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

1 Μ2

 

w2

,

(12.5.4)

 

 

 

 

 

2

 

 

òð

 

 

 

 

 

 

 

 

[ – коэффициент потери энергии, тогда

 

ãäå 1 Μ2

 

 

 

 

E

 

 

 

[

w2

.

 

 

 

(12.5.5)

 

 

òð

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

2 h1 h2

 

 

 

и, естественно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ä

2

1

 

 

ä

 

,

 

 

 

 

w

 

h

h

 

 

 

(12.5.6)

величина E может быть выражена как

 

 

 

 

 

òð

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

h

h

 

.

 

 

(12.5.7)

 

 

òð

 

 

 

 

 

2

 

 

 

 

 

Hаконец, из (12.5.5) получаем

165

'Eòð [ h1 h2 .

(12.5.8)

Приравнивая между собой правые части уравнений (12.5.7) и (12.5.8), имеем

h

h

[ h

h

.

(12.5.9)

2

1

2

 

 

Коэффициенты M и [ для реальных каналов, конечно, не могут быть найдены (определены) термодинамическим путем.

12.6. Общие закономерности течения.

Закон обращения воздействий

Рассмотрим теперь уравнение, описывающее наиболее общий случай течения:

wdw vdP gdz dl

dl .

(*)

òåõí

òð

 

Анализ этого уравнения позволяет сделать интересные заключе- ния о возможных способах ускорения потока.

Значение dP может быть выражено следующим образом:

§ wP · dP ¨ ¸ © wv ¹S

§ wP ·

dv ¨ ¸ dS. (12.6.1)

© wS ¹v

§ wP ·

Частную производную ¨ ¸ © ws ¹v

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

§ wx ·

§ wy ·

§ wz ·

1

с учетом ¨

 

¸

¨

 

¸

¨

 

¸

 

 

 

© wy ¹z

© wz ¹x © wx ¹y

 

§ wP ·

¨ ¸ © wS ¹v

поскольку

§ wv ·

¨ ¸ © wS ¹P

§ wP ·

§ wv ·

,

 

¨

 

 

¸

¨

 

 

 

¸

(12.6.2)

 

 

 

 

© wv ¹S © wS ¹P

 

 

§ wv ·

§ wT ·

,

 

¨

 

 

¸

¨

 

 

¸

(12.6.3)

 

 

 

 

© wT ¹P © wS ¹P

 

 

à

§ wT ·

¨ ¸ © wS ¹P

§ wT ·

T

,

T ¨

 

¸

 

 

CP

© wq ¹P

 

166

то из (12.6.2) получаем

§ wP ·

¨ ¸ © wS ¹v

§ wP ·

§ wv ·

T

 

 

¨

 

¸

¨

 

¸

 

.

(12.6.4)

 

 

 

© wv ¹S ©wT ¹P

CP

 

 

С учетом этого соотношения находим из (12.6.1)

 

§ wP · ª

 

§ wv ·

T

º

 

dP

¨

 

¸ «dv

¨

 

 

 

 

¸

 

dS ».

(12.6.5)

 

 

 

 

 

 

 

© wv

¹S ¬

 

© wT ¹P

CP

¼

 

Из уравнения Лапласа c

 

v

2

§ wP ·

 

 

 

 

 

¨

 

 

 

 

 

¸

следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

© wv

¹S

 

 

 

 

 

 

 

wP

 

 

 

 

c2

.

 

 

 

 

(12.6.6)

 

 

 

 

wv

 

 

v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Входящий в (12.6.5) дифференциал dv выразим из уравнения не-

разрывности в дифференциальной форме

df

 

dv

 

dw

, тогда

f

 

 

 

 

 

 

 

 

 

 

 

 

v w

 

§ df

 

 

dw ·

,

 

 

 

 

 

 

 

dv v¨

 

 

 

¸

 

 

(12.6.7)

 

 

 

©

f

 

 

w ¹

 

 

 

 

 

 

 

 

а дифференциал энтропии

 

 

 

 

 

 

 

 

 

 

 

 

 

dS

dq

 

dq

 

.

 

 

 

 

 

 

âíåø

òð

(12.6.8)

 

 

 

 

T

С учетом (12.6.6)–(12.6.8) получаем из (12.6.5)

vdp c

2

ªdf

 

dw

 

1

§ wv ·

dq dq

òð º

 

âíåø

 

 

«

 

 

 

¨

 

¸

 

 

».

(12.6.9)

 

 

 

 

 

 

 

 

 

¬

f

 

w v

© wT ¹P

CP

¼

 

Подставляя это соотношение для vdP в исходное уравнение (*) и группируя члены, содержащие dw, в левой части уравнения, получим

167

M 2 1

dw

 

df

 

 

1

§

wv

·

 

dq

 

 

1

dl

 

 

 

 

¨

 

 

 

¸

 

 

 

2

 

 

w

 

 

f

 

 

vCP © wT

¹P

 

g

c

 

 

 

 

1

§ wv ·

 

 

 

1

 

 

 

 

(12.6.10)

 

 

 

¨

 

 

¸

 

dqòð

 

 

 

 

dlòð

 

 

dz.

 

 

 

 

 

c

2

c

2

 

 

 

vCP © wT ¹P

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что dq

= dl , имеем

òð

òð

M 2 1

dw

 

df

 

 

1

§

wv

·

dq

 

 

 

 

 

 

¨

 

 

 

 

 

¸

 

 

 

w

f

vCP © wT

¹P

 

 

 

ª 1

 

§ wv ·

 

 

1

º

 

 

g

«

 

 

 

¨

 

 

 

¸

 

 

 

 

»dlòð

 

 

 

 

 

 

 

c

2

c

2

 

¬vCP © wT ¹P

 

 

 

¼

 

 

 

1

 

 

 

dlòåõí

c2

 

 

(12.6.11)

dz.

 

Проанализируем полученное соотношение. В случае, когда

q= 0, l = 0, dz = 0, уравнение (12.6.11) превращается в уже изве-

âíåø òð

стное нам уравнение

M 2 1

dw

 

df

.

w

 

 

 

f

Для случая течения в трубе постоянного сечения (df = 0) при

l= 0, l = 0 и dz = 0, но при наличии подвода (или отвода) тепла к

òåõí òð

потоку получаем из (12.6.11)

M 2 1

dw

 

1

 

§

wv

·

dq .

(12.6.12)

w

 

 

¨ ¸

 

 

 

 

 

âíåø

 

 

 

vCP © wT ¹P

 

 

P

 

 

§ wv ·

 

 

 

 

 

¨

 

 

¸

 

 

 

 

Поскольку всегда C > 0 и обычно

© wT ¹P

>0, òî èç (12.6.12) ñëå-

 

 

 

 

 

 

 

дует, что в дозвуковом потоке (М < 1) подвод тепла (dq

> 0) ïðè-

 

 

 

 

 

 

 

 

 

 

 

 

âíåø

водит к ускорению потока (dw > 0), отвод тепла – к его торможению. При подводе тепла газ в потоке расширяется, и его скорость увеличивается. Любопытно отметить, что дозвуковой поток жидко-

§ wv ·

сти, у которой ¨ ¸ < 0 (вода при t < 4°С), будет замедляться при

© wT ¹P

подводе к нему тепла.

Соответственно в сверхзвуковом потоке (М > 1) подвод тепла, как видно из (12.6.12), приводит к торможению потока, а отвод тепла – к

168

его ускорению.

 

 

 

Hа этих выводах основан

 

 

 

принцип устройства так называ-

 

 

 

 

 

емого теплового сопла. Тепло-

 

 

 

 

 

вым соплом называется труба

 

 

 

постоянного сечения, поток в

 

 

 

которой ускоряется за счет под-

Q

 

Q

âíåø

 

âíåø

вода или отвода тепла через

 

 

 

Критическое

стенки трубы (рис. 12.6.1).

 

 

 

сечение

До тех пор, пока скорость по-

 

Ðèñ. 12.6.1

 

 

 

 

тока не достигнет скорости зву-

 

 

 

 

 

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

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

ходящей в потоке газа.

Рассмотрим процесс течения в негоризонтальной (dz z 0) трубе

постоянного сечения при dq = 0, dl

= 0, dl = 0. Для этого слу-

âíåø

 

 

òð

 

 

òåõí

чая из (12.6.11) получаем

 

 

 

 

M 2 1

dw

 

 

g

dz.

(12.6.13)

w

 

 

 

c2

 

Из соотношения следует, что дозвуковой (М < 1) поток газа, движущийся вверх (dz > 0), ускоряется (dw > 0), а сверхзвуковой (М > 1) поток, движущийся вверх, замедляется (dw < 0).

Эти выводы представляют интерес для анализа процессов исте- чения природного газа из скважин (сечение которых постоянно по высоте).

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

dl = 0 и dz = 0, но при наличии потерь энергии на трение получим

òåõí

M

2

1

dw

ª

1

§ wv ·

 

 

 

«

 

¨

 

¸

 

 

w

 

 

 

 

 

¬vCP © wT ¹P

 

1 º

»dl . (12.6.14)

c2 ¼ òð

Из (12.6.14) следует, что при М<1 поток с трением ускоряется (dw > 0) до звуковой скорости, но перейти через скорость звука он не сможет, поскольку для этого нужно было бы отводить тепло от потока, а тепло трения всегда подводится к потоку. Hевозможность в

169

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

Упомянем еще об одном типе сопла – о так называемом расходном сопле, принцип действия которого состоит в следующем. Если ввести представление о плотности потока в канале как о величи- не расхода газа через единицу площади поперечного сечения канала

j

G

,

(12.6.15)

 

 

f

то из рассмотрения обычного сопла Лаваля следует, что в дозвуковой (сужающейся) части сопла j растет (df < 0 при G = const), достигает максимума в критическом сечении сопла и затем в сверхзвуковой (расширяющейся) части сопла уменьшается (df > 0 при

G= const).

Âтрубе постоянного сечения ( f = const) этого же эффекта можно добиться изменением расхода газа G путем вдувания или отсоса этого газа через отверстия в боковой поверхности трубы (так называемое расходное сопло, рис. 12.6.2).

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

Уравнение (12.6.11), позволяющее установить знак воздействия в зависимости от величины М, носит название закона обращения

воздействий.

Критическое

сечение

Ðèñ. 12.6.2

170