ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 343
Скачиваний: 6
§ 3. Модели плазмы, с малым ß |
283 |
Тогда уравнение (78) принимает вид |
|
|||
d2i|) |
1 dip |
»m*ъ2 IP: |
|
|
dr2 |
dr |
r 2 |
|
|
|
|
r-\-a |
re+ (Я, г) R dR |
|
|
- 4леС[т J |
- лг_ (r, o ] - (84) |
||
|
[4Д2г2__(Д2_(_г2 _ а2)2] |
|||
|
|
|
i/2 |
Магнитное и электрическое поля того типа, который мы рассматриваем, приводят к появлению компонент дрейфовых ско ростей в направлении z, однако эти компоненты не включаются в уравнения. Из формулы (81) можно получить выражения для /•-компоненты дрейфовой скорости ионов и ее дивергенции:
== - -Т Ж Ь 5+ 4 |
Ѵ2Ф] *ітф. |
(85) |
|
divv+= — |
(-§§•) + |
Ѵ2ф] |
(86) |
Дрейфовую скорость нулевого порядка для ионов можно запи |
|||
сать в виде Ѵ+ = фо Ѵ+(г), |
где |
|
|
У+ (г) |
= г (Qm + |
Qe). |
(87) |
Предполагается, что функция распределения энергии ионов имеет вид 6-функции. Частота прецессии из-за градиента магнитного ноля равна
о _ сГ B z dB |
,q„. |
Ü M - ~ ^ W ~ d F ’ |
(88> |
где Т — энергия ионов. Электрическая составляющая частоты пре цессии, вызываемой электрическим полем нулевого порядка, равна
Г |
|
4лec (1 — Г) j r\+(r')r' dr', |
(89) |
где ц+ (г) — невозмущенная плотность ионов. |
В формуле (89) |
использовано предположение, что ионная и электронная плотности имеют одинаковую пространственную зависимость, т. е. г]_ (г) =
= Гг]+ (г), где |
Г — постоянная |
величина. С |
учетом |
соотноше |
|||||
ний (85)—(89) уравнение (79) принимает вид |
|
|
|||||||
|
дп_і_ |
ІСП: |
i f f |
^ |
+ |
~ Ѵ3ФJ = |
О, |
(90) |
|
|
dt |
||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G(r) = ^ - V +(r) |
■m |
( & м Ң- £2_g), |
|
|
||||
|
|
|
B z |
dN+ |
-N+ |
d ( BzY] |
|
||
|
|
|
В 2 |
dr |
|
dr \ B 2 } \ ' |
|
||
При решении |
уравнения (90) |
мы подставляем Д2ф из уравнения |
|||||||
|
Ѵ2ф = — 4лес I~п+— п_ -f- |
V2/i+J |
|
286 Гл. 7. Модели плазмы без столкновений
k = 4яес, а — ионный ларморовский радиус. Интегралы в форму
лах находятся численно в точках с использованием значений |
р"+ 1 |
ау ”+1, / = 0 , 1 , 2 Подынтегральные функции имеют |
осо |
бенности в граничных точках, но интегрируемы. В окрестности
граничных точек |
использовалась |
теорема |
о |
главном |
значении, |
а в оставшейся |
области, г — а + |
Д ^ 7і |
^ |
г Jr а — |
А, приме |
нялась формула трапеций. Для обеспечения нужной точности интегрирования использовалось не менее десяти точек, так что
на отрезке а умещалось не менее пяти точек г;-. Для малых значений г, когда область интегрирования содержит меньше десяти точек гj, требуется интерполяция.
Для разностной аппроксимации уравнения (96) используется следующее неявное разностное уравнение:
р П + і ----р п |
. |
- J ж = Т ъ (ѵГ1+ чЪ + HjW] +
щн, [уГ-бГ1+-т(-алм+ълГ1-с№!)] +
+ka)Hj [у? — 6? + -^- ( — ajy% 4+ bjyf - crf$-i)] -
Аналогично аппроксимируется и уравнение (97):
_у7-1 |
— jG j (рГ 1 + Pi) - |
HjV] - |
4 |
kaW i [ p r 1 - ѳ г 1 + |
|
At |
|||||
|
+ 4 |
( - ^ р м + ^ рГ |
- ^ р"-!1)] - |
||
|
- |
ka)Hj- i [ p |
? |
Ц- - |
( aj- pѲ ? ? + bjp71 + ]+ - cjp^) ] |
Написанные выше разностные уравнения можно представить в виде системы линейных алгебраических уравнений
|
Ajyr;: i + р |
г - B j y ^ + |
Cjy]_+i = |
d -, |
(102) |
||
- |
A j Р - 1 + |
B jpJ+1+ |
тГ 1 - |
Cjp™ = К Ь |
(1оз> |
||
где |
|
|
|
|
|
|
|
|
Aj = 22 ka*HjdjAt, |
|
|
||||
Bj ==y GjAt -\-^-ka2jHjAt -f-^ kajHjbjAt, |
|
||||||
|
Cj = — kajHjCjAt, |
|
|
||||
Dj = И;-у™+1 -f- p" -j- Bjy1} — |
+ |
HjW’jAt |
g kaj (öf -(- 8j+1) At,. |
||||
K nj = Ajpl, - |
Bjp] -f y? + Cjp?„, - |
HjVjAt + 4 |
kä) (&) + |
Ѳ"+1) At. |