CC BY
23
УДК530.17:536.76
Методы прогнозирования псевдокритических параметров растворов хладагент/масло
Д-р техн. наук В. П. ЖЕЛЕЗНЫЙ, канд. техн. наук Ю. В. СЕМЕНЮК, Т. Л. ЛОЗОВСКИЙ, А. С. МАРКВАРТ
Одесская государственная академия холода 65082, Украина, г. Одесса, ул. Дворянская 1/3
In present paper a new method for the prediction of the pseudocritical parameters for the refrigerant/oil solution (ROS) is reported. The method is based on using limited experimental information. In proposed models for determination of the pseudocritical parameters of the ROS a few-constant scaling equations with universally changing indexes have been applied. The prediction results for the pseudocritical parameters of the ROS R600a/Azmol, R600a/Reniso WF 15A and R245fa/Planetelf ACD 100FY are presented.
Key words: solution refrigerant/oil, pseudocritical temperature, scaling equation, capillary constant, surface tension, saturated vapor pressure.
Ключевые слова: раствор хладагент/масло, псевдокритическая температура, уравнения скейлинга, капиллярная постоянная, поверхностное натяжение, давление насыщенных паров.
Введение
Информация о критических и псевдокритических параметрах веществ и растворов необходима для моделирования свойств сложных термодинамических систем. Используя критические параметры, можно выразить соотношения между различными свойствами в инвариантной форме и придать им универсальный характер.
Однако, вопросы прогнозирования критических параметров растворов хладагент/масло (РХМ) в настоящее время недостаточно изучены.
Компрессорные масла являются сложными или, как предпочитают называть некоторые авторы, «нетрадиционными» термодинамическими системами. Под этим термином следует понимать объекты исследования, для которых отсутствует информация о составе, псевдокритических параметрах, средней молекулярной массе.
Очевидно, что широко применяемые структурноаддитивные методы расчета [I] из-за неопределенности состава компрессорных масел не могут быть использованы для прогнозирования псевдокритических параметров РХМ. Единственным достоверным источником информации о псевдокритических параметрах таких систем могут являться данные, полученные импульсным методом [2]. По этим причинам информация о псевдокритических параметрах сложных термодинамических систем до сих пор остается очень ограниченной.
Опубликованные в литературе данные о плотности, поверхностном натяжении, фазовых равновесиях систем жидкость—пар и жидкость—жидкость, вязкости, теплоемкости хладагентов и растворов хладагент/масло позволяют приступить к обобщению полученной эмпирической информации [3].
Следует отметить, что к настоящему времени опубликовано лишь несколько работ, в которых обсуждаются способы определения псевдокритических параметров РХМ 12—6].
В настоящей работе для определения псевдокритических параметров РХМ предлагается использовать наиболее доступную информацию (о давлении насыщенных паров, плотности и капиллярной постоянной РХМ)
с привлечением малоконстантных уравнении скеилинга [3,7, 8].
Уравнения, которые используются для аппроксимации экспериментальных данных по термодинамическим свойствам веществ на линии насыщения в широком интервале температур, включая окрестность критической точки, должны отвечать нескольким требованиям:
— незначительное количество параметров, имеющих определенный физический смысл;
— высокие экстраполяционные возможности применяемых корреляций;
— инвариантность значений параметров корреляций относительно диапазона используемой эмпирической информации.
В полной мере указанным требованиям отвечают уравнения скейлинга с универсально изменяющимися эффективными показателями степени, предложенные в работах [3, 7, 9].
Определение псевдокритической плотности РХМ
В основу предлагаемой методики определения псевдокритической плотности РХМ положена пропорциональная зависимость между скейлинговым значением парахора и критическим мольным объемом [9]:
[Р]с = 0,750К,
(1)
Скейлинговое значение парахора является комплексом критических индексов и амплитуд [9]:
Ро
(2)
где М — молекулярная масса;
Ро, ст0, /3, м — критические амплитуды и критические показатели для разности ортобарических плотностей Ар и поверхностного натяжения а, соответственно.
Как показано в работе [9], для расчета скейлингового значения парахора необходима информация о критических индексах и амплитудах. Эту информацию можно
получить при аппроксимации данных о разности плотностей на линиях кипения и конденсации и поверхностного натяжения уравнениями скейлинга с универсально изменяющимися эффективными показателями степени
13, 7, 9|
Ар = р' - р" = Pot0F(t\ а = сго^^1\
(3)
(4)
где £ = (1 — Т/Тс) — приведенная температура;
Р(£) и /(£) — универсальные кроссоверные функции, которые могут быть рассчитаны по формулам, представленным в литературе [7]:
F« = 1-I-2278RI) + 1'32s2.»(t)-
(5)
t1'5 t f{t) = 1 - 0,03534^-7^ - 0,31656
in (t)
+0,34246
in (t)
+
ln(£)'
(6)
Следует отметить, что в предложенных формах ап-проксимационных зависимостей критические показатели, наряду с амплитудами, приобретают индивидуальные для каждого вещества значения. Однако в условиях ограниченной эмпирической информации (в задачах прогнозирования) их можно принять равными универсальным теоретическим величинам /3 = 0,3245, д = 1,2586 (см. 110]).
В настоящей работе скейлинговое значение парахора для РХМ предлагается рассчитывать по правилу аддитивности
[P]c,rtiix — [P]c,refCref “Ь [-^>]с,ог/(1 £ге/),
(7)
^ІГПІХ -- ^^ref^ref "Ь А ( 1 £ге/)
(8)
и выполнив переход к массовым концентрациям хге/, можно определить значение псевдокритической плотности РХМ
0,750Мт,х
Рс —
[Р] с,mi-,
0,750
£ref{&о /Ро)ге/ (1 xref){a0 /po)oii
(9)
Таким образом, в отличие от опубликованных методов определения псевдокритической плотности [4—6], точность предложенного метода прогнозирования не зависит от погрешности экспериментальной информации о свойствах РХМ.
На рис. 1 приведены концентрационные зависимости псевдокритической плотности РХМ, рассчитанные из данных [11 — 13] по описанной выше методике.
хгф кг/кг
где [Р]с,ге/ и [Р]с,оіі — скейлинговые значения парахора для хладагента и масла, соответственно;
£ге; — мольная доля хладагента в растворе.
Таким образом, располагая значением скейлингового парахора [Р]с,шгж из(1), рассчитав среднюю (кажущуюся) молекулярную массу
Рис. 1. Зависимость псевдокритической плотности от концентрации для растворов Я600а/А1то1, ЯбООа/Яепш \VF15A и Я245/а/Р1апе1е/АСО \00FY
Определение псевдокритической температуры РХМ по плотности на линии кипения
В основе предлагаемой методики лежит уравнение для аппроксимации плотности РХМ на линии кипения р„пХ, которое было предложено в работе [7]:
Ьхртгх = 1п (рс{Хге{)) + В[ {хге{)90Р^) , (10)
где 5;(а:Ге/) — критическая амплитуда;
в = 1п (Тс(хге1)/Т);
~~ универсальная для неассоциированных веществ кроссоверная функция
в0,4
m = 1 — 1.Ш-.
(11)
Отметим два важных достоинства формулы (9). Во-первых, для определения псевдокритической плотности РХМ не требуется информация о молекулярной массе масла. Во-вторых, для расчета псевдокритической плотности РХМ достаточно лишь данных о свойствах компонентов раствора.
В предложенном уравнении (10) псевдокритическая температура может определяться как подгоночный параметр при описании данных о плотности РХМ. Выбор этого уравнения продиктован несколькими обстоятельствами. Во-первых, в уравнении содержится всего два подгоночных коэффициента — Bi(xref) и Tc(xref). Критический индекс для растворов определяется по правилу аддитивности. Во-вторых, коэффициент Bi(xref) при xref = const практически не зависит от интервала температур, в котором содержится исходная информация [7]. В-третьих, уравнение позволяет описывать плотность жидкости на линии кипения практически во всем интервале температур существования жидкой фазы [7].
Определение псевдокритической температуры РХМ по капиллярной постоянной
Температурная зависимость капиллярной постоянной является практически прямолинейной. Поэтому данное термодинамическое свойство может быть использовано для прогнозирования псевдокритической температуры термически нестабильных веществ, включая компрессорные масла и РХМ. Другим достоинством применения капиллярной постоянной для прогнозирования псевдокритической температуры является простота ее измерения в эксперименте без использования информации о других термодинамических величинах.
Однако при интерпретации экспериментальных данных и определении псевдокритических параметров РХМ в рамках термодинамического моделирования поверхностных свойств необходимо учитывать одно важное обстоятельство. Анализ результатов выполненных исследований показывает, что адсорбция летучего компонента (хладагента) на поверхности жидкой фазы раствора может привести к существенному изменению состава поверхностного слоя РХМ по сравнению с составом во всем объеме жидкости [1, 11, 14, 15}. Поэтому поверхностные свойства раствора, такие как капиллярная постоянная, поверхностное натяжение и давление насыщенных паров зависят, прежде всего, от эффективной концентрации его поверхностного слоя [1, 3, 11, 14, 15].
Очевидно, что измеренные в эксперименте значения объемных свойств (плотность, вязкость) и поверхностных свойств должны быть отнесены к соответствующим (отличающимся друг от друга) составам раствора. Если состав объемной фазы раствора хге/у может быть надежно определен при экспериментальном исследовании, то вопрос об эффективных значениях концентрации xrefts и толщины поверхностного слоя остается недостаточно изученным.
В работе [3] предлагается методика определения эффективного состава поверхностного слоя, основанная на применении парахорных корреляций. В рамках данной методики для трех экспериментально исследован-
ных РХМ были установлены зависимости между объемными xrefy и поверхностными xrej s концентрациями [11—13]. На рис. 2 показан пример такой зависимости.
Полученная информация о составе поверхностного слоя РХМ позволила произвести коррекцию экспериментальных данных по поверхностным свойствам. Суть коррекции состоит в отнесении экспериментальных значений капиллярной постоянной к эффективному составу поверхностного слоя (вместо концентрации хладагента в жидкой фазе РХМ). Процедура коррекции экспериментальных данных показана на рис. 3. По мнению авторов, именно эти скорректированные экспериментальные данные должны использоваться при термодинамическом моделировании фазовых равновесий и поверхностного натяжения РХМ.
В рамках предлагаемой методики псевдокритическая температура определяется как подгоночный параметр при описании данных по капиллярной постоянной РХМ скейлинговым уравнением [7]
°ynix{xref) = al(xTef)lni,{t), (12)
где а§, п — критические амплитуда и индекс, который для РХМ рассчитывается с использованием правила аддитивности для компонентов;
?/>(£) — универсальная для нормальных жидкостей кроссоверная функция [7, 9]
+ 1,5
ip(t) = 1 - 0,04762-—. (13)
Ini
Поскольку концентрация хладагента в жидкой фазе раствора и эффективная концентрация хладагента в поверхностном слое различны, значения псевдокритических параметров, получаемых при описании объемных свойств (плотность, вязкость) и поверхностных свойств, должны отличаться. Однако при х = idem должно соблюдаться условие равенства псевдокритических параметров, найденных для поверхностного слоя и жидкой фазы раствора.
xrefv, кг/кг
Рис. 2. Изотермы концентрации хладагента в поверхностном слое растворов Я600а/А1то1
Рис. 3. Коррекция экспериментальных данных для капиллярной постоянной растворов КбООа/Агто!
Проведенный анализ показывает, что значения псев-докритических температур, выделенных при описании плотности уравнением (10) и при описании капиллярной постоянной уравнением (12) согласуются в пределах 10—15 К (рис. 4). Этот результат косвенно подтверждает корректность методики определения эффективной концентрации поверхностного слоя [3].
Концентрационные зависимости псевдокритических температур для трех экспериментально исследованных РХМ [ 11 — 131 показаны на рис. 5.
где ад — критерий Риделя;
Ь и с — подгоночные параметры.
Выбор уравнения (14) для прогнозирования псевдо-критического давления РХМ продиктован несколькими обстоятельствами. Во-первых, в этом уравнении содержится всего три параметра: ад, Ь и с, причем показатель с очень мало зависит от химической природы вещества. Поэтому его значение без особого ущерба для точности аппроксимации экспериментальных данных можно принять равным с = 2,64 (см. [8]). Во-вторых, если экспериментальные данные достоверны, то значения коэффициентов уравнения (14) практически не зависят от интервала параметров, в котором содержится исходная информация. В-третьих, уравнение (14) позволяет описывать давление насыщенных паров веществ практически во всем интервале температур существования жидкой фазы.
В предлагаемой методике определения псевдокрити-ческого давления используется псевдокритическая температура, определенная по изложенной выше методике. Критерий Риделя расчитывался по известному в литературе (см. [1]) уравнению
ад = 4,926а; + 5,818, (15)
■ из данных пор' ---- из данных по сг
где ш — критерий Питцера, который расчитывается по формуле
X 800-г 750-
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
хгф кг/кг
Рис. 4. Концентрационная зависимость псевдокритической температуры растворов ЯбООа/Яет&о И7/715А, полученная из экспериментальных данных по плотности и капиллярной постоянной
Хгф кг/кг
Рис. 5. Концентрационная зависимость псевдокритических температур растворов R600a/Azmol, R600a/Reniso WF 15А и R245fa/Planetelf ACD 100FY
Определение псевдокрнтического давления РХМ
Для прогнозирования значений псевдокрнтического давления РХМ предлагается использовать предложенную в работе [8] корреляцию для описания давления насыщенных паров РХМ на линии кипения Ps,mix
lii Ps.mix In (Pc(xrej)) aR(xref)d b(xref)6 , (14)
где Т/Тс = 0,7.
После подстановки формул (15) и (16) в уравнение (14) и выполнения алгебраических преобразований, получаем:
1п Рз,шх = Ь (Ре(хге/)) - ИМ + в2’64) - М, (17)
где кх = 0,593767; к2 = 3,76446.
Проведенное исследование показывает, что критерий Риделя является линейной функцией коэффициента Ь:
ад = к\Ь + &2. (18)
Таким образом, в уравнении (17) остается два неизвестных параметра — Рс и Ь, которые могут быть выделены из ограниченных по объему экспериментальных данных о давлении насыщенных паров РХМ.
Изложенная методика была применена для прогнозирования значений псевдокритического давления трех ранее изученных РХМ (см. [11—13]). При этом предварительно была проведена корректировка экспериментальных данныходавлении насыщенных паров РХМ. Приведенные в работах [11 — 13] значения давления Р^тгх были отнесены к рассчитанным значениям концентрации поверхностного слоя. Зависимости псевдокритического давления от концентрации хладагента в поверхностном слое показаны на рис. 6, а критические кривые изученных растворов хладагент/масло в координатах Р—Т приведены на рис. 7.
■£, 200 160 120 80 40 0
----- R600a/Azmol
-----R600a/Reniso WF 15A
-----R245fa/Planetelf ACD 100FY
---R600a/Azmol
-----R600a/Reniso WF 15A
.-----R245fa/Planetelf ACD 100FY
—I-----------------
0,8 1
xref, 5. Kr/Kr
Рис. 6. Зависимость псевдокритического давления растворов ЯбООа/Аипо!, ЯбООа/Четю № 15А и Я245/а/Р1апе1е1/АСО ЮОРУ от концентрации поверхностного слоя
Выводы
Данные о псевдокритических параметрах растворов хладагент/масло имеют важное значение для моделирования термоди намических свойств этих растворов в условиях ограниченного объема экспериментальной информации, неопределенного состава и термической нестабильности компрессорных масел. В статье приведена новая методика прогнозирования псевдокритических параметров для РХМ. Предложенная модель основана на применении малоконстантных уравнений расширенного скейлинга с универсально изменяющимися эффективными показателями степени для описания плотности жидкой фазы и капиллярной постоянной РХМ. Приведены результаты прогнозирования псевдокритических параметров для растворов R600a/Azmol, R600a/Reniso WF I5A и R245fa/Planetelf ACD 100FY. Достоверность рассчитанных значений псевдокритических параметров РХМ подтверждается как согласованием с результатами экспериментального измерения псевдокритической температуры компрессорных масел импульсными методами [2], так и результатами моделирования фазовых равновесий и поверхностного натяжения изученных РХМ [3, 11-13].
Проведенные исследования показывают, что с использованием полученной информации о псевдокритических параметрах возможно термодинамическое моделирования фазовых равновесий в диапазоне приведенных температур до т = Т/Тс = 0,9. Указанный интервал температур в полной мере охватывает интервалы параметров эксплуатации холодильной техники.
Список литературы
1. Рид Р., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей. — JI.: Химия, 1982.
2. Skripov P. V., Starostin A. A., Volosnikov D. V, Zhelezny V. P. Comparison of thermophysical properties for oil/refrigerant mixtures by use of pulse heating method // Int. J. of Refrigeration. 2003. N 26.
3. Лозовский T. JL, Семенюк Ю. В., Ивченко Д. А., При-ходченко И. А., Железный В. П. Поверхностное натяжение смесевых хладагентов и растворов хладагент/масло. Эксперимент методы прогнозирования. Ч. 3: Поверхностное натяжение растворов хладагент/масло // Холодильная техника и технология. 2009. № 6 (122).
Рис. 7. Критические локусы Р—Т растворов Я600а/А1то1, ЯбООа/Яетво И^ 15А и Я245/а/Р1апе1е1/АСО \00FY, находящихся в поверхностном слое
4. Рувинский Г. Я. Разработка кубических единых уравнений состояний для описания термодинамических свойств веществ и систем компрессорное масло—фреон // Теплофизические свойства веществ и материалов (ГСССД. Сер. «Физические константы и свойства веществ»), 1989. №28.
5. Yokozeki A. Solubility of Refrigerants in Various Lubricants// Int. J. ofThermophys. 2001. Vol. 22. N 4.
6. Mermond Y, Feidt М., Marvillet C. Proprietes thermodynamiques et physiques des melanges de Fluides frigorigemes et d'huiles // International Journal of Refrigeration. 1999. Vol. 22. Issue 7.
7. Zhelezny V. P. The Methods of Prediction of the Properties for Substances on the Coexistence Curve Including Vicinity of the Critical Point // Proceedings of the NATO Advanced Research Workshop on Nonlinear Dielectric Phenomena in Complex Liquids. Jaszowiec-Ustron, Poland, 10—14 May 2003, edited by S. J. Rzoska, V. P. Zhelezny.
8. Zhelezny V. P., Cheremnykh Y. V., Bisko E. G. New vapor pressure equation for non-associated substances // High Temperatures — High Pressures. 1999. Vol. 31. N 2.
9. Zhelezny V. P., Semenyuk Yu. V., Ancherbak S. N., Emel'yanenko N. V. The Temperature Dependence of Parachor // Russian J. of Phys. Chem. A. 2009. Vol. 83. N 2.
10. Анисимов М. А. Критические явления в жидкостях и жидких кристаллах. — М.: Наука, 1987.
11. Zhelezny P. V., Zhelezny V. P., Procenko D. А., Ancherbak S. N. An experimental investigation and modelling of the thermodynamic properties of isobutane-compressor oil solutions. Some aspects of experimental methodology //1 nt. J. of Rerigeration. 2007. Vol. 30.
12. Ниченко С. В., Семенюк Ю. В., Железный В. П. Экспериментальные исследования энтальпии растворов хладагента R600a с маслом Reniso WF 15А // Холодильная техника и технология. 2008. № 1 (111).
13. Zhelezny V. P., Semenyuk Yu. V., Ancherbak S. N., Skripov P. V. Investigacion experimental у modelo de las propiedades termodinamicas de las soluciones R245fa / lubricante del compresor Localizacion // Frio-calor у aire acondicionado. 2008. N 399.
14. Нугачевич П. П., Бегляров Э. М., Лавыгин И. А. Поверхностные явления в полимерах. — М.: Химия, 1982.
15. Русанов А. И. Фазовые равновесия и поверхностные явления. — Л.: Химия, 1967.
