Добавил:
okley
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Диплом / 111
.pdf
Министерство образования è íàóêè
Ðоссийскîé Федерации
Федеральное агенство ïî образованию
ÌÎÑÊОВСКИЙ ИНЖЕНЕРНО-ФИЗИЧЕСКИЙ ИНСТИТУТ
(ÃÎÑÓДАРСТВЕННЫЙ УНИВЕРСИТЕТ)
Ôàêóëüòåò ýêспериментальной è теоретическîé физики
Кафедра теоретическîé ядерной физики 32
Диплом íà соискàíèå степени áàêалавра
Áåæенцев Борис Владимирович
Ðасчет ионизационного состоÿíèÿ
плазмы аргона
íà основе параметрическîé системы
кинетических уравнений
Москва 2007
Íàучный ðóêîâîдитель:
ä.ô.-ì.í., профессор
Â.Ã. Новикîâ
Аннотàöèÿ
Íà основе аргоновой плазмы проведено óсреднение äåòальных
уровней энергии â эффективные группы ñ íåêоторой средней
энергией è ñòатистическим весом. Íà основании элементарных
процессов, протекающих â плазме, áûëè написаны è численно
решены уравнения поуровневой кинетики. Проведен расчет
степени ионизации вещества â случае прозра÷íîé è оптически
плотной плазмы. Òàêæå проведен расчет âûõîäà излучения èç
плазмы â приближåíèè òîíêîãî ñëîÿ. Ðасчеты проводились äëÿ
диапазона температур T = 1ö1000 ýÂ è плотностей свобоäíûõ
электронов
Ne= 1010ö
10
24
1/ñì3.
Автор благодарит çà íàучное ðóêîâîдство Новикîâà Владимира
Григорьевича , à òàêæå çà помощь â работе ñòаршего íàучного
сотруäíèêà èç îòäåëà 6 ÈÏÌ èì. Ì.Â. Келдыша À .Ä .
Соломянную.
Ñîäåðæàíèå
Введение 2
1 Ïîñòановкà çàäà÷è 3
2 Основные процессы â плазме 5
2.1 Столкновительные è радиационные процессы . . . . . . . . 5
2.2 Î ñêоростÿõ столкновительных è радиационных процессах 6
2.2.1 Ó÷åò òîíêîé структуры уровней . . . . . . . . . . . . 8
3 Система уравнений поуровневой кинетики 10
3.1 Система уравнений äëÿ определения êонцентраций ионов . 10
3.2 Квазистационарный случай . . . . . . . . . . . . . . . . . . . 11
4 Ðåçóëüòàòû расчетов 12
4.1 Ðåçóëüòàò óсреднения äåòальных уровней энергии . . . . . 12
4.2 Ионизационное равновесие . . . . . . . . . . . . . . . . . . . 12
4.3 Излучение оптически òîíêîé плазмы . . . . . . . . . . . . . 22
Заключение 25
Список литературы 26
1
Введение
Исследование аргоновой плазмы представляет большой интерес,
благоäàðÿ åå применению â современных отраслях техники.
Àêòóальность òåìû связана ñ òåì, ÷òî термодинамические è оптические
свойства неравновесной плазмы плазмы, íåîáõîäèìû ïðè ìîделировании
физических процессов, êоторые используютсÿ äëÿ создания источникîâ
излучения, ïðè плазменной обработкå материалов, äëÿ исследования
è оптимизации термоядерных мишеней. Â изучаемых современной
физикîé процессах реализуютсÿ самые разнообразные óсловия: îò
температур ïîðÿäêà абсолютного íóëÿ äî температур â недрах звезд,
îò ãазовых плотностей äî плотностей, â äåñÿòêè ðàç превышающих
плотность òâåðäîãî òåëà. Â связи ñ ýòèì решение многих àêòóальных
проблем физики требуåò âñå большего объема информации î свойствах
вещества â различных, â òîì числе ýêстремальных, óсловиях.
Äëÿ описания свойств вещества â äîñòаточном широкîì диапазоне
температур è плотностей, êàê правило, неприменимы сравнительно
простые ìîäåëè. Ñ другой стороны, проведение ýêспериментов ñ целью
получения данных î свойствах вещества â ýêстремальных óсловиях
обычно связано ñî çíàчительными техническими òðóдностÿìè è â ðÿäå
случаев îêазываютсÿ чрезвычайно дорогостоящими. Поэтому требуåòñÿ
создание быстрой è эффективной ìåòîäèêè расчетà свойств плазмы äëÿ
использования â ãàçîдинамических программах.
 данной работе произвоäèòñÿ расчет ионизационного состоÿíèÿ
плазмы аргона. Â íåé рассматриваютсÿ основные элементарные
процессы, происõîдящие â ãîðÿ÷åé плазме, òàêèå, например, êàê
ионизация электронным óдаром, фотоионизация, возбуждение, ãашение
è другие. Íà основании ýòèõ процессов пишутсÿ уравнения поуровневой
кинетики, à затем, после èõ решения, получаютсÿ êонцентрации ионов
â различных состояниях â зависимости îò температуры è плотности
электронов. Òàê êàê энергетический спектр многоэлектронного атома
состоит èç ìíîæентва уровней, òî система кинетических уравнений
îêазалась áû громадной. Поэтому провоäèòñÿ óсреднения реальной
ñõåìû уровней. Ñ помощью ýòèõ êонцентраций рассчитываетсÿ средняя
степень ионизации вещества.
2
1 Ïîñòановкà çàäà÷è
Ìû áóäåì изучать õарактеристики плазмы â äâóõ предельных
случаях: â приближåíèè оптически òîíêîé плазмы ( прозра÷íàÿ плазма )
è â приближåíèè оптически плотной плазмы. Â общем случае äëÿ
решения çàäà÷è возникàåò íåîáõîдимость совместного рассмотрения
процессов переноса фотонов è поуровневой кинетики электронов [1].
 квазистационарном приближåíèè è áåç ó÷åòà процессов рассеяния
уравнение переноса фотонов имеет âèä
(−→Ω ∇)Iω= (ω, Iω) − κ(ω, Iω)I
ãäå
Iω интенсивность излучения фотонов ñ энергией~ω,
êоэффициент ïîãлощения â 1/ñì,
способность,−→Ω
единичный вектор â направлении распространения
= (ω, Iω)
объемная излучательная
ω
κ = κ(ω, Iω)
(1)
излучения.
Энергетический спектр реального многоэлектронного атома состоит
èç ìíîæества уровней, õарактеризуåìûõ соответствующим набором
квантовых чисел. Система кинетических уравнений, записанных äëÿ
êаждого èç ýòèõ уровней, âûãлядела áû совершенно необозримой.
Поэтому ïðè решении кинетических çàäà÷ произвоäèòñÿ íåêоторое
¾огрубление¿ реальной ñõåìû уровней. Â первую очередь îíî связано
ñ объединением групп близкî располоæенных уровней â íåêоторые
эффективные уровни. Следуåò отметить, ÷òî, êонечно, íå существуåò
универсального способа òàêîãî объединения è, более òîãî, способ
пригоäíûé ïðè решении îäíèõ çàäà÷, ìîæåò áûòü непригоäíûì äëÿ
других. Ïðè получении огрубленных энергетических ñõåì различных
атомов è ионов использованы следующие ðåêомендации.
1. Â îäèí уровень объединяютсÿ группы энергетически близких
состоÿíèé. Этому уровню приписываютсÿ эффективный ñòатистический
âåñ è эффективная энергия, вычисляемые ïî формуëàì
gks=
X
γJ ∈ks
g
, Eks=
ksγJ
X
1
g
ks
γJ ∈ks
g
ksγJEksγJ
(2)
ãäå k çàðÿä èîíà, s номер объединенного уровня, J полный момент,
γ
дополнительный набор квантовых чисел, êоторые определяют
состоÿíèå ñ энергией
E
è ñòатистический âåñ
ksγJ
g
ksγJ
= 2J + 1
.
Коэффициенты ïîãлощения è излучательная способность плазмы
определяютсÿ населенностÿìè основных è возбужденных состоÿíèé
3
ионов, òî åñòü êоличеством ионовk-îé кратности, íàõîдящихñÿ â
íåêотором состоÿíèès. Обозна÷èì ýòè êонцентрации через
nks= nks(−→r ).
Òàêæå áóäåì использовать относительные êонцентрации:
xks= nks/n
i
ãäå плотность ионов
ni= ρNA/A,
A
атомный âåñ (äëÿ аргона A = 39,948 ã/ìîëü. ),
NA= 6.022·10231/
ìîëü
постоянная Àâîãàäðî.
Äëÿ êонцентрации частиц
ñ êîîðдинатой−→r
â состоÿíèè
nks, íàõîдящихñÿ â малом îáú¼ìå
k s
, выполняетсÿ óсловие балланса
(уравнение ñòатистическîãî равновесия ) [2]:
∂n
ks
= −∇ · (nks~v) +
∂t
X
k0s
0
(n
w
k0s
k0s0→ks
0
− nksw
ks→k0s
0
),
òî åñòü изменение å¼ ñî временем áóäåò определятьсÿ потокîì частиц
через ограничивающую ýòîò îáú¼ì поверхность è темпом поступления
èç других состоÿíèé
процессов. Здесь
âåðîятность перехîäà èîíà èç состоÿíèÿksâ состоÿíèå
k0s0çà ñ÷åò радиационных è столкновительных
~v
макроскопическàÿ ñêорость,
w
ks→k0s
0s0
k
0
.
полная
Благоäàðÿ òîìó, ÷òî â нашем случае ñêорости ãàçîдинамических
процессов существенно меньше ñêоростей элементарных процессов,
влиянием макропереноса ìîæíî пренебречь. Òîãäà уравнение примет
âèä
∂n
∂t
ks
=
X
k0s
0
(n
w
k0s
k0s0→ks
0
− nksw
ks→k0s
0
).
Âåðîятности столкновительных процессов определяютсÿ плотностью
свобоäíûõ электронов è èõ распределением ïî энергиям. Âåðîятности
радиационных процессов зависÿò îò функции распределения
фотонов, точнее îò òîãî, êàêàÿ часть излучения îñòà¼òñÿ â веществе
(реабсорбируåòñÿ). Åñëè заметная часть излучения реабсорбируåòñÿ,
òî уравнения 1 è 2 íåîáõîäèìî решать совместно. Â ýòîì случае
êонцентрации различных ионов сильно зависÿò îò интенсивности
4
излучения, à îíà, â ñâîþ очередь, зависит îò êонцентрации ионов через
êоэффициент ïîãлощения è излучательную способность плазмы [1].
 точной ïîñòановкå çàäà÷è â формуëû äëÿ âåðîятностей
элементарных процессов âõîäèò спектральная плотность излучения
W (ω/T )
, êоторая определяетсÿ решением уравнения переноса â заданной
геометрии.  работе рассмотрен предельный случай прозра÷íîé
плазмы, òîãäà
W (ω/T ) = 0
, è случай êîãäà выполняетсÿ ëîêальное
термодинамическîå равновесие. Òîãäà ïîëå излучения определяетсÿ
формуëîé Планкà ñ температурой ïîëÿ равной электронной температуре
W (x) =
1
, x =
ex− 1
ω
θ
2 Основные процессы â плазме
2.1 Столкновительные è радиационные процессы
Áóäåì рассматривать плазму, состоÿùóþ èç электроновeè ионов
Nksразличного çàðÿäàk, íàõîдящихñÿ â íåêотором состоÿíèès, è
фотонов~ω. Â плазме èäóò процессы взаимодействия различных частиц
äðóã ñ другом è ñ излучением. Соответственно процессы называютсÿ
столкновительными è радиационными.
 работе рассматривались неупругие столкновения:
1. Возбуждение электронным óдаром è обратный процесс ãашения:
k
N
+ e À N
s
2. Ионизация электронным óдаром è тройная ðåêомбинация:
k
N
+ e À e + e0+ N
s
Òàêæå рассматривались процессы взаимодействия ионов ñ фотономи:
1. Ïîãлощение è излучение â спектральных линиях:
k
N
+ ~ω À N
s
2. Фотоионизация è фоторекомбинация:
k
N
+ ~ω À N
s
k
0
s
k+1
s
0
+ e
k
s
0
s
0
+ e
k+1
0
Ðадиационные процессы делятсÿ íà спонтанные è вынужденные.
Первые зависÿò толькî îò свойств излучающей системы. Вторые åùå
è îò плотности излучения â интервале частот, соответствующем
рассматриваемому перехîäó. Èõ çíàчение возрастàåò ïî ìåðå
уменьшения плотности плазмы è ñíèæåíèÿ частоты столкновительных
процессов.
5
2.2 Î ñêоростÿõ столкновительных è радиационных
процессах
Основную ðîëü â столкновительных è радиационных процессах
â плазме играют электроны, à возбуждение òÿæелыми частицами
(ионами ) существенно толькî äëÿ перехîäîâ между очень близкими
уровнями ( ìû пренебрегàåì ýòèì процессом). Â данной работе
учитываютсÿ íå âñå элементарные процессы. Â частности íå
рассматриваетсÿ êомптон-эффект è тормозные процессы, òàê
êàê ïðè протекàíèè ýòèõ реакций íå изменяютсÿ интересующие
íàñ êонцентрации ионов. Íåò автоионизационных процессов è
диэлектрическîé ðåêомбинации. Â ðàìêàõ рассматриваемой ìîäåëè
ñ óсреднением ýòîò вопрос íå äîñòаточно проработàí. Â частности,
использование ñòандартного ïîäõîäà â нашем приближåíèè привоäèò
ñ ðåçêîìó ñíèæåíèþ степени ионизации, ÷òî связано ñ преувеличением
ðîëè ýòèõ процессов â ðàìêàõ параметрическîé ìîäåëè. Поэтому â
основных расчетàõ диэлектронная ðåêомбинация è автоионизация íå
учитывались.
Обозна÷èì çà
 случае реакций ñ фотонами ýòî число спонтанных èëè вынужденных
перехîäîâ â единицу времени, ïðèõîдящихñÿ íà îäèí àòîì. Â
случае столкновительных процессов ýòî число актов возбуждения â
единицу времени íà îäèí àòîì, êоторое пропорционально êонцентрации
электронов
Ne[3]:
ωifâåðîятность перехîäà èç состоÿíèÿiâ состоÿíèåf.
ωif= Ne< vσif>,
ãäå
σif эффективное сечение процесса возбуждения â ñì2,v ñêорость
налетающих электронов.
Çíàÿ функцию распределения электронов ïî энергиям
писать
ãäå
4ε
< vσif>=
энергия возбуждения.
∞
Z
vσif(ε)F (ε) dε ,
4ε
F (ε)
, ìîæåì
Заметим, ÷òî â случае термодинамическîãî равновесия
распределение свобоäíûõ электронов ïî энергиям
ìàêсвелловским:
F (ε) =
√
2
πT
√
−ε/T
εe
3
,
F (ε)
являетсÿ
ãäåT температура â ýÂ.
6
Âåðîятность процесса фотоионизации равна:
∞
Z
ãäå
Ei энергия ионизации,
wiν=
σiν сечение фотоионизации,
Nωcσiνdω,
E
i
Nω плотность
квантов ñ энергиейω.
Åñëè â íàчальном состоÿíèè имеютсÿ òðè частицы, òî выражåíèÿ
äëÿ âåðîятностей âûãлядят èíà÷å. Например, äëÿ 3-¼õ частичной
ðåêомбинации имеем
Z Z
wir= N
2
e
v1v2σrF (ε1)F (ε2) dε1dε2.
Здесь
σrимеет размерность ñì4·
c, íî òàêæå называетсÿ сечением.
Òàêèì образом, получаем зависимости âåðîятностей
рассматриваемых процессов îò êонцентрации частиц (фотонов èëè
электронов):
ионизация:
3-¼õ частичная ðåêомбинация:
фотоионизация:
фоторекомбинация:
ïîãлощение:
излучение:
возбуждение:
ãашение:
wii∼ N
wir∼ N
phi
w
phr
w
abs
w
e
2
e
∼ W (ω/T )
∼ Ne(1 + W (ω/T ))
∼ W (ω/T )
wem∼ (1 + W (ω/T ))
wex∼ N
dex
w
∼ N
e
e
Выпишем приближ¼ííûå формуëû äëÿ âåðîятностей основных
элементарных процессов (èõ è áóäåì использовать â дальнейшем )
[1, 4, 5].
1. Столкновительная ионизация электрона ñ уровняnläëÿ состоÿíèÿ
ks
â приближåíèè Лотца [6]:
ãäå
wii(ks → k + 1, s) = 1.28 · 1016N
N
1
N
e
, E1(x) =
1/2
θ
A
C0=
7
C
nl
0
E1(ω/θ)
ω
∞
Z
e
x
(3)
−t
dt
t
θ = 36.75 · T · 10
ω = E
k+1,s
0
− Eks энергия перехîäà.
−3
температура â атомных единицах,
2. Тр¼хчастичная ðåêомбинация íà уровеньnlâ состоÿíèåks:
wir(k + 1, s0→ ks) = C
ks→k+1, s
0
wii(ks → k + 1, s0)
ãäå
g
N
C
ks→k+1, s
0
= 0.704
g
k+1, s
ks
0
N
1
e
e
3/2
θ
A
3. Фотоионизация ñ уровняnläëÿ состоÿíèÿks:
∞
Z
phi
w
(ks → k + 1, s0) = Nnlu
ks→k+1, s
0
ω/θ
ãäå â приближåíèè Крамерса
3/2
u
ks→k+1, s
l
- орбитальный момент.
0
= 4.45 · 1010Z
|ω|
ef
2(2l + 1)
, Zef= k + 1,
4. Фоторекомбинация íà уровеньnlâ состоÿíèåks:
∞
Z
−t
phr
w
(k + 1, s0→ ks) = C
ks→k+1, s
0
u
ks→k+1, s
0
ω/θ
e
t
ω/θ
W (t)
dt
t
[1 + W (t)] dt
(4)
(5)
(6)
2.2.1 Ó÷åò òîíêîé структуры уровней
Äëÿ ó÷åòà òîíêîé структуры уровней ñêорости перехîäîâ между
óсредненными состояниями ñ îäíîé степенью ионизации имеем формуëû
(учитываем, ÷òî
X
ks→ks
0
=
γJ ∈ks
w
X
w
ks0→ks
=
γ0J0∈ks
E
ksγJ
X
γ0J0∈ks
X
0
γJ ∈ks
0
< E
g
ksγJ
g
g
ks0γ0J
ks
g
ks
ks0γ0J
ex
[w
ks→ks
0
[w
0
0
)
(γJ → γ0J0) + w
0
dex
(γ0J0→ γJ) + w
ks0→ks
abs
(γJ → γ0J0)],
0
ks→ks
em
ks0→ks
(7)
(γ0J0→ γJ)].
(8)
Здесь учтено, ÷òî состоÿíèå
посредством îдноэлектронного перехîäà
ks0получаетсÿ èç состоÿíèÿ
nl → n0l0.
ks
8