Теплофизика высоких температур, 2021, T. 59, № 5, стр. 679-685

Плотность жидкого калия на линии насыщения: краткая история длиною в 50 лет

Л. Р. Фокин 1*, Е. Ю. Кулямина 1**

1 Объединенный институт высоких температур РАН
Москва, Россия

* E-mail: lfokin@mail.ru
** E-mail: kulyamina.elena@gmail.com

Поступила в редакцию 17.11.2020
После доработки 17.11.2020
Принята к публикации 19.05.2021

Полный текст (PDF)

Аннотация

В статье рассматривается созданная в Объединенном институте высоких температур РАН установка уникальной конструкции, реализующая методику пикнометра постоянного объема и переменной массы. Получены данные о плотности жидкого калия на линии насыщения и на изобаре 100 атм в интервале температур 336−2030 К. Отмечается связь результатов реализованной в ОИВТ РАН методики и измерений сжимаемости воды и водяного пара, проводимых в МЭИ, начиная с работы В.А. Кириллина 1953 г. Анализируется достоверность экспериментальных данных о плотности жидкого калия 1976–1977 гг. в сравнении с результатами новых независимых измерений 1990 г. Найдено положение идеального, как у разреженного газа, коэффициента термического расширения жидкости на линии насыщения и на изобаре калия. Экспериментальные данные о плотности, сжимаемости, идеальных коэффициентах термического расширения использованы для оценки положения линии единичного фактора сжимаемости и линии инверсии эффекта Джоуля–Томсона для плотного калия в интервале температур $\sim {\kern 1pt} 1600{\kern 1pt} - {\kern 1pt} 2300$ К и давлений до $\sim {\kern 1pt} 2$ кбар. Рассмотренные данные на термодинамической поверхности жидкого калия нуждаются в проведении новых экспериментальных и теоретических работ и расширении перечня анализируемых свойств.

ВВЕДЕНИЕ

В связи с юбилеем ОИВТ РАН есть все основания вспомнить сотрудников и работы, которые проводились на первых шагах становления института, проследить примеры использования результатов тех лет в последующие годы.

Вспоминается отдел жидкостей и газов начала 1970-х гг., помещения которого располагались в основном на пятом этаже института на Красноказарменной ул., 17а. Это было время, когда руководитель отдела Шпильрайн Эвальд Эмильевич уже защитил докторскую диссертацию (1963 г.) и была опубликована монография “Теплофизические свойства щелочных металлов”, подготовленная коллективом авторов МЭИ и ИВТАН [1]. Основные направления работ отдела были связаны с исследованиями теплофизических свойств рабочих тел и теплоносителей атомной энергетики и МГД-процессов на жидких металлах. Можно вспомнить, что, когда осенью 1956 г. один из авторов пришел на кафедру инженерной теплофизики МЭИ для выполнения дипломной работы, Шпильрайн Э.Э. взял амбарную книгу и на обложке крупно написал ЩЕЛОЧНЫЕ МЕТАЛЛЫ.

С левой стороны зала Л-501 на пятом этаже институтского корпуса располагался ряд стендов для исследования свойств щелочных металлов, на которых работали молодые сотрудники, в будущем известные специалисты по теплофизическим свойствам веществ:

− Мозговой А.Г.: капающий пикнометр с жидким металлом, расположенный в автоклаве с чувствительными весами. Управление весами под давлением инертного газа осуществлялось через стенку камеры с помощью сельсинов, ${{T}_{{\max }}} = 2030$ К, ${{p}_{{\max }}} = 250$ атм;

− Каган Д.Н.: метод смешения, декремент энтальпии жидких металлов, блок: высокотемпературная печь – кипящий водяной калориметр, ${{T}_{{\max }}} = 1350$ К, герметичная ампула с магнитным подвесом, позволяющим проводить эксперимент при новых параметрах без разборки установки;

− Савченко В.А.: электрическое сопротивление столба чистого металла и сплавов, ${{T}_{{\max }}} = 1300$ К;

− Крайнова И.Ф.: теплопроводность, метод осевого теплового потока для столба жидкости, ${{T}_{{\max }}} = 1300$ К;

− Фомин В.А. и Сковородько С.Н.: вязкость жидких металлов, метод затухающих крутильных колебаний стаканчика, целиком заполненного жидким металлом без свободных поверхностей, ${{T}_{{\max }}} = 1300$ К.

В правом крыле зала Л-501 находились, в частности, два стенда для измерения давления насыщения пара, здесь работали:

− Белова А.М.: метод сброса давления, ${{T}_{{\max }}} = 1200$ К;

− Никанорова Э.В.: метод точек кипения, тарировка установки по давлению насыщения ртути, ${{T}_{{\max }}} = 1300$ К.

Результаты проведенных исследований свойств щелочных металлов представлены публикациями Шпильрайна Э.Э. и соавторов в основном в журнале “Теплофизика высоких температур”, например [2], и в монографии [1]. По этой тематике известна единственная статья отдельно Эвальда Эмильевича [3].

Каждая из перечисленных установок, методика эксперимента, результаты содержат много интересного. Но здесь остановимся на работах Мозгового А.Г., посвященных измерениям плотности жидкого калия [46], так как, во-первых, описание установки и основные результаты не публиковались после издания монографии 1970 г. [1] и, во-вторых, в работах 2012 г. [7, 8] по идеальному коэффициенту термического расширения (КТР) на линии насыщенной жидкости были использованы экспериментальные данные [46].

УСТАНОВКА И ДАННЫЕ А.Г. МОЗГОВОГО ДЛЯ КАЛИЯ

Установка представляет собой оригинальный вариант взвешенного пьезометра переменной массы. Ранее в 1963 г. Шпильрайн Э.Э. и Якимович К.А. на установке такого типа провели измерения плотности жидких лития и рубидия до $\sim {\kern 1pt} 1150$ К, но при атмосферном давлении [9]. Мозговой А.Г., сохраняя высокую точность методики, создал аппаратуру и провел измерения плотности жидкого калия в интервале температур 347−2030 К на линии насыщения и на изобаре 100 атм [46]. Аналогов данного эксперимента неизвестно.

Ампула пикнометра, в которой находится жидкий металл при высокой температуре, подвешена на нити к аналитическим весам, расположенным в зоне комнатных температур. Пикнометр и весы находятся в кожухе под давлением инертного газа. При нагреве основного объема происходит расширение жидкого металла, и избыток жидкости удаляется через капилляр, приваренный снизу к ампуле. Равновесие на коромыслах весов, нагрузка гирь от грамм до миллиграмм устанавливается снаружи экспериментатором с помощью двух сельсинов: сельсина датчика (снаружи) и сельсина приемника (внутри). Здесь вспоминается бессальниковый насос в циркуляционном контуре установки А.Е. Шейндлина по измерению теплоемкости водяного пара.

Для пикнометра переменной массы расчетное выражение плотности исследуемой жидкости при параметрах эксперимента имеет вид

(1)
$\rho = {{\left( {m - \Delta m - {{m}_{{{\text{кап}}}}}} \right)} \mathord{\left/ {\vphantom {{\left( {m - \Delta m - {{m}_{{{\text{кап}}}}}} \right)} {\left( {{{V}_{0}}\left( {1 + \alpha \Delta T} \right)} \right)}}} \right. \kern-0em} {\left( {{{V}_{0}}\left( {1 + \alpha \Delta T} \right)} \right)}},$
где m – начальная масса жидкости в ампуле; Δm, mкап – массы удаленной жидкости и жидкости в капилляре; ${{V}_{0}}\left( {1 + \alpha \Delta T} \right)$ – начальный объем и поправка на термическое расширение ампулы. Весы, естественно, фиксируют общий вес жидкости и ампулы пикнометра. В работах [46] много внимания уделено деталям эксперимента, определяющим высокий уровень результатов, таким как чистота исходного металла до и после проведения экспериментов, очистке инертных газов, тарировке весов при работе под давлением, специальным исследованиям коэффициентов термического расширения конкретных образцов стали и сплавов молибдена, использованных при изготовлении ампул пьезометров, и т.п. Чтобы получить значение плотности на бинодали жидкой фазы, в установке при заданной температуре создавалось давление на несколько атмосфер выше расчетного давления насыщения, повышалась температура, через капилляр сбрасывался избыток жидкости и фиксировались показания весов. В каждой серии экспериментов постепенно повышалась температура и, если нужно, давление, что позволяло получать массив данных на изобаре или на бинодали жидкости без разборки установки.

Результаты экспериментов для калия в интервале температур 347−1785 К опубликованы в [4] и затем позже в интервале температур 1642–2027 К в [5, 6]. В статье [4] приведены экспериментальные данные ${{\rho }}\left( T \right)$ для пяти серий экспериментов на линии насыщения и двух серий на изобаре 100 атм. Это на примере калия дает возможность, в частности, проследить, как лежит изобара плотности по отношению к линии плотности насыщенной жидкости. На основании анализа случайных и систематических погрешностей эксперимента здесь приводятся оценки доверительных погрешностей экспериментальных данных: 0.15% ниже 1000 К, 0.2% до 1500 К и 0.25% при более высокой температуре. Краткие публикации [5, 6] включают 8 экспериментальных точек ${{\rho }}\left( T \right)$, полученных с применением ампул из разных сплавов. Для измерения температуры в интервале 1300−2100 К применялись термопары PtRh 6-30. В статьях [5, 6] не сообщается, проверялась ли тарировка термопар после высокотемпературного эксперимента. Авторская оценка погрешности этих данных 0.3%. Все рассмотренные выше данные на линии насыщения калия аппроксимированы единым полиномом шестой степени по температурам.

Работы подобного рода проводились в 1970-е гг. в МЭИ и в ИВТАН и по другим свойствам газов, жидкостей, твердых тел, и именно они обеспечили признание и авторитет теплофизической школы института. Надо отметить преемственность, связь этих работ с исследованиями свойств веществ на кафедрах МЭИ. В частности, пьезометры переменной массы широко использовались на кафедрах ТОТ (теоретических основ теплотехники) и ИТФ (инженерной теплофизики) для определения термодинамических свойств воды и водяного пара, диоксида углерода и других рабочих веществ при высоких параметрах. Если для определения плотности жидкого калия измерялся переменный вес капающего пикнометра (см. уравнение (1)), то на установках для определения удельного объема воды и водяного пара проводился на изотерме выпуск вещества из пьезометра в кюветы, масса выпусков по ходу эксперимента суммировалась и тем самым определялся удельный объем от низких до высоких давлений. Описание конструкции и работы пьезометра переменной массы в работе Кириллина В.А. и Румянцева Л.И. приводится в Трудах МЭИ за 1953 г. [10], а также в учебнике [11]. Этот экспериментальный метод и полученные результаты послужили основой докторской диссертации Кириллина В.А. Область параметров при исследовании свойств воды и водяного пара по этой методике постепенно расширялась, одновременно улучшалась приборная база и повышалась точность измерений удельного объема. Зубарев В.Н. провел новые измерения и опубликовал в 1955 г. экспериментальные данные до 500°С и 950 атм [12]. Улыбин С.А. расширил массив экспериментальных данных в околокритической области [13]. Александров А.А., используя методику пьезометра переменной массы, подготовил и защитил диссертацию на основе результатов измерений удельных объемов воды и водяного пара до 1200 атм и 900°С [14]. Через 10 лет эти исследования были расширены до 2000 атм [15]. Целое поколение экспериментаторов МЭИ при исследовании сжимаемости газов и газовых смесей использовало метод присоединяемых баллончиков с адсорбентом для отбора вещества из пьезометра, который применял Алтунин В.В. в 1959 г. при работе с диоксидом углерода [16]. На примере рассмотренных работ для щелочных металлов [16], воды и водяного пара [1015], диоксида углерода [16] можно видеть перспективы применения пикнометров и пьезометров переменной массы при исследовании плотности жидкостей и газов в зависимости от температуры и давления.

Возвращаясь к работам по исследованию плотности жидких щелочных металлов, надо отметить, что в значительной степени по инициативе Мозгового А.Г. был подготовлен и опубликован в 1983 г. обзор экспериментальных данных. Также были предложены рекомендуемые значения плотности и КТР с оценками погрешностей на линии насыщения всех щелочных металлов, кроме франция, при температурах, на 10−15 К не доходящих до принятых в те годы оценок критических параметров [17], а именно для Li до 3660 К, для Na до 2500 К, для K до 2270 К, для Rb до 2080 К и для Cs до 2030 К [18]. В этой работе процедура анализа данных включала два этапа. Вначале было проведено обобщение экспериментальных данных о плотности всех металлов на линии насыщения в приведенных координатах ${{{\rho }} \mathord{\left/ {\vphantom {{{\rho }} {{{{{\rho }}}_{c}}\left( {{T \mathord{\left/ {\vphantom {T {{{T}_{c}}}}} \right. \kern-0em} {{{T}_{c}}}}} \right)}}} \right. \kern-0em} {{{{{\rho }}}_{c}}\left( {{T \mathord{\left/ {\vphantom {T {{{T}_{c}}}}} \right. \kern-0em} {{{T}_{c}}}}} \right)}}$ в рамках метода наименьших квадратов с весами с помощью полинома шестой степени. Затем для каждого металла отдельно результаты этого обобщения и индивидуальные данные о плотности обрабатывались совместно в координатах $\rho \left( T \right)$. Результаты повторной обработки для каждого элемента были приняты в качестве рекомендуемых значений о плотности и КТР. Обзор [17] включает оценки погрешностей этих свойств, в частности для плотности калия на линии насыщения 0.25% до 1300 К, 0.5% в интервале 1300−1500 К, 1% при Т > 1500 К. Позже этот обзор был воспроизведен в издании IUPAC по свойствам щелочных металлов [19], а рекомендуемые таблицы плотности этих элементов на линии насыщения включены в справочник по жидкометаллическим теплоносителям 1988 г., подготовленный Каганом Д.Н. и коллегами [20].

Что касается рекомендуемых значений плотности жидкого калия, то они до 1600 К практически воспроизводят экспериментальные данные Мозгового А.Г. 1976 г. [4]. Однако при более высоких температурах рекомендуемые значения лежат систематически выше, а экспериментальные данные ниже. В интервале температур 1730−2026 К отклонение экспериментальных данных 1977 г. [5, 6] увеличивается с 3 до 15%. В 1985 г. Мозговой А.Г. защищал кандидатскую диссертацию в форме научного доклада (руководитель Якимович К.А.) [21] через 8 лет после публикации своих экспериментальных данных по калию и через два года после издания обзора [18]. В качестве результата исследований автореферат включает таблицу данных о плотности, КТР и коэффициенте сжимаемости на линии насыщения калия до 2280 К. Однако выше 1600 К для плотности автор приводит не свои экспериментальные данные, а рекомендуемые значения из обзора [18]. Какие-либо соображения по этому поводу не приводятся. Подробное описание установки ИВТАН с капающим пикнометром, но без какого-либо анализа полученных экспериментальных данных о плотности жидкого калия приводится в международном справочном издании 1992 г. [22].

Ампула, капилляр и сам жидкий металл в рассматриваемом эксперименте находятся под давлением инертного газа, а именно гелия. Известно, что Мозгового А.Г. постоянно беспокоил вопрос, не мог ли газ из объема пикнометра в заметных количествах диффундировать через холодный капилляр в горячий объем ампулы и в жидкий калий, плотность которого на линии насыщения при высоких параметрах резко падает. Эта проблема растворимости газов в металлах нашла отражение в серии публикаций Шпильрайна Э.Э., Сковородько С.Н., Мозгового А.Г. [2325]. Была разработана “дырочная” модель растворимости инертных газов в металлических расплавах, и, например, в работе [24] предложены оценки растворимости гелия во всех щелочных металлах при температурах 600−1500 К и давлениях до 10 МПа. Однако попытки перенести эти оценки на анализ диффузии и растворимости гелия в жидком калии для реальной установки Мозгового А.Г. отсутствуют. Вопрос о том, почему высокотемпературные данные [5, 6] оказались систематически ниже результатов обобщений [18, 19], остается открытым.

В 2000 г. Хензел с соавторами опубликовали новые результаты измерений плотности пара и жидкости на линии насыщения калия от 600 К до температуры, близкой к критической [26], которые совместно с данными [5, 6] образуют представительную базу данных о свойствах жидкого калия. В работе [26] была использована та же экспериментальная методика фиксации изменения наклона изохор при прохождении линии насыщения рубидия и цезия [27]. Давление на калий в пьезометре постоянного объема через капилляр волосяного сечения передавалось от внешней среды (аргон), что позволяло наблюдать термограмму процесса нагрева–охлаждения. Измерения температуры проводились с помощью тарированных W–Re-термопар. По оценкам авторов, погрешность определения температуры фазового перехода составляла 2−3 градуса. В работе [26] на рисунке в координатах ${{\rho --}}T$ в интервале 600−2160 К показаны экспериментальные точки на обеих ветвях бинодали, а также их диаметр ${{\left( {{{\rho }}{\kern 1pt} '\,\, + {{\rho }}{\kern 1pt} ''} \right)} \mathord{\left/ {\vphantom {{\left( {{{\rho }}{\kern 1pt} '\,\, + {{\rho }}{\kern 1pt} ''} \right)} 2}} \right. \kern-0em} 2}$. Видно, что эти линии при экстраполяции сходятся к температуре $T\sim 2178$ К и к значению плотности ${{\rho }}\sim 0.17$ г/см3, которые объявлены критическими параметрами калия. Данные оценки отличаются от рекомендованных ранее [18, 19] величин, на которые ориентировался Мозговой А.Г. в своей диссертации [21]. Независимые результаты измерений критических параметров цезия и рубидия [27, 28] также дают новые по сравнению с [18, 19] оценки критических параметров. Есть основания именно новые экспериментальные данные для критических параметров щелочных металлов использовать в дальнейшем при разработке различных корреляций термодинамических свойств этих металлов. Отметим, что исследования свойств щелочных металлов на линии насыщения проводились также в ЭНИН им. Г.М. Кржижановского в 1970 г. В частности, Е.И. Чернеев проводил измерения давлений насыщения и критических параметров калия [29]. В эксперименте был использован метод сброса давления и фиксации температурного отклика. Известно, что вблизи критической точки этот метод теряет чувствительность, и к этим результатам относились весьма осторожно. Тем не менее данные Чернеева Е.И. для калия ${{T}_{c}} = 2140 \pm 13$ К, ${{p}_{c}} = 14.9 \pm 0.8$ МПа весьма близки к результатам 2000 г. [26]. В то же время в работах [26, 29] не приводятся важные детали температурных измерений в районе 2000 К, такие как температурное поле и контроль тарировки термопар после цикла измерений. Поэтому весьма вероятно, что для критической температуры, предложенной в [26], надо принимать осторожную оценку погрешности $\delta {{T}_{c}} = 10$ К, пока не будут получены новые экспериментальные данные.

На рисунке в работе [26] для жидкого калия на линии насыщения показаны 23 экспериментальные точки в интервале 600–2175 К. Современная компьютерная техника позволяет с высокой точностью провести оцифровку таких данных, и на рис. 1 дано сравнение старых и новых данных на бинодали жидкой фазы.

Рис. 1.

Значения плотности жидкого калия на линии насыщения: 1 – данные [18, 21], 2 – высокотемпературные экспериментальные данные [5, 6], 3 – экспериментальные данные [26].

Можно видеть, что приблизительно до 1600−1650 К экспериментальные данные ИВТАН 1976 г. и новые данные 2000 г. [26] в границах оценок их погрешностей ~0.2 и 1.25% практически совпадают. Этот факт является важным результатом проведенного анализа свойств жидкого калия. От 1650 К новые экспериментальные данные в интервале температур 1780–2030 К лежат заметно выше данных [5, 6].

В свою очередь, результаты обобщения экспериментальных данных о плотности щелочных металлов на линии насыщения [18], в частности для калия, по отношению к данным работы [26] лежат выше и по температурам до 2270 К, т.е. попадают в область температур, которые выше новой оценки критической температуры калия.

Для обобщения экспериментальных данных в рамках весового МНК-метода был создан массив ${{\rho }}\left( T \right)$ значений ${{n}_{i}}$ с оценками их погрешностей: в интервале 336−1600 К ${{n}_{1}} = 14$, $\delta {{\rho }} = 0.3\% $ [21]; при 1700 К ${{n}_{2}} = 1$, $\delta {{\rho }} = 0.5\% $ [21]; 600–2130 К ${{n}_{3}} = 20$, $\delta {{\rho }} = 1.5\% $ [26]. Для аппроксимации этих данных в интервале от температуры плавления до критической было рассмотрено несколько параметрических зависимостей и в результате с учетом небольшого числа данных в околокритической области было принято простое выражение, берущее начало от функции Гуггенхейма [30]:

(2)
$\begin{gathered} {{\rho }}\left( T \right) = \\ = {{{{\rho }}}_{c}}\left[ {1 + {{a}_{1}}\left( {1 - {{\tau }}} \right) + {{a}_{2}}{{{\left( {1 - {{\tau }}} \right)}}^{2}}... + {{a}_{6}}{{{{\tau }}}^{{2.5}}}{{{\left( {1 - {{\tau }}} \right)}}^{{{1 \mathord{\left/ {\vphantom {1 3}} \right. \kern-0em} 3}}}}} \right], \\ \end{gathered} $
где ${{{{\rho }}}_{c}} = 0.175$ г/см3, ${{\tau }} = {T \mathord{\left/ {\vphantom {T {{{T}_{c}}}}} \right. \kern-0em} {{{T}_{c}}}}$, ${{T}_{c}} = 2178$ К. Данное уравнение обеспечивает правильную асимптотику поведения плотности и КТР при $T \to {{T}_{c}}$. Показатель степени 1/3 в (2) принят в соответствии с логикой статистической обработки на конкретном массиве данных. В то же время надо иметь в виду, что в ряде работ для щелочных металлов авторы для критического индекса принимают значение 0.5. Это значение встречается, например, в работах Мартынюка М.М. [31], связанных с построением бинодали Cs, Rb, K и использованием трехпараметрического уравнения Ван-дер-Ваальса с переменным показателем n в выражении ${a \mathord{\left/ {\vphantom {a {{{V}^{n}}}}} \right. \kern-0em} {{{V}^{n}}}}$. Однако это предмет специального анализа, для проведения которого нужно рассматривать данные на обеих ветвях бинодали. Демпфирующий множитель ${{{{\tau }}}^{{2.5}}}$ в последнем члене (2) введен “от руки”, чтобы ослабить влияние околокритической области на регулярную часть уравнения.

Результаты аппроксимации: параметры a1 = = 5.68240, ${{a}_{2}} = - 1.51739$, ${{a}_{6}} = 1.86491$; отклонения точек массива данных от расчетной зависимости (2) меньше заданных погрешностей, кроме точки при 1700 К, где вклад в МНК-функционал – 1.3; в безразмерном МКН-функционале в три раза меньше точек. Отклонения экспериментальных данных от расчетной зависимости показаны на рис. 2. Разброс экспериментальных данных по плотности жидкого калия при $T > 2000$ К увеличивается, что отражает трудности проведения эксперимента в околокритической области.

Рис. 2.

Отклонение экспериментальных данных о плотности насыщенного калия ρэкс от результатов расчетной зависимости (2): 1 – [4], 2 – [26].

ИДЕАЛЬНЫЙ КОЭФФИЦИЕНТ ТЕРМИЧЕСКОГО РАСШИРЕНИЯ

На линии насыщенной жидкости для всех известных веществ с ростом температуры плотность падает, и произведение ${{\rho }}T$ проходит через максимум. При температуре экстремума ${{T}_{{\text{э}}}}$ на изолинии КТР ${{\alpha }} = - {{{{\rho }}}^{{ - 1}}}{{dp} \mathord{\left/ {\vphantom {{dp} {dT = {1 \mathord{\left/ {\vphantom {1 {{{T}_{{\text{э}}}}}}} \right. \kern-0em} {{{T}_{{\text{э}}}}}}}}} \right. \kern-0em} {dT = {1 \mathord{\left/ {\vphantom {1 {{{T}_{{\text{э}}}}}}} \right. \kern-0em} {{{T}_{{\text{э}}}}}}}}$ равен обратной температуре, как у идеального газа. Ранее это обстоятельство было рассмотрено для большой группы веществ металлов и неметаллов [7] на основе экспериментальных данных и обобщенных зависимостей ${{\rho }}\left( T \right)$ для жидкостей на линии насыщения. В этой работе из металлов рассматривался только цезий, для которого с учетом новых экспериментальных данных [27] ${{T}_{{\text{э}}}} = 1510$ К, значение критической температуры ${{T}_{c}} = 1938$ К и отношение ${{\theta }} = {{{{T}_{{\text{э}}}}} \mathord{\left/ {\vphantom {{{{T}_{{\text{э}}}}} {{{T}_{c}} = 0.78}}} \right. \kern-0em} {{{T}_{c}} = 0.78}}$. Устойчивость этого отношения для других щелочных металлов имеет важное значение при анализе термодинамического подобия их свойств [8].

Проведенный выше анализ данных о плотности жидкого калия на линии насыщения показывает, что температура экстремума ${{\rho }}T$ для функции (2) равна ${{T}_{{\text{э}}}} = 1645$ К, т.е. попадает в область согласованных данных о плотности [4, 26]. При этом для калия отношение ${{\theta }} = {{{{T}_{{\text{э}}}}} \mathord{\left/ {\vphantom {{{{T}_{{\text{э}}}}} {{{T}_{c}} = 0.75}}} \right. \kern-0em} {{{T}_{c}} = 0.75}}$ несколько ниже, чем приведенная выше оценка для цезия, однако эти оценки не противоречат друг другу в пределах погрешностей их определения. Значение отношения θ на практике связано, в частности, с тем, что оно позволяет оценивать критическую температуру веществ на основе ограниченного массива экспериментальных данных о плотности жидкости на линии насыщения. В работах [7, 8] показано, что для большой группы неорганических веществ, от инертных газов до ассоциированной воды, отношение θ меняется в узких пределах и равно 0.8 ± 0.02. Переход щелочных металлов в другую группу веществ требует дальнейшего анализа, в первую очередь для натрия.

В [4] представлены также две серии экспериментальных данных о плотности жидкого калия на изобаре 100 атм в интервале температур 876−1775 К. Эксперимент на изобаре в методе капающего пикнометра проще, и экспериментальные данные надежнее, чем на линии насыщения, так как на изобаре не проводится дополнительная операция тщательной подгонки состояния вещества к бинодали. По этим данным произведение ${{\left( {{{\rho }}T} \right)}_{p}}$ имеет максимум при температуре Tэ ≈ 1580 К, которой соответствует идеальный изобарный КТР (здесь ${{{{\alpha }}}_{{\text{p}}}} = {1 \mathord{\left/ {\vphantom {1 {{{T}_{{\text{э}}}}}}} \right. \kern-0em} {{{T}_{{\text{э}}}}}}$). Этому состоянию соответствует температура на линии инверсии эффекта Джоуля–Томсона (эффекта дросселирования при постоянной энтальпии). Известно, что эта линия инверсии при высоких параметрах  в p–T-координатах проходит через максимум линии единичного фактора сжимаемости zσ = 1, где  z = p/(ρRT), R – газовая постоянная. Для уравнения Ван-дер-Ваальса линия zσ = 1 (линия Бачинского, или Zeno line) в координатах ${{\rho }}{\kern 1pt} - {\kern 1pt} T$ ‒ это прямая ${{{\rho }} \mathord{\left/ {\vphantom {{{\rho }} {{{{{\rho }}}_{B}} + {T \mathord{\left/ {\vphantom {T {{{T}_{B}} = 1}}} \right. \kern-0em} {{{T}_{B}} = 1}}}}} \right. \kern-0em} {{{{{\rho }}}_{B}} + {T \mathord{\left/ {\vphantom {T {{{T}_{B}} = 1}}} \right. \kern-0em} {{{T}_{B}} = 1}}}}$. Здесь параметр ${{{{\rho }}}_{B}}$ – оценка плотности, которая получается при экстраполяции линии zσ = 1 на нулевую температуру, а параметр ${{T}_{B}}$ – обычно фиктивная температура, которая получается при экстраполяции zσ = 1 из области высоких плотностей на нулевую плотность.

В свое время на основе экспериментальных данных ИВТАН о плотности жидкого калия на линии насыщения и на изобаре 100 атм [4] были сделаны оценки среднего изотермического коэффициента сжимаемости β(Т) в широкой области температур [32]. Для плотной жидкости этот коэффициент меняется незначительно. Это позволяет провести пересчет данных о плотности насыщенной жидкости к давлению и плотности на линии единичного фактора сжимаемости. В статье [33] сделаны такие оценки для калия в интервале температур 400–1600 К, и при экстраполяции полученного отрезка Zeno line на нулевую температуру и нулевую плотность получены оценки параметров ${{{{\rho }}}_{В}} = 0.92$ г/см3 и ${{T}_{B}} = 4670$ К, близкие к рекомендуемым значениям. Значения плотности и температуры для максимума функции zσ = 1 модели Ван-дер-Ваальса в p–T-координатах вдвое меньше указанных выше величин, и оценка давления p = 2.3 × 103 бар, т.е. приблизительно в 15 раз выше оценки критического давления калия. Это дает возможность очень приближенно прямо нарисовать линию инверсии Джоуля–Томсона для калия от 100 атм и $\sim {\kern 1pt} 1580$ К к высоким параметрам ${{{{\rho }}}_{{\text{э}}}} = 0.92$ г/см3 и ${{Т}_{{\text{э}}}} = 2300$ К, естественно, в предположении монотонного изменения функции. Тем самым, объединяя результаты исследований за много лет, получаем общую картину термодинамической поверхности плотного калия в широкой области температур и давлений.

ЗАКЛЮЧЕНИЕ

Исследования свойств щелочных металлов в широкой области температур и давлений − одна из ярких страниц в истории института. В настоящей статье в качестве примера рассмотрена работа по определению плотности жидкого калия на линии насыщения: создание уникальной установки, анализ достоверности опытных данных при сравнении с результатами новых независимых измерений, использование специфических особенностей данных на линии насыщения и на изобаре в виде идеального КТР для зондирования линии инверсии эффекта Джоуля–Томсона плотного калия от 100 до 2000 атм.

Применение пьезометров постоянного объема и переменной массы для измерений плотности жидкостей и газов при высоких температурах и давлениях, используемых в течение многих лет в МЭИ и ИВТАН, начиная с работ Кириллина В.А. 1953 г., показало свою высокую работоспособность. В начале 1970-х гг. Мозговой А.Г. в ИВТАН под руководством Шпильрайна Э.Э. и Якимовича К.А. создал оригинальную установку для измерения плотности жидкого калия высоких параметров, реализующую метод капающего пикнометра. При этом вместо выпуска вещества в кювету и ее взвешивания при обычных условиях на установке ИВТАН проводилось взвешивание самой ампулы с металлом на аналитических весах, находящихся внутри пьезометра под давлением. Управление весами осуществлялось “через стенку” с помощью сельсинов. По сведениям авторов, эта установка до сих пор не имеет аналогов.

Экспериментальные данные ИВТАН о плотности жидкого калия на линии насыщения, опубликованные в 1976 г., в интервале температур 347–1650 К согласуются в пределах ошибок экспериментов с результатами новых независимых наблюдений 2000 г. Однако в свое время Мозговой А.Г. при защите диссертации отказался от результатов своих измерений плотности калия при более высоких температурах (до 2030 К), заменив их некоторыми обобщенными данными. Как оказалось, новые экспериментальные данные подтвердили справедливость такого решения.

Данные о плотности жидкого калия на линии насыщения имеют характерную особенность в виде идеального КТР в районе 1645 К. Отношение температуры этой особенности ${{T}_{{\text{э}}}}$ к критической температуре калия равно приблизительно ${{{{T}_{{\text{э}}}}} \mathord{\left/ {\vphantom {{{{T}_{{\text{э}}}}} {{{T}_{c}} = 0.75}}} \right. \kern-0em} {{{T}_{c}} = 0.75}}$ и близко к подобной величине для цезия. Это значение отличается от аналогичной оценки ${{{{T}_{{\text{э}}}}} \mathord{\left/ {\vphantom {{{{T}_{{\text{э}}}}} {{{T}_{c}} = 0.80 \pm 0.02}}} \right. \kern-0em} {{{T}_{c}} = 0.80 \pm 0.02}}$, сделанной ранее для большой группы неорганических веществ [7]. Дальнейшего анализа заслуживают вопросы принципиальности характера отличия оценок и специфики щелочных металлов.

Результаты измерений плотности жидкого калия, полученные в ИВТАН в 1976–1977 гг., являются базой для анализа других особенностей поведения этого вещества при высоких температурах и давлениях. Информация об идеальном КТР на изобаре 100 атм и данные о сжимаемости жидкого калия позволяют наметить поведение линии инверсии эффекта Джоуля–Томсона от $\sim {\kern 1pt} 1600$ К до ее пересечения с линией единичного фактора сжимаемости ${{z}_{{{\sigma }}}} = 1$ при температуре $\sim {\kern 1pt} 2300$ К и давлении ~2 × 103 атм. Объединением результатов прошлых лет и новых данных получается, хотя и очень приближенная, картина термодинамической поверхности плотного калия в широком интервале параметров, и дело за новыми экспериментальными и теоретическими работами, которые смогут уточнить эти оценки.

Авторы выражают благодарность сотрудницам библиотеки ОИВТ РАН и в первую очередь Л.М. Яковлевой за библиографическое обеспечение данной работы.

Список литературы

  1. Шпильрайн Э.Э., Якимович К.А., Тоцкий Е.Е., Тимрот Д.Л., Фомин В.А. Теплофизические свойства щелочных металлов / Под ред. Кириллина В.А. М.: Изд-во стандартов, 1970. 488 с.

  2. Шпильрайн Э.Э., Солдатенко Ю.А., Якимович К.А. и др. Экспериментальное исследование теплофизических и электрофизических свойств жидких щелочных металлов при высоких температурах // ТВТ. 1965. Т. 3. № 6. С. 930.

  3. Shpilrain E.E. Thermodynamic Properties of Fluid Metals at Medium and Low Pressures // Experimental Thermodynamics. V. II. Experimental Thermodynamics of Non-reacting Fluids. London: Butterworth, 1968. P. 953.

  4. Шпильрайн Э.Э., Якимович К.А., Мозговой А.Г. Экспериментальное исследование плотности жидкого калия при температурах до 1800 К и давлениях до 100 атм // ТВТ. 1976. Т. 14. № 3. С. 511.

  5. Шпильрайн Э.Э., Якимович К.А., Мозговой А.Г. Экспериментальное исследование плотности жидкого калия при температурах до 1740–2030 К // ТВТ. 1977. Т. 15. № 5. С. 1104.

  6. Yakimovich K.A., Mozgovoi A.G., Dubinin V.V. The Density of Liquid Potassium at Near-critical Temperatures // High Temp.‒High Press. 1979. V. 11. P. 543.

  7. Фокин Л.Р. Идеальный коэффициент термического расширения на линии насыщенной жидкости // Мониторинг: наука и технологии. 2012. № 2. С. 84.

  8. Фокин Л.Р. Соответственное состояние для коэффициентов термического расширения насыщенных жидкостей // ТВТ. 2012. Т. 50. № 3. С. 467.

  9. Шпильрайн Э.Э., Якимович К.А. Экспериментальная установка для исследования плотности жидких металлов // ТВТ. 1963. Т. 1. № 2. С. 173.

  10. Кириллин В.А., Румянцев Л.И. Экспериментальное исследование термодинамических свойств воды и водяного пара при высоких давлениях и температурах // Тр. МЭИ. 1953. Вып. XI. С. 20.

  11. Кириллин В.А., Шейндлин А.Е. Исследования термодинамических свойств веществ. М.–Л.: Госэнергоиздат, 1963. 510 с.

  12. Кириллин В.А., Зубарев В.Н. Экспериментальное исследование удельных объемов воды и водяного пара при сверхвысоких давлениях // Теплоэнергетика. 1955. № 11. С. 19.

  13. Кириллин В.А., Улыбин С.А. Экспериментальное исследование сжимаемости воды и водяного пара при температурах, близких к критическим // Теплоэнергетика. 1958. № 4. С. 53.

  14. Вукалович М.П., Зубарев В.Н., Александров А.А. Экспериментальное исследование удельных объемов водяного пара до температур 700–900°С при давлениях до 1200 кГ/см2 // Теплоэнергетика. 1962. № 1. С. 49.

  15. Зубарев В.Н., Прусаков П.Г., Барковский В.В. Экспериментальное определение удельных объемов водяного пара при температурах 673−873 К и давлениях до 200 МПа // Теплоэнергетика. 1977. № 8. С. 77.

  16. Вукалович М.П., Алтунин В.В. Экспериментальное исследование зависимости p–v–t углекислоты // Теплоэнергетика. 1959. № 11. С. 58.

  17. Покрасин М.А., Рощупкин В.В., Фокин Л.Р., Хандомирова Н.Э. Интерполяционные уравнения и таблицы давлений насыщения натрия, калия, рубидия и цезия // Теплофизические свойства веществ и материалов. 1983. № 19. С. 33.

  18. Шпильрайн Э.Э., Якимович К.А., Сковородько С.Н., Мозговой А.Г. Плотность и тепловое расширение жидких щелочных металлов // Обзоры по теплофизическим свойствам веществ. 1983. № 6(44). 89 с.

  19. Shpilrain E.E., Yakimovich K.A., Fomin V.A., Skovorodjko S.N., Mozgovoi A.G. Density and Thermal Expansion of Liquid Alkali Metals // Handbook of Thermodynamic and Transport Properties of Alkali Metals / Ed. Ohse R.W. IUPAC. Boston: Blackwell SP, 1985. Ch. 6.3.

  20. Быстров П.И., Каган Д.Н., Кречетова Г.А., Шпильрайн Э.Э. Жидкометаллические теплоносители тепловых труб и энергетических установок. М.: Наука, 1988. 264 с.

  21. Мозговой А.Г. Термодинамические свойства жидкого калия при высоких температурах. Автореф. дис. … канд. техн. наук. М.: ИВТ АН СССР, 1985. 21 с.

  22. Shpilrain E.E., Yakimovich K.A., Mozgovoi A.G. Apparatus for Continuos Measurement of Temperature Dependence of Density of Molten Melts by Method of a Suspended Pycnometer at High Temperatures and Pressures // Compendium on Thermophysical Property Measurement Methods. V. 2. Recommended Techniques and Practices / Eds. Maglich K.D., Cezairliyan A., Peletsky V.E. N.Y., USA: Plenum Press, 1992. P. 601.

  23. Мозговой А.Г., Сковородько С.Н., Шпильрайн Э.Э., Якимович К.А. Растворимость инертных газов в жидких щелочных металлах. Эксперимент // Обзоры по теплофизическим свойствам веществ. 1987. № 6(68). 80 с.

  24. Шпильрайн Э.Э., Сковородько С.Н., Мозговой А.Г. Новые данные о растворимости инертных газов в жидких щелочных металлах // ТВТ. 2002. Т. 40. № 6. С. 891.

  25. Сковородько С.Н., Мозговой А.Г., Резенова Е.В. Растворимость газообразного гелия в жидкой ртути при температурах 600−1500 К // ТВТ. 2011. Т. 49. № 2. С. 138.

  26. Hensel F., Hohl G.F., Schaumleffer D., Pilgrim W.C. Empirical Regularities in the Behavior of the Critical Constants of Fluid Alkali Metals // Z. Phys. Chem. 2000. V. 214. P. 823.

  27. Junst S., Knuth B., Hensel F. Observation of Singular Diameters in the Coexsistanse Curves of Metals // Phys. Rev. Lett. 1985. V. 55(20). P. 2160.

  28. Vargaftik N.B., Kozhevnikov V.F., Ermilov P.N. Experimental p-v-T Data end Equation of State of Liquid Cesium up to 2000 K and 600 bar // High Temp.‒High Press. 1984. V. 16. P. 233.

  29. Чернеев Е.И. Измерение параметров критической точки калия // Тр. Энерг. ин-та им. Г.М. Кржижановского. 1977. Вып. 61. С. 58.

  30. Guggenheim E.A. The Principle of Corresponding States // J. Chem. Phys. 1945. V. 13. P. 253.

  31. Martynyuk M.M., Balasubramanian R. Equation of State for Fluid Alkali Metals: Binodal // Int. J. Thermophys. 1995. V. 16(2). P. 533.

  32. Шпильрайн Э.Э., Якимович К.А., Мозговой А.Г. Изотермическая сжимаемость жидкого калия при высоких температурах // ТВТ. 1976. Т. 14. С. 1099.

  33. Apfelbaum E.M., Vorob’ev V.S. The Wide-range Method to Constract the Entire Coexsistence Liquid-gas Curve and to Determine the Critical Parameters of Metals // J. Phys. Chem. B. 2015. V. 119. P. 11825.

Дополнительные материалы отсутствуют.