Для всех энергетических установок с двухфазным циклом (пар и конденсат) необходимы двухпараметрические таблицы термодина мических свойств рабочего агента, так как в процессах цикла у них имеются фазовые переходы и протекают процессы цикла вблизи кри тической точки рабочего агента. Особенно важно упорядочить теп ловые расчеты аналитическими методами для установок, работающих на водяном паре. Если для энергетических установок с газовым ра бочим агентом в их современном исполнении еще можно довольство ваться уравнением состояния в форме (204), то для водопаровых установок давно назрела необходимость уточнить уравнение состоя ния жидкости и пара, включая и линии фазового перехода в соот ветствующих областях диаграмм состояния этого рабочего агента. Работы в указанном направлении ведутся.
Основным результатом рассмотрения таблиц [22] в § 10 явился вывод о необходимости уравнений состояния для каждой однофазо вой области диаграммы состояния рассматриваемого рабочего агента, для двухфазовой области и для разделяющих их линий фазовых пере ходов. Сейчас уже нельзя удовлетвориться наличием одного уравне ния состояния водяного перегретого пара, которое считалось доста точным до настоящего времени, а приходится искать взаимосвязь параметров Т, р и v и в других фазовых состояниях Н 20.
Так как все такие формы уравнения состояния получают только экспериментальным путем, то необходимо, чтобы точки на линиях, разделяющих соседние области, определялись одним и тем же зна чением параметров Т, р и п, независимо от подхода к ним из той или другой области.
Как ни парадоксально это выглядит на первый взгляд, надо сказать, что основательное знакомство с теплофизикой воды и водя ного пара особенно важно для специалистов по газотурбинным уста новкам. Водяной пар давно применяется в технике как рабочий агент энергетических установок, его теплофизические свойства хорошо изучены и экспериментально и теоретически. Это может слу жить подспорьем при изучении паров других технических жидкостей, применение которых в энергетике возможно. Большое количество закономерностей, установленных для воды и водяного пара, несом ненно, останется принципиально применимым и для других жидко стей и их паров.
Возвращаясь к последним работам по расчетам показателя адиа
баты к, |
остановимся на областях перегретого пара, |
обозначенных |
на рис. |
13 цифрами 2 и 4. В работе |
[25] в качестве уравнения состоя |
ния для рассматриваемого случая |
принято уравнение |
|
|
a = - g r = l + Bp + Cp2 + Dp3. |
(385) |
Для расчета к используется уравнение (382), где значение произ
водной ^ берется из уравнения (385). После подстановки по-
Лучаем выражение для |
показателя а д и а б а т ы перегретого водяного |
пара: |
|
|
|
|
|
ср 1 4- Вр+Ср2+ |
Рр3 |
(386) |
|
cv 2Dp3+ |
Ср2 — 1 |
|
|
Значение изобарной |
теплоемкости |
ср |
берется из таблиц [22], |
а значение изохорной определяется из дифференциального соотно шения
Производные, входящие в это выражение, определяются по урав нению (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 в зависимости от давления
й температуры. Эти данные приводятся в табл. 31. По ним на рис. 48 построены изотермы функциональной зависимости k = k (р, Т). Пользуясь табл. 31 или кривыми на рис. 48, можно при расчетах термодинамических процессов по известным параметрам р и Т при нимать значение показателя адиабаты k [25 и 26].
Рис. 48. Зависимость |
показателя адиабаты k перегретого водяного |
пара |
от давления и температуры. |
Термодинамические процессы, совершаемые рабочим агентом — водой и водяным паром — часто заходят из однофазной области диа граммы состояния в двухфазную или обратно. В таких случаях линия процесса пересекает граничную кривую, разделяющую эти области. Естественно, что фазовые переходы и разные фазовые об ласти состояния рабочего агента требуют специальных исследова ний [82].
Исходной зависимостью в этой области остается формула (380),
но отыскание значения частной производной Должно базиро
ваться на других физических свойствах, отличающихся от свойств
Таблица 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" — полные производные вдоль пограничной линии.
Степень сухости х может быть определена через энтропию:
откуда
?М£-<,-'>+£Ф <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) может быть преобразовано следующим образом:
и уравнение (392) примет вид
d v |
\"дФ |
(395) |
|
|
Таким образом, аддитивность величины |
да s в двухфазной |
области установлена, что существенно упрощает задачу получения простого аналитического уравнения для зависимости k = k (Т, х).
Действительно, подставив (389) и (395) в уравнение (380), получим следующее уравнение для показателя адиабаты в двухфазной области параметров состояния чистого вещества:
|
Ц Т , х ) = — - |
v ' ( 1 — х ) + v "x |
(396) |
|
(£)>-*>+ ( Э Т |
|
|
■ |
Для расчетов по формуле (396) надо иметь значения производ-
ных ^ да удФ |
^ ^ d v у'ДФ |
которые можно совместно написать в виде |
д р Л |
dp ) s |
|
/ 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,