CC BY
23
УДК 536.71
Метод построения фундаментального уравнения состояния, учитывающего особенности критической области
Рыков С.В., Кудрявцева И.В., Рыков А.В., Курова Л.В.
Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики
Полученное в работе уравнение учитывает асимметрию реальной жидкости относительно критической изохоры. Количественная оценка его точности и размера рабочей области проведена на примере описания равновесных свойств хладона R218. Погрешность описания плотности и изо-хорной теплоемкости в однофазной области, теплоты парообразования, давления и плотности на линии фазового равновесия предложенным неаналитическим уравнением состояния сравнима с погрешностью соответствующих экспериментальных данных.
Ключевые слова: R218, асимметричное уравнение состояния, фундаментальное уравнение состояния, равенство химических потенциалов, термодинамические свойства.
Рассмотренный в работе метод является развитием метода построения широкодиапазонного уравнения состояния, предложенного в работе [7]. Уравнения состояния, разработанные в рамках подхода [7], качественно удовлетворяют всем требованиям, которые обычно предъявляют к единым аналитическим уравнениям состояния и обеспечивают, во-первых, предельный переход в области малых плотностей Р и давлений к уравнению состояния идеального газа:
(1)
Во-вторых, равенство химических потенциалов т на обеих ветвях линии насыщения:
(2)
В-третьих, выполнение правила Планка-Гиббса:
В-четвертых, критических условий:
Однако, при описании широкой окрестности критической точки автор [7] ограничился передачей особенностей поведения вещества в соответствии с моделью решеточного газа и не учел асимметрию реальной жидкости относительно критической изохоры. Пути устранения этого недостатка подхода, развитого в [7], и обсуждаются в данной работе.
С этой целью были проанализированы уравнения состояния, в которых учитывалась асимметрия системы жидкость-газ относительно критической изохоры [4, 6]. Обратим внимание, что в соответствии с современной теорией критических явлений [5], равновесные свойства реальной системы жидкость-газ на изолиниях в окрестности критической точки, должны удовлетворять следующим зависимостям:
- на критической изотерме:
1&4Л/Р |Л
4Р1« = ± А, |ЛРГ + А |ЛРГ ± °2 1дРПф + М'
с\ ^ = с® |д рГ"р + с® |Др|1"и+Аур ± С® |Др|
(-а+дур
(5)
на критической изохоре:
■ I „ | | _ | |2-01 | |
л^Ц=о = £?-1М+^М + АМ +^2 Iх
|1+Д,
1Др=0
(6)
ар^ = Н + а Iх Г31 + 1хН™+°2 Iх
|2-а+Д
|1+Ди
1Др=0
- на линии равновесия:
= ^ к I+^ к га+А к г*+°2 к
(7)
Подходы к решению поставленной проблемы можно разбить на две группы: соответственно, основанные на параметрическом подходе [2, 3] и основанные на методе построения масштабных асимметричных уравнений состояния в физических переменных [4, 6, 7]. Область применения уравнений состояния, основанных на параметрической модели, в основном, ограничена широкой окрестностью критической точки. В тоже время, в работах [4, 6] показано, что предложенные в них уравнения состояния ни качественно, ни количественно не уступают параметрическим уравнениям состояния. В силу того, что регулярная
часть термодинамической поверхности описывается и количественно и качественно верно едиными уравнениями состояния в физических переменных, представляется целесообразным для решения обозначенной в данной работе проблемы в дальнейшем использовать физические переменные плотность-температура.
Фундаментальное уравнение состояния, положенное в основу рассматриваемого в данной работе метода, имеет вид:
(8)
где ^ (Т) - идеально-газовая составляющая свободной энергии Гельмгольца;
Т = Т (ш ) - уравнение линии насыщения; аг (х) - масштабные функции свободной энергии [3]:
(9)
(10)
(12)
+ П31(х + х34)у+А2 +С3
Здесь х - обобщенная масштабная переменная Рыкова; а , Р , у , § , А г -критические индексы.
В формуле (8) функции 1%0 (Т) и А§ (Т) - это аналитические функции температуры, которые представлены в виде:
Щ т2
г=3 >0
Функция £ (Т, Т) представлена в виде произведения двух функций:
(13)
Если в термическое уравнение, полученное на основе (8) подставить
Т = Т.
, то в результате получим выражение для давления на линии упругости:
Функция $(Т,Т) выбрана в соответствии с [7]:
(14)
(15)
в которой показатель степени ^ I 2.
Обратим внимание, что частная производная $ ш = (3$/3ш) т обращается в
нуль в каждой точке линии фазового равновесия.
Термическое уравнение состояния с учетом уточнений, сделанных о функциях £ (Т, Т) и $ (Т, Т), принимает следующий вид:
Рс
1-сх+Дз
г=0
И / " I
рГЙ
+
г=3
X
т2 г=1
А
-ДГ^
5ЯЗ
ДГГ
Ш2
■т0и0со1п со
7-; («Я [Чи
г
уио-1'
7
+
(16
«з Л СО
Е Е
сог; (со)
Г
ДГ^
-1
Рс 1=2 >0
В силу того, что т = ^ + р / р , используя формулы (8) и (16), можно показать, что равенство (1) в рамках рассмотренного в работе подхода выполняется строго в каждой точке линии фазового равновесия.
Предложенное в данной работе фундаментальное асимметричное уравнение состояния также верно, как и асимметричное масштабное уравнение состояния [10], передает как требования (1) ё (4), так и зависимости (5) ё (7).
Количественная оценка его точности и размера рабочей области проведена на примере описания равновесных свойств хладона R218. Погрешность описания плотности и изохорной теплоемкости в однофазной области, теплоты парообразования, давления и плотности на линии фазового равновесия предложенным фундаментальным асимметричным уравнением состояния сравнима с погрешностью соответствующих экспериментальных данных. Полученные результаты (см. рис. 1 и рис. 2) позволяют сделать вывод о том, что рабочая область фундаментального асимметричного уравнения состояния (8) ограничена по температуре диапазоном 190^430 ^ а по давлению 50 МПа.
)
8 р, % 1,2 1
0,8 0,6 0,4 0,2 0 -0,2 -0,4
320 340 360 380 400 420 440 Т, К
о
о 1 □ 2 А 3 X 4
X □
X □ о О X Лл □
X □ 1 к
□ А о
Рис. 1. Отклонения значений плотности, рассчитанных по уравнению состояния (16) данной работы, от экспериментальных данных [1]: 1 - 564,6 кг/м3; 2 - 479,9 кг/м3; 3 - 375,3
кг/м3; 4 - 187,8 кг/м3.
S Cv, % 1,5 1
0,5 0
-0,5 -1 -1,5
330 340 350 360 370 380 T, K
Рис.2. Отклонения значений Cv хладона 218, рассчитанных на основе фундаментального асимметричного уравнения состояния (8) данной работы, от опытных данных Рыкова В.А. [6] и Рябушевой Т.И. [9]. Изохоры: 1 - 806 кг/м3; 2 - 503 кг/м3; 3 - 685 кг/м3; 4 - 377
кг/м3; 5 - 995 кг/м3.
Список литературы:
1. Brown I. A. Physical properties of perfluoropropane // J. Chem. Eng. Data. -1963. Vol. 8, № 11. P. 106-108.
2. Sengers J.V., Levelt Sengers J.M.H. A universal representation of the thermodynamic properties of fluids in the critical region // Intern. J. Thermoph. - 1984. -V.5, № 2. - P. 195-208.
3. Адамов Ш.П., Анисимов М.А., Киселев С.Б., Смирнов В.А. Уравнение состояния и методы расчета теплофизических свойств в окрестности критической точки // ИФЖ. - 1981, - Т. 40, № 2, - С. 163-164.
4. Кудрявцева И.В., Рыков В.А., Рыков С.В. Асимметричное единое уравнение состояния R134a // Вестник Международной академии холода. - 2008. -№ 2. - С.36-39.
5. Ма Ш. Современная теория критических явлений. - М.: Мир. - 1980. -298 с.
6. Рыков В.А. Анализ закономерностей изменения термодинамических свойств веществ в широком диапазоне параметров состояния, включая окрестность критической точки и метастабильную область // Дис. на соискание уч. ст. канд. техн. наук. - Л.: ЛТИХП, - 1988. - 275 с.
7. Рыков В.А. Единое неаналитическое уравнение состояния газа и жидкости и таблицы термодинамических свойств аргона и хладагентов R134а, R218, R134а // Дис. на соискание уч.ст. докт. техн. наук. СПб.: СПбГУНиПТ, 2000. -456 с.
8. Рыков С.В., Кудрявцева И.В. Выбор структуры асимметричных масштабных функций свободной энергии в физических переменных // Вестник Международной академии холода. - 2009. - № 1. - С. 43-45.
X
х X о
о
О 1
о □ jid А □ 2 Д3 X 4
€
о 5
9. Рябушева Т.И. Исследование изохорной теплоемкости холодильных агентов // Дис. на соискание уч. ст. канд. техн. наук. - Л.: ЛТИХП, - 1979. -189 с.
10. Рыков А.В., Кудрявцева И.В., Рыков В.А. Асимметричное масштабное уравнение состояния хладона R23 // Вестник Международной академии холода. - 2012. - № 4. - С. 26-28.
Method of construction of the fundamental equation of state considering features of critical area
Rykov S.V., Kudryavtseva I.V., Rykov A.V., Kurova L.V.
National Research University of Information Technologies, Mechanics and Optics
The equation received in work considers asymmetry of a real liquid concerning critical isochoruses. The quantitative estimation of its accuracy and the size of working area is spent on an example of the description of equilibrium properties R218. The error of the description of density and isochoric thermal capacities in single-phase area, warmth of steam formation, pressure and density on a line of phase balance the offered not analytical equation of state is comparable with a margin error corresponding experimental data.
Key words: R218, asymmetrical equation of state, fundamental equation of state, equality of the chemical potentials, thermodynamic properties.
