Файл: Васильев В.К. Термодинамические основы исследовательского проектирования судовых энергетических установок.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 270

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Соответственно этому значению я 02 находим строку в таблице (непосредственно или путем интерполяции) и берем из нее значе­

ния s°, 12 и г2. Таким образом определим параметры точки 2. Если требуется знать параметры других точек процесса, то надо задаться давлениями в этих точках и для каждой из них проделать указанные операции.

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

Из уравнения (6) можно видеть, что при постоянной темпера­ туре (АТ = 0) при любом значении подынтегральной функции ин­ теграл, определяющий величину s°, будет равен нулю. Тогда фор­

мула (17) для изменения энтропии в и з о т е р м и ч е с к о м

п р о ­

ц е с с е будет

 

ASl_2 = — Rln-g-.

(20)

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

изотермам.

и In p 2lpi = 0.

Тогда урав-

Для

изобарного процесса p j p 1 = 1

ние (17)

дает

 

 

 

Asi_2 = Asi_2*

(21)

Следовательно, в и з о б а р н о м

п р о ц е с с е

приращение

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

Изучая энергетические циклы с рабочим агентом в виде идеаль­ ных газов, обычно задаются предельными изотермами с температу­ рами Т з (верхняя изотерма) и Ту (нижняя изотерма). «Карнотизация» действительного цикла установки показывает, что при идеализации всех процессов к. п. д. цикла зависит от отношения указанных тем­ ператур Т 3/Ту и непрерывно повышается с увеличением этого отно­ шения. Из этого следует, что надо стремиться к повышению Т 3 и снижению Ту. Обе эти тенденции ограничены внешними условиями работы установки: первая из них — уровнем развития техники и условиями эксплуатации, а вторая — свойствами природных охла­ дителей рабочего агента перед началом его процесса сжатия. То и другое проектировщики установки должны тщательно .учитывать,

24


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

§4. ВЫБОР МАКСИМАЛЬНОГО ДАВЛЕНИЯ р2

ВПРОСТОМ ЦИКЛЕ ГТУ

Что касается максимального давления р 2, то его назначение яв­ ляется делом проектировщика; следует выяснить целесообразность и возможности повышения этого давления. Учитывая сказанное о зна­

чении давления р ъ можно

отметить, что выбор отношения

давле­

ний p jp x и (давления р х)

дает и значение другого. Прежде

всего

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

Состав воздуха, принятый также

при составлении таблиц [72],

включает следующие компоненты (объемн. доли):

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

■ - 0,7803

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

- 0,2099

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

• 0,0094

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

■0,0001

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

-

- 0,0003

 

Ит о г о - •

-1,0000

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

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

Рассмотрим влияние отношения давлений р^/рх на идеальный простой цикл ГТУ с использованием табличных данных теплофизи-

25


 

 

 

 

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

установки

 

 

 

 

 

 

 

 

 

 

 

Расчет

Отношения давлений

 

Расчет изоэнтропийного

процесса

процесса

 

внешнего

 

 

(задаю тся)

 

 

 

 

сжатия 1—2

 

 

теплооб­

 

 

 

 

 

 

 

 

 

 

 

мена

1

 

2

 

3

4

5

6

7

8

 

9

 

 

 

 

 

 

 

 

Я

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

Я

 

 

 

ай

|

 

■i

 

 

 

 

 

 

£

-«Ч.

 

 

 

 

 

 

 

 

 

 

 

 

1 (-

 

 

 

*

С4 и

 

I U

 

 

 

о.|

 

 

и

d

сС

 

 

as

 

 

 

 

 

*

] 1

 

 

 

t i | Cl

 

 

 

о

as

X

 

II *

 

 

С

 

 

й о

 

5Й(М

 

 

-ЕС

ъ\ £

 

 

 

о сч

< К

 

СУа5

 

 

 

К

 

 

 

1

 

0,0000

0,00000

1,2014

15,00

288,30

6,6617

0,00

 

795,50

2

 

0,6931

0,19891

2,4082

78,06

351,66

6,8607

63,36

 

732,14

3

 

1,0986

0,31529

3,6123

120,96

394,96

6,9770

106,66

 

688,84

4

 

1,3863

0,39786

4,8164

154,39

428,89

7,0596

140,59

 

654,91

5

 

1,6094

0,46189

6,0205

182,13

457,13

7,1236

168,83

 

626,67

6

 

1,7918

0,51423

7,2246

206,01

481,61

7,1760

193,31

 

602,19

7

 

1,9459

0,55832

8,4287

227,09

503,19

7,2202

214,89

 

580,61

8

 

2,0794

0,59677

9,6328

246,04

522,74

7,2585

234,44

 

561,06

9

 

2,1972

0,63058

10,8369

263,27

540,57

7,2923

252,27

 

543,26

10

 

2,3026

0,66083

12,0410

279,13

557,03

7,3225

268,73

 

526,77

1 1

 

2,3979

0,68818

13,2451

293,82

572,32

7,3499

284,02

 

511,48

1 2

 

2,4849

0,71315

14,4492

307,54

586,74

7,3767

298,74

 

497,06

13

 

2,5649

0,73611

15,6533

320,41

600,21

7,3978

311,91

 

483,59

14

 

2,6391

0,75740

16,8574

332,54

612,94

7,4191

324,64

 

470,86

15

 

2,7081

0,77721

18,0615

344,04

625,04

7,4390

336,74

 

458,76

16

 

2,7726

0,79572

19,2656

354,96

636,56

7,4574

348,26

 

447,24

17

 

2,8332

0,81311

20,4697

365,36

647,56

7,4748

359,26

 

436,24

18

 

2,8904

0,82952

21,6738

375,31

658,11

7,4912

369,81

 

425,69

19

 

2,9444

0,84502

22,8779

384,85

668,25

7,5068

379,95

 

415,55

20

 

2,9957

0,85975

24,0820

394,01

678,01

7,5215

389,71

 

405,79

 

П р и м е ч а н и е .

Параметры заданных точек 3 и / (см. рис. 1):

т о ч к а

3

760° С

к а

 

1 t l = 15°С;

Т 1 =

288,15

К; Sj =

6,6617 кД ж /(кг- К) = sj;

il = 288,30 кД ж /кг;

( s j =

Газовая постоянная для атмосферного воздуха R = 0,2870 кД ж /(кг-К ).

 

 

 

 

 

 

 

 

Таблица 1

при заданном значении предельных температур Т3/Тг (см. рис.

1)

 

 

 

 

 

 

 

Расчет

 

 

 

 

Расчет изоэнтропийного процесса

процесса

Расчет

полезной

работы

внешнего

 

расширения 3—4

 

 

 

теплооб­

 

 

 

 

 

 

 

 

 

мена

 

 

 

 

10

И

12

13

14

15

16

17

18

 

 

 

Я

 

 

 

а: м

айм

 

 

и

ай

 

 

 

<

-

<

 

 

as

h

7со и

Г

V)

1 б*

1

 

 

£

 

f-

Н «я

 

О

 

п

•-» *

 

< 'й

<

 

< £

 

as

as

 

 

С-1

о

н

* *

II *

Л *

1

||

5

неч

f-w

 

 

 

СУ*

< 3

 

Q

ТГ

 

 

 

 

 

 

129,85

760,00

1083,8

8,0051

0,0000

795,50

0,0000

0,00000

с о

64,925

593,33

895.33

7,8062

188.47

607.03

125,11

0,17088

7,9930

43,382

506.53

799.53

7,6897

284.27

511.23

177,61

0,25784

5,6303

32,462

449,48

737.53

7,6072

346.27

449.23

205,68

0,31406

4,8619

25,970

407,68

692,64

7,5432

391,16

404,34

222,33

0,35478

4,4978

21,642

375,06

657,86

7,4908

425,94

369,56

232,63

0,38631

4,3042

18,550

348.53

629,78

7,4466

454,02

341,48

239,13

0,41186

4,1818

16,231

326.31

606.41

7,4083

477.39

318.11

242.95

0,43302

4,1161

14,428

307.31

586,51

7,3745

497,29

298,21

245,02

0,45104

4,0813

12,985

290,74

569,14

7,3442

514,66

280,84

245,93

0,46686

4,0662

11,805

276,12

553,91

7,3169

529,89

265,61

245,87

0,48070

4,0672

10,821

263,05

540,35

7,2919

543,45

252,05

245.01

0,49292

4,0815

9,9885

251,29

528,19

7,2690

555,61

239,89

243.70

0,50394

4,1034

9,2750

240.60

517,16

7,2477

566,64

228,86

242,00

0,51395

4,1322

8,6567

230.83

507,11

7,2279

576,69

218,81

239.95

0,52304

4,1675

8,1156

221.84

497,84

7,2094

585,96

209,54

237.70

0,53148

4,2070

7,6382

213.53

489.33

7,1920

594.47

201.03

235,21

0,53918

4,2515

7,2140

205,81

481.41

7,1756

602.39

193.11

232,58

0,54636

4,2996

6,8342

198.61

474,01

7,1600

609,79

185,71

229,84

0,55310

4,3508

6,4925

191,87

467,07

7,1453

616,73

178,77

227.02

0,55945

4,4049

7\j = 1033,15

К; s3 =

7,4879 кД ж /(кг-К );

ig = 1083,8

кД ж /кг;

Sg = 8,0051

кД ж /(кг-К ); т о ч -

6,6617 к Д ж /(к г-К ).).

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

Исследование будет носить расчетный характер; отношение тем­ ператур Т 31Т1.считается заданным и неизменным, варьируется же отношение давлений p j p 1 с единичным интервалом, начиная с 1 и кончая 20, — пределы, которые ограничивают в настоящее время применение ГТУ простого цикла.

Результаты расчетов представлены в табл. 1. В первом и втором столбцах таблицы даны значения с интервалами, равными единице, отношений p j p 1 и их натуральных логарифмов. В третьем столбце даются значения произведения R In р г!рх при R = 0,2870 кДж/(кг-К) для атмосферного воздуха. Далее идет расчет простого цикла, изо­ браженного в диаграмме Т —s на рис. 1. Исходной точкой выбрана точка 1 начала процесса сжатия. От нее последовательно для каждого единичного интервала — отношения p j p i — с использованием таблиц 2 [72] рассчитывается изоэнтропийный процесс сжатия. Рас­ чет ведется по строкам таблицы. Цифры первой строки определяются по табл. 2 [72] при известной температуре Ц и являются исходными

26

27

 


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

Рис. 2. График изменяемости параметров простого идеального цикла ГТУ.

вого столбца таблицы. В результате этих расчетов заполняются 4, 5, 6 и 7-й столбцы таблицы. Далее рассчитываем столбец 8, определяю­

щий величину изоэнтропийного теплоперепада ДгТ_2 =

г2 — П, при­

чем величина

занесенная в первую строку столбца 6,

принимается

одинаковой для всех остальных строк этого столбца. Величина же г2 берется из строк столбца 6 как конечный параметр, отсчитываемый для каждого интервала от ее значения в первой строке. Этим закан­ чивается расчет процесса сжатия (точка 2 на рис. 1).

Следующим процессом цикла является изобарный процесс внеш­ него нагрева сжатого рабочего агента (23). Столбец 9 дает последо­ вательный расчет по интервалам внешнего теплообмена Qx = is — t2, причем i3 принимается по заданной температуре точки 3 из табл. 2

28

[72 ], ' оставаясь неизменной величиной при расчете всех строк столбца 9. Величина i2 (энтальпия конца процесса сжатия) берется из строк столбца 6 рассматриваемой таблицы. Цифры столбца 9 показывают, что по мере увеличения давления конца сжатия сни­ жается внешний теплообмен в нагреве рабочего агента. Это проис­ ходит из-за повышения температуры конца процесса сжатия (цифры столбца 5) при неизменности заданной температуры конца изобар­ ного нагрева (точка 3).

Характер зависимости Q = Q (p2/pi) показан на рис. 2, иллю­ стрирующем результаты расчетов в табл. 1. Там же показана линия роста изоэнтропийной работы компрессора по мере увеличения дав­

ления конца сжатия Ai$ = г2 — Ц.

Следующим процессом цикла является процесс расширения, иду­ щий изоэнтропийно от точки 3 до точки 4, определяемой конечной изобарой р х этого процесса. По ходу рассчитываемого процесса, те­ кущего в направлении, обратном направлению процесса сжатия, параметры рабочего агента от самых высоких их значений в точке 3 (первая строка рассматриваемой таблицы) снижаются до самых низких p2/Pi = 20 (предельной цифры расчетов). Параметры про­ цесса расширения рассчитываются для тех же интервалов отноше­ ния p 2/pi, что и параметры процесса сжатия (столбец 1 таблицы). Цифровые значения рассчитываемых параметров берутся из табл. 2 [72], по тому же методу, как при расчете процесса сжатия. Резуль­ таты расчетов занесены в столбцы 10—13 табл. 1.

Столбец 14 определяет зависимость изоэнтропийного теплопере-

пада турбины Ai] = is — t4 от отношения давлений p 2/pi в преде­ лах значения этого отношения от 1 до 20 через интервалы, равные единице. При этом величина i3остается во всех строках таблицы оди­ наковой, равной энтальпии заданной точки 3, а г4 берется из строк таблицы в столбце 12. Результат этих расчетов иллюстрируется ли­

нией Ai] = Ai] (p2/pi) (рис. 2). .

В столбце 15 рассчитывается внешний теплообмен

Qs = Qa (-§-) = *«— *‘i-

Значение г4 берут из столбца 12 таблицы, а г4 — по заданным параметрам точки 1 (см. рис. 1) начала процесса сжатия. Эта точка является вместе с тем точкой конца внешнего теплообмена Q2. На рис. 2 результат табличного расчета показан линией Q2 = Q2 (р2/рх)- Последние столбцы таблицы (16—18) содержат результаты опре­ деления энергетических показателей идеального цикла рассчитывае­

мой ГТУ. Полезная удельная работа Ai] Ai* f (pjpi) дана

в столбце 16, причем At] и г'з берутся соответственно из столбцов 14 и 8. К. п. д. теоретического цикла ГТУ:

 

Ai]

Ai] _ Qt — qz

Ъ ~

Qi

Qi

29


дан в столбце 17, а удельный массовый расход рабочего агента —• в столбце 18. Этот показатель является величиной, обратной полез­

ной удельной работе, и

определяется формулой

D

г/кДж.

Табличные результаты расчетов проиллюстрированы рис. 2. По­ лученные данные совпадают с результатами расчетов идеальных газо­

турбинных циклов, рассматриваемых, например, в § 4 труда

[70].

На основании расчетов можно сделать некоторые выводы.

р 21р1

1. Полезная

работа с

увеличением отношения давлений

растет от нуля

при p j p i

= 1 до некоторого максимума при

опре­

деленном значении отношения pJPi- Этот рост в указанных пределах увеличения pjp-i идет сначала быстро, затем, перейдя максимум, по­ лезная работа начинает монотонно и медленно снижаться. Отмечен­

ный здесь максимум зависит от отношения температур

Т а1Тг, и по­

этому, изменив

заданную величину отношения Т 3/Т г,

следует за­

ново выполнить

показанный здесь табличный расчет.

 

2. При том же отношении давлений, при котором получен макси­ мум удельной полезной механической энергии, имеет место и минимум массового удельного расхода рабочего агента D г/кДж. Рассматри­ вая линию D = D {pjpi) на рис. 2, можно видеть, что вблизи мини­ мума эта линия идет весьма полого. Минимум получается при p j p х =

=9,6, причем видно, что отклонения от этого отношения давлений

вту и другую сторону в пределах от 6 до 18 дают малые увеличения D, которыми можно пренебречь. Трудности конструирования компрес­

соров с высоким отношением давлений вынуждают не придержи­ ваться оптимального отношения р 2/Ръ а принять такое его значение, которое позволяет получить наиболее мощный компрессор. Можно, например, принять отношение давлений pJP i = 6. Однако при от­ ношении температур Т 31Т1, при котором рассчитывалась табл. 1, такое отступление от оптимального значения p%tpx снижает к. п. д. цикла с 45,9 до 38,6%, с чем приходится считаться.

3. Особенностью расчетов простого идеального цикла ГТУ яв­ ляется одинаковость отношения наивысшего давления к наинизшему в процессах расширения и сжатия. Это объясняется полным отсутствием энергетических потерь в цикле. В действительности же вследствие наличия потерь давления в газовых потоках отноше­ ние p 2lpi в процессе сжатия приходится принимать более высоким, чем в процессе расширения. Влияние сопротивлений в потоке рабо­ чего агента на совершение полезной работы в цикле весьма суще­ ственно и требует специального рассмотрения.

4. При анализе графиков рис. 2 обращают на себя внимание линии внешнего теплообмена Q1 и Q2- По мере увеличения p j p i обе эти величины падают. Особенно существенно снижение количества теп­ лоты Q1; сообщаемой рабочему агенту в камере сгорания. Это обстоя­ тельство при фиксированной температуре в точке 3 (начала процесса

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

р 2

в указанной точке неизбежно будет уменьшаться энтропия s3.

Так

30