Самосогласованная микрополевая модель неидеальной плазмы




Скачать 182.53 Kb.
Название Самосогласованная микрополевая модель неидеальной плазмы
Дата публикации 17.05.2014
Размер 182.53 Kb.
Тип Автореферат
literature-edu.ru > Математика > Автореферат


ИНСТИТУТ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ

РОССИЙСКОЙ АКАДЕМИИ НАУК
На правах рукописи


Павлов

Алексей Сергеевич

Самосогласованная микрополевая модель неидеальной плазмы

Специальность 05.13.18 – математическое моделирование,

численные методы и комплексы программ
Автореферат

диссертации на соискание ученой степени

кандидата физико-математических наук

Москва

2005

Работа выполнена на механико-математическом факультете Московского государственного университета им. М.В.Ломоносова и в Институте математического моделирования РАН

Научный руководитель – член-корреспондент РАН Калиткин Н.Н.
Официальные оппоненты – доктор физико-математических наук,

профессор Змитренко Н.В.

– доктор физико-математических наук,

профессор Лебо И.Г.

Ведущая организация – Московский физико-технический институт (Государственный университет).

Защита состоится «____» _______________________________2005 года в _______часов на заседании диссертационного совета К002.058.01 при институте Институт математического моделирования РАН по адресу: 125047, Москва, Миусская пл., 4а.

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


Автореферат разослан «_____»________________2005г.

Ученый секретарь диссертационного совета

кандидат физико-математических наук Прончева Н.Г.

Общая характеристика работы
Актуальность темы. Многократная ионизация разреженных газов описывается системой классических уравнений Саха. Для нее доказано существование и единственность решения, а в [1] предложен эффективный метод решения, пригодный для реализации на ЭВМ.

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

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

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

Сначала строились классические дебаевские поправки ([4]). Однако большинство физиков считало, что эти поправки сильно завышают эффект неидеальности при высоких плотностях. (Экспериментов, с которыми можно было бы сравниться, не существовало.) К тому же в классической модели Дебая-Хюккеля нарушается электронейтральность среды в целом, а также наблюдается не единственность решений уравнений ионизационного равновесия. Столь серьезные недостатки вообще ставят под сомнение правомерность использования данной модели. Позднее стала популярной модель БДХ (Дебай-Хюккель в большом каноническом ансамбле,[5]). Она приводит к гораздо меньшим эффектам неидеальности; на наш взгляд она их сильно занижает. Эти и многие другие работы обосновывались различными приближениями теории возмущений с помощью фейнмановской диаграммной техники. Обзор это направления есть в [6].

Другим подходом, который мы считаем наиболее перспективным в настоящий момент, является микрополевая модель неидеальности ([7]). В ней использована идея плазменного микрополя, давно применявшаяся в спектроскопии. Было предложено рассчитывать термодинамические поправки на неидеальность с помощью этого микрополя. Эта модель существенно лучше согласовывалась с результатами экспериментов для очень плотной плазмы, проведенных в 1980-е годы, чем [4]-[6].

Однако оказалось, что микрополевая модель [7] требует доработки. Во-первых, использовавшееся там выражение для сдвига потенциала ионизации не позволяло написать простые выражения для поправок на неидеальность к термодинамическим

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

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

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

Практическая ценность работы. Разработан численный алгоритм для решения обобщенных уравнений Саха с микрополевой неидеальностью на основе которого построен программный комплекс «Плазма-элемент», реализованный на языке С++. Он позволяет производить расчет ионизации и термодинамических характеристик многократно ионизированной однокомпонентной плазмы для любого элемента периодической системы от водорода () до лоуренсия (). В качестве примеров расчетов в работе представлены таблицы относительной электронной концентрации, давления и энергии для плазмы лития, алюминия и меди.

Апробация работы. Полученные результаты докладывались и обсуждались на международной конференции «VII тематические харитоновские чтения» (Саров, 2005). По материалам диссертации сделаны доклады на совместном семинаре Института математического моделирования РАН и кафедры математического моделирования Московского физико-технического института (октябрь 2005), на кафедре общих проблем управления механико-математического факультета МГУ им. М.В.Ломоносова (октябрь 2005).

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения, списка литературы и приложения. Общий объем диссертации 76 страниц, рисунков 11, таблиц 11. Список литературы включает 47 наименований.

Публикации. По теме диссертации опубликовано 4 работы. Их список представлен в конце автореферата.
Краткое содержание работы.
Введение включает обоснование актуальности темы, содержит формулировку основных целей работы и краткое содержание глав.

В первой главе приводятся некоторые результаты, необходимые для дальнейших изысканий. В начале обсуждается модель ионизационного равновесия для плазмы. В рамках модели рассматривается газо-плазменная смесь i-х элементов с порядковыми номерами . Она считается состоящей из ионов разной кратности k этих элементов (k=0 означает нейтральные атомы, которые присутствует при невысоких температурах) и свободных электронов. Атомы и ионы рассматриваются как частицы с жесткой внутренней структурой, т.е. с системой возбужденных уровней, имеющих энергии относительно основного состояния частицы () и статистические веса . Для отрыва электрона от (k-1) -кратного иона надо затратить энергию, равную k-му потенциалу ионизации . Обозначим относительные концентрации электронов через , относительные концентрации ионов i-го сорта и кратности ионизации k - через . Они удовлетворяет следующим балансным соотношениям:

, , , (1)

здесь - доли различных элементов в смеси (по числу атомов). Они заданы начальным составом плазмы.

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

, , (2)

где есть поправка на неидеальность, а атомно-ионные и электронные члены соответственно определяются идеальногазовыми выражениями [Б4]

, . (3)

При этом атомы и ионы считаются классическими, для электронов учитывается частичное фермиевское вырождение. Здесь V есть средний объем атомной ячейки, T – температура, - потенциал k-ой ионизации i-го элемента, - статистическая сумма по всем возбужденным состояниям k-кратного иона i-го элемента, - функции Ферми-Дирака, - масса иона, измеренная в электронных массах и - химический потенциал электронов. Последний связан с параметрами вещества соотношением

. (4)

Минимизируя функционал (2)-(3) по всем концентрациям с учетом соотношений (1) и (4), получим систему уравнений ионизационного равновесия для всех допустимых индексов i, k

, (5)

где (6)

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

(7)

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

Для неидеальной плазмы поправка становится существенной; величины же даже при заметной неидеальности могут оказаться пренебрежимо малыми в некоторых диапазонах температур и плотностей вещества.

Далее в данной главе предлагается обзор наиболее известных моделей неидеальности, в частности широко известные и упомянутые ранее модели МДХ и БДХ. Приводится выражение для микрополевой поправки на неидеальность из [7]. В ней вклад ионов предлагается аппроксимировать следующим достаточно громоздким выражением:

, (8)

где , , , .

При низких температурах и больших плотностях значительна кон­центрация нейтральных частиц. Их влияние предлагается учитывать следующим образом:

(9)

Заметим, что сдвиг (9) построен в довольно грубом приближении, когда нейтральные частицы моделируются твер­дыми сферами. (Жесткая механическая система - отрыв электронов с внешней орбиты обусловлен столкновением с электроном другого атома). Окончательно поправка на взаимодействие, учитывающая оба вклада, ионный и электронный, выглядит следующим образом:

(10)

Вторая глава посвящена построению микрополевой поправки на неидеальность и численного алгоритма для решения системы уравнений ионизационного равновесия с микрополевым сдвигом потенциалов.

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

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

Рис. 1 Атом (ион) во внешнем однородном поле; разрез по направлению поля. Сплошные линии – невозмущенные уровни, пунктирные – возмущенные.



Значение зависит от функции распределения микрополя. Например, если взять распределение Хольцмарка, то

,. (11)

Тогда получим следующее выражение для сдвига потенциала ионизации:

(12)

Удалось подобрать несложное выражение для , чтобы соответствующий ему сдвиг потенциала ионизации достаточно точно аппроксимировал (12):

. (13)

Для него сдвиг потенциала к-й ионизации равен

(14)

Здесь и постоянные.

Из соображений более точной аппроксимации (12) было выбрано . В этом случае выражение (14) будет почти точно совпадать с хольцмаркским сдвигом (12), хотя внешне оно от него несколько отличается. Однако упомянутый ранее микрополевой сдвиг Волокитина-Голосного-Калиткина из [7] фактически не успевает выйти на хольцмаркскую ассимптотику даже для газов невысоких плотностей; оказывается примерно в три - четыре меньше. Поэтому и в данном выражении для сдвига потенциала ионизации (14) было принято .

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

, (15)

где - статистический вес электронного уровня, - энергия электронного уровня, отсчитанная от основного состояния.

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

, (16)
Простейшие оценки формфактора основаны на следующем соображении. Пусть потенциал ионизации k-кратного иона равен , а энергия некоторого уровня есть (она отсчитана от основного состояния). Тогда для ионизации с данного возбужденного уровня достаточно энергии . Если эта величина много меньше, чем средний сдвиг потенциала , то нахождение электрона на данном уровне маловероятно: . В противном случае уровень практически заселен: . Таким образом этот эффект имеет вероятностный характер, и вероятность заселения должна зависеть от отношения . На основании всего выше сказанного в качестве формфактора нами было выбрано следующее выражение:

, (17)

Поскольку микрополевой сдвиг (14) , то , как и должно быть.

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

Систему уравнений ионизационного равновесия (5) можно, пользуясь балансными соотношениями (1), свести к одному уравнению:

, (18)

где , а x – относительная концентрация электронов. Тогда относительные концентрации ионов вычисляются по формуле:

. (19)

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

В третьей главе описывается программный комплекс, который реализует описанный во второй главе численный алгоритм для решения системы уравнений ионизационного равновесия (5) с предложенным в работе микрополевым сдвигом потенциала ионизации (14). Он реализован на языке C++ и состоит из непосредственно программного ядра и прилагающейся к нему таблице потенциалов ионизации элементов, взятой из [8]. В этой таблице содержатся потенциалы ионизации всех элементов периодической системы от (водород) до (лоуренсий), причем кратностей от до z. Потенциалы низших ионизаций элементов взяты из экспериментальных данных, собранных в [9]. Для более высоких кратностей ионизации потенциалы построены в [8] на основе расчетов по релятивисткой модели Хартри-Фока-Слэтера с учетом специальных интерполяционных поправок.

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

При расчете таблицы для выбранных диапазонов T и V программа организована таким образом, что на каждом шаге в качестве начального приближения и принимаются значения в предыдущей точке расчетов. Это уменьшает число итераций и экономит время. Поэтому рассчитать всю таблицу данных сразу получается быстрее, нежели по отдельности в каждой точке. Пока скорость работы данного программного комплекса не очень высока. Она сильно зависит от атомного номера элемента. Например, время, затраченное на расчет приведенных в диссертационной работе таблиц (каждая содержит 273 точки) для Cu составляет около 1 часа, в то время как аналогичные таблицы для Al вычисляются за пару минут. Причин здесь несколько. Во-первых, сам алгоритм является достаточно трудоемким (порядка операций сложения, умножения и возведения в степень), во-вторых, программный комплекс написан на языке С++, который сам по себе не приспособлен для таких расчетов. И, наконец, пока еще мы не ставили перед собой задачу оптимизации вычислительного процесса.

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

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

Для демонстрации работоспособности программного комплекса были проведены расчеты плазмы лития, алюминия и меди. Расчеты проводились в диапазоне температур (тыс.К) и объемов с шагом и соответственно (T и V в атомных единицах). Для каждого элемента рассчитаны три таблицы: электронная концентрация x, давление P и энергия E. В таблицах приведены значения десятичных логарифмов указанных величин. Таблицы представлены в диссертационной работе.

Проведены также сравнительные расчеты ионизации на примере лития с микрополевой поправкой из [7] и предложенным в данной работе сдвигом потенциалов ионизации. Результаты для () и (тыс.К) представлены на рис.11. Видно, что некоторое количественное различие наблюдается лишь в области первой ионизации, в то время как качественно модели идентичны. Причем количественная разница уменьшается с увеличением объема.
Рис.11 Расчет ионизации для лития по микрополевой модели [7] и модели, предложенной в данной работе.

На рисунке представлено изменение относительной электронной концентрации х при V=10000 ()



В заключении сформулированы основные результаты работы.

В приложении приводится таблица потенциалов ионизации, используемая в программном комплексе и текст программы.

Основные результаты



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

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

3. С использованием программного комплекса проведены численные расчеты, которые показали:

- модель разумно воспроизводит явление металлизации;

- аккуратный учет возбужденных уровней в статистических суммах слабо влияет на ионизацию и термодинамику; доминирующий вклад вносит основное состояние и его обрезание;

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

Литература



1. Калиткин Н.Н., Царева Л.С. Метод расчета ионизации на ЭВМ.// ЖВММФ, 1971, т.11, №3, С.772-773.

2. Баско М.М. Уравнение состояния металлов в приближении среднего иона. М.: Препр. ИТЭФ, 1982, №57, 44с.

3. Калиткин Н.Н., Ритус И.В., Миронов А.М. Ионизационное равновесие с учетом вырождения электронов. М.: Препр. ИПМ, 1983, №46, 27с.

4. Тиман Л.Б. Влияние взаимодействия ионов на их равновесные концентрации в случае многократной термической ионизации газа.// ЖЭТФ, 1954, т.27, вып. 6(12), с.708-711.

5. Ликальтер А.А. Взаимодействие атомов с электронами и ионами в плазме. // ЖЭТФ, 1969, т.56, №1, с.240-245.

6. Эбелинг В., Крефт В., Кремп Д. Теория связанных состояний и ионизационного равновесия в плазме и твердом теле. М.: Мир, 1979, 262с.

7. Волокитин В.С., Голосной И.О., Калиткин Н.Н. Широкодиапазонное уравнение состояния вещества. 1.Анализ моделей неидеальности. // Известия вузов; физика, 1994, №11, с.23-43.

Волокитин В.С., Голосной И.О., Калиткин Н.Н. Широкодиапазонное уравнение состояния вещества. 2.Микрополевая модель. // Известия вузов; физика, 1995, №4, с.11-31.

8. Калиткин Н.Н., Кузьмина Л.В., Ритус И.В. Потенциалы ионизации и ударные адиабаты плазмы. М.: Препр. ИПМ, 1986, №120, 24с.

9. Физические величины: Справочник/ А.П. Бабичев, Н.А. Бабушкина, А.М. Братковский и др. Под ред. И.С. Григорьева, Е.З. Мейлихова – М.; Энергоатомиздат, 1991.-1232с.

Основные результаты опубликованы в работах



Калиткин Н.Н., Павлов А.С. Метод расчета состава неидеальной плазмы. Математическое моделирование, 2004, т.16, №12, с.61-68.

Калиткин Н.Н., Павлов А.С. Микрополевая неидеальность и металлизация плазмы. Математическое моделирование, 2005, т.17, №6, с.21-32.

Калиткин Н.Н., Павлов А.С. Программный комплекс «Плазма-элемент». Математическое моделирование, 2005, т.17, №10, с.14-30.

Калиткин Н.Н., Павлов А.С. Микрополевая неидеальность и металлизация плазмы. Экстремальные состояния вещества, детонация, ударные волны. Труды международной конференции VII харитоновские тематические научные чтения, Саров, 2005, под ред. А.Л.Михайлова, с.312-317.

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

кандидата физико-математических наук

Отпечатано в Институте математического моделирования РАН

Тираж 100 экз.


Добавить документ в свой блог или на сайт

Похожие:

Самосогласованная микрополевая модель неидеальной плазмы icon Модели данных, использующих в гис, можно разделить на гибридную модель,...
Енной и атрибутивной информации. Пространственная (координаты) и топологическая информация хранится в собственных файлах гис-приложения....
Самосогласованная микрополевая модель неидеальной плазмы icon Литература Александров А. Ф., Рухадзе А. А. Основы электродинамики...
Термодинамический подход к описанию молекулярных явлений. Термодинамическое равновесие. Температура
Самосогласованная микрополевая модель неидеальной плазмы icon 2 Особенности эконометрических методов прогнозирования
Таким образом, в самом общем виде эконометрическая модель – это дескриптивная модель исследуемого объекта, в структуре которой можно...
Самосогласованная микрополевая модель неидеальной плазмы icon План 1 июня 2011 г. 1: 35 Структура приложения(mvc) База данных Использованные технологии
«Модель-представление-контроллер» — шаблон проектирования, в котором модель данных приложения,пользовательский интерфейс и управляющая...
Самосогласованная микрополевая модель неидеальной плазмы icon Воздействие холодной плазмы на структуру суставного гиалинового хряща
Кроме этого принимается во внимание необходимость точного восстановления анатомической формы суставных поверхностей, обеспечивающего...
Самосогласованная микрополевая модель неидеальной плазмы icon Интегрированная модель аттестации педагогических работников на квалификационную категорию
Ю комплексного установления теоретической подготовленности педагога и результативности его практической деятельности разработана...
Самосогласованная микрополевая модель неидеальной плазмы icon Программа Государственного экзамена по подготовке магистра по направлению...
Влияние объемного заряда на движение заряженных частиц. Самосогласованная система уравнений электронного пучка
Самосогласованная микрополевая модель неидеальной плазмы icon Модель «Портфолио +» как одно из звеньев мониторинга качества воспитания
Наименование проекта Модель «Портфолио +» как одно из звеньев мониторинга качества воспитания
Самосогласованная микрополевая модель неидеальной плазмы icon Модель профессионально значимых качеств педагога дополнительного...
Модель профессионально значимых качеств педагога дополнительного образования можно предоставить
Самосогласованная микрополевая модель неидеальной плазмы icon Проекта
Модель многоуровневого многопрофильного профессионального образовательного комплекса 66
Литература


При копировании материала укажите ссылку © 2015
контакты
literature-edu.ru
Поиск на сайте

Главная страница  Литература  Доклады  Рефераты  Курсовая работа  Лекции