Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы

.pdf
Скачиваний:
17
Добавлен:
22.10.2023
Размер:
12.01 Mб
Скачать

364Гл. 6. Взаимодействие электромагн. полей с сист. многих частиц

3.Магнитное взаимодействие частиц пренебрежимо мало, поскольку оно по порядку величины в c2/v2 раз меньше величины взаимодействия, определяемого элек­ трическим полем.

4.Движение электронов достаточно хорошо описы­ вается уравнением Власова. Поэтому на рассматриваемые частоты налагается требование

.V

2 л

In Л

.г, . I

®

 

(З Л )

где те — время между

столкновениями.

 

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

Запишем уравнения Власова для электронов (т_ —

— т, е_ — е)

3 / 1 1 ’ ,

д1<У

ду

О,

(3.2)

д% ^

дт

АФ =

8"1j / ? > d v - * 2 6 ( r - r , ) ) .

 

(3.3)

з

Здесь через гj обозначены координаты ионов. Преобразуем теперь эти уравнения, перейдя к системе

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

R = r -----~ Е 0е~ш ,

V = v + —

Е0е -’“‘.

(3.4)

т о ) 2 и

таз

1

'

Соответственно соотношения для преобразования произ­ водных запишутся следующим образом:

§ 3. Решения с учетом электрон-ионных корреляций 365

( л ) г ,1 - ( е т ) в , г

(3,5)

a i M - k l ,

В результате уравнения Власова преобразуются к системе

 

 

<5/1

a / i 11

,

е

<5Ф

а /! 11

= 0,

 

 

 

dt

 

aR

 

т

aR

 

a v

(3.6)

 

АФ =

е~хе [ j /I1' d\ -

2

б (R— %е~ш -

 

г/)] .

 

 

 

 

 

 

 

з

 

 

 

 

 

 

Здесь использована сокращенная запись

 

 

 

 

 

 

1 =

 

 

 

Е„.

 

 

 

(3.7)

Линеаризация уравнений

(3.6) сводит эту систему к виду

[см. (2.15)

и (2.17)]

 

 

 

 

 

 

 

 

 

 

 

 

df1-

V-

a/j_

 

 

е

a<x>!

а /0_

= 0,

 

 

dt

 

aR

 

 

 

aR

 

av

 

Aa>1= -k

 

 

 

 

 

 

 

 

 

 

 

 

(3.8)

'+ J A-

 

2

6 (R -

 

6e_to< -*>)]■

Предполагается,

 

что

 

 

3

 

/0_

имеет вид

 

функция

 

 

 

 

 

п_

 

 

 

I

 

V*

 

 

(3.9)

 

 

 

 

(2я)3/2

6ХР ^

 

2Vi

 

 

 

 

 

 

 

 

где Fr —среднеквадратичная

скорость

 

электронов.

Применим теперь к уравнениям (3.8) преобразование

Фурье

по пространству R:

 

 

 

 

 

 

 

 

 

 

ah-

 

ik• V/i_+ * ~

®ik • - =0,

 

 

dt

 

 

'*

'

 

т

 

 

 

 

 

 

*2Ф1 =

— ^

{ (2я)3/2 и_6 (к) +

 

 

 

 

 

 

(З.Ю)

+ [

dV — (2я)_3/а 2

ехр [ гк(%e~ш + г,)] } .

з

(Вследствие предположения о нейтральности невозму­ щенной плазмы первый и третий члены в преобразованном

366 Гл. 6. ЁзаимодействиёэлеКтромйен. пдлейссибт. Многихчастиц

по Фурье уравнении Пуассона для к == 0 взаимно сокра­ щаются.)

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

L = i \ - к, Л = — | k - ( ^ ) O i [ / i - ] ; (3.11)

тогда рассматриваемое дифференциальное уравнение мож­ но записать в общем виде

 

 

( 4 i + L )t‘- = R -

 

<ЗЛ2>

Однородное

уравнение имеет

решение

St = exp (— Lt).

Примем теперь, что fi- = Stg. Тогда получим

 

 

 

 

t

 

 

 

 

? i_ = j St_vR{t')dt'.

(3.13)

Потребовав,

чтобы

f1_ ( t = — oо) = 0, и

использовав под­

становку x —t — t',

найдем

 

 

 

 

U- = j SXR Ifi-

T)]dr-

(3.14)

 

 

о

 

 

 

Подставив это выражение во второе уравнение системы (3.10), получим, что

 

00

/<2Ф ,=

0 1( t - x ) d x ^ e - ^ k . - % d \ ~

 

о

 

— (2л)_3/г 2 ехр[ — гк- (1е_,’“‘ + г ^ )]|. (3.15)

3

Интегрирование по скоростям дает

- fo o

кп^

(2я)1/2 V2

j в lXkV dV еХР (

2У| )

-

ixk2

n_

 

 

 

(2я)1/2VT —оо

- гт/с2и_ exp [ — (Ут^т)2] • (3.16)

§ 3. Решения 6 учетом элёктрон-ионных корреляций

36?

Отсюда сразу получаем

оо

Ф 1 = — ft»p_ j exp [

j

( F t A:t) 2J

(t — x) d t +

 

o

 

 

 

 

+ (2я)"3/2(

^ )

2 ex p [ - ik .( le - ‘“‘+ rj)].

(3.17)

 

 

i

 

 

Строгое определение проводимости существует лишь для предельного случая Е 0 -*■ 0, т. е. для | ->■ 0. Это позволяет разложить последний член в правой части выражения

(3.17):

2 exp [—ik‘(%e~iat + г;)] =

7

 

= (1 — ik.-\e~m ) 2 exp (— гк-r,).

(3.18)

i

 

(Вообще говоря, правомерность такого приближения нуж­ дается в доказательстве, так как вследствие сингулярно­ стей, возникающих в ионных членах, оказывается суще­ ственным вклад от очень больших значений к . Детальное исследование данного вопроса показывает, что предполо­ жение Ik < 1 дает правильные результаты в низших поряд­

ках разложения по

справедливые для

любых зна­

чений к . )

(3.18) состоит из

постоянной

Поскольку выражение

и переменной по времени части, предположим для функ­ ции Ф4 аналогичное представление в виде Ф* = Фс +

+ Ф/ш (0- Подстановка этого выражения в (3.17)

приво­

дит к выражению для постоянного слагаемого в виде

Фс=

е

2 exp (— ik-rj).

(3.19)

(2я)3/2е0 (fc2 + £0|,_/F|)

3

 

 

 

 

Из структуры уравнения (3.17) также очевидно, что функ­

ция Фйш = С exp ( ~ ia t) должна удовлетворять осцил­ лирующей части этого уравнения. Используя это обстоя­ тельство, найдем

k'l е - ш 2 ехр(— г'к-Г;), (3.20)

Ф*ш (2л)3/2 г 0к2 Ако

368 Гл. 6. Взаимодействие электромагн. полей с сист. многих частиц

где

AftU) = l + to^_ j exp^tcor—- ^ p ^ - jrd T .

(3.21)

и

 

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

ния Пуассона. В результате выражение для Ф ^ изме­ нится на множитель /с2/ (к2 — со2/с2), а выражение для Фс останется неизменным. Однако асимптотические формулы для А к(Лпоказывают, что основной вклад в общее выраже­

ние для Ф4 дают значения к « a>/VT со/с (см. первое предположение), так что величина корректирующего мно­ жителя близка к единице.

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

Потенциал, создаваемый в точке г в лабораторной системе координат, совпадает с потенциалом в точке R

в введенной нами движущейся системе координат, и по­ этому он определяется формулой для обратного преобра­ зования Фурье:

®i(r) = (2я)-3/г jeik,rexp (ik -|e-itot) (Фс-)-Фйт) dk. (3.22)

Подставляя в (3.22) выражения (3.19) и (3.20) и опуская

члены

высших

порядков по

имеем

 

Ф‘ W =

(2я)3 80 Х

 

 

х 1

{ № +

+ гк-&~Ш ( *2 + a>a_/Ff, ~ Ш

^ ) } Х

 

 

X 2

ехР [&*(г— г./)]

(3.23)

з

§ 3. Решения с учетом электрон-ионных

корреляций 369

Результирующее среднее поле, действующее на ионы,

тогда равно

 

<Е> = —

(з-24)

X

 

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

Используя «приближение случайной фазы», которое

вытекает из предположения о случайном характере рас­

пределения ионов,

 

2 е х р [гк>(гг г,-)] = и+

(3.25)

i,3

 

и выражение (3.23), из (3.24) получим следующее выра­

жение

для

средней

напряженности

электрического поля:

® =

^

1

к ^

) ‘ - “ ( и ч ^

- - и

Ь

) Л - (3-26)

Здесь

также

учтено, что первый

член

под

интегралом

в выражении (3.23) не дает вклада в интеграл (3.26).

Интегрирование

по угловым координатам в (3.26)

является тривиальным и сводит тензор к)(к к единичному тензору, умноженному на (4я/3) к2.Таким образом, резуль­ тирующее выражение для средней напряженности поля можно записать в виде

(Е> = 6л2еп -iat \ **( fc2+ co“ !V\ №Ahc -)d&. (3.27)

Обычный метод обрезания при &макс = 2ле0тУт1е2 устра­ няет расходимость для малых значений прицельных пара­ метров.

В лабораторной системе координат плотность тока, который в данной модели создается только движением электронов, определяется соотношением

j = —e j v / t_dv.

(3.28)

2 4 - 0 1 2 9 1

370 Гл. 6. Взаимодействие электромагн. полей с сист. многих Частиц

Дифференцируя его по времени и используя уравнение Власова, получаем

W = e W - J v) (

v

ф 1) X

 

х

j ^ ) ( vdv- (3-29)

В соответствии с предположением (1) здесь можно прене­ бречь первым членом в правой части. Тогда интегрирова­ ние по частям оставшихся членов дает

3j

е2

аф4

\

(3.30)

-rf = П- —

dr

j*

dt

т

 

Усредняя по всем возможным ионным конфигурациям, как это было сделано при выводе формулы (3.27), полу­ чаем

^ = » - - £ ( Е0е - “Ч ( Е » .

(3.31)

Вспомним теперь, что коэффициент электропроводно­ сти а (со) связывает электрическое поле и плотность тока

j (со) = о (со) Е0 (со).

(3.32)

Сравнивая это соотношение с преобразованием Фурье уравнения (3.31), получим выражение для комплексной проводимости

ст =

г'е°

f1 — 7 И ) ’

(3-33)

где

 

 

 

Ь

 

 

 

"мако

 

 

1 № = 6я2е0тш2 |

(

к2+ ш2_/F^, ~ №Aka> ) к* d k

(3 -34)

— малая в рассматриваемой области частот величина по­

рядка 1/(ОТс.

Удельное сопротивление плазмы описывается соотно­ шением

р = т = - г й т 1 1+ / М1-

<3-35>

В предельных случаях, когда со/сор_ <С 1 или соответствен­ но (о/(ор_ 1, величина р (3.35) может быть вычислена

 

<ji 3. Решения с учётом элёктрон-иЬнких корреляций

371

аналитически.

В случае

to

сор_

 

 

 

Р =

{ t

1

( 1

8 " )

2 m 0mV»T ]

+

 

 

 

2 \V2

 

'е2со?,_

 

2Fy/cMaKc )2- l ] }

.

(3.36)

 

* ( n )

24леот^1'«>

0)n_

 

 

 

Для (о > ®p-

 

 

e2Ci>|_

 

 

 

 

 

 

V*

 

 

 

p = - = r { [ 1+ '2")

24ne0mV^U) ]

 

 

 

+г(4Г

e2(Op_

 

 

 

(3.37)

24яе0т1/ |.сй

 

 

 

 

 

 

 

 

 

где С

0,577

— постоянная Эйлера.

 

 

 

Заметим, что авторы работ [6, 7] рассчитали числен­

ными

методами

 

величину р

во всей

области

изменения

Ф н г. 27. Удельное сопротивление плазмы Re [р (со)] в области частот со > vc [6].

Значения

Re[p(co)] нормированы на

величину р0 =

ЕУ32 ег/3т~ (4ле0)*] X

 

X (т./в)*/* яг 5,8- lO-seV'Om,

где значения в

даны в эВ.

частот..

Их

результаты

для зависимости

действительной

и мнимой

(реактанса)

частей р приведены на фиг. 27

и 28. Из фиг. 27 видно, что действительная часть удель­ ного сопротивления на частотах ниже плазменной яв­ ляется постоянной величиной. Это также непосредственно

24*

372 Гл. 6. Взаимодействие электромазн. полей с сист. многих частиц

следует из формулы (3.36). Сравнение данных резуль­ татов с результатами, которые мы получим в следующем параграфе для низких частот со сор_, показывает, что на верхнем пределе частот они совпадают. Это вполне

понятно,

так как

наши результаты верны и в области

1/тс <

со

< ш р_.

 

Приведенная функциональная зависимость обнаружи­

вает

небольшое

возрастание удельного сопротивления

на частотах чуть

выше плазменной. Даусон и Оберман

Ф и г . 28. Реактанс плазмы 1ш [р (ш)] в области частот

ы > vc

(в тех же единицах, что и величина Re [р (со)] на фиг. 27)

[6].

В приведенном графике ив величины Im [р (со)] вычитается реактанс свобод­ ных электронов—ш/еоИр-.

показали, что это возрастание связано с возбуждением продольных плазменных колебаний. Убывание удельного сопротивления в области высоких частот можно объяснить уменьшением эффективного поперечного сечения столкно­ вений: столкновения с большими значениями прицельного параметра представляют собой медленный процесс и не дают вклада в высокочастотные поля. Таким образом, максимальное эффективное значение прицельного пара­ метра составляет величину порядка VT/ со.

Формула (3.37) и соответствующая ей зависимость на фиг. 28 показывают, что изменение реактанса можно рас-

§ 3. Решения с учетом электрон-ионных корреляций

373

сматривать как результат изменения эффективной массы электронов вследствие взаимодействия с ионами. Заме­ тим, что зависимость реактанса также обнаруживает небольшой пик вблизи плазменной частоты, что объясняется теми же причинами, что и возрастание действи­ тельной части удельного сопротивления.

Теперь, после того как мы вычислили проводимость и удельное сопротивление плазмы с учетом корреля­ ционных эффектов, можно также решить дисперсионное уравнение (2.39)

det ( тг21_1_- 1-----— а ) = 0

(3.38)

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

 

Подставим

выражение

(3.33)

 

для

проводимости а

в дисперсионное уравнение.

 

 

 

 

 

получаем

 

Для

продольных электростатических волн

 

 

 

со2 =

со2_ [ 1 - / ( с о ) ] ,

 

(3.39)

а для поперечных электромагнитных волн

 

 

 

со2 = fc2c2

+

со2_ (1

-

/

(со)).

(3.40)

Принимая во внимание неравенство

 

| I

(со) | <С

кото­

рое

уже

было

использовано

выше,

находим

из (3.39)

для

потенциальных электростатических

волн

 

 

 

 

co = cop_ [l

— |/( с о )]

,

 

(3.41)

а для поперечных волн, согласно (3.40), имеем

 

 

 

со - (А*С*+

K - f 2- 4

 

1 Н -

(3-42)

Таким образом получаем выражение для коэффициентов

поглощения в виде

 

 

 

 

 

 

 

 

 

 

 

а ц =

со;,_ 1 ш

[

/

(с о )],

 

(

(3.44)

Соседние файлы в папке книги из ГПНТБ