CC BY
45
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
УДК 536.71
Особенности расчета термодинамических характеристик сверхкритических флюидов, используемых в скф-технологиях
Канд. техн. наук Рыков С.В., канд. техн. наук Кудрявцева И.В., 165627@niuitmo.ru д-р техн. наук Рыков В.А. rikov.zav.kaf@irbt-itmo.ru Университет ИТМО Институт холода и биотехнологий 921002, Санкт-Петербург, ул. Ломоносова, 9
Дано физическое обоснование новому уравнению состояния, используемому для расчета равновесных характеристик сверхкритических флюидов (СКФ). Особенностью данного метода является, во-первых, возможность использовать его во всей рабочей области, включая область сильно развитых флуктуаций плотности, во-вторых, простоту расчетной схемы. Последнее обстоятельство обусловлено тем, что уравнение состояния, лежащее в основе расчетной схемы, в отличие от известных уравнений состояния, удовлетворяющих масштабной гипотезе, не содержит нелинейных индивидуальных параметров.
Ключевые слова: сверхкритические флюиды, диоксид углерода, уравнение состояния,
кроссоверная функция, изохорная теплоемкость.
Features of calculation of thermodynamic characteristics of the supercritical
fluids used in scf-technologies
Rykov S.V., Kudryavtseva I.V., Rykov V.A.
University ITMO
Institute of Refrigeration and Biotechnologies 191002, Russia, St. Petersburg, Lomonosov str., 9
The physical substantiation is given the new equation of state used for calculation of equilibrium characteristics of supercritical fluids. Feature of the given method is, first, possibility to use it in all working area, including area of strongly developed fluctuations of density, secondly, simplicity of the settlement scheme. Last circumstance is caused by that the equation of a condition underlying the settlement scheme, unlike the known equations of a condition satisfying to a scale hypothesis, does not contain nonlinear individual parametres.
Key words: supercritical fluids, carbon dioxide, equation of state, crossover function, isochoric heat capacity.
В работах [1-8] рассмотрены возможные пути решения проблемы расчета равновесных свойств сверхкритических флюидов в широкой области параметров состояния, в том числе в области сильно развитых флуктуаций плотности. Актуальность данной проблемы обусловлена тем, что в последние три десятилетия в пищевой и косметической промышленности все большее развитие получают технологии, в основе которых лежит использование растворяющей способности сверхкритических флюидов. Эти технологии дают возможность извлекать из растительного сырья вкусовые и ароматические компоненты, которые используются в качестве добавок в конечной продукции [9].
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
Отличительной особенность критической области является существенное, в разы, изменение ряда теплофизических характеристик, таких как изохорная Cv и изобарная C
теплоемкости, изотермическая сжимаемость Кг, что, с одной стороны, делает эту область параметров состояния перспективной для разработки инновационных технологий в пищевой, косметической и фармацевтической промышленности, с другой стороны, предъявляет повышенные требования к методам расчета равновесных свойств рабочих тел в этой области.
Как показано в [1-8] на примере расчета термодинамических характеристик диоксида углерода, нашедшего широкое применение при реализации СКФ-технологии в пищевой и косметической промышленности, перспективным является подход, в основе которого лежит идея использовать в качестве исходной термодинамической функции свободную энергию Гельмгольца F, которая в переменных плотность р и температура T является характеристической функцией [10-15].
К недостаткам подхода, развитого и реализованного в [1-8] в настоящее время относят только то, что в отличие от методов расчета равновесных свойств индивидуальных веществ на основе параметрических уравнений состояния [16], уравнение состояния, разработанные в рамках метода псевдокритических точек [16-22]:
F р,Т =ЯТсщ/ t,со |Др|8+1<т0 х +
i i о)
+F0 T + ДПпю + ДГ(о£ Z CijTl ДР >
7=1 7=0
лежащее в основе метода [1-8], не является физически обоснованным. Последнее утверждение основано на том, что масштабная функция свободной aQ x , входящая в
(1), фактически «сконструирована» таким образом, чтобы удовлетворять масштабной гипотезе.
В формуле (1) использованы обозначения: р - плотность; Т - абсолютная температура; ю = р/рс; Др = со-1; t=T/Tc; т-t-l; Т]_ =1/^ — 1; Тс - критическая температура; рс - критическая плотность; ft, со - кроссоверная функция; а и у -критические индексы изохорной теплоемкости Cv и изотермической сжимаемости Кт, соответственно; [3 и 8 - критические индексы кривой сосуществования и критической изотермы, соответственно; R - газовая постоянная; F0 T - идеально-газовая
составляющая свободной энергии; х = т /1 Др^ - масштабная переменная.
Покажем, что выражение (1) с масштабной функцией ao x в виде [23-26]:
а(х) = А(х + х1 )2~а + В(х + х1 У + С (2)
может быть строго рассчитано на основе феноменологической теории критических явлений [27].
Рассчитаем на основе свободной энергии Гельмгольца (1) и термодинамического равенства S = - 8F/дТ р выражение для энтропии S:
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
S р,Т =-Rutf ю |Ар|
(1-а)/р ,
«з h i
а0 х
щ h 1
-Ft Т - Д1пш-Дш£ £ Cjjx{ Ар ' +ДГ‘юХ £ уС^т/"1 Ар
/'=1 у=0 /'=1 7=0
Выделим в (3) сингулярную составляющую энтропии:
AS р,Т =-|Др|(
(3)
|(1-а)/р ,
6/q Л
где
AS р,Т = S-Sreg / Rf(со)
п3 J3 1 . п3 J3 1 , .
* ■ ^ У1 г
(4)
(5)
р,Г =-F(j Т -Flnco-Fco^ ^ Qyij7 Ар +7F ^ jQj^i Ар • (6)
/'=1 7=0 /=1 7=0
Рассчитаем функцию а'д х на основе феноменологической теории критических
явлений, согласно которой масштабную гипотезу, по аналогии с [28], можно представить в виде уравнения:
1-га
AS- Cv
Представим параметр m в виде:
ф0+ф2-Щ .
m = Ар • Cv
(7)
(8)
Здесь AS = pTcl рс s р,Т ~ sq Р-,Т ; s - энтропия; s0 р ,Т -регулярная функция температуры; рс - критическое давление; (3 - критический индекс кривой
сосуществования; Ар = р/рс -1; ф0 и ф2 - постоянные коэффициенты.
В асимптотической окрестности критической точки имеем:
С„(р,Г) = ^|ДрГа/|3 ДГ + ЛГ! Л (9)
где х = т/| Ар|1р - масштабная переменная; т = Г / Г -1; х1 - параметр линии псевдокритических точек [17], значение которого определяется из уравнения:
TCV (р) ~Тс | Л |1/р
—^---------= ~Х11ЛР| '
Tc
Подставим (8) в уравнение (7) и получим:
(10)
— ^ Р,Т -S0 Р,Т = Cv “ ( ф0 + ф2 Ар " С,
Рс \
2 ^
(11)
Учитывая зависимость (9), найдем
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
s р,Т -s0 Т =Аа |Ар|
PL.
Pc
Из (12) найдем выражение для энтропии:
р,Т = Т|4р| Рс
1-а /р 1-а * у-1 Р Т
JC + JCj +ф2 JC + JCj +I—^s0 1 ,
Pc
(13)
а-1 2(3
где А* = А а ф0; ср*=ф2 -Аа / ср0.
Из выражения (13) следует, что с точностью до постоянных коэффициентов А и В мы приходим, если воспользоваться термодинамическим равенством F = — J*S” р,Т dl', к
масштабной функции свободной энергии Гельмгольца (2).
Таким образом, метод описания равновесных свойств сверхкритических флюидов, предложенный в [1] и апробированный на примере описания термодинамических свойств диоксида углерода, получил физическое обоснование в рамках феноменологической теории критических явлений и гипотезы Бенедека [29]. Следует отметить, что фундаментальное единое уравнение состояния (1) находит широкое применение не только при расчете равновесных свойств сверхкритических флюидов, используемых в СКФ-технологиях, но и для расчета термодинамических характеристик технически важных веществ, используемых в холодильной и криогенной техники [3038]. Дальнейшее совершенствование рассмотренного в работе метода расчета равновесных свойств СКФ связано с введением в структуру уравнения (1) вместо масштабной переменной x «обобщенной» масштабной переменной x [39] и использованием предложенных в [40-51] моделей линии фазового равновесия, удовлетворяющих масштабной теории.
Литература
1. Кудрявцева И.В., Рыков А.В., Рыков В.А. Метод расчета равновесных свойств сверхкритических флюидов, используемых в СКФ-технологиях // Научный журнал НИУ ИТМО. Серия: Процессы и аппараты пищевых производств. 2013. № 2. С. 29.
2. Рыков С.В., Кудрявцева И.В., Рыков В.А. Асимметричное масштабное уравнение состояния аргона в переменных плотность-температура // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2008. № 2. С. 6-11.
3. Кудрявцева И.В., Рыков В.А., Рыков С.В., Селина Е.Г., Курова Л.В. Метод расчета плотности и теплоты парообразования двуокиси углерода // Научный журнал НИУ ИТМО. Серия: Процессы и аппараты пищевых производств. 2013. № 1. С. 25.
4. Кудрявцева И.В., Рыков А.В., Рыков В.А. Новое уравнение для «кажущейся» теплоты парообразования // Научный журнал НИУ ИТМО. Серия: Процессы и аппараты пищевых производств. 2013. № 2. С. 30.
5. Рыков С.В. Выбор структуры масштабных функций асимметричного уравнения состояния // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2009. № 2. С. 1-6.
6. Рыков А.В., Кудрявцева И.В., Рыков В.А. К вопросу описания термодинамической поверхности, включая критическую область, уравнениями состояния в физических
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
переменных // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 1. С. 4.
7. Кудрявцева И.В., Рыков А.В., Рыков С.В. Анализ структуры непараметрического уравнения состояния скейлингового вида // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 2. С. 2.
8. Kozlov A.D., Lysenkov V.F., Popov P.V., Rykov V.A. Single non-analytic equation of R218 chladon state // Инженерно-физический журнал. 1992. Т. 62. № 6. С. 840-847.
9. Алиев А.М., Раджабов Г.К., Степанов Г.В. Компонентный состав экстракта шишкоягод можжевельника продолговатого, полученного экстракцией сверхкритическим диоксидом углерода // Сверхкритические флюиды: Теория и практика. 2012. Т. 7. № 3. С. 20-29.
10. Кудрявцева И.В., Демина Л.Ю. О структуре фундаментального уравнения состояния, учитывающего асимметрию жидкости и пара // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2009. № 1. С. 35-38.
11. Рыков В.А. Масштабное уравнение состояния в физических переменных // Теплофизика высоких температур. 1986. Т. 25. № 2. С. 345.
12. Кудрявцева И.В., Рыков А.В., Рыков В.А. Непараметрическое уравнение состояния скейлингового вида и метод псевдокритических точек // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 2. С. 4.
13. Рыков В.А Уравнение состояния в критической области, построенное в рамках метода нескольких «псевдоспинодальных» кривых // Журнал физической химии. 1985. Т. 59. № 10. С. 2605.
14. Рыков А.В., Кудрявцева И.В., Рыков С.В. Непараметрическое масштабное уравнение состояния, не содержащее дифференциальных биномов // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 2. С. 7.
15. Лысенков В.Ф., Попов П.В., Рыков В.А. Параметрические масштабные уравнения состояния для асимптотической окрестности критической точки. Обзоры по теплофизическим свойствам веществ // ТФЦ - М.: ИВТАН. 1992. № 1 (93). С.3-80.
16. Рыков В.А. Определение «псевдоспинодальной» кривой на основе термодинамических равенств dT/ds v = 0 и dv/dp т = 0 // Журнал физической
химии. 1985. Т. 59. № 11. С. 2905.
17. Рыков В.А. О гипотезе «псевдоспинодальной» кривой // Журнал физической химии. 1986. Т. 60. № 3.
18. Rykov V.A. Method of constructing a single equation of state satisfying the requirements of the scaling hypothesis // Journal of Engineering Physics and Thermophysics. 1985. T. 48. № 4. C. 476-481.
19. Рыков В.А. Масштабное уравнение состояния в р-Г переменных с учетом неасимптотических членов // Журнал физической химии. 1985. Т. 59. № 8. С. 2069.
20. Rykov V.A., Varfolomeeva G.B. Method of determining a structural form of the free energy satisfying the requirements of the scaling hypothesis // Journal of Engineering Physics and Thermophysics. 1985. Т. 48. № 3. С. 341-345.
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
21. Рыков В.А. Структура сингулярных членов свободной энергии, верно воспроизводящих неасимптотические поправки термодинамических функций // Инженерно-физический журнал. 1985. Т. 49. № 6. С. 1027-1033.
22. Rykov V.A. Structure of the singular terms in the free energy correctly reproducing the nonasymptotic corrections to the thermodynamic functions // Journal of Engineering Physics and Thermophysics. 1986. Т. 49. № 6. С. 1502-1508.
23. Лысенков В.Ф., Рыков В.А., Яковлева М.В. Рабочая область асимптотических масштабных уравнений состояния // Теплофизика высоких температур. 1990. Т. 28. № 5. С. 1034.
24. Рыков А.В., Кудрявцев Д.А., Рыков В.А. Метод расчета параметров масштабной функции свободной энергии // Научно-технический вестник Поволжья. 2013. № 5. С. 50-53.
25. Рыков В.А.Уравнение «псевдоспинодальной» кривой // Журнал физической химии. 1985. Т. 59. № 10.
26. Рыков В.А. Анализ масштабного уравнения состояния, основанного на гипотезе «псевдоспинодальной» кривой // Журнал физической химии. 1985. Т. 59. № 9. С. 2354.
27. Мигдал А.А. Уравнение состояния вблизи критической точки // ЖЭТФ. 1972. Т. 62. № 4. С. 1559-1573.
28. Берестов А.Т. Исследование уравнения состояния в широкой окрестности критической точки // Дис. на соискание уч.ст. канд. физ.-мат.. наук. -М.: - 1977. -103 с.
29. Бенедек Дж. Спектроскопия оптического смещения и ее применения к задачам физики, химии, биологии и техники // УФН. 1972. Т. 106. С. 481-504.
30. Носков А.Н., Петухов В.В. Изменение параметров состояния пара хладагента в элементарном рабочем процессе маслозаполненного холодильного винтового компрессора // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2007. № 1. С. 10-14.
31. Докукин В.Н., Емельянов А.Л., Носков А.Н. Результаты испытаний маслозаполненного винтового компрессора малой производительности в высокотемпературных режимах // Вестник Международной академии холода. 2009. № 1. С. 6-8.
32. Носков А.Н., Петухов В.В., Чернов Н.П. Результаты испытаний маслозаполненного винтового компрессора малой производительности // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2007. № 1. С. 15-18.
33. Носков А.Н., Зимков А.А. Расчет процесса сжатия маслозаполненного холодильного винтового компрессора // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 1. С. 3.
34. Носков А.Н., Зимков А.А. Расчет процесса всасывания маслозаполненного холодильного винтового компрессора // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2012. № 1. С. 2.
35. Бараненко А.В., Кириллов В.В., Сивачев А.Е. О выборе хладоносителя для систем косвенного охлаждения // Вестник Международной академии холода. 2010. № 2. С. 22-24.
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
36. Ховалыг Д., Бараненко А.В. Методы расчета градиента давления двухфазного потока при течении в малых каналах // Вестник Международной академии холода. 2012. № 1. С. 3-10.
37. Бараненко А.В., Кириллов В.В. Разработка хладоносителей на основе электролитных водно-пропиленгликолевых растворов // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2007. № 3. С. 38-41.
38. Цветков О.Б. Исследования теплофизических свойств холодильных агентов в санкт-петербургском государственном университете низкотемпературных и пищевых технологий // Вестник Международной академии холода. 2011. № 1. С. 8-9.
39. Рыков С.В. Метод построения асимметричного масштабного уравнения состояния в физических переменных // Диссертация на соискание ученой степени кандидата технических наук / Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий. Санкт-Петербург, 2009.
40. Кудрявцева И.В., Камоцкий В.И., Рыков С.В., Рыков В.А. Расчет линии фазового равновесия диоксида углерода // Научный журнал НИУ ИТМО. Серия: Процессы и аппараты пищевых производств. 2013. № 2. С. 31.
41. Устюжанин Е.Е., Шишаков В.В., Абдулагатов И.М., Попов П.В., Рыков В.А., Френкель М.Л. Скейлинговые модели для описания термодинамических свойств на линии насыщения: проблемы и некоторые решения // Сверхкритические флюиды: Теория и практика. 2012. Т. 7. № 3. С. 30-55.
42. 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 // Russian Journal of Physical Chemistry B. 2012. Vol. 6. N 8. P. 912-931.
43. Устюжанин Е.Е., Абдулагатов И.М., Попов П.В., Шишаков В.В., Рыков В.А. Скейлинговые модели для описания термодинамических свойств на пограничной кривой: характеристики и критерии // Ультразвук и термодинамические свойства вещества. 2008. № 34-35. С. 159-171.
44. Кудрявцева И.В., Рыков А.В., Рыков В.А. Модифицированное уравнение линии насыщения, удовлетворяющее требованиям масштабной теории // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 2. С. 3.
45. Рыков В.А. Термодинамические свойства R23 на линии насыщения в диапазоне температур от 180 до 298 К // Вестник Международной академии холода. 2000. № 4.
46. Кудрявцева И.В., Рыков А.В., Рыков С.В. Описание метастабильной области непараметрическими уравнениями состояния скейлингового вида // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 2. С. 6.
47. Кудрявцева И.В., Рыков С.В., Рыков В.А. Описание метастабильной области в рамках параметрического представления масштабной теории // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 2. С. 5.
48. Кудрявцева И.В., Рыков С.В., Рыков В.А. Непараметрическое уравнение состояния скейлингового вида и расчет равновесных свойств сверхкритических флюидов //
Научный журнал НИУ ИТМО.
Серия «Процессы и аппараты пищевых производств»
№ 1, 2014
Научный журнал НИУ ИТМО. Серия: Процессы и аппараты пищевых производств. 2013. № 2. С. 28.
49. Рыков А.В., Кудрявцева И.В., Рыков В.А. Асимметричное масштабное уравнение состояния хладона R23 // Вестник Международной академии холода. 2012. № 4. С. 26-28.
50. Кудрявцева И.В., Рыков А.В., Рыков В.А., Рыков С.В. Единое неаналитическое уравнение состояния перфторпропана, удовлетворяющее масштабной теории критических явлений // Вестник Международной академии холода. 2013. № 3. С. 22-26.
51. Рыков А.В., Кудрявцева И.В., Рыков С.В. Уравнения линии насыщения и упругости хладона R218 // Вестник Международной академии холода. 2013. № 4. С. 54-57.
