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

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

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

§ 1. Приближение Власова

133

<Е)= - 2

j dy' '

J

 

х

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

X /(1> (vr

, vp ; ^)-)-Евнсшн- (1-4)

Здесь

через

Евнешн обозначено

внешнее

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

поле.

Уравнение

(1.3) в отсутствие корреляций, когда

 

 

 

 

S/lV

О,

 

 

(1.5)

совпадает

с

уравнением Власова [1]

 

 

 

 

5/<1>

3/<1)

 

5/d)

=0,

( 1. 6)

 

 

 

dt

д»г

■еа(Е>• dip

 

 

 

 

причем (Е)

 

дается выражением (1.4).

 

 

Условие (1.5) накладывает ограничения на примени­ мость уравнения Власова. В реальных системах это усло­ вие выполняется только приближенно. Для того чтобы обсудить, насколько существенный вклад дает член с пар­ ными корреляциями, приведем здесь некоторые резуль­ таты, которые будут получены ниже, при исследовании уравнений цепочки более высокого порядка. Мы покажем, что при определенных условиях вклад от члена, стоящего в правой части уравнения (1-3), порядка величины /(1>/тс, где тс — время между столкновениями электронов г)

(1.7)

Юр_ In Л

Здесь

СОт

/

4 ппе2 \ V* t д

зе 3/г

= 96,

(1.8)

\

т / ’

— (4лп)1^2<

 

 

 

где о»р_ — хорошо известная электронная плазменная частота, а Л — плазменный параметр, который пропорци­ онален числу частиц 6 в дебаевской сфере.

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

т < т с.

(1.9)

!) См. для сравнения определения, приведенные в приложении.

134 Гл. 3. Неравновесные состояния кулоновской системы

Если плотность плазмы значительно ниже критической, т. е. когда Л^> 1, мы можем пренебречь корреляциями вплоть до характерных времен порядка времени электрон­ ных плазменных колебаний тр_. Это очень важное заме­ чание, поскольку большинство коллективных явлений, описываемых уравнением Власова, имеют характерные времена порядка тр_.

Чтобы оценить возможности применения приближения Власова к реальным физическим системам, в табл. 1 при-

Таблица 1

Характерные параметры для некоторых плазменных сред

Плазма (полностью

П, СМ-3

е/хв, к

ионизованная)

Межзвездный газ

1

102

Солнечная корона

106

10«

Разряд

низкого

10»

5-10*

давления . . .

Разряд

высокого

Ю15

5-103

давления . . .

Термоядерная

1019

108

плазма . . . .

Л

V е

Тс> с

 

 

10?

10-4

10

1010

10-2

10

4-105

4-10-ю

2-10-в

10

4С 1 О

3-10-12

108

Ю-12

10-6

ведены характерные параметры для некоторых типичных плазменных сред.

Заслуживает внимания тот факт, что Власов сам пред­ ложил некоторую гипотетическую систему, которая точно описывается уравнением Власова, — так называемую дисперсионную модель [2]. Эта модель основана на беспре­ дельном разбиении каждой частицы на равные части, так что га-»- оо. Процесс такого разбиения подчиняется следующим условиям:

lim (гае) = const,

lim (nm) = const.

(1.10)

71-*-оо

71—►оо

 

Из выражений (1.10) следует, что elm = const, и посколь­ ку тепловая скорость — постоянная величина, 9 -»- О,

§ 1. Приближение Власова

135

Совершенно очевидно, что при условиях (1.10) урав­ нение Власова и время плазменных колебаний тр_ оста­ ются неизменными. То же самое справедливо и для дебаев­ ской длины, которая является произведением времени плазменных колебаний и тепловой скорости. Из всего ска­ занного следует, что

lim

б == lim Л = оо.

(1 .11)

П -* - о о

П ~+ СО

 

Из этого условия и выражения (1.7) получаем тс ->оо, что означает отсутствие корреляций между частицами.

1.2.Общие свойства

Вдальнейшем для простоты будем рассматривать одно­

компонентную систему,

а поэтому

в

уравнениях (1.4)

и (1.6) опустим индексы р и г .

При

этом, чтобы избежать

расходимости электрического

поля

в

пределе N —>- оо

и V —*- оо, предположим,

что в системе существует ком­

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

У Евнешн =4япсе = J /(1>(г\v'; t)dv'dv'. (1.12)

Отсюда для однокомпонентной системы имеем следующие соотношения:

5/<i>

,

a/<i)

^ < е > - 4 ^ = о,

(!-!3)

d t

1 '

дт

г m ч" / ' д \

<Е)= — J

 

 

{ j /аЧ-\ V ; О * ' - * } *'•

 

 

 

 

(1.14)

1. Уравнение Власова является уравнением типа Лиувилля в р-пространстве:

Т =

+

=

(1Л5)

!) Здесь мы заменили зависимость от импульса на зависимость от скорости, как это общепринято в литературе.

136 Гл. 3. Неравновесные состояния кулоновской системы

Это весьма удивительный факт, поскольку при выводе соотношений (1.13) и (1.14) мы пренебрегали только взаимодействием между отдельными частицами, а не кол­ лективным взаимодействием всех частиц. Уравнение (1.15) означает, что функция распределения /(1>остается посто­ янной вдоль траектории частицы.

2. Без ограничения общности можно представить функ­ цию распределения /(1) в виде функционала /(1)(^4V. Под­ ставляя этот функционал в уравнение Власова, полу­ чаем

V

 

Отсюда видно, что любая функция

от интегралов дви­

жения A v есть решение уравнения

Власова. Например,

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

3. Легко найти производную функции Н по времени, если предположить, что /(1) удовлетворяет периодическим граничным условиям в произвольно большом, но ограни­ ченном объеме р-пространства:

^ = т г “ М /',’1п,° ’* ' гу=

= ji^ L in /< 1>drdv+ j ^ - d r d \ = 0. (1.17)

Отсюда следует, что функция Н, а вместе с ней и энтропия S — хвН являются постоянными. Уравнение Власова и его решения обратимы во времени, что легко показать, заменяя в соотношениях (1.13) и (1.14) время t на —t

искорость v на — v.

4.Если распределение /О является положительно

определенным в начальный момент времени, оно остается положительно определенным и для всех остальных момен­ тов времени.

Допустим, что /(И принимает отрицательные значения, тогда должна существовать точка (г0, v0), в которой в мо­ мент времени t0функция /(1)в первый раз становится нулем.

§ 1. Приближение Власова

137

Поскольку в окрестности этой точки функция /(1>положи­ тельна, в самой точке должны выполняться следующие соотношения:

dV/<!>

О для всех v < a ,

dqt ... dqj

(1.18)

( A q . ^ y y ^ A Y 1’ > 0 ,

где

ч = ( ; ) .

о м )

С другой стороны, уравнение Власова может быть записа­ но в операторной форме

о>/<1)

t)

d / < i >

( 1. 20)

dt — с (q,

dq

откуда получаем производные по времени более высокого порядка

d v / a >

(C( 4 . 0 - ^ ) V/(l)-)-члены,

содержащие производные

atv

 

—— более низкого порядка, чем v.

( 1. 21)

Учитывая (1.18), получаем

 

 

 

 

о

для v < a ,

( 1.22)

d v / < D

 

для

v = а (положи­

~ dtv

* ( c ( q , l ) . ^ - ) ° / a '

 

тельно определена),

 

 

 

 

причем следует отметить, что оператор d/dq действует так­ же на с (q, t). Однако соответствующие слагаемые не дают вклада благодаря первому соотношению (1.18). Из соотно­ шений (1.21) и (1.18) следует, что daf lVdta — положитель­ но определенная функция.

5. Любая функция вида

/(1) (г, v, t) = h (у)

(1-23)

является решением уравнения Власова при условии, что выполняется соотношение (Е) = 0. Это соотношение безу­ словно справедливо в линейной теории Ландау, где имеет

138 Гл . 3. Неравновесные состояния кулоновской системы

место следующее равенство:

j/i(v )d v = nc.

(1-24)

1.3. Линейное приближение

Нелинейное уравнение Власова можно решить только в нескольких частных случаях. Поэтому до сих пор основ­ ные усилия были направлены на решение линеаризован­ ного уравнения Власова. Представим функцию / (1) в виде

/<*> = (0)/<1) + <D/«\

(1.25)

предполагая, что

 

(1)/<i) < <°)/(1,.

(1.26)

Будем также считать, что функции (0>/(1) и (1>/ а> являются гладкими, так что дифференцирование и интегрирование не меняет порядок их величины. Функция распределения нулевого порядка <0>/(1) предполагается однородной и ста­ ционарной, а внешние поля — равными нулю. В осталь­ ной части этой главы будем пользоваться упрощенными обозначениями:

<0,/ а, =

/0>

а>/ (1, = /.

(1.27)

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

^ + v

df_ + — Е. 9fо

= 0,

dt ^

dr

1 m

д \

(1.28)

Е = —Уф,

 

 

 

 

У2ф = 4ле j fdv.

Эти уравнения образуют систему интегро-дифференциаль- ных уравнений для /. Для решения такой системы уравне­ ний применяются различные методы.

§2. РЕШЕНИЯ ЛИНЕЙНОГО УРАВНЕНИЯ ВЛАСОВА

2.1.Метод собственных решений

Для того чтобы найти собственные решения системы (1.28), отделим зависимость от пространственных коорди­ нат и зависимость от скорости и времени, т. е. представим /

§ 2. Решения

линейного уравнения Власова

139

в виде

= / (г) / (v, t).

 

/

(2.1)

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

Подставляя (2.1) в (1.28) и применяя хорошо известную процедуру вычислений (метод разделения переменных), получаем функцию, характеризующую пространственную зависимость функции распределения

/ (г) — A exp (X - г).

(2.2)

Чтобы обеспечить ограниченность / (г) при всех значениях г, необходимо выбрать постоянную разделения X чисто мнимой величиной, т. е.

/ (г) = A exp (fk-r),

(2.3)

где к — действительная величина. Известно, что решения в виде нормальных мод позволяют к любому явлению в про­ странстве координат г применить фурье-анализ:

й ( к , v, t) =

(2я)-я/а [ й (г, v, t)e~ik"dr,

Q (г, v, t) =

,

(2.4)

(2п) - 3/2 j

й (k, v, t) eik T dk.

Фурье-преобразование в (1.28) дает следующие диффе­

ренциальные

уравнения:

 

 

( 4 - +

ik .v)/< k, V, 0 =

О к-9- ^

(2.5)

И

 

 

 

 

Ё (к, t)= -ik<5(k, t),

 

 

Ф(к, t) = ^g- { /(к , v',

t)d\'.

(2‘6)

Уравнение для поля Е можно также получить непосредСтвенщ» из (1.14) и (1.25), используя теорему свертки

140 Гл. 3. Неравновесные состояния кулоновской системы

из фурье-анализа и соотношение

/

1

\ ~

1

(2.7)

\

| г —г' |

/

&2

(2я)3/2

 

Очевидно, что составляющая скорости, перпендику­ лярная волновому вектору к,

vj_ = v к (v • к), к = -|-

(2.8)

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

Q (к, и, t) = j Q (к, v, t) 6 к • v) dv.

(2.9)

Здесь через и обозначена компонента скорости, параллель­ ная вектору к. В результате получаем

[ ж + Л и ) fa .

[ j fa .

(>**']

 

 

 

( 2. 10)

Чтобы найти собственные значения задачи, восполь­ зуемся соотношением

Г(к, и, t) = f(k, и)е~™\

(2.11)

которое приводит к

=

J f ( u')du'>

(2.12)

где ир — комплексная фазовая скорость, равная

ир = ~к '

(2-13)

Здесь и в дальнейшем для простоты записи мы не всегда будем отмечать зависимость от величины к. Фазовая ско­ рость, вообще говоря, является комплексной величиной вследствие комплексности частоты со. Причина этого заключается в том, что уравнение (2.12) не является само­ сопряженным. Из уравнения (2.12) непосредственно нахо­

дим функцию распределения / (и), поскольку интеграл в правой части этого соотношения представляет собой постоянную величину, определяемую нормировкой самой

§ 2. Решения линейного уравнения Власова

141

функции распределения, которая в линейной теории про­ водится обычным способом.

В результате интегрирования уравнения (2.12) полу­

чаем

дисперсионное

соотношение

 

 

fc2 = ig fl Г dfo {и)!ди

2

 

 

т J (и—ир)

'

которое определяет

фазовую

скорость

ир (к) и частоту

<й (к)

как функции волнового

вектора к.

В общем случае

интегрирование (2.14) представляется сложной пробле­ мой, поскольку значение этого интеграла не является одно­ значно определенным из-за наличия полюса при и = ир. В дальнейшем мы увидим, как решается эта проблема.

Рассмотрим сначала два примера, в которых такой проблемы не возникает. Это имеет место либо в случае, когда df0 (и)/ди = 0 вблизи и = ир, либо когда ир — комплексная величина. Такое положение имеет место во многих практически интересных случаях.

Прежде чем перейти к обсуждению этих задач, опреде­

лим функцию

распределения F (и) следующим образом:

/„

(и) = F (и) j /о (u') du' = F (и) п.

(2.15)

Предполагая, что полная функция распределения вклю­ чает различные сорта частиц у (вследствие многокомпо­ нентное™ плазмы), запишем

/ov (и) = Fv (и) j /ov (и') du' = Fv (и) nv.

(2.16)

Плазменная частота каждой из компонент равна

®2,v= mv

Тогда дисперсионное соотношение (2.14), обобщенное на случай многокомпонентной плазмы, запишется в виде

dF: ~ l j Ldu = e (k,<o) = 0. (2.18)

V = 1

- о о

Конечно, теперь имеется столько кинетических уравнений, сколько компонент частиц плазмы, но все они приводят к одному и тому же дисперсионному соотношению (2.18).

142Гл. 3. Неравновесные собтОяниА Кулондвской системы

Впоследующих параграфах мы обсудим физичес­ кий смысл соотношения (2.18) для различных важных примеров.

2.2. Моноэнергетические пучки частиц с заданным направ­ лением скорости

Рассмотрим систему пучков заряженных частиц, кото­

рые движутся в направлении вектора к с различными, но строго определенными скоростями:

Fy (и) = б щ).

(2.19)

Будем считать, что ни один из этих пучков не обладает тепловым разбросом скоростей. Подставляя (2.19) в (2.18) и используя соотношения (2.1.3) *), получаем

*! = 2 к ^ Ь

г -

<2-20>

V

 

 

Можно показать, что условие /

(и)

/ 0 (и) выполняется

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

1. Uy = 0. Этот случай, когда скорости пучков всех заряженных частиц равны нулю, является особенно про­ стым, но до некоторой степени представляет лишь теорети­ ческий интерес. Дисперсионное соотношение в этом случае принимает вид

(O' = 2 < = сор 0 '

(2.21)

v

 

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

_ Зсо

(2.22)

х) Для ссылок на формулы другой главы здесь и далее употреб­ ляется запись тремя числами: первое обозначает номер главы, второе — параграф и третье — порядковый номер формулы. —

Прим. ред.

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