УДК 536.71
Анализ экспериментальной информации о равновесных свойствах r218 на основе неаналитического уравнения состояния
Рыков А.В., Кудрявцева И.В., Рыков С.В.
Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики
В статье проведен анализ новых экспериментальных данных о теплофизиче-ских свойствах R218 и их сравнение с опубликованными ранее. Сравнение проведено с помощью неаналитического уравнения состояния. Приведена структура уравнения состояния.
Ключевые слова: уравнение состояния, термодинамические свойства, R218.
Появление новой экспериментальной информации дает возможность уточнять имеющиеся уравнения состояния. Однако важным является проверка новых экспериментальных данных на согласованность с имеющимся массивом опытных данных. В работе [1] приводятся новые термические для хладона R218. При этом авторы [1] показали, что эти данные хорошо согласуются с расчетными значениями, полученными на основе уравнения состояния Леммона и Спана [2]. С другой стороны, уравнения состояния Барышева [3] и Рыкова [4] в области плотной жидкости существенно расходятся с результатами, полученными на основе [2].
Для анализа экспериментальной информации о термических и калорических свойствах R218 воспользуемся уравнением состояния, которое качественно верно описывает как регулярную часть термодинамической поверхности, так и область сильно развитых флуктуаций:
Здесь Р - плотность; T - абсолютная температура; t=T / Tc ; t = t - 1; t1 = 1/1 - 1; w = p / p c ; А p = w - 1; Tc - критическая температура; Р c - критическая плотность; R - газовая постоянная; a и у - критические индексы изохорной теплоемкости Cv и изотермической сжимаемости КТ, соответственно; ß и d - критические индексы кривой сосуществования и критической изотермы, соответственно; f (w) и f j (t ) - сглаживающие функции; F0 ( T ) - идеально-газовая составляющая свободной энергии; x = t /1А р |1ß - масштабная переменная.
Термическое уравнение состояния, рассчитанное на базе выражения свободной энергии Гельмгольца (1) имеет следующий вид:
где
Здесь % (р ,Т) - сжимаемость; Ь (х) - масштабные функции химического потенциала.
Сглаживающие функции f (®) и /у() в (1), (2) описываются следующими зависимостями:
(3)
В данной работе использованы следующие масштабные функции ai (х):
£23(х) = Д (х + + С4
(4)
+ С
где = Д3=Д4=у+М;
Коэффициенты Д и Д находятся из равенств:
Здесь
Постоянные С1, входящие в масштабные функции (6), находятся из равен-
ства (2 - а + А 1) а (х = - х0) + х^ (х = - х0) = 0.
Для того чтобы термическое уравнение состояния (2) передавало равновесные свойства R218 в соответствии с требованиями МТ, коэффициенты с¿,^^ выбирались таким образом, чтобы выполнялись равенства:
где п = 1,2,3,4.
С учетом (5) уравнение состояния (2) принимает вид:
+оо^Сг.0 (Лр)! (¿со+ Ар)+ сот;1 (2со-3)Си +
1=6
+со2т1 (Зсо- 4 )С21 + оот^Сд (Ар)! 1 (¡'оо+ Др) +
"3 ъ
(О
х
¡'-О 3=2
(6)
"1 "2
г=0 ;'=0
где функции у1, у2 и у3 описываются зависимостями:
=-15>4/12+5,8/12Др-1Д/6(Др)2+0,05(Др)'
72=5-4Др + 3(Др)2-2(Др)3 + (Др)4
>-3 = 4-ЗДр+2(Др)2-(Др) + (Др)"
Для хладона Я218 имеется экспериментальная информация об изохорной теплоемкости [5, 6], полученная как в критической области, так и в регулярной.
Коэффициенты уравнения (1) определяются из условия минимума функционала
Для расчета идеально-газовых составляющих термодинамических функций использовано уравнение для идеально-газовой составляющей изобарной теплоемкости, предложенное в [6]:
С° = 0Д170+ 2,908 Ю-3 Г- 2,226 -10^Г2, кДж/СкгЕ).
Были рассмотрены варианты, в которых вес для данных [1] на квазиизохорах менялся от 1,7 до 0,7 (соответственно, вариант I и вариант II). На рис. 1 - рис. 6. представлены полученные результаты. При этом вес для данных Барышева на изотермах 133,15 К - 313,15 К выбран равным 1. На основе анализа рис. 1 видно, что увеличение веса экспериментальных данных [1] привело к существенному улучшению (до 0,28 %) описания этих данных уравнением (1) на квазиизохорах 1746 кг/м3 и 1795 кг/м3.
При увеличении веса данных [1] не произошло существенного изменения описания термических [5] и калорических экспериментальных данных [4, 6] относительно варианта I (см. рис. 2, 5, 6).
Из рис. 4. видно, что в случае варианта I экспериментальные данные Барышева описаны значительно хуже, чем при малом весе данных [1] в интервале температур 213,15 К - 253,15 К.
На основе полученной информации можно сделать вывод, что новая экспериментальная информация [1] не противоречит ранее полученной. И ее можно включать в массив экспериментальных данных при построении уравнения состояния. Однако, следует иметь в виду, что в случае включения в массив исходной информации экспериментальных данных [1] с весом больше 1 приводит к существенному снижению точности описания линии фазового равновесия в окрестности тройной точки. Это обусловлено тем, что расхождение между данными [1] и [3] при в диапазоне температур 213,15 К - 253,15 К достигает 0,4% по плотности и носит систематический характер, а это приводит в области плотной жидкости к значительному расхождению по давлению.
8р, % 0,4
0,3
0,2
0,1
0
-0,1
-0,2
-0,3
-0,4
-0,5
200 220 240 260 280 300 320 340 Т, К
Рис. 1. Отклонения расчетных значений от экспериментальных данных [1]:
1 - варианта I; 2 - варианта II.
8р, % 4
3
2
1
0
-1
-2
320 340 360 380 400 420 440 Т, К
$ 01
02
о
о о о о
О $ о я % 6 А 8 , •
о О В л О т« <*> СО с
8
Рис. 2. Отклонения расчетных значений от экспериментальных данных [5]:
1 - варианта I; 2 - варианта II.
8р, %
1
0,5
0
-0,5
-1
-1,5
280 300 320 340 360 380 400 420 Т, К
01 02 О о о о О Л
о о $ О е
о Л О ЛО л о° А, Л О
Ч> в 8 3 М4' Л ><» о о 8 о^ о О < оо о о о о 0
; о о о о < ' < о о
< < о > к
Рис. 3. Отклонения расчетных значений от экспериментальных данных [3]:
1 - варианта I; 2 - варианта II.
8р, % 0,2
0,1
0
-0,1
-0,2
-0,3
-0,4
о & о
Й 0 о
8 $ с о Й* 0
ос о о Я о о Л Л о о° о о ^о ЛО
о о о о < о V
01 3>о о о О О г о
02 о О оо
1,3
1,4
1,5
1,6
1,7
1,9 р, г/см3
Рис. 4. Отклонения расчетных значений от экспериментальных данных [3]:
1 - варианта I; 2 - варианта II.
§ Cv, % 2
1,5
1
0,5 0
-0,5 -1 -1,5 -2
$
О 1
в ^ о 2
о«
<> *
о ^о^ ft toi < 8 $ ov о
о >v о о О о л
О о * > фО * 0
#
344
346
348
350
352
354
T,, K
Рис. 5. Отклонения расчетных значений от экспериментальных данных [4]:
1 - варианта I; 2 - варианта II.
§ Cv, % 2,5
2
1,5
1
0,5 0
-0,5 -1 -1,5
о
о О 1
о о 2
о о о
о * 0 <р>
о V
8
330
340
350
360
370
380
T, K
Рис. 6. Отклонения расчетных значений от экспериментальных данных [6]:
1 - варианта I; 2 - варианта II.
Список литературы:
1. Klomfar J., SoucC kovaa M, Paatek J. Experimental study of p-p-T relationship of compressed liquid phase for octafluoropropane and two near azeotropic ternary HFC/HC mixtures // J. Chem. Eng. Data. 2012. Vol. 57. P. 1627-1634.
2. Lemmon E. W., Span R. Short Fundamental Equations of State for 20 Industrial Fluids // J. Chem. Eng. Data 2006, 51, 785-850.
3. Барышев В.П. Комплексное исследование теплофизических свойств фреона -218 //Дис. ... канд. техн. наук. Л., 1981. С. 204.
4. Рыков В.А. Единое неаналитическое уравнение состояния газа и жидкости и таблицы термодинамических свойств аргона и хладагентов R134а, R218, R134а // Дис. на соискание уч.ст. докт. техн. наук. СПб.: СПбГУНиПТ, 2000. -456 с.
5. Brown I.A. Physical properties of perfluoropropane//J. Chem. Eng. Data. 1963. Vol. 8, № 11. P.106-108.
6. Рябушева Т.И. Исследование изохорной теплоемкости холодильных агентов // Дис. ... канд. техн. наук. Л., 1979. С. 189.
The analysis of the experimental information on equilibrium properties r218 on the basis of not analytical equation of state
Rykov A.V., Kudryavtseva I.V., Rykov S.V.
National Research University of Information Technologies, Mechanics and Optics
In article the analysis of new experimental data about thermalphysic properties R218 and their comparison with published earlier is carried out. Comparison is spent by means of not analytical equation of state. The structure of an equation of state is given.
Key words: equation of state, thermodynamic properties, R218.