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

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

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

динамики для приращения энтропии ds имеем

.

d i

vdp

d s = ----

(206)

По уравнениям (204) и (205) можно определить значение энтро­ пии

(207)

То

где Т о и р 0 являются параметрами, при которых произвольно выбрано s = 0.

Особое значение имеет изоэнтропийный процесс, подчиняющийся уравнению изоэнтропы

pvk — const,

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

Для

идеального газа имеем:

 

 

 

 

 

ср cv — R.

 

 

Из

последних

уравнений следует

 

 

 

 

к — \ R и cv =

R.

Дифференцируя

уравнение (205),

найдем

 

 

 

di — cpdT-

k

R d T

 

 

 

1

 

или, с

учетом (204),

 

 

 

 

 

di

T d (pv).

 

(208)

(209)

(210)

(211)

(212)

(213)

Если считать теплоемкость ср постоянной, то постоянны также cv и к. Тогда

i =

cp ( T - Т0)

(214)

s = c „ m

( ^ - ) - R i n ( i ) .

(215)

Дифференцируя при указанных условиях уравнение (208), полу­ чаем

— +'& — = 0.

(216)

р 1 V

 

166

Если в изоэнтропийном процессе изменения состояния идеального газа в начальный момент его параметры были р х, Т г и а параметры в текущий момент р, Т и и, то из уравнений (204) и (208) можно по­ лучить зависимости для расчетов изоэнтропийного процесса идеаль­ ного газа:

 

_ р _ _ А i_\*.

JT_

/ j ' A * - 1.

 

А-1

 

J L — ( J L \ к

(217)

 

Pi

V V ) '

7 \

\ v )

тг

\ Pi /

 

Обозначив

изоэнтропийное

изменение

энтальпии

идеального

газа Ats,

из формул (214) и (217) будем иметь

 

к—11

 

 

 

 

 

 

 

 

Дц

с, ( Л - Т ) =

 

 

 

с „ п

 

к

 

(1 -

Т ) =

[ 1 - ( - £ - )

(218)

Подставив значение ср из формулы (211) и учитывая уравнение (204), найдем при изоэнтропийном расширении идеального газа

 

 

 

-1 1

М.

■ РА

(■ к)

(219)

 

§ 22. СВОЙСТВА И УРАВНЕНИЕ

СОСТОЯНИЯ РЕАЛЬНЫХ ГАЗОВ И ПАРОВ

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

Для идеального газа формула (204) дает

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

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

в виде

а = ^ = а(р,Т).

(220)

Установлено, что а зависит от природы реального газа или пара.

Внастоящее время, когда все новые и новые вещества исследуются

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

167

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

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

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

Обычно теплофизические лаборатории, занимающиеся указан­ ными исследованиями, обрабатывают и публикуют свои работы в виде таблиц физических свойств рабочих агентов. В частности, широко известны и используются в расчетах таблицы теплофизических свойств воды и водяного пара М. П. Вукаловича. Последнее их изда­ ние вышло в 1969 г. в Издательстве стандартов [22]. К ним прило­ жены диаграммы i — s для водяного пара с предельными парамет­ рами ра =1000 ат и t = 800° С. Как таблицы, так и диаграммы в рассматриваемом издании даны в системе МКГСС и в междуна­ родной системе СИ. Такие же таблицы и диаграммы i — s для воды и водяного пара имеются и в других изданиях у нас и за границей.

Следует отметить, что основным уравнением, характеризующим физические свойства рабочего агента, в таких таблицах и диаграммах

служит уравнение состояния вещества.

В указанных выше таблицах

Вукаловича это уравнение дано

в виде

 

pv =

RT

+ 2

Рг+1

(221)

Vi

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

вириальные

коэффициенты,

вычисляемые

по теории ассоциации, разработанной В. П. Вукаловичем и И. И. Но­ виковым.' Эта теория изложена в текстовой части таблиц. Наш труд базируется на таблицах Вукаловича.

Обозначив заключенный в квадратные скобки множитель в пра­

вой части (221) через а, можно

выразить уравнение состояния в виде

pv

= a R71,

(222)

где а есть коэффициент сжимаемости уравнения (220). Если найти зависимость а от давления и температуры водяного пара, пользуясь таблицами [22], то формула (222) даст те же результаты, что и фор­ мула (221).

Значения двухпараметрической зависимости величины а от дав­ ления и температуры можно определить из таблиц теплофизических свойств воды и водяного пара [22] в пределах давлений и темпера­

тур

этих таблиц.

Был выбран ряд давлений р (в барах): р = 20;

40;

60;

80;

100;

150;

200;

220;

250;

300;

400;

600;

800;

1000

и для

этого

ряда

составлено 14

вспомогательных

таблиц, причем

каждая из них содержит необходимые термодинамические параметры

168

 

 

 

 

 

 

 

 

 

 

Таблица 11

 

Термодинамические параметры

и энергетические функции

воды и ее пара в изобарных

процессах

 

 

 

 

 

 

(давление

100-106, Па)

 

 

 

Т емпература

 

 

 

 

 

Энергетические функции

 

 

 

 

Удельный

Энтропия S,

pv

Ts

и

- f

S

t , °с

г, К

объем и*102,

к Д ж / (кг-К )

м3/кг

 

 

 

 

кД ж /к г

 

 

 

 

 

 

 

 

 

 

 

50

323,15

0,10077

 

0,6989

217,8

10,077

225,849

207,72

18,13

8,05

1 0 0

373,15

0,10386

 

1,2992

426,5

10,386

484,796

416,11

68,69

58,30

150

423,15

0,10843

 

1,8312

638,1

10,843

774,872

627,26

147,61

136,77

2 0 0

473,15

0,11480

 

2,3176

855,9

11,480

1096,57

844,42

252,15

240,67

250

523,15

0,12406

 

2,7794

1085,9

12,406

1454,04

1073,49

380,55

368,14

300

573,15

0,13978

 

3,2494

1343,7

13,978

1862,39

1329,72

_ 532,67

518,69

320

593,15

1,924

 

5,7120

2782,0

192,4

3388,07

2589,60

798,47

606,07

340

613,15

2,147

 

5,8782

2882,1

214,7

3604,22

2667,40

936,82

722,12

350

623,15

2,242

 

5,9464

2924,2

224,2

3705,50

2700,00

1005,5

781,30

375

648,15

2,453

 

6,0930

3017,4

245,3

3949,18

2772,10

1177,1

931,78

400

673,15

2,641

 

6,2158

3098,5

264,1

4184,16

2834,40

1349,8

1085,66

425

698,15

2,812

 

6,3238

3172,5

281,2

4414,96

2891,30

1523,7

1242,46

450

723,15

2,974

 

6,4220

3242,2

297,4

4644,07

2944,80

1699,3

1401,87

500

773,15

3,277

 

6,5984

3374,1

327,7

5101,55

3046,40

2055,1

1727,45

550

823,15

3,561

,

6,7568

3500,4

356,1

5561,86

3144,30

2417,6

2061,46

600

873,15

3,833

 

6,9025

3624,0

383,3

6026,92

3240,70

2786,2

2402,92

700

973,15

4,356

 

7,1656

3866,6

435,6

6973,20

3431,00

3542,2

3106,60

800

1073,15

4,860

 

7,4050

4111,3

486,0

7946,68

3625,30

4321,4

3835,38

и энергетические функции рабочего агента для ряда температур °С:

t =

50;

100;

150; 200; 250; 300; 320; 340; 350; 375 400; 425; 450;

500;

600;

700;

800.

При постоянном давлении, равном одному из членов указанного ряда давлений, при всех значениях температур ряда t °С из таблиц [22 ] были выбраны термодинамические параметры v, s и i и рассчи­ таны энергетические параметры pv и Ts путем перемножения их сомножителей. Рассчитывались также функции и, / и g по формулам:

и = i pv, f = и Ts\ g — i Ts.

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

Для примера здесь дана одна из серии таблиц для каждого члена ряда давлений. Табл. 11 рассчитана при постоянном значении р = =■ 100 бар (10~5 Па), причем в столбцах приведены цифровые зна­ чения параметров и функций, а в строках — их значения для ряда температур, указанного выше. Четырнадцать таких таблиц дают необходимый материал для термодинамических расчетов. Пользуясь ими, можно выполнять расчеты процессов энергетических циклов и производить различные исследования, и не располагая табли­ цами [22].

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

х = х (Т, р).

Для примера на рис. 16 построены такие изотермы для двух­ параметрической зависимости энтальпии от давления и температуры. На этой же диаграмме можно построить такие же изотермы значений любой другой функции х = х (Т, р), измеряемой в кДж/кг.

Комплекс таблиц типа табл. 11 может быть использован и для исследования двухпараметрической функциональной зависимости

а = а(Т, р) =

.

Удобно в координатах р, а построить изотермы значений а для всего принятого температурного ряда или для части его. Для каж­ дого постоянного значения температуры надо найти лишь значение функции pv, кДж/кг, в зависимости от давлений. Это можно сделать, используя серию таблиц типа табл. 11, не производя никаких до­ полнительных расчетов. Тогда оказывается удобным зависимость а = а (Т, р) представить в виде единой таблицы (табл. 12). По ней

170

 

 

 

 

 

 

 

 

Таблица 12

 

Значения

коэффициента а

= а ( Т , р )

=

DV

для воды

и пара

 

 

к/

 

 

 

 

 

 

 

 

 

 

Температура

Энергети­

 

Давление р, 10s Па

 

 

ческая

 

 

 

 

 

функция

 

 

 

 

 

 

 

t , °С

г, К

R7\

20

40

 

60

80

100

150

к Д ж /к г

 

50

323.15

149,137

0,01356

0,02709

0,04061

0,05410

0,06757

0,10716

100

373.15

172,212

0,01211

0,02420

0,03625

0,04829

0,06631

0,09025

150

423.15

■195,288

0,01116

0,02229

0,03340

0,04447

0,05552

0,08304

200

473.15

218,363

0,01059

0,02114

0,03165

0,04213

0,05257

0,07853

250

523.15

241,439

0,02363

0,02073

0,03100

0,04122

0,05138

0,07657

300

573.15

264,514

0,94891

0,88993

0,82022

0,73342

0,05284

0,07814

320

593.15

273,745

0,95564

0,90595

0,84954

0,78380

0,70284

0,08075

340

613.15

282,975

0,96122

0,91867

0,87167

0,81901

0,75872

0,08652

350

623.15

287,590

0,96387

0,92423

0,88105

0,83313

0,77958

0,59877

375

648.15

299,128

0,96882

0,93565

0,90002

0,86171

0,82005

0,69652

400

673.15

310,665

0,97339

0,94494

0,91507

0,88352

0,85011

0,75612

425

698.15

322,203

0,97702

0,95244

0,92699

0,90055

0,87274

0,79748

450

723.15

333,741

0,97980

0,95871

0,93701

0,91448

0,89111

0,82924

500

773.15

356,816

0,98426

0,96834

0,95209

0,93538

0,91834

0,87398

550

823.15

379,892

0,98767

0,97543

0,96280

0,95016

0,93737

0,90460

600

873.15

402,968

0,99016

0,98062

0,97095

0,96107

0,95119

0,92650

700

973.15

449,119

0,99395

0,98771

0,98179

0,97578

0,96990

0,95487

800

1073,15

495,271

0,99622

0,99259

0,98867

0,98500

0,98128

0,97220

 

П р и м е ч а н и е . В расчетах

взято значение R =

0,46151 кД ж /(кг- К).

Ф унк­

ция ри от Т и р берется из таблиц изобар, форма которых показана на примере табл. 11

при р — 100-10~6, Па.

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

лений логарифмический

масштаб.

Из рис. 26 следует,

что изотермы значений а сходятся к точке

с координатами а = 1

и р = 1, проходя в значительной степени

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

а,

соответствующими

началу и концу

фазовых превращений (х =

=

1

и х = 0). Удельные объемы в этих точках даются в таблицах

[22],

давление ps в них известно, и легко рассчитываются значения

а '

в точке х = 0 и а"

в точке х = 1.

Форма такого расчета дана

в табл. 13.

 

 

171

■ч»

t o

800° 700 600 550 500 450 400°С

Рис. 26. Изотермы коэффициента а = а ( Т , р ) для воды и пара.

 

 

 

 

 

Таблица 13

 

Расчет значений

коэффициента а на линии

насыщения

 

р, 10s Па

t , °с

т, К

RТ , кД ж /кг

м8/к г

К".10ь,

м*/кг

180

356,96

630,11

290,802

0,18380

0,7534

90

303,31

576,46

266,042

0,14179

2,046

45

257,41

530,56

244,859

0,12691

4.402

27

228,06

501,21

231,313

0,12050

7.402

18

207,10

480,25

221,640

0,11678

11,031

9

175,36 '

448,51

206,992

0,11213

21,484

4,5

147,91

421,06

194,323

0,10885

41,410

2.7

• 129,98

403,13

186,026 .

0,10700

66,878

1.8

116,93

390,08

180,026

0,10579

97,775

0,9

96,71

369,86

170,694

0,10412

187,01

р, 105 Па

p v ' , кД ж /кг

p v ", кД ж /кг

а '

а"

180-

33,0840

135,774

0,11377

0,46689

90

12,7611

184,140

0,04797

0,69215

45

5,71095

198,090

0,02332

0,79058

27

3,25350

199,854

0,01400

0,86400

18

2,10204

198,558

0,00948

0,89586

9

1,00917

193,356

0,00529

0,93412

4,5

0,48982

186,345

0,00252

0,95894

2.7

0,28890

180,571

0,00155

0,97067

1.8

0,19042

175,995

0.00106

0,97761

0,9

0,09371

168,309

0,00052

0,98603

П р и м е ч а н и е . По данным этих расчетов на рис. 26 и 27 нанесены линии насыщения в координатах р —а.

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

а = a (s).

Это не противоречит рассмотренной выше функциональной зави­ симости а — а (Т , р), так как энтропия s сама является двухпара­ метрической функцией состояния рабочего агента:

s = s (Т, р).

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

173

4^

 

 

Расчеты изоэнтропийного процесса в координатах р

— а

 

Таблица 14

 

 

 

 

 

 

 

S =

6,0000 кДж/(кг-К)

 

 

 

 

р , 105 Па

t , ° С

Т , К

R Т, кД ж /к г

/, кД ж /кг

и*102, м3/к г

i-pv, кД ж /кг

а

X

900

766,75

1039,90

479,924

3628,74

0,45680

411,120

0,856635

450

620,22

893,37

412,299

3373,41

0,73591

331,160

0,803202

270

521,69

794,84

366,827

3215,26

1,07420

290,034

0,790657

180

452,25

725,40

334,779

3110,31

1,47335

265,203

0,792173

90

343,74

616,89

284,701

2983,78

2,51955

226,760

0,796483

45

257,41

530,26

244,859

2787,43

4,37884

197,048

0,804741

0,9946

27

228,06

501,21

231,313

2689,18

6,95269

187,723

0,811550

0,9383

18

207,10

480,25

221,640

2614,58

9,99973

179,995

0,812105

0,9055

9

175,36

448,51

206,992

2494,40

18,5514

166,963

0,806614

0,8628

4,5

147,91

421,06

194,323

2382,92

34,3817

154,718

0,796186

0,8298

2,7

129,98

403,13

186,026

2306,07

54,1462

146,195

0,785788

0,8093

1,8

116,93

390,08

180,026

2248,15

77,7257

139,906

0,777145

0,7947

0,9

96,71

369,86

170,694

2154,69

144,414

129,972

0,761432

П р и м е ч а н и е . Изоэнтропы в координатах

р

— а нанесены

на рис. 27.

Там же нанесены

изоэнтропы

s =

6,1142; 5,5000; 5,0000;

4,5000 кД ж /(кг- К), рассчитанные по таблицам, аналогичным данной.

На рис.

27

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

К ,

параметры ^которой за ­

фиксированы

по таблицам [22], рекомендуемым для

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

Эти параметры;

647,30 К;

Рс = 221,2* 10 5 Па, t c

= 374,15° С;

R = 0,46151 кД ж /(кг. К). Рассчитанная

по ним величина а с составляет

 

 

 

 

 

 

ре°с

221,2.0,00317-105

„ оои, ол

 

 

 

 

- с

 

Кгс

0,46151-647,3.103

’’ д '

 

 

 

энтропийных процессов расширения современных энергетических

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

циклы.

В пароводяных циклах значение постоянной энтропии процессов

расширения лежит в пределах 4,000 ■— 8,000

кДж/кг. Поэтому для

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

ной энтропии процессов

расширения были

приняты

значения

s, кдж/(кг-К): s = 4,500; 5,000; 5,500; 6,000; 6,114.

энтропии

По таблицам [22]

при

этих значениях

постоянной

составлялись таблицы

изменяемости давлений от 900 до

0,9 бар

Рис. 27. Линии постоянной энтропии в координатах а и р для воды

иее пара.

синтервалами: р = 900; 450; 270; 180; 90; 45; 27; 18; 9; 4,5; 2,7; 1,8; 0,9. Этот ряд выбран по изменяемости отношения давлений в про­

цессе расширения к начальному давлению

=

900 бар (10~5 Па):

— = 1; 0,5; 0,3; 0,2; 0,1; 0,05; 0,03; 0,02;

0,01;

0,005; 0,003; 0,002;

о д а .

По таблицам [22 ] при этих условиях определялись температура t, удельный объем v и удельная энтальпия i. По найденным данным рас­ считывалось значение произведений pv, кДж/кг, и RT, кДж/кг, и по

ним вычислялось значение а = - £ ~ . Примером таких табличных КI

расчетов может служить табл. 14 (при значении постоянной энтро­ пии s= 6,0000 кДж/(кг-К). Аналогично получены результаты рас­ четов и при других значениях постоянной энтропии.

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

175

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