Разработка справочных данных о коэффициентах переноса паров щелочных и щелочноземельных металлов на основе теоретических расчетов и анализа результатов экспериментов тема автореферата и диссертации по физике, 01.04.14 ВАК РФ

Мухтаров, Эркин Саидович АВТОР
кандидата технических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
1992 ГОД ЗАЩИТЫ
   
01.04.14 КОД ВАК РФ
Автореферат по физике на тему «Разработка справочных данных о коэффициентах переноса паров щелочных и щелочноземельных металлов на основе теоретических расчетов и анализа результатов экспериментов»
 
Автореферат диссертации на тему "Разработка справочных данных о коэффициентах переноса паров щелочных и щелочноземельных металлов на основе теоретических расчетов и анализа результатов экспериментов"

МОСКОВСКИЙ ордена ЛЕНИНА и ордена ОКТЯБРЬСКОЙ РЕВОЛЮЦИИ ЭНЕРГЕТИЧЕСКИЙ ИНСТИТУТ

РАЗРАБОТКА СПРАВОЧНЫХ ДАННЫХ О КОЭФФИЦИЕНТАХ ПЕРЕНОСА ПАРОВ ЩЕЛОЧНЫХ И ЦЩЯОЧНОЗЕЛЕЯЬШХ МЕТАЛЛОВ НА ОСНОВЕ ТЕОРЕТИЧЕСКИХ РАСЧЕТОВ И АНАЛИЗА РЕЗУЛЬТАТОВ ЭКСПЕРИМЕНТОВ

Специальность 01.04.14 Теплофизика и молекулярная физика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

На нравах рукописи

\

ЫУХТАРОВ ЭРКИН 0АВД03ИЧ

n

МОСКВА - 1992

Работа выполнена на кафедре Инженерной теплофизики Московского енергетического института

Научный руководитель: доктор {ызико-математичесюи наук А.Ы. Сел,=на$

Официальные оппоненты: доктор технических наук, профессор К.Ы. Лр&рьев

доктор тзхшгческих наук, ведущий научный сотрудник М-Р. Фошн

Ведущая организация: Московский авиационный институт

Зкцзтв диссертации состоится • "_" октября 1992 года в

_ часов но заседании специализированного совета К 053.16.02 при

¡Московском энергетическом институте по адресу: Москва, Красноказарменная ул., д.17, корпус Т, 2 этаж, к.206. (каф. ИТФ).

С диссертацией можно ознакомиться в библиотеке МЭИ.

Отзывы на автореферат в двух екземлярах просим направлять по адрес}": 105835 ГСП, Москва Е-250, Красноказарменная ул., д. 14, Учений совет Ю'Л.

Автореферат разослан "_" _ 1992 г.

Ученый секретарь специализированного совета К 053-16.02 к.ф.-м.н., с.и.о.

В.И. Ниш

/. ,' Г~ёTrn \

—1 ' ' j.wes-rnrwSi . _

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность теми. В последние годы в связи с развитием различного рода специальных установок растет интерес к применешно паров щелочных и щелочноземельных металлов, которые могут использоваться как в качестве теплоносителей (тепловые труби в теплонапря-жешшх двигательных установках ), так и в качестве , рабочих тел (стабилизирующие кошоненты в элементах конструкции ядерных реакторов; лазеры на парах металлов ). Интерес к названным веществам определяется их высокой теплопроводаостью, относительно низким давлением паров при высоких тешературах, больной теплотой парообразования, наличием в газовой фазе обратимых реакций типа диссоциации (щелочные металлы), слабой связью валентных электронов в атомах и молекулах. Литий к тому же является замедлителем быстрых нейтронов и источником получения трития, что определяет перспективы его использования в термоядерных реакторах.

Проектирование подобних устройств делает одним из важных вопросов современной энергетики получение надежных данных о теплофи-зических свойствах указанных веществ и, в частости, коэффициентах вязкости и теплопроводности. Расширение банка данных о коэффициентах переноса является семо по себе вашей задачей.

Среда наиболее известных таблиц коэффициентов вязкости и теплопроводности паров лития и натрия »ложно отметить данные Э.Э.Шпильрайяа, К.А. Якимовича и др., а таете таблицы РСД группы Н.Б. Варгафтика, B.C. Ярпша с сотрудниками, составленные в 1970 и 1984 г.г. с учетом имевшихся на то время опытных данных и теоретических оценок. N

К настоящему времени помимо новых опытных данных о вязкости и теплопроводности для лития и натрия опубликованы результаты кван-товомехаяических расчетов потенциальных криых взаимодействия атомов. Имеется также информация о потенциалах взаимодействия, полученная из данных спектроскопических экспериментов. Наличие точных потенциалов позволяет рассчитать коэффициенты переноса паров лития и натрия с погрешностью, не выходящей за рамки экспериментальных ошибок. Поэтому для паров лития и натрия на сегодня .имеется воз-мояюсть получить на основе анализа экспериментальных и расчетных данных надежные таблицы коэффициентов вязкости и теплопроводности.

ьпелне отгдчйкгу.е требованиям, предъявляемым к справочным данный.

Д.;:я Ейров щелочноземельных металлов в литературе опубликованы оксш-.рнмепталь.ние данные только о вязкости к теплопроводности пара магния, а тав;:е вязкости пара кальция, Имеетиеся в литературе дан-кие о цзаимодейстыш атомов металлов второй группы позволяют провести расчет коэффициентов вязкости и теплопроводности паров указанных веществ по формулам строгой кинетической теории газов.

"Работы. 1.Выполнить теоретически;! расчет коэффициентов переноса паров лития и натрия на основе имевдихся в литературе кьантоБомеханкческлх и спектроскопических данных о взаимодействии атомов.

2.Сопоставить результаты расчетов и имеющиеся экспериментальные данные на основе согласования расчетных и опытных даншх равработьть справочные данные о коэффициентах переноса паров лития к натрия в широком диапазоне параметров состояния.

3..Выполнить теоретический расчет ко&фХкциентов переноса паров металлов второй группы периодической системы элементов (Бе, Е^, Са, Бг, Ва, 1п, Сй).

Основные результаты и их научная новизна.

1.В диапазона температур 800-2000 К на основе современных литературных данных о потенциалах взаимодействия атомов для лития и натрия рассчитаны все усредненные сечения столкновения мономер-мономер и нокомер-дпмер, необходимые для расчета коэффициентов переноса во втором приближении теории Энекога-Чепмеиа.

2. Ка основании теоретического расчета вязкости н теплопроводности паров лития и натрия и сопоставления его результатов с имеющимися експериментальнкми данными делается вывод о целесообразности использоваш!я квантовомеханичеких и спектроскопических данных о взаимодействии атомов и возможности получения на етоЯ основе значений коэффициентов переноса указанных веществ, практически не уступающих по точности теплсфизкческшу вксперименту.

3. На основе согласования вкспериментальных и расчетных данных о коэффициентах вязкости и теплопроводности разработаны уточненные таблицы справочных дашшх об указанных коэффициентах переноса паров лития и натрия в диапазоне температур 800-2000 К и давлений 0,1-1000 кПа.

4. Впервые получены расчетные данные об усредненных сечениях

столкновений атомов и коэф|ядиентах вязкости а теплопроводности паров металлов второй группы: Ве, Са, Бг, Рл, 7.п к С<Д.

Практическая значимость работы.

1. По материалам диссертации.. подготовлены 'стандартные справочные данные о теплопроводности и вязкости паров лития и натрия, аттестованные ВНИЦ МВ Госстандарта РФ. Данные о койф5ищ;ентах переноса паров металлов второй группы используются ■ предприятием-заказчиком. Таблицы теплопроводности и вязкости перечисленных веществ переданы в банк данных Научно-информационного центра тепло-физических свойств чистых веществ (топлофизичеекий центр) КВТРАН.

2. Разработана программа расчета коэффициентов вязкости и теплопроводности химически реагирующих газовых смесей, которая может быть использована, в частности, в расчетах указанных величин для паров других щелочных металлов, а такие эффективная программа статистической оценки экспериментальных данных нелинейным методом наименьших квадратов.

Апробация работа. Основные результаты работы докладывались на на заседании Секции теплсфизических свойств вецеетв Научного Совета "Теплофизика и теплоэнергетика" АН СССР (1990 г.); на заседай®! рабочей группы по свойствам диссоцилруивдх газов Советской комиссии по теплофизическим таблицам'газов и жидкостей технически важных веществ (Минск, 1991 г.); на семинаре по методам оценки достоверности физико-химических данных, проводимого секцией теплофизических свойств Научного Совета "Теплофизика и теплоенэргегшеа" РАН совместно с секцией "Банк данных "Металлургия" Научного Совета по физико-химическим основам металлургических процессов" (1991 г.): а также на научных семинарах в МЭИ, >ДАИ и ИВГРАН.

Публикации. По теме диссертации опубликованы 4 статьи [1-4].

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения, списка литературы и приложения. Работа изложена на 127 страницах, из них 12 с рисунками. 25 с таблицами. Библиографический список содержит 118 наименований.

СОДЕРЖАНИЕ РАБОТЫ

Е<2 введении обосновывается необходимость исследований, проведенных в диссертации, даются краткий обзор состояния рассматриваемых вопросов, общая постановка решаемых в диссертации задач, о

также перечисляются основные результаты, выносимые на защиту.

Б первой главе анализируются литературные данные о потенциалах взаимодействия атомов и атомов с молекулами в парах лития и натрия. Отмечается, что при температурах Т£ 2000 К атомы щелочных металлов находятся в основном состоянии 2Б. В этом состоянии два атома могут взаимодействовать по двум потенциальным кривым: син-глетной Ф^'Чг), описывающей притяжение между атомами, которое • обусловливает воеможность образования димера Ме , и триплетной

»■п. 3

Фд (г), описнващей отталкивание между атомами на малых и средних расстояниях г и притяжение на больших.

Проводится сравнение параметров синглетной потенциальной кривой, полученных "в результате вариационного решения уравнения Шре-дангера и из спектроскопических экспериментальных данных об елект-ронно-колеОательно-вращательных спектрах, излучаемых димерами при переходах из возбужденных электронных состояний в основное. Срав-некие показывает, что погрешность теоретических потенциальных кривых по межатомным расстояния^ находится на уровне 0,1-0,01 1, а по энергиям (в температурных единицах) не превышает 100-200 К (1-3& от онергии диссоциации Р ). Хотя экспериментальная информация о потенциалах отсутствует, можно ожидать, что точность теоретических триплетных потенциальных кривых, рассчитанных тем же методом, что и «шглетних, хотя и на основе более грубых пробных волновых функций, не намного хуже.

Потенциалы взаимодействия атомов лития и натрия, описанные в разделе 1.1 используются в разделе 1.2 при вычислении потенциалов взаимодействия 4>1а димера Ые2 с межатомным расстоянием г и мономера Ме, находящегося на расстояниях г и тгд от атомов 1 и 2, "связанных" в димер. Расчет потенциала Ф проводится,исходя из потенциальных поверхностей трех атомов Ф^21\ определяемых методом двухатомных комплексов в молекулах на основе комбинирования двухатомных потенциалов Ф{1> и ©(э> э з

Ф1В(г.Э;г12) = - Ф»>(г1а) , (1)

где г - расстояние от мономера 3 до центра димера 12; 3 - угол между осью димера и прямой, соединяющей мономер и центр димера.

Результаты главы 1 используются во второй главе для расчета усредненных сечений столкновений атом-атом и атом-молекула (для

лития и натрия) 0[*,И>(Т). фигурирующих в соотношениях кинетической теории газов, основанных на решении уравнения Больцмана методом Энскога-Чепмена.

Сечения столкновений атомов щелочных металлов могут бить найдены усреднением вкладов от с:шглетного и триллетного взаимодействий:

1 3

q"-'1 = _ q<«,.m> + 1 g<i..)(3> (2)

4 ■ 4

Вычисление сечений выполнялось по программе, алгоритм которой основан на квадратурной формуле Гаусса-Чебышева. Он описан в работе A.M. Семенова, А.О. Еркимбаева (1986 г.) и дополнен в настоящей работа возможностью интерполяции таблично заданных потенциальных кривых кубическими сплайнами.

Выбор различных исходных данных сравнительно слабо влияет на рассчитанные значения сечений î расхождения не превышают 4%, то есть находятся в пределах погрешности обсуждаемых расчетов. Результаты расчетов находятся в приемлемом согласии с данными, полученными другими авторами: В.Л. Белов , H.H. Ключников (1965 г.): Б.М. Смирнов, М.11. Чибисов (1971 г.); B.C. Ярган B.C., H.A. Вани-чева, В.И. Долгов (1980 г.); Э.Э. "Шпильрайн Э.Э., А.Я. Полшцук

(1979 г.): В.И. Долгов (1984 г.); A.M. Семенов, А.О, Еркимбаев (1986 г.): Holland P.M., Blolsi Ь. (1936, 1987 гг.), Niestro de Castro С.A., Parelerla J.M.N.A. (1990 г.).

При рассмотрении столкновений мономера с димером главной трудностью на пути определения сечений столкновений можно считать несферичность их взаимодействия ( частота неупругих' столкновений с изменением колебательно-вращательных состояний обычно на порядок меньше частоты упругих столкновений). Можно показать, что отношение времени г , за которое димер успевает сделать один оборот, ко

Эр

времени столкновения тс составляет »1.4+-2, т.е. димер за время одного столкновения не успевает сделать один оборот. И для учета несферичности взаимодействия предлагается следующая процедура: полагая, что во время столкновения ориентация димера не изменяется, вычисляем Q^' "'(Г;!* для потенциала (1), в котором rJ2 и 9 рассматриваются как фиксированные параметры. Затем производим усреднение результата по всем ориентациям, считая их равновероятными, и по всем возможным межатомным расстояниям в димере г , ис-

пользуя статистическую сумму по связанным состояниям двух атомов. Погрешность рассчитанных величин ' ° 1 как ожидается, порядка 20Ж, хотя ота оценка скорее интуитивна. Проведено сопоставление этих величин с результатами предшествующих работ (Э.Э. Шгшльрайн, А.Я. Полищук (1979,1980 г.)).

Третья глава посвящена теоретическому расчету коэффициентов переноса паров лития и натрия, а также разработке на основе анализа и обработки результатов расчетов и экспериментальных данных справочных таблиц о коэффициентах вязкости и теплопроводности указанных веществ.

Расчет выполнен в предположении, что пары указашшх веществ

представляют собой бинарную смесь мономеров Ы., На и димеров Ыг,

На„ , реагирующих по схема 2

2Пв ^ Ме2 (3)

Предполагается, что неупругио столкновения Не + Не . Не, + Ме„ и

3 2 3

реакция (3) то влияют на распределение частиц по скоростям, так что коэффициент вязкости и вклад в коэффициент теплопроводности от переноса тепла поступательными степенями свобода молекул Хпост мокно рассчитать по формулам для смесей разреженных одноатомннх газов, а вклады в коэффициент теплопроводности от переноса тепла внутренними степенями свобода молекул Хг,тт и от переноса теплового аффекта реакции Хр - по формулам Гиршфельдера и Гиршфельдера-Ерокау: 1=г?пост, *=*поот + >вн + • Неточности в результатах расчета, связанные с этими допущениями, по-видимому, не превосходят онибок, вызвашшх неопределенностью исходных дашшх.

Хотя в парах лития и натрия, помимо димеров Не2, имеются три-меры Не3 и тетрамеры Ме4, их концентрации столь малы, что влиянием вклада стих частиц в величины ч п 1 мозно пренебречь в пределах

погрешности расчета. Константа равновесия реакции (3) К ,

р

необходимая для расчета состава смеси, вычислялось на основе синглетного потенциала ыезкзтошюго взаимодействия (г).

Выполненные на основе наилучших квантовомеханических потенциалов межатомного взаимодействия и внчислешшх в главе 2 "квантозо-кехыглческих" сечений столкновений частиц расчеты коэффициентов вязкости и теплопроводности паров лития и натрия в интервалах температур 800зТг2000 К и давлений 102«р*10в Па (М) и 10э*р510в Па (Ыа) приведены в (33 и представлены на рис.1,2.

т,к

т,к

Рис.1. Вяпкость t) и теплопроводность к пара лития.

3 2.5 2 1.5

В00 1000 1200 1400 1600 1800 2000

.. т,к

0.05 0.04

0.03 0.02

800 1000 1200 1400 1600 1800 2000

т,к

Рис.2. Вязкость г) и теплопроводность X пара натрия.

77-ю , Па-С

I I I I

/уЛ '•7 \

-50 1С|° ЮООкПа ) I I I

• На основе оценок основных составляющих погрешности теоретического расчета, проведенных в работе, погрешности результатов оцениваются следующим образом: в наилучших условиях (низкие давления, высокие температуры) в наихудших условиях (высокие

давления, низкие температуры) бг}«10-15Й, 6Х«20-30%.

В разделе Э.1 показано,что результаты расчета коэффициентов в пределах погрешности согласуются с данными большинства (11 из 13) опубликованных работ, в которых измерялись коэффициенты переноса паров лития и натрия. Еще лучшее согласие достигается в случае использования в расчетах "спектроскопических" сечшшй столкновений

и К , вычисленной на основе синглетной потенциальной кривой, кото-р

рая получена из данных спектроскопических экспериментов. Результаты этого раздела позволяют надеяться, что путем некоторой корректировки параметров теоретического расчета удастся достигнуть оптимального согласия вычисленных на его основе значений коэффициентов переноса с имеющимися экспериментальными данными.

Вопрос о правильной корректировке параметров теоретического расчета решается при сопоставлен;« полученных в главе 2 расчетных данных о сечениях столкновений частиц с сечениями, определенными путем раздельной обработки экспериментальных данных каждого -из экспериментов о ковффициентах переноса нелинейным методом наименьших квадратов. Такое сопоставление проведено в разделе 3.2 и представлено на рис.3,4, которые наводят на мысль, что разумный компромисс между результатами цитированных выше исследований, требуемый для создания широкодиапвзоиных справочных данных о коэффициентах, переноса паров лития и натрия, может быть достигнут, если выполнить расчет этих коэффициентов переноса по упомянутой вше методике [4] и на основе тех же исходных данных, лишь немного скорректировав вычисленные в главе 2 температурив зависимости сечений атом-атом (спектроскопических) и атом-молекула 11:

где и чг- константы, близкие к 1, а комбинации сечений "атом-молекула" А)а и В1г принимаются в соответствии с описанным вше теоретическим расчетом, равно как и сечение "молекула-молекула" ^аз'а)• Беличина которого ( й(2,г)/П = 50 А2 ) практически не вли-

21 19

17 10

13

800 1000 1200 1400 1600 1800 2000

г\и, 2) -г

N а

а

* ! 1 т.к

80

60

40

20

Г" 5Г 1 2!-

3

1 е

*

Т.К

6i-

boo 1000 1200 1400 1боо 1800 2000

Рис.3. Сечения столкновений мономер-мономер q'^'31

'------ - и мономер-дамер 5{g'n в паре лития

1-данная работа. Эксперимент! 2 - В.И. Долгов (г) ,1984 г.)! 3 -И.Ф. Степаненко, Н.И. Сидоров, Ю.В. Тарлаков, B.C. Яргин ( г),1906 г.); А- В.М. Капитонов (Х.1983 г.): 5- Д.Л. Тимрот. В.В. Махров, Ф.И. Пильненъский (Х.1983 г.); 6- Ю.К. Якимович (X, 1987 г.).

Сплошная линия - расчет ио кгантовомеханическим данным; штриховая-по спектроскопическим донным

Рпс.4. Сечения столкновений мономер-мономер

и мономер-димер Q J ^'1' в паре натрия

1- данная работа. Эксперимент: 2 - И.Ф. Стегтненко, H.H. Сидоров,

Й.В. Торлаков, B.C. Ярган (т),193б г.): 3- H.H. Сидоров, Ю.П. Тар-

лаков, B.C. Ярпш (г?,1976 г.): 4- Д.Л. Тимрот, А.Н. Варава (п, .1977

г.); 5- Н.Б. Варгафтик, A.A. Вощшпш (Х,19б7 Г.): б- Н.Б. Bopra-Jh-

тик, A.A. Вощишш, В.В. Керженцев (А,1973 г.); 7- Д.Л. Тимрот,

В.В. Махров, В.И. Свириденко(Х,1976 г.); 8-В.В. Махров, Ф.И. Ппль-

непьский (X, 1983 г.); Д.Л. Тимрот, Е.Е ТощсиЯ (Х.1967 г.)-9

Сплошная линия - расчет по квантовомехапическим данным! ттриховая-по спектроскопическим дагашм

яет на результаты расчетов вязкости и теплопроводности рассматриваемых паров. Теплота диссоциации дачера АН принимается по результатам наиболее тщательно выполненных спектроскопических измерений: ЛН(Ыа)=102,89±0,2, ДН№2)=70,9±1,3 кДж/моль.

Методика расчет, таблиц справочных ванных. Для оценки опти-л л

мальных параметров д( и цг использована следующая процедура. Сначала методом наименьших квадратов раздельно для каждой из экспериментальных работ оценивались температурные зависимости сечений столкновений и 1 '(Т). Затем оптимальные параметры ч

и находились усреднением этих "експериментальных" (рис.3.4) сечений (рассматриваемых как случайные функции) путем минимизации по параметру функции

Т<к)

| £ | ^"Чт) - ^-»(т)]3®. (5)

а 1 т(к>

Здесь диапазон температур, в котором измерены коэффици-

енты переноса к-го эксперимента; оценка сечений, соответст-

вующая етим данным. В таком подходе все опытные данные имеют равное представительство в обработке( в традиционно применяемой совместной обработке методом наименьших квадратов при наличии систематических погрешностей данные большего объема приобретают пропорционально большее влияние на результаты обработки), а также появляется возможность учесть априорную информацию о расчетных сечениях столкновений путем включергая их в число усредняемых "экспериментальных" сечений столкновений. В качестве теоретических использовались вычисленные в данной работе "спектроскопические" сечения, как согласукщиеся с расчетами других авторов.

"Веса" в (5) или полагались равными 1, или уточнялись итерациями [4] по формуле

V ( 1 + вк / [£У(»-1>} г1 •

но оценки параметров при этом слабо менялись.

Обработка данных, основанная на усреднении (5), проводилась в несколько шагов. На каждом шаге после корректировки вычисленными множителяами <з и дз по формуле (4) сечений столкновений на их основе рассчитывались коиЗфициенты переноса. По распределению откло-

нений экспериментальных значений г) и X от вычисленной сглаживающей зависимости проверялась взаимная согласованность данных о вязкости и теплопроводности. При этом использовался критерий Фишера.> При необходимости из обработки исключались наиболее сильно уклоняющиеся от усредняющей зависимости опытные данные. Но следующем шаге процедура обработки и анализ согласованности повторялись для оставшихся данных.

Л А А

Результаты такого подхода Ш- 99. 4=1,49; Иа- 4=0,95.

А 13-1

4=0,62) использованы для расчета справочных данных о коеффщиентах переноса паров лития и натрия [4].

Проведенные .в диссертации анализ результатов"обработки и сравнение с другими справочными данными позволяют надеяться, что получении е в работе таблицы значений коэф1ади9нтов переноса пара лития и натрия наиболее надежны по сравнению с имеющимися в настоящее время литературшдм данными.

В четвертой главе описиваотся расчет коэффициентов вязкости и теплопроводности паров щелочноземельных металлов, а также полученная корреляция между усредненными сечениями столкновений атомов второй группы периодической системы элементов.

Атомы щелочноземельных металлов Ве, Mg, Са, Бг, Ва, находящиеся в основном электронном состоянии взаимодействуют друг с другом по единственной потенциальной кривой Ф(г). У всех рассмат-ргваемых веществ на этой потенциальной кривой имеется неглубокий минимум , отвечающий возмокности образования димеров Моа. Однако концентрация димеров в парах металлов второй группы (в отличие от щелочных) ничтожна, и при проведении теллофизических расчетов эти пары можно рассматривать как одноатомные.

В настоящее время наиболее надежным источником информации о взаимодействии атомов щелочноземельных металлов являются экспериментальные данные об влектронно-колебательно-вращателышх спектрах. Сечения столкновений (Т) атомов щелочноземельных металлов, необходимые для определения коэффициентов вязкости и теплопроводности, вычислялись в интервале температур 800-2ООО К, исходя из Ф(г), по той же методике, что применялась для расчета сечений столкновений атомов щелочных металлов. Рассчитывались также значения второго группового интеграла Ьа, позволяющие оценить константу равновесия К реакции и тем самым - доля димеров х„ в парах щелоч-

Р я ,

поземельных металлов. При Т=ЮОО К и р=103 Па (что заведомо выше давления насыщениях паров) х оказываются равными 0.03& для Ве, 0,1,4 для М£, 0,6Я для Са и 0,6^5 для Зг, а при меньших давлениях и более высоких температурах и подпвно пренебрежимо малыми.

Результаты выполненных расчетов вязкости и теплопроводности паров щелочноземельных металлов [1] в настоящее время удается сравнить о экспериментом только для магния и кальция. Погрешность опубликованных для пара магния результатов измерений вязкости и теплопроводности, по оценкам авторов, ответственно 4 и 10%. Данные хорошо согласуются и езду собой: соотношение Х/г)=15к/4га (к -постоя!иая Больцмана, т - масса атома), справедливое для одноатомных газов, выполняется с погрешностью «ЗЙ. В то же время измеренные значения г) и X вдвое превышают величины, приведенные в [1], а сечения ц'3,2', извлекаемые из экспериментальным данных, соответственно вдвое шике теоретических сечений.

В работе показывается, что объяснить эти расхождения неточностью исходных данных невозможно, аце меньше оснований искап причины этих расхождений в оиибках численной процедуры расчета сечений. Наконец, об отсутствии в использованной схеме расчета значительных ошибок как вычислительного, так и методического характера убедительно свидетельствует вполне приемлемое согласие с ¡экспериментом результатов выполненного в главе 3 по той же схеме расчета для более сложного случая пароз щелочных металлов.

Между тем, опубликованные после выполненных нами расчетов результаты измерений реактивным методом в интервале температур 1710-2020 К вязкости пара кальция, погрешность которых авторы (Д.Л. Тимрот, Б.Ф. Реутов, О.В.Кречетов) оценили в 656, почти совпадают с расчетными значениями г): различие 3-4%.

В разделе 4.2 показывается., чсо в литературе отсутствуют надежные данные о потенциалах взаимодействия атомов Ва, Ъп и Сй. Малая продуктивность попыток оценить коэффициенты переноса, используя простые, но грубые квантовомеханические модели потенциала взаимодействия, обусловливает поиск каких-либо корреляций, связывающих коэффицентн переноса паров металлов второй группы периодической системы и позволяющих рассчитать вязкость и теплопроводность Ва, гп и са.

' Корреляция сечений спомоюЯенхй а-.о.юб мшусб (торсй С целью найти указадпу» корреляцию и уточнить на ее осксве рассчя-ташше в работе на основе пстеидозлов ьэогалодеАстевя Морзе (аппроксииируиадх имеющиеся в литературе данлн» о ьз&имоя'.'Ястйии атомов Ва, СЛ и 2я) значения сочкш.1 столоторовий отомов Вч, со. и Яп будем исходить из того, что естеотшяшт нзеитг-Зом, оародоля»-ккм величину сечений в рлду ьтомов со скс.ным:: гаотгожк эл-'ктрок-вими оболочками, является квчдр'зт слетерэр-ского рш-тус-ч атома г., который характеризует эффективный размер электронно.'*, оболочки. Этот параметр позволяет ввести безроз-^-рнчо сечения /ла2.

Простейшее предположение состоит г> тог.;, что д.-я" всех "подо:них" веществ (а необходим™ условием подобия является, очевидно, сходство внеиних электронных оболочек) зависимость 0* от безразмерной температуры Т*= кТ/1 (I - поюгак&л копизэщга атсмо), учитывавшей различие энергетических характеристик атомов, окзкется одно.1 и той же.

В действительности одинаковость структуры впеглнх ?лектрс;л:ых оболочек всех рассматриваемых атомов являетел необходимым, но отнюдь не достаточным условием подобия меа;ато"ного взаимодействия и, как следствие, теплофигического подобия соответствуют«: волости.

Рис.Корреляция безразмерных сечений столкновений атомов металлов вторсЯ группы 0* и их атошшх номеров (зарядов ядра) Ъ при различных приведенных температурах Т*.

Существенное влияние на характер межатомного взаимодействия, очевидно, должно оказывать строение внутренних электронных оболочек. Последовательное заполнение внутренних электронных оболочок дол-хно, по-видимому, приводить к закономерному изменению относительной интенсивности сил межатомного притяжения и к соответствующему изменению величин 0*.

Сказанное наводит на мысль попытаться проанализировать зависи-_ »

мость 0 от атомного номера I как параметра, характеризующего относительное положение соответствующих атомов в данной группе периодической системы, при одинаковых безразмерных температурах Т*. Из рис.5 видно, что таким способом действительно удается построить единую корреляцию <3*(Т* позволившую рассчитать усредненные сечения столкновений 5'3'а>(рио.б) и коэффициенты вязкости и теплопроводности в парах бария, цинка и кадмия при температурах 800-

2000 К.

й(2'2) 7Г , А2

б

зо

35

25

20

15

10

800

1200

1600

2000

т,к

Рис.6. Сечения столкновений атомов О второй группы

П<2,2>

в парах металлов

ОСНОВНЫЕ ВЫВОДЫ

1. На основе имеющихся в литературе кванговомехаяических и спектроскопических данных о взаимодействии атомов лития и натрия рассчитаны сечения столшовениП мономеров в парах указанных метал-

лов. Показано, что неопределенность исходник данных, вносимая неточностью потенциальных кривых взаимодействия, мало сказывается на значениях сечений столкновений.

С использованием метода "фиксированных ориентация" вычислены сечения столкновений атома с молекулой в парах лития и натрия.

2. Без привлечения какой бы то ни было эмпирической информации на основе "квантовомеханических" сечений столкновений атомных частиц вычислены коэффициенты вязкости и теплопроводности паров лития и натрия. Результаты настоящего расчета приводят к выводу о том, что в настоящее время для ряда веществ, круг которых бистро расширяется, чисто теоретический расчет по точности практически уже не уступает теплофйзачеекому эксперименту, а по широте охвата области параметров состояния, полноте и детальности доставляемой информации и, разумеется, по экономическим соображениям такой подход и вовсе вне конкуренции.

3. Проведены анализ взаимной согласованности экспериментальных-и расчетных датгых о вязкости и теплопроводности паров лития и пзтртя и их совместная обработка на основе известных формул строгой кинетической теорет газов. Предлог:езгшй в даосертацшг :тодход позволяет в тактильно ¡гирсксм дпзяззогге параметров состояния достигнуть оптимального согласия Г'ЭТхеленяих теоретическим путем коэффициентов переноса с кмсеезижся експер-иментальнимп данни/и и получить согласоважше на базе единой неделя оценки коэффициентов вязкости и теплопроводности иэрез литил v, натрия, вполне удовлетворяющие требованиям, предъявляете;?* к справочным данным. Указанные оценки с анализом основная еоставлящих обдах погрешностей предлагаемых данных положены в основу стандартных сггравоч1Шх данных о коэффициентах переноса исследуемых веществ.

4. Исходя из спектроскопически дашшх о взаимодействии атомов вычислены усредаенпие сечения столкновения атомов Q<3'я1 (?) и второй групповой интеграл для паров да-ючноэемельных металлов: Бе, Mg, Са, Sr. Нз основе сечений стаххновевяЯ по формулам для разреженных однозтомяых газов вьгчяслекн вязкость и теплопроводность указанных веществ. Для пзров беряялкя, кальция и стронция эти данные получены впервые.

5. Полученная sevascssiocn безразмерна! сечений столкновений от атомного номера (зерядэ) при одинаковых приведенных температурах

позволила рассчитать ксаф$ициеити вяслости и теплопроводности паров бария, цинка и кадмия.

Основное содержание диссертации опубликовано в следующих работах:

1.Мухтаров З.С., Семенов A.M. Расчет теплофизических свойств паров щелочноземельных металлов // Теплофизика высоких температур.- 1937.- Т.25.- Кб,- С.1094-1099.

2.Семенов A.M., Мухтаров Э.С. Расчет коэффициентов переноса паров щелочных металлов "из первых принципов"// Тр. ин-та / Московский энергетический институт.- 1990.- *236.- С.70-74.

3.Мухтаров Э.С., Семенов A.M. Невмпирический расчет коэффициентов переноса паров лития и нагрия // Теплофизика высоких температур.- 1990.- Т.28.- N1.- С,56-63.

4.Мухтаров Э.С.«Семенов A.M. К разработке справочных данных о коэффициентах переноса паров лития и натрия U Теплофизика высоких температур,- 1992,- Т.30.- N3.- С. ^-55.

WfJjQQ аим fail

™ Типшрафня MSII. Кр,и-»о|;аиряг№гя, U.