![](/user_photo/_userpic.png)
книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы
.pdf§ 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)
И*
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, если начальное возмущение является максвелловским
с температурой не ниже, чем температура основной плазмы. Однако для произвольного начального возмущения эта задача еще пол ностью не рещена,
![](/html/65386/283/html_27H_8AcMss.O9nw/htmlconvd-R5p8_T168x1.jpg)
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и |
—т — |
Зи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