Экспериментально-расчетное исследование термодинамических свойств октафторпропана и декафторбутана тема автореферата и диссертации по физике, 01.04.14 ВАК РФ

ООЗ^ычаг-" На правах рукописи

Кузнецов Кирилл Игоревич

Экспериментально-расчетное исследование термодинамических свойств октафторпропана и декафторбутана

Специальность - 01.04.14 «Теплофизика и теоретическая теплотехника»

2 6 НОЯ 2009

АВТОРЕФЕРАТ диссертации па соискание ученой степени кандидата технических наук

г. Москва 2009 год

003484920

Работа выполнена на кафедре теоретических основ теплотехники им. М.П. Вукаловича Московского энергетического института (технического университета)

Научный руководитель — кандидат технических наук,

старший научный сотрудник Сухих Андрей Анатольевич

Официальные оппоненты — доктор технических наук, профессор

Махров Виктор Владимирович

— кандидат технических наук, доцент Ларкин Дмитрий Константинович

Ведущая организация - Федеральное государственное унитарное

предприятие «Российский научно-технический центр информации по стандартизации, метрологии и оценке соответствия»

Защита состоится «04» декабря 2009 г. в 11:30 часов на заседании диссертационного совета Д 212.157.04 при Московском энергетическом институте (техническом университете) по адресу: г. Москва, Красноказарменная ул., д. 17., корп. Т, ауд. Т-206

С диссертацией можно ознакомиться в библиотеке МЭИ (ТУ),

Отзыв на автореферат, заверенный печатью учреждения, просим направлять по адресу: 111250, Москва, Красноказарменная ул., 14, Ученый совет МЭИ (ТУ).

Автореферат разослан « -3 » И0 Я2009 г.

Ученый секретарь

диссертационного совета Д 212.157.04, к.ф.-м.н. доцент

В.И. Мика

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы

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

Рабочие вещества фторуглеродного класса инертны, не горючи, широко используются в качестве низко- и среднетемпературных теплоносителей в атомных, химических и космических технологиях. Например, смесь перфторпропана (С3Р8) с элегазом (БРб) применяется в качестве хладагента (11510а) в холодильных и теплонасосных установках. Другим примером является использование фторуглеродов в климатической системе космической станции «Мир» и т.д. Исследования последних лет показали, что рабочие вещества фторуглеродного класса имеют термическую стойкость, достаточную для реализации и высокотемпературного термодинамического цикла.

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

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

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

Для ряда рабочих тел фторуглеродного состава, отвечающих требованиям, предъявляемым к современным энергоустановкам, таким как энергетическая эффективность, низкая либо нулевая токсичность, совместимость с конструкционными материалами и отвечающими экологическим

требованиям Киотского и Монреальского протоколов к настоящему времени данные, необходимые для расчетов циклов, уже существуют. Например, для тетрафторэтана (CF4), циклофторбутана (C4F8) были выполнены оценочные расчеты термодинамических циклов. При температуре на входе в турбину 530 °С в сверхкритическом цикле с регенерацией теплоты от температуры выхода из питающего насоса с использованием циклофторбутана (С4Р8) термический КПД составил 46%. Высокая термодинамическая эффективность применения в качестве рабочих тел фторуглеродного состава (октафторпро-пана C3F8, циклофторбутана C4FS, декафторбутана C4F|0) для энергетических установок с бинарными циклами впервые была отмечена в работах Гохштейна Д.П. и соавторов ещё в 60-ые годы прошлого века. Однако расчеты циклов с приемлемой точностью могут быть выполнены только по результатам настоящей работы.

Цель работы

Модернизация экспериментальной установки для измерения плотности фторуглеродов в широком диапазоне параметров состояния, охватывающем все расчетные точки проектируемых циклов энергетических установок. Экспериментальное исследование плотности октафторпропана (CjFg) и декафторбутана (C4Fio), как наиболее перспективных рабочих тел в широком диапазоне температур и давлений. Получение термического уравнения состояния каждого из исследуемых веществ на основе полученных экспериментальных данных и экспериментальных данных других авторов из ранее опубликованных источников. Расчет термодинамических свойств исследуемых веществ по полученным уравнениям состояния.

Основные результаты работы:

1) Проведена модернизация экспериментальной установки для измерения плотности чистых газообразных и жидких веществ в интервале температур 300.. .800 К и давлений 1... 10 МПа.

2) Разработана методика исследований с электронным протоколированием результатов измерений в процессе опытов.

3) Получены новые экспериментальные данные по плотности C3F8 и C4Fio, в диапазоне давлений 1...10 МПа и температур 393,15 К. ..773,15 К объемом в 158 экспериментальных точек.

4) Построены уравнения состояния октафторпропана (C3F8) и декафторбутана (C4F|0) на основе собственных экспериментальных данных и данных других авторов, доступных на момент проведения

5) рабочшганы таблицы термодинамических свойств октафторпропана (C3FS) и декафторбутана (C4F10), охватывающие интервал давлений 0,1... 10 МПа и температур 293Д5...823Д5 К в газовой области.

6) Проведен анализ полученных результатов и сравнение их с имеющимися данными других авторов, в том числе и по теплоемкости С3Р8.

Научная новизна

1) Расширен диапазон экспериментальных данных по термическим свойствам C3F8 и C4F10 при повышенных температурах в газовой фазе, так и при средних температурах в области жидкости.

2) Данные о термодинамических свойствах декафторбутана (C4FI0) в газовой фазе и в области сверхкритических состояний получены впервые.

3) Разработаны экспериментально обоснованные многоконстантные уравнения состояния вириального типа для расчета термодинамических свойств (C3FS) и (C4F10) в диапазоне, необходимом для проектирования энергетических циклов на рабочих телах фторуглеродного состава.

Практическая ценность работы.

Диссертационная работа выполнена в соответствие с хозяйственным договором № 280-07/2216070-МЭИ от 02.06.07 между ГОУВПО МЭИ (ТУ) и НИКИЭТ им H.A. Доллежаля «Расчетно-экспериментальное исследование теплофизических и теплотехнических свойств фторуглеродов с целью оценки их использования качестве рабочего тела 2-го контура реакторной установки на быстрых нейтронах».

Работа поддержана грантом РФФИ № 09-08-00741-а «Исследование термодинамических свойств рабочих тел фторуглеродного состава для теплосиловых циклов».

На защиту выносятся:

Результаты экспериментальных исследований плотности октафторпропа-на (C3F8; Тс = 345,02 К, рс = 2,670 МПа) и декафторбутана (C4F,0; Тс = 386,35 К, рс = 2,324 МПа).

Таблицы термодинамических свойствах октафторпропана (C;,F8) и декафторбутана (C4F10) в широкой области параметров состояния.

Публикации

Основные положения диссертационной работы представлены в публикациях [1 - 4].

Структура и объем работы

Диссертация общим объемом 116 страниц состоит из 8 разделов, включающих в себя введение, главы основной части, заключение и таблицы термодинамических свойств рабочих тел фторуглеродного состава. Диссертация включает в себя 32 иллюстрации, 16 таблиц, 63 формулы, ссылки на 52 литературных источника.

Содержание работы

Во введении (глава I) обосновывается актуальность темы диссертации. Сформулирована цель и основные задачи работы. Выбраны объекты исследований.

Вторая глава посвящена обзору опубликованных литературных данных о теплофизических свойствах веществ, подлежащих исследованию в настоящей работе. По результатам изучения опубликованных данных сделан важный вывод о невозможности расчетов паросиловых циклов энергетических установок, использующих октафторпропан (С3Р8) и де-кафторбуган (С4Рю) без пополнения базы экспериментальных данных по плотности в газовой области при повышенных температурах для СзР8, а для С4Рю и еще и в жидкой фазе, на кривых кипения и конденсации. Этот вывод особенно важен, так как данных о термических свойствах декаф-торбутана очень мало, что подтверждается представленной в диссертаци-ии диаграммой исследованных областей в р-Т координатах. Диаграмма

проводились исследования, представленные соответствующими линиями. Данные, полученные автором, представлены значками (в), другие данные МЭИ, полученные Закопыриным М.А. представлены значком (А), остальные данные представлены значком (•). Только на такой базе собственных опытных данных и привлеченных литературных становится возможным построение достоверных широкодкапазонных уравнений состояния.

Третья глава посвящена обзору основных методов исследования плотности веществ. В главе выявлены достоинства и недостатки различных методов и способов их реализации. Выбран метод пьезометра постоянного объема, позволяющий провести измерения с привлечением метрологического оборудования самого высокого уровня научно-исследовательской лаборатории кафедры ТОТ МЭИ.

В четвертой главе описывается экспериментальная установка для проведения измерений плотности веществ в широком диапазоне температур и давлений. Схема экспериментальной установки представлена на рисунке 2.

Вентиляция

Система измерения и регулирования температуры

ш

вода

Рисунок 2. Схема экспериментальной установки

I - пьезометр постоянного объема; 2 - жидкостный термостат; 3 -мембранный дифференциальный манометр; 4 - капилляр; 5 - система выпуска исследуемого вещества; 6 - аналитические весы; 7 - вторичный прибор мембранного дифференциального манометра; 8 - мешалка жидкостного термостата; 9 - двигатель мешалки; 10 - вакуумный насос;

II - горизонтальный разделитель масляных систем; 12 - грузо-поршневые манометры МП-6, МП-60, МП-600.

В экспериментальной установке для проведения основных измерений использовались измерительные приборы, имеющие прецизионную точность. Для измерения давления применялись грузо-поршневые манометры МП-6 класса 0,02, МП-60 класса 0,02 и МП-600 класса 0,02. Для измерения температуры опыта использовался платиновый термометр сопротивления, имеющий погрешность 0,01°С в комплекте с полуавтоматическим потенциометром Р-3 63/1 класса 0,001.

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

При подготовке к проведению эксперимента необходимо установить в систему выпуска вещества предварительно взвешенный на аналитических весах съемный баллончик. Линия, соединяющая съемный баллон с мембранным разделителем, предварительно вакуумируется перед каждым выпуском. Под баллон, присоединенный к выпускной линии, подводится сосуд Дыоара. который далее заполняется жидким азотом. При открытии выпускного вентиля на дифманометре происходит вымораживание выпущенного исследуемого вещества в съемный баллон. После закрытия выпускного вентиля - вымораживание продолжается, и давление в линии становится равным давлению СуиЛНмацйИ Прй ТсмпсраТурс КИПсНИЯ ЖИДКОГО аЗОТа. ТаКОС ДаБЛСНИС вполне соизмеримо с форвакуумом при откачке. Поэтому остаточная масса вещества в линии пренебрежимо мала. После отсечения линии закрытием вентиля на съемном баллончике последний отсоединяется, нагревается до комнатной температуры и взвешивается на аналитических весах. Изменение массьт вещества в съемном баллоне в результате описанных выше операций -является массой выпущенного вещества между двумя стационарными состояниями. Последний выпуск производится до полного равновесия между давлением в пьезометре и давлением в съемном баллоне при температуре жидкого азота. Так как данных о давлении сублимации для исследуемых веществ нет, то автор считает допустимым считать, что давление в системе ниже давления в тройной точке потому что температура тройной точки выше температуры кипения жидкого азота при атмосферном давлении. Данные о давлении в тройной точке найдены в литературе только для СзР8, оно равно 20 Па. Остаточная масса в пьезометре, найденная по идеально-газовому уравнению для минимальной изотермы (Т= 373,15 К) составила 0,4 мг, что пренебрежимо мало и может не учитываться в опыте. Масса в пьезометре в каждом состоянии может быть определена, как сумма масс выпусков от данного состояния до конечного. Плотность вещества в каждом состоянии может быть найдена, как

где V - объем пьезометра при данных условиях опыта, т - масса исследуемого вещества в пьезометре в данном состоянии. Объем пьезометра, приведенный к 20°С, разгруженный от давления составляет 204,53±0,01 см3. При проведении опытов учитывалась масса вещества, находящаяся в верхней полости мембранного разделителя, а также масса вещества в капилляре, соединяющем пьезометр и мембранный разделитель. Для нахождения этих масс в конструкции экспериментальной установки предусмотрено измерение температуры в полости дифманометра и в каждой точке по длине капилляра. Расчет этих масс проводился с ис-

пользованием таблиц ГССД. Данные о материале пьезометра для внесения необходимых поправок к его объему в зависимости от температуры и давления опыта были взяты из литературных источников для стандартного состава использованного сплава ЭИ-617.

Для нахождения объема пьезометра и балластного объема в полости мембранного разделителя были проведены специальные опыты с использованием хорошо изученных веществ (воды и углекислоты). Методика их проведения описана в пятой главе. Кроме того, в пятой главе описана методика проведения опытов по измерению плотности исследуемых веществ октаф-торпропана (СзР8) и декафторбутана (С4Р;о) и аттестационных опытов по определению плотности хорошо изученных веществ (диоксида углерода и шес-тифтористой серы). В главе описана методика расчета экспериментальных данных по результатам проведенных измерений, анализ погрешностей измерений и полученных экспериментальных значений плотности, сравнение результатов аттестационных измерений с данными источников, имеющих статуе стандартных. Сравнение результатов проведенных аттестационных опытов со стандартными справочными данными приведены в таблице 1. Таблица 1 - Результаты аттестационных опытов

т, к р, МПа рЯ1СТ1, кг/м3 Р лит, кг/м3 5,%

Шестифтористая сера

448,24 9,3840 460,055099 459,28 0,17

448,24 7,9433 380,573313 379,91 0,17

448,25 5,9949 275,613614 274,95 0,24

448,25 3,9799 173,762923 173,45 0,18

448,27 1,1762 47,5492757 47,55 -0,0015

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

333,10 | 5,9193 122,86 122,79 0,057

333,12 | 9,8767 283,14 282,78 0,13

По результатам анализа погрешностей измерений установлено, что в области давлений 2... 10 МПа во всем температурном диапазоне измерений погрешности без учета вероятных ошибок отнесения не превышают 0,25 %. С учетом вероятных ошибок отнесения погрешность не превышает 0,27%. На высоких изотермах (выше 523К), при давлениях ниже 2 МПа, погрешность возрастает до 0,31%. В надкритической области, где величина производной ^^ф ] существенно возрастает, полная погрешность

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

главы. В настоящем автореферате таблицы экспериментальных данных приведены в приложениях А1, А2, АЗ.

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

Г в!

¡=1 1=о

где,

г= pv / RT ; е> = р/рк

приведенная плотность; т

Т/ 71,

■ приве-

денная температура. Для получения коэффициентов Ь^ уравнения использован метод наименьших квадратов на основе минимизации функционала (Б).

k=l

г-1

i=l j=0

ю т

-J

(3)

J

где весовая функция W = l/(8kz')2, 5К - предельная относительная погрешность экспериментальных значений сжимаемости 7?. Для описания использовалась компьютерная программа написанная на языке FORTRAN в пакете компьютерных программ FORTRAN POWER STATION 4.0. Для массива экспериментальных Р,р,Т-данных величина среднеквадра-тического отклонения аппроксимации вычислялась по соотношению:

02=1ОО

±{l-zp/z3J/(n-m)

(4)

где г, 7/ - соответственно экспериментальное и расчетное значение сжимаемости в каждой точке;

п - число экспериментальных точек; т - число эмпирических констант уравнения. По известным, описанным в различных литературных источниках соотношениям, были рассчитаны термодинамические свойства С^Чо и С3Р8. Данные о идеально-газовой энтальпии, энтропии и теплоемкости были взяты из таблиц

идеально-газовых функций, рассчитанных в Техасском университете. Для построения уравнения состояния, способного устойчиво описывать термодинамические свойства исследуемых веществ с возможностью монотонной экстраполяции на параметры состояния до 823,15 К и 12 МПа были получены 15 эквивалентных уравнений для С4Р|о и 10 эквивалентных уравнений для СзР8, имеющие различные по величине и даже по знаку коэффициенты, при этом описывающие массив экспериментальных данных с одинаковым (либо близким) среднеквадратическим отклонением. По усредненным коэффициентам

было построено уравнение для каждого из рассматриваемых фторуглеродов.

Проведенное сравнение отклонений рассчитанных термодинамических свойств декафторбутана по каждому из эквивалентных уравнений со свойствами, рассчитанными по усредненному уравнению показало их хорошее взаимное согласование. Для С4Рю среднеквадратические отклонения энтропии не превысили 0,03%, максимальные не превысили 0,61 %; среднеквадратические отклонения энтальпии не превысили 0,1%, максимальные не превысили 0,62%; среднеквадратические отклонения изобарной теплоемкости не превысили 0,06%, максимальные не превысили 2,41%. Усредненные коэффициенты уравнения (2), описывающего декафторбутан представлены в таблице 2.

Таблица 2 - Коэффициенты усредненного уравнения С4Р)0

Ью .9576033100+01 Ь32 -.1849829060+03

Ь„ -.8928028750+01 Ьзз .4977634090+03

Ъп -.6096927620+02 Ь34 -.2312576680+03

Ь,з .1368991180+03 Ь4о -.1260729980+03

Ь)4 -.1053456220+03 Ь4; .1988294520+03

Ь,5 .2758479650+02 Ь42 .2403330890+03

Ьго -.7430784640+02 Ь43 -.5776038160+03

Ьц .2231880920+03 ъ44 .2639042070+03

Ь22 -.2069015080+03 Ь5о .8113591210+02

Ь23 .4817723880+02 Ь51 -.2331574570+03

Ь24 .2361993250+01 Ь52 .1875471330+03

Ьц .7648560750+01 ь53 .1467010330+01

Ьзо .1185352450+03 ъ54 -.3683574460+02

Ь3! -.19931142Ю+03

Для С3Р8 специальное исследование взаимных отклонений в эквивалентных уравнениях не проводилось, так как существует возможность сравнений с данными других авторов по изобарной и изохорной теплоемкостях в наиболее трудно измеряемых и трудно описываемых областях. Усредненные коэффициенты уравнения (2), описывающего октафторпропан представлены в таблице 3.

Таблица 3 - Коэффициенты уравнения состояния октофторпропана

Ью 10,0991666 Ь32 376,6919530

Ьи -69,1750574 Ь33 -359,6544650

Ь,2 188,5019770 Ьз4 127,1500710

Ь,з -252,6102350 ь40 -14,7356347

Ь|4 162,5383900 Ь41 101,6199720

Ь15 -40,4872052 ь42 -238,3964660

Ь20 -23,1360654 Ь4з 234,8300040

Ь21 157,2898870 Ъа4 -84,1432884

Ь22 -403,9661080 Ь5о -8,8830671

Ьгз 487,8757060 Ь5, 29,5663888

Ь24 -272,6943450 ь52 -30,9287941

ь25 54,6259367 Ь5з 8,0789663

Ьзо 28,99.16996 ь54 2,3821803

Ьз1 -172,1530880 |

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

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

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

Заключение

1. Модернизирована экспериментальная установка для исследования плотности веществ в жидком и газообразном агрегатном состоянии в диапазоне температур 25...500 °С и давлений от атмосферного до 10 МПа.

2. Проведено экспериментальное измерение плотности декафторбутана ДО,о) в широком диапазоне параметров состояния с погрешностью измерений без вероятных ошибок отнесения не хуже 0,24 % во всей области исследований. Погрешности измерений с учетом вероятных ошибок отнесения при давлениях ниже 2 МПа не превышают 0,35 %, в надкритической области (Т=293,15...413,15К, р=2...3,5 МПа) не превышают 0,45 %. В остальной, рабочей области параметров состояния, используемых при расчетах энергетических установок, полная погрешность измерений не превышает 0,3 %. Результаты измерений представлены в таблицах А2, АЗ приложения.

3. Проведено экспериментальное измерение плотности октафторпропана (СзРз) в широком диапазоне параметров состояния с погрешностью измерений без вероятных ошибок отнесения не хуже 0,24 % во всей области исследований. Погрешности измерений с учетом вероятных ошибок отнесения при давлениях ниже 2 МПа не превышают 0,35 %, в надкритической области (1=248,15...353,15К, р=2,8...3,8 МПа) не превышают 0,45 %. В остальной, рабочей области параметров состояния, используемых при расчетах энергетических установок, полная погрешность измерений не превышает 0,3 %. Результаты измерений представлены в таблице А1 приложения.

4. Общее число опытных точек, включая измерения плотности в жидкой фазе, составило 162.

5. Проведено расчетно-теоретическое исследование основных термодинамических свойств декафторбутана (СдРш) и октафторпропана (С3Р8) с использованием экспериментальных данных настоящей работы и экспериментальных данных других авторов. На этой базе разработаны многоконстантные уравнения состояния вириального типа, достоверно описывающие свойства исследуемых веществ в широком диапазоне параметров состояния (р = 0Д...10 МПа, Т = 293,15...823,15К). Среднее квадратическое отклонение такого описания составило 0,24 % для октафторпропана и 0,27 % - для декафторбутана.

6. Построены таблицы термодинамических свойств исследуемых веществ для обеспечения инженерных расчетов как циклов, так и основных аппаратов энергетических установок на фторуглеродных рабочих телах.

Приложения:

Таблица А1 - Экспериментальные значения плотности октафторпропана в однофазной области

453,29 2,9896 171,06 673,69 4,8333 162,60

453,17 1,9685 106,92 673,45 3,8909 131,32

493,01 9,8694 543,15 673,37 2,9140 98,164

493,03 7,9292 438,23 673,04 2,0339 68,363

493,07 5,9165 319,35 773,23 10,0214 273,17

493,05 4,9494 261,58 773,16 8,8729 245,63

493,12 3,9537 203,62 773,22 7,5579 211,98

AQl 17 • »» ч» j ■ • 3,0261 151,58 772.96 5,9207 168.22

493,15 1,9635 95,295 773,15 4,9167 140,92

493,26 1,0066 47,506 773,22 3,9952 115,08

532,94 9,8349 462,04 773,24 2,9948 86,83

533,15 7,8540 370,92 773,23 1,1382 33,190 I

533,15 5,8966 275,71 !

Таблица А2 - Экспериментальные значения плотности декафторбутана в однофазной области

473,19 2,9088 221,60 672,99 2,0774 86,864

| 473,15 1,9570 137,35 672,99 2,0766 87,210

473,14 1,0762 70,101 672,93 1,1425 47,499

523,15 9,8815 685,40 772,86 9,8859 337,48

523,14 8,9354 634,92 772,96 7,9112 277,82

523,20 7,9117 556,47 773,05 6,0063 212,02

523,26 6,0402 417,69 772,96 5,0418 177,771

523,26 5.0087 340.56 773,04 4,0276 141,72

523,29 3,9232 257,18 | 773,14 3,1128 108,59

523,27 2,9366 184,96 773,16 2,0977 72,494

523,16 1,3519 117,93 773,16 1,0846 38,170

523,08 0,9878 56,410 |

Таблица АЗ - Экспериментальные значения плотности декафторбутана в жидкости и на кривой насыщения

т. к Р, МПа р, кг/мЗ Состояние

292,04 0,31328 1531,6 (рО

292,24 0,34872 1531,6 Жидкость

292,27 0,34969 1531,6 Жидкость

296,54 2,4888 1531,2 Жидкость

303,27 5,7143 1530,6 Жидкость

307,94 8,2047 1530,2 Жидкость

310,97 9,7439 1529,9 Жидкость

310,96 6,2021 1506,9 Жидкость

310,93 3,4307 1483,8 Жидкость

308,93 0,53296 1463,1 (Р')

311,57 1,50158 1462,4 Жидкость

323,21 6,3252 1461,5 Жидкость

323,21 4,1884 1441,0 Жидкость

317,61 0,68318 1422,6 <р')

323,19 2,5175 1422,6 Жидкость

333,31 6,2949 1422.0 Жидкость

343,64 10,1596 1421,0 Жидкость

343,70 8,1137 1400,8 Жидкость

343,62 5,2851 1366,1 Жидкость

328,77 0,92121 1366,4 (рО

339,36 3,9222 1366,4 Жидкость

364,92 12,0846 1365,1 Жидкость

364,97 6,4303 1287,02 Жидкость

364,96 3,7756 1227,85 Жидкость

355,37 1,74297 1227,85 (рО

356,72 2,0054 1227,9 Жидкость

375,81 6,1611 1226,8 Жидкость

391,58 9,6770 1225,8 Сверх крит. сост.

391,63 3,4320 1139,3 Сверх крит. сост.

391,66 3,0294 982,51 Сверх крит. сост.

391,77 2,6249 747,72 Сверх крит. сост.

391,62 2,4846 475,69 Сверх крит. сост.

Публикации по теме диссертации

1. Кузнецов К.И., Скородумов C.B., Сухих A.A. Экспериментальное исследование PvT поверхности рабочих тел фторуг-леродного состава. // Энергосбережение и водоиодготовка, 2009, №2, с. 28-31.

2. Сухих A.A., Кузнецов К.И., Закопырин М.А., Утенков В.Ф., Скородумов C.B. Исследование плотности декафторбутана в широком диапазоне температур и давлений // Известия ВУЗов, 2009, №7,8, с. 27-36.

3. Александров A.A., Хасаншин Т.С., Кузнецов К.й. мембранный дифференциальный манометр для работы при высоких давлениях. // Известия ВУЗов. Энергетика. Минск 1976, №5, с. 146 - 149.

4. Кузнецов К.И., Скородумов C.B., Сухих A.A.. Методика экспериментального исследования удельных объемов на кривой насыщения и в паровой фазе новых рабочих тел энергетических установок методом пьезометра постоянного объема: методическое пособие - М.: Издательский дом МЭИ, 2008. - 20 с.

Подписано в печатьЯлЛ За к. MX Тир fCC Пл Ш' Полиграфический центр МЭИ(ТУ) Красноказарменная ул.,д.13

1. Введение.

2. Обзор литературных данных по теплофизическим свойствам октафторпропана и декафторбтана.

3. Обзор методов измерения плотности.

3.1 Метод гидростатического взвешивания.

3.2 Метод пикнометра.

3.3 Метод пьезометра постоянного объема.

3.4 Метод последовательных расширений (Барнетта).

3.5 Методы косвенных измерений плотности.

3.6 Обоснование выбранного метода измерений плотности.

4. Конструкция экспериментальной установки.

4.1 Устройство жидкостного термостата.

4.2 Конструкция мембранного дифманометра.

4.3 Узлы системы измерения и регулирования температуры.

4.4 Схема измерения температуры.

4.5 Конструкция заправочного узла экспериментальной установки.

4.6 Пьезометр постоянного объема.

5. Проведение экспериментальных исследований.

5.1 Определение объема пьезометра.

5.2 Проведение аттестационных опытов.

5.3 Методика проведения опытов.

5.4 Порядок проведения измерений.

5.5 Обработка экспериментальных данных.

5.6 Пример обработки измерений плотности.

5.7 Погрешности измерений.

5.8 Чистота исследуемых веществ.

5.9 Результаты экспериментальных исследований.

6. Расчет термодинамических свойств рабочих тел фторуглеродного состава.

6.1 Построение полиномиальной зависимости.

6.2 Построение уравнения состояния для расчета термодинамических свойств декафторбутана.

6.3 Расчет термодинамических свойств октафторпропана и сравнение полученных результатов с данными литературных источников.

7. Выводы о результатах экспериментальных и расчетно-теоретических иссдований термодинамических свойств рабочих тел фторуглеродного состава.

8. Таблицы термодинамических свойств рабочих тел фторуглеродного состава.

8.1 Таблица термодинамических свойств декафторбутана.

8.2 Таблица термодинамических свойств октафторпропана.

9. Литературные источники.

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

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

Возможным становится применение фторуглеродных соединений в качестве рабочих тел второго контура атомных энергоустановках на быстрых нейтронах при использовании в первом контуре жидкометаллического теплоносителя. В настоящее время в атомной энергетике США и Франции разрабатываются и строятся двухконтурные ядерные энергетические установки, в которых преобразование энергии во втором контуре основано на цикле Брайтона. В качестве рабочих тел используются гелий и диоксид углерода. В России разрабатывается реакторная установка (РУ) «БРЕСТ-ОД-ЗОО» на быстрых нейтронах, в которой в качестве рабочего вещества (РВ) второго контура рассматривается также газ.

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

Термодинамические расчеты теплосиловых циклов показывают возможность достижения более высоких кпд по сравнению с водопаровыми циклами (Ренки-на) на рабочих телах фторуглеродного состава. Высокая молекулярная масса фторуглеродов позволяет повысить энергонапряженность оборудования и, следовательно, сделать его более компактным по сравнению с водопаровым и, тем более, газовым аналогами. Рабочие вещества (РВ) фторуглеродного класса инертны, не горючи, широко используются в качестве низко- и среднетемпера-турных теплоносителей в атомных, химических и космических технологиях. Например, смесь перфторпропана (C3F8) с элегазом (SF6) применяется в качестве хладагента (R510a) в холодильных и теплонасосных установках. Другим примером является использование фторуглеродов в климатической системе космической станции «Мир» и т.д.

Исследования последних лет показали, что рабочие вещества фторуглеродного класса имеют термическую стойкость, достаточную для реализации и высокотемпературного термодинамического цикла. Предварительный анализ показал, что применение фторуглеродов позволяет реализовать прямой сверхкритический цикл с изобарами нагрева не более 10 МПа и температурой острого пара до 600 °С, и конденсацией со сбросом теплоты в окружающую среду. Известно, что сверхкритические (треугольные циклы) имеют более высокий эксер-гетический кпд по сравнению с циклом Ренкина. Сжатие рабочего тела может осуществляться с помощью насосов (как и в водо-паровых циклах), вместо многоступенчатых компрессоров (в случае гелия, оксида углерода и смеси аргона с воздухом). Предварительный термодинамический анализ энергетических характеристик цикла на вышеперечисленных рабочих веществах (РВ) и на предлагаемых фторуглеродных РВ показывает преимущество последних. Применение в энергоустановках в качестве РВ фторуглеродов является новаторским и требует выполнения комплекса научных исследований, опытных и расчетных проработок, которые невозможно сделать без решения основных теплофизических проблем.

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

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

Для ряда рабочих тел фторуглеродного состава, отвечающих требованиям, предъявляемым к современным энергоустановкам, таким как энергетическая эффективность, низкая либо нулевая токсичность, совместимость с конструкционными материалами и отвечающими экологическим требованиям Киотского и Монреальского протоколов к настоящему времени уже существуют данные, необходимые для расчетов энергетических циклов. Например, для тетрафторэтана (CF4), циклофторбутана (C4F8) были выполнены оценочные расчеты термодинамических циклов. Основные константы фторуглеродов даны в таблице 1. Вариант такого цикла представлен на рисунке 1. Например, при 773,15 температуре на входе в турбину 530 °С в сверхкритическом цикле <5"з,15 с регенерацией теплоты от ^ температуры выхода из н питающего насоса с ис-*'3'15 пользованием циклофторбутана (C4F8) термический КПД составил 46%.!

373,15

273,15

Рис. 1.1 - Сверхкритический теплосиловой цикл на циклофторбутане (C4F8)

Таблица 1.1— Основные константы фторуглеродов

Вещество Обознач. Мол. масса, кг/кмоль Критические параметры тнк> к Тпл, к

Тс, к Рс, МПа Р с, кг/м3 cf4 r14 88,005 227,5 3,745 630,0 145,12 89,53 c2f6 r116 138,01 293,03 3,040 620,0 195,00 173,10 c3f8 r218 188,02 345,05 2,670 615,0 236,5 125,0 c4f8 rc318 200,04 388,47 2,780 547,9 267,18 232,95 c4f10 r31(10) 238,02 386,35 2,324 600,0 271,15 145,0 c5f12 r41(12) 288,03 421,85 2,048 619,4 302,45 148,2

Высокая термодинамическая эффективность применения в качестве рабочих тел фторуглеродного состава (октафторпропана СзГ8, циклофторбутана c4f8, декафторбутана c4f10) для энергетических установок с бинарными циклами впервые была отмечена в работах Гохштейна д.П. и соавторов [2] ещё в 60-ые годы прошлого века. Однако, расчеты циклов с приемлемой точностью могут быть выполнены только по результатам настоящей работы.

7. ВЫВОДЫ О РЕЗУЛЬТАТАХ ЭКСПЕРИМЕНТАЛЬНЫХ И РАСЧЕТНО-ТЕОРЕТИЧЕСКИХ ИССДОВАНИЙ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ РАБОЧИХ ТЕЛ ФТОРУГЛЕРОДНОГО СОСТАВА

1. Модернизирована экспериментальная установка для исследования плотности веществ в жидком и газообразном агрегатном состоянии в диапазоне температур 25.500 °С и давлений от атмосферного до 10 МПа.

2. Проведено экспериментальное измерение плотности декафторбутана (C4F10) в широком диапазоне параметров состояния с погрешностью измерений без вероятных ошибок отнесения не хуже 0,24% во всей области исследований. Погрешности измерений с учетом вероятных ошибок отнесения при давлениях ниже 2 МПа не превышают 0,35%, в надкритической области (Т=293,15.413,15К, р=2.3,5 МПа) не превышают 0,45%. В остальной, рабочей области параметров состояния, используемых при расчетах энергетических установок, полная погрешность измерений не превышает 0,3%. Полная погрешность измерения плотности в жидкой фазе не превышает 0,25%.

3. Проведено экспериментальное измерение плотности октафторпропана (СзР8) в широком диапазоне параметров состояния с погрешностью измерений без вероятных ошибок отнесения не хуже 0,24% во всей области исследований. Погрешности измерений с учетом вероятных ошибок отнесения при давлениях ниже 2 МПа не превышают 0,35%, в надкритической области (Т=248,15.353,15К, р=2,8.3,8 МПа) не превышают 0,45%. В остальной, рабочей области параметров состояния, используемых при расчетах энергетических установок, полная погрешность измерений не превышает 0,3%.

4. Разработаны уравнения для расчета термодинамических свойств октафr S;

7 = 1 + W Ъ юЧ"-1 торпропана и декафторбутана вида ** Zj Zj и ' работающие в

5=1 j=o диапазоне перегретого пара и сверхкритического состояния.

5. Построено устойчивое уравнение для описания термодинамических свойств декафторбутана с использованием методики эквивалентных урав нений состояния (15 эквивалентных уравнений).

6. Проведено расчетно-теоретическое исследование основных термодинамических свойств декафторбутана (C4F10) и октафторпропана (C3F8) с использованием экспериментальных данных настоящей работы и экспериментальных данных других авторов. Что позволило создать уравнения, хорошо описывающие свойства исследуемых веществ в широком диапазоне параметров состояния.

7. ТАБЛИЦЫ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ РАБОЧИХ ТЕЛ

ФТОРУГЛЕРОДНОГО СОСТАВА Обозначения табличных величин: Р - давление; р — плотность; Z - сжимаемость; h - энтальпия, S — энтропия; Ср — изобарная теплоемкость; Cv - изо-хорная теплоемкость; W - скорость звука; Т - температура.

1. Некоторые особенности парогазовых схем с неводяными парами / Гох-штейн Д.П., Смирнов Г.Ф., Киров B.C.// Теплоэнергетика. 1966, №1, С. 20-24.

2. Расчетно-экспериметальные данные о теплофизических свойствах пер-фторпропана и смесей пяти низкокипящих фторсоединений: / Отчет МЭИ; Руководитель работы Алтунин В.В. № ГР 0191.0016142; №0293.0001614. - Москва, 1992. -70с.

3. Промышленные фторорганические продукты: Справочник / Б.Н. Максимов, В.Г. Баранов, И.Л. Серушкин и др. Л.: Химия, 1990

4. Барышев В.П. Комплексное исследование теплофизических свойств фреона-218: автореф. дис. канд.техн. наук.ОТИХП- Одесса, 1982.-20с.

5. Единое неаналитическое уравнение состояния хладона R218 / Козлов А.Д., Лысенков В.Ф., Рыков В.А. // ИФЖ 1992.Т.62,№6, с.840-847.

6. Перельштейн И.И., Парушин Е.Б. Термодинамические и теплофизиче-ские свойства рабочих веществ холодильных машин и тепловых насосов. М: Легкая и пищевая промышленность. 1984. - 232 с.

7. Комплексные исследования теплофизических свойств озонобезопасных хладагентов / Кессельман П.М., Железный В.П. // Холод, техника. 1992.№ 11/12. с.16-18.

8. Physical and chemical properties of pure fluorocarbons / R.D. Fowler, J.M. Hamilton, J.S. Kasper et al // Ind. Eng. Chem. 1947. Vol. 39, N 3, P. 375378.

9. The properties of n-butforane and its mixtures with n-butane / Simons J.H., Mausteller J.W. // J.Chem. Phus. 1952. Vol. 20, N 10. p.l 516-1519.

10. Physical properties of n-perfluorobutane / Brown J.A., Mears W.H. // J. Phys. 1958. Vol. 62, N 10. P.960-962.

11. Second viral coefficients for the system n-butane + perfluoro-n-butane / Tripp T.B., Dunlap R.D. // J. Phus. Chem. 1962. Vol. 66, N 4. P.635-639.

12. Вторые вириальные коэффициенты органических соединений и их смесей / Маркузин Н.П. // Химия и термодинамика растворов 1968. Вып. 2. с. 212-238.

13. Vapor pressures of pure substances. Amsterdam / Boublik Т., Fried V., Ha-la E. // Elsevier sci.publ.corp. 1973. P.146-147.

14. Surface tension of perfluoropropane, perfluorobutane, . / McClure I.A., Soares V.A., Edmonds B. // J. Phus. Chem., 1960. Vol. 51, N 1. P.535-639

15. Таблицы стандартных справочных данных. Хладон R218. Плотность, энтальпия, энтропия, изобарная и изохорная теплоемкости / Рыков В.А,

16. Устюжанин Е.Е., Попов П.В., Кудрявцева И.В., Рыков С.В. ГСССД 211-05 // ГССД, М., 2005

17. TRC Thermodynamic Tables, Non-Hydrocarbons. JANAF Thermochemical Tables 4th floor stacks College Station,TX: Thermodynamics Research Center, Texas A&M University System, QD511 N57 1998, 2985- (QD305 H5 T45) 9 volumes, loose leaf Page s-6520

18. Тепло и массообмен. Теплотехнический эксперимент. Справочник. Е.В.Аметистов, В.А.Григорьев, Б.Т.Емцев и др.; под общей редакцией В.А.Григорьева и В.М.Зорина. - М.: энергоиздат, 1982г 512с.

19. Гидростатический метод определения плотности жидкостей при давлениях до 5000 кгс/см" / В.Н.Разумихин // Труды институтов Комитета стандартов, мер и измерительных приборов, М., 1960, т. 46 (106), с. 96 -106.

20. Установка для измерения плотности жидкостей гидростатическим методом при давлениях до 10000 кгс/см" / В.А.Борзунов, В.Н.Разумихин // Труды институтов Комитета стандартов, мер и измерительных приборов, М., 1964, т. 75 (135), с.143 142.

21. Исследование термодинамических свойств веществ / В.А.Кирилин, А.Е.Шейндлин М. - Л.: Госэнергоиздат, 1963, 559 с.

22. Теплофизические свойства органических теплоносителей, / М.П.Вукалович, Ю.М.Бабиков, Д.С.Рассказов, М.: Атомиздат, 1970, 236 с.

23. С.Л.Ривкин, Е.А.Кременевская. Исследование плотности фреона 11. -Теплофизические свойства веществ и материалов // ГСССД, 1971, вып.4, с. 5 —17.

24. Экспериментальное исследование P-V-T зависимости бутил- бутилового эфира в широком интервале параметров состояния / К.Д.Гусейнов, Н.М.Байрамов // Журнал физической химии, вып. «Химическая термодинамика и термохимия»М., 1979, с. 64-65.

25. Экспериментальное определение удельных объемов тяжелой воды в интервале температур от 200 до 425°С при давлениях до ЮООбар / Н.В.Цедерберг, А.А.Александров, Т.С.Хасаншин, Д.К.Ларкин. // Теплоэнергетика, М., 1973, №8, с. 13-15.

26. Экспериментальное исследование удельных объемов тяжелой воды в интервале температур от 20 до 200°С и при давлениях до 1000 бар / Н.В.Цедерберг, А.А.Александров, Т.С.Хасаншин // Теплоэнергетика, М., 1972, №10, с.65-69.

27. Утенков В.Ф., Филатов Н.Я. Экспериментальное и расчетно-теоретическое исследование удельных объемов хладона R151A.: Химия и компьютерное моделирование. Бутлеровские сообщения 2002. Приложение к спецвыпуску №10.

28. Техническая термодинамика: учебник для вузов. В.А.Кирилин,

29. B.В.Сычев, А.Е.Шейндлин. 4-е изд., - М.: Энергоатомиздат, 1983г. 416с.

30. Исследование скорости распространения звука в тяжелой воде при высоких давлениях / А.А.Александров, А.И.Кочетков. — Теплофизические свойства веществ и материалов // ГСССД, 1979, вып. 14, с. 60 — 62.

31. К оценке погрешности измерения скорости звука в зависимости от точности регистрации сложения сигналов / Кочетков А.И., Кузнецов К.И., Кузнецов Г.Г. //Термодинамика и теплофизические свойства веществ: Сб. науч. Трудов № 131. М.: Моск.энерг.ин-т. 1987.

32. Физические свойства сталей и сплавов, применяемых в энергетике. Под ред. Неймарк Б.Е. М.: Изд-во «Энергия», М.: 1967,120 стр.

33. Теплофизические свойства материалов ядерной техники Чиркин B.C. -М.: Изд-во «Атомиздат», 1968, 484 стр.

34. Исследование термодинамических свойств веществ / В.А.Кирилин, А.Е.Шейндлин. М. Л., Госэнергоиздат, 1963, 559 с

35. Мембранный дифференциальный манометр для работы при высоких давлениях / Александров А.А., Хасаншин Т.С., Кузнецов К.И. // Известия ВУЗов. Энергетика. Минск 1976, №5, с. 146 149.

36. Зарубежные интегральные микросхемы для промышленной электронной аппаратуры / Нефедов А.В. М.,Издательство: Энергоатомиздат, 1989, 288 стр.

37. Лурье А.И. Теория упругости. Изд-во «Наука», М., 1970, 675 стр.

38. Экспериментальное исследование PVT поверхности рабочих тел фто-руглеродного состава / Кузнецов К.И., Скородумов С.В., Сухих // А.А. Энергосбережение и водоподготовка. М., 2009, №2, с. 28 31.

39. В.В. Алтунин. Теплофизические свойства двуокиси углерода. М., Издательство стандартов, 1975, с. 546.

40. Термическое уравненме состояния шестифтористой серы в широкой области параметров состояния / Присяжный А.П., Тоцкий Е.Е., Устюжанин Е.Е. // Теплофизика высоких температур. 1989. - Т. 27, №21. C.400 403.

41. Таблицы теплофизических свойств воды и водяного пара / Александров А.А., Григорьев Б.А.: Справочник. Рек. Гос. Службой стандартных справочных данных. ГСССД Р-776-98. — 2-е изд., стереот. М.: Издательский дом МЭИ, 2006.

42. Теплофизические свойства ртути: Рек. Гос. Службой стандартных справочных данных. ГСССД / Вукалович М.П., Иванов А.И., Фокин Л.Р. и др. М., Издательство стандартов, 1971.

43. Исследование плотности декафторбутана в широком диапазоне температур и давлений / Сухих А.А., Кузнецов К.И., Закопырин М.А., Утенков В.Ф., Скородумов С.В. // Известия ВУЗов. Проблемы энергетики, 2009, №7,8, стр. 27-36.

44. Brown I.A Physical properties of perfluoropropane//J. Chem. Eng. Dta.-1963. Vol.8, №11. P.106-108.

45. Термодинамические свойства азота. / B.B. Сычев, А.А. Вассерман, А.Д. Козлов, Г.А. Спиридонов В.А. Цымарный М., Издательство стандартов, 1977, с. 352.

46. Heat Capacity of Gaseous Perfluoropropane / Masi J.F., Flieger H.W., Wick-lund J.S. Journal of Research of the National Bureau of Standards. Vol.52, No.5, 1954.

47. Lemmon E.W. and Span R., «Short Fundamental Equations of State for 20 Industrial Fluids» E.W. Lemmon, R., J. Span Chemical Engineering. Data, 51: p. 785-850, 2006.

48. Плотность и уравнение состояния фреона- 218 / Геллер В.З., Поричан-ский Е.Г., Барышев В.П. // Известия ВУЗов. Серия «Энергетика».-1980. №6.-с.119-123.