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

Батыгин&co

.pdf
Скачиваний:
374
Добавлен:
10.02.2015
Размер:
7.66 Mб
Скачать

620

ГлаваXIV

где \{t) — скорость частицы в точке г = H(t). Для вычисления полного тока j(r, £) нужно умножить (5) на число частиц в объеме (dRo), обладающих начальной скоростью vo, и проинтегрировать по всем возможным значениям Ro, v0 :

j(r,t) = enJv(t)S(r-K(t))f(v0)(dv0)(dRo).

(6)

Начнем с интегрирования по координатам. Аргумент «^-функции зависит от Ro сложным образом, поэтому перейдем кновой переменной интегрирования R = Ro + vof +Ri(Ro,£). Вычисляя якобиан преобразования с точностью до членов, линейных по Ео , получим

O(Hx,Ky,Hz)

После этого интегрирование по (dR) не представляет труда ипроводится с помощью формулы типа (П1.4). Подставляя под интеграл (6) выражение (7) и пренебрегая снова слагаемыми к • Ri в показателях экспонент, получим

j(M) =en/{vo +

*g + fVo(tE\2}/(vo)^vo), (8)

J {

т(и - k • vo)

m(w - k • v0) )

где E = Eoexp[i(k • г — wt)]. Точка k • vo = и не является особой точ-

кой

подынтегрального выражения, так как предполагается, что /(vo) = 0

при

к • vo =и>. Поэтому можно произвести

разложение по отноше-

нию

vo/vy = kvo/ijj, предполагая характерные

скорости частиц малыми

по сравнению с фазовой скоростью волны. Это позволяет представить (8) в виде

-/{•

ге(к • E)v0

+ тпш

Предполагая, что f(vo) не зависит от углов, получим

§ 2. Коллективные движения в плазме

621

Здесь

Е ( к ' Е ) к Е Е Е

v2 =2n I vjf(v)v±dv±dvll.

(11)

В случае распределения Максвелла v2 = Т/т.

Из выражения (10) находим тензор проводимости:

те2

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

', (13)

будем иметь

еа0(ш,к) = е±а0--^-\+е\\-^-,

(14)

где

,2 _ 4тгпе2

т

_

По сравнению сослучаем отсутствия теплового движения (v?, =0)возник новый важный эффект — зависимость е от к,пространственная дисперсия. В связи сэтим диэлектрическая проницаемость стала тензорной величиной. Зависимость е от к объясняется тем, что ток в некоторой точке создается частицами, приходящими из соседних областей, поле в которых не равно полю врассматриваемой точке. Пространственная неоднородность поля вместе степловым движением частиц и приводят к пространственной дисперсии диэлектрической проницаемости.

Поведение резонансных (к • v0 = ш) частиц рассмотрено в этой задаче недостаточно корректно, в связи с чем при расчете утеряна малая мнимая часть £ц, которая описывает передачу энергии от волны кчастицам (затухание Ландау, которое существует даже вбесстолкновительной плазме).

622

Глава Ш

 

Wv

3kvj

Wv

3v$

870. Vy, = -Z

+ -5—^ «

k'1"9 =

~HT^' ® отсутствие теплового

движения vg = 0. Таким образом, плазменные колебания распространяются

врезультате переноса электромагнитной энергии частицами.

871.а) р(г, t) = р(г,{0) cos wpt,

где 0 = %гг. В случае б), кроме плазменных колебаний плотности за-

ряда с частотой wp, происходит его релаксация из-за теплового движения частиц, р(х, t) = 0 при t > оо.

ПРИЛОЖЕНИЯ

1. (5-ФУНКЦИЯ

(5-функция Дирака определяется равенствами:1

8{х -а) =

оо при х = а;

6(х — a)dx = 1.

Интегрирование в(П1.2) выполняется попромежутку Дпроизвольной длины, содержащему внутри себя точку о.

(5-функция удовлетворяет следующим соотношениям:

5(х) = 6(-х),

(П1.3)

J

 

Jf(x)5(x-a)dx=

 

А

 

8(ах) = -гт8(х),

(П1.5

|а|

 

где f(x) непрерывная функция.

Трехмерная (5-функция определяется аналогичными соотношениями:

6(т - а) = 5(х - ах)5(у - ay)S(z - az)=

//(r)J(r - a)^=(^ ( a ) ' -ли а внутри объема К,

J

\0,

если а внеобъема V.

А

 

 

Здесь /(г) —непрерывная функция.

1 Математически корректное определение 5-функции требует обобщения обычного понятия функции (см.:И.М. Гельфанд, Г.Е.Шилов «Обобщенные функции и действия над ними», т. I,Физматгиз, 1958;В.С.Владимиров «Уравнения математической физики», «Наука», 1967). {-функция относится к классу сингулярных обобщенных функций.

624

Приложения

6(x —a,a)

dS(x —a,a)

С помощью J-функции можно описывать распределение в пространстве заряда точечной частицы. Объемная плотность такого распределения

р(г) = е6(т - а),

(П1.8)

где е —заряд частицы, а —радиус-вектор точки,

в которой находится частица.

Можно определить

также производную

от J-функции. Точный

ее смысл содержится

в формуле

 

/я*) дх

ах = — -да

 

д6(х - a)

J

df(a)

(П1.9)

 

 

 

которая получается интегрированием по частям. Аналогично определяются производные высших порядков:

 

 

f f(x)6(n\x-a)dx=

(-l)n/n(a).

H|

 

Сама J-функция

может

рассматриваться

4

 

как производная от функции, испытывающей

 

 

Рис.

135

в некоторой точке а конечный скачок 6. Ес-

ли /(а + 0) - /(а - 0) = 6, то

 

 

 

 

 

 

тг" = Ь6(х —а) + ограниченная функция.

(П 1.11)

 

 

ох

 

 

Наглядное представление о J-функции и ее производных можно получить, рассматривая график некоторой непрерывной функции 6(х —а, а), такой что / 6(х — a, a)dx = 1 (рис. 135). Параметр а характеризует ши-

А

рину промежутка, в котором 5(х — а, а) отлична от нуля. J-функция и ее производные определяются как пределы:

Six —а) = lim 6(x —а, а),

а—>0

д6(х - а)

=

,.

д6(х - а, а)

 

lim

 

ОХ

 

а->0

ОХ

и т. д.

1.8-функция

625

Свойства ^-функции приобретают многие несингулярные функции, зависящие от параметра, при определенных предельных значениях этого параметра. Например:

8(х) = limД • - a ^ - j ,

 

(П1.12)

 

<*->0 тг

ОТ + X2

 

 

 

6 ( х ) = lim к

.

й

 

 

 

fc-»oo 7Г

 

 

 

 

 

При А;и п целых

 

 

 

 

 

+ОО

 

 

 

 

цт±МШ=

у

5(х-пж).

(П1.14')

 

fc—oo ^ S1I1X

-^-^

V

'

V

'

п = — о о

Легко убедиться в том, что любое из представлений (П 1.12-П 1.14') согласуется со всеми определениями (П 1.1—П 1.5), а также с определением (П 1.9) производной от <5-функции. При вычислении интегралов вида f f(x)S(x —

— a)dxc помощью представлений типа (П 1.12-П 1.14) нужно иметь в виду, что предельный переход должен производиться после выполнения интегрирования, например при использовании (П 1.12):

J/ f(x)S(x -a)dx= limа-оУ/ f(x)S(x - a) dx.

Путем рассмотрения интеграла Фурье (или с помощью представления (П 1.13))можно получить еще одно полезное представление ^-функции:

оооо

S(x) = ^ j

eikxdk

= \

j cos kxdk.

(П1.15)

—оо

 

 

О

 

К J-функции близки две другие обобщенные функции 6+(х)

и 6-(х).

Они определяются равенствами, сходными с (П 1.15):

 

 

 

оо

 

 

6±(х)

= ^

/ e ± i k x

dk.

(П 1.16)

 

о

 

 

Функции 8+ и 6- связаны с ^-функцией соотношениями

±Р±,

626 Приложения

так что 6(х) = 6+(х) + 6-(х). Символ Р 1 в формуле (П 1.17) представляет главное значение интеграла:

О2

0 1

=У(а) ± ± що \7Ж

„,+

/ M J

2l / v

' 27Ге-о[1/ х - a

J х-а

у

 

ai

о+е

 

где ai < a < О2, £ > 0.

Если аргумент ^-функции является однозначной функцией независи-

мой переменной х, то имеет место формула

 

^ * - * ) .

( П 1 Л 8 )

где а» — корни уравнения tp(x) = 0.

 

2. СФЕРИЧЕСКИЕ ФУНКЦИИ ЛЕЖАНДРА

Сферическая функция порядка I, т, зависящая от полярных углов д, а, определяется формулой:

в которой целые числа I, т

удовлетворяют условиям I ^ 0, —I ^ т ^ I;

5т = (—1)т при т ^ 0, 6т

= 1 при т < 0. Через Pj m обозначен присо-

единенный полином Лежандра

 

(П2.2)

ь

ь

'Часто вместо обозначения Р f используют обозначение •£.

2. Сферические функции лежандра

627

где Pj(x) — обычный полином Лежандра. Онсовпадает с Р/т (х) прит = 0:

I. (П2.3)

Присоединенные полиномы Лежандра удовлетворяют дифференциальному уравнению

Приведем некоторые формулы, полезные при работе со сферическими функциями:

(0,0) =

= (-1)"

(П2.5)

+ l)P/+i(x) = (2* + l)xfl(x) - ifl-i(x),'

(П2.6)

Сферические функции с Z = 0, 1,2 имеют вид:

* v

/5

(П2.7)

'20=V4^

Y2,±1 =Ti/gsi

, У2,±2 = «/^si

Сферические функции образуют на поверхности сферы полную ортонормированную систему функций от д, а. Это означает, что

(П2.8)

'Символом п\\обозначено произведение всех последовательных целых чисел, имеющих ту же четность, что и п, от 2 до п при п четном и от 1 до п при п нечетном.

628

Приложения

где du =sin •ddfida элемент телесного угла, и что произвольная функция от д, а с интегрируемым квадратом модуля может быть разложена вряд

0?,a). (П2.9)

Коэффициенты a/m определяются формулами

aim = / ljm(0, <*)/(*,a) du.

(П 2.10)

Функции вида r~l~1Yim('d,a) иr'yjm(i9, a) называют шаровыми гармониками. Легко проверить с помощью (П2.4), что шаровые гармоники являются частными решениями уравнения Лапласа вовсех точках, кроме г = 0:

[ Г ^ а ) ] = 0 , A[rlYlm(#,a)} = 0.

Решением уравнения Лапласа является также суперпозиция шаровых гармоник с произвольнымикоэффициентами

оо

I

 

V = Е

Е («Ы"'"1 + blmrl)Ylm(d, а).

(П 2.12)

J=0 m=-l

 

Если г(г, i9,a) и г'(г', i9', a') — радиусы-векторы двух точек пространства, причем г > г' (см. рис. 7), то

Б = Т Ц т =

Х

= У -ТГГР/(СО87)- (П2.13)

 

= У ТГГ

R k - r ' | ^ - 2 1 ^ 0 0 8 7 + ^» ^ г

Функция 1/Д называется производящей функцией для полиномов Лежандра. Имеет место следующая теорема сложения для сферических функций:

5Z Yi{Vcc)Yf{d'cc')

(П2.14)

m=—l

Углы i9, a и д', а! входят в (П2.14) вполне симметричным образом. Подстановка (П2.14) в (П2.13) приводит к разложению:

1

Е , 0 , / ^ t + 1 ^ m ( ^ a ) ^ ( ^ , a 0 .

(П2.15)

3. Цилиндрические функции

629

Из формулы (П2.13) следует (если положить у = е~^) разложение:

1

^ i ^

(П2.16)

3. ЦИЛИНДРИЧЕСКИЕ ФУНКЦИИ

Цилиндрические функции Zp(kx) являются решениями уравнения Бес-

селя

xr dx

Решение, которое при p ^ 0 ограничено в точке г = 0, называется цилиндрической функцией первого рода (или функцией Бесселя):

< п з - 2 )

Так как в уравнение (П3.1) входит р2, то J _ p также является решением этого уравнения. То же относится к любой линейной комбинации Jp и J_p . Цилиндрическая функция второго рода (функция Неймана)1 определяется следующим образом:

.

(113.3)

Sinp7T

V

Часто употребляются также цилиндрические функции третьего рода (функции Ханкеля):

(П3.4)

В качестве общего решения уравнения Бесселя может быть взята линейная комбинация с произвольными коэффициентами любых двух линейно независимых цилиндрических функций. Такими функциями являются, в частности, Jp(x) и J_p (x), если р не равно целому числу. При р = п,

'Эту функцию называют иногда функцией Вебера и обозначают Yp(x).