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

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

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

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

Добавлен: 11.04.2024

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

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

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

Для всех энергетических установок с двухфазным циклом (пар и конденсат) необходимы двухпараметрические таблицы термодина­ мических свойств рабочего агента, так как в процессах цикла у них имеются фазовые переходы и протекают процессы цикла вблизи кри­ тической точки рабочего агента. Особенно важно упорядочить теп­ ловые расчеты аналитическими методами для установок, работающих на водяном паре. Если для энергетических установок с газовым ра­ бочим агентом в их современном исполнении еще можно довольство­ ваться уравнением состояния в форме (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