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

книги из ГПНТБ / Гульельми А.В. Геомагнитные пульсации и диагностика магнитосферы

.pdf
Скачиваний:
15
Добавлен:
25.10.2023
Размер:
15.82 Mб
Скачать

В плазме низкого давления ( Л Т 4

В2/8п) затухание

гидро­

магнитных волн на частотах со <<J й 4

обусловлено главным

обра­

зом черепковским поглощением энергии волнового поля части­

цами плазмы.

Е с л и в

плоской

волне считать заданным

волновой

вектор

к

= Re к,

то

соответствующее диспе2:>сионное

уравнение

определит комплексную частоту со - j - іу. Декремент

з а т у х а н и я

альвеновской волны в

электронном газе [ И З ]

 

т =

 

'»„ LT

ш 3

..

2 „

ctg2 0).

 

(7.15)

4

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

v\

=

2TJme\

 

 

Te

— температура

электронов.

Предпо­

лагается, что

А <^

ѵтс

с> ig

0 <|J Q{/cù,

co/ßj <С sin 2

0/2 cos Ѳ.

Черепковское поглощение альвеновской волны ионами холодной

плазмы

{ѵт{

<^SA) экспоненциально

мало. Если и электронный газ

холодный в том смысле, что А

ѵу , то затухание и на

электро­

нах также экспоненциально

мало. Тепловых

поправок

к

фазовой

скорости альвеновских волн не возникает.

 

 

 

 

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

в плазме

не очень высокого

давления

имеет

вид

 

 

 

 

 

АѢ*

+

fcj_i4. ( 1

 

 

Аѵтк-±к

X

 

 

 

со-

2 /\

 

к

 

 

 

X

 

 

 

 

 

 

AW-

 

T

m

 

1 +

ft-

+

i 7 f - ) ' ^ " e

x P

ѵ~г_к'„

2

77

m.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.16)

Это выражение справедливо при А

< ^ ѵТе\

cos Ѳ|.

При

обратном

неравенстве

имеем

 

 

 

 

 

 

 

 

со2 == А%къ

+

ѵт.к\

(1 +

Y~

 

 

 

 

 

 

 

- UTTT

 

 

A4"-

\

 

 

(7.17)

 

A v T i k

 

E X P :

 

 

 

 

l

 

 

 

 

 

 

Затухание воли описывается третьим членом в правых частях

уравнений (7.16) и (7.17). Фазова я скорость определяется

выра­

жением

 

 

 

 

 

4 = 4 3 + S 2 s i r ï 2 0 ,

 

 

 

где S2

=

(Т,

+ 2ТЛ/ГПІ при А

<< ѵТе cos Ѳ и 5 2 = 2

( Г е +

Т^щ

при А

^>

ѵТе

cosO. (Напомним,

что рассматривается

случай

плаз ­

мы низкого давления A J§> і>тѵ)

Медленный магнитный звук (акустические волны) сильно за­ тухает (у ~ со), если ТЕ ^ ТІ. Затухание обусловлено черенков-

5L

ским поглощением волн

ионами.

Если же

плазма

неизотермична

( ^ е

 

то затухание

мало

и

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

из

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

уравнения

 

 

 

 

 

 

 

 

 

 

 

m

I т

 

 

 

 

 

 

1

~ ѵ

\ г г )

е х

' М - Т

Декремент

затухания при Те

 

Т\ равен

 

 

 

Г =

] /

g ^ & S f c o s G I ,

 

 

 

 

 

(7.18)

где S'2

=

TJnii.

 

 

 

 

 

 

Циклотронное затухание в плазме низкого давлении велико лишь д л я волн альвеновского типа на частотах со ^ Q,h Распрост­ ранение становится невозможным в узкой полосе вблизи частоты ионно-циклотропного резонанса

 

 

о.

 

 

 

7.19)

i

_

- Z i

 

 

 

 

 

 

 

 

 

Здесь

у

—•

со •— к (А*ѵтУ'>.

 

 

 

В

плазме

высокого давления

г 7'( - > • 5'2 /8

л.)

альвеповскпе

волны

остаются слабозатухающими

даже п])п Те

^

Г,-. Поглоще­

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

мент (у

со) — (со Q,)'2 < ^ 1 мал. Лпгаь

при

/1

ѵт. (ш/й/)*»

погло­

щение альвеновских волн становится

аномально

большим (у-— со).

Что

касается

быстрой

и медленной

магпитозвуковых

волн,

то при углах 0, не близких

к 0 и я / 2 в изотермической

плазме вы­

сокого давления они распространяться не могут. Только в

сильно

непзотермической

плазме

с горячими

электронами и

холодными

ионами возможно распространение слабозатухающих волн этого типа. Довольно громоздкие формулы для декрементов в этом слу­

чае

приведены в монографии 11141.

§

8. Неустойчивость плазмы

Задачу об устойчивости к бесконечно малым возмущениям можно рассматривать как частный случай задачи о спектре собст­ венных колебаний плазмы. В однородной безграничной плазме спектр собственных частот находится путем исследования корней дисперсионного уравнения

Det {ft*:-ар -

-

~~ вдр (to, fc)} - U.

(8.1)

При

заданном

и вещественном fc частоты

со (А"), определяемые из

(8.1),

будут в

общем случае комплексными. Плазма неустойчива,

если

хотя бы

одна из

собственных частот

имеет положительную

мнимую

часть у =

I m

со ^> 0.

Нарастающие волны черпают энер­

гию

из

кинетической

энергии

частиц плазмы.

 

 

 

Тензор

диэлектрической проницаемости. Б у д у т волны

нарас­

тать

или

затухать

— полностью

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

видом

тензора

еа р( со, />"),который, в свою очередь, зависит от исходного

распределе­

ния заряженных частиц по скоростям / ( « ) . Конкретное

выражение

для е а р может быть получено путем интегрирования

линеаризован ­

ного

кинетического

уравнения с самосогласованным

полем [113];

 

 

 

 

 

 

 

 

 

 

 

(8.2)

где

 

 

 

 

 

 

 

 

 

 

 

G a ß [/ (о)]

^-{AsL

 

±Bl^I}f{v),

 

 

 

 

L

 

ЫѴ I

т:

 

д

 

д

 

 

 

 

со

- dvj_

arz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= av.

 

 

д

vx

a

 

 

 

M

\ ^ - +

і'т ^ —

-4-

 

 

 

S. = S . ( » - ^ - s Q ) .

Непулевые элементы матриц Л'ар и 5 „ з равны:

 

 

Ахх =

 

 

 

 

л;.

 

л 2 ,

 

 

 

 

 

 

 

B\Z

=

л ,

где

/ ,

(я.) — функция

Бесселя, a =k±Vj_/Q.

Ф у н к ц и я £(со)

 

опре­

делена

следующим

образом:

 

 

 

 

 

 

го

 

 

 

 

 

 

 

 

£ ( о) = i jj е ш

df

=

— іяб

( о),

 

 

 

где

Р

означает взятие

интеграла

в смысле главлого значения.

Ф у н к ц и я распределения

/„ частиц

сорта а

зависит только от

про­

дольной v. и поперечной ѵ± компонент скорости ѵ и нормирована так, что

\fadV 1,

где dv — 2 лѵ і dv \ dvz.

Формула (8.2) дает наиболее общее выражение для тензора диэлектрической проницаемости однородной иерелятивистской

плазмы.

Неустойчивость холодной плазмы с малой примесью энергич­ ных частиц. Выбирая конкретные выражения для / (v±, vz) и под-

53

с т а в л яя (8.2) в (8.1), получим трансцендентное уравнение для отыскания со при заданном к. В общем случае строгое решение со­ пряжено с большими вычислительными трудностями. Поэтому целесообразно воспользоваться приближенным методом расчета.

Матрицу

в

фигурных

скобках (8.1) обозначим через

Мар.

 

=

^2<5аз — kjbß

— 8 a ß .

(8.3)

Выделим

в

е а р

эрмитову и аитиэрмитову части:

 

Saß

=

ß a ß

+

i&aß,

 

 

где

 

 

 

 

 

 

Saß

=

(«aß

+ 6 р а ) ,

 

 

fc'aß =

~2f

( e a ß

— E ß a ) .

 

 

Будем считать волновой вектор вещественным, а частоту комплекс­ ной (со —> со + іу). В предположении, что у <СС со и e"aß <С ea ß про­ изведем разложение:

de g

Baß (Ш + ІТ) ~ 6aß (Ш) + J S a ß (eu) +

Іу

Подставляя это в (8.3), получим

Л / а 3 = М а р + і М а р ,

где

Л / а Р Ä ft2Oaß ÄaA:9 — -^- 8 a ß ,

(8/1)

,1 г

СО- "

Т

2

e.

c / U J

"

f c

 

ш

 

 

 

Л/aß Ä

с- - В .* р

с 2

<Эсо

Чтобы найти ы(/с), полностью пренебрежем антиэрмитовой частью eœ p. Собственные частоты при этом вещественны, а урав ­ нения Максвелла приобретают вид

Л/а р(со,/с) Eß = 0.

Условие существования нетривиальных решений определяет coa(fc)

Det {Maß 1/с,со° (&)]} = 0,

(8.5)

где a — тип нормальной волны. Можно ввести векторы поляри ­ зации а", удовлетворяющие условию ортоиормировки

а* (oaß

a? = ô a v

(8.6)

54

и соотношению

 

кГОА,> =

daEaßdl,

(8.7)

в котором

со = соа (А;). Поле Е

представляется в виде суммы нор­

мальных

волн

 

Е = 2

Еааа,

 

а

 

 

где Е" — амплитуда нормальной волны типа а (см. [113]).

Д л я отыскания у

(к)

используем

следующее

приближение

по

еа р. Умножим векторное

уравнение

Mapftîj

=

0 на

а"а

В силу

эрмитовости

матриц

М а Р

и М а р имеем паа

Ма$а\

= 0,

откуда

со

f ôw/i I

- 1 "

 

 

 

 

 

,п

O N

T = — — {» - ^ r j

е а 5 я а я з -

 

 

 

 

( 8 - 8 )

Значок a опущен. При написании (8.8) использованы формулы

(8.4)

и

(8.7).

 

 

 

 

 

 

 

 

 

Выражения (8.5)—(8.8) решают вопрос о спектре

колебаний

плазмы, но они применимы лишь при у

со и еа ( з.<^; еа р. Эти

не­

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

основную

массу плазмы

составляют

«холодные» частицы с концентрацией

N

и имеется небольшая примесь «горячих» частиц с

концентрацией

N'.

В силу

аддитивности

тензора диэлектрической

проницаемости

e a ß

=

Bap''1

+

3

Вдали

от полюсов

показателя

преломления

можно

пренебречь

антиэрмитовой частью

е а р л . В

силу

же

нера­

венства N' <^ N можно

пренебречь эрмитовой частью

е а р р .

Та­

ким образом, показатель преломления, поляризация,

групповая

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

Энергич­

ными частицами в данном случае

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

лишь

иикре^

мент.

 

 

 

 

 

 

 

 

 

 

 

Выразим в явном виде инкремент через показатель преломле­ ния холодной плазмы (6.21) и функцию распределения горячих частиц:

сс

 

ï = 2 V,

(8-9)

где

Ф* = - sec2 Ѳ

+ К) е - +

' Si I* Ѳ

[ ^ - К ) ^

F = 4?I-Ö— л ,

9CÜ

55

= - 2 л 2 -± \ dv±dv.ô (со - кхѵ: - sQ) E^L (/),

*

i i =

[ - f Л

( « ) f ,

^ 2

=

[Л (fl)l2 ,

Яl s

=

-

Л

(а)] ІЛ (а)].

 

Расчет инкрементов. Пусть функция распределения

энергичных

частиц

имеет

вид

[1341

 

 

 

 

 

 

 

 

 

 

 

 

.

 

(«• х /и - і ) г "

 

 

 

ехр

( " z - ' V

 

(8.10)

 

/,,( 1 ' - Ы' г ) =

-

Г

^

Г

е х р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"'я

 

 

 

где q =

0, 1, 2

 

 

 

Функция (8.10) довольно гибко описывает раз­

личные ситуации. Так . при «7 =

0, »|| = 0

имеем

двухтемператур-

ное максвеллово

распределение;

при q = 0. =

ы'ц и Н| =/= 0' —

продольпып пучок

частиц;

при q >

1 распределение

имеет мак­

симум при

-

 

 

= q'nv^c

относительной полушириной

поряд­

ка

< Ï'J_>'

 

9/,г/; наконец,

при г/—> ос,

u>j_, ц —> 0 и ш_|_ )/ г/ =

п±

const имеем ô-образігое

распределение/(wj_, уц) =

(2 лѵ±)

— «_!_) ô ( ^

 

— нц).

 

 

 

 

 

 

 

 

 

 

Рассмотрим

случай

продольного

распространения

(0 =

0)

поперечных волн. В к л а д

в (8.2) дают лишь

члены

с s

— +

1. Пе-

исчезающими

компонентами

safi

являются

гѵх

и

е ѵ | / .

Подставляя

(8.10) в

(8.9),

получим

 

 

 

 

 

 

 

 

 

 

 

Т ^

- Y я

 

со

ды->ѵ ~1

Г

со - f 9-а — кги !

 

X

 

 

 

2

 

\kt\w

dm-

ехр !

 

 

 

 

 

 

 

 

 

І,2и

„ \

w

 

 

 

 

 

 

 

 

 

 

 

X I 11

 

\

ici

 

 

 

 

 

 

 

 

(8.11)

 

со

/

!t"T,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К а к обычно, верхний и нижний знаки относятся к Х- и J# -волнам соответственно. Показатель преломления п = с/с/со в (8. Л ) опре­ деляется параметрами холодной плазмы.

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

56

в о л ны распространяются почти перпендикулярно к внешнему маг­ нитному полю [135].

В результате черепковского взаимодействия с резонансными частицами (протонами) магнитозвуковые волны нарастают с ин­ крементом

T = nQ-p w

^

ô ( о -

M , ) ,

(ь.12)

где

 

 

 

 

 

 

•о

 

 

 

 

Д л я функции

распределения вида / оо б (гг — н2 ) имеем

[135

T == -

пОр

/ V

- ^ у.— J l ( н) / 0

(я),

(8.13)

 

 

N

I fr,I »

 

 

где

 

 

 

 

 

ßp

I

Л-2

 

 

Неустойчивость

возникает при

2,4 Q p / » .

 

Характер неустойчивости. До сих пор речь шла об устойчивос­ ти плазмы по отношению к возмущениям, имеющим вид плоской монохроматической волны. Реальные же возмущения так или ина­ че локализованы в пространстве. Положим, в начальный момент t = 0 имеется возмущенно

E(x,0)=^Ekë^dk,

причем спектральная плотность Ек имеет максимум при к = к0. Полуширина спектральной плотности Ак связана с размерами области локализации возмущения Ах соотношением AJiАх — 1. В последующие моменты (1^> 0) возмущение представляет собой суперпозицию нормальных воли

Е (х, 1-) = \ Е><е х Р

+ 1 ( к х ~ с о " г ) ] f № <

( 8 - 1 4 )

где комплексная частота сод- + іу/%- связана с волновым вектором к дисперсионным уравнением (8.1). Образуем квадратичную комби­ нацию величии (8.14) и проинтегрируем ее по всему пространству. Получим

I Е I 2 dx = ^ I Е„ |2 e2 Y *' dk.

(8.15)

Очевидно, что при yh ^> 0 энергия поля экспоненциально

нараста­

ет со временем.

 

57

Выражение (8.15) ничего не говорит относительно области ло­ кализации растущего возмущения. Здесь возможны две ситуации [114, 116]: а) при t~> энергия поля неограниченно нарастает во всем пространстве (абсолютная неустойчивость); б) в любой фик­ сированной точке пространства возмущение все время остается ограниченным (конвективная неустойчивость). При конвективной неустойчивости в однородной безграничной среде возмущение экс­ поненциально нарастает лишь в системе отсчета, движущейся с групповой скоростью пакета. Если же область, где имеет место конвективная неустойчивость, ограничена и если волны могут проникать через границу этой области без заметного отражения, то истинной неустойчивости не возникает вообще, т. е. полная энер­ гия ноля останется конечной. Подобная область плазмы может служить усилителем, но не генератором воли. Естественно, что это обстоятельство придает особую актуальность вопросу о характере неустойчивости плазмы в магнитосфере.

Приближенный критерий абсолютной неустойчивости в случае

одномерных волновых

пакетов

дан в [116. 1231. Пусть

инкремент

у (к)

имеет

максимум

при А- =

кт. Раскладываем у (Л) и со (к) в

ряд по степеням к — кт

с точностью

до квадратичных

членов:

со (к) =

м т + Ѵгр (к — к,„) -і- а (к —

кт)2,

 

 

 

 

 

 

 

 

 

 

 

 

(8.16)

T(Ä)

=

T m - ß ( f t - f c m ) a .

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

i n

=

а(кт),

Тт = ï(A-'m),

 

 

 

а

=

д*л/Ш*

],=,,„,

 

ß =

-

3l T /25fca | , = , m .

 

Так как п р и больших

t эволюция пакета слабо зависит от фор­

мы начального

возмущения,

то примем для простоты,

что

Ек

 

= £ 0 е х р

[ - 4~Zo(fe -

fco)2] •

 

(8.17)

Используя выражепия (8.16), (8.17), преобразуем (8.14) к виду

E(z,t)

=

Ev

(ß + ta) t + ^ 4

X

 

(8.18)

 

 

 

J

 

 

[

 

 

 

4 (ß + ta) t +

4

]

Отсюда видно, что при t —> эо амплитуда неограниченно

нарастает

в данной

точке

пространства,

если

 

 

^ < »

!

= 4

ï m ^ .

 

 

(8.19)

58

Е с ли

выполняется это условие, то неустойчивость абсолютная;

при

обратном неравенстве — конвективная . Обобщение критерия

(8.19) па случай трехмерных волновых пакетов дано в работе [136.. § 9. Нелинейное взаимодействие гидромагиитпых волн

В неустойчивой плазме экспоненциальный рост первоначально малых возмущений не может продолжаться долго. За время поряд ­ ка •— ІІу с момента перехода к неустойчивому состоянию включа­ ются нелинейные процессы, ограничивающие амплитуду колеба­ ний.

Роль нелинейных процессов в динамике магнитосферы весьма велика. Укажем например, что питч-угловая диффузия приводит к самоограничению потока захваченных частиц во внешнем радиа­ ционном поясе [14, 16, 138, 146]. В последнее время, однако, ста­ новится все более ясным, что характер нелинейных процессов в магнитосфере не так прост, как это представлялось ранее: для узкополосиых колебаний в диапазоне Pel (жемчужины) более важным, чем питч-угловая диффузия, по-видимому, является за­ хват частиц в потенциальные я м ы волны [11, 141, 142]; д л я широ­ кополосных колебаний (гидромагнитные шипения, Ipdp) при опре­ деленных условиях основную роль играют распадные процессы [143, 1441; наконец, вполне очевидно, что перенос резонансных частиц вследствие дрейфа поперек узких волновых пучков должен приводить к специфическим нелинейным эффектам, не возникаю­ щим в одпородиой плазме (§ 14, 15).

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

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

вости (со) ^> 0)

компенсируется утечкой энергии в

области по­

глощения

(со) <

0).

 

В линейном приближении различные волны распространяются

независимо

друг от друга (принцип суперпозиции).

Фактически,

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

В квадратичном по амплитуде приближении нелинейное взаимо­ действие воли сводится к распаду одной волны на две и к слия-

59

ниіо двух волн в одну. Резонансные условия взаимодействия имеют вид [1471

 

ci) шI -|- и.,,

к --- к\ -\- к.,,

 

 

где

оз, А' частота и волновой

вектор исходной волны;

co] j 2 ,

Л!,) 2

частоты и

волновые векторы продуктов распада.

 

 

Следует различать две предельные ситуации: а) фазы взаимо­

действующих волн

хаотичны; б)

фазы взаимодействующих

волн

фиксированы. Характерное время сбоя фазы порядка обратной ширины спектра пакета воли Доз. Поэтому, если t — характерное

время нелинейного процесса, то при ^Дсо

1 фазы волн в ходе

взаимодействия не успевают заметно измениться и их можно счи­

тать фиксированными. Если же ^Дю !îs> 1,

то следует пользоваться

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

Рассмотрим взаимодействие гидромагнитпых волн, имеющих случайные фазы. Будем при этом использовать общий метод расче­ та слабонелнпеппых эффектов, изложенный в монографии [120] (см. также [116, 1241).

Исходными будут служить уравнения магнитной гидродинами­ ки (5.22)—(5.28). Считаем невозмущенную плазму неподвижной, стационарной и однородной. Разложим и. (.«. /) в интеграл Фурье:

tt (.«, О =--=

^/^^""-'""ашак

 

 

 

 

 

 

 

п выразим Фурье-компоненты возмущении р,

В и т. д. через

пШк

с помощью (5.22), (5.24)—(5.28). В линейном

приближении произ­

вольное

поле ііи,к

 

можно

представить

и

виде суперпозиции

нор­

мальных

волн:

 

 

 

 

 

 

 

 

 

 

ііи>к

= Sа'і<хРі<8 (м — i»l

)•

 

 

 

 

 

 

(9.2)

Здесь

 

(à'It u n i собственные частоты и нормированные

собствен­

ные векторы (векторы поляризации) волны типа ѵ;

— амплитуда

волны

типа v.

 

 

 

 

 

 

 

 

 

 

Будем предполагать, что в линейном приближении фазы волн

иекоррелированы

и введем спектральную

плотность

 

 

•<ч^|ѵ> = | ^ l 2

ô ( A - -

к').

 

 

 

 

 

 

 

Угловые

скобки

здесь означают

усреднение

по фазам.

Величина

j \|)fc j 1

имеет смысл спектральной плотности

энергии

колебаний

типа

V

на единицу массы плазмы.

 

 

 

 

 

 

 

 

 

 

 

 

 

Нелинейное

взаимодействие

приводит

к

слабой

корреляции

между

волнами,

в результате чего |

| будет медленно

нарастать

60

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