Например, Бобцов

К вопросу описания термодинамической поверхности, включая критическую область, уравнениями состояния в физических переменных

УДК 536.71
К вопросу описания термодинамической поверхности, включая критическую область, уравнениями состояния в физических пере­
менных
Рыков А.В., Кудрявцева И.В., Рыков В.А.
Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики
В работе проведен подробный анализ одного из современных подходов к по­ строению масштабного уравнения состояния вещества. Отмечены положи­ тельные и отрицательные стороны такого подхода. Проведено сравнение с масштабным уравнением авторов. Ключевые слова: масштабное уравнение, линия псевдокритических точек, крити­ ческие индексы.
До настоящего времени является актуальной задача построения в физических переменных уравнения состояния жидкости и газа, которое с одинаковой точно­ стью описывало бы как регулярную часть термодинамической поверхности, так и окрестность критической точки [1–6]. И это несмотря на то, что решению этой проблемы – построению единого уравнения состояния, описывающего всю термо­ динамическую поверхность системы жидкость–газа с учетом требований совре­ менной теории критичеких явлений за последние полвека посвящено большое число работ, среди них [7–10]. Особенно непростым является выбор нерегулярной составляющей единого уравнения состояния в физических переменных. В работе [6] эта задача решается в переменных плотность-температура на основе предло­ женного авторами неаналитического уравнения состояния

(1)

где p – давление; pc – критическое давление; p = pclas ( ω ,t ) – аналитическое
уравнение состояния; p = pscal ( ω ,t ) – масштабное уравнение состояния; ω = ρ / ρ c
и t = T / Tc – приведенные плотность и температура, соответственно; ρ c – критиче­ ская плотность; Tc – критическая температура; ∆ ρ = ω − 1; τ = t − 1; λ и µ – посто­ янные.

Масштабное уравнение

в (1) имеет следующий вид [6]:

1

(2)
. Неаналитическое уравнение (2) получено из предложенного в [6] уравнения состояния
(3)
путем вычисления определенного интеграла, входящего (3), в приближении
. (4)
Из (2) следует, что уже первая частная производная ( ∂ p ∂ T ) ρ расходится в
каждой точке критической изотермы. Однако авторы [6] справедливо указывают, что выражение (1) «не пригодно в пределе τ → 0 , ∆ ρ → − 1». С нашей точки зре­ ния выражения (1) и (3) некорректны в более широкой области, т.е. если τ → 0 , ∆ ρ № 0 , т. е., по крайней мере, практически в любой точке критической изотермы за исключением критической точки.
В целом структура уравнения (2) представляется довольно сложной. Дей­ ствительно, вычисляя интеграл (алгоритм вычисления интеграла подробно изло­ жен в [5]), входящий в уравнение состояния (2), получим выражение, трудоемкое для практического применения:
2

. Наиболее простое уравнение состояния, в основном, удовлетворяющее тре­ бованиям МТ, согласно (5), имеет вид:
. С другой стороны, в работах [7–10] для асимптотической окрестности крити­ ческой точки предложено масштабное уравнение состояния в физических пере­ менных:
. (5)
Где F ( ρ ,T ) – свободная энергия Гельмгольца; F0 ( T ) и A0 ( T ) – аналитиче­ ские функции температуры, а масштабная функция свободной энергии a0 ( x) ,
предложенная Рыковым В.А., определяется выражением:
. (6)
Как показано в работах Рыкова В.А. [8–12], такой выбор масштабной функ­ ции свободной энергии позволяет удовлетворить всем требованиям масштабной теории и качественно верно описать область метастабильных состояний – вос­ произвести на термодинамической поверхности термическую спинодаль в соот­ ветствии с равенством:
.
Казалось бы, масштабная функция (6) имеет много неопределенных парамет­ ров. Если считать, что критические индексы суть универсальные величины, то их
3

семь: четыре линейных и три нелинейных. На самом деле их только два и те ли­ нейные. Убедимся в этом.
Как известно, масштабные функции химического потенциала h(x) , изохор­ ной теплоемкости f (x) и изотермической сжимаемости fz (x) связаны с масштаб­ ной функцией a0(x) равенствами:
, (7)
, (8)
. (9) Согласно формуле (8), для того, чтобы обеспечить выполнение равенства хи­ мических потенциалов в каждой точке линии насыщения x = − x0 достаточно по­ требовать, чтобы постоянная C0 находилась из равенства:
. На основе совместного анализа масштабных функции
, (10)
(11) рассчитанных на основе (6)–(9), и соответствующих им масштабных функций ли­ нейной модели (ЛМ):
(12)
4

, (13)
. (14) В формулах (12)–(14) масштабная переменная определяется равенством
получим систему уравнений, устанавливающую связь между параметрами масштабных функций в физических переменных и линейной модели:
(15)
(16)
, (17)
, (18)
( )где ε = φ 1 φ 2 ; φ i = xi / x0 (i =1 или 2); φ i0 = φ i − 1; k1 = b2 − 1 β .
Параметры φ 1, φ 2 , φ 3 являются универсальными с точностью до универсаль­ ности критических индексов. Следовательно, для данного набора критических ин­ дексов достаточно один раз решить систему уравнений (19), найти значения x1 , x2 и x3 по формуле
5

, (19)

и

из (3.49),

(3.59) определить

A1 ,

A2

и

B3 . В формуле (3.62)

φ

* i

– решения

системы

(3.61) для заданного набора критических индексов.

В силу универсальности x1 , x2 и x3 уравнение

, (20)
описывает линию особых точек изохорной теплоемкости, т.е. (20) – это уравнение линии псевдокритических точек.
Разработанная на основе анализа степенных функционалов и метода псевдо­ критических точек масштабная функция свободной энергии (6) в физических переменных не уступает в точности описания однофазной области параметриче­ ским масштабным функциям (12–14) уравнения ЛМ.
Согласно масштабной теории критических явлений, поведение изотермиче­
ской сжимаемости KT = ρ − 1 ( ∂ ρ / ∂ p) T определяется следующим выражением. Найдем для функции (2) значение частной производной ( ∂ p / ∂ ρ ) T :

. На линии τ = − qp ∆ ρ 1/β , которую авторы [6] позиционируют как линию син­ гулярности изохорной теплоемкости, из (20) получим
. (21) Но автором работ [13,14] впервые показано, что если на некоторой линии расходится изохорная теплоемкость, то в каждой точке этой линии, за исключени­ ем критической точки выполняются равенства:
(22) 6

Выражение (21) ни одному равенству (22) не удовлетворяет, следовательно линия τ = − qp ∆ ρ 1/β не является линией сингулярности изохорной теплоемкости.
В заключении заметим, что не смотря на то, что масштабное уравнение со­ стояния в физических переменных (3) имеет сложную структуру, предложенный авторами подход представляется интересным и перспективным.
Список литературы:
1. Рыков А.В., Кудрявцева И.В., Рыков В.А. Асимметричное масштабное уравнение состояния хладона R23 // Вестник Международной академии холода. – 2012. – № 4.
2. Рыков С.В., Кудрявцева И.В., Демина Л.Ю. Единое уравнение состояния R717, учитывающее особенности критической области // Вестник Международ­ ной академии холода. – 2009. – № 4. – С. 29–32.
3. Кудрявцева И.В., Рыков В.А., Рыков С.В. Асимметричное единое уравне­ ние состояния R134a // Вестник Международной академии холода. 2008 г. Выпуск № 2. С. 36–39.
4. Рыков С.В., Багаутдинова А.Ш., Кудрявцева И.В., Рыков В.А. Асиммет­ ричное масштабное уравнение состояния // Вестник Международной академии холода. 2008 г. Выпуск № 3. С. 30–32.
5. Безверхий П.П., Мартынец В.Г., Матизен В.Г. Скейлинговское уравнение состояния в реальных переменных для флюидов с учётом асимметрии // Вестник СибГУТИ. 2009. № 3. С. 105–116.
6. Безверхий П.П., Мартынец В.Г., Матизен В.Г. Объединенное уравнение со­ стояния флюидов, включающее регулярную и скейлинговую часть // Сверхкрити­ ческие флюиды: Теория и Практика. Т. 3. № 3. 2008. С. 13–29.
7. Рыков В.А. Анализ масштабного уравнения состояния, основанного на ги­ потезе «псевдоспинодальной» кривой // Журнал физической химии. – 1985. – Т. 59, № 9. – С. 2354–2356.
8. Рыков В.А. Масштабное уравнение состояния в физических переменных // Теплофизика высоких температур. – 1986. – Т.25, № 2. – С. 345.
9. Рыков В.А. Масштабное уравнение состояния, верно воспроизводящее ме­ тастабильную область // Инженерно-физический журнал. – 1985. – Т. 49, № 3. – С. 506–507.
10. Рыков В.А. Масштабные функции свободной энергии Ar, C2 H6, CO2, Xe, N2, O2. // Журнал физической химии. – 1985. – Т. 59, Вып. 3. – С. 792.
11. Рыков В.А. Метод расчета ρ-Т- параметров границы устойчивости одно­ родного состояния вещества // Журнал физической химии. – 1985. – Т. 59, № 8. – С. 2070–2072.
12. Рыков В.А. Методика выбора масштабной функции свободной энергии // Журнал физической химии. – 1984. – Т. 58, № 11. – С. 2852–2853.
13. Рыков В.А. О гипотезе «псевдоспинодальной» кривой // Журнал физиче­ ской химии. –1986. – Т. 60, № 3. – С. 789–793.
7

14. Рыков В.А. Определение «псевдоспинодальной» кривой на основе термо­ динамических равенств (∂T/∂S)v = 0 и (∂V/∂p)T = 0 // Журнал физической химии. – 1985. – Т. 59, № 11. – С. 2905–2906.
15. Лысенков В.Ф., Шустров А.В. О возможности построения масштабного уравнения состояния газа и жидкости в физических переменных // Инженерно-фи­ зический журнал. – 1984. – Т. 47, № 4. – С. 602–608.
To a question of the description of a thermodynamic surface, including critical region, by equations of state in physical variables
Rykov A.V., Kudryavtseva I.V., Rykov V.A.
National Research University of Information Technologies, Mechanics and Optics
In article the detailed analysis of one of the up-to-date approaches to build-up of a scale equation of state of substance is spent. Plus are scored and negative sides of such approach. Comparison with the scale equation of authors is spent. Key words: scale equation, pseudocritical points line, critical index.
8