Правительство Российской Федерации
Федеральное государственное бюджетное образовательное учреждение
высшего профессионального образования
«Санкт-Петербургский государственный университет»
Кафедра физической механики
Клюшников Георгий Николаевич
Динамика заряженных частиц в геомагнитном поле
Выпускная квалификационная работа
Допущен к защите.
Зав. кафедрой:
проф., д. ф.-м. н. Морозов В. А.
Научный руководитель:
проф., д. ф.-м. н. Колесников Е. К.
Рецензент:
к. ф.-м. н. Чернов С. В.
Санкт-Петербург
2016
SAINT-PETERSBURG STATE UNIVERSITY
Department of physical mechanics
Georgiy Klyushnicov
Dynamics of charged particles in a geomagnetic field
Graduation Qualification Paper
Admitted for defence.
Head of the chair
doctor of sc., prof. V. Morozov
Scientific supervisor:
doctor of sc., prof. E. Kolesnikov
Reviewer:
ass. of prof, PhD S. Chernov
Saint-Petersburg
2016
Содержание
1 РАЗРЕШЁННЫЕ ОБЛАСТИ В ЗАДАЧЕ О ДИНАМИКЕ
ЗАРЯЖЕННОЙ ЧАСТИЦЫ В СУПЕРПОЗИЦИИ ДИПОЛЬНОГО И ОДНОРОДНОГО МАГНИТНЫХ ПОЛЕЙ
7
1.1 Каноническая формулировка задачи о движении заряженной
частицы в суперпозиционном поле. Общее соотношение, задающее разрешённые области. . . . . . . . . . . . . . . . . .
7
1.2 Уравнения для границ разрешённой области. Качественные
особенности разрешённых областей. . . . . . . . . . . . . . .
10
2 ЗАДАЧА О ДИНАМИКЕ ВЕДУЩЕГО ЦЕНТРА ЗАРЯЖЕННОЙ ЧАСТИЦЫ В СУПЕРПОЗИЦИИ ПОЛЯ МАГНИТНОГО ДИПОЛЯ И ОДНОРОДНОГО МАГНИТНОГО ПОЛЯ
13
2.1 Вывод квадратуры для дрейфа ведущего центра . . . . . . .
13
2.2 Условия применимости дрейфовых уравнений движения . .
21
3 ОБЛАСТИ ВЫСЫПАНИЯ ЗАРЯЖЕННЫХ ЧАСТИЦ В
ГЕОМАГНИТНОМ ПОЛЕ, ПРИБЛИЖЕННОМ ПЕРВЫМИ СФЕРИЧЕСКИМИ ГАРМОНИКАМИ РЯДА ГАУССА
24
3.1 Алгоритм численного моделирования . . . . . . . . . . . . .
27
3.2 Основные результаты расчётов . . . . . . . . . . . . . . . . .
28
3
ВВЕДЕНИЕ
Задачи о движении высокоэнергетичных заряженных частиц в магнитном поле Земли давно привлекают внимание исследователей. Работы Штермера [13], направленные на математическое обоснование гипотезы Гольдштейна-Биркеланда о корпускулярной природе полярных сияний, и исследования Леметра и Валларта [12], [14], посвященные построению разрешенного конуса космических лучей, явились основой классической теории,
изучающей, главным образом, особенности распространения в магнитном
поле Земли заряженных частиц, приходящих с больших расстояний от центра Земли.
Вместе с тем в последние годы на передний план вышли задачи динамики в геомагнитном поле заряженных частиц высокой энергии, запускаемых в околоземное космическое пространство на конечном расстоянии
от Земли. К числу возможных источников указанных частиц относится
альбедное корпускулярное излучение, возбуждаемое в материале крупных
объектов космической инфраструктуры (многофункциональных долговременных орбитальных станций, солнечных электростанций, стартовых орбитальных платформ), появление которых в околоземном пространстве возможно в ближайшем будущем. Источниками частиц высокой энергии могут быть и космические ускорители высоких энергий, концепции которых в
настоящее время разрабатываются как в нашей стране, так и за рубежом.
Заметим, что отмеченные источники частиц высокой энергии могут явиться дополнительным фактором антропогенного загрязнения ближнего космоса, значительно более опасным для космической техники в силу высокой проникающей способности частиц высокой энергии по сравнению, например, с загрязнением техногенными микрочастицами, воздействующими
только на внешние поверхности космического аппарата и его аппаратуры.
Настоящая работа состоит из трёх глав. Первая глава посвящена иссле4
дованию областей, являющихся разрешёнными для движения заряженной
частицы в суперпозиции дипольного и однородного магнитных полей. При
этом основное внимание уделено построению разрешённых областей в ранее
не исследованном случае движения в суперпозиционном поле с магнитным
моментом диполя, сонаправленным внешнему магнитному полю, для положительных значений постоянной Штермера. Вид разрешённых областей
для случая, когда магнитный момент диполя направлен противоположно
индукции однородного поля, рассмотрен в работах [6, 7], а для случая, когда магнитный момент диполя сонаправлен индукции однородного поля и
постоянная Штермера отрицательна – в работе [8].
Во второй главе рассмотрена задача о динамике заряженной частицы в
суперпозиции поля магнитного диполя и однородного магнитного поля, параллельного магнитному моменту диполя. Указанная модель может быть
использована для описания магнитного поля Земли на достаточно больших
расстояниях от центра Земли (порядка десяти земных радиусов). На указанных расстояниях основной вклад в поле внутренних источников даёт
первая (дипольная) гармоника геомагнитного потенциала, а поле внешних
гармоник описывается моделью однородного магнитного поля. Известно,
что точное решение системы дифференциальных уравнений динамики заряженной частицы в магнитном поле, представляющем из себя суперпозицию дипольного и однородного магнитного полей, в общем случае не может
быть получено. В то же время в случае дипольного магнитного поля было
найдено (см. [1]) аналитическое решение динамической задачи с использованием так называемых дрейфовых уравнений движения, описывающих
движение ведущего центра заряженной частицы по соответствующей силовой поверхности дипольного магнитного поля. Во второй главе показано,
что использование дрейфовых уравнений движения позволяет получить
аналитическое решение и задачи о дрейфе ведущего центра заряженной
частицы в суперпозиционном поле.
5
Наконец, в третьей главе исследуются области высыпания на земную поверхность высокоэнергетических заряженных частиц для различных значений энергий и углов, задающих начальную скорость частицы. Число частиц
в ансамбле принимается равным 100000. При этом магнитное поле Земли
моделируется суперпозицией одной первой либо первых четырёх сферических гармоник геомагнитного потенциала.
6
1
РАЗРЕШЁННЫЕ ОБЛАСТИ В ЗАДАЧЕ
О ДИНАМИКЕ ЗАРЯЖЕННОЙ ЧАСТИЦЫ В СУПЕРПОЗИЦИИ ДИПОЛЬНОГО
И ОДНОРОДНОГО МАГНИТНЫХ ПОЛЕЙ
1.1
Каноническая формулировка задачи о движении
заряженной частицы в суперпозиционном поле. Общее соотношение, задающее разрешённые области.
В цилиндрической системе координат (ρ, ϕ, z) с осью z, сонаправленной
магнитному моменту M, релятивистский гамильтониан заряженной частицы с зарядом e имеет вид
s
HR = c Pρ2 + Pz2 +
Pϕ eAϕ
−
ρ
c
2
+ m20 c2 ,
(1.1)
где Pρ = mρ̇, Pϕ = mρ2 ϕ̇ + eρAϕ /c и Pz = mż – обобщённые импульсы,
соответствующие обобщённым координатам ρ, ϕ и z (точка обозначает проp
изводную по времени t), m0 – масса покоя частицы, m = m0 / 1 − v 2 /c2 ,
v – скорость частицы, c – скорость света (см. [4]).
Из (1.1) следует, что релятивистский гамильтониан HR не зависит от
времени явно. Значит, HR является интегралом движения:
HR = mc2 = const.
(1.2)
(1.2) означает, что релятивистская масса m остаётся неизменной во всё
время движения частицы. Следовательно, движение релятивистской частицы с массой покоя m0 можно рассматривать как движение нереляти7
вистской частицы массы m с классической функцией Гамильтона
2 !
1
Pϕ eAϕ
H=
Pρ2 + Pz2 +
−
,
2m
ρ
c
(1.3)
Из (1.3) видно что, гамильтониан H не зависит от координаты ϕ, поэтому
∂H
P˙ϕ = −
= 0.
∂ϕ
(1.4)
Таким образом, обобщённый импульс Pϕ является интегралом движения.
Представим его в виде
emΓ
,
c
Γ – постоянная размерности обратной длины.
Pϕ =
(1.5)
Нетрудно проверить, что в рассматриваемом случае суперпозиционного
магнитного поля векторный потенциал
M ρ B0 ρ
A=
+
eϕ ,
r3
2
(1.6)
где B0 – проекция однородного магнитного поля на ось z (см. также [11]).
С учётом (1.6) гамильтониан (1.3) можно записать в виде
H=
1
Pρ2 + Pz2 + V,
2m
(1.7)
где
2
ρ
B0 ρ
e2 M 2 Γ
V =
− −
.
(1.8)
2mc2 ρ r3 2M
Гамильтониан (1.7) описывает движение частицы в так называемой ве-
дущей плоскости, проходящей в любой момент времени через частицу и
магнитную ось и движущейся вместе с частицей. Такое движение эквивалентно движению в потенциальном поле с эффективным потенциалом
(1.8). Уравнения движения частицы в указанном поле могут быть записаны в виде
∂V
,
mρ̈ = −
∂ρ
∂V
mz̈ = −
.
∂z
8
(1.9)
Эффективный потенциал V не зависит от времени, поэтому полная энергия частицы E в ведущей плоскости является интегралом движения:
E=
m 2
ρ̇ + ż 2 + V = const.
2
(1.10)
Перейдём к безразмерным переменным, зависящим от штермеровской
p
единицы длины Cst = eM/mvc:
ρ̃ =
ρ
z
vt
, z̃ =
, t̃ =
.
Cst
Cst
Cst
(1.11)
Будем обозначать дифференцирование по t̃ штрихом. Тогда в новых переменных (1.11) уравнения (1.9) примут вид
∂ Ṽ
,
ρ̃00 = −
∂ ρ̃
z̃ 00 = − ∂ Ṽ ,
∂ z̃
где новый эффективный потенциал
2
3
ρ̃
ρ̃
1 ΓCst B0 Cst
−
− 3 .
Ṽ =
2
ρ̃
M
r̃
(1.12)
(1.13)
С учётом (1.11) запишем интеграл (1.10) следующим образом:
1 0 2
1
0 2
(r̃ ) + (z̃ ) + Ṽ (ρ̃, z̃) = .
2
2
(1.14)
Условимся в дальнейшем опускать значок ∼“. Как следует из (1.14),
”
разрешённые области движения частицы определяются неравенством V ≤
1/2, которое может быть представлено в виде
(1.15)
3
где γ = −Cst Γ/2 – постоянная Штермера, α = −B0 Cst
/2M .
Заметим, что уравнения (1.12) и (1.14) эквивалентны уравнениям движения точки единичной массы в эффективном потенциальном поле с потенциалом (1.13).
9
1.2
Уравнения для границ разрешённой области. Качественные особенности разрешённых областей.
Итак, в любой момент времени безразмерные координаты частицы должны удовлетворять неравенству (1.15), задающему так называемую разрешённую область движения. Точки, координаты которых не удовлетворяют
(1.15), принадлежат запрещённой области.
Полагая ρ = r cos ψ, где ψ – геомагнитная широта, получим из (1.15), что
В рассматриваемом случае индукция однородного поля параллельна магнитному моменту диполя, поэтому α < 0. В дальнейшем ограничимся рассмотрением ранее не исследованного случая положительных значений постоянной Штермера γ (см. также [3]).
Границы области (1.16) определяются уравнениями
2
4γ
2
r2 +
r
+
= 0,
µ cos ψ
µ cos2 ψ
µ
2
4γ
2
r3 −
r2 +
r
+
= 0,
µ cos ψ
µ cos2 ψ
µ
r3 +
(1.17)
(1.18)
где µ = −2α > 0 (см. [8]).
В уравнении (1.18) все коэффициенты при степенях r положительны.
Указанное уравнение не имеет вещественных положительных корней и,
следовательно, не может задавать границу области (1.16). Таким образом,
в рассматриваемом случае граница области (1.16) задаётся вещественными
положительными корнями уравнения (1.18), существующими при условии
D ≤ 0:
3
2 !
16
4
1
1
2
D=
γ
−
+
γ
−
+
9 cos6 ψ 3µ3
3µ
µ4
9µ
8
2
1
γ
−
+
≤ 0.
+ 3
3µ cos3 ψ
9µ
µ2
10
(1.19)
При условии (1.19) уравнение (1.18) имеет два положительных корня
r
r
p
p
ϕ + 4π
ϕ
r1 = 2 − cos
,
r2 = 2 − cos ,
(1.20)
3
3
3
3
где угол ϕ определяется соотношением
√
3 3q
;
cos ϕ = − p
2 −p3
4
1
2
4γ
8
p=
γ−
,
q=
1+
−
.
µ cos2 ψ
3µ
µ
3µ cos3 ψ 27µ2 cos3 ψ
Полагая cos3 ψ = u, запишем (1.19) в виде
A + Bu + Cu2 ≤ 0,
(1.21)
где
3
2 !
1
1
2
+ 4 γ−
,
γ−
3µ
µ
9µ
2
8
,
B= 3 γ−
µ
9µ
1
C = 2.
µ
16
A=
9
4
3µ3
Значения u должны удовлетворять условию
0 ≤ u ≤ 1.
(1.22)
Анализ показывает, что не все пары значений (µ, γ) (а значит, не все тройки коэффициентов (A, B, C)) имеют физический смысл, т. е. соответствуют
каким-либо начальным данным. При выполнении системы неравенств
B 2 > 4AC,
(1.23)
p
B 2 − 4AC > B.
пара (µ, γ) соответствует некоторым начальным данным, а решение (1.21)
в области (1.22) существует и может быть записано в виде
0 ≤ u ≤ min(u2 , 1),
11
(1.24)
где u2 =
√
B 2 − 4AC − B /2C.
Как видно из (1.24), качественный вид разрешённой области зависит от
значения корня u2 .
В случае u2 < 1 решение (1.21) в области (1.22) имеет вид 0 ≤ u ≤ u2 . При
таких значениях u граница разрешённой области существует для значений
магнитной широты
π
i h
√
√
π
3
3
ψ ∈ − , − arccos u2 ∪ arccos u2 ,
.
2
2
(1.25)
Как следует из (1.25), разрешённая область в случае u2 < 1 состоит из
двух компонент, границы каждой из которых не пересекают плоскость магнитного экватора. Каждому значению широты (1.25) соответствуют значения r1 и r2 , задаваемые формулами (1.20). Качественный вид разрешённой области в рассматриваемом случае показан на Рис. 1, для которого
√
γ = 1.17, µ = 0.2, u2 = 0.812, arccos 3 u2 = 21.078°, D(ψ = 0) = 12.182 > 0.
При u2 ≥ 1 решение (1.21) в области (1.22) имеет вид 0 ≤ u ≤ 1; граница
разрешённой области существует для любых значений магнитной широты
ψ. Если u2 = 1, то разрешённая область, как и при u2 < 1, состоит из
двух компонент, соединяющихся друг с другом в плоскости магнитного
экватора. Качественный вид разрешённой области для указанных значений
u2 показан на Рис. 2, для которого γ = 1.15186, µ = 0.2, D(ψ = 0) = 0.
Наконец, при u2 > 1 разрешённая область является однокомпонентной, что
видно из Рис. 3, соответствующего значениям параметров γ = 0.1, µ = 0.2,
D(ψ = 0) = −315.444 < 0.
12
2
ЗАДАЧА О ДИНАМИКЕ ВЕДУЩЕГО ЦЕНТРА ЗАРЯЖЕННОЙ ЧАСТИЦЫ В СУПЕРПОЗИЦИИ ПОЛЯ МАГНИТНОГО ДИПОЛЯ И ОДНОРОДНОГО МАГНИТНОГО ПОЛЯ
2.1
Вывод квадратуры для дрейфа ведущего центра
Дадим общую формулировку дрейфовых уравнений движения.
Рассмотрим частицу массы m и заряда q, движущуюся в магнитном поле B и поле сил немагнитной природы F. Будем предполагать, что выполнены условия слабой неоднородности поля B (см. [1]):
ρL |(∇B)⊥ | B,
v T |(∇B) | B,
k g
k
(2.1)
где ρL = mv⊥ /qB – ларморовский радиус, Tg = 2πm/qB – период ларморовского вращения, v⊥ – проекция скорости частицы на плоскость, перпендикулярную B, vk – составляющая скорости частицы в направлении B,
|(∇B)⊥ | – модуль проекции градиента B на плоскость, перпендикулярную
B, |(∇B)k | – модуль проекции градиента на направление B. В однородном поле, как известно, движение можно рассматривать как наложение
вращения по ларморовской окружности и движения центра ларморовской
окружности вдоль силовой линии. Если обозначить через r радиус-вектор
точки, в которой находится частица, то в каждый момент времени
r = a + c,
(2.2)
где a – вектор, проведённый из центра вращения (ведущего центра) в
точку, в которой находится частица; c – текущее положение ведущего центра.
13
При выполнении условий (2.1) движение частицы сохранит вид (2.2), но,
в отличие от движения в однородном поле, величины a и c медленно меняются со временем: a = a(t), c = c(t), а ведущий центр может двигаться
поперёк геомагнитного поля, перемещаясь с одной силовой линии на другую.
Скорость частицы v будет состоять из двух частей: скорости вращения
w и скорости дрейфа u:
v=
dr da dc
=
+
= w + u.
dt
dt
dt
(2.3)
С помощью метода усреднения из точного уравнения движения частицы
могут быть получены уравнения динамики ведущего центра (см. [5]):
du
mb
= Fb − µb∇B,
dt
!
!
2
(2.4)
2u
2cµu2k
du⊥
cB
k
u⊥ = F − µ 1 + 2 ∇B − m
×
+
(rotB)⊥ ,
w
dt
qB 2 w2 qB
где b = B/B, величины с индексами ⊥ обозначают проекцию на плоскость ϕ = const, движущуюся вместе с частицей (будем называть такую
плоскость ведущей), µ = p2⊥ /2γmB – эквивалентный магнитный момент
частицы, который при выполнении условий (8) является приближённым
интегралом движения (см. [1]).
При этом продольное движение ведущего центра описывается первым
уравнением (2.4), а поперечное – вторым уравнением (2.4).
В рассматриваемом случае сила F равна нулю (берутся в расчёт только силы магнитной природы). Также равно нулю и последнее слагаемое в
правой части второго уравнения (2.4) (вследствие безвихревого характера
поля). Кроме того, будем во втором уравнении (2.4) пренебрегать слагаемыми, содержащими производную
du⊥
dt
(вследствие малости u⊥ по сравнению
с v).
Пусть (r, λ, ϕ) – сферическая система координат, где r –расстояние до
центра Земли, λ – широта, ϕ – долгота. Тогда из (2.4) следует, что u⊥ имеет
14
одну отличную от 0 азимутальную составляющую uϕ :
dϕ
cm w2
uϕ = r cos λ
=− 2
+ u2k (∇B)⊥ .
dt
qB
2
(2.5)
Угол ϕ будем отсчитывать в направлении дрейфа частицы. Поскольку в
рассматриваемом случае проекция u⊥ на ведущую плоскость равна 0, в
этой плоскости ведущий центр в процессе движения остаётся на фиксированной силовой линии.
Пусть в начальный момент времени частица обладает скоростью v, направленной под углом θ к направлению местного магнитного поля. Тогда
w0 = v sin θ, uk0 = v cos θ – начальные значения w и uk , где величины с
индексом 0 соответствуют их значениям при t = 0. Учитывая, что в магнитостатическом поле модуль скорости v остаётся постоянным и используя
сохранение магнитного момента µ, можно показать (см. [1]), что текущее
значение продольной компоненты скорости ведущего центра
r
2µ(B1 − B)
uk =
,
m
(2.6)
где B1 = B0 / sin2 θ – значение индукции в «зеркальной» точке, где происходит отражение частицы от магнитного зеркала.
Введём безразмерную индукцию η, связанную с индукцией B и параметром re следующим образом:
B = M η/re3 .
(2.7)
Аналогично (2.7) запишем B1 :
B1 = M η0 /re3 = m(u2k + w2 )/2µ,
(2.8)
где η0 – значение η в точке, где происходит отражение от магнитного зеркала.
Скорость движения вдоль силовых линий поля равна
15
√
uk =
dr2
+
dt
r2 dλ2
s
=
dr
dλ
2
+ r2
dλ
.
dt
(2.9)
Функцию r(λ, α) найдём из уравнения силовых линий, которое удобно представить в виде (см. [2])
(2.10)
r3 + Cr/(α cos2 λ) − 1/α = 0,
где C = (1 − αre3 )/re .
Определим, какое значение re соответствует точке с радиус-вектором r .
Для этого решим кубическое уравнение
(2.11)
re3 + Cre /α − 1/α = 0.
Физический смысл имеют только положительные корни уравнения (2.11).
В зависимости от знаков D = 1/(4α2 ) + C 3 /(27α3 ) и α возможны 4 случая.
1. Если D > 0 и α > 0 (замкнутые кривые на рис. 1), то по формуле
Кардано получаем, что единственный положительный корень –
r
re =
3
1
2α
+
q
r
1
4α2
+
C3
27α3
+
3
1
2α
−
q
1
4α2
+
C3
27α3 .
Остальные два корня являются комплексными.
2. Если D < 0 и α > 0 (разомкнутые силовые линии), то уравнение
(2.11) имеет один положительный корень, определяющийся по формуле
p
re = 2 −C/(3α) cos(ψ/3), где угол ψ удовлетворяет условию cos ψ =
p
−27α/C 3 /2.
3. Если D 6 0 и α < 0 (разомкнутые силовые линии, не пересекающие экваториальную плоскость), то уравнение (2.11) не имеет положительных
корней. В этом случае один корень re вещественный и отрицательный,
а два другие комплексные.
16
4. Если D < 0 и α < 0 (кривые, пересекающие экваториальную плоскость, на рис. 2), то существуют два положительных корня re1 и re2 ,
выражения для которых имеют вид
p
re1 = 2 −C/(3α) cos(ψ/3),
где cos ψ = −
p
re2 = 2 −C/(3α) cos(ψ/3 + 4π/3),
p
−27α/C 3 /2.
Таким образом, в случаях 1 и 2 параметр re определён единственным
образом, в случае 3 уравнение для определения re не имеет положительных корней, а в случае 4 одному значению соответствуют 2 значения
re . Замкнутой силовой линии соответствует меньшее значение re , равное
p
re2 = 2 −C/(3α) cos(ψ/3 + 4π/3); силовой линии, уходящей на бесконечp
ность (разомкнутой) – большее значение re1 = 2 −C/(3α) cos(ψ/3).
Ограничимся движением ведущего центра вдоль замкнутых силовых поверхностей. С математической точки зрения это значит, что если уравнение
(2.11) при заданном C будет иметь два корня, то из них следует выбрать
наименьший.
Введём безразмерный параметр
δ = B0 re3 /M.
(2.12)
В (2.12) все величины кроме δ выражены в единицах СГС.
Тогда постоянную C можно записать в следующем виде:
C = (2 − δ)/(2re ).
(2.13)
Решение кубического уравнения (2.10), в котором и α выражены через
δ, зависит от знака соответствующего ему дискриминанта
1
(2 − δ)3
D= 2+
δ
27δ 3 cos6 λ
.
17
Будем представлять решение уравнения (2.10) в виде
r = F1 (λ, δ)re .
Если D > 0, то
v
v
s
s
u
u
u
u
3
3
3
(2 − δ)
(2 − δ)3
1
1
t1
t1
F1 (λ, δ) =
+
+
+
−
+
;
δ
δ 2 27δ 3 cos6 λ
δ
δ 2 27δ 3 cos6 λ
(2.14)
(2.15)
если же D < 0, то
r
F1 (λ, δ) = 2
где cos ψ =
1
δ
ψ + 4π
δ−2
cos
,
2
3δ cos λ
3
(2.16)
q
Производную
27δ 3 cos6 λ
(δ−2)3 , sin ψ > 0.
dr
dλ представим в виде
dr
= F2 (λ, δ) re .
dλ
(2.17)
При D > 0
s
3
(2 − δ) sin λ
(2 − δ)3
1
F2 (λ, δ) =
×
+
27δ 3 cos7 λ
27δ 3 cos6 λ δ 2
s
s
− 23
− 23
(2 − δ)3
1
1
(2 − δ)3
1
1
×
+
+
−
−
+
;
δ
27δ 3 cos6 λ δ 2
δ
27δ 3 cos6 λ δ 2
(2.18)
при D < 0
r
δ − 2 tg λ
ψ + 4π
cos
−
3δ cos λ
3
p
sin 2λ
ψ + 4π
− 3 δ (δ − 2) q
sin
,
3
3
δ δ (δ − 2) − 27 cos6 λ
F2 (λ, δ) = 2
где, как и ранее, cos ψ =
1
δ
q
27δ 3 cos6 λ
(δ−2)3 ,
sin ψ > 0.
Из (2.8) и определения штермеровской единицы длины
18
(2.19)
2
Cst
=
|q| re
c
s
re M
.
2mµη0
(2.20)
С учётом (2.6) – (2.9) получаем
dλ
1
=q
dt
r2 +
s
dr 2
dλ
2µM (η0 − η)
.
mre3
(2.21)
Теперь преобразуем вторую формулу (2.4). Сначала посчитаем модуль
вектора B:
s
B=M
2
δ
−
r3 re3
2
sin2 λ +
1
δ
+
r3 re3
2
cos2 λ.
(2.22)
Определим безразмерную индукцию η по формуле
q
η(r̂, λ, δ) = (2 − δr̂3 )2 sin2 λ + (1 + δr̂3 )2 cos2 λ,
где r̂ = r/re .
По формулам для градиента в сферической системе координат
3M 3 cos λ
1
1
4
(∇B)⊥ = 10 4 2
δ+ 3 +
2δ − 3 sin2 λ − δ + 3 cos2 λ
re F1 B
F1
F1
F1
1
2
+ sin2 λ
−
2δ
−δ
.
F13
F13
(2.23)
Выделим функцию F3 (λ, δ), не зависящую от r:
(∇B)⊥ =
M
F3 (λ, δ),
F14 re4 η 2
(2.24)
где
4
F3 (λ, δ) = 3 cos λ
2δ − 3 sin2 λ−
F1
1
1
1
2
2
2
− δ + 3 cos λ
δ + 3 + sin λ
− 2δ
−δ
.
F1
F1
F13
F13
^⊥ , поперечный полю B:
Определим безразмерный градиент (∇B)
19
F3 (λ, δ)
re4 (∇B)⊥
^
= 4
.
(∇B)⊥ =
M
F1 (λ, δ)η 2
(2.25)
Из (2.5), (2.8), (2.14) и определения µ следует, что
dϕ
cµM (η − 2η0 )
=
(∇B)⊥ .
dt
qB 2 re4 F1 (λ, δ) cos λ
(2.26)
Из (2.14), (2.20), (2.21) и (2.26) получаем производную
dϕ
=
dλ
M 1−
η
2η0
q
(∇B)⊥
r2 +
q
1−
2 B 2 cos λ
F1 re Cst
η
η0
dr 2
dλ
· sign(−q),
(2.27)
где sign(q) – знак q. Не умаляя общности рассуждения будем считать, что
sign(−q) = 1.
Наконец, (2.14), (2.17), (2.25) и (2.27) дают зависимость широты от долготы ведущего центра:
Zλ
ϕ − ϕ0 =
0
e=
где Φ
p
2 1 − η (∇B)
^Φ
e
⊥
2η0
re
q
dλ,
η
Cst
2
F1 η cos λ 1 −
(2.28)
η0
F12 + F22 .
На основе формулы (2.28) были построены графики зависимостей λ(ϕ).
Рассматривались только пространственно ограниченные движения, каковыми заведомо являются движения по замкнутым силовым линиям (2.10).
Таким силовым линиям соответствуют значения r(λ, α) в форме (2.14) –
(2.15) при 0 < δ < 2 и в форме ((2.14), (2.16)) при δ < 0. Соответствующие
графики для значений δ = 0.1, δ = −0.1 представлены на Рис. 4 – 5.
На каждом из графиков представлены кривые 1 – 4, соответствующие
следующим значениям угла θ между v и B: 10°, 30°, 40°, 80°. Как видно из
рисунков, при возрастании θ максимальная широта λ убывает.
20
2.2
Условия применимости дрейфовых уравнений движения
Вернёмся к условиям слабой неоднородности (2.1). Иначе можно записать их в виде
|(∇B)⊥ |
ρL
1,
B
(2.29)
Tg ∂B
1,
B ∂t
За условия применимости дрейфового приближения, как и в [4], примем
неравенства
ρL |(∇B)⊥ | 6 β 2 B,
1
v T |(∇B) | 6 β 2 B,
k g
k
2
(2.30)
где β1 и β2 – максимальные значения отношений ρL |(∇B)⊥ |/B и vk Tg |(∇B)k |/B.
Будем считать β12 и β22 такими же, как и для дипольного поля. По данным [4], β1 ≈ 0.46, β2 ≈ 0.66.
Для преобразования неравенств (2.30) вычислим теперь (∇B)⊥ и (∇B)k .
Если re = r/Cst , то из (2.22) следует, что
2M sin λ 1
M cos λ 1
Br =
− α , Bλ = −
+ 2α .
3
3
Cst
re3
Cst
re3
(2.31)
Из (2.31) найдём модуль вектора магнитной индукции:
s
2
2
q
M
1
1
2
B = Br2 + r2 Bλ2 = 3 4 3 − α sin λ +
+ 2α cos2 λ. (2.32)
Cst
re
re3
Удобно записать (2.32) в виде
B=
M B̃
3 ,
r̃Cst
(2.33)
q
где B̃(e
r, λ, α) = 4 (1 − αr̃3 )2 sin2 λ + (1 + 2αr̃3 )2 cos2 λ.
Проекция градиента B на направление поля B равна
1 ∂B
1 ∂B
3M 3 sin λ
1
(∇B)k =
Br +
Bλ = 2 10 4
4 α − 3 sin2 λ−
B ∂r
r ∂λ
B Cst re
re
(2.34)
1
2
1
1
− 3 + 2α cos2 λ
− 2α + cos2 λ 4α − 3
2α + 3
.
3
re
re
re
re
21
Введём обозначение
Fk (r̃, λ, α) = 6π sin λ 2 4 αr̃3 − 1 sin2 λ−
− 1 + 2αr̃3 cos2 λ 1 − αr̃3 + cos2 λ 4αr̃3 − 1 2αr̃3 + 1 .
Тогда
r, λ, α)
M 3 Fk (e
r, λ, α) M Fk (e
=
(∇B)k =
.
10 r
4r
e 2 Cst
2πB 2 Cst
e10
2π B
e4
Таким образом, первое неравенство (2.30) равносильно
e2
β2
B
q
k r , λ, α)
(2.35)
(2.36)
Таким же образом можно представить (∇B)⊥ в виде
(∇B)⊥ =
M F⊥ (e
r, λ, α)
,
4
4
2
e
Cst re B
(2.37)
где
4 αr̃3 − 1 sin2 λ − 2αr3 + 1 cos2 λ 2αr̃3 + 1 +
+2 sin2 λ cos λ 1 − 4αr̃3 1 − αr̃3 .
F⊥ (r̃, λ, α) = 3
Значит, второе неравенство (2.30) равносильно
re 6 p
β1
v⊥ /v
e2
B
p
.
|F⊥ (e
r, λ, α)|
(2.38)
Вместо неравенств (2.36) и (2.38) будем далее использовать достаточные
условия
e2
B
re 6 β1 p
,
|F⊥ (e
r, λ, α)|
e2
B
re 6 β2 q
(2.39)
Обозначим через Ω область, задаваемую неравенствами (2.39), а дополнение к ней – CΩ. Установим, как меняется вид областей Ω и CΩ при изме3
нении модуля индукции однородного магнитного поля B0 = 2M/Cst
. В ка-
честве заряженной частицы рассмотрим электрон с энергией E = 100 МэВ
22
[10]. Сначала разберём случай B0 ↑↓ M. Результаты расчётов показывают, что при переходе от дипольного поля (Рис. 6) к суперпозиционному
область CΩ (закрашена серым) становится ограниченной. Для значений
B0 < 17.39γ (1γ = 10−5 Гс) область Ω является двухкомпонентной (Рис.
7), а для B0 > 17.39γ – однокомпонентной (Рис. 8 – 9). При увеличении B0
площадь области CΩ уменьшается и в предельном случае B0 → ∞ область
Ω превращается во всю полуплоскость (см. также [9]).
Обратимся теперь к случаю B0 ↑↑ M. В этом случае, как и в рассмотренном ранее, для конечных значений B0 область CΩ является ограниченной.
Для значений B0 , меньших критического значения, равного 65.5γ, область
Ω является двухкомпонентной (Рис. 10), а для значений B0 , больших указанного критического значения, область Ω становится однокомпонентной
(Рис. 11 – 12). При увеличении B0 , как и в случае B0 ↑↓ M, площадь
области CΩ уменьшается, и в предельном случае B0 → ∞ область Ω превращается во всю полуплоскость.
23
3
ОБЛАСТИ ВЫСЫПАНИЯ ЗАРЯЖЕННЫХ
ЧАСТИЦ В ГЕОМАГНИТНОМ ПОЛЕ, ПРИБЛИЖЕННОМ ПЕРВЫМИ СФЕРИЧЕСКИМИ ГАРМОНИКАМИ РЯДА ГАУССА
Запишем уравнение движения заряженной частицы с массой покоя m0 и
зарядом e в геомагнитном поле с потенциалом V запишем в виде
e
d2 r
=
−
v × ∇V,
dt2
m0 γc
p
где c – скорость света, γ = 1 − v 2 /c2 – лоренц-фактор.
Потенциал V представляется в виде ряда Гаусса
n+1 X
∞
n
X
RE
m
V = RE
(gnm cos mϕ + hm
n sin mϕ)Pn (cos θ),
r
n=1
m=0
(3.1)
(3.2)
где RE – радиус Земли, gnm и hm
n – коэффициенты, определяемые в результате измерений геомагнитного поля, Pnm (cos θ) – многочлены Шмидта.
24
Таблица 1. Значения коэффициентов gnm и hm
n
n
m
gnm , 10−4 Гс
−4
hm
Гс
n , 10
1
0
-2944.2
0
1
-150.1
479.7
0
-244.5
0
1
301.3
-284.6
2
167.7
-64.2
0
135.1
0
1
-235.2
-11.5
2
122.6
24.5
3
58.2
-53.8
0
90.8
0
1
81.4
28.3
2
12.0
-18.9
3
-33.5
18.1
4
7.0
-33.0
2
3
4
Представим (3.1) в виде системы дифференциальных уравнений первого
порядка
dr
= v,
dt
e
dv
=−
v × ∇V.
dt
m0 γc
(3.3)
Индукция магнитного поля B выражается через потенциал V : B = −∇V .
Из (3.2) следует, что компоненты вектора B в сферической системе координат выражаются рядами
25
Br = −RE
n+2 X
∞
n
X
RE
n=1
∞
X
r
m
(gnm cos mϕ + hm
n sin mϕ)Pn (cos θ),
m=0
n+2 X
n
dPnm (cos θ)
RE
,
(gnm cos mϕ + hm
sin
mϕ)
n
r
dθ
n=1
m=0
∞
n
n+2
X RE
X
Pnm (cos θ)
m
m
Bϕ = RE
,
(mgn sin mϕ − mhn cos mϕ)
r
sin
θ
n=1
m=0
Bθ = −RE
(3.4)
В сферической системе координат с центром в центре Земли систему (3.3)
удобно представить в следующем удобном для численного интегрирования
виде
dr
dt
dθ
dt
dϕ
dt
dvr
dt
dvθ
dt
dv
ϕ
dt
= vr ,
=
=
=
=
=
vθ
,
r
vϕ
,
r sin θ
eωr vθ2 + vϕ2
+
,
mc
r
eωθ vr vθ vϕ2 ctg θ
−
+
,
mc
r
r
eωϕ vr vϕ vϕ vθ ctg θ
−
−
,
mc
r
r
(3.5)
где ωr = vθ Bϕ − vϕ Bθ , ωθ = Br vϕ − vr Bϕ , ωϕ = vr Bθ − vθ Br .
Положим B0 =
M
3
RE
(M – модуль магнитного момента Земли). Перейдём
в (3.5) к безразмерным переменным
Br
Bθ
Bϕ
, B̂θ =
, B̂ϕ =
,
B0
B0
B0
vr
vθ
vϕ
r
v̂r = , v̂θ = , v̂ϕ = , r̂ =
,
c
c
c
RE
B̂r =
ct
τ=
.
RE
В новых переменных (3.6) уравнения (3.5) принимают вид
26
(3.6)
dr̂
dt
dθ
dt
dϕ
dt
dv̂r
dt
dv̂θ
dt
dv̂
ϕ
dt
где
= v̂r ,
v̂θ
,
r̂
v̂ϕ
=
,
r̂ sin θ
=
v̂θ2 + v̂ϕ2
= Kωr +
,
r̂
v̂r v̂θ v̂ϕ2 ctg θ
= Kωθ −
+
,
r̂
r̂
v̂r v̂ϕ v̂ϕ v̂θ ctg θ
= Kωϕ −
−
,
r̂
r̂
(3.7)
ω̂r = v̂θ B̂ϕ − v̂ϕ B̂θ , ω̂θ = B̂r v̂ϕ − v̂r B̂ϕ , ω̂ϕ = v̂r B̂θ − v̂θ B̂r ,
2
K = sign(q) CRst2 – безразмерный коэффициент, Cst – штермеровская едиE
ница длины.
3.1
Алгоритм численного моделирования
Пусть заряженные частицы запускаются с околоземное космическое пространство из точечного источника, находящегося на геостационарной орбите, с фиксированной кинетической энергией и произвольными направлениями скорости. При этом некоторые частицы ансамбля будут падать на
поверхность Земли, а некоторые будут улетать от центра Земли на бесконечно большое расстояние. Будем называть областью высыпания множество точек вида (ϕf , λf ), где ϕf и λf – соответственно долгота и широта
точки земной поверхности, в которую попала частица.
В качестве заряженной частицы будем рассматривать электрон. Направление начальной скорости частицы будем характеризовать углом с местной вертикалью i и углом между проекцией вектора скорости на плоскость
местного горизонта и восточным направлением q. Предположим, что i и q
являются случайными величинами с равномерным законом распределения
в промежутках [2π/3, π] и [−π, π] соответственно. Для каждого электро27
на ансамбля производится численное интегрирование уравнений движения
(3.7) разностным методом Рунге-Кутта-Мерсона четвёртого порядка с автоматическим выбором длины шага. Вычисления проводились до тех пор,
пока не оказывалось выполненным одно из следующих условий:
1. Траектория электрона пересекает сферу радиусом RE ;
2. Длина траектории превышает установленное предельное значение 20RE ;
3. Электрон удаляется от центра Земли на столь большое расстояние, что
его дальнейшее движение сопровождается неограниченным монотонным ростом координаты r (указанное расстояние принимается равным
10RE ).
3.2
Основные результаты расчётов
С использованием описанного алгоритма были построены области
высыпания на земную поверхность электронов с энергиями 15, 30 и 60
МэВ. Считалось, что источник излучения находится в плоскости географического экватора на высоте 5.63 RE на нулевом меридиане. Для
сравнения областей высыпания при n = 1и n = 4 были получены изображения областей высыпания электронов рассматриваемых энергий на
поверхность Земли с учётом только первой гармоники геомагнитного
поля, а также с учетом первых четырёх сферических гармоник. Полученные результаты представлены на Рис. 13 – 18.
Также для каждого значения энергии электронов было определено
число точек высыпания. Результаты расчётов приведены в Табл. 2.
28
Таблица 2. Число точек высыпания
на земную поверхность при n = 1 и n = 4
Ek , ГэВ
n=1
n=4
15
5263
4288
30
7414
7479
60
12813
13172
Результаты исследований, приведённые в данной главе, были получены
с использованием вычислительных ресурсов Ресурсного Центра «Вычислительный центр СПбГУ».
29
Заключение
Кратко сформулируем основные результаты, полученные в главах 1 – 3.
1. Как показывают результаты приведённого анализа, в случае, когда
B0 ↑↑ M и γ > 0, разрешённая область Ω всегда является неограниченной. При этом в зависимости от параметра u2 Ω может состоять
как из одной, так и из двух компонент. При 0 < u2 < 1 в структуре
разрешённой области можно выделить две неограниченные компоненты, лежащие в верхней (z > 0) и нижней (z < 0) полуплоскости. При
u2 = 1 компоненты Ω соединяются друг с другом в плоскости магнитного экватора. Для значений u2 > 1 Ω однокомпонентна, причём её
ближняя и дальняя границы существуют для любых значений магнитной широты ψ.
2. Доказано, что приближение ведущего центра применимо к описанию
динамики заряженной частицы, движущейся в поле, представляющем
собой суперпозицию дипольного и однородного магнитных полей. Выведена квадратура, описывающая дрейф ведущего центра по силовой
поверхности суперпозиционного поля. С использованием этой квадратуры построены зависимости широты от долготы, соответствующие
различным направлениям начальной скорости частицы. Определены
области применимости дрейфового приближения для суперпозиционного поля. Полученные результаты позволяют сформулировать основные качественные особенности областей Ω и CΩ:
(a) В отличие от дипольного случая для конечного значения B0 запрещённая область Ω является ограниченной.
(b) Как при B0 ↑↑ M, так и при B0 ↑↓ M для значений B0 , меньших
некоторого критического значения, разрешённая область Ω являет30
ся двухкомпонентной, а при значениях B0 , больших критического,
область Ω является однокомпонентной.
(c) C ростом величины внешнего поля B0 площадь области CΩ уменьшается и в пределе B0 → ∞ разрешённой областью становится вся
полуплоскость.
3. Установлено, что независимо от модели геомагнитного поля области
высыпания электронов на поверхность Земли качественно не отличаются. При низких энергиях область остается разделенной недостижимой для инжектируемых частиц областью, обладая при этом симметрией. При этом вторая, третья и четвёртая гармоники в совокупности
оказывают заметное влияние на форму области высыпания.
31
Список литературы
1. Альвен Г., Фельтхаммар К.-Г.: Космическая электродинамика. Основные принципы М., Мир, 1967, 260 с.
2. Колесников Е. К., Клюшников Г. Н. Задача о движении ведущего центра заряженной частицы в суперпозиции магнитного диполя и однородного магнитного поля// Вестник СПбГУ, серия 1. 2015, №1, с. 123
– 134.
3. Колесников Е. К., Клюшников Г. Н. Разрешённые области в задаче о
динамике заряженной частицы в суперпозиции дипольного и однородного магнитных полей// Вестник СПбГУ, серия 1. 2016, №2, с. 293 –
299.
4. Колесников Е. К., Филиппов Б. В. Некоторые задачи эволюции заряженных частиц в поле магнитного диполя. Изд. ЛГУ, 1974. 72 с.
5. Ленерт Б.: Динамика заряженных частиц. М.: Атомиздат, 1967, 351 c.
6. Шалимов В. П., Швачунов И. Н. Изучение движения заряженных частиц в поле магнитного диполя, находящегося во внешнем магнитном
поле, методом Штермера. I //Космические исследования. 1966. Т. IV.
Вып. 2. C. 208 – 220.
7. Шалимов В. П., Швачунов И. Н. Изучение движения заряженных частиц в поле магнитного диполя, находящегося во внешнем магнитном
поле, методом Штермера. II //Космические исследования. 1966. Т. IV.
Вып. 3. C. 395 – 403.
8. Katsiaris G. A., Psillakis Z. M. Allowed regions for the motion of charged
particles in superposed dipole and uniform magnetic fields//J. Astrophysics
and Space Science. Vol. 126. P. 69 – 87. 1986.
32
9. Klyushnikov G. Precipitation regions in the problem of charged particles
dynamics in the Earth’s magnetic field//International Student Conference
«Science and Progress». Abstracts. Nov. 2015.
10. Kolesnicov E., Kluishnicov G. About drift approach applicability regions
in the problem for dynamics of a charged particle in a superposition
of dipole and uniform magnetic fields //International Conference
on Mechanics - Seventh Polyakhov’s Reading, 2015/Feb. 2015. DOI:
10.1109/POLYAKHOV.2015.7106736.
11. Lemaire J. F. The effect of a southward interplanetary magnetic field on
Stormer’s allowed regions//Space Res. Vol. 31, №5. P. 1131 – 1153. 2003.
12. Lemaitre G., Vallatra M. S. On Compton’s latitude effect of cosmic
radiation//Phys. Rev. V. 43, № 2. P. 87-91. 1933.
13. Stormer C. The polar Aurora. London-New-York: Oxford university press,
1955.
14. Vallatra M. S. On the allowed cone of cosmic radiation. Toronto: University
press, 1938.
33
Рис. 1: γ = 1.17, µ = 0.2
Рис. 2: γ = 1.15186, µ = 0.2
Рис. 3: γ = 0.1, µ = 0.2
Рис. 1 – 3. Разрешённые области движения(запрещённые области светло-серые, разрешённые
тёмно-серые).
34
Рис. 4: Графики λ(ϕ) при δ = 0.1 (L – широта λ, f i – долгота ϕ)
Рис. 5: Графики λ(ϕ) при δ = −0.1 (L – широта λ, f i – долгота ϕ)
Рис. 6: B0 = 0γ, rm = 60 Rz
35
Рис. 7: B0 ↑↓ M, B0 = 16γ, rm = 60 Rz
Рис. 8: B0 ↑↓ M, B0 = 17.39γ, rm = 60 Rz
Рис. 9: B0 ↑↓ M, B0 = 18γ, rm = 60 Rz
36
Рис. 10: B0 ↑↑ M, B0 = 60γ, rm = 60 Rz
Рис. 11: B0 ↑↑ M, B0 = 65.5γ, rm = 60 Rz
Рис. 12: B0 ↑↑ M, B0 = 70γ, rm = 60 Rz
37
Рис. 13: Ek = 15GeV , n = 1
Рис. 14: Ek = 30GeV , n = 1
38
Рис. 15: Ek = 60GeV , n = 1
Рис. 16: Ek = 15GeV , n = 4
39
Рис. 17: Ek = 30GeV , n = 4
Рис. 18: Ek = 60GeV , n = 4
40
Отзывы:
Авторизуйтесь, чтобы оставить отзыв