Научная статья на тему 'Анализ экспериментальной информации о равновесных свойствах r218 на основе неаналитического уравнения состояния'

Анализ экспериментальной информации о равновесных свойствах r218 на основе неаналитического уравнения состояния Текст научной статьи по специальности «Физика»

CC BY
70
26
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УРАВНЕНИЕ СОСТОЯНИЯ / EQUATION OF STATE / ТЕРМОДИНАМИЧЕСКИЕ СВОЙСТВА / THERMODYNAMIC PROPERTIES / R218

Аннотация научной статьи по физике, автор научной работы — Рыков А. В., Кудрявцева И. В., Рыков С. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Рыков А. В., Кудрявцева И. В., Рыков С. В.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

n 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.

Текст научной работы на тему «Анализ экспериментальной информации о равновесных свойствах r218 на основе неаналитического уравнения состояния»

УДК 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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

о & о

Й 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.

i Надоели баннеры? Вы всегда можете отключить рекламу.