CC BY
29
УДК 621.039.534:621.574.3
Диэлектрическая проницаемость аммиака как рабочего тела холодильных установок
Д-р техн. наук К. М. АРЕФЬЕВ
Санкт-Петербургский государственный политехнический университет 195251, Санкт-Петербург, Политехническая ул., 29 Д-р техн. наук Ю. В. МУЛЕВ НПО «ЮМАС»
121552, Москва, Ярцевская ул., 29—2 Канд. техн. наук О. В. БЕЛЯЕВА, Т. А. ЗАЯЦ ГНУ «ОИЭЯИ — Сосны» НАН Беларуси 220109, Республика Беларусь, г. Минск, ул. академика А. К. Красина, 99
The actual experimental findings array on the ammonia dielectric permeability in the wide range of parameters of thermodynamic states was analyzed. An equation was deduced which generalizes the experimental findings array according to the Kirkwood-Onsager model. Universal one-parameter dependence of the ammonia dielectric permeability on the ratio of the density of a medium in this state to the absolute temperature was elaborated. The ratios proposed maybe used for diagnostics of ammonia state parameters during phase transitions in refrigerating plants.
Keywords: permeability, ammonia, parameters of thermodynamic states, diagnostics, the refrigerating plants. Ключевые слова: диэлектрическая проницаемость, аммиак, параметры термодинамического состояния, диагностика, холодильные установки.
Электрофизические свойства полярных веществ являются важными характеристиками, позволяющими анализировать строение молекул вещества. Методы электрофизической диагностики, в основу которых положена диэлектрическая проницаемость, дают также возможность оценивать термодинамическое состояние рабочего тела при фазовых переходах, в двухфазном состоянии, при определении концентраций различных веществ в смесях. Контроль термодинамических параметров аммиака в условиях работы холодильных установок актуален для обеспечения стабильной работы аммиачных насосов, оптимизации энергетических затрат на привод современных винтовых компрессорных агрегатов, повышения энергетической эффективности систем хладоснабжения.
Цель настоящей работы — определение вида зависимости диэлектрической проницаемости, как измеряемого параметра (электрическая емкость), от термодинамических свойств на основе опытных данных с целью разработки систем диагностики аммиака, как рабочего тела холодильных установок.
Обзор экспериментальных данных по диэлектрической проницаемости аммиака
Опытные данные по диэлектрической проницаемости аммиака е относительно вакуума, приводимые в работах [1—11], охватывают температурный диапазон от 198 до 483 К, диапазон давлений от 0,02 до 203 М Па. В большинстве опытов величина е определялась методом колебательного контура на частотах порядка 1 МГц. Заявленная в эксперименталь-
ных работах погрешность величин е в газовой фазе при температуре Тс ниже критической составляет около 0,1%, в жидкой фазе — в пределах 0,1—0,5%, включая жидкость при давлении насыщения, и достигает 1,5% в однофазной области при Т> Т [5]. В наиболее ранних работах [1—2], а также в справочных изданиях [8—10] погрешности опытных измерений не указаны. В табл. 1 даны ссылки на литературные источники, диапазоны параметров состояния и число опытных точек по диэлектрической проницаемости в различных областях термодинамического состояния, которые использованы при построении обобщающих уравнений.
Предварительный анализ массива опытных данных проведен с использованием существующих уравнений для расчета диэлектрических характеристик. Наши оценки показывают, что разброс опытных данных в газовой фазе при Т< Т составляет около 0,1%. Опытные данные по диэлектрической проницаемости жидкого аммиака согласуются со значениями, полученными в условиях, отвечающих пограничной кривой со стороны жидкости. Разброс опытных точек в жидкой фазе, а также в однофазной области при сверхкритических параметрах достигает +5% и определяется неточностью измерений. Из всего массива рассмотренных опытных данных диэлектрической проницаемости аммиака около половины (98 точек) приведены в работе [5]. Следует отметить, что в литературе не встречаются опытные данные, полученные в условиях, отвечающих пограничной кривой со стороны пара, а экспериментальное исследование насыщенной жидкости ограничено температурой 323 К. Области
Таблица 1
Интервалы параметров состояния и число опытных точек по диэлектрической проницаемости аммиака в газе, жидкости и однофазной области при сверхкритических параметрах
Термодинамическое состояние Число точек Интервал
температура, К давление, МПа
Газовая фаза при Т<Т [1,2,9-11] 48 227-398 0,045-7
Жидкая фаза [5] 32 213-403 0,02-203
Жидкость при давлении насыщения [3—8] 41 198-323 -
Однофазная область при Т> Т [1,2,5] ‘ 91 413-483 2-119
200
300
400
Г, К
нения является модель Кирквуда—Онзагера, доработанная Harris и Alder, в которой использован так называемый g-фактор Кирквуда. Этот фактор определяет взаимодействие молекул вещества при различных термодинамических параметрах. В интерполяционном уравнении используется соотношение для g-фактора Кирквуда в форме, предложенной Harris и Alder
(в — 1) _ Яф Г а, , gH2 3 Is»
(£ + 2)
9г
ЗkTh, (2є+1)(є + 2)
(1)
где к — константа Больцмана; к = 1,380658 Дж-К-1;
уУа—число Авогадро, Л^А = 6,022 136 7-1023моль'; а — средняя молекулярная поляризуемость изолированной молекулы аммиака, а = 2,5146-1040 С2Дж_| м2;
е0 — диэлектрическая проницаемость вакуума, е0 = 0,885419-10-" СЩж-'м-1;
р — дипольный момент изолированной молекулы аммиака, р = 1,480 = 4,93675-1030 С-м; р — молярная плотность; рс — критическая плотность аммиака, рс= 0,13212- 105моль/м3;
Т— критическая температура аммиака, Т=405,4 К; Т= Т/Т, р= р/р — приведенные параметры состояния.
Значения дипольного момента молекулы и критических параметров аммиака взяты из работы [13].
Уравнение (1) перепишем как
е - 1 е + 2
= А-
Рис. 1. Экспериментальные данные по статической диэлектрической проницаемости аммиака: о — газовая фаза при Т < Т; ш — жидкая фаза, включающая данные при давлении насыщения;
А — однофазная область при Т> Т
параметров состояния, в которых проведены экспериментальные исследования диэлектрической проницаемости аммиака, показаны на рис. 1.
Обобщение массива опытных данных по диэлектрической проницаемости аммиака уравнением Кирквуда—Онзагера
Рассмотрим согласованность набора экспериментальных данных по диэлектрическим константам аммиака (относительной диэлектрической проницаемости) е при различных параметрах состояния. Для решения этой задачи как прототип использована модель для полярного вещества с водородной связью, которая была принята в качестве стандарта при обобщении статической диэлектрической проницаемости воды и пара Международной Ассоциацией по свойствам воды и водяного пара — МАСВВП (The International Association for the Properties of Water and Steam — IAPWS) [12]. Физической основой формулировки этого урав-
(2е+ 1)(е + 2) где А и В имеют вид
А^д2 pg Л & к ' Т’
п- МАа В ~ Зє„ Р-
+ 5,
(2)
(3)
(4)
Физически корректный корень уравнения (2) для диэлектрической константы равен
\ + А + 5В + уі9+2А + \8В+А2 + ЮАВ + 9В2 ,г,
е=---------------------------------------Р)
4-4 В
Для значений фактора g в нашей задаче по аналогии с работой [12] выбрано следующее соотношение:
g = 1 + А рР Г"-25 + Аг рг Тр + А, р2 Т'р5 +
+ А4 р2 Гр 25 + As р3Р Г/ + А, р3„ г2'5
(6)
При обобщении опытных значений е аммиака при различных параметрах состояния использован массив данных по диэлектрической проницаемости в интервалах температур 198—483 К, давлений 0,02—203 МПа. Значения постоянных коэффициентов А. получены методом нелинейной регрессии (табл. 2), исходя из условия, что сумма квадратов отклонений опытных значений е от величин, рассчитанных по уравнению (1) минимальна.
На рис. 2 представлены значения корреля-
Таблица 2
Коэффициенты А. уравнения (6) для ^-фактора
і А.
1 -0,587376
2 1,44497
3 -0,423389
4 0,0299065
5 0,0501455
6 0,00465624
ционного параметра Кирквуда g, рассчитанные по уравнению (1) для опытных данных по диэлектрической проницаемости в трех областях термодинамического состояния от параметра р/Т.
Изменение корреляционного фактора g, рассчитанного по уравнению (1), находится в пределах от 1 в области газовой фазы и сверхкритичес-ких параметров до 1,75 в области жидкой фазы. Из рис. 2 видно, что фактор g является однозначной функцией параметра р/ Т. только в интервале (р/Т)<(рс/Т) = 0,555.
При р/Т> 0,555 на рис. 2 можно видеть изолинии, которые отвечают различным изохорам. Полученные результаты позволяют делать вывод, что во всех трех областях термодинамического состояния появляются ассоциаты, наблюдаемый дипольный момент которых ц* = #|/2ц больше дипольного момента изолированной молекулы аммиака ц, т. е. фактор g больше единицы.
Ддя рассмотренного массива опытных данных среднеквадратическое отклонение (СКО) экспериментальных значений е от величин, рассчитанных посоотношению( 1 )составило2,5%снадежностью0,95. В каждой отдельной области параметров состояния СКО составляет: в газовой фазе — 0,94%, в жидкой фазе — 3,2%, в однофазной области при сверхкри-тических параметрах — 2,4% с надежностью 0,95. Таким образом, результаты обобщения опытных данных для трех областей термодинамического состояния позволяют сделать вывод, что разброс опытных точек в жидкой фазе и в однофазной области при сверхкритических параметрах количественно согласуется с результатом анализа, выполненного отдельно для этих двух областей параметров состояния.
Для расчета диэлектрической проницаемости жидкости и газа при давлении насыщения для величин давления насыщения аммиака р на основе
‘ s
данных [13] нами получено соотношение в форме уравнения Вагнера:
Ps = Рсexp [(TJT) (—7,128002т + (Г)
+ 0,962513х15-2,992915т3'5)],
где рс = 11,333 МПа — критическое давление [13]; т= 1-Т/Т.
Значения диэлектрической проницаемости жидкости и пара при давлении насыщения, г' и е", рассчитаны по уравнению (1), диэлектрическая проницаемость аммиака при параметрах, отвечающих критической точке по модели Кирквуда—Онза-
Рис. 2. Зависимость значений фактора g, рассчитанных по уравнению (1) для опытных данных по є в трех областях термодинамического состояния, от параметра р/Т: о — газовая фаза при Т < Т; ш — жидкая фаза, включающая данные при давлении насыщения;
Д — однофазная область при Т > Т
200 300 400 500 Т, К
Рис. 3. Диэлектрическая проницаемость аммиака по модели Кирквуда—Онзагера (1) в интервале температур от 200—485 К и давлений 0,05—120 МПа:
1 и 2 — диэлектрическая проницаемость аммиака пара и жидкости при давлении насыщения, г" и г’, соответственно; 3—9— изобары, отвечающие давлениям 0,1; 5; 10; 20; 40; 80 и 120 МПа
гера (1) равна 4,1.
На рис. 3 представлены расчетные значения диэлектрической проницаемости аммиака в газовой и жидкой фазах, а также в однофазной области при сверхкритических параметрах, рассчитанные по модели (1) в интервале температур от 200—485 К и давлений 0,05-120 МПа.
Следует отметить, что часть расчетных значений диэлектрической проницаемости аммиака выходит за пределы параметров измерений.
Однопараметрическая зависимость диэлектрической проницаемости аммиака от отношения плотности к температуре
Уравнение, рассмотренное выше сложно по форме и неудобно для практических расчетов. В [ 14] Ю. В. Мулевым и С. Н. Смирновым основе экспериментальной базы данных предположена методика расчета диэлектрической проницаемое-
ти на основе зависимости е от одного параметра — отношения плотности среды в данном состоянии к абсолютной температуре р/Т как для однофазных состояний, так и для жидкости и пара при давлениях насыщения. Однако не проработанность в [14] всего массива экспериментальных данных, недостаточность математической обработки затруднило применение предложенных зависимостей в разработке диагностических систем.
В настоящей работе выполнено обобщение всего существующего массива имеющихся опытных данных по е аммиака в трех областях термодинамического состояния: газообразном при Т < Г, жидком, включая данные при давлении насыщения и однофазном состоянии при Т > Г (рис. 4), проведена их коррекция с соответствующими температурными шкалами и эти данные проанализированы на основе существующих статистических теорий погрешностей. В настоящей работе предлагается единое уравнение расчета диэлектрической проницаемости аммиака при различных экспериментально исследованных параметрах состояния.
В координатах е-р/Т все опытные значения е укладываются на единую кривую. Имеются диапазоны параметра р/Т, в которых одновременно располагаются опытные данные, полученные в разных областях термодинамического состояния. В интервале р/Т < 0,555 расположены точки, отвечающие газообразному состоянию при Т < Т и однофазной области при Т> Т. В интервале р/Т> 0,555 находятся точки, отвечающие жидкости, включая данные при давлении насыщения, а также однофазной области при Т > Т. Данные, представленные на рис. 4, обобщены зависимостью вида
8=1+1 £ (р/77, (8)
1
где коэффициенты Кх = 2,3915; К2 = 6,8567; К} = —3,5171; К4 = 0,8685; К5 = —0,0863 найдены методом линейной регрессии. Величина СКО экспе-
Рис. 4. Зависимость экспериментальных значений е аммиака от параметра р/Т в трех областях термодинамического состояния: о — газовая фаза при Т < Тс;
• — жидкая фаза, включающая данные при давлении насыщения; Д — однофазная область при Т > Т
риментальных значений е от сглаженных для всего массива точек составила 2,9% с надежностью 0,95. Для различных областей термодинамического состояния эти отклонения составляют: в газовой фазе (0,0015 < р/Т < 0,133) при Т< Т — 0,60%, в жидкой фазе (1,12 <р/Т< 3,68) — 3,0% и в однофазном области (0,02 < р/Т< 1,30) при Т> Т — 3,5%, с надежностью 0,95. При параметрах, отвечающих критической точке, величина диэлектрической проницаемости, рассчитанная по уравнению (8) равна 3,9 и в пределах погрешности расчета согласуется со значением, рассчитанным по уравнению Кирквуда—Онзагера (1).
Наибольший разброс имеют опытные данные, отвечающие однофазной области при температуре выше критической в окрестности критической точки в интервале 0,1 > р/Т < 1, что может быть связано со сложностью проведения эксперимента в указанной области параметров.
Таким образом, выполненное нами обобщение существующего массива опытных данных по диэлектрической проницаемости аммиака в широком диапазоне параметров состояния показало возможность применения упрощенной модели представления е как однопараметрической функции отношения плотности к температуре. Предложенная методика расчета диэлектрической проницаемости отличается простотой и наглядностью. На зависимости е-р/Тдля всего массива экспериментальных данных расслоение по фазовой составляющей не наблюдается. Это позволяет утверждать, что для описания опытных данных по е аммиака, как в однофазных областях, так и на линиях фазового перехода можно использовать однопараметрическую зависимость, т. е. говорить об универсальности этой зависимости.
Однопараметрическая зависимость диэлектрической проницаемости теплоносителя может быть использована для диагностики параметров аммиака при фазовых переходах в холодильных установках. Используя предложенную зависимость (8) и измеряемую полезную (в вакууме) электрическую емкость первичного преобразователя при учете температуры (для двухфазных состояний давление и температура взаимосвязаны), по значению электрической емкости можно определить плотность. Вследствие значительной разности диэлектрической проницаемости аммиака жидкости и пара при давлении насыщения (до 25 раз) измерения диэлектрической константы двухфазной смеси будут однозначно характеризовать соотношение жидкой и паровой составляющих в этом веществе. Возможность распространения однопараметрической зависимости диэлектрической проницаемости на двухфазную область состояния анализировалась для водного теплоносителя в [ 15].
Заключение
Выполненная аппроксимация существующего массива опытных данных по диэлектрической проницаемости аммиака в трех областях термодинамического состояния по модели Кирквуда-Он-
загера, доработанной Harris и Alder с применением g-фактора Кирквуда, позволяет описывать е аммиака со среднеквадратическим отклонением 2,5% при надежности 0,95.
Предложенная универсальная однопараметрическая зависимость диэлектрической проницаемости от отношения плотности к температуре в однофазных областях и при параметрах фазового перехода более удобна в применении и позволяет проводить расчеты е аммиака со среднеквадратическим отклонением 2,9% при надежности 0,95. Кроме того, предложенная зависимость может быть использована для разработки систем контроля термодинамических параметров аммиака, как рабочего тела холодильных установок.
Для уточнения полученных обобщений необходимо более детальное экспериментальное определение диэлектрической проницаемости аммиака при различных параметрах термодинамического состояния.
Список литературы
1. Keyes F. G., Kirkwood J. G. The dielectric constant of ammonia as a function of temperature and density// Phys. Rev. 1930. Vol. 36, № 15.
2. Uhlig H. H., Kirkwood J. G., Keyes F. G. The dependence of the dielectric constant of gases on temperature and density // J. Chem. Phys. 1932, Vol. 54. №4.
3. Hooper G. S., Kraus C. A. A. Properties of electrolytic solutions. XIV. The dielectric constant of some solutions of electrolytes in benzene // J. Amer. Chem. Soc. 1934. Vol. 56. №1.
4. Grubb H. М., Chittum J. F., Hunt H. Liquid ammonia as a solvent. VI. The dielectric constant of liquid ammonia//J. Amer. Chem. Soc. 1936. Vol. 58.
5. ВиЬаск М., Harder W. D. The static dielectric constant ofammonia to high pressures and temperatures.
1. Computer assisted data evaluation // Ber. Bunsenges. Phys. Chem. 1977. №81.
6. Billaud G., Demortier A. A. Dielectric constant of Liquid Ammonia from —35 to +50°C and its influence on the Association between Solvated Electrons and Cation. //J. Phys. Chem. 1975. Vol. 79. №26.
7. Brady М. M. Loss Measurements ofWet Textiles at 9 GHz//Journal of Microwave. 1968. №3 (4).
8. Ахадов Ю. Я. Диэлектрические параметры чистых жидкостей. — М.: МАИ, 1999.
9. Никольский Б. П., Григоров О. Н., По-зин М. Е. Общие сведения строения вещества. Свойства важнейших веществ. Лабораторная техника. Справочник химика — М. —Л.: Химия. 2-е изд., перераб. и доп., 1966. Т. 1.
10. Чайлдс У. Физические постоянные. Справочное пособие. — М.: Физматгиз. 2-е изд. ис-правл., 1962.
11. Watson Н. Е. The Dielectric Constants of Ammonia, Phosphine and Arsine // Proceedings of the Royal Society of London. 1927. Series A. Vol. 117. №776.
12. Fernandez D. P., Goodwin A. R. H., Lemmon E. Wet al. A Formulation for the Static Permittivity of Water and Steam at Temperatures from 238 К to 873 К at Pressures up to 1200 MPa, Including Derivatives and Debye-Huckel Coefficients. //J. Phys. Chem. Ref. Data. 1997. Vol. 26. №4.
13. Lemmon E. W, McLinden М. O., Hubner M. L. NIST Standard Reference Database 23, NIST thermodynamic and transport properties of refrigerants and refrigerant mixtures — REFPROP, version 7.1, (2003) // Standard Reference Data Program/National Institute of Standards and Technology, 2003.
14. Мулев Ю. В., Смирнов С. Н. Методика расчета диэлектрической проницаемости воды до 923 К и 1 ГПА//ТВТ. 1992. №1.
15. Арефьев К. М., Беляева О. В., Мулев Ю. В., Мулев М. Ю., Заяц Т. А. Обобщение данных по диэлектрической проницаемости жидкой и паровой фаз воды с целью определения паросодержа-ния двухфазных состояний // ИФЖ. 2011. Т. 84. №5.
