Фазовое равновесие фреона R245fa от 242 к до критической точки Текст научной статьи по специальности «Физика»

УДК 536.423.1

Фазовое равновесие фреона R245fa

от 242 К до критической точки

А. Ж. ГРЕБЕНЬКОВ, В. П. ЦУРБЕЛЕВ, О. В. БЕЛЯЕВА, Т. А. ЗАЯЦ, Б. Д.ТИМОФЕЕВ

Объединенный институт энергетических и ядерных исследований — Сосны Национальной академии наук Беларуси г. Минск, Республика Беларусь

Wide use of R245fa is restrained by the lack of consistent data on thermophysical properties. This paper reviews available different sources, and it is concluded that the most crucial inconsistency exist between measurements of the phase equilibrium parameters. In order to provide accurate data on properties of R245fa, the vapor-liquid equilibrium parameters over a wide temperature range relevant to commercial applications of R245fa were measured. The information obtained, along with other reliable and verified data on the phase equilibrium parameters for R245fa, was fitted with appropriate property correlation models.

Анализ информации из доступных оригинальных источников о термодинамических свойствах R245fa показал, что имеется несогласованность параметров равновесия жидкость — пар (VLE), представленных различными исследователями, особенно в низкотемпературной области, где расхождение между данными по давлению насыщенного пара превышает 2 %. К/,£-данные вблизи критической точки ограничиваются только несколькими измеренными величинами. PVT-ц,анные в однофазных областях, позволяющие рассчитывать величины плотности на линии насыщения, представлены пятью оригинальными работами, которые приемлемо согласуются друг с другом, но имеется широкая область параметров состояния, где все еще отсутствуют прямые измерения.

При создании Стандартной базы данных N 1ST 23 (последней версии REFPROP7.1) использована расширенная модель соответственных состояний (ECS), которая базируется наэксперимен-тальных данных, и, если последние недостаточно проверены, точны и последовательны (что справедливо для имеющегося набора свойств R245fa), параметры модели могут быть также неточны. Главная цель настоящего исследования — представить точные данные о параметрах фазового равновесия R245fa в широком диапазоне температур, важном для коммерческих систем кондиционирования воздуха.

Литературные данные

Давление насыщенных паров

Опубликованы три работы, представляющие экспериментальные данные по давлению насыщенного пара R245fa (табл. 1). В сравниваемых температурных диапазонах только данные, полученные

Nicola [8] и Bobbo [4], согласуются в пределах общей погрешности. Расхождение между этими данными и экспериментальными результатами, полученными Sotani et al. [101, также не превышает 5 кПа, но имеет систематический характер.

Для параметров фазового перехода жидкость — пар в области критической точки R245fa в литературе имеются только величины критической температуры и плотности, рассчитанные на основе измерений критического индекса рефракции и представленные в работе Schmidt et al. |9]. Эти параметры использованы в REFPROP 7.1 [7j, где критическая точка R245fa имеет следующие параметры: Г = 427,2 К, рк = 3640 кПа, рк= 517 кг/м3.

Плотность на линии насыщения

Известны три источника экспериментальных данных о плотности R245fa в жидкой фазе, которые позволяют рассчитать плотность насыщенной жидкости (табл. 1). Информация, представленная в [7], основана на PVT -данных R245fa, полученных на изотермах Defibauch и Moldover [6|. Погрешности измерения температуры, давления и плотности в этой работе составили ±0,01 К, ±0,5 кПа и ±0,5 кг/м^ соответственно, плотность жидкости при давлении насыщения представлена в [6J в виде полинома. Sotani et al. [10] измерили плотность жидкой фазы с погрешностью менее ±0,09 %. Результаты обобщены уравнением состояния типа уравнения Тейта с максимальным отклонением 0,07 %, по которому авторы рассчитали плотность насыщенной жидкости. Bobbo et al. [5] экспериментально определили плотность жидкости также на изотермах и полученные данные аппроксимировали уравнением состояния типа уравнения Тейта с максимальным отклонением 0,21 %.

Обзор данных по параметрам фазового равновесия R245fa

Автор и дата Температурный диапазон, К Диапазон давлений, кПа Количество экспериментальных точек Погрешность*

Давление насыщенных паров Погрешность, кПа

Nicola G.D., 2001 (8] 265 - 352 Давление насыщения 51 0,5

Sotani Т., et al., 1999 110] 293 - 426 34 5

BobboS., et al., 2001 [41 293 - 313 3 1

Плотность жидкой фазы Погрешность, кг/м3

Defibauch et al. 1997 [6] 250 - 370 макс. 6500 1041 < 0,5

Sotani et al. 1999110] 298 - 323 макс. 200000 113 0,09%

Bobboet al. 2003 [5] 283 - 343 макс. 25000 200 0,1 %

Плотность газовой фазы Погрешность, кПа

Nicola 2001 18) 313 - 367 216-610 32 1,5

Sotani et al. 1999 110] 313 -343 151 - 507 98 3

* Взяты из оценок оригинальных работ.

Имеющиеся данные о плотности насыщенной жидкости хорошо согласуются между собой, максимальное отклонение не превышает общую погрешность в сравнимом диапазоне температур и давлений. Исключение составляют результаты [5], полученные на изотерме 343 К, которые выпадают из общего набора данных.

Экспериментальная информация о плотности R245fa в газовой фазе имеется в двух работах (табл. 1). Эти данные хорошо согласуются друг с другом в сравниваемом диапазоне параметров состояния. Nicola [8] и Sotani at al. [10] измерили плотность газа в сравнительно небольшом температурном диапазоне около 40 К. Плотности насыщенного пара [8] и [ 10] отличаются на 2 — 3 %, а данные по плотности пара на пограничной кривой REFPROP V. 7.1 [4] расположены выше данных [8] на 0,8 — 0,9 %.

Имеется одна теоретическая работа, в которой Schmidt et al. [9] представили скейлинговые уравнения для расчета плотности насыщенных жидкости и пара, которые получены на основе измерений рефракции.

Аппаратура и методика эксперимента

Подготовка образца Образцы хладагента R245fa в количестве 20 кг приобретены в фирме Honeywell Fluorine Products Europe B.V. Поставщик указал чистоту вещества 99 %. Тем не менее содержание возможных примесей было проверено посредством хроматографического ана-

лиза двумя методами. Один из них основан на химическом разделении компонентов; другой использует разработанный нами метод определения примесей [2], которые трудно отделимы, но могут иметь разные точки кипения (например, изомеры R245cb или R245ca). Чувствительность аппаратуры составляет около 0,01 массового процента.

Чтобы устранить неконденсируемые газы, каждый образец проходил несколько циклов замораживания жидким азотом с последующей эвакуацией этих примесей. Результаты показали высокую чистоту вещества (не менее 99,9 %).

Аппаратура

Параметры равновесия жидкость — пар R245fa были исследованы с помощью пьезометра постоянного объема в варианте «адиабатического калориметра постоянного объема», который был успешно использован в ряде исследований [1,3, И, 12]. Пьезометр состоит из цилиндрического сосуда из нержавеющей стали с толщиной стенки 1 мм, окруженного двумя изотермическими защитными экранами, и установлен в вакуумной камере. В верхней части имеется клапан для наполнения пьезометра образцом хладагента, а также мембранный разделитель с индикатором «нуля», чувствительность которого составляет 0,005 кПа. Для компенсации давления на мембрану и передачи импульса давления в поршневой манометр или ртутный барометр высокой точности используется азот. Вдоль

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

Автоматическая система поддерживает в пьезометре заданную температуру с точностью ±510-4 К. Вакуумная камера откачивается до давления

1,3-10-7 кПа. Когда температурный диапазон измерений ниже комнатной температуры, вакуумная камера охлаждается жидким азотом. Объем пьезометра при комнатной температуре (0,0729478 м3) был определен с ошибкой 0,06 % калибровкой на дистиллированной воде. Во время калибровки объема был также определен коэффициент изотермического расширения [к= АУ/(Ур) = 3,26-10_3 кПа-1] .

Отдельные серии измерений были выполнены на установке с пьезометром большего объема. Пьезометр оснащен системой отбора части исследуемого

Рис. 1. Схема установки для измерения параметров фазового’ равновесия:

1 — контрольно-измерительная аппаратура;

2 — магнитный сердечник; 3 - гильза с платиновым термометром сопротивления; 4 - пьезометр;

5 — изотермические защитные экраны; 6 — мембранный разделитель; 7 — баллон с азотом; 8 - поршневой манометр (или высокоточный ртутный барометр для измерения низких давлений); 9 — промежуточный (дозирующий) вентиль с подготовленным образцом;

10 — вакуумный насос; 11 — омметрические нагреватели

вещества из паровой и жидкой фаз, что позволяет определять состав вещества или возможные примеси. Экспериментальная ячейка состоит из цилиндрического сосуда длиной около 0,5 м, установленного в термостате, температура которого регулируется с помощью электронной схемы и термометра сопротивления. Изменения температуры в пределах ячейки не превышают ± 0,01 К в течение периода измерения на каждой изотерме. Давление измеряется кварцевым преобразователем давления, откалиброванным с помощью поршневого манометра, а температура — платиновым термометром сопротивления. Погрешности измерения давления и температуры составляют 0,1 % и 0,02 К соответственно.

Процедура РУТ-измерений

РУТ-измерения как в жидкой, так и газовой фазах проводились на квазиизохорах (кривых постоянной загрузки). Образец исследуемого хладагента загружали из промежуточного сосуда в предварительно от-вакуумированный пьезометр через вентиль и капилляр с калиброванным объемом. До и после загрузки (выгрузки) промежуточный сосуд взвешивали и массу образца в пьезометре определяли как разницу масс с ошибкой 0,5 мг. Затем пьезометр разогревали и измеряли температуру и давление. Переход на следующую квазиизохору выполняли посредством удаления определенной части образца из пьезометра в промежуточный сосуд. Погрешность измерений плотности обычно была меньше 0,1 % с доверительной вероятностью 0,95. Температуру измеряли с погрешностью 0,02 К, а ошибка измерения давления составляла около 0,05 %.

Процедура УЬЕ-измерений

При измерении р — Г-параметров кривой насыщения жидкость — пар масса образца К245Га, загруженного в пьезометр в основной серии, примерно соответствовала значению критической плотности и составила 39,6702 г. В критической области такое количество соответствовало плотности 537 г/м3. Это позволило измерить давление пара в температурном диапазоне от 242 К до критической точки в течение одной серии эксперимента.

Каждое отдельное измерение занимало около часа, что гарантировало установившийся режим и условия равновесия. Температурное приращение между измерениями составляло около 5 К в диапазоне от 240 К до 390 К, а затем при приближении к критической точке приращение уменьшалось примерно до 3 К.

При измерении плотности на линии насыщения температуру насыщения определяли в р — Т-диаг-

426----1--------1-------1-------1-------1----

О 10 20 30 40 Ь мин

Рис. 2. Пример термограммы (зависимость температуры Т от времени 0 в окрестности критической точки:

1 - нагрев; 2 — охлаждение

рамме с точностью 0,1 К как точка пересечения двух экстраполированных отрезков квазиизохор в однофазной и двухфазной областях. В окрестности критической точки производная (др/дТ)у при переходе из двухфазной области в однофазную меняется не так сильно. Поэтому температура фазового перехода в этой области параметров определялась методом квазистатических диаграмм [3], который основан на явлении скачкообразного изменения изохорной теплоемкости при фазовом переходе со стороны жидкости. Чтобы реализовать этот метод, пьезометр с образцом хладагента нагревался в адиабатических условиях с помощью нагревателя, расположенного в пьезометре. Скорость нагрева обычно составляла около 1 К/ч. В результате была получена зависимость изменения температуры (температурный ход) от времени (термограмма, рис. 2). Значение температуры фазового перехода соответствует точке пересечения экстраполированных отрезков термограмм в однофазной и двухфазной областях. Затем нагреватель пьезометра отключался, пьезометр медленно охлаждался, и термограммы снимали снова. Усредненное в этих двух режимах значение температуры фазового перехода, полученное с погрешностью 0,1 К, относилось к плотности на линии насыщения для данной изохоры.

Результаты измерений и сравнение данных

Параметры кривой насыщения жидкость - пар С помощью пьезометра постоянного объема получено 61 значение давлений насыщенного пара Я245^а в температурном диапазоне 242...427 К в одной специальной серии вдоль сверхкритической квазиизохоры. Другой набор данных по давлению пара составил 36 экспериментальных точек, полученных для некоторых других квазиизохор в условиях насыщения (табл. 2). В табл. 3 приведены дан -

ные, которые получены с использованием VLE-установки. Это независимая серия, которая состоит из 10 измерений, не дает какую-то новую дополнительную информацию по отношению к результатам основных серий, но подтверждает их точность.

Наши измерения величин давления насыщенных паров хорошо согласуются с данными, представленными Nicola [8] и Bobbo etal. [4]. Максимальное отклонение в сравниваемом температурном интервале не превышает 0,3 %. В области низких температур отклонение находится в пределах 0,13 кПа. В температурном диапазоне от 330 до 350 К, где наибольшее число измерений выполнил Nicola, отклонения достигают 1,9 кПа. Согласие наших данных с экспериментальными результатами, представленными Sotaní et al. [10], менее очевидно. В низкотемпературной области данные отклоняются на величину до 4 кПа, что соответствует 2,5 % от измеренного давления. В высокотемпературном интервале отклонение достигает 17 кПа. В интервале температур до 340 К наши данные согласуются с корреляцией, использованной в REFPROP 7.1 [7], но с повышением температуры расхождение увеличивается, так что в критической области достигает 13 кПа.

Используя описанные выше метод квазистатических термограмм и метод аппроксимации в р — Т-диаграмме отрезков кривых опытных квазиизохор, мы определили 32 значения плотности насыщенной жидкости и 27 значений плотности насыщенного пара (табл. 4). Наши экспериментальные данные согласуются с результатами, представленными в литературе, в температурном диапазоне, где существует экспериментальная информация о плотности R245fa в однофазных состояниях. Для других температур и, что особенно важно, для области вблизи критической точки наши данные являются единственной экспериментальной информацией о плотности насыщенных фаз. Тем не менее наши измерения в области критической точки хорошо согласуются (среднеквадратичное отклонение 0,2 % для насыщенной жидкости и 3 % для насыщенного пара) с расчетами, выполненными в рамках масштабной теории [9]. При температурах до 410 К имеется удовлетворительное совпадение с моделью, использованной в REFPROP 7.1 [7].

Корреляция и анализ данных

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

Экспериментальные значения давления насыщенных паров (метод пьезометра постоянного объема)

Температура X К(1) 241,826 247,923 250,316 250,361 251,656 251,898 257,261 258,715 261,823 262,003

Давление насыщеных паров рз, кПа (II) 9,0502 14,030 16,426 16,545 17,597 18,117 23,891 25,974 30,111 30,396

Температура Т, К (I) 264,344 270,504 276,364 277,230 279,360 284,179 284,535 285,946 291,375 295,943

Давление насыщеных паров ра, кПа (II) 34,619 46,766 61,386 63,466 69,683 85,622 86,998 91,945 114,46 137,63

Температура Т, К (I) 297,109 299,619 300,450 301,650 303,568 304,533 305,046 306,988 308,245 309,997

Давление насыщеных паров р8, кПа (II) 142,66 157,86 160,70 168,55 181,11 187,85 190,92 202,78 212,89 224,86

Температура % К (I) 310,659 314,607 316,308 320,399 321,876 326,884 328,464 333,450 338,968 341,757

Давление насыщеных паров р,, кПа (II) 230,84 261,56 277,69 314,50 331,34 385,55 403,49 466,52 544,23 585,80

Температура Т, К (I) 345,808 352,760 357,924 364,828 367,908 371,188 374,175 375,553 377,338 380,023

Давление насыщеных паров рз, кПа (II) 653,69 780,66 886,57 1044,6 1121,6 1208,2 1291,2 1329,5 1383,7 1468,2

Температура Т, К (I) 380,362 380,790 381,349 382,194 384,166 385,589 387,09 387,32 390,034 393,485

Давление насыщеных паров р5, кПа (II) 1477,8 1490,4 1509,1 1536,9 1600,4 1650,3 1701,3 1708,7 1806,6 1938,4

Температура X К (I) 393,507 394,604 394,615 396,911 400,45 401,412 403,318 406,063 406,361 406,474

Давление насыщеных паров р4, кПа (II) 1937,9 1983,1 1982,1 2073,7 2222,2 2265,2 2348,2 2476,1 2486,9 2493,9

Температура Т, К (I) 408,836 409,581 410,645 411,441 412,549 413,675 414,264 416,128 416,577 418,646

Давление насыщеных паров р4, кПа (II) 2606,0 2643,6 2695,7 2737,4 2794,1 2852 2883,1 2984,8 3009,9 3123,6

Температура Т, К (I) 419,461 419,685 420,334 421,189 421,705 421,894 422,789 423,599 423,889 423,948

Давление насыщеных паров р4, кПа (II) 3171,5 3184,1 3223,1 3271,8 3301,6 3314,2 3369,9 3419,4 3435,3 3439,7

Температура Т, К (I) 424,403 424,463 424,641 424,886 425,477 426,096 427,101 ' - - -

Давление насыщеных паров р8, кПа (II) 3467,4 3473,3 3484,8 3498,8 3538,6 3576,3 3643,7 - - -

термограмм вдоль околокритических квазиизохор. Критическая плотность рк определена путем наи-Это значение равно 427,158 К (ІТ8-90) с погрет- лучшего подбора параметров масштабного уравне-ностыо 0,02 К, которая в основном определяется ния [3]: точностью измерения температуры платиновым термометром сопротивления.

*0-7,7

Экспериментальные значения давления насыщенных паров (УЬЕ-установка)

2 = пП-т^8&1(р -л)*В* +В1^-Тг) +..., (3)

Таблица 3

Температура Т, К (I) 293,11 293,25 298,06 303,24 313,12 323,04 333,22 343,56 353,40 354,04

Давление насыщеных паров р5, кПа (II) 123,43 123,05 147,85 178,08 250,37 343,40 463,26 618,95 794,79 803,55

Таблица 4

Экспериментальные значения плотности насыщенной жидкости (р') и пара (р")

Т, К

255,998

263,500

277,251

286,471

314,911

353,176

371,352

394,518

402,760

413,026

420,388

425,363

426,880

427,076

427,151

427,153

Р . кг/м3

1445,5

1427,3

1390,3

1367,6

1290,4

1170,6

1101,6

991,79

936,28

859,87

777,71

686,08

597,90

576,19

539,31

539,27

Т, К

255,998

263,500

263,500

277,251

277,251

286,471

286,471

314,911

314,911

353,176

353,176

371,352

371,352

427,155

427,157

427,158

Р,

кг/м5

1445,2

1427,3

1427,5

1390,3

1392,5

1367,6

1368,0

1290,4

1289,3

1170,6

1171,3

1101,6

1099,8

536,06

528,06

520,42

Т, К

266,593

293,601

305,138

323,261

329ДОЗ

337,136

346,230

359,254

382,404

386,805

409,417

419,712

423,197

426,409

427,046

427,116

Р , кг/м5

2,6070

7,2240

10,806

19,286

23,702

29,244

37,199

52,794

92,699

100,50

181,19

262,26

309,38

402,18

456,14

470,99

Т, К

427,149

266,593

305,138

323,261

329,203

337,136

346,230

359,254

362,810

382,404

386,805

Р , кг/м3

490,37

2,5650

11,071

19,432

23,125

28,972

37,180

52,294

57,213

92,058

102,20

которое может быть получено из известного степенного разложения:

/"’-ft

р.

= В,(\-Тгу +В&-Т,)'-' + в2(1-гу+...Д4)

где Тг-Т/Тк - приведенная температура;

В0 - критическая амплитуда;

В[У В2, и т.д. - другие амплитуды масштабного уравнения, р = 0,34, а = 0,12 и А = 0,49.

Функция Z=J{\ - Т)°'\ построенная на базе наших измерений плотностей насыщения (рис. 3), имеет две ветви, соответствующие насыщенной фазе, которые должны пересекаться друг с другом при Z=B0 в точке (/ -Т)°'5= 0. Если критическая плотность определена неправильно, жидкостная и паровая ветви Z-фун-кции не пересекаются в этой точке (см. рис. 3, би в). Наилучшим образом подобранное значение критической плотности составляет рк =516,68 кг/м3 с оцененной погрешностью 0,5 кг/м3.

В температурном диапазоне от 240 К до критической точки результаты настоящего исследования по давлению насыщенных паров ps, а также данные, представленные Nicola [8] и Bobbo et al. [4], аппроксимированы выражением типа уравнения Вагнера:

inA=^¿4(i-rrr,

Pm *r 1

(5)

Рис. 3 Масштабная функция Z-f(l- Т)°-$, построенная на основе измерений: а - рк = 516,68 кг/м3; б - рх= 518 кг/м3; в — рк =515 кг/м3

где а, и А{ — параметры регрессии (табл. 5).

Критическое давление рк было также получено как параметр регрессии и по оценкам составляет 3647 кПа с погрешностью 0,05 %. Стандартное отклонение всех данных, включенных в это обобщение, составляет 0,67 кПа. На рис. 4 показаны отклонения имеющихся ps от уравнения (5).

В температурном диапазоне от 250 К до критической точки результаты нашего исследования плотности жидкости и пара в насыщенном состоянии аппроксимированы следующими выражениями:

Р ~Рк

: А>(1 - TrfM+- Trf.

(6)

Параметры регрессии уравнения (5)

i 1 2 3 4 5 6 7

ai 1 1,5 2,5 3 4 5 6

Л -7,87353351 1,96083166 0,20906769 -14,9888978 47,51861087 -82,92485701 43,09954343

V 1 1 • о

Л Л V -з □ —4

\ Лд Д -5

□ ь д д

Щ 1 . . 1 1 |

______________i___________________i___________________i___________________i---------

250 300 350 400 Т, К

Рис. 4. Отклонения 5 величин давления насыщения ps от уравнения (5):

1 — наше исследование; 2 — наше исследование на VJ.E-установке; 3 — Nicola [8]; 4 - ВоЪЪо et al. [4]; 5 - Sotani et al. [10]; 6 - REFPROP 7.1 [7]

-£^=-в0(1-г,Г34+£д"(1-гУ > (7)

Рк ¿=1

где В0 - критическая амплитуда;

Ь'п Ь”у В/ и В" - другие параметры регрессии (табл. 6).

Стандартное отклонение данных по плотности, включенных в обобщение, составляет 0,27 % для плотности насыщенной жидкости и 0,99 % для плотности насыщенного пара. На рис. 5 и 6 пока-

Таблица 6

Параметры регрессии уравнений (6) и (7)

Рис. 5. Отклонения 6 плотности насыщенной жидкости р'от полинома (6):

1 — наше исследование; 2 - Defibaugh [6]; 3 - Bobbo [4]; 4 - Sotani [10]; 5 - REFPROP 7.1 [7]

250 300 350 400 T, К

Рис. 6. Отклонения 8 плотности насыщенного пара р" от полинома (7):

1 - наше исследование; 2 - Nicola [8]; 3 - Sotani [10]; 4- REFPROP 7.1 [7]

i ь: i в; V в;

0 0,34 2,041145 0,34 2,041145

1 1 -0,24512768 1 1,09378121

2 1,5 8,99699246 2 1,57779004

3 2 -25,45105251 3 -7,55413064

4 2,5 28,16685893 4 18,26381319

5 3 -10,21987307 5 -17,18161668

заны отклонения имеющихся данных по насыщенным плотностям от уравнений (6) и (7).

Данная работа выполнена в рамках проекта № 1256-RP корпорации «Американское общество инженеров теплоснабжения, холодильных устройств и кондиционирования воздуха» (ASHRAE), и проекта N° 05К-001 Белорусского республиканского фонда фундаментальных исследований.

Список литературы

1. Вер минская А. Б., 1'ребеньков А. Ж., Цурбелев В. II. Исследование термодинамических свойств диоксида азота в широкой области критического состояния // Теплофизические свойства веществ и материалов, 1989. №. 27.-М.: ГССД.

2. ТимофеевВ. Д., Беляева О. В., НехайчикН.А., Николаев В А. Газохроматографический анализ азотсодержащей газовой смеси // Заводская лаборатория. 1981. Т. 47, № 7.

3. Цурбелев В. //., 1'ребеньков А. Ж. Параметры фазового равновесия тетраоксида азота и их асимптотическое поведение в критической области: II. Плотность// Известия Академии наук Беларуси. Сер. физ.- энерг. наук. 1985, №. 1.

4. Bobbo S., Fedele L., Scattolini М., Camporese R. Isothermal VLE measurements for the binary mixtures H FC-134a+ H FC-245fa and HC-600a + HFC-245fa// Fluid Phase Equilibria . 2001. V. 185.

5. Bobbo S., Fedele I., Scattolini M. Camporese R. Compressed Liquid Densities and Saturated Liquid Densities OF HFC-245fa International congress of refrigeration 2003. Washington. DC (ICR402).

6. Defibaugh I). R., and Michael R. Moldover. Compres-sed and Saturated Liquid Densities for 18 Halogenated Organic Compounds.//Journal of Chemical and Engineering Data. 1997. Vol. 42. N 1.

7. McLinden. M. ()., HuberM. L., lemmon. E. W, andPeskin. A. P. G. 2003. NIST Standard Reference Database 23. NIST Thermodynamic Properties of Refrigerants and Refrigerants Mixtures Database (REFPROP). Version 7.1 (Gaithersburg: National Institute of Standard and Technology).

8. Nicola G. I). Vapor Pressure and Gas Phase P-V-T Data for 1.1.1.3.3-pentafluoropropane (R245fa)//JournaI of Chemical and Engineering Data. 2001. Vol. 46, N 6.

9. Schmidt J. W., Carrio-Nava E., MoldoverM. R. Partially halogenated hydrocarbons CHFCL-CF3, CF3- CH3 CF3- CHF- CHF2, CF3- CH2 - CF3, CHFr CH2-CHF2, CF3- CH2 — CHF2; CF3-0- CHF2: critical temperature, refractive indices, surface tension and estimates of liquid, vapor, and critical densities //Fluid Phase Equilibria. 1996. V. 122.

10. Sotani Kubota H. Vapor pressures and PVT properties of. I 1.1.1.3.3-pentafluoropropane HFC-245fa// Fluid Phase Equilibria. 1999. V. 161.

11. Tsurbelev VP, GrebenkovAJ., andKlepatsky P.M. 1998. Volumetric Behavior of Refrigerant Mixture R32/R125 // Proc. of the Fifth Asian Thermophysical Properties Conference, vol. 2, M.S. Kim and S.T. Ro, ed., Seoul National University, Seoul, Korea.

12. Tsurbelev V. P., GrebenkovAJ., and Klepatsky P.M. 2000. P-V-T and phase equilibrium data of a R32-R125 refrigerant mixture // Proc. of the 14th Symposium on Thermophysical Properties, Boulder, Colorado, USA, June 25-30, 2000. University of Colorado at Boulder. Collection of Preprints.