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

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

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

§ 3. Преобразование Фурье по. времени

163

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

Рассмотрим функцию г (к, со) в пределе низких и высо­ ких частот для частиц с максвелловским распределе­ нием Fv.

В пределе низких частот (сог 0) выражение (3.10)

принимает вид

 

 

 

е (к, со = 0) = 1 2

j

К k) [<9*

+ шб (u)] du.

V

 

 

 

 

 

 

(3.13)

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

Fv (u, k) = cve x p [—

.

(3.14)

Подставляя эту функцию в (3.13), окончательно полу­ чаем

 

 

 

 

4 + 2

PV

кЧI

 

 

 

 

 

 

 

 

кЩу

 

 

 

 

где

 

 

 

 

 

 

 

 

(3.15)

 

 

 

 

1

^ 4ппче2

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

Для

фурье-компоненты потенциала мы имеем

 

 

 

Ф ( к

с о __q \

4 я е

гевнешн (fc,

со — 0 )

4 яегавн еш н ( к ,

со — 0 )

16)

'

'

кг

l-M/dn)-2

№ + №

'

 

С помощью теоремы свертки обратное преобразование

выражения (3.16)

дает

 

 

Ф (г, со — 0) — е

^ пвнеШ (г , со — 0)

е х р

Яв

 

 

-г' I

 

 

 

(3-17)

И*

(3.21)
Умножим эти уравнения на exp(iat) и проинтегриру­ ем по t от —оо до +оо. Функция Ai связана с функцией Дирихле следующим образом:
Ф ВНешн-

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

При получении формулы (3.17) было учтено, что

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

В пределе высоких частот сог оо знаменатель в выра­ жении (3.10) можно разложить в ряд

ОО

(3.19)

Затем, подставив функцию распределения (3.14) в (3.10), найдем

(3.20)

Как и следовало ожидать, электрическая восприимчи­ вость, пропорциональная е — 1, из-за наличия инерции частиц в пределе очень высоких частот стремится к нулю как сор/со2.

3.2. Начальная задача

Чтобы решить начальную задачу, рассмотрим опять уравнения (3.1), но без внешнего возмущения При этом

{ 4 t + iku ) Т{и, k, t) = i- ^ -kn[ 9F°^uu' ] Ф (к, t),

 

 

§

3.

Преобразование Фурье

по

времени

165

Согласно (2.39),

производная функции

равна

 

 

 

 

 

Т Г = «(*).

 

 

(3-23)

Отсюда следует,

что

 

 

 

-{-00

^

 

+00

 

 

 

j

 

&1вш

dt =

- J 7(u, к, t) [8 (t) + гсоД,] <*»t dt

(3.24)

— 00

 

 

—00

 

 

 

и

соответственно

 

 

 

- | - o

o

___

 

—|—00

 

 

 

j

 

Aiei(0<4 f

dt=

7(u>к, 0) — гео j

^ е ш ]{и, k, t) dt.

Используя вместо (3.3) разложение1)

 

(3.25)

 

 

 

 

 

 

+ 0 0

 

 

 

 

 

Q (и, '•ш)=тт JQ(u, r, t) A±eiaitdt,

(3.26)

а также выражение (3.25), можно уравнение (3.21) пере­ писать в следующем виде:

г (ки — со) J (и, Л, ю) — i +

Лге [ -- ° ^ ’ ^ ] Ф (к, со) =

^

= w * 1{u' k' t= 0 ) •

+°° _

(3.27)

Ф(/с, (х>)=~-

j f(u, /5,

(o)dw.

— оо

Разделив первое уравнение на i (ки — со) и проинтегри­ ровав по и от —оо до +оо, а также используя второе

уравнение (3.27), мы получим потенциал Ф (к, со), выраженный через заданное его начальное значение. При

этом можно также найти / (и, к, со) из первого уравнения (3.27). Обратное преобразование по пространству и вре­ мени дает искомые функции. Таким образом в принципе задача решена.

!) Знаком •— сверху здесь и в дальнейшем обозначается изобра­ жение Лапласа.

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

Прежде чем перейти к оценкам, сделаем еще одно заме­ чание. Вводя обозначения

ш = ip, рг =

toi, Pi = —con

(3.28)

имеем

oo

 

 

 

£2 {и, г, р) = ^

f e~pt £2 (и, г, t) dt.

(3.29)

~\/2n J

 

Здесь £2 (и, г, р) является теперь изображением Лапла­ са для функции £2 (и, г, <). Конечно, мы должны отме­ тить, что контур интегрирования обратного преобра­ зования в p-плоскости отличается от соответствующего контура в ш-плоскости. Действительно, р = рТ + ipt, где Pi изменяется от —оо до +°°> а рг — постоянная величина, которая больше максимального инкремента нарастания rim.

§ 4. РЕШЕНИЕ УРАВНЕНИЯ ВЛАСОВА МЕТОДОМ ЛАНДАУ

4.1. Решение Ландау Подставляя (3.28) и (3.29) в (3.27) и исключая функ­

цию f(u, к, р), находим выражение для потенциала

+оо

 

_ 1

[f(u,

к, t =

0)/(iku-{-p)]du

*

 

4л« V 2.

+оо

 

(4.1)

Ф{к,р) = ^ е

 

 

 

 

1 — г—g-

l

[F^(u)/(iku-j-p)]du

 

 

 

 

— 00

 

 

 

Отсюда имеем

 

 

 

 

 

 

a+ioo

Г +5°

[7(u, k , t = 0)/(iku +

'i

 

|

£

p)]du j

 

 

 

 

 

-----------------------

)<**•

 

 

| l

-

? j

[P^(u)/(iku +

p)]du

l

(4.2)

§ 4. Решение уравнения Власова методом Ландау

167

Отметим, что контур интегрирования представляет собой вертикальную линию в комплексной p-плоскости с о > цт , т. е. находящуюся справа от

всех

сингулярностей

подынте­

Im(“)

 

грального выражения.

 

 

 

 

Ограничимся теперь рассмо­

 

1 (р/к)

трением функций распределения

 

вида

Fо (и), которые,

согласно

 

 

Re(и)

критерию Пенроуза,

являются

 

 

 

 

 

устойчивыми и поэтому не при­

 

 

 

водят к комплексным корням ир

 

 

 

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

соотношения.

1гт' П

 

Следовательно,

подынтеграль­

 

ное

выражение

(4.2)

представ­

 

 

 

ляет собой регулярную функ­

 

i (р/к)

цию как в полуплоскости при

 

рг >

0, так и при рТ<

0, поэ­

 

-- -\*г

- УRe (и')

тому

величину

а >

0

можно

 

 

 

выбрать произвольным образом.

 

 

 

Непосредственно

вычислить

 

 

 

интеграл (4.2) невозможно, по­

1ш(и)

 

скольку выражение под инте­

 

 

 

гралом имеет разрыв на мни­

 

 

 

мой оси р, что легко видеть,

 

 

 

если

применить

формулу Пле-

 

 

Ко (и}

меля (2.74). Поэтому мы будем

 

 

 

ф 1(р/к)

следовать методу, предложен­

 

ному Ландау и основанному на

 

 

 

аналитическом

продолжении

Фиг.

14.

Деформация

функций. Помимо тех условий,

контура

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

которые были уже наложены на

при аналитическом продол­

Fо (и), потребуем, чтобы функ­

 

жении.

ции

F0 (и) и

/ (и,

к,

t = 0)

 

 

 

имели аналитическое продолжение, а функция Ф (к, р) достаточно быстро стремилась к нулю при | р \ —*- оо. Накладывая эти требования, мы стремимся к тому, чтобы можно было применить теорию вычетов к вычислению интеграла (4.2), в котором путь интегрирования замы­ кается по большой полуокружности в отрицательной полуплоскости. В общем случае этого сделать нельзя

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

из-за наличия разрывов в подынтегральном выраженииОднако значение интеграла вдоль указанного контура интегрирования (от а — ioo до о + too) не зависит от значений подынтегрального выражения вне этого кон­ тура. Следовательно, если можно найти аналитическое продолжение подынтегрального выражения в области, расположенной слева от пути интегрирования, то для вычисления интеграла (4.2) можно применить теорию вычетов.

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

D(k, р) = 1 — ^

2 + ° °

du для

рг > 0, (4.3)

f

 

-» и— i -т-

 

 

D(k,p) = 1

 

du + inF'0

j -)

для рт= О,

 

 

 

 

(4.4)

D (к, р) =

Fq(и)

du + 2niF'0 (-у-)

для Р г< 0 .

и i ук-

 

 

(4.5)

 

 

 

 

Аналогично строится аналитическое продолжение числи­ теля подынтегрального выражения (4.2).

При

таком аналитическом продолжении функции

D (к, р)

ее нули могут лежать только в открытой левой

р-полуплоскости х). Эти нули являются изолированными,

х) Заметим, что при таком аналитическом продолжении функция D (к, р) может иметь нули, хотя предполагалось, что знаменатель в первоначальном виде никаких нулей не имеет,

§ 4. Решение уравнения Власова методом Ландау

169

поскольку существование неизолированных нулей (точек сгущения) означало бы, что D (к, р) = О1).

Применяя теорему вычетов, сдвинем путь интегриро­ вания в обратном преобразовании Лапласа на величину

а = — а и значения интеграла

в полюсах представим

в виде суммы вычетов:

 

 

 

—a-j-ioo

Ф ( М ) = = 2 ^ + ^

J

ф ' (к, p)evtdp. (4.6)

v

— a — гео

Здесь Ф' (к, р) — аналитическое

продолжение функции

Ф (к, р), ас , является вычетом функции Ф' (к, р) в полю­ се pvСуммирование распространяется на все нули/) (к, р). Предполагая, что интеграл в (4.6) стремится к нулю при a - > оо, получаем

Ф {к, t) = 2

cvepv4 ~

Coepot.

(4.7)

V

t-*oo

 

Таким образом, поскольку

p vr <

0, начальное возмуще­

ние затухает. Это затухание называется затуханием

Ландау [5].

сю остается только

Как видно из (4.7), в пределе t

слагаемое, обладающее наименьшей по модулю действи­

тельной частью | рг |. Поэтому среди

корней

уравнения

*) Для максвелловской функции распределения F0(u) не выпол­

няется предположение Ф (к, р) -*■ 0 при | р |

оо. В

этом случае

D (к , р) имеет бесконечное число нулей в левой р-полуплоскости [4], которые сгущаются на бесконечности и располагаются вблизи линий arg р = ± Зя /4. Однако любая вертикальная узкая полоса содержит только конечное число нулей D (к, р) [14]. Отсюда следует, что обратное преобразование определяется следующим образом:

Ф ( М ) = 2

cvePv<+

0(ef><),

Pv r>P

 

где p v — полюса в полосе р <

pvr < 0.

Кроме того, можно пока­

зать, что вклад от пути интегрирования, сдвинутого в сторону беско­ нечности, становится исчезающе малым для произвольного момента времени t > 0, если начальное возмущение является максвелловским

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

170

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

D (к,

р) = 0 найдем корень р 0 с наименьшей действитель­

ной

частью

| pr |. С этой целью

разложим функцию

D (к,

р) в степенной ряд

по мнимым значениям p t.

Это можно

сделать с

помощью

обобщенной формулы

Племеля [12]. Допустим, что функция g (и) удовлетворяет всем требованиям, которые предъявлялись к функ­ ции F0 (и)\ тогда интеграл

' w = +J -Ш-**

<«>

—QO

 

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

оо

/(*) = 2 4 - / ,n,(*o)(*-*o)n-

(4.9)

п = 0

 

Из (4.8) для значения z0, лежащего в верхней полуплоско­ сти, находим

—j— ОО

- j - 0 0

00 — оо

Отсюда производная /г-го порядка может быть записана по индукции в виде

/ <л,ы = \ I^ - du-

(4.11)

— оо

 

На основании формулы Племеля (2.74) для z0 = х0 + iy0 можно получить

+0 0

lim I im(z0) = ^ ‘

f

.g.JAu) . du |_ ingln>(x0).

(4.12)

»<Г+0

Л

U~ X°

 

Отсюда следует, что при х0 = х и у0 = О

ОО

+ 0 0

I(x + iy)= 2 ■ ¥ ! >Su-JUJ-duA

( д ) ] •

 

(4.13)

Tf=0

-09

§ 4. Решение уравнения Власова методом Ландау

171

Положим g (и) =

F'0 (и), х =

p t/k

и у = рТ/к]

тогда

разложение D (к,

р) по мнимым значениям р примет вид

 

9

 

 

 

 

cod

 

 

 

D(k,p) = 1— ^ 2

 

 

 

 

n=0

 

 

X & J

Р(°П+1У du

| - m . F

„ n + 1 ) (j-)~

(4.14)

Поскольку мы ищем нуль функции D (к, р) — 0 с наимень­ шим значением | рг |, можно пренебречь членами с более высокими степенями рт/к. Тогда, приравнивая правую часть (4.14) нулю, после разделения мнимой и действитель­ ной частей получаем

“+ _г

 

 

и

 

 

 

 

o= kf; ( ~ ^ - ) + ^ ^

+(

- Ш - с 1 и + о ( - 1|-). (4.16)

 

 

u + ~

t

 

Дифференцируя (4.15) по к,

находим

 

P i — k

dPi

-{-00

П (и ) du.

 

2к

d k

f

(4.17)

W

 

3

« + %

 

 

Исключая интегралы в (4.16) и (4.17), приходим к сле­ дующему соотношению:

P i _ \

/

к

d p t

(4.18)

к )

\

p i

d k

 

В это соотношение мы должны подставить p t из (4.15). Для простоты ограничимся рассмотрением малых значе­ ний к. Тогда существенный вклад в интеграл дает только область | и | pilk (случай низких температур), поэтому можно разложить зцаменатель в ряд по степеням ик/рj.

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

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

к2

 

 

 

 

 

 

 

ku

k2u2

k3u3

+

j du.

(4.19)

 

 

 

 

 

 

 

 

~Pi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда после интегрирования по частям получаем

 

 

к2

fe

+°°

р , .

 

 

к . п к2 о 2 /с3

 

 

 

 

f

 

)

/

.. j du.

—Г ~

------ \

J

^0

( U

---------- Н

—т —

Зи2

—г +

<*1

Pt

 

 

 

w

 

'

Pi

Pi

 

Pi

 

(4.

0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что (и) = 0,

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<^>

В нулевом приближении мы получаем

решение р\ =

сор,

и в первом приближении

 

 

 

 

 

 

 

 

 

 

Pi «

±

(Op [ l + 3 (и2) -^ -]1/2 «

 

 

 

 

 

 

»

 

±

 

Юр

 

+

у (ц2)

 

+ (Ор.

 

(4.22)

Это решение совпадает с формулой (2.58). В самом деле, приближение (4.15) приводит нас обратно к мнимой оси р (т. е. к действительной оси ир). Более того, пренебрежение сингулярностью сводит уравнение (4.15) к дисперсионно­ му уравнению, полученному для случая кусочно-гладкой функции распределения с провалами. Так же как и ранее, здесь мы ограничились рассмотрением малых значений к. Для того чтобы найти соответствующее выражение для рг, продифференцируем выражение (4.22) по к. Простые вычис­ ления приводят к следующему результату:

1

к dpi Pi dk

2 ( } to2

Jr2

(4.23)

 

1 — 3 (u2) ~ « 1.

^ 2 v CO2

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