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

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

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

Рассмотрение полученного графика приводит к следующим заклю­ чениям.

1.Изоэнтропы с большим значением s (две верхние линии на рис. 27) идут почти горизонтально, т. е. сохраняют в ходе изоэнтропийного процесса расширения довольно высокое и почти постоян­ ное, причем однозначно определяемое значение а я» 0,8.

2.Изоэнтропы изгибаются по мере снижения s и приближения этой линии в области критической точки К ■Изгиб увеличивается по

мере приближения изоэнтропы к точке К. - 3. Изгиб изоэнтропы, наблюдаемый в области, близкой к крити­

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

4. Значениям энтропии выбранных изоэнтропийных процессов, по данным диаграммы рис. 27, примерно соответствуют постоянные значения а, определяемые начальной и конечной точками процесса расширения (если эти точки достаточно удалены от критической):

при s =

5,5

a *==>0,7;

»

s ==

5,0

а

0,6;

»

s = 4,5

а ^

0,5.

Таким образом, можно считать, что эксперимент подтверждает закономерность прямой связи постоянной энтропии процесса рас­ ширения с коэффициентом а.

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

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

J76

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

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

(222)а = 1, т. е. случай, когда его уравнение состояния может быть принято в форме (204). Обычно такие условия в действительности соблюдаются при достаточно малых давлениях р или достаточно

высоких температурах Т.

В технической литературе 124], [23], 172] и других источниках, публикующих результаты экспериментальных исследований тепло­ физических свойств газов и паров, предлагается расчетный метод учета влияния давления на параметры и функции состояния допол­ нительно к учету влияния температуры. В случае газов, подчиня­ ющихся уравнению состояния (204), теплоемкости ср и cv являются функциями только температуры Т и не зависят от давления р. Реаль­ ные газы и пары в отличие от идеальных имеют двупараметричес­ кую зависимость теплоемкостей от температуры и давления:

 

ср= ср(Т,р) и cv = cY(T,p).

 

 

(223)

Сравнивая эти значения теплоемкостей с их значениями при а =

= 1, т. е. с величинами ср0

= срТ и cvm (Т), получаем

 

 

 

 

СрСро —|—Дер и Су —- Сусо Ф- Асу

 

 

(224)

В этих уравнениях через ср0 и cVa} обозначены теплоемкости

при

нулевом значении р и бесконечно большом значении

V, т.

е.

при

условиях,

обеспечивающих реальным газам и парам свойства идеаль­

ного газа,

у которого а =

1.

члены

Аср и

Из изложенного вытекает, что корректирующие

АСу будут только функциями давления р. Это значит, что зависимость теплоемкостей от температуры у реальных газов такая же, как у соот­ ветствующих им идеальных газов. Такое предположение, касающееся теплоемкостей, разбивает их на два слагаемых, из которых одно зависит только от Т, а другое только от р. Предлагаемый метод корректирующих членов требует по аналогии и вириальной формы уравнения состояния:

pv = RТ + Вр,

(225)

где второй вириальный коэффициент В — является только функцией Т [23 ]. Как показано в [23 ], можно с достаточной степенью точности выразить этот коэффициент через параметры критической точки ве­ щества рабочего агента. Обозначая их подстрочным значком «с», получаем

В-

9R Т с

54RT);

128р .

(226)

 

№ р сТ 2

12 В. К . Васильев

 

177

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

( дсу\

_

т (]Рр_\ ( dv_\

 

\ д р ) т

 

1 \ д Т 2 J v \ д р ) т

 

( fcp \

_

т / d*v \

 

\ д р ) т ~

\ д Т 2 ) р ’

 

откуда

 

 

 

Дс» =

-

( 7' ( 1 р ) / р

<227>

и

 

О

 

 

 

 

&c« = l T ( T w U i ) T dP-

!228>

О

 

 

 

Первые и вторые производные в (227) и (228) можно найти из урав­ нения состояния, пользуясь значением В из (226). Получим:

( д 2р \

2R

( д В \

1

Rт

/ д 2В

\

(229)

\ д Т 2 ) v ~

(у В ) 2 \ д Т ) v + (у — В )2 V, д Т 2

) у

 

и

 

 

 

 

 

 

 

\ д р ) т ~

RТ

'

V д Т 2 ) Р ~

V д Т 2 ) Р ■

 

Подставив эти значения производных в (227) и (228) и взяв зна­

чения производных

вириального коэффициента В из

(226),

найдем:

( д В \

2.54.R7^_

/ а » В \

6.54.R7^_

 

\ d T J v ~

128-рсГ3

\ д Т 2

) у

128- р сТ*

'

^ ^

 

А

J

128

р с

Т 3

И

 

 

 

 

 

 

 

 

Ас =

f

324

 

dP

 

 

(232)

 

Р

J

128

Рс

Т 3

 

 

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

178

агента. Экспериментально установлено значение коэффициента сжи­ маемости а для этой точки нормальных рабочих газообразных аген­ тов а = 9/32 = 0,2812, что дает

отсюда находим

(233)

Подставляя это значение R в (232), будем иметь

Acv= 3vcT2cpT 3 и Acp = 9vcT2cpT 3.

(234}

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

В табл. 15 приведены данные по параметрам критической точки газообразных рабочих агентов [23].

Совершенно так же можно, пользуясь уравнением состояния в форме (225), вычислить поправки At и As к значениям i'„ и s0, полу­ чающимся при использовании уравнения состояния (204) для а = 1. Имеем:

i = i о + At и s = s0 + As.

(235)

Корректирующие члены можно найти по зависимостям:

(236)

и

(237)

из которых получаем

р

(238)

о

12

179

 

 

 

 

 

 

 

 

 

Т а б л и ц а 15

 

 

 

Критические значения параметров

состояния некоторых веществ

 

Наименование вещества

Химическая

n ,

105 Па

 

t c, °C

 

формула

 

(кгс/см 2)

 

 

 

 

 

I. Н е о р г а н и ч е с к и е в е щ е с т в а

 

 

 

Азот

 

азота

N2

 

33,9

(34,6)

 

— 147,2

Окись

NO

 

65,4

(66,7)

 

—92,9

Закись

азота

N20

 

72,7

(74,1)

 

36,5

Четырехокись азота

n 2o 4

 

1 0 0 , 0

(1 0 2 )

 

158

Аммиак

 

NH3

 

113,0(115,2)

 

132,4

Аргон

 

 

Ar

 

48,6 (49,6)

 

1 2 2 , 0

Бор

треххлористый

BC13

 

38,7

(39,5)

 

178,8

Бор

трехфтористый

BF3

 

49,8

(50,8)

 

1 2 , 2

Бром

 

 

B r2

 

122,6 (125,0)

 

302,2

Водород

 

H 2

 

12,9 (13,2)

 

—240,0

Водород бромистый

H Br

 

85,5

(87,2)

 

90,0

Водород йодистый

H J

 

83,1

(84,7)

 

151,0

Водород хлористый

HC1

 

82,7

(84,3)

 

51,4

Водород цианистый

HCN

 

50,7

(51,7)

 

183,5

Водяной

пар

H20

 

221,28 (225,65)

 

374,2

Воздух

 

 

 

37,76 (38,5)

 

— 140,7

Гелий

 

 

He

 

2,29 (2,34)

 

—267,9

Кислород

четырехфтористый

o 2

 

50,4

(51,4)

 

— 118,9

Кремний

SiF4

 

37,2

(37,9)

 

— 14,2

Криптон

 

Kr

 

55 (56)

 

—63

Неон

 

 

 

Ne

 

27,3

(27,8)

 

—228,3

Олово четыреххлористое

SnCl4

 

38,4

(39,2)

 

388,7

Двуокись

серы

s o 2

 

78,7

(80,3)

 

157,2

Сероводород

H2s

 

90,1

(91,9)

 

100,3

Трехокись серы

s o 3

 

84,7

(86,4)

 

218,3

Хлор

 

 

 

Cl

 

77,1

(78,6)

 

144,0

 

 

 

I I . У г л е р о д и с т ы е

с о е д и н е н и я

 

 

Метан

 

 

с н 4

 

46,4

(47,3)

 

—82,1

Цианистоводородная кислота

CHN

 

50,8

(51,7)

 

183,5

Метиловый спирт

CH3 OH

 

79,7

(81,3)

 

240,0

Окись углерода

СО

 

35,0 (35,7)

 

— 138,7

Двуокись

углерода

c o 2

 

73,9

(75,4)

 

31,1

Сероуглерод

c s 2

 

73,8

(75,3)

 

273,0

Фосген

 

 

CC120

 

56,8

(57,9)

.

181,7

Четыреххлористый углерод

CC14

 

45,6 (46,5)

283,1

Ацетилен

 

C2H2

 

62,9

(64,1)

 

36,0

Этилен

 

 

c 2H 4

 

51,4(52,4) .

 

9,6

Уксусная

кислота

C2H4 0 2

 

58,0 (59,1)

 

321,6

Этан

 

 

 

C2H 0

 

48,8

(49,8)

 

32,3

Этиловый спирт

C2H 5OH

63,9

(65,2)

 

243,5

Дициан

 

 

c 2n 2

 

58,9

(60,1)

 

126,6

Пропадиен

c 3H 4

 

52,5

(53,5)

 

120,7

Пропин

 

 

C3 H4

 

53,5

(54,6)

 

128

Пропилен

 

СзН6

 

46,0

(46,9)

 

91,4

Ацетон

 

 

c 3H 6o

 

47,7 (48,6)

 

235,0

Пропан

 

 

СзНд

 

42,6

(43,4)

 

96,8

1,3-бутадиен

C4 H e

 

43,1

(44,0)

 

161,8

Бутан

 

 

C4 H10

 

36,5

(37,2)

 

152,8

2 -метилпропан

c 4H 10

 

37,5

(38,2)

 

134,0

180

П р о д о л ж ен и е т абл . 15

Наименование вещества

Химическая

Р с > 10е Па

 

 

формула

(кгс/см2)

*е’ °с

 

Бутиловый спирт

С4 Н 9ОН

49,0

(50,0)

287

 

Пентан

С5Й12

33,4

(34,1)

197,2

 

2 -метилбутан

с5н12

33,2

(33,9)

187,8

"

Бензол

Q H ,

50,8

(51,8)

290,5

Фенол

С6Н вО

61,3

(62,5)

419

 

Анилин

C6H,N

53,0

(54,1)

426

 

Циклогексан

СсН12

41,3

(42,1)

279,9

 

Гексан

QHl

30,0

(30,6)

234,8

 

4

 

 

Гептан

С7н1в

27,3

(27,8)

266,8

 

Октан

QHl8

25,0

(25,5)

296,2

 

; - 0s +

j

[ —

( -d Pf- f ) J

(239)

Следовательно, корректирующие члены уравнений (235) будут:

M = = \ [ v - T { - w ) P \ d P

(240)

 

AS =

1

{~w~)pdp-

(241)

 

 

 

Частную производную

 

можно

получить

из уравнения

состояния (225), откуда

имеем:

 

 

 

 

 

, = ^

+

В,

 

 

причем величина В берется по формуле (226).

 

Тогда

 

 

 

 

 

 

/ 5 о \

R , 2-54 R т с

 

\ д Т ) Р

р

+

128

р сТ 3

 

Л(

9

ЯТС

 

3-54

р сТ 2 р;

(242)

128

р с

 

128

As = R In р -

2.54

 

(243)

128 р сТ 3 Р-

181

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

Найдя по формулам (224) и (235) значения энергетических пара­ метров газа с учетом влияния давления, можно получить значение а

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

аот р при разных постоянных значениях t (как это сделано на рис. 26 для водяного пара).

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

применяемого

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

§ 23.

РАСЧЕТЫ КОЭФФИЦИЕНТА СЖИМАЕМОСТИ

 

ДЛЯ АТМОСФЕРНОГО ВОЗДУХА

Посмотрим, как можно определить термодинамические параметры при условии возможности использования уравнения состояния идеального газа (204). С этой целью удобно воспользоваться табли­ цами [72].

При составлении этих таблиц воздух был принят в таком объем­ ном составе:

Атмосферный а з о т ..................................

0,7901

К и с л о р о д ...................................................

0,2099

причем состав атмосферного азота принимался следующий:

А з о т ............................................................

0,9876

А р г о н ............................................................

0,0119

Двуокись углерода ..............................

0,0004

В одород.......................................................

0,0001

Соответствующий этому объемному весовой (массовый) состав таков:

Атмосферный а з о т ..................................

0,7680

К и с л о р о д ...................................................

0,2320

Воздух принимался сухим (без содержания водяных паров).

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

Тогда

p v

1 ,

9

Т с

Р_

 

R Т

" г

128

Г

Рс

'

Положив

 

 

 

 

 

а =

1

9

Т с

_Р_

(244)

128

 

 

 

Т

Рс '

 

182

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

Т с

р

(245)

т = -=- и я =

— ,

ТРс

после чего формула (244) запишется так:

а = 1 + ^ 8 т(1 —6т2) я.

(246)

Следует обратить внимание на то, что относительные параметры обозначены не так, как в табл. 6, поскольку в таком именно виде они использованы в таблицах [72], рекомендуемых нами для рас­ четов.

В формуле (246) коэффициент при я по изотермам постоянный

-y^g-T (1 — 6т2) = С.

(247)

Тогда

(248)

а = 1 + Ся.

Пользуясь приведенными выше соотношениями и, в частности, формулой (244), можно расчетным путем найти для газов зависимость а от Т и р. Для воздуха эти расчеты выполнены в табл. 16, где исполь-

Т а б л и ц а 16

Результаты расчетов для воздуха а = а ( Т , р ) по формулам (246)—(248)

Температу­

 

т —

6т2

1 - 6 т 2

С

а = 1 4- Ся

 

132,45

ра Т , К

Т С

Т

Т

 

 

 

 

 

 

 

 

 

273.15

0,48490

1,4108

—0,4108

—0,01397

0,9624

323.15

0,40987

1,0080

—0,0080

—0,00023

0,9994

373.15

0,35495

0,75594

0,24406

0,00609

1,0164

423.15

0,31301

0,58875

0,41215

0,00907

1,0244

473.15

0,27993

0,47017

0,52983

0,01043

1,0280

зована по изотермам прямолинейная зависимость а от я. Это при­ ближенный результат, и весьма желательно сопоставить его с экспе­ риментальными данными по теплофизическим параметрам газов. В труде [1001 опубликованы данные Гольборна и Отто о значениях а для воздуха и водорода (табл. 17). Примем их для сравнения с циф­ рами табл. 16.

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

из них в координатах р, а (рис. 28)

имеет р = 0 и а = 1,0000 для

всех изотерм. Другую точку прямой

можно выбрать произвольно,

причем

в табл. 16 она имеет для

всех изотерм общую абсциссу:

р = 100

1 0 '5 Па (см. табл. 16, столбец 6). На графике рис. 28 нане­

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

183

Таблица 17

 

 

Значения

а =

p v

 

 

 

 

~RT

 

 

 

 

 

 

 

 

 

 

1. Д л я

в о з д у х а

 

 

 

 

Температура t, °С

 

 

~ I05 П а

' ’

50

 

100

ISO

200

 

0

 

0

1 , 0 0 0 0

1 , 0 0 0 0

 

1 , 0 0 0 0

1 , 0 0 0 0

1 , 0 0 0 0

1 0

0,9945

0,9990

 

1 , 0 0 1 2

1,0025

1,0031

2 0

0,9895

0,9984

 

1,0027

1,0051

1,0064

30

0,9851

0,9981

 

1,0045

1,0078

1,0097

40

0,9812

0,9982

 

1,0065

1,0108

1,0132

50

0,9779

0,9986

 

1,0087

1,0139

1,0168

60

0,9751

0,9993

 

1 , 0 1 1 2

1,0172

1,0205

70

0,9730

1,0004

 

1,0139

1,0206

1,0243

80

0,9714

1,0018

 

1,0169

1,0242

1,0282

90

0,9704

1,0036

 

1 , 0 2 0 1

1,0279

1,0322

1 0 0

0,9699

1,0057

 

1,0235

1,0319

1,0364

2. Д л я в о д о р о д а

 

 

 

Температура t,

°С

 

 

105

п а

—100

- 5 0

 

50

100

200

 

—150

0

0

1 , 0 0 0 0

1 , 0 0 0 0

1 , 0 0 0 0

1 , 0 0 0 0

1 , 0 0 0 0

1 , 0 0 0 0

1 , 0 0 0 0

1 0

1,0032

1,0064

1,0064

1,0061

1,0055

1,0049

1,0039

2 0

1,0073

1,0130

1,0130

1 , 0 1 2 2

1 , 0 1 1 1

1,0098

1,0078

30

1 , 0 1 2 2

1,0199

1,0197

1,0183

1,0166

1,0148

1,0118

40

1,0180

1,0271

1,0265

1,0245

1 , 0 2 2 2

1,0197

1,0157

50

1,0245

1,0345

1,0334

1,0307

1,0277

1,0246

1,0196

60

1,0319

1,0422

1,0404

1,0370

1,0332

1,0295

1,0235

70

1,0402

1,0501

1,0476

1,0433

1,0388

1,0345

1,0274

80

1,0492

1,0584

1,0584

1,0496

1,0443

1,0394

1,0313

90

1,0591

1,0668

1,0622

1,0560

1,0498

1,0443

1,0353

1 0 0

1,0699

1,0756

1,0697

1,0625

1,0554

1,0492

1,0392

рассчитанные по табл. 16.

На том же графике сплошными линиями

нанесены изотермы,

рассчитанные по табл.

17 (по данным Гольборна

и Отто). Сопоставление штриховых и сплошных линий изотерм выяв­ ляет результаты допущений, сделанных при выводе формул (246)— (248) для воздуха.

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

184

газа следует находить по известному Значению универсальной газо­ вой постоянной pR = 8314,2 Дж/(кмоль • К):

р_н- R

и'

Молекулярные (атомные) веса, удельные плотности и абсолютные значения энтальпии и энтропии при t = 0° С для одноатомных, двухатомных и трехатомных неорганических газов даны в табл. 18,

Рис. 28.

Графики табличных данных значений а для воздуха.

заимствованной

из труда

[23]. Для воздуха принимаем р = 28,97

и получаем значение газовой

постоянной

 

R = - - -

g ^ =

° .2870 кДж/(кг-К).

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

Продолжим расчеты коэффициента а для атмосферного воздуха. Взяв из табл. 15 соответствующие параметры критической точки для воздуха:

Тс = 132,45 К и рс = 37,76-К Г 5 Па,

185

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