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

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

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

или убывать. Удобно ввести плотность числа «квантов» типа ѵ по­ средством соотношения

M

(9.3)

Кинетическое уравнение, описывающее изменение УѴ^, можно по­ строить феноменологически путем расчета баланса прямых (—>) и обратных (<•—) процессов 1120]. Стрелки обозначают соответствен­ но распад одной волны на две и слияние двух волн в одну.

 

Рассмотрим,

6

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

fci

н а

полны л

V

с волновыми в е к т о р а м и / ^ и к:

ц.

> ѵлПусть

 

и

 

 

вероятность этого процесса, в единицу времени отнесенная к ин­

тервалу dk^kJÇln) ,

 

равна

Р* * (к,

 

7«^,/.',.,).

Полная

вероятность

пропорциональна Nk

для поглощаемых и Ni.-

-+- 1 для

излучаемых

волн. (Подробности

см. в

[120].)

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

случае р.

поглощается,

a ѵ и

т) излучаются,

так что полная вероятность в

интервале dkldkj(2nf

 

равна

 

 

 

 

 

 

P^(Nl\-i)(Nl

 

+

l)Nl.

 

 

 

 

 

 

(9.4)

В обратном

процессе

ц. •-— vi] волны

ѵ и і] поглощаются,

что

дает

P^NlNKNl^-

і).

 

 

 

 

 

 

(9.5)

Согласно принципу

детального равновесия вероятности

P'f'

пря ­

мого и обратного

процессов

одинаковы. Чтобы найти изменение

Ni- в единицу времени в результате процессов

р ~± ѵп следует вы­

честь

(9.5) из (9.4), умножить на с№ч іі/ѵ.2 ;'(2л)е

и проинтегрировать

по

Л",

п Л\2 . Получим

 

 

 

 

 

 

а #

- ) ^ г -

{ММ

-

л"» (.vî, - A t , ) } .

(9-6)

Здесь опущен член, соответствующий спонтанному излучению

волн

V

волнами (.1.

 

 

 

 

 

 

Аналогичным

образом можно

рассмотреть

процессы ѵ ri

pn,

в

которых также появляются

и исчезают волны типа ѵ. Наконец,

в

(9.0) следует произвести суммирование по возможным типам волн

и 11-

Пространственная неоднородность усредненного по фазам волнового поля учитывается заменой в (9.6) частной производной

dldt

на didt =

д/dt -|- -wi'p д/дх,

где « г р =

даі/дк.

Если волны сла­

бо затухают или нарастают в результате

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

с части­

цами плазмы, то в левой части уравнения (9.6)

нужно

добавить

член

y'kNh,

где у'к — инкремент нарастания

колебаний. (Сле­

дует

заметить,

что инкремент

нарастания

числа

квантов

Л"/,., т. е.

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

Gl

Сказанного достаточно для того, чтобы попять, как строится кинетическое уравнение для волн. До сих пор. однако, оставалась неопределенной вероятность Р^п- У к а ж е м без вывода, что в слу­ чае гидромагнитных волн она равна

РГ

(Л, ки

/ , у == (2л ) 7 б

- fcj

- fcä)6

( ol - со1,:, - u>2,)

X

X

-^ѴтГ

I ° ^ / Ѵ « 3 (Ä *. * і ,

**.

. <"*., »2.) a & Ä

I2.

(9.7)

Здесь

Л^,-а

определяется следующим

образом.

Если

произвести

Фурье-траисформацшо исходной системы уравнений (5.22)—(5.28),

то с точностью

до

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

членов

будем

иметь

AijUj-cfc =\NiaßiuiUßdK.

 

 

 

 

 

 

 

 

(9.8)

Здесь Л £ j (ш, А') =

Lu

(/>•) — co2ô, j

;

линейный

оператор L,j опре­

деляется формулой

(7.7); di=

S (со —

о), — oj2) ô

(Л; — кл

— Л; 2 ) X

X dtö1d'.o.,dArldA,2-,

 

« t

=

 

w2 =

 

и ^ .

Нелинейный

оператор

i V i c c p ( « ,

л- ,,

AM,

0 ,

(ÜJ.

ш2 ) имеет вид

 

 

 

 

Л ^ 3

=

-

^

 

№*»

+

S'Fw

+ A*Gw),

 

 

 

(9.9)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

іар =

Ач(А.+

 

Tfci»)A-.,3(cü,.'tu),

 

 

 

 

 

 

 

 

— П а т

П ( З п

{rïmnfejife^fcg;

-f- Пі„(/і1 т /і"2 . — k2mkl:)

2.

-\-

+

( C Ü 2 / M ) U - ] - f

Uimknkzku

(W2 /CÖ) +

г^2пкки}

 

 

— ПірПатЗзАіг

(Мо /м) (&l m fcj, — &гПт р).

 

 

 

 

Здесь

П а р

=

ô a

P

 

— z a z p ; внешнее

магнитное

поле

направлено

вдоль

оси

2.

 

 

 

 

 

 

 

 

 

 

 

 

 

Е с ли учесть анизотропию давления плазмы, то исходное урав ­ нение в квадратичном приблюкепии будет иметь вид (9.8). Линей ­ ный оператор ЬХ] определяется формулой (7.13). Выражение для Niaß. имеющее весьма громоздкий вид, здесь не приводится (см. [143]).

Значительное упрощение возникает при рассмотрении строго

продольного распространения

взаимодействующих воли ± 0)

в плазме низкого давления (S

<gÇ А). В этом случае уравнения

для Фурье-компонент волігового поля в квадратичном приближе ­

нии имеют

вид

 

 

(SW-a*)ux

= - a y U l [ ^ } ( u 1 ± u a ± ) ,

(9.10)

( АѢ^ -

со2) и ±

= ЛаІ M k l U u u , x .

(9.11)

Здесь к

= lcz. При выводе (9.10) мы пренебрегли членами,

описы­

вающими взаимодействия с участием трех

продольных

волн.

(При к±

сколь угодно малом, но не равном

нулю,

эти процессы

запрещены в силу законов сохранения (9.1).)

П р а в

а я часть у р а в ­

нения (9.10) описывает колебания магнитного давления в поле двух поперечных волн (член типа V 2/8л) в уравнении движения). Правая часть (9.11) описывает колебания плотности в продольной

волне и скорости в поперечной (член типа pda/dt

в уравнении дви­

жения) .

= 0 и S <^ А нелинейный оператор

ІѴ,-а р определяется

При

формулой

 

 

Лriaß = -

l h + -g- fc2) zAß

- ~ AlZ A* -

MA*} .

 

 

 

(9.12)

Распадные

процессы и связанная

с ними нелинейная релакса­

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

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

случае к изменению как питч-углового, так и энергетического

спектра частиц. Если, однако,

ш < ^

то можно пренебречь

из ­

менением энергии по сравнению

с изменением питч-углов. В

этом

случае эволюция функции распределения резонансных частиц

описывается

 

уравнением [138]

 

£

L e

J _

a

J

L

W ,

i a | ü ,

(9,13)

dt

 

sin

 

ôa

\

dec

 

где коэффициент

питч-угловой диффузии

равен

D

=

 

- 2

ь2

.

 

(9.14)

— 1

 

Здесь а = arctg (ѵ±/ѵ^) — питч-угол частицы; bl — спектральная плотность гидромагнитных шумов; ѵ\\ Qp/k. При Ак к коэф­ фициент диффузии по порядку величины равен

I

плотность энергии гпдромапштішх

шумов

] Q

\

плотность энергии геомагнитного

поля

J " р '

63

Исследованию питч-углового рассеяния частиц в .магнитосфере посвящены работы [137—139, 145, 1461. Здесь мы отметим лишь следующее. Рассеяние приводит к гибели частиц в плотной атмо­ сфере в результате попадания в конус потерь. При так называемой

слабой диффузии (}/"D* Тц

а0 ) среднее

время жизни частицы

[138]

 

 

 

 

 

О

\ П а

I

 

 

 

где D* =

(v\\lv)

D\

тц ~ Lrjv\\

— время пробега частицы от вер­

шины силовой линии до Земли;

а 0 — раствор конуса потерь. При

сильной

диффузии

((Z)*t||)'« 2§> а 0 ) время

жизни

Нелинейное взаимодействие резонансных частиц и волн, имею­ щих достаточно узкий спектр, приводит к тому, что появляется группа частиц, запертых в потенциальных ямах волнового поля. Число захваченных частиц пропорционально корню-квадратному из амплитуды, т. е. может быть большим даже при относительно невысокой амплитуде волны [1211.

Захваченные частицы совершают колебания около положения равновесия Ц41, 142, 1481. Сделаем элементарную оценку частоты колебаний резонансного прототта в поле монохроматической альве-

новскон волны. Д л я

этого перейдем в систему отсчета, движущую ­

ся с фазовой скоростью волны. Относительная

фаза волны и части­

цы равна

ср = kz{t)

Qpt

- j -

const. Дважды

продифференцируем

Ф (t)

по времени: <| — кѵг.

В выбранной системе отсчета

электри­

ческое поле волны отсутствует,

так что па частицу действует лишь

сила

Лоренца (е/с) [vb\.

С учетом

этого получаем

ѵ, —

(е/трс)

X

X ѵ^Ь^ sin ср.

 

 

 

 

 

 

 

 

 

 

 

Такилі

образом,

 

ср .—• kvj_Qp

sin

ср,

где

Ц , =

eb±/mpc.

Частота

колебаний,

очевидно,

равна

со — (kvjQ')1-''.

Отсюда

при

ѵ± — vz

с

учетом резонансного

условия kvz

~~ Qp

имеем

 

 

с о ~ К Б ^ р -

 

 

 

 

 

 

 

 

 

(9.15)

Захваченными оказываются частицы, продольная скорость кото­ рых отличается от резонансной скорости vz — Q,,/A: на величину àvz порядка

Avz. ^ -, /И'~

Если амплитуда волны экспоненциально нарастает, то увели ­ чивается (также экспоненциально) и число захваченных частиц. Это приводит к осцилляциям инкремента с частотой порядка со.

64

Предполагается, что за время нескольких осцилляции инкремент может существенно уменьшиться 9 .

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

но

узкий спектр, то

t со ІІЬ в

случае

распадпых процессов [147]

и

і м I/ft1.'* в случае

захвата

частиц

волной.

0 Более подробно эффект захвата частпц волпой проаналпзпрован в работах [121, 141, 142, 148-150].

А. В. Гульельми, В. А. Троицкая

Г л а в а I I I

КОЛЕБАНИЯ II ВОЛНЫ В МАГНИТОСФЕРЕ

§10. Волновое уравнение

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

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

Такой подход впервые использован в работе [151]. Были выве­ дены два одномерных волновых уравнения, описывающих коле­ бания альвеновского типа («тороидальный» и «гофрированный» моды) и двумерное уравнение для колебаний магиитозвукового типа («полоидальный» мод). Большое число работ посвящено ис­ следованию этих уравнений и интерпретации на их основе неко­ торых типов пульсаций [152—173].

Данный и следующий параграфы посвящены собственным ко­ лебаниям магнитосферы 1 . Основное внимание уделено анализу спектра альвеновских колебаний. Альвеиовские колебания инте­ ресны тем, что с ними .можно отождествить локализованные по широте пульсации Рс4, 5, имеющие высокую корреляцию в сопря­ женных точках. Д л я нас альвеиовские колебания интересны еще и потому, что они являются полезным диагностическим средством.

Прежде всего запишем исходные уравнения в форме, удобной для анализа. Введем систему криволинейных коордииат ( I 1 , £2 , £3 )

1Мы оставим без рассмотрения интересную проблему колебаний геомагнит­ ного хвоста. Ей посвящены работы [174—180].

66

с метрическим тензором giH. Квадрат

элемента длины будет равен

dl* =

gihdlldlk.

 

К а к

известно, при использовании

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

следует делать различие между контравариаптными (индексы свер­ ху), ковариантными (индексы снизу) и физическими (индексы в

скобках) компонентами векторов и тензоров. Связь

между ними

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

метрическим тензором.

 

Уравнения электромагнитных колебаний малой амплитуды в

криволинейных координатах имеют вид

 

(rot Е)а =

~ Ьа,

(10.1)

(rot Ь)а = -

eapii0 .

(10.2)

Здесь принято,что среда стационарна, Е ос ехр (— iat), компонен­ ты тензора диэлектрической проницаемости е а р (со) определяются в каждой точке пространства локальными значениями параметров

среды. В

общем случае

тензор е а р (со) неэрмитов.

 

Используя

тензорную

плотность

о'1"'"', можно

записать

(rot Е)л

= g^(vot

Ef = І 2 | &ИЕІЛ.

 

 

(10.3)

Здесь Eu

j =

dEildl?,

g =

Del { g „ p } . Компоненты 6 i t m

равны пу­

лю, если

хотя бы два любых индекса

одинаковы;

б 1 2 3

= + 1, а

остальные иеисчезагощие компоненты меняют знак при перестанов­ ке любых двух индексов [181].

Исключая из (10.2) поле b с помощью (10.1) и используя (10.3),

запишем

волновое

уравнение

в

виде

 

 

 

 

 

 

 

 

=.-£„,.„

= ~ е , , ( с о ) £ Л ' .

 

 

(10.4)

К уравнению

(10.4)

следует

добавить

уравнение

 

 

 

(УІЧіЕЪ

= 0,

 

 

 

 

 

(10.5)

 

которое

вытекает

из

divJD = 0 ,

Ö ; =

гиЕ! и

формулы

div а

=

= (IlYІ)

(Yêai)>i-

Далее следует задать граничные условия на по­

верхности Земли, на нижней

границе

ионосферы, на поверхности

плазмосферы

и на внешней границе магнитосферы. Задача

будет

заключаться в определении поля

Е из уравнений (10.4), (10.5) пр и

соответствующих граничных

условиях . Поле b затем определяется

из уравнения

(10.1).

 

 

 

 

 

 

 

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

Введем

локальную

декартову

систему координат (ж1 , хг,

х3)

с

осью X3 , направленной по касательной к силовой линии геомаг­

нитного

поля . Тогда

физические компоненты

тензора

проницае-

 

 

 

 

 

 

 

 

3*

 

«

мости

e

в

приближении холодной

плазмы

даются

формулой

(6.6).

 

 

 

 

 

 

 

 

 

 

Ковариаитиые компоненты тензора проницаемости выражаются

через

физические

следующим

о б р а з о м 2 :

 

 

 

_

дхі

дѵ';

 

 

 

 

 

 

 

Нетрудно

убедиться, что

в холоднонлазмеином

приближении

e-aß =

/

ВМ

р£J -

. ,

Y7&

,

в*вр

 

n n r s

е{ga ß

i £ Ô a Y 3

 

h il —j^--.

(10.6)

Здесь

В

= У

В*Вa.

 

 

 

 

 

 

Уравнения

(10.4), (10.5) с тензором

проницаемости (10.6) будут

исходными при рассмотрении свободных колебаний магнитосферы.

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

при

которых были написаны эти уравнения .

 

а) Амплитуда колебаний мала настолько, что связь D с Е

мож­

но считать линейной. П р и рассмотрении пульсаций с амплитудой порядка нескольких гамм линейное приближение вполне допусти­

мо. Но его, вероятно, нельзя использовать даже в качестве

первого

шага при анализе пульсаций с амплитудой —

50у, иногда

наблю­

даемых на далеких магнитных оболочках.

 

 

б) Связь D с Е предполагается локальной,

т. е. не учитывается

пространственная дисперсия. Приближение локальной связи при­ менимо, если размер неоднородности волнового поля значительно превышает гирорадиус ионов, а период колебаний намного меньше времени пробега электронов на дистанцию порядка длины волны [114]. В диапазоне пульсаций первое условие легко может быть

нарушено,

особенно на периферии .магнитосферы.

в) Среда

стационарна, а положение и форма

внешней поверх­

ности границы магнитосферы — фиксированы.

Предположение о

стационарности параметров плазмы внутри магнитосферы довольно

естественно. Что же

касается внешней поверхности, то считать

ее фиксированной,

строго говоря, нельзя . Фактически должна

рассматриваться задача о колебаниях объема замагничеииой плаз ­

мы со свободной внешней поверхностью. Если,

однако,

волновое

поле локализовано вдали от внешней границы

(как, например,

в случае альвеновских колебаний достаточно

глубоких

магнит­

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

г) Плазма предполагается холодной. Этим самым исключаются из рассмотрения волны акустического типа, дрейфовые волны, эффекты анизотропии давления плазмы и т. п.

Если система (Ç1, Ç2, £3 )

ортогональна, то е а ( 3 = (gagß)'he.at?,, где ga — диа­

гональные компоненты

метрического тензора.

68

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

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

сравнительно

простых

вычислительных средствах.

 

 

В диапазоне

частот

со

Qe

компоненты диэлектрического тен­

зора без учета

соударений

между частицами имеют

вид

 

- - • " * *

 

- - -

" в

, Ч « - 4 - .

(10.7)

Отсюда

видно,

 

что в гидромагнитном диапазоне ( с о < ^ £2,-) можно

пренебречь педиагональными

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

 

 

газ =

/

-

ВВ„ \

 

B„BR

(10.8)

 

е Ц

- ^ )

+ Л -^г->

где

е ^

CÛQ І /й?, и ~

—cùU e / со2.

 

 

Гирочастота

протонов Qp

=

3 - 10 3 /L 3 уменьшается с удалением

от

Земли. Н а

оболочке

L ^

10 гиропериод протонов

около двух

секунд,

так что аппроксимация

(10.8) является вполне

удовлетво­

рительной при рассмотрении пульсаций с периодами десятки — сотни секунд.

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

Как обычно делается при анализе первых гармоник собственных колебаний магнитосферы, положим Е В = 0, т. е. пренебрежем параллельной геомагнитному полю составляющей электрического

вектора. Основанием для этого

может

служить следующее. Со­

гласно

[114],

продольная составляющая

электрического

вектора

Е у ~

(е/п)

Л и ) 2 Е±,

где к±

и к~1 — характерные размеры не­

однородности волнового

поля поперек и вдоль магнитных

силовых

линий.

В типичных условиях е/т] х co2 /Qe Q; — Ю - 8 ч - Ю - 7 , т. е.

весьма

малая величина.

Предположим теперь, что выбрана такая ортогональная систе­ ма координат, в которой силовые линии геомагнитного поля — прямые, совпадающие с координатными линиями t,1. Тогда уравне-

69

иие (10.4) с учетом (10.8) при условии Е В = 0 примет вид

JL

^

*

*

* а _ £

£ «

) .

(ІО.9>

Здесь a =

2, 3,

ß = 5 —a; gj

— диагональные

компоненты

мет­

рического

тензора.

 

 

 

 

 

 

Система

(10.9) без правой

части

представляет

собой

пару

укороченных уравнений, которые в однородной среде описывают одномерные (альвеновскне) колебания. Возникает вопрос, при каких условиях решение укороченного уравнения в неоднородной среде является хорошим приближением к одному из решений урав­ нения (10.9). Очевидно, достаточным условием будет тождествен­

ное выполнение

равенства

 

 

 

ж

=

ж

-

 

 

 

 

( 1 0 Л 0 )

Условие (10.10) может быть выполнено в аксиально симметрич­

ной магнитосфере. Будем считать, что

£3 ср у г л о в а я

коорди­

ната,

меняющаяся от

0 до 2л, и что

ö#a p/dcp

— 0, ô/J/dcp - 0,

de/dcp = 0.

Согласно

условию (10.10),

линейно

поляризованные

колебания

Е

=

(0, Ег,

0) также должны быть аксиально

симмет­

ричными: дЕ-2 /дср=0. Колебания такого типа называют тороидаль­ ными.

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

поляризацией Е — (0, 0, Е3).

Чтобы при этом выполнялось условие

(10.10), в разложении

поля по азимутальным гармоникам

должны

преобладать большие

азимутальные

числа

 

 

Ей = —-^Щ—*0

при

w - > со.

 

 

(10.11)

Подобные колебания

будем

условно

называть

гофрированными.

Уравнения тороидальных

и гофрированных

колебаний

выведе­

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

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

Поле

Е

в окрестности данной спловой линии удобно рассмат­

ривать

в

координатах

 

?=1,

 

l* = q,

£ 3 = Ѳ ,

70

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