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

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

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

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

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

Обратимся теперь к дисперсионному уравнению (2.39)

и формуле

для

тензора

диэлектрической

проницаемости

(2.40). Пусть v3

— компонента

скорости v

в направлении

волнового вектора к, а

и

— компоненты, перпендику­

лярные к.

Тогда дисперсионное уравнение запишется

в виде

 

 

 

 

 

 

C2fe2

det {;

IjL +

С02 l2±-

йр0

(

О н

dfо

 

О)2

 

 

dv

 

 

w

df°

b

 

v) (v - a;. ■-

k

 

 

v3k

— j d v | = d e t A = 0.

(2.55)

 

 

 

 

Обозначим компоненты тензора А через A^v. Тогда можно написать, что

— 6^v ^1-

°р0

ир0 Г

^U3^V3

c2A2

2

0)2

2 > (2-56)

 

 

 

где

 

 

 

 

с,

 

dfo

 

 

= ‘ 1

v»t>v~d^

■dv.

(2.57)

 

v3k 0)

 

 

Заметим, что тензоры А и С являются симметричными:

2 2

 

Г

U V

— Г

А

— А

(2.58)

 

°

 

'-'VLU

-^Liv

A V U ‘

Напомним теперь, что мы исследуем случай

ur « vT

с

и

поэтому можем

пренебречь

единицей по

сравнению

с

величиной с2к2/ со2,

что

физически означает

пренебреже­

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

Anv— — ~~z [fyiv (& 2 4- ^с) + к\С^у\, г Де ке —-— (2.59)

§ 2. Решение

Нри наличии многочаст,

коллект. корреляций 355

для ц и v, не равных

3,

и

 

 

 

 

 

 

^4зз =

~2 { —

 

(1 + ^зз) } •

(2.60)

Связь между электромагнитными и электростатиче­

скими волнами отсутствует при условии

 

 

А 13 =

А23 =

0

или

С13 =

С23 =

0.

(2.61)

Это условие можно также переписать в форме

 

 

 

 

 

d [(t>!> F0 (у3)]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2 -62)

где используются обозначения

 

 

 

 

 

F0 (v3) = j

dvi dv2,

<Hj> =

-

 

J

vJQdvi dv2.

(2.63)

Тогда достаточным условием отсутствия связи между

волнами, очевидно,

является

тождество

 

 

 

 

(у р

=

0

для всех

у3.

 

(2.64)

В этом случае дисперсионное уравнение для электроста­

тических волн принимает

вид

 

 

 

 

 

 

 

 

л\2

 

 

 

 

 

 

 

 

 

1 — J ^ ( 1 +

C33) = 0,

 

(2.65)

откуда с учетом

тождества

 

 

 

 

 

 

1 = j f„dv=

-

 

 

 

 

 

Jf

 

v3— u

dv

(2.66)

следует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

dfo

 

 

 

 

 

1-

Wpo

(*

3

dv3 dv = 0.

 

(2.67)

 

 

G>2 U

J

V3 —

U

 

 

 

 

Интегрируя по частям, можно показать, что этот резуль­ тат совпадает с ранее полученным дисперсионным соотно­ шением для потенциальных электростатических волн [см.

формулу (3.2.18)].

23*

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

При выполнении условия (2.64) дисперсионное соотно­ шение для электромагнитных волн может быть записано в виде

( | - + 1 + Си ) ( - J + l + Caa) —C J,= 0.

(2.68)

Чтобы решить это уравнение, следует сначала записать интегралы Сп , С22 и С12 следующим образом:

dXuv

 

(v3 — u r

dXllV

dvз

 

d%nv

 

 

dvз

,.

—j-------d v з

J v3— u

i?3—

+

 

'*“* 5

(Уз- и г)2+и? ’ (2‘69)

где %nv

=

(Vv-Vv) T0 (^3)1

затем выделить действительную

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

ряда оценок по порядку величины.

 

 

Перепишем выражение для

коэффициентов

 

'\IV '

((i^ v ) Fq (v3)) dv3

 

 

 

v3— u

 

 

 

 

 

<!Vv> Fo(v3) dv3

где p = 1,2, v = l , 2,

(2.70)

- l

(V3 — u)2

 

 

 

 

и учтем,

что | и |2 л: (у2).

Тогда

 

 

 

 

\и\2ж | v3 — u |2,

j

F0(v3)dv3= 1,

 

 

С ^ = 0 ( 1),

ц = 1 ,

2,

v = l , 2.

(2.71)

Из (2.67) видно, что

 

 

 

 

 

fc!„=0(*2c) < * ! c= o

( | | ) .

(2.72)

Абсолютное значение волнового вектора /сэм рассматри­ ваемой электромагнитной волны много меньше величины волнового вектора к эс соответствующей электростатиче­ ской волны. Формула (2.72) дает весьма полезную инфор­ мацию для оценки величины связи электростатических

§ 2. Решение при наличии многочаст,, коллект. корреляций

357

и электромагнитных волн. Для этого запишем дисперсион­ ное соотношение, сохранив в нем все члены, в виде

l + ' p " ( 1 + С ц )

*

2

г

* 2

г

Ж

 

° 12

Ж

0 *3

* 2

п

1 + Ж

( 1 +22

* 2

п

23

= 0. (2.73)

Ж

с ‘ 2

¥

°

 

kl

si

* 2

г

 

 

 

 

W

Ll3

Ж

° 23

 

 

 

 

Рассмотрим вначале те волны, волновые векторы кото­ рых близки к кдс, где-А:1с » o)|,0/(i>2). Используя оценку коэффициентов по порядку величины, детерминант (2.73) можно переписать в виде

- ^ + § ( 1 + С з з ) = 0 ( 1 т ) ,

(2‘74)

или

 

1 + Си (* , В) - £ = о ( $ ) .

(2.75)

С помощью разложения к — кэс + Ак, где кэс представ­ ляет собой решение уравнения (2.75) с правой частью, строго равной нулю, мы получим

т " = ° ( ! - ) -

<2-76>

Полагая дС33/дк = О (С33/к) = О (1/к) и принимая во внимание (2.72), найдем, что в соответствии с исходным предположением

(2-77)

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

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

ккс, т. е. к2 » к\. Тогда в дисперсионном соотношении (2.73) все коэффициенты по порядку величины будут близки к единице и, следовательно, членами связи в этом случае пренебрегать нельзя.

2.4.Электромагнитные неустойчивости

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

/о (v) = / 0 (— v).

(2.78)

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

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

оо

 

 

j vfo {v, #,

cp)dv = const,

(2.79)

j п2/о (v, #,

ф) dv = const,

(2.80)

о

 

 

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

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

ОО

 

j vnfo (v, ft, ф) dv = const

(2.81)

0

§ 2. Решение при наличии многочаст, коллект. корреляций

359

для всех значений п. В частности, для изотропии давле­ ния нужно было бы потребовать выполнения условия (2.81) при п — 4. Тем не менее во всех практических слу­ чаях при выполнении условий (2.79) и (2.80) распределе­ ние по скоростям в плазме вряд ли будет анизотропным, если она не приготовлена каким-либо особым способом. В практических случаях трудно ожидать, чтобы функция распределения обладала таким свойством.

Исследования Кана не дают ответа на все вопросы. Из них не следует, что плазма заведомо будет устойчива, если ее функция распределения удовлетворяет условиям (2.79) и (2.80). Эти условия являются лишь необходимыми. Кроме того, они не дают возможности получить резуль­ таты, если распределение по скоростям не является цен­ трально-симметричным. Однако Кан утверждает, что все подобные распределения неустойчивы.

Чтобы проанализировать указанные результаты, вспом­ ним, что мы пренебрегли корреляцией. Следовательно, приведенные выше утверждения справедливы, когда рас­ сматриваемые интервалы времени малы по сравнению с временем между столкновениями, т. е. когда [см. фор­ мулу (П.З) в приложении]

а > «го • (2.82)

Поскольку мы приняли, исходя из соотношения (2.72), что | со| = О {(йр0/с), и предположили, что | и | = = О (vT), то рассматриваемые явления обусловлены вол­ нами низких частот. Поэтому с учетом неравенства (2.82) получим

^ = ° ( - !r ) » J¥ -

(2-83>

Отсюда вытекает, что описанные выше неустойчивости могут проявляться только при достаточно низкой плот­ ности и высокой (но нерелятивистской!) температуре

плазмы.

Фрид [3] и Фюрт [4] исследовали физический меха­ низм усиления электромагнитных волн. Поскольку авто­ ров интересовал лишь основной механизм, а не проявле­ ние частных деталей, они изучали поведение системы электронов в простых экспериментальных условиях.

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

Начнем с рассмотрения начального возмущения маг­ нитного поля типа

Bz = Beihx, Вх = Ву = 0,

(2.84)

которое представляет собой волновое возмущение поля

внаправлении координаты х. Вектор магнитного поля направлен по оси z.

Вначальный момент времени нет движения электронов

внаправлениях z или х. В направлении у все электроны обладают одинаковой скоростью. Тогда функция распре­

деления электронов запишется следующим образом:

/о (v) = ab(vx) b(vz) b(vl aa) при a > 0.

(2.85)

Электрон, имеющий начальную скорость vy, в магнитном поле испытывает ускорение

vx =

— VyBz= co^Vy,

(2.86)

где coj, = eBjm — переменная циклотронная частота. Вследствие этого в направлении оси х возникает поток у-компоненты импульса. Производная по времени от этого потока определяется соотношением

mvyvx = — mvl<x)L.

(2.87)

Проводя усреднение с помощью функции распреде­ ления / = / 0 + Л (г, v; t) и объединяя члены первого порядка, получаем

д (vyvx)

a2aiL.

(2. 88)

Ft

 

 

В силу закона сохранения импульса этот поток (или дав­ ление) вызывает в соответствии с уравнением непрерыв­ ности изменение компоненты скорости (vy):

д {vy)

d{i>xvy)

(2.89)

dt

дх

 

Уравнение Максвелла, содержащее VxB, в пренебреже­ нии током смещения сводится к соотношению

Щ В ^ = po/j, = + ро е п _ ( V y ) .

(2.90)

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

корреляций 361

Исключая с помощью этого соотношения

(у ^ и з

уравне­

ний (2.88) и (2.89), найдем

 

 

 

 

ik

d*Bz

(и0ен_

9 (а2а>ь)

ik

С 0 р _

а2В ,

(2.91)

dfl

дх

С2

Отсюда следует, что начальное возмущение Вг нарастает по экспоненциальному закону с инкрементом (ор_а/с. Это согласуется с выводами, полученными Каном.

Выражение (2.91) не совпадает с результатами более общего исследования Вайбеля [5]. Тем не менее рассмо­ тренная грубая модель с ее частным случаем начального распределения все же может дать некоторое качественное объяснение неустойчивости.

Пренебрежение током смещения d E l d t в уравнении (2.90) эквивалентно пренебрежению единицей по сравне­ нию с величиной с2А:2/(о2 в дисперсионном соотношении. Это показывает, что токи смещения не дают вклада в раз­ витие рассматриваемой неустойчивости. Ее механизм заключается просто в том, что начальное возмущение магнитного поля приводит к такому ускорению электро­ нов, что возникающий поток импульса воздействует на начальное распределение токов в направлении усиления флуктуации поля. С физической точки зрения усиление начального возмущения поля можно рассматривать как сжатие (пинч-эффект) плазмы в токовые слои. Если при этом предположить, что компонента скорости vx не равна нулю, то должен возникнуть эффект суперпозиции, обусловленный движением частиц системы в направлении оси х. Этот эффект противодействует развитию неустой­ чивости, полученной в отсутствие движения в направле­ нии оси х.

§3. РЕШЕНИЯ С УЧЕТОМ ЭЛЕКТРОН-ИОННЫХ КОРРЕЛЯЦИЙ

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

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

Выше мы уже обсуждали вопрос об излучении, свя­ занном с одночастичными взаимодействиями. В данном параграфе мы исследуем процессы излучения и поглоще­ ния, обусловленные той частью ускорения частиц, кото­ рая приобретается при взаимодействии одновременно со всеми частицами плазмы. Эффективное поле частиц плаз­ мы обычно представляется в виде суммы коллективного поля, входящего в уравнение Власова, и поля, обуслов­ ленного корреляциями отдельных частиц и представляе­ мого остаточным членом правой части цепочки уравне­ ния ББГКИ. В принципе эффекты от обоих этих сла­ гаемых могут быть названы «тормозным излучением». Разумеется, вклад в тормозное излучение, обусловленный коллективным полем, представляет собой явление, кото­ рое не может встретиться при одночастичном рассмотре­ нии. С другой стороны, вклад в тормозное излучение за счет полей, обусловленных корреляциями, должен при­ сутствовать в расчетах, относящихся к одночастичным процессам взаимодействия.

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

В принципе мы должны были бы вновь обратиться к об­ щим уравнениям, рассмотренным в § 1, чтобы составить цепочку уравнений для величин /2, (Е)2 и (В)2, сделав при этом ряд предположений относительно корреляций третьего порядка / 3, (В ) 3 и (Е )3. Однако громоздкие вы­ кладки, приведенные в § 1, показывают, что подобная процедура малопривлекательна.

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

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

достаточную точность для многих практических прило­ жений.

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

когда выполнено условие со ^

мр_.

К сожалению,

как мы видели выше, в области и > с,

представляющей

практический

интерес, приходится

иметь дело с частотами со > сор_, что не позволяет исполь­ зовать уравнение Фоккера — Планка. Поэтому в ука­ занной области мы будем следовать методу, предложен­ ному Даусоном и Оберманом [6, 7].

3.1. Модель Даусона — Обермана

Применимость

модели

Даусона — Обермана ограни­

чена следующей областью частот и скоростей:

^

2п

In Л

st.

ю >

— « ®р- —д - ,

u^>vT.

Приведем ряд предположений, которые лежат в основе этой модели.

1. Исследуемая система, как всегда, представляет собой неограниченную пространственно однородную среду, со­ стоящую из электронов и равного числа ионов противо­ положного знака заряда. В среде присутствует слабое однородное электрическое поле Еоб- *"*. Предположение относительно однородности электрического поля означает, что область рассмотрения ограничена достаточно большими

длинами волн

или, строго говоря, рассматривается об­

ласть фазовых

скоростей,

превосходящих

среднеквадра­

тичную скорость частиц.

(т_/т+ — 0), а

распределение

2. Ионы неподвижны

их случайно.

 

 

 

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