Жданов Явления переноса в газах и плазме 2008
.pdfj |
|
= σ |
|
E |
|
, |
σ |
|
|
|
= σ0 , |
(7.57) |
|
|
|
|
|
где E = k ( E k). Таким образом, плотность тока вдоль на-
правления магнитного поля определяется тем же значением электропроводности σ , что и в отсутствие магнитного поля.
Полная плотность тока определяется как j = j + j , где |
|
|
j = k (j k) , |
j = j − j = k ×(j×k) . |
(7.58) |
Решение уравнения (7.54) для составляющих плотности тока, перпендикулярных магнитному полю, можно найти в результате следующих простых операций. Образуем векторное произведение
уравнения (7.54) с вектором k .В результате имеем
σ(E ×k)= (1+ s)(j×k)−β0 j .
Выражая отсюда (j×k)
(7.57), находим
σ(E +E )= j
и подставляя его в исходное уравнение
+ j + |
β02 |
j + |
σβ0 |
(E ×k)+ s j . |
|
1+ s |
|||||
|
|
1+ s |
|
Если исключить из этого уравнения соотношение для продольных
величин (7.57), то для j получаем |
|
|
|
|
|
|||
|
σ(1+ s) |
|
σβ |
0 |
(k ×E |
|
) , |
(7.59) |
j = |
|
E + |
|
|
||||
(1+ s)2 +β02 |
(1+ s)2 |
+β02 |
|
где E = [k ×(E ×k)].
Окончательное выражение для плотности тока проводимости в плазме j = j + j можно представить тогда в виде
j = σ |
|
|
|
E |
|
|
|
+σ E +σH (k ×E ) . |
(7.60) |
|
|
|
|||||||
|
|
201
Здесь σ|| , σ и σH – продольная, поперечная и холловская элек-
тропроводности плазмы, определяемые соответственно выражениями
σ |
|
= σ , |
σ = σ |
|
1 + s |
|
, |
σH |
= σ |
|
β0 |
|
|
|
(7.61) |
||||||||
|| |
(1 + s)2 +β02 |
(1 + s)2 |
+ β02 |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
В полностью ионизованной плазме |
δ0 = 0 |
|
и τ0 = τei . В этом |
||||||||||||||||||||
случае |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
n |
e |
2 |
|
|
|
|
|
1 |
|
|
|
|
|
|
ω |
τ |
|
|
|
|
σ|| = σ = |
e |
|
|
τei |
, |
σ = σ |
|
|
|
|
|
, σH |
= σ |
|
e |
|
ei |
. (7.62) |
|||||
me |
|
|
+ ωe2 |
τei2 |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
1 |
|
|
1+ ωe2 τei2 |
|
Из полученных выше выражений следует, что магнитное поле не влияет на электропроводность плазмы вдоль силовых линий магнитного поля. Если электронный параметр Холла мал (β0 <<1 ), то
магнитное поле почти не влияет на токи поперек магнитного поля.
Однако при βo ≥1 это влияние оказывается заметным. При этом магнитное поле может заметно ослаблять ток в направлении попе-
речной компоненты электрического поля E , а также создавать ток Холла в направлении, перпендикулярном как электричеcкому, так и магнитному полю. Заметим, что для частично ионизованной
плазмы, удовлетворяющей условию τen−1 ≥ τei−1 , имеем |
|
δ = (1−α)2 ββ0 ~ ε <<1 , |
(7.63) |
i |
|
поэтому влияние параметра проскальзывания ионов s = δ0β02 на-
чинает сказываться на значениях σ и σH гораздо позднее, чем
влияние параметра Холла для электронов, а для полностью ионизованной плазмы эффект проскальзывания вообще отсутствует. В
последнем случае при ωe τei >>1 получаем
σ = |
σ |
|
, |
σH = |
σ |
|
. |
(7.64) |
||
ω2 |
τ2 |
|
|
|||||||
ω |
τ |
|
||||||||
|
|
|
ei |
|
||||||
|
e |
|
ei |
|
|
e |
|
|
202
7.4. Диффузия в слабоионизованном газе
Рассмотрим слабо ионизованный газ, образованный из электронов, ионов и нейтральных атомов одного сорта. Напомним, что условие слабой ионизации определяется как
νei << νen . |
(7.65) |
Поскольку при этом α <<1, по отношению к независимой диффузии заряженных частиц нейтральный компонент выступает как
среда, скорость которой u n фактичеcки совпадает со среднемас-
совой или среднемолярной скоростью смеси. Молярные диффузионные потоки заряженных частиц определены при этом как
Jm = n |
α |
(u |
α |
− u |
n |
), |
α = e,i . |
(7.66) |
α |
|
|
|
|
|
Рассмотрим вначале случай, когда магнитное поле отсутствует, а электрическое поле действует в направлении оси X . Будем полагать, что в том же направлении изменяются плотности заряженных компонентов. Тогда уравнения (7.17) для электронного и ионного компонента запишутся в виде
n m |
|
ν |
|
|
(u |
|
|
−u |
|
|
)= −n eE − |
dpe |
|
, |
|||||||||||||||||
|
|
|
|
|
|
|
|
dx |
|||||||||||||||||||||||
e |
e |
|
|
en |
|
|
|
e |
|
|
|
|
n |
|
|
|
e |
|
|
|
|||||||||||
n |
μ |
|
|
ν |
|
|
|
(u |
|
−u |
|
|
)= n ZeE − |
dpi |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
dx . |
(7.67) |
||||||||||||||||||||||
i |
|
in |
|
|
in |
|
|
i |
|
|
|
|
|
n |
|
|
i |
|
|
|
|||||||||||
Для молярных диффузионных потоков отсюда получаем |
|||||||||||||||||||||||||||||||
|
|
|
|
J m = −D |
|
dne |
−μ n E |
|
, |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
e |
|
|
|
|
|
|
|
e |
|
|
dx |
e e |
|
|
|
|
|
|
|
|
|
||||
|
|
|
J m |
= −D |
|
dni |
+ Zμ n E |
|
, |
|
|
|
|
(7.68) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
i |
|
|
|
|
|
|
|
i |
|
dx |
i i |
|
|
|
|
|
|
|
|
|
|||||||
где использованы соотношения |
|
|
|
pe = nekT , |
|
|
pi = ni kT . |
||||||||||||||||||||||||
Коэффициенты диффузии электронов и ионов определены как |
|||||||||||||||||||||||||||||||
|
De = |
|
kTe |
|
|
|
, |
|
|
|
Di = |
|
kTi |
|
. |
(7.69) |
|||||||||||||||
|
m ν |
en |
|
|
|
|
|
|
μ |
in |
ν |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
e |
|
|
|
|
|
|
|
|
|
|
|
|
|
in |
|
203
При этом коэффициенты подвижности μα |
и диффузии Dα связа- |
|||||
ны соотношением Эйнштейна |
e |
|
|
|||
μ |
α |
= |
D |
(7.70) |
||
|
||||||
|
|
|
α . |
|||
|
|
|
kTα |
|
||
Заметим, что в соответствии с оценкой (7.51) |
||||||
De >> Di |
, |
μe >>μi , |
(7.71) |
поэтому скорость диффузии электронов намного превышает скорость диффузии ионов.
Полученные выше выражения легко обобщить и на случай диффузии заряженных частиц поперек магнитного поля. Уравнения (7.17) для электронного и ионного компонента в рассматриваемом нами случае слабо ионизованной плазмы принимают вид
Jem +βe (Jem ×k)= −De ne −μeneE , |
|
Jim +βi (Jim ×k)= −Di ni + Ziμi niE , |
(7.72) |
где βe = ωe τen и βi = ωi τin . Рассмотрим в качестве примера ци-
линдрический столб плазмы, помещенный в продольное (направленное по оси цилиндра) магнитное поле. Предположим, что электрическое поле и все градиенты направлены лишь в радиальном направлении. Тогда из решения уравнений (7.72) для радиальных потоков поперек магнитного поля получаем
Jem = −De dndre −μe neEr ,
|
J m |
= −D |
|
dni |
+ Zμ |
i |
n E |
r |
, |
(7.73) |
|||
|
|
|
|||||||||||
|
i |
|
i dx |
|
i |
|
|
||||||
где |
|
|
|
|
|
|
|
|
|
|
|
||
Dα = |
kTα |
bα = |
|
Dα |
|
, |
|
(α = e,i) . |
(7.74) |
||||
|
1 +βα2 |
|
|||||||||||
|
e |
|
|
|
|
|
|
|
Заметим, что в сильном магнитном поле, когда βα2 |
>>1, имеют ме- |
|
сто неравенства |
μe << μi , |
|
De << Di , |
(7.75) |
т.е. условия, обратные условиям (7.71). Это означает, что скорость радиальной диффузии ионов поперек магнитного поля может зна-
204
чительно превышать скорость радиальной диффузии электронов. Мы оставляем здесь в стороне анализ холловских диффузионных потоков, перпендикулярных как электрическому, так и магнитному полям.
Рассмотрим теперь вопрос о так называемой амбиполярной диффузии в слабоионизованной плазме. Пусть магнитное поле отсутствует. При диффузии заряженных частиц в поле градиентов их собственного парциального давления из-за высокой подвижности электронов по сравнению с ионами может возникать разделение зарядов. Возникающее при этом электрическое поле (поле объемного заряда) начинает препятствовать диффузии электронов и, наоборот, способствовать диффузии ионов. Это поле должно расти до тех пор, пока не станут одинаковыми диффузионные потоки электронов и ионов, что соответствует установлению квазистационарного состояния. Совместная диффузия электронов и ионов носит название амбиполярной диффузии.
В качестве конкретного примера рассмотрим радиальную диффузию заряженных частиц в длинной цилиндрической трубе в условиях, когда на ее стенках происходит интенсивная рекомбинация частиц, благодаря чему возникают радиальные градиенты плотности заряженных частиц. Пусть в начальный момент времени квазинейтральность плазмы точно осуществляется во всем объеме. В отсутствие внешнего электрического поля из-за большого различия
вкоэффициентах диффузии электронов и ионов (см. условие (7.71)) вблизи стенок будет нарастать объемный отрицательный заряд и возникнет радиальное электрическое поле объемного заряда. В соответствии с соотношениями (7.68) это поле, так же, как и градиенты плотности заряженных частиц, будет воздействовать на их диффузию. Квазистационарное состояние наступает, когда диффузионные потоки сравниваются по своей величине, что соответствует отсутствию электрического тока в плазме. (Последнее возможно, например, при отсутствии контакта плазмы с внешними проводниками). Действительно, из определения тока проводимости
вплазме при выполнении условия ne = Zni получаем
j = −e(Jem − Jim ). |
(7.76) |
205
Условие j = 0 можно использовать для того, чтобы найти величину напряженности электрического поля, возникающего в плазме в квазистационарном состоянии. Приравнивая выражения для J em и
J im (7.68) и полагая |
dne |
dr = Zdni |
dr , находим |
|
|
||||||||||||
|
|
ne Er = |
Di − ZDe |
|
|
dni |
|
|
|
|
(7.77) |
||||||
|
|
μe +μi |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
dr |
|
|
|
|
|
|||
(Мы заменили здесь координату x радиальной координатой r ). |
|||||||||||||||||
Подставляя поле Er |
(7.77) в выражение для |
диффузионного пото- |
|||||||||||||||
ка ионов (7.68), находим |
|
dni |
|
|
|
|
|
|
|||||||||
|
|
Jim = Jem = −DA |
, |
|
|
|
|
(7.78) |
|||||||||
где |
|
|
dr |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
μe |
|
|
|
|
|
||
D |
A |
= μeDi + Zμi De = D |
|
|
|
1+ ZTe |
. |
(7.79) |
|||||||||
|
|
|
|
||||||||||||||
|
μ |
e |
+μ |
i |
i μ |
e |
+μ |
|
|
T |
|
|
|||||
|
|
|
|
|
|
|
|
i |
i |
|
|
Коэффициент DA носит название коэффициента амбиполярной диффузии. При выполнении условий (7.71) получаем
|
|
|
|
Te |
(7.80) |
DA = Di 1 |
+ Z T . |
|
|
i |
|
Если Te =T и Z =1 , имеем |
DA = 2Di , т. е. коэффициент ам- |
биполярной диффузии всего лишь в два раза больше коэффициента диффузии ионов, но значительно меньше коэффициента свободной (униполярной) диффузии электронов. Таким образом, амбиполярное электрическое поле существенно уменьшает направленную скорость электронов.
При наличии продольного магнитного поля в результате аналогичных преобразований вместо соотношений (7.78) получаем соотношения для радиальных диффузионных потоков поперек магнитного поля в виде
m |
m |
dni |
|
|
|
Ji = Je = −DA |
|
, |
(7.81) |
||
dr |
|||||
|
|
|
|
206
где DA получается из (7.79) заменой коэффициентов Dα и μα на соответствующие поперечные коэффициенты Dα и μα . В сильном магнитном поле при выполнении условий (7.75) имеем
|
|
|
|
|
Ti |
|
|
|
D |
A |
= D |
1 |
+ |
. |
(7.82) |
||
ZT |
||||||||
|
e |
|
|
|
||||
|
|
|
|
|
e |
|
При Te =T и Z =1 получаем, что коэффициент амбиполярной
диффузии поперек сильного магнитного поля равен удвоенному коэффициенту поперечной диффузии для электронов
DA = 2De = |
2De |
. |
(7.83) |
|
|||
|
β2 |
|
|
|
e |
|
При наличии градиента температуры в плазме в выражениях для диффузионных потоков в общем случае необходимо учитывать термодиффузионные члены. Анализ показывает, что для диффузионного потока электронов вклад термодиффузии связан с учетом
лишь парциального приведенного потока электронов hTe , а в уравнениях диффузионных потоков для ионов и нейтралов фигурируют
лишь парциальные тепловые потоки hTi и hTn . При учете термо-
диффузии следует иметь в виду, что в исходных уравнениях (7.17), на основе которых получаются соотношения для потоков (7.68), фигурируют градиенты парциального давления компонентов. С
учетом соотношения pα = nαkTα эти градиенты можно представить как
dp |
α |
|
|
dn |
α |
|
|
dT |
|
|
|
= k T |
|
+ n |
α |
α |
|
. |
|||
|
|
|
|
|
||||||
dx |
|
α |
dx |
dx |
|
|||||
|
|
|
В результате молярные диффузионные потоки электронов и ионов в слабо ионизованном газе запишутся в виде
J m = −D |
|
dne |
−μ |
n |
E − n D |
(1 |
+ αe |
) |
1 |
|
|
dTe |
, |
|
||||||
|
|
|
|
|
|
|
|
|
||||||||||||
e |
e |
|
dx |
|
e |
e |
e e |
|
T |
|
|
T |
|
dx |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
|
|
|
|
||
m |
|
dni |
|
|
|
|
|
i |
|
|
1 dTi |
|
|
|||||||
Ji |
= −Di |
|
|
+ Zμi ni E −ni Di (1+αT |
) |
|
|
|
|
. |
(7.84) |
|||||||||
|
dx |
|
T |
|
dx |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
207
Постоянные термодиффузии αTe и αTi в этом случае, как правило,
заметно меньше единицы, а для некоторых потенциалов взаимодействия частиц могут вообще обращаться в нуль. Выражения для
αTi могут рассчитываться по известным формулам для бинарной смеси газов, компонентами которой служат ионы и атомы. Выражения для постоянной термодиффузии электронов αTe будут обсуждаться в следующей главе.
208
ГЛАВА 8. КОЭФФИЦИЕНТЫ ПЕРЕНОСА В ЧАСТИЧНО ИОНИЗОВАННОЙ ПЛАЗМЕ
В настоящей главе мы применим методы вычисления коэффициентов переноса, использованные ранее для газовых смесей, для того чтобы оценить их значения в частично ионизованном газе. При этом для вязкости и теплопроводности плазмы вполне удовлетворительные результаты могут быть получены на основе метода, использующего концепцию средней длины свободного пробега частиц. Для определения электропроводности и вклада термодиффузии в плотность тока проводимости и тепловые потоки в плазме достаточно эффективным оказывается подход, основанный на использовании метода баланса импульса. Разумеется, существенное упрощение подходов дает в ряде случаев лишь качественное описание рассматриваемых явлений, поэтому там, где это оказывается возможным, наряду с приближенными приводятся точные результаты кинетической теории.
8.1. Электропроводность частично ионизованного газа
Согласно результатам предыдущей главы, плотность тока проводимости в плазме в отсутствие магнитного поля определяется выражением
j = σE = σ (E + |
1 |
pe ) , |
(8.1) |
|
ne e |
||||
|
|
|
где σ– коэффициент электропроводности плазмы
σ = |
nee2 |
= |
nee2 |
τ0 . |
(8.2) |
|
meν0 |
me |
|||||
|
|
|
|
Средняя эффективная частота столкновений с передачей импульса представляет собой сумму соответствующих частот электроннейтральных и электрон-ионных столкновений
ν0 = νei + νen . |
(8.3) |
По точности вычислений этот результат соответствует первому приближению Чепмена-Каулинга.
Общее выражение для эффективной частоты νek имеет вид
209
ν |
|
= |
4 |
n |
|
v |
|
Q(1) |
|
|
|
8kTe |
1 2 |
|
ek |
|
k |
e |
, |
ve |
= |
|
(8.4) |
||||||
|
||||||||||||||
|
|
3 |
|
|
ek |
|
|
. |
||||||
|
|
|
|
|
|
|
|
|
|
|
πme |
|
Среднее эффективное сечение столкновений Qek(1) определяется
при этом выражением (7.24). Конкретный расчет его требует знания дифференциального сечения рассеяния электрона на тяжелой частице либо проинтегрированного по углам рассеяния эффектив-
ного сечения столкновений с передачей импульса Qek(1)(v). Для
электрон-нейтральных столкновений из эксперимента определяется обычно зависящая от скорости (или энергии) электрона эффектив-
ная частота столкновений с передачей импульса νek (v) = v Qek(1)(v). В этом случае для определения νek можно воспользоваться выра-
жением (7.25). |
|
|
|
Q(1) |
|
Для электрон-ионных столкновений |
ν |
ei |
и |
могут быть |
|
|
|
|
ei |
|
рассчитаны до конца и определяются выражениями (см. формулы
(7.34) и (7.36))
|
4 |
|
|
(1) |
|
(1) |
|
|
|
Ze |
2 |
|
2 |
1 |
|
|
|
|
|
|
|
π |
|
|
|
|
|
||||||
νei = |
|
ni |
ve |
Qei |
, |
Qei |
= |
|
|
|
|
|
|
|
ln Λ . (8.5) |
|
3 |
2 |
|
|
|||||||||||||
4πε0 |
|
(kT |
)2 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
В единицах СИ имеем
νei = 3,64 10 |
−6 |
ni |
ln Λ |
|
|
|
|
|
(1) |
|
|
−10 |
ln Λ |
|
|
|
T 3 2 |
с-1 |
, |
|
Qei |
= 5,85 |
10 |
|
|
м2 |
, |
||||
|
|
|
T 2 |
||||||||||||
|
|
|
e |
|
|
|
|
|
|
|
|
|
e |
|
|
|
|
|
|
|
|
|
T 3 |
1 2 |
|
|
|
|
|
|
|
|
|
Λ =1,24 10 |
7 |
|
e |
|
|
|
|
|
(8.6) |
||||
|
|
|
|
ne |
. |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для слабо ионизованного газа можно пренебречь электронионным взаимодействием по сравнению с взаимодействием электронов с нейтральными атомами. Этому условию отвечает, очевидно, требование
210