Научная статья на тему 'Динамика заряженных частиц в расщепленных тонких токовых слоях'

Динамика заряженных частиц в расщепленных тонких токовых слоях Текст научной статьи по специальности «Физика»

CC BY
47
8
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по физике, автор научной работы — Оводков Д. А., Попов В. Ю., Малова X. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Динамика заряженных частиц в расщепленных тонких токовых слоях»

УДК 533.92

ДИНАМИКА ЗАРЯЖЕННЫХ ЧАСТИЦ В РАСЩЕПЛЕННЫХ ТОНКИХ ТОКОВЫХ СЛОЯХ

Д. А. Оводков, В. Ю. Попов, X. В. Малова

(.кафедра математики)

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

Введение

Динамика заряженных частиц в «колоколообразных» (или «одногорбых») тонких токовых слоях (с толщиной L порядка ионного ларморовского радиуса) является предметом интенсивного изучения на протяжении нескольких последних десятилетий [1-3]. В последние годы благодаря исследованиям спутников GEOTÂIL и CLUSTER были обнаружены двойные токовые слои, профили плотности тока которых имеют максимумы на краях и минимум в центре [4, 5]. Структура и эволюция TTC не могут быть описаны в рамках МГД-теории, так как гирора-диус ионов, носителей тока, порядка толщины слоя. Более адекватной является кинетическая теория, где детальное описание движения разных групп частиц очень важно для структуры и динамики TTC как целого. Как удалось показать ранее [6], процессы хаотического рассеяния частиц в TTC могут играть ключевую роль как в образовании двойных токовых слоев, так и в их разрушении. Однако этот вопрос пока недостаточно хорошо изучен, и требуется более детальное рассмотрение динамики частиц в двойных токовых слоях. В настоящей работе проводится исследование влияния общей структуры токового слоя (одногорбый или двугорбый профиль плотности тока) на структуру фазового пространства, в котором движутся заряженные частицы. С этой целью исследуется зависимость фазовых областей от величины sjLjи от формы профиля плотности тока, одногорбой или двугорбой, где L — полутолщина области обращения поля, ро максимальный ларморовский радиус. Авторы попытались подойти к обратной задаче, т. е. по структуре фазового пространства сделать предположения о структуре и возможной эволюции TTC.

1. Динамика заряженных частиц в TTC и сечения Пуанкаре

В тонких токовых слоях движение частиц можно разделить на три основных типа: пролетные частицы на разомкнутых (или спейсеровских [7]) орбитах, захваченные частицы на круговых орбитах и ква-зизахваченные частицы на так называемых «огур-

цовых» орбитах [8] (рис. 1). Спейсеровские ионы, замагниченные вне слоя, приходят из бесконечности и уходят на бесконечность. В толще слоя спейсеровские ионы размагничиваются, совершая «меанд-ровое», петляющее движение с пересечением нейтральной плоскости. «Огурцовые» орбиты в центре слоя также совершают меандровое вращение, но вне слоя их орбиты почти замкнуты, и такие частицы многократно возвращаются в слой в процессе крупномасштабного вращения. Они могут захватываться в слой из пролетных орбит или в процессе рассеяния уходить из слоя на бесконечность. «Круговые» орбиты полностью замкнуты, они не выходят из слоя и не пересекают сепаратрису, разделяющую два типа движения — с пересечением токового слоя и без пересечения. Основными носителями тока в TTC являются спейсеровские ионы. Частицы на «огурцовых» орбитах не дают вклада в полный ток, так же как и ионы на круговых траекториях. Однако локальный ток этих частиц ненулевой, он может изменять профиль основной плотности тока так, что в центре слоя появляется локальный минимум, т.е. слой становится расщепленным [6]. Последние спут-

никовые исследования демонстрируют, что такие слои могут играть существенную роль в динамике суббурь [9, 10].

Движение пролетных и квазизахваченных частиц в одногорбых TTC носит квазиадиабатический характер. Это означает, что при пересечении сепаратрисы магнитный момент частицы (или, что практически то же самое, приближенный инвариант движения Iz = j^§vzdz) претерпевает скачок, величина которого много меньше величины самого магнитного момента. В работе [3] динамика частиц была изучена при помощи сечений Пуанкаре [11], которые представляют собой двумерное изображение множества точек пересечения частицами плоскости Z = 0 в фазовых координатах. На рис. 2 показано, что фазовое пространство частиц состоит из трех непересекающихся разных областей: область интегрируемого движения (А, круговые орбиты), динамического хаоса (В, квазизахваченные орбиты) и квазирегулярного движения (С, спейсеровские орбиты). В работе [8] было показано, что квазиадиабатическое приближение выполняется, когда параметр bn\/Ljра, описывающий динамику частиц, много меньше единицы.

Рис. 2. Сечение Пуанкаре для «колоколообразного» профиля плотности тока TTC, zc = 0, а = 1.78, Ъ„ = 0.1. Сечение построено в координатах {<тЬ„х,рх}, причем безразмерная координата аЪ„х и импульс рх таковы, что все точки сечения Пуанкаре отображаются на единичный круг

Существует два основных механизма етохаети-зации движения частиц при рассеянии в сильно искривленных магнитных полях. Первый связан со скачками Iz при пересечении сепаратрисы, второй механизм проявляется при перекрытии нелинейных резонансов (возмущенных инвариантных торов движения) в окрестности регулярного движения [12, 13]), что позволяет объяснить хаотическое движение частиц в отсутствие сепаратрисы, например в трилинейной модели [14]. Рассматриваемая в данной работе модель магнитного поля позволяет с единых позиций подходить к изучению динамики частиц в созданных ранее моделях: трилинейной, клиновид-

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

2. Модель

Рассматривалась принятая в геофизике еолнеч-но-магнитосферная система координат, в которой ось X направлена от Земли к Солнцу, Z — перпендикулярно плоскости эклиптики, а ось Y направлена с утренней на вечернюю сторону Земли. Предполагалось, что магнитное поле имеет две компоненты: Вх = BqF(Z) и Bz = Вп. Исследовались две модели TTC — колоколообразная и двугорбая. Профиль плотности тока «колоколообразного» токового слоя имеет вид

Jy=Ldh^fy (1)

где Jq = ^, L — полутолщина области обращения магнитного поля, Bq — магнитное поле на бесконечности, с — скорость света.

Постоянное электрическое поле, всегда присутствующее в токовом слое магнитосферного хвоста, не принималось во внимание, так как движение частиц рассматривалось в системе координат де Хофф-манна-Теллера.

В настоящей работе модель двойного (или расщепленного) слоя представляет собой суперпозицию двух «одногорбых TTC» (рис. 3). Расстояние между горбами контролируется параметром Zc, а ширина горбов — параметром L:

Можно показать, что при таком задании магнитного поля полный ток в системе будет сохраняться при изменении параметров Zc я L.

Jy/Jf) 0.5

0

Рис. 3. Профиль плотности тока для расщепленного TTC

В предельных случаях при Zc ^ 0 получается гиперболическая модель магнитного поля, при L ^ 0 — трилинейная модель, а при Zc, Ь —> 0 — клиновидная модель магнитного поля, изученные ранее в [8, 14, 15].

3. Анализ движения частиц в TTC

Гамильтониан движения в исследуемом поле ра-

вен

1

PÎ) + *=[P,--Mx,Z)]2. (з)

Поскольку гамильтониан системы не зависит явно от У, то канонический момент импульса Ру является интегралом движения. Это позволяет понизить размерность исследуемой задачи и рассматривать двумерные уравнения движения в переменных X, Z. Гамильтониан движения в безразмерных переменных имеет следующий вид:

1 1 . 2 =

Pz) + 2 Mz)

bnax)z

(4)

где

ay{z) = -

a2 In ch (72 In ch

(z - Zç) a

(z + zc) a

2a2 In ch I ^ I

(5)

и использованы безразмерные переменные px,y,z =

— 1 х,у,

, _ {Y,Z,Zc\

X =

X

bn =

a =

ipo

VpoL

ларморовский радиус).

Ру

Ь„а '

mvo Bn

Во ' " \/ Po

Получая уравнения движения из гамильтониана (4), имеем

Рх = bncr(ay(z) - Ъпох), Pz = -bx(z)cr(ay(z) - bnax), ® = Рх ! Z = Pz,

где

1 /

bx(z) = - ( tanh-

tanh ■

а

а

(6)

(7)

(8)

(9)

Начальные условия выбирались следующим образом: координата го =0 фиксирована, координата х задается в диапазоне — 1 < Ъпахо < 1 с постоянным шагом, аналогично ^1<рха<1, рхо вычисляется из закона сохранения энергии (Рго)2 = 1 - (ЬпСГХо)2 - (рх0)2 .

Для численного интегрирования уравнений движения с заданными начальными координатами и скоростями был использован метод Рунге-Кутта четвертого порядка с переменным шагом.

Траектория частиц отслеживалась таким образом, что при пересечении частицей плоскости г = 0 на множестве рх, Ъпах отмечалась соответствующая точка. Совокупность всех таких точек на множестве рх,Ьпсгх представляет собой сечение Пуанкаре.

4. Результаты численного эксперимента

Для описания динамики частиц в «колоколооб-разных» токовых слоях часто используется пара-

минимальный радиус мак-

метр к = Яс/рь, где Нс кривизны магнитной силовой линии, р1 еимальный ларморовский радиус, в частности для простой параболической модели магнитного поля {Вх = -Вот! Ву = О' В* = Вп}, параметр к = Ьпа [3, 8]. В работе [16] параметр к был обобщен на случай магнитных полей, где полутолщина области обращения Ь меньше, чем гирорадиус вращения ро:

к =

Ъпа, а > 1,

Ъп{ 1 + ст), ег < 1.

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

Keff —

{KCeiï, Oeff > 1,

Ml+Ceff)) °eff < 1,

где параметр aes = \JL+p^c определяет отношение

области обращения поля к гирорадиусу вращения ро в случае двугорбой конфигурации. Здесь приводятся типичные сечения Пуанкаре, отображающие основные свойства рассматриваемых систем. На рис. 2 показано сечение Пуанкаре для колоколообразно-го слоя, которое соответствует значению параметров zc = 0, о = 1.78 или к = 0.178. Аналогичные сечения были построены для расщепленного TTC (рис. 4) с гс = 0.356, о = 1.78, параметр кев = 3.08. Сравнивая рисунки для колоколообразной и двугорбой конфигураций слоя, можно видеть, что общая структура фазовой плоскости остается неизменной, но при увеличении параметра zc (а фиксировано) происходит увеличение стохастической области

0.5

0

-0.5

■и

h „ах

о

0.5

Рис. 4. Сечение Пуанкаре для расщепленного TTC,

0.5

Рх

0

-0.5

О

0.5

Ъпах

1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

0.5

О

-0.5

Рх

7.

^ а -

v.-";-:-;;-.-/. J ......•»•-•• - '

' Л - Д. ; i.'-ftiil'^iV-,Г' • '—- - ' -

- w-.'

•N . -,««% -,

\Ч\ '

.'.4 •/< ijj-r.f . -ïv.v.v'i /

' 4.Т'4; -'¿о.?, 'i

' " 'л' . .. "

.' 'Л î V -v

. *Î 5-.'/: ¡V}' " • * ' ' Г

0

0.5

bnax

Pue. 5. Сечения Пуанкаре для «колоколообразного» (а) и расщепленного (б) TTC. Показана лишь часть фазового круга рх2 + (кх)2 ^ 1. Рис. 5, а соответствует значению параметров zc = 0, а = 0.86, рис. б — значению параметров

гс = 1, а = 0.5

на фазовой плоскости, при этом количество вытянутых «пальцеобразных» структур уменьшается. На обоих рисунках четко видны три области фазового пространства, соответствующие трем типам орбит. Пунктирная линия соответствует сепаратрисе, которую пересекают пролетные и квазизахваченные частицы. На рис. 5 изображены сечения Пуанкаре для колоколообразного (а) и двойного (б) слоев при очень малых cres = 0.86. Как видно из рисунка, фазовые области мало отличаются друг от друга, что свидетельствует о слабой зависимости траекторий частиц от структуры токового слоя. Поскольку при малых значениях параметра cres частицы обладают большой энергией и соответственно ларморовским радиусом много больше толщины слоя, то такие частицы не чувствуют тонкой структуры TTC и их рассеяние на двугорбом слое происходит так же, как и на одногорбом слое.

Выводы

Численное исследование динамики заряженных частиц в нераещепленных и расщепленных TTC показало, что при умеренно малых значениях к существуют отличия в объемах фазовых областей квазизахваченной плазмы. Было продемонстрировано, что расщепление слоя вызывает увеличение относительного объема рассеянной плазмы, что может играть критическую роль для эволюции и структуры тонкого токового слоя [6]. Управляющими параметрами, определяющими размеры фазовых областей, могут служить как расстояние между максимумами плотности тока Zc, так и отношение их ширины к ларморовскому радиусу ионов а. Увеличение фазовой области захваченных частиц сопровождается общим увеличением времени жизни популяции плазменных частиц в слое.

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

Работа выполнена при частичной финансовой поддержке РФФИ (грант 05-02-17003), ГФЕН (04-02-39021), программой научных школ (НШ-1739.2003.2) и JNTAS (03-51-3738).

Литература

1. Sonnerup B.U.O. // J. Geophys. Res. 1971. 76. P. 8211.

2. Harris E.G. 11 Nuovo Chimento. 1962. 23. P. 115.

3. Chen J., Palmadesso P.J. // J. Geophys. Res. 1986. 91.

P. 1499.

4. Hoshino M. Nishida A. Mukai T. et al. // J. Geophys. Res.

1996. 101. P. 24775.

5. Asatio Y.T., Mukai M., Hoshino Y. et al. // J. Geophys. Res.

2003. 108. P. 1019.

6. Zelenyi L.M., Delcourt D.C., Maloisa H.V. et al. // Geophys.

Res. Lett. 2002. 29. P. 49-1.

7. Speiser T. W. 11 J. Geophys. Res. 1965. 70. P. 4219.

8. Büchner J., Zelenyi L.M.// J. Geophys. Res. 1989. 94.

P. 11821.

9. Runois A., Nakamura R., Baumjohann W. et al. // Geophys.

Res. Lett. 2003. 30. P. 8-1.

10. Runov A., Nakamura R., Baumjohann W. et al. // Geophys. Res. Lett. 2003. 30. P. 33-1.

11. Лихтенберг А., Либерман М. Регулярная и стохастическая динамика. М., 1984.

12. Чириков Б.В. Нелинейный резонанс. Новосибирск, 1977.

13. Заславский Г.М. Стохастическая динамика систем. М., 1984.

14. Chen J, Mitchell H.G, Palmadesso P.J. // J. Geophys. Res. 1990. 95. P. 15141.

15. Alexeev /./., Malova H. V. // Advances in Space Research. 1995. 16. P. 205.

16. Савенков Б.В., Зеленый Л.М., Зогин Д.В. // Физика плазмы. 1997. 23. С. 436.

Поступила в редакцию 21.04.05

i Надоели баннеры? Вы всегда можете отключить рекламу.