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

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

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

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

Основным результатом рассмотрения таблиц [22] в § 10 явился вывод о необходимости уравнений состояния для каждой однофазо­ вой области диаграммы состояния рассматриваемого рабочего агента, для двухфазовой области и для разделяющих их линий фазовых пере­ ходов. Сейчас уже нельзя удовлетвориться наличием одного уравне­ ния состояния водяного перегретого пара, которое считалось доста­ точным до настоящего времени, а приходится искать взаимосвязь параметров Т, р и v и в других фазовых состояниях Н 20.

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

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

Возвращаясь к последним работам по расчетам показателя адиа­

баты к,

остановимся на областях перегретого пара,

обозначенных

на рис.

13 цифрами 2 и 4. В работе

[25] в качестве уравнения состоя­

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

принято уравнение

 

 

a = - g r = l + Bp + Cp2 + Dp3.

(385)

Для расчета к используется уравнение (382), где значение произ­

водной ^ берется из уравнения (385). После подстановки по-

267

Лучаем выражение для

показателя а д и а б а т ы перегретого водяного

пара:

 

 

 

 

 

ср 1 4- Вр+Ср2+

Рр3

(386)

 

cv 2Dp3+

Ср2 — 1

 

 

Значение изобарной

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

ср

берется из таблиц [22],

а значение изохорной определяется из дифференциального соотно­ шения

cv ср -\-Т

(387)

Производные, входящие в это выражение, определяются по урав­ нению (385), и подстановка их в уравнение (387) дает

.

R[(B + T- f ) + (c+ 7'-3r)>’+ (D+ r f - ) ' ,S+ 7 ]2

Су Ср -\----------------------------------------------------------

j

C + 2 D p — ^

(388)

Значения коэффициентов в уравнениях (385)—(387) подобраны как функции температуры Т и даны в табл. 30. Там же приведены

Таблица 30

Значения коэффициентов и их производных по температуре в уравнении состояния водяного пара (385)

t , °с

ВЛО4,

4 ? - 107’

С-10е, бар-2

+ + Л0>'

£>•10»,

d T

бар*”1

бар-3

 

 

бар-1• К-1

 

бар-2.К-‘

 

бар"3*К~3

50

—290,92

+4502,9

—3670,4

112 907

—424,22

5562,75

1 0 0

— 144,19

1831,8

—844,04

23 718,2

— 192,86

3383,53

150

—80,890

849,94

—220,49

5 638,75

—76,621

1469,40

2 0 0

—50,012

440,13

—65,721

1 497,77

—28,632

577,05

250

—33,320

250,09

—22,648

449,26

1 0 , 2 2 0

215,31

300

—23,480

153,66

—8,9688

155,74

—3,4721

76,893

350

— 17,244

100,73

—3,8316

62,851

-1 ,1 1 5 0

26,004

375

— 14,953

83,241

—2,6217

41,893

—0,61843

14,727

400

— 13,050

69,602

1,7530

28,557

—0,33971

8,1617

425

— 11,450

58,812

1,1566

19,735

—0,18711

4,4125

450

10,091

50,160

—0,74317

13,703

—0,10542

2,3274

475

—8,9279

43,135

—0,45653

9,4722

—0,06276

1,2075

500

—7,9240

37,366

—0,25950

6,4564

—0,04047

0,63474

550

—6,2887

28,565

—0,039942

2,7344

—0,02130

0,24612

600

—5,0259

22,279

+0,044198

0,85350

—0,01060

0,20044

650

—4,0334

17,647

+ 0,062112

— 0,01406

—0,00047

0,20640

700

—3,2425

14,144

+0,051905

—0,32907

+0,00985

0,19811

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

268

й температуры. Эти данные приводятся в табл. 31. По ним на рис. 48 построены изотермы функциональной зависимости k = k (р, Т). Пользуясь табл. 31 или кривыми на рис. 48, можно при расчетах термодинамических процессов по известным параметрам р и Т при­ нимать значение показателя адиабаты k [25 и 26].

Рис. 48. Зависимость

показателя адиабаты k перегретого водяного

пара

от давления и температуры.

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

Исходной зависимостью в этой области остается формула (380),

но отыскание значения частной производной Должно базиро­

ваться на других физических свойствах, отличающихся от свойств

269

Таблица 31

П оказатель ад и аб аты k для перегретого водяного п ара

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

р ,

 

 

 

 

 

 

 

 

бар

425

450

475

500

550

600

650

700

 

0

1,286

1,285

1,281

1,276

1,271

1,265

1,260

1,255

1

1,285

1,282

1,280

1,276

1,271

1,266

1,261

1,255

5

1,285

1,282

1,279

1,276

1,271

1,266

1,261

1,256

10

1,284

1,282

1,279

1,276

1,271

1,266

1,261

1,256

25

1,282

1,281

1,278

1,276

1,271

1,267'

1,262

1,257

50

1,281

1,279

1,278

1,276

1,272

1,267

1,263

1,259

75

1,281

1,280

1,279

1,277

1,273

1,269

1,265

1,261

100

1,282

1,282

1,280

1,278

1,275

1,271

1,267

1,263

125

1,284

1,285

1,283

1,281

1,277

1,273

1,269

1,267

150

1,291

1,289

1,286

1,284

1,280

1,276

1,273

1,271

175

1,294

1,290

1,287

1,282

1,279

1,277

1,276

200

1,301

1,295

1,290

1,285

1,283

1,282

1,282

225

1,310

1,300

1,294

1,289

1,288

1,288

1,290

250

1,307

1,298

1,294

1,293

1,296

1,300

275

1,304

1,298

1,300

1,305

1,312

300

1,310

1,303

1,308

1,316

1,328

350

 

1,329

1,347

1,369

400

1,390

1,427

450

1,513

Р ,

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

бар

100

150

200

250

300

350

375

400

 

0

1,323

1,318

1,312

1,306

1,300

1,294

1,291

1,287

1

1,320

1,316

1,311

1,305

1,300

1,294

1,291

1,287

5

1,305

1,301

1,298

1,293

1,290

1,287

10

1,299

1,297

1,295

1,291

1 288

1,286

25

1,287

1,285

1,286

1,285

1,284

50

1,269

1,280

1,280

1,281

75

1,276

1,279

1,280

100

1,275

1,278

1,281

125

1,283

150

 

 

 

 

 

 

 

 

перегретого пара. В двухфазной области удельный объем рабочего агента суммируется наличием обеих фаз:

v = v' (1 — х)

+ v"x.

(389)

Покажем, что производная

в двухфазной области является

тоже аддитивной. Дифференцируя равенство (389)

по давлению при

s = const, получаем

 

 

+

 

<390>

где do' и do" — полные производные вдоль пограничной линии.

270

Степень сухости х может быть определена через энтропию:

откуда

?М£-<,-'>+£Ф <391)

Далее из (390) и (391) получим

( £ ) . - ( 3 - - £ £ ) о - М - £ - £ £ Ь <зэ2>

Полную производную от v по р вдоль линии насыщения напишем в виде

 

d v n

/ d v \ л ' н

/ d v \ л

d s n

(393)

 

d p

V d p ) s

\ ds ) р

d p

 

 

I d v \ л . Н

/ d v \л. н

 

 

 

Здесь Ы

и 1 ^ 1

— частные производные, взятые в точках

пересечения изобары с

линией насыщения, и надстрочный ин­

декс «л. н» относится к этим точкам.

 

 

 

Производные в формуле (393) могут быть

взяты со стороны одно­

фазной или двухфазной области. Условимся

брать их, идя от двух­

фазной области. Так как

 

 

 

( д о _ \

_

/ _ З Г \ _

d T

\ ds )

р ~ \ d p ) s

\ d p ) s ~

d p ’

то уравнение (393) может быть преобразовано следующим образом:

dv*

d T

d p

d p

< М 4 >

и уравнение (392) примет вид

d v

\"дФ

(395)

 

 

Таким образом, аддитивность величины

да s в двухфазной

области установлена, что существенно упрощает задачу получения простого аналитического уравнения для зависимости k = k (Т, х).

Действительно, подставив (389) и (395) в уравнение (380), получим следующее уравнение для показателя адиабаты в двухфазной области параметров состояния чистого вещества:

Ц Т , х ) = — -

v ' ( 1 х ) + v "x

(396)

(£)>-*>+ ( Э Т

 

Для расчетов по формуле (396) надо иметь значения производ-

ных ^ да удФ

^ ^ d v у'ДФ

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

д р Л

dp ) s

 

271

/ d v \ л .

н дф

\~ д р )

' '•^ТУ пРоизводную можно определить уравнением

р \ л-ндф

Су'нлф/ ^ Г \ а

(397)

V d p ) s

Т \ d p )

 

Вывод этой зависимости можно найти в [82]. Здесь через с £ ндф обо­ значена теплоемкость cv на линии йасыщения со стороны двухфаз­ ной области. Уравнение (397) позволяет довольно точно вычислять

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

/ d v \ 'дФ

( d v \"дФ

 

 

 

и \~gfi)

 

 

По изотерме в двухфазной области имеем

 

скдф — С1/ДФ— (и" v')

Т .

(398)

Вывод этого равенства найдем в [83].

Используя (398), из (397)

 

 

 

/ д о \ Л . Н д ф

получаем уравнение, связывающее между собой значения

(

на правой и левой частях линии насыщения:

 

( £ Г - ( £ Г = - ^ - » '> ( £ ) !£-

(339)

Это

уравнение

может служить

для взаимной увязки

величин

( d v \"дФ

/ d v \'дФ

> вычисленных

независимо одна от

другой.

 

и y-fyjs

Следует также отметить, что при переходе через линию насыще­ ния показатель k , подобно ряду других термодинамических величин, меняется скачком от величины кл- Н°Ф (однофазная область) до вели­ чины АлнДФ(двухфазная область). Скачкообразное изменение k при переходе линии процесса через линию насыщения обусловлено на­

личием разрыва производной

в точке перехода.

Величина скачка, который имеет в точке перехода производ-

ная

d v ^

может быть получена

по формуле

 

 

 

 

 

 

Ч £ Г Ч £ Г

*

«

°

 

°

>

 

 

 

 

 

 

г л . н оф

 

с р

Здесь сд ноф— изобарная теплоемкость в точке пересечения изо­ бары с линией насыщения со стороны однофазной области (надстроч­ ный индекс «оф»). Вывод формулы (400) можно найти в труде [83].

На базе полученных зависимостей (397)—(400) были выполнены расчеты в интервалах температур от 0 до 340° С, причем величины

( 7 ^ ) дФи

дф вычислялись по формуле (397). Производ­

ная dT/dp при этом определялась из соотношений для функциональ­ ной зависимости р = р (Т), которую можно найти в таблицах [22]

и других соответствующих источниках. Значение теплоемкости Су' н дф принималось по данным [83]. Полученные расчетами значения част­ ных производных даны в табл. 32,

272

Таблица 32

 

 

 

 

 

d v \ л.

н дф

 

 

 

 

 

З н ач ен и е

производн ой WJs

показателя ади абаты k воды и пара

 

 

 

 

 

на линии насыщения

 

 

 

 

 

 

 

 

Вода на линии насыщения

Сухой насыщенный пар

 

t , °с

/

dv

\'ДФ

к дФ

А/оФ-НГ3

/ dv \ ”ДФ

-

й"ДФ

*"°Ф

 

*

др

'&

- ( - а Д

 

 

мб/кг

 

 

м5/кг

 

 

 

0

0,7692-10-1

0,0002088

3602

2,9968

 

1,105

1,331

1 0

0,2149-10-1

0,0003720

1818

0,7704

 

1,104

1,332

2 0

0,6676-10“ 2

0,0006296

967

0,2190

 

1,108

1,330

30

0,2280

-10-2

0,001018

538

0,6845-10 '1

1 , 1 1 2

1,326

40

0,8467-10"3

0,001582

312

0,2328-10 "1

1,116

1,323

50

0,3388-Ю "3

0,002376

187

0,8547-10-2

1 , 1 2 0

1,321

60

0,1449- Щ- 3

0,003456

116

0,3364-10- 2

1,124

1,320

70

0,6592-10-4

0,004885

74,2

0,1410

-10"2

1,127

1,319

80

0,3176-10-4

0,006730

48,6

0,6244-10“ 3

1,130

1,317

90

0,1600-10"4

0,009056

32,6

0,2915- 10-з

1,133

1,314

1 0 0

0,8467-10“ 6

0,01193

22,3

0,1426

-10-3

1,136

1,311

ПО

0,4669-10“5

0,01542

15,5

0,7273-10“ 4

1,139

1,309

1 2 0

0,2674

-10-6

0,01959

1 1 , 0

0,3861-1 0 -4

1,141

1,306

130

0,1585

-ТО"6

0,02450

•7,88

0,2125

-10-4

1,142

1,303

140

0,9701 - 10~ 6

0,03019

5,71

0,1209

-10“ 4

1,142

1,300

150

0,6116-10“ °

0,03674

4,19

0,708210“ Б

1,142

1,297

160

0,3959

-10- 6

0,04418

3,11

0,4266

-10-5

1,141

1,294

170

0,2625- Ю- 6

0,05255

2,32

0,2632- Ю" 6

1,141

1,291

180

0,1781 • 10“ в

0,06192

1,74

0,1663- ю - 6

1,140

1,288

190

0,1234- Щ- 6

0,07230

1,31

0,1074-10-5

1,138

1,286

2 0 0

0,8711 - 10- 7

0,08373

0,997

0,7067- 10-е

1,135

1,284

2 1 0

0,6262-10“7

0,09625

0,761

0,473710- 6

1,131

1,282

2 2 0

0,4577-10- 7

0,1099

0,583

0,3229-10-6

1,126

1,279

230

0,3398

-10"7

0,1247

0,448

0,2235-10“ 6

1 , 1 2 0

1,275

240

0,2558

-10-7

0,1407

0,345

0,1569

-10-«

1,114

1,270

250

0,1952

-10-7

0,1580

0,266

0,1116-10-®

1,106

1,266

260

0,1509

-10“7

0,1766

0,205

0,8032

-10-7

1,096

1,262

270

0,1180

-10"7

0,1966

0,158

0,5845-10“ 7

1,085

1,257

290

0,7452- Ю' 8

0,2414

0,0924

0,3182-10-7

1,057

1,248

300

0,6013- ю - 8

0,2664

0,0697

0,2375

-10-7

1,040

1,243

310

0,4900-1.0-8

0,2934

0,0522

0,1784

-10“7

1 , 0 2 0

1,23

320

0,4034- Щ- 8

0,3228

0,0385

0,1348-10-7

0,9955

1 , 2 2

330

0,3359-10-8

0,3546

0,0277

0 , 1 0 2 0 - 1 0 " 7

0,9692

1 , 2 1

340

0,2826IQ' 8

0,3893

0,0197

0,7698-10

 

0,9402

 

Для расчетов k использовалась формула (396). Величины произ­

водных ДФи ( |~ ) ДФ> входящих в эту формулу, даны

в табл. 32. В той же таблице приведены значения показателя адиа­ баты &л- ндФна линии насыщения (надстрочный индекс «л. н») со стороны двухфазной области (надстрочный индекс «дф»). Этот пока­ затель рассчитывался при х = 0. При значении х = 1 рассчитан другой показатель, так же приведенный в табл. 32.

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

273

Кроме того, были рассчитаны для сравнения показатели адиа­ баты k’0^ и k"°Ф на линии насыщения со стороны однофазной области

(надстрочный индекс

«оф») в тех

же точках

по

формуле

 

 

Ф \ л'

 

ф

121

£Л. Н оф __

 

г л . н оф

(401)

 

К

d v ) т

д Т ■X

 

 

 

 

 

Рис. 49. Номограмма для определения показателя адиабаты влажного пара k ДФ с подходом от двухфазной области.

Эта формула была получена из (380) и (397), рассмотренных при пересечении линии насыщения линией процесса, идущей из одно­ фазной области.

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

Н°Фи

Н°Ф> необходимые

для расчетов по формуле (401), взяты из [98], а значение с£'ноф—

из [83]. Полученные в результате расчетов значения &л- н°4 также помещены в табл. 32. Следует отметить, что весьма большие значе­ ния k' объясняются малой сжимаемостью воды.

На рис. 49 даны изотермы зависимости k = k (Т, х) на основе расчетов, выполненных по формуле (396).

ГЛАВА III

ПРОИЗВОЛЬНЫЕ ПРОЦЕССЫ РАСШИРЕНИЯ* ОПТИМАЛЬНЫЙ ПРОЦЕСС И ЕГО

ТЕРМОДИНАМИЧЕСКИЙ РАСЧЕТ

§ 34. ТЕРМОДИНАМИЧЕСКИЕ СВОЙСТВА ПРОИЗВОЛЬНОГО ПРОЦЕССА РАСШИРЕНИЯ. ХАРАКТЕРНЫЕ ТОЧКИ ПРОЦЕССА

В процессе развития методов получения работы из теплоты теория турбин развивалась исключительно на базе адиабатного процесса расширения. Ход этого процесса и его свойства были заложены не только в основу теории и тепловых расчетов турбин, но и в конструк­ тивные формы «адиабатных» машин. Большие скорости потока рабо­ чего агента, текущего через проточную часть турбины, позволяли пренебрегать при надлежащей изоляции теплообменом потока с окру­ жающей средой и считать поток теплоизолированным в процессе расширения. Это значило, что кинетическая энергия в потоке созда­ валась только за счет соответствующего уменьшения внутренней энергии U!, которой поток обладал в начальной точке процесса расширения. Конструкция проточной части турбины обеспечивала непрерывный отвод кинетической (механической) энергии потока через лопаточные венцы вращающихся дисков ротора и через его вал на привод полезного сопротивления (генератор электрического тока, гребной винт и т. п.).

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

18*

275

Таким же образом развивалась теория И отыскивались оптималь­ ные решения задач выбора принципиальных схем турбинных энер­ гетических установок. Действительно, оптимальный цикл Карно в диаграмме Т— s изображался прямоугольником со сторонами, па­ раллельными координатным осям. Это было наглядно удобно и при­ вычно, пока не возник вопрос: всегда ли и при всяких ли внешних условиях адиабатный процесс расширения является единственным и бесспорно наилучшим?

На этот вопрос могут быть два ответа:

1.В энергетических установках адиабатный процесс оптимален,

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

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

позволяет

увеличить энерговооруженность

потока рабочего агента

в процессе

его расширения путем подвода

не только теплоты, но

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

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

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

276

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