![](/user_photo/_userpic.png)
книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы
.pdf§ 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 |
\ ~ |
4л |
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)
х) Для ссылок на формулы другой главы здесь и далее употреб ляется запись тремя числами: первое обозначает номер главы, второе — параграф и третье — порядковый номер формулы. —
Прим. ред.