CC BY
6
УДК 536.71
Анализ различных моделей среднего диаметра линии фазового равновесия R236ea
Канд. техн. наук С. В. РЫКОВ1, канд. техн. наук И. В. КУДРЯВЦЕВА, д-р техн. наук В. А. РЫКОВ, А. В. СВЕРДЛОВ, М. НУРЫШЕВА
1togg1@yandex.ru Университет ИТМО
Выполнен анализ различных вариантов линии фазового равновесия холодильного агента R236ea, разработанных в рамках модели среднего диаметра Вегнера ([1-а] -модели) и [2ß] -модели. Особенность предложенного подхода заключается в том, что плотность насыщенного пара р- и во всех трех моделях среднего диаметра fd описывается одним и тем же уравнением Клапейрона - Клаузиуса, в которое вместо теплоты парообразования r введена «кажущаяся» теплота парообразования r*: r*=r/[1 - р -/р+], в то время как плотность насыщенной жидкости р+ описывается в этих моделях разными уравнениями, но имеющими один и тот же набор слагаемых с одинаковыми показателями степени. Таким образом, предложенная система уравнений для линии фазового равновесия имеет физически обоснованную структуру и позволяет в модели среднего диаметра fdучесть все компоненты, включая линейный компонент т=1 - T/T. Перечисленные модели среднего диаметра апробированы при расчете линии фазового равновесия R236ea в диапазоне температур от Тг (тройная точка) до Tc (критическая точка). Ключевые слова: линия упругости, «кажущаяся» теплота парообразования, хладон R236ea, диаметр порядка, средний диаметр, линия фазового равновесия, критические индексы.
Информация о статье:
Поступила в редакцию 12.05.2019, принята к печати 20.08.2019 DOI: 10.17586/1606-4313-2019-18-3-87-93 Язык статьи — русский Для цитирования:
Рыков С. В., Кудрявцева И. В., Рыков В. А., Свердлов А. В., Нурышева М. Анализ различных моделей среднего диаметра линии фазового равновесия R236ea. // Вестник Международной академии холода. 2019. № 3. С. 87-93.
Analysis of various models of the average diameter of phase equilibrium line R236ea
Ph. D. S. V. RYKOV1, Ph. D. I. V. KUDRYAVTSEVA, D. Sc. V. A. RYKOV, A. V. SVERDLOV, M. NURYSHEVA
'toggl@yandex.ru
ITMO University
The analysis of various variants of the phase equilibrium line of the refrigerant R236ea, developed in the framework of the model of the average diameter of the Wegner or [1 - a] -model, [2p] -model, is performed. The peculiarity of the proposed approach is that the density of saturated steam p- and in all three models of average diameter fd is described by the same Clapeyron-Clausius equation, in which instead ofthe heat of vaporization the "apparent" heat of vaporization r* is introduced: r*=r/ [1 - p -/p+] while the density of the saturated liquid p+ is described in these models with different equations, but having the same set of terms with the same exponents. Thus, the proposed system of equations for the phase equilibrium line has a physically justified structure and makes it possible to take into account all components in the average diameter fd model, including the linear component t=1 - T/Tc. The listed models of average diameter were tested when calculating the phase equilibrium line R236ea in the temperature range from Tr (triplepoint) to Tc (criticalpoint).
Keywords: line of elasticity, "apparent" heat of vaporization, R236ea, diameter of the order, average diameter, line of phase equilibrium, critical indices.
Article info:
Received 12/05/2019, accepted 20/08/2019 DOI: 10.17586/1606-4313-2019-18-3-87-93 Article in Russian For citation:
Rykov S. V, Kudryavtseva I. V., Rykov V. A., Sverdlov A. V., Nurysheva M. Analysis of various models of the average diameter of phase equilibrium line R236ea. VestnikMezhdunarodnoi akademii kholoda. 2019. No 3. p. 87-93.
Введение
Информация о равновесных свойствах холодильных агентов и других технически важных веществ на линии фазового равновесия, используемых в энергетике, имеет важное практическое значение. Получение этой информации во всем диапазоне линии насыщения, от температуры Тг тройной точки до температуры в критической точке Тс является сложной и дорогостоящей задачей. Поэтому для расчета параметров линии насыщения, линии упругости и теплоты парообразования используются расчетные методы [1-5], которые не требуют наличия информации о плотности р- и давлении р, насыщенного пара, плотности р+ насыщенной жидкости и теплоте парообразования г во всем диапазоне температур от Т,г до Тс.
Чтобы неопределенность расчетных методов была минимальной, необходимо использовать такие уравнения для описания р - , р+, р, и г, которые качественно верно, то есть в соответствии с требованиями современной физики, описывают поведение линии фазового равновесия от тройной точки до критической точки. В частности, к таким требованиям относится качественное и количественно верное описание, в рамках используемой модели линии насыщения, поведения среднего диаметра
d = Р-+£~1, 2рс
(1)
(T ) = T • dPs-,
v ' r dT
где р! — давление насыщенного пара; г* — «кажущаяся» теплота парообразования, которая связана с теплотой парообразования г зависимостью:
r = r
(
1 -Р-
Р+,
(5)
Привлекательность описания плотности р- (Т), на основе формулы (4), обусловлена тем, что и теплота г*, и давление р, измеряются экспериментальным путем [14]. При этом следует учитывать, что характер поведения функции рх (Т) теоретически и экспериментально обоснован и в окрестности критической точки [10]:
ps(t) = pc 1 +a1t + a2t2 a+o(t2 a)
(6)
и в непосредственной близости от тройной точки [15]:
Ps (T) = Рс exp I - a0
(7)
где а0 > 0 — индивидуальный параметр; р , — критическое давление.
В работе [14] предложена зависимость р, (Т), которая удовлетворяет требованиям (6) и (7):
Ps(T)=Pc^x2/t
где р+ — плотность насыщенной жидкости; р- — плотность насыщенного пара; рс — критическая плотность.
При описании линии фазового равновесия до 2000 г. использовались фактически две модели среднего диаметр. Во-первых, это линейная модель [9, 10], в рамках которой функция fd описывается зависимостью:
fd (T ®TC ) = At + o (t), (2)
где т= 1 - T/T; T — абсолютная температура; о — символ Ландау.
Теория критических явлений, бурный расцвет которой пришелся на вторую половину прошлого века, на основе ренорм-группового анализа и e -разложения [11] предсказала другую модель поведения криволинейного диаметра [12, 13]:
fd (T ® Tc) = At1-a + o(t1-a), (3)
где а — критический индекс изохорной теплоемкости, а -0,1 [11].
Модель среднего диаметра (3) была предложена Вег-нером и в дальнейшем мы будем обозначать ее как [1 - а]-модель.
Модель среднего диаметра [2ß ]
Модель [1 - а] до 2000 г. была признана единственной теоретически обоснованной моделью среднего диаметра линии насыщения и, вследствие этого, получила широкое распространение. Видимо поэтому в то время не привлекла внимание исследователей работа [14], выполненная в ЛТИХП (г. Ленинград, 1988 г.), в которой было предложено описывать плотность паровой ветви линии насыщения на основе уравнения:
l + qT + a2W2 а +
М
+ fl3k |2"а+А + 1> г=4
МО
(8)
где а1 — постоянные коэффициенты, i е{1,2,3,...,М}; п (г) е N ; t = Т / Тс; А — первый неасимптотический критический индекс: А - 0.5 .
Выражение для г* в [14] выбрано таким образом, чтобы удовлетворить требованиям масштабной теории и позволить количественно с малой неопределенностью описать опытные данные о «кажущейся» теплоте парообразования R218:
r*(T) = ^ (d0 + d1tb + d3 tb-D +d4 t1-a + Xd-1 ^
рс
i=4
(9)
где di — постоянные коэффициенты; у (г) е N.
Заметим, что выражение (9) использовалось в работах [16-18] при описании свойств различных холодильных агентов.
Система уравнений (4), (8) и (9) позволяет передать плотность р- на паровой ветви линии насыщения с малой неопределенностью, не превышающей неопределенности опытных данных [19] (см. [20]), и в соответствии с моделью среднего диаметра [2р]:
fd (T ® Tc ) = At2b+ B t1-a+ o (t1-a)
(10)
описать поведение р
Р-1 Т =РсЦ +Dltb + D2t2b + Dзtb+А+DAtl-a + Л. (11)
Начиная с 2000 г., благодаря теоретическим работам [21, 22], модели (10) и (11) получили признание и начали исследоваться многими авторами [2-4, 23-25].
В работе [26] показано, что точность описания р - , на основе уравнения (4), можно повысить, если в функцию (9) ввести дополнительное слагаемое, пропорциональное т2в:
V
Pc
+ </3к1Р+А+^|тГ" +
1_а ■
i=5
Р+(Т) = РС
r2ß.
d
d3
do = a ; Di =-i, D2 =-3, D3 =
ж dL _*2L
do2 do
D4 =-
Ж "l3
3
di_ -2 did2
do
do
Таким образом, выбирая в (13) коэффициенты D*
*
и D4 в соответствии с (14):
^ ¿2 * у3 = ^ = -~г и А
D* = D3 = 4 - ^ и d4 = D4 = - 2 dld2
do2 do
^ d 3 d0
do
D* =-D3 =-
Ж di2
ц
и d4 = d4 = d3 - 2 ä
ioo
фА=е
n=1
ioo
ф =е p
n=i
Q ( pex - pras *ips,n\rs n is n
Q (p -,ex -p-,ras
(17)
(18)
где Q
(12)
Ps n
и Qp - n — вес соответствующей n -ой опытной
где 5 (0=1+(5 - 0 р.
Заметим, что введение дополнительного слагаемого в (12) не приводит к изменению структуры выражения (11) вблизи критической точки.
Теперь для описания линии фазового равновесия в соответствии с моделью (10) необходимо использовать уравнение для жидкостной ветви линии насыщения, имеющего ту же структуру при Т ® Тс, что и уравнение (11). Поэтому мы выбрали выражение для функции р+ (Т) в виде:
г
(13)
10
+ В4х1~а + £5т3(3 + £ £гт1+( )а
1=6 , Наш анализ выражений (4), (8) и (11) показал, что коэффициенты этих уравнений связаны между собой следующими равенствами:
ех _ех га ? _га ?
точки р8 п или р п' ; р5п и г„ — значения давления и плотности, рассчитанные соответственно по уравнениям (7) и (4), (12).
При этом критические индексы выбраны в соответствии с масштабной теорией [11]: а=0,11, р = 0,325 и Д=0,51. В ходе поиска коэффициентов уравнения (7) на основе минимизации функционалов (17) и (18), мы установили значения критических параметров R236eа: р5=3,41695 МПа, р5=563 кг/м3, Гс=412,3801 К.
В результате получены значения коэффициентов и параметров уравнений (7) и (4), (12), которые приведены в табл. 1 и 2.
Результаты сравнительного анализа опытных данных [27, 29, 32, 34, 35] и расчетных значений р™ и р_га" представлены на рис. 1, 2.
Таблица 1
Коэффициенты и параметры (7)
Table 1
Coefficients and parameters (7)
(14)
i a, n (,) i a, n (i)
0 13,8 — 5 41,56405515467 3
1 8,691970045447 — 6 -33,60392555014 5
2 207,8636296092 — 7 35,90836752245 7
3 92,62472654246 — 8 -36,44499800073 9
4 -262,2190347318 2 — — —
(15)
Коэффициенты и параметры (12) Coefficients and parameters (12)
Таблица 2
Table 2
получим модель среднего диаметра (10).
Если в (13) коэффициенты и D4* удовлетворяют равенствам:
i d, i d,
0 8,691970045447 3 -42,8041137641
1 15,0479824374 4 29,2631979938
2 17,6143631770 5 17,8226277957
(16)
Коэффициенты уравнения (13) устанавливались на основе данных [27, 28, 34], путем поиска минимума следующего функционала:
то получим для среднего диаметра ^ модель [1 - а] (3).
Линия фазового равновесия R236ea
На основе уравнений (4), (7), (12) и (13) мы описали для R236ea данные о давлениир5 [27-32], плотности р+ и р- [27-29, 34-35] для двух вариантов
— вариант I: выполняются равенства (15), модель [2р];
— вариант II: выполняются равенства (16), модель [1 - а].
Коэффициенты уравнений (4), (7) и (12) устанавливались на основе данных работ [27, 29, 32, 34, 35], путем поиска минимума функционалов:
Ф
100 =х
Qr + n (p+ex-
p+ras n
(19)
где Qp +n
вес соответствующей п-й опытной точки
+ ех + га ?
Р п ; Рп' — значения плотности, рассчитанные по уравнению (13).
В результате мы получили следующие значения коэффициентов и параметров уравнения (13), которые приведены в табл. 3 (вариант [2р] (15)) и табл. 4 (вариант [1 - а] (16)). На рис. 3 мы ограничились представлением информации о точности описания р+'ех [27, 28, 33, 34] в рамках модели [2р], так как модель [1 - а] опи-
o
o
8p-, % 5
0
-5
-10
-15
-20
-25
-30
-35
-40
150 200 250 300 350 400 T, K
pmr-'g
□ 0
0
01 □2 Д3 Х4 Ж5 Об О
Рис. 1. Отклонения 5ps = ( pf - prsas ) / pf^ 100% давления prsas на линии упругости, рассчитанные поуравнению (8), от известныгх значенийpesx [27-33]: 1 — [27]; 2 — [28]; 3 — [28)/; 4 — [29]; 5 — [30]; 6 — [31]; 7 —[32]; 8 — [33]; 9 — [28]
Fig. 1. Pt-essure 5ps = 0pf-prsas ) / pf ■ 100%% deviations p™ on the elasticity line calculated by equation (88)) from pei known values [27-33]: 1 — [272; 2 — [28]; 3 — [28]; 4 — [29]; 5 — [30]; 6 — [31 ]; 7 — [32]; 8 — [33r;9 — [28]
Puc. 3. OmHocumejibHbie omKJioHeHun 5p+ = (p+,ex -p+'ras ) / p+,ex ■ 100% momnocmu p+/as HanuHuu Hacu^eHun, paccuumatmbe noypaBHemmo (13),om u3BecmHux SHaietiuu p+'ex [27, 28, 33, 34]: 1 — [34]; 2 — [27]; 3 — [28]; 4-[33]
Fig. 3. Relative 5p+ = ( p+^ - p+'ras ) / p+ ^ ■ 100% density p+'ras
deviations on the saturation line calculated by equation (13) from p+'ex known values [27, 28, 33, 34]: 1 — [34]; 2 — [27]; 3 — [28]; 4 — [33]
Puc. 2. OmmoHemn Sp- = (p~,£x -p-'ras) / p~'ex ■ 100%nnomHocmu p-'ras Ha nmuu Hacbiuf.eHUM, pacciumaHHbie no ypaeHeHuw (13), om u3eecmHbix 3HaieHuu p-,£% [27, 29, 33-35]: 1 — [34]; 2 —
[27]; 3 -7 [27]; 4 — [35]; 5 — [29]; 6 — [33] Fig. 2. Density Sp- = (p-ex - p-'ras) / p-ex ■ 100% deviations p-'ras
on the saturation line calculated by equation (13) from p-'ex known
values [27, 29, 33-35]: 1 — [34]; 2 — [27]; 3 — [27]; 4 — [35]; 5 — [29]; 6 — [33]
Ta6m^ 3
Коэ(|)(|)IIЦIIентI>I H napaMeTpti (13). EapiiaHT I, Mog7e.n1> [2p]
Table 3
(Coefficients and parameters (133). Option I, model [2p]
i D, г D,
1 1,73 1151069517248 6 -0,12488741255639
2 -1,924558361 371886 7 869,949176054
3 0,970720385332132 8 -2136,78780101
4 -1,827836692473682 9 3948,40761253
5 -48,5649922569975- 10 —626,224-0064601
Таблица 4
Коэффициенты и параметры (13). Вариант II, модель [1 - а]
Table 4
Coefficients and parameters (13). Option II, model [1 - а]
i D, i D,
1 1,731251069517248 6 -0,1248741255639
2 ^,924558361371886 7 1048,08344232
3 -0,970720385332132 8 -2661,46796047
4 -1,827836692473682 9 2481,26263996
5 ^8,56499225699752 10 -810,761915946
сывает эти данные с несколько меньшей точностью (но в пределах их экспериментальной неопределенности).
Выводы
Предлагаемая модель линии фазового равновесия позволяет: во-первых, описать опытные данные о плотности р - , р+ и давлении р, R236ea в пределах их экспериментальной неопределенности (рис. 1-3), во-вторых, обеспечить согласованность значений р-г, р+и р™,
рассчитанных по предлагаемой методике (уравнения (4), (7), (12), (13)) в диапазоне от тройной точки до критической точки, в-третьих, провести анализ различных моделей среднего диаметра и выбрать оптимальный для исследуемого вещества.
Для линии фазового равновесия холодильного агента R236ea мы рекомендуем использовать модель среднего диаметра [2р], как обеспечивающую описание плотности на линии насыщения с большей точностью, чем модель среднего диаметра [1 - а].
Литература
1. Vorobev V. S., Ochkov V. F., Rykov V. A., Rykov S. V., Ustyuzhanin E. E., Pokholchenko V. A. Development of combined scaling models for liquid and gas densities at the saturation line: Structures and numerical data for SF6 // J. Phys.: Conf. Ser. 2019. V. 1147. P. 012017.
2. Rykov S. V., Kudryavtseva I. V., Rykov V. A., Ustyuzhanin E. E., Ochkov V. F. Analysis of the saturation line on the basis of Clapeyron — Clausius and Gibbs — Duhem equations // J. Phys.: Conf. Ser. 2019. V. 1147. P. 012017.
3. Apfelbaum E. M., Vorob'ev V. S. The Similarity Relations Set on the Basis of Symmetrization of the Liquid-Vapor Phase Diagram // J. Phys. Chem. B. 2015. V. 119. No 26. P. 8419-8424.
4. Apfelbaum E. M, Vorob 'ev V. S. The Wide-Range Method to Construct the Entire Coexistence Liquid — Gas Curve and to Determine the Critical Parameters of Metals // J. Phys. Chem.
B. 2015. V. 119. No 35. P. 11825-11832.
5. Шишаков В. В. Комбинированные скейлинговые модели для инженерных расчетов термодинамических свойств на линии насыщения // Дис. на соискание уч. ст. канд. техн. наук. М.: МЭИ, 2014. 230 с.
6. Рыков С. В., Кудрявцева И. В., Рыков В. А. Методика расчета линии фазового равновесия хладагентов от тройной до критической точки // Холодильная техника. 2017. № 3.
C. 26-30.
7. Ustyuzhanin E. E., Shishakov V. V., Abdulagatov I. M., Popov P. V., Rykov V. A., Frenkel M. L. Scaling models of thermodynamic properties on the coexistence curve: problems and some solutions // Russ. J. Phys. Chem. B. 2012. V. 6. No 8. P. 912-931.
8. Ochkov V. F., Rykov V. A., Rykov S. V., Ustyuzhanin E. E., Znamensky B. E. Extrapolation of IAPWS-IF97 data: The liquid and gas densities on the saturation line near the critical point of H2O // J. Phys.: Conf. Ser. 2018. V. 946. P. 012119.
9. Cailletet L., Mathias E. Recherches sur les densités des gaz liquéfiés et de leurs vapeurs saturées // J. Phys. Theor. Appl. 1886. V. 5 (1). P. 549-564.
10. Goldstein R., Arola A. Liquid-vapor asymmetry at the critical point // Acc. Chem. Res. 1989. V. 22. P. 77-82.
11. Ма Ш. Современная теория критических явлений. M.: Мир. 1980. 298 c.
12. WidomB., Rowlinson J. S. New Model for the Study of Liquid — Vapor Phase Transitions // J. Chem. Phys. 1970 V. 52. P. 1670.
13. Hemmer P. C., Stell G. Fluids with Several Phase Transitions // Phys. Rev. Lett. 1970. V. 24. P. 1284.
14. Рыков В. А. Анализ закономерностей изменения термодинамических свойств веществ в широком диапазоне параметров состояния, включая окрестность критической точки и метастабильную область // Дис. на соискание уч. ст. канд. техн. наук. Л.: ЛТИХП,1988. 275 с.
15. Алтунин В. В. Теплофизические свойства двуокиси углерода. М.: Изд-во стандартов, 1975. 546 с.
16. Козлов А. Д., Лысенков В. Ф., Попов П. В., Рыков В. А. Единое неаналитическое уравнение состояния хладона 218 // ИФЖ. 1992. Т. 62. № 6. С. 840-847.
17. Рыков С. В. Метод построения асимметричного масштабного уравнения состояния в физических переменных // Дис. на соискание уч. ст. канд. техн. наук. СПб.: СПбГУ-НиПТ. 2009. 198 с.
18. Полторацкий М. И. Метод построения фундаментального уравнения состояния и термодинамические таблицы гек-
References
1. Vorobev V. S., Ochkov V. F., Rykov V. A., Rykov S. V., Ustyuzhanin E. E., Pokholchenko V. A. Development of combined scaling models for liquid and gas densities at the saturation line: Structures and numerical data for SF6. J. Phys.: Conf. Ser. 2019. V. 1147. P. 012017.
2. Rykov S. V., Kudryavtseva I. V., Rykov V. A., Ustyuzhanin E. E., Ochkov V. F. Analysis of the saturation line on the basis of Clapeyron — Clausius and Gibbs — Duhem equations. J. Phys.: Conf. Ser. 2019. V. 1147. P. 012017.
3. Apfelbaum E. M., Vorob'ev V. S. The Similarity Relations Set on the Basis of Symmetrization of the Liquid-Vapor Phase Diagram. J. Phys. Chem. B. 2015. V. 119. No 26. P. 8419-8424.
4. Apfelbaum E. M., Vorob'ev V. S. The Wide-Range Method to Construct the Entire Coexistence Liquid — Gas Curve and to Determine the Critical Parameters of Metals. J. Phys. Chem. B. 2015. V. 119. No 35. P. 11825-11832.
5. Shishakov V. V. Combined scaling models for engineering calculations of thermodynamic properties on the saturation curve. Moscow. 2014. 230 pp. (in Russian)
6. Rykov S. V., Kudryavtseva I. V., Rykov V. A. Method of calculating the phase equilibrium line of refrigerants from the triple to the critical point. Kholodilnaya Tekhnika. 2017. No 3. P. 2630. (in Russian)
7. Ustyuzhanin E. E., Shishakov V. V., Abdulagatov I. M., Popov P. V., Rykov V. A., Frenkel M. L. Scaling models of thermodynamic properties on the coexistence curve: problems and some solutions. Russ. J. Phys. Chem. B. 2012. V. 6. No 8. P. 912931.
8. Ochkov V. F., Rykov V. A., Rykov S. V., Ustyuzhanin E. E., Znamensky B. E. Extrapolation of IAPWS-IF97 data: The liquid and gas densities on the saturation line near the critical point of H2O. J. Phys.: Conf. Ser. 2018. V. 946. P. 012119.
9. Cailletet L., Mathias E. Recherches sur les densités des gaz liquéfiés et de leurs vapeurs saturées. J. Phys. Theor. Appl. 1886. V. 5 (1). P. 549-564.
10. Goldstein R., Arola A. Liquid-vapor asymmetry at the critical point. Acc. Chem. Res. 1989. V. 22. P. 77-82.
11. Ma Sh. Modern Theory of Critical Phenomena. New York, NY: Roudedge. 1976.
12. Widom B., Rowlinson J. S. New Model for the Study of Liquid — Vapor Phase Transitions. J. Chem. Phys. 1970. V. 52. P. 1670.
13. Hemmer P. C., Stell G. Fluids with Several Phase Transitions. Phys. Rev. Lett. 1970. V. 24. P. 1284.
14. Rykov V. A. Analysis of regularities of change of substance thermodynamic properties in a wide range of state parameters including the critical point neighborhood and metastable region. Leningrad. 1988. 275 pp. (in Russian)
15. Altunin V. V. Teplofizicheskiye svoystva dvuokisi ugleroda (Thermophysical properties of carbon dioxide). Moscow: Standards Publishing. 1975. 546 pp. (in Russian)
16. Kozlov A. D., Lysenkov V. F., Popov P. V., Rykov V. A. Unified non-analytical equation of state of freon 218. J. Eng. Phys. Ther-mophys. 1992. V. 62. P. 611-617. (in Russian)
17. Rykov S. V. A method for constructing an asymmetric scale equation of state in physical variables. 2009. C. Sc. thesis (Saint-Petersburg: St. Petersburg State University of Refrigeration and Food Engineering). 198 pp. (in Russian)
18. Poltoratskiy M. I. The method of constructing the fundamental equation of state and the thermodynamic tables of hexafluoro-
сафторпропана (R236ea) // Дис. на соискание уч. ст. канд. техн. наук. СПб.: Университет ИТМО, 2018. -165 с.
19. Duschek W., Kleinrahm R., Wagner W. Measurement and correlation of the (pressure, density, temperature) relation of carbon dioxide. II. Saturated-liquid and saturated-vapor densities and vapor pressure along the entire coexistence curve // J. Chem. Thermodyn. 1990. V. 22. P. 841-864.
20. Кудрявцева И. В., Камоцкий В. И., Рыков С. В., Рыков В. А. Расчет линии фазового равновесия диоксида углерода // Научный журнал НИУ ИТМО. Серия. Процессы и аппараты пищевых производств. 2013. № 4. С. 11.
21. Fisher M. E., Orkoulas G. The Yang-Yang anomaly in fluid criticality: experiment and scaling theory // Phys. Rev. Lett. 2000. V. 85. No 4-24. P. 696-699.
22. AnisimovM. A. Universality versus nonuniversality in asymmetric fluid criticality // Condens. Matter Phys. 2013. V. 16. No 2. P. 23603: 1-10.
23. Станкус С. В., Хайрулин Р. А., Мартынец В. Г., Молодое Ю. И. Плотность перфторгексана в окрестности критической точки испарения // Вестник НГУ Серия: Физика. 2013. Т. 8. № 1. С. 73-77.
24. Polikhronidi N. G., Abdulagatov I. M, Batyrova R. G., Stepa-nov G. V., Ustuzhanin E. E., Wu J. T. Experimental study of the thermodynamic properties of diethyl ether (DEE) at saturation // Int. J. Thermophys. 2011. V. 32. P. 559-595.
25. Ustjuzhanin E. E., Ochkov V. F., Znamensky V E., Shishakov V V., Znamensky V. E., Rykov S. V. Investigation of gas and liquid densities on the saturation line: some scaling models and numerical data on H2O example // J. Phys.: Conf. Ser. 2017. V. 891. P. 012346.
26. Рыков С. В., Кудрявцева И. В., Рыков В. А. Анализ линии фазового равновесия, разработанной на основе уравнения Клапейрона-Клаузиуса и различных моделей среднего диаметра // Вестник Международной академии холода. 2018. № 4. С. 28-34.
27. Gruzdev V. A., Khairulin R. A., Komarov S. G., Stankus S. V. Thermodynamic Properties of HFC-236ea // Int. J. Thermophys. 2008. V. 29. P. 546-556.
28. Defibaugh D. R., Gillis K. A., Moldover M. R., Schmidt J. W., Weber L. A. Thermodynamic properties of CF3-CF-CHF2, 1,1,1, 2,3,3-hexafluoropropane // Fluid Phase Equilib. 1996. V. 122. No 1-2. P. 131-155.
29. Zhang H., Sato H., Watanabe K. Vapor Pressure of 1,1,1,2,3,3-hexafluorpropane (R-236ea) from 300 to 410 K // J. Chem. Eng. Data. 1995. V. 20. P. 1281-1284.
30. Di Nicola G, Giuliani G. Vapor Pressure and PVT Measurements for 1,1,1,2,3,3-hexafluorpropane (R-236ea) // J. Chem. Eng. Data. 2000. V. 45. P. 1075-1079.
31. Bobbo S., Fedele L., Scattolini M., Camporese R. Vapor+Liquid Equilibrium Measurements and Correlation of the Binary Refrigerant Mixtures Difluoromethane (HFC-32) +1,1,1,2,3,3-Hexa-fluoropropane (HFC-236ea) and Pentafluoroethane (HFC-125) +1,1,1,2,3,3-Hexafluoropropane (HFC-236ea) at 288.6, 303.2, and 318.2 // Int. J. Thermophys. 2000. V. 21. P. 781-790.
32. Bobbo S., Fedele L., Camporese R., Stryjek R. Vapor-liquid equilibrium for the three binary systems 1,1,1,2,3,3-hexafluoro-propane with dimethyl ether or propane, and 1,1,1,3,3,3-hexafluoro-propane with dimethyl ether // Fluid Phase Equilib. 2000. V. 174. P. 3-12.
33. Lemmon E. W., Huber M. L., McLinden M. O. Reference Fluid Thermodynamic and Transport Properties (REFPROP), version 9.0// National Institute of Standards and Technology, 2010.
propane (R236ea). 2018. C. Sc. thesis (Saint-Petersburg: ITMO University). 165 pp. (in Russian)
19. Duschek W., Kleinrahm R., Wagner W. Measurement and correlation of the (pressure, density, temperature) relation of carbon dioxide. II. Saturated-liquid and saturated-vapor densities and vapor pressure along the entire coexistence curve. J. Chem. Thermodyn. 1990. V. 22. P. 841-864.
20. Kudryavtseva I. V., Kamotskii V. I., Rykov S. V., Rykov V. A. Calculation carbon dioxide line of phase equilibrium. Scientific Journal NRU ITMO. Processes and Food Production Equipment. 2013. No 4. P. 11. (in Russian)
21. Fisher M. E., Orkoulas G. The Yang-Yang anomaly in fluid criticality: experiment and scaling theoryro Phys. Rev. Lett. 2000. V. 85. No 4-24. P. 696-699.
22. Anisimov M. A. Universality versus nonuniversality in asymmetric fluid criticality. Condens. Matter Phys. 2013. V. 16. No 2. P. 23603: 1-10.
23. Stankus S. V., Khairulin R. A., Martynets V. G., Moloro-dov Yu. I. Density of Perfluorohexane near the Evaporation Critical Point. VestnikNSU. Series: Physics. 2013. V. 8. No 1. P. 73-77. (in Russian)
24. Polikhronidi N. G., Abdulagatov I. M., Batyrova R. G., Stepa-nov G. V., Ustuzhanin E. E., Wu J. T. Experimental study of the thermodynamic properties of diethyl ether (DEE) at saturation. Int. J. Thermophys. 2011. V. 32. P. 559-595.
25. Ustjuzhanin E. E., Ochkov V. F., Znamensky V. E., Shishakov V. V., Znamensky V. E., Rykov S. V. Investigation of gas and liquid densities on the saturation line: some scaling models and numerical data on H2O example. J. Phys.: Conf. Ser. 2017. V. 891. P. 012346.
26. Rykov S. V., Kudryavtseva I. V., Rykov V. A. Analysis of the phase equilibrium line developed on the basis of the Clapeyron-Clau-sius equation and various models of average diameter. Journal of International Academy of Refrigeration. 2018. No 4. P. 28-34. (in Russian)
27. Gruzdev V. A., Khairulin R. A., Komarov S. G., Stankus S. V. Thermodynamic Properties of HFC-236ea. Int. J. Ther-mophys. 2008. V. 29. P. 546-556.
28. Defibaugh D. R., Gillis K. A., Moldover M. R., Schmidt J. W., Weber L. A. Thermodynamic properties of CF3-CF-CHF2, 1,1,1,2,3,3-hexafluoropropane. Fluid Phase Equilib. 1996. V. 122. No 1-2. P. 131-155.
29. Zhang H., Sato H., Watanabe K. Vapor Pressure of 1,1,1,2,3,3-hex-afluorpropane (R-236ea) from 300 to 410 K. J. Chem. Eng. Data. 1995. V. 20. P. 1281-1284.
30. Di Nicola G., Giuliani G. Vapor Pressure and PVT Measurements for 1,1,1,2,3,3-hexafluorpropane (R-236ea). J. Chem. Eng. Data. 2000. V. 45. P. 1075-1079.
31. Bobbo S., Fedele L., Scattolini M., Camporese R. Vapor+Liquid Equilibrium Measurements and Correlation of the Binary Refrigerant Mixtures Difluoromethane (HFC-32) +1,1,1,2,3,3-Hex-afluoropropane (HFC-236ea) and Pentafluoroethane (HFC-125) +1,1,1,2,3,3 -Hexafluoropropane (HFC-236ea) at 288.6, 303.2, and 318.2. Int. J. Thermophys. 2000. V. 21. P. 781-790.
32. Bobbo S., Fedele L., Camporese R., Stryjek R. Vapor-liquid equilibrium for the three binary systems 1,1,1,2,3,3-hexafluoropropane with dimethyl ether or propane, and 1,1,1,3,3,3-hexafluoropropane with dimethyl ether. Fluid Phase Equilib. 2000. V. 174. P. 3-12.
33. Lemmon E. W., Huber M. L, McLinden M. O Reference Fluid Thermodynamic and Transport Properties (REFPROP), version 9.0. National Institute of Standards and Technology, 2010.
34. Aoyama H., Kishizawa G., Sato H., Watanabe K. Liquid Coexistence Curves in the Critical Region and the Critical Temperatures and Densities of 1,1,1,2-Tetrafluoroethane (R-134a), 1,1,1-Trifluoroethane (R-143a), and 1,1,1,2,3,3-Hexafluoropropane (R-236ea) // J. Chem. Eng. Data. 1996. V. 41. No. 5. P. 1046-1051.
35. Груздев В. А., Комаров С. Г. Экспериментальное исследование давления и плотности насыщенного и перегретого пара фреона R-236ea от 20 до 150 °С // Теплофизика и аэромеханика. 2006. Т. 13. № 3. С. 443-451.
34. Aoyama H., Kishizawa G., Sato H., Watanabe K. Liquid Coexistence Curves in the Critical Region and the Critical Temperatures and Densities of 1,1,1,2-Tetrafluoroethane (R-134a), 1,1,1-Trifluoroethane (R-143a), and 1,1,1,2,3,3-Hexafluoropropane (R-236eaj. J. Chem. Eng. Data. 1996. V. 41. No. 5. P. 1046-1051.
35. Gruzdev V. A., Komarov S. G. Experimental study of pressure and density of saturated and superheated vapor of refrigerant R-236ea from 20 to 150°C. Thermophys. Aeromech. 2006. V. 13. No 3. P. 411-418. (in Russian)
Сведения об авторах
Рыков Сергей Владимирович
к. т. н., доцент факультета низкотемпературной энергетики Университета ИТМО, 191002, Санкт-Петербург, ул. Ломоносова, 9, togg1@yandex.ru
Кудрявцева Ирина Владимировна
к. т. н., доцент факультета систем управления и робототехники Университета ИТМО, 197101, Санкт-Петербург, Кронверкский пр., 49, togg1@yandex.ru
Рыков Владимир Алексеевич
д. т. н., профессор, доцент факультета низкотемпературной энергетики Университета ИТМО, 191002, Санкт-Петербург, ул. Ломоносова, 9, togg1@yandex.ru
Свердлов Александр Викторович
магистрант факультета низкотемпературной энергетики Университета ИТМО, 191002, Санкт-Петербург, ул. Ломоносова, 9, togg1@yandex.ru
Нурышева Меруерт
магистрант факультета низкотемпературной энергетики Университета ИТМО, 191002, Санкт-Петербург, ул. Ломоносова, 9, togg1@yandex.ru
Information about authors
Rykov Sergey Vladimirovich
Ph. D., associate professor of Faculty of Cryogenic Engineering of ITMO University, 191002, Russia, St. Petersburg, Lomonosov str., 9, togg1@yandex.ru
Kudryavtseva Irina Vladimirovna
Ph. D., associate professor of department of Faculty of Control Systems and Robotics of ITMO University, 49 Kronverksky Pr., St. Petersburg, 197101 Russia, togg1@yandex.ru
Rykov Vladimir Alekseyevich
D. Sc., professor, associate professor of Faculty of Cryogenic Engineering of ITMO University, 191002, Russia, St. Petersburg, Lomonosov str., 9, togg1@yandex.ru
Sverdlov Aleksandr Viktorovich
Undergraduate of Faculty of Cryogenic Engineering of ITMO University, 191002, Russia, St. Petersburg, Lomonosov str., 9, togg1@yandex.ru
Nurysheva Meruert
Undergraduate of Faculty of Cryogenic Engineering of ITMO University, 191002, Russia, St. Petersburg, Lomonosov str., 9, togg1@yandex.ru
РОД экспо
/
27-я международная выставка продуктов питания, напитков и сырья для их производства
/ ^0-14 февраля 2020
«Продэкспо» - самая крупная выставка России по размеру выставочных профессиональному интересу и охвату рынка, по данным Общероссийского рейтинга. «Продэкспо» представляет все отрасли продовольственного рынка.
Организатор выставки:
АО «Экспоцентр»
123100, Москва, Краснопресненская наб., 14 E-mail: mezvist@expocentr.ru www.prod-expo.ru
площадей,
