Мазур В. А. МГД-резонатор в лобовой части магнитосферы и его

advertisement
Солнечно-земная физика. Вып. 14. (2009) С. 52–66
УДК 533.951
МГД-РЕЗОНАТОР В ЛОБОВОЙ ЧАСТИ МАГНИТОСФЕРЫ И ЕГО ВОЗБУЖДЕНИЕ
ВНЕМАГНИТОСФЕРНЫМИ ГИДРОМАГНИТНЫМИ ВОЛНАМИ
(МОНОХРОМАТИЧЕСКИЕ КОЛЕБАНИЯ)
В.А. Мазур
MHD-RESONATOR IN THE FRONT PART OF THE MAGNETOSPHERE AND ITS EXCITATION BY
EXTRAMAGNETOSPHERIC HYDROMAGNETIC WAVES (MONOCHROMATIC OSCILLATIONS)
V.А. Mazur
В рамках одномерно-неоднородной модели среды теоретически исследованы свойства МГД-резонатора в лобовой
части магнитосферы (его собственные частоты и собственные моды) и его возбуждение монохроматическими гидромагнитными волнами, падающими на магнитосферу из солнечного ветра. Декремент затухания собственных мод обусловлен поглощением в точках альфвеновского резонанса и излучением в солнечный ветер через неполностью отражающую
магнитопаузу. Показано, что влияние магнитосферного резонатора приводит к резонансному характеру уровня накачки
магнитосферы падающими волнами. В магнитосферу проникают волны только в узких спектральных диапазонах вблизи
собственных частот резонатора, ширина которых порядка декремента затухания собственных мод. Волны с другими
частотами отражаются от магнитопаузы.
The one-dimensional inhomogeneous model of medium was applied to study theoretically the MHD-resonator properties in
the front part of the magnetosphere (its eigenfrequencies and eigenmodes) and its excitation by monochromatic hydromagnetic
waves incident on the magnetosphere from the solar wind. Damping factor of eigenmodes is determined by the absorption at
Alfvén resonance points and by the radiation into the solar wind through the partially reflecting magnetopause. The magnetospheric resonator effect was demonstrated to result in the resonant behavior of the level of pumping of the magnetosphere by
incident waves. Waves penetrate into the magnetosphere only in narrow spectral ranges near the resonator’s eigenfrequencies
whose width is of order of eigenmode damping factor. Waves with other frequencies are reflected from the magnetopause.
Введение
Проникновение гидромагнитных волн из солнечного ветра в магнитосферу является одним из
основных механизмов возбуждения в ней гидромагнитных колебаний. Достаточно сказать, что геомагнитные пульсации Рс3, практически постоянно существующие в дневной полусфере, энергия которых
превосходит энергию всех остальных пульсаций
вместе взятых [Гульельми, Троицкая, 1973], возбуждаются именно таким механизмом [Guglielmi, 1974;
Гульелми и др., 1976; Potapov, Mazur, 1984]. К этому
же классу можно отнести и пульсации, возникающие в результате событий SC (sudden commencement)
[Baumjohann, et al., 1984; Lysak, Lee, 1992; Fujita,
Tanaka, 2006].
В теоретических работах, посвященных такому
механизму, основное внимание уделяется возбуждению проникающей в магнитосферу волной альфвеновского резонанса [Kivelson, Southwood, 1986;
Mann, et al., 1995; Walker, 2002]. Но в этих работах
не учитывается тот факт, что магнитосфера может
обладать свойствами МГД-резонатора. Между тем,
возможность того, что роль резонатора может играть как магнитосфера в целом, так и различные ее
области, давно обсуждается в литературе [Kivelson,
Southwood, 1985; Walker, 1998; Leonovich, Mazur,
2006]. Совершенно ясно, что наличие в магнитосфере МГД-резонатора должно решающим образом
влиять на эффективность проникновения в нее гидромагнитных волн, амплитуду колебаний в ней и
уровень возбуждения альфвеновского резонанса.
Настоящая работа посвящена теоретическому исследованию этого круга вопросов.
Мы ограничиваемся рассмотрением лобовой части магнитосферы и прилегающей к ней области сол52
нечного ветра (magnetosheath). Представляется, что
наиболее естественно падение гидромагнитных волн
из солнечного ветра именно на лобовую часть магнитосферы. Кроме того, к лобовой части магнитосферы примыкает застойная область солнечного
ветра, поэтому исследование не осложняется сдвиговым течением на магнитопаузе, не проявляется
неустойчивость Кельвина–Гельмгольца, столь важная на флангах магнитосферы.
Для описания указанных областей магнитосферы
и солнечного ветра мы используем одномернонеоднородную модель среды, в которой учитывается только радиальная неоднородность. Именно эта
неоднородность приводит к запиранию МГД-волны
между магнитопаузой и точкой поворота в глубине
магнитосферы.
В настоящей статье излагаются результаты
первой части работы, посвященной монохроматическим колебаниям (т. е. колебаниям с фиксированной частотой). Результаты такого исследования представляют самостоятельный интерес, кроме того, они являются основой для изучения нестационарных и стохастических колебаний, которым посвящены вторая и третья части работы,
публикуемые в этом же номере, непосредственно
за настоящей статьей.
Модель среды
Схематическое изображение лобовой части магнитосферы и прилегающей к ней части солнечного
ветра, а также используемая нами одномерная модель этих областей представлены на рис. 1. Из этого
рисунка видно соответствие между различными
элементами модели и реальной системы.
Координата x декартовой системы координат, направленная вдоль неоднородности, меняется в модели
МГД-резонатор в лобовой части магнитосферы и его возбуждение…
рии отношение сSW/cAM будет считаться малым параметром. Оно связано с большим скачком плотности на магнитопаузе. Из условия равновесия нетрудно получить, что
cSW / cAM ~ (ρ M / ρ W )1/ 2 .
(1)
Исходные уравнения
Пространственно-временную зависимость возмущенной величины ξx, проекции вектора смещения
r
плазмы ξ, представим в виде
ξ x ( x, y , z , t ) = ξ% x ( x) U N ( z ) e
i ( k y y − ωt )
.
(2)
U N ( z ) = sin k z z , k z = k zN ≡ πN / L,
(3)
Здесь
где N=1, 2, 3, … – продольное (вдоль магнитного
поля) волновое число. Зависимость от z вида (3) следует из граничного условия
ξ x |z = 0, L = 0,
Рис. 1. Соответствие между элементами реальной среды и ее одномерной модели.
которое в приближении идеальной гидродинамики
вытекает из предполагаемой нами идеальной проводимости ионосферных плоскостей. Учитывая циклический характер азимутальной координаты y в
реальной магнитосфере, будем считать, что ky=m/a,
где m=0, 1, 2, … – азимутальное волновое число, a –
характерное значение экваториального радиуса, отвечающего силовой линии длины L. Таким образом,
решение (2) определяется волновыми числами m и N.
При исследовании возбуждения магнитосферных
колебаний внемагнитосферными волнами наибольшее значение имеют гармоники с наименьшими
значениями N, поскольку они более эффективно
возбуждаются внешней волной. Поэтому для оценок
мы будем считать, что N=1. Для азимутального волнового числа будем считать характерными значения
m=3÷5. Меньшие значения m неоправданны из-за
реально существующей азимутальной асимметрии
магнитосферы и еще более из-за азимутальной
асимметрии падающей волны, а для больших значений резко уменьшается эффективность возбуждения.
Предполагая применимость приближения идеальной гидродинамики, можно получить уравнение, описывающее структуру колебания по координате x [Duhau, Gratton, 1975]:
в пределах (–∞, ∞). Точка x=xM соответствует магнитопаузе. По оси z плазма в модели ограничена
проводящими плоскостями z=0, L, изображающими
магнитосопряженные области ионосферы. Чтобы
сохранить одномерный характер модели, мы примиряемся с наличием этих плоскостей и в области солнечного ветра x>xM. Ось y моделирует азимутальную координату магнитосферы.
r
r
Магнитное поле B0 = B0 ( x) будем считать направленным по оси z как в магнитосфере, так и в
солнечном ветре. Впрочем, как будет видно ниже,
направление магнитного поля в солнечном ветре не
играет большой роли.
На рис. 2 представлены схематические графики
зависимостей от x существенных для нашей задачи
параметров среды. Отметим, что в солнечном ветре
для нас будет важен только тот факт, что cS >> cA ,
конкретные же зависимости cS(x) и cA(x) в этой области существенного значения не имеют.
Для оценок по порядку величины примем следующие значения: L=1.5·105 км (что соответствует
силовой линии геомагнитного поля, экваториальный
радиус которой равен a=7RE=5·104 км), сAM=103 км/с,
сSW=2·102 км/с, сAW=50 км/с. В развиваемой ниже тео-
d ρ 0 Ω 2A dξ% x
+ ρ 0 Ω A2 ξ% x = 0.
dx K 2 dx
(4)
Здесь ρ0=ρ0(x) – плотность плазмы,
Ω 2A = ω2 − k z2 cA2 ,
K2 =
ω4
− kt2 ,
ω (c + cS2 ) − k z2 cA2 cS2
2
2
A
kt2 = k y2 + k z2 .
В дальнейшем мы предпочтем вместо уравнения
(4) для ξx использовать уравнение для функции
ρ Ω 2 dξ%
% x) = − 0 2 A x ,
ς(
(5)
dx
K
Рис. 2. Схематические графики скорости Альфвена
сA(x) и скорости звука cS(x) в магнитосфере и солнечном
ветре.
которое имеет вид
53
В.А. Мазур
d 2 ς% d ln ρ 0 Ω A2 dς%
−
+ K 2 ς% = 0.
dx
dx
dx 2
где величина
(6)
ik y
1 dς% %
%
, ξy =
ς.
2
ρ 0 Ω A dx
ρ 0 Ω 2A
точке xM. Из формул (7) и (11) видно, что условия
(12) означают непрерывность на магнитопаузе смещения плазмы ξx и полного возмущенного давления P.
Общее выражение для закона сохранения энергии МГД-волн имеет вид [Ландау, Лифшиц, 1982]
r
∂W
+ div S + Q = 0.
(13)
∂t
Здесь
r
r
r
B 2 ρ 0 v 2 pρ r c r r
 E B  + pv
+
+
,S =
W=


8π
2
2ρ0
4π
– плотность энергии и плотность потока энергии
волны, Q – плотность диссипируемой мощности, а
r
v – возмущенная скорость плазмы. В рассматриваемом одномерном случае уравнение (13) принимает
вид
(7)
r
Компоненты возмущенных электрического E и
r
магнитного B полей более простым образом выраr
жаются через компоненты вектора смещения ξ , но
при необходимости, используя соотношения (7), их
можно выразить и через ς% . Имеем (далее для кратi ( k y − ωt )
∂W ∂S
+
+ Q = 0,
∂t
∂x
кости опускаем множитель e y
)
Ex = E% x ( x) U N ( z ), E y = E% y ( x) U N ( z ), Ez = 0,
Bx = B% x ( x) VN ( z ), By =
B% x = k z B0 ξ% x , B% y = k z B0 ξ% y ,
(9)
×( f g * + f * g ). Усреднение по z функций U N2 ( z ) и
(10)
VN2 ( z ) дает множитель 1/2, а среднее от произведения U N ( z )VN ( z ) равно нулю.
Учитывая сказанное, имеем
r
r
1  | B% |2 ρ 0 ω2 | ξ% |2 p% • ρ% + p% ρ% • 
,
W= 
+
+
(15)

4  8π
2
4ρ0


B% z = −( B0 ξ% x )′ − ik y B0 ξ% y .
Для возмущенного давления
p = p% ( x) U N ( z ), p% =
= − p0′ ξ% x +
ω2 cS2
%
ς.
ω (c + cA2 ) − k z2 cS2 cA2
2
2
S
Sx =
Наконец, для полного возмущенного давления:
P = p+
B0 Bz
% x) U N ( z ).
= ς(
4π
{ς%}x
M
c % %* %*
ω
( E y Bz + E y Bz ) + i ( p% ξ% *x − p% *ξ% x ).
32π
8
(16)
Что касается диссипируемой мощности, то в
рамках приближения идеальной МГД она равна нулю. Однако в точке альфвеновского резонанса диссипацией нельзя пренебрегать в принципе – она остается конечной даже при стремлении диссипативных коэффициентов к нулю. Подробнее этот вопрос
рассмотрен ниже, где и приведена формула для Q .
Используя соотношения (7–10), выражение (16)
можно привести к виду
(11)
Из последнего соотношения ясен физический
% x).
смысл используемой нами функции ς(
Граничное условие для уравнения (6) при x→–∞,
где имеет место область непрозрачности, сводится к
требованию ограниченности решения, что фактически приводит к его экспоненциальному падению.
При x→∞ в области прозрачности (распространения
волны) граничное условие зависит от постановки
задачи. В задаче о собственных модах резонатора
таким условием является требование убегания волны, а в задаче о падении волны на магнитосферу
задается амплитуда падающей волны.
На магнитопаузе, при x=xM, коэффициенты уравнения (6) испытывают скачки. Из этого уравнения
следуют условия сшивки
 1 d ln ς% 
= 0, 
 = 0,
2
 ρ 0 Ω A dx  xM
(14)
где черта над буквой означает усреднение по координате z и по фазе φ=kyy–ωt (последнее равносильно
усреднению по координате y или по времени t).
Использование комплексной записи для возмущенных величин F=feiφ, где f=f1+if2, подразумевает, что
в действительности F=Refeiφ=f1cos φ–f2sin φ. Поэтому
для усредненных билинейных комбинаций имеем
FG → Re F ⋅ Re G = (1/ 2)( f1 g1 + f 2 g 2 ) = (1/ 4) ×
(8)
ωB
ωB
E% x = i 0 ξ% y , E% y = −i 0 ξ% x .
c
c
Компоненты Bx и By имеют другую зависимость
от координаты z, описываемую функцией
VN ( z ) = cos k z z :
= B% y ( x) VN ( z ), Bz = B% z ( x) U N ( z ),
означает скачок величины f в
M
Некоторое преимущество этого уравнения состоит
в том, что в нем эффекты альфвеновского резонанса и
точек поворота описываются разными членами.
Через функцию ς% можно выразить все другие
возмущенные величины. Для компонент вектора
r
смещения ξ имеем
ξ% x =
{ f }x
Sx = i
ω  dς% *
dς% 
ς% − ς% *  .
2 
8ρ0 Ω A  dx
dx 
(17)
Аналогичное выражение для W слишком громоздко, и мы приведем его в следующем разделе
только для магнитосферной части модели, где оно
значительно упрощается.
Приближение ВКБ
Аналитическое решение основного уравнения (6)
в общем случае, конечно, невозможно. Поэтому для
его решения будем использовать приближение ВКБ.
(12)
54
МГД-резонатор в лобовой части магнитосферы и его возбуждение…
Решения (20) могут быть написаны для каждой
из областей, в которых применимо приближение
ВКБ, причем с различными значениями констант Ci.
В разделяющих эти области окрестностях особых
точек основное уравнение (6) приводится к определенным стандартным формам, решение которых
хорошо известно. После этого методом сращивания
асимптотических разложений получают единое решение для всех x.
Решение (20) в области прозрачности можно
трактовать как суперпозицию падающей и отраженной волн. Запишем его следующим образом:
Как будет видно ниже, для интересующих нас колебаний это приближение, безусловно, выполняется в
области солнечного ветра, x>xM. В магнитосферной
области, x<xM, оно выполняется для высоких гармоник магнитосферного резонатора, но находится на
границе своей применимости для наиболее важных
первых гармоник. Однако приближение ВКБ обладает замечательным свойством – оно дает качественно правильные результаты даже на границе своей применимости. Поскольку характер используемой нами модели среды позволяет рассчитывать
только на результаты качественного характера, мы
без колебаний будем использовать формулы приближения ВКБ.
Согласно методу ВКБ решение ищется в виде
ς% ( x) = exp[iΨ ( x)] . Для применимости метода необходимо, чтобы фаза Ψ ( x) быстро менялась, т. е.
чтобы большим был параметр
S = l | d Ψ / dx |>> 1,
1/ 2
 ρ Ω2 
% x) = C0  0 A 
ς(
 kx 

 x

 x
 
× exp  −i ∫ k x ( x′)dx′ + R exp i ∫ k x ( x′)dx′  .
 xR

 xR
 

где l – характерный масштаб изменения коэффициентов уравнения. Решение для фазы ищется в виде асимптотического разложения по степеням 1/S:
Sx =
(19)
Ψ 0 ( x) = ± ∫ k x ( x)dx, k ( x) = K ( x),
2
ρ0 Ω 2A =
k x ( x)
1
Ψ1 ( x) = −i ln C + i ln
,
2
ρ 0 ( x)Ω A2 ( x)
K2 =
где C – произвольная постоянная. С учетом этих двух
порядков общее решение представляется в виде
% x) =
ς(
ρ0 Ω
kx
1
2 2
A
(
)
(
)
ς% ′′ −
(20)
Неравенство (18) при подстановке в него нулевого приближения (19) дает условие
ω2
− kt2 = q 2 − k y2 .
cA2
(25)
( q 2 )′
ς% ′ + (q 2 − k y2 )ς% = 0.
q2
(26)
Применяя к уравнению (26) приближение ВКБ,
имеем
(21)
k x2 ( x) = q 2 ( x) − k y2 ≡
Оно означает, что длина волны возмущения 1/kx
должна быть много меньше масштаба неоднородности. Однако это условие необходимо, но недостаточно. Полная формулировка условия применимости приближения ВКБ имеет вид | Ψ ′0 ( x ) |>>| Ψ1′ | .
Отсюда вытекают два неравенства:
| k x |3 >>| dk x2 / dx |, | k x Ω A2 |>>| d Ω A2 / dx | .
B02 2 2 ω2
q , q = 2 − k z2 ,
4π
cA
Учитывая также, что B0=const, основное уравнение (6) приводим к виду
C0 exp −i k x dx + C1 exp i k x dx  .
∫
∫


| k x | l >> 1.
1
1
ω ( − | C0 |2 + | C1 |2 ) = − ω | C0 |2 (1− | R |2 ) . (24)
4
4
Решение для магнитосферы
В этой области cS=0, поэтому
Два первых члена разложения имеют вид
2
x
(23)
Величину R будем называть комплексным коэффициентом отражения. Для решений (20) и (23) плотность потока энергии (17) принимает простой вид:
(18)
Ψ ( x) = Ψ 0 ( x) + Ψ1 ( x) + ...
×
ω2
− kt2 ,
c ( x)
2
A
где обозначено kt2 = k y2 + k z2 . Так как мы предполагаем, что функция cA(x) монотонно падает на интервале (–∞, xM), то уравнение для точки поворота
k x2 ( xR ) = 0 имеет в этом интервале решение, если
ω > kt cAM .
(22)
(27)
При этом интервал (–∞, xR) – это область непрозрачности (k x2, < 0), а интервал (xR, xM) – область
Эти условия очевидным образом нарушаются в
некоторых окрестностях точки поворота xR, в которой k x2 ( xR ) = 0 , и точки альфвеновского резонанса
прозрачности (k x2 > 0). Для существования собственных мод в магнитосфере необходимо наличие в
ней области прозрачности, поэтому мы будем предполагать выполнение условия (27). Точка альфвеновского резонанса определяется уравнением
q2(xA)=0. Она лежит в области непрозрачности левее
точки поворота: xA<xR.
Вблизи этих точек будем использовать разложения
xA, в которой Ω 2A ( xA ) = 0. Вдали от этих особых
точек для производной применима оценка d/dx~1/l,
и тогда оба неравенства (22) сводятся к условию
(21). Нетрудно убедиться, что выполнение условия
(21) вдали от особых точек обеспечивает тот факт,
что размеры окрестностей этих точек, в которых
приближение ВКБ неприменимо, малы по сравнению с масштабом неоднородности. Конкретные оценки размеров этих окрестностей будут даны ниже.
cA ( x) =
55
ω  x − xA 
ω  x − xR 
1 −
 , cA ( x) = 1 −
 . (28)
2lA 
2lR 
kz 
kt 
В.А. Мазур
нологически оно сводится к замене в выражении (25)
для q2 величины ω на ω+iγ, где γ>0 – бесконечно
малая. Это приводит к сдвигу особой точки в нижнюю полуплоскость: xA→xA–iε, где ε=2γ/ω, т. е. к
обходу ее сверху. При таком обходе из (33) получается (34).
При | x − xA | << ∆ A два главных члена разложения выражений (33) и (34) совпадают:
Далее будем считать, что lA~lR~l – величины одного порядка. Используя (28), условия применимости приближения ВКБ (22) можно привести к виду
| x − xA | >> ∆ A ≡ 1/ k y , | x − xR | >> ∆ R ≡ (lR / kt2 )1/ 3 . (29)
Размеры окрестностей особых точек, в которых
неприменимо приближение ВКБ ∆A и ∆R, малы по
сравнению с соответствующими масштабами неоднородности lA и lR, если
k y l >> 1.
 ( x − xA ) 2
x − xA 
% x) = ς A 1 +
(35)
ς(
ln
.
2
2∆ A
∆A 

% x) в точке x=xA конечна: ς(
% xA ) = ς A .
Функция ς(
Но из соотношений (7)–(9) видно, что величины
ξ% y , E% x , B% y имеют особенность вида (x–xA)–1, а вели-
(30)
Условие (30) можно считать конкретизацией общего соотношения (23) на области, удовлетворяющие неравенствам (29). В реальной магнитосфере, с
учетом того, что ky=m/a и l~a, где a – характерный радиус магнитной оболочки, на которой рассматривается
возмущение, неравенство (30) сводится к условию
m >> 1.
(31)
чины ξ% x , E% y , B% x – особенность вида ln(x–xA).
В области применимости метода ВКБ между
точками xA и xR, т. е. при x − xA >> ∆ A , xR − x >> ∆ R ,
Напомним, что мы считаем характерными значения m=3÷5, так что с некоторой натяжкой условие
(31) можно считать выполненным.
При выполнении условия (30), (31) магнитосферная область x<xM разбивается на пять интервалов: три интервала, в которых применимо приближение ВКБ, лежащих в областях x<xA, xA<x<xR, x>xR,
и две окрестности особых точек – | x − xA | ≤ ∆ A и
| x − xR | ≤ ∆ R . В первых трех интервалах можно написать соответствующие выражения приближения
ВКБ. В разделяющих их окрестностях особых точек,
используя разложения (28), основное уравнение (26)
можно упростить и привести к определенным стандартным формам: вблизи точки xA – к уравнению
Бесселя, а вблизи точки xR – к уравнению Эйри, для
которых также можно написать соответствующие
решения. После этого, учитывая граничное условие
при x→–∞ и используя метод сращивания асимптотических разложений, можно все произвольные
константы, задающие решения в разных интервалах,
выразить через одну константу – амплитуду волны в
магнитосфере. Не приводя детали достаточно очевидных выкладок, кратко опишем получающееся
таким образом решение.
В первой области применимости приближения
ВКБ x<xA, | x − xA | >> ∆ A ,
x

(36)
exp  ∫ | k x ( x ′)dx ′  .
x

 A

Это выражение представляет собой растущую
экспоненту. Удерживать на ее фоне падающую экспоненту, т. е. второе линейно независимое решение
в выражении (20), не позволяет точность метода
ВКБ. Выражение (36) можно переписать в виде:
% x ) = CM
ς(
q 2 ( x)
k x ( x)
1/ 2
1/ 2
x

exp  ∫ | k x ( x ′) | dx ′  ,
x

 R

(37)
(38)
xR
Γ AR =
∫ | k ( x) | dx.
(39)
x
xA
Экспоненциальный множитель exp ГAR характеризует возрастание амплитуды колебания от точки
xA до точки xR. Фактор ГAR при выполнении условия
(30) велик: ГAR>>1. Интеграл (39) можно вычислить
аналитически, если предположить, что второе из
выражений (28) справедливо вплоть до точки x=xA:
2
2
2 ky
(40)
Γ AR = (k y ∆ R )3 =
k y lR .
3
3 kt2
В общем случае это выражение можно использовать как оценку по порядку величины.
Вблизи точки поворота, при | x − xR | << lR ,
1/ 2
 x − x   xA − x 
% x) = ς A  A
ς(
 K1 
 при x<xA,
 ∆A   ∆A 
q 2 ( x)
k x ( x)
где
CM = −iC exp Γ AR ,
x

exp  ∫ | k x ( x′) | dx′  .
(32)
x

 A

Вблизи точки альфвеновского резонанса, при
| x − xA | << lA ,
q 2 ( x)
% x) = C
ς(
k x ( x)
% x) = −iC
ς(
x −x
ς% ( x) = ς R Ai  R
.
 ∆R 
(41)
Здесь Ai(η) – функция Эйри [Справочник по специальным функциям, 1979]. Точка поворота x=xR,
будучи особой точкой асимптотического разложения
ВКБ, не является особой для самого уравнения (26).
Поэтому решение (41) аналитично в точке xR.
Наконец, в области прозрачности x>xR и при условии применимости приближения ВКБ x–xR>>∆R
имеем
(33)
 x − xA    x − xA 
 x − xA  
% x) = ς A 
ς(
  K1 
 − −iπ I1 
  при
 ∆ A    ∆ A 
 ∆ A  
x>xA.
(34)
Здесь K1(ξ) и I1(ξ) – модифицированные функции
Бесселя [Справочник по специальным функциям,
1979]. Точка x=xA есть точка ветвления решения (33),
(34). Правило ее обхода хорошо известно [Ландау,
1949; Линь, 1958;Арцимович, Сагдеев, 1979]. Феноме-
 q 2 ( x) 
ς% ( x) = 2CM 

 k x ( x) 
56
1/ 2
π

cos  Ψ ( x ) −  .
4

(42)
МГД-резонатор в лобовой части магнитосферы и его возбуждение…
Здесь выражение для фазы Ψ(x) конкретизировано следующим образом:
2
 ω2 − kt2 cAM

2
Ψ M ( ω) = k t l M 

2 2
3
 kt cAM 
x
Ψ ( x) =
∫k
x
( x ′) dx ′.
(43)
Плотность энергии колебаний (15) в магнитосферной области x<xM можно привести к виду
W=
(44)
ς R = 2 π k y ∆1/R 2 CM .
W=
Из приведенных выражений ясно видно, что рассматриваемое колебание в области прозрачности
представляет собой быстрый магнитный звук (БМЗ)
с дисперсионным уравнением
+
а в точке альфвеновского резонанса – альфвеновскую волну с дисперсионным уравнением
Фаза приближения ВКБ (43) зависит как от параметра от частоты колебания: Ψ=Ψ(x, ω). Важную
роль в дальнейшем изложении будет играть значение фазы в точке x=xM, которое является функцией
одной только частоты:
∫
xR ( ω )
W ( x) =
x
(46)
1
.
v g ( x, ω)
∫
W ( x ) dx =
2πωτ g (ω) CM
xR
(51)
B02
2
.
(52)
В дальнейшем изложении важную роль будут
играть решения в магнитосфере, удовлетворяющие
граничному условию на магнитопаузе
dς%
|x = 0.
(53)
dx M
Оно равносильно условию ξ% | = 0, означаю-
(47)
x xM
где
щему, что нормальное к магнитопаузе смещение
плазмы равно нулю. Ниже мы увидим, что в задаче
о собственных колебаниях магнитосферы это условие является нулевым приближением по малому
параметру cSW/cAM для условий сшивки на магнитопаузе. Поэтому соответствующие решения будем
называть собственными модами.
Подстановка в (53) решения (42) дает условие
−1
k c2 k
 ∂k 
v g =  x  = x A = x cA
ω
k
 ∂ω 
есть проекция полной групповой скорости быстрого
магнитного звука cA на ось x.
Функции (46) и (47) определены только для
2
ω2 > kt2 cAM
. Их легко вычислить, если ω2–kt2×
×c2AM<<kt2c2AM. Тогда на всем интервале интегрирования в (46) можно использовать линейное разложение
 x − xM 
cA ( x) = cAM  1 −
,
2lM 

B02
2
xM
WR =
x
M
M
d
dx
k
x
dx
=
(
,
ω
)
,
x
∫
∫
dω xR ( ω )
v ( x, ω)
xR ( ω ) g
2 πω CM
Полная энергия, заключенная в области прозрачности, дается выражением
Подчеркнем, что нижний предел в этом интеграле также является функцией частоты: xR=xR(ω).
Производная τg=dΨM(ω)/dω есть время пробега с
групповой скоростью волны частотой ω от точки
поворота до магнитопаузы. Действительно,
τ g (ω) =
k x (k 2 + k z2 ) 2 
π 
sin  Ψ ( x ) −  
2
q
4 

Усредним это выражение по мелкомасштабным
пространственным осцилляциям, длина волны которых порядка 1/kx<<l. Сохраняя для такого среднего
прежнее обозначение, имеем
ω2 = cA2 k z2 .
Ψ M (ω) ≡ Ψ ( xM , ω) =
2π
2
CM ×
B02
4
2
2
2
π
 q + k y (k + k z )

×
cos 2  Ψ ( x) −  +
2
4
kx q


ω2 = cA2 k 2 , k 2 = k x2 + k y2 + k z2 ,
ω2
− kt2 dx.
2
cA ( x )
4
2
2
2

π  q + k y (k + k z ) 2 k 2 + k z2
| ς% | +
| ς% ′ |2  . (50)

2
4
4
2 B0 
q
q

В области прозрачности, подставляя решение
(42), получаем отсюда
(45)
xM
(49)
1/ 2
Все константы, определяющие амплитуду колебания в различных областях, можно выразить через
одну – любую из них. Выразим их через CM – амплитуду в области прозрачности. Выражение для C
следует из (38). Для констант ςA и ςR имеем
2 kz
exp(−Γ AR )CM ,
π k y lA1/ 2
,
2

l  ω2 − kt2 cAM
τ g ( ω) = 2 M 
 .
2 2
cAM  kt cAM 
xR
ςA = i
3/ 2
Ψ M (ω) = π(n + 1/ 4),
(54)
где n = 0, 1, 2, ... – волновое число по координате x.
Уравнение (54) вполне аналогично условию Бора–
Зоммерфельда в квантовой механике, оно определяет собственные частоты резонатора ω = ωn . Если в
нем использовать выражение для фазы (49), то
(48)
справедливое при xM–x<<lM. Будем считать, что величина lM того же порядка, что lA и lM, т. е. порядка l.
Имеем
57
2/3
 1  3 

1 
1
ωn = kt cAM 1 +   n +  
.
2/3 
4   (kt lM ) 
 2  2 
(55)
В.А. Мазур
Это выражение применимо, если n<<ktlM (тогда
размер локализации моды по оси x мал по сравнению с lM, и на этом размере применимо разложение
(48)). С другой стороны, приближение ВКБ требует
n>>1. Эти условия совместимы, поскольку ktlM>>1.
Описанное выше решение ς% ( x) также зависит от
ω как от параметра, т. е. ς% = ς% ( x, ω). При ω=ωn это
решение представляет собой собственную функцию:
ς n = ς n ( x) ≡ ς% ( x, ωn ).
Вблизи частоты ω=ωn функцию ΨM(ω) можно
представить в виде
Ψ M (ω) = π(n + 1/ 4) + (ω − ωn )τ gn ,
где τ gn = τ g (ωn ). Из этой формулы видно, что характерное расстояние между соседними собственными частотами равно τ −gn1 .
Поглощение в окрестности альфвеновского
резонанса
Общее решение в области прозрачности (23) для
магнитосферы можно записать в виде
Для значений n<<ktlM, когда мода локализована
на размере, много меньшем lM, для нее можно использовать выражение (41) во всей области прозрачности. Действительно, в этом случае можно
считать, что разложения (28) в точке поворота и (48)
на магнитопаузе совпадают, так что lR=lM и
1/ 2
 q 2 ( x) 
% x ) = C M e iπ / 4 
ς(

 k x ( x) 
Сравнение с полученным выше решением (42)
дает для коэффициента отражения значение
RM = e− iπ / 2 . Его модуль RM = 1 , т. е. потоки энергии
Если подставить эти выражения в (41), то при
ω=ωn получаем
Здесь ηn – нули функции
падающей на точку отражения и отраженной от нее
волн равны. Это означает, что эффект поглощения
энергии в окрестности альфвеновского резонанса в
решении (42) утерян. Причина этого заключается в
отмеченном выше обстоятельстве – точность метода
ВКБ не позволяет определить в интервале (xA, xR)
падающую экспоненту на фоне растущей. В результате при переходе в область x>xR теряется экспоненциально малый эффект поглощения.
Учтем его следующим образом. Используя най% xA ) = ς A
денное решение, в частности связь между ς(
и CM как нулевое приближение, вычислим мощность,
диссипируемую в окрестности альфвеновского резонанса. После этого с помощью уравнения сохранения энергии (14) найдем поправку к коэффициенту отражения RM.
Вблизи точки альфвеновского резонанса выражение для плотности энергии (50) принимает вид
(56)
Ai ′(−η),
т. е.
Ai ′(−ηn ) = 0. Нули ηn (n=0, 1, 2, …) представляют
собой возрастающую последовательность положительных чисел, характерное значение которых порядка единицы. При n>>1
3 
1 
ηn =   n +  
4 
2 
2/3
.
Собственные частоты выражаются через эти
числа:
 1 η

n
ωn = kt cAM 1 +
.
1/ 3
 2 ( kt lM ) 
(58)
× {exp [ −i Ψ ( x) ] + RM exp [i Ψ ( x ) ]} .
 ω

− 1 .
xR ≡ xR (ω) = xM − 2lM 
k
c
 t AM


x −x
ς n ( x) = ς R Ai  − ηn + M
.
∆R 

×
(57)
При n>>1 выражение (57) для частоты переходит
в (55), а выражение (56) для собственной функции в
области прозрачности – в (42). Подчеркнем, что в
практическом плане важны моды с n~1, которые
описываются выражениями (56) и (57). Графики
первых трех собственных мод ςn(x) изображены на
рис. 3.
W ( x) =
2 2
π k y kz
2
ςA .
B02 q 4 ( x)
(59)
Это выражение в точке x=xA обращается в бесконечность. Для его регуляризации необходимо учесть
эффекты, выходящие за рамки идеальной МГД.
В условиях земной магнитосферы важнейшим из
таких эффектов является затухание колебаний в
ионосфере, обусловленное ее конечной проводимостью. Вблизи точки xA колебание представляет собой альфвеновскую волну. Декремент затухания
альфвеновских волн в ионосфере вычислялся во
многих работах [Newton, et al., 1978; Southwood,
Hughes, 1983]. Применительно к нашей модели среды такой расчет дан в приложении. В дальнейшем
для нас будет важен только тот факт, что этот декремент γA много меньше частоты волны ω. Характерные значения для стоячих альфвеновских волн в
магнитосфере γA/ω~10–1÷10–3 [Newton, et al., 1978].
Согласно результатам приложения, конечная
проводимость ионосферы приводит в уравнении (26)
к замене q 2 → q% 2 = q 2 + ik z2 ε P , что равносильно замене ω → ω + iγ A , где εP связано с γA соотношением
Рис. 3. Графики собственных функций магнитосферного резонатора ςn(x).
58
МГД-резонатор в лобовой части магнитосферы и его возбуждение…
εP=2γA/ω. Такая замена согласуется с использованным выше правилом обхода особой точки. Вблизи
точки x=xA, там, где применимо разложение (28),
имеем
RM = TM eiφM .
Для модуля получаем
TM = 1 − 2δ A .
x − xA + iε P lA
q 2 = k z2
.
lA
Фаза в нулевом порядке, как мы видели, имеет
значение φM=–π/2. В следующем порядке положим
Для комплексных значений q2 множитель
q =(q2)2 в знаменателе выражения (59) заменяется на
π
+ α M , α M << 1.
2
Примененный нами метод определяет поправку к
TM, но не позволяет найти поправку αM. Ниже будет
видно, что первая поправка определяет декремент
затухания собственной моды, а вторая – малую поправку к вещественной части частоты. Последняя
же не представляет существенного интереса.
φM = −
4
2
q 2 . Таким образом,
2 2
π k y kz
ςA
W ( x) = 2
B0 q 2 2
2
=
4 δ A lA CM
2
B02  ( x − xA ) 2 + ε 2P lA2 
.
(60)
Здесь использовано соотношение (44) и введен
малый параметр
Решение для солнечного ветра
Для интересующих нас колебаний в магнитосфере имеется область прозрачности. Это означает, что
их частота ω > kt cAM . Поскольку в солнечном ветре
cA << cAM , cS << cAM , то в этой области выполняются неравенства
δ A = (1/ 2) exp ( −2Γ AR ) ,
характеризующий степень проникновения возмущения из области прозрачности в точку альфвеновского резонанса через барьер непрозрачности между
точками xR и xA.
Полная энергия, сосредоточенная в окрестности
альфвеновского резонанса,
∞
WA =
∫ W ( x)dx =
2 π ω δ A CM
2
(61)
.
B02 γ A
−∞
ω >> kt cA , ω >> k z cA , ω >> k z
Q ( x) = 2 γ AW ( x) =
B02
2
ε P lA
( x − xA )
2
+ ε P2 lA2
k x2 =
QA =
∫ Q( x)dx = 2γ
−∞
WA =
A
4 πωδ A CM
B02
∫k
,
(62)
dx′
∫ c ( x′)
xM
(64)
– время пробега от магнито-
S
паузы до точки x с локальной скоростью cS(x). Из
приведенных формул видно, что рассматриваемое
колебание представляет собой в солнечном ветре
быстрый магнитный звук, который в силу неравенства cS>>cA можно считать просто звуком. Из соотношений (63), (64) следует, что kx>>kt, т. е. рассматриваемые волны падают на магнитопаузу почти нормально. Волны, не удовлетворяющие этому условию,
полностью отражаются от магнитопаузы, так как магнитосфера является для них областью непрозрачности.
Общее решение приближения ВКБ (23) в солнечном ветре с учетом соотношений (64) можно
записать в виде
1/ 2
% x) = CW [ρ0 ( x) cS ( x) ] ×
ς(
Используем для S x в области прозрачности
выражение (24). Сравнивая (23) и (58), с учетом
первого из соотношений (25) видим, что
× {exp [ −iωτS ( x) ] + RW exp [iωτS ( x) ]} .
С0 = ( 4π / B02 ) CM . Вместе с (62) это дает
2
2
(63)
( x′)dx′ = ωτS ( x).
x
S x ( x) + QA = 0.
RM
x
Здесь τS ( x) =
Замечательная особенность этих выражений состоит в том, что они не зависят от γA и остаются конечными в пределе γA→0.
Для стационарных колебаний с вещественной
частотой в уравнении (14), очевидно, ∂W / ∂t = 0.
Если два оставшиеся члена уравнения проинтегрировать по x от –∞ до некоторого значения x внутри
области прозрачности и учесть, что S x (−∞) = 0, то
получим
2
.
xM
2
.
ρ0 Ω A2
ω2
,
= ωρ0 ( x)cS ( x ),
kx
cS2 ( x)
x
а полная диссипируемая мощность – выражением
∞
cA2 + cS2
Из них следует, что в солнечном ветре нет точек
поворота, точек альфвеновского резонанса и точек
резонанса медленного магнитного звука, т. е. приближение ВКБ применимо в этой области везде.
В силу неравенств (63) имеем
Интеграл здесь формально распространен на
бесконечные пределы, так как при использовании
выражения (60) он фактически набирается в малом
интервале шириной порядка εPlA вблизи точки x=xA.
Плотность диссипируемой мощности дается выражением
4 ω δ A CM
cA cS
(65)
Величину CW = C0 ω будем называть амплитудой падающей волны. Поток энергии падающей
волны есть
= 1 − 4δ A .
2
Представим коэффициент отражения в виде
S W = CW / 4.
59
В.А. Мазур
При использовании этих обозначений формула
(24) для полного потока энергии принимает вид
(
S x = − S W 1 − RW
2
го интереса. Эта поправка далее будет игнорироваться. По этой же причине будем игнорировать
% .С
малое слагаемое αM/2 в выражении (69) для Ψ
M
учетом сказанного имеем:
).
Условия сшивки на магнитопаузе (12) при подстановке в них решения в магнитосфере (58) и решения в солнечном ветре (65) принимают вид
(
ω = ωn − iν n , ν n = ν Wn + ν An , ν Wn =
= δ W / τ gn , ν An = δ A / τ gn .
Физические причины затухания собственных
мод становятся особенно ясны при рассмотрении
энергетических соотношений. Амплитуда убегающей волны определяется уравнением (66):
)
CM eiπ / 4 e −iΨ M + RM eiΨ M =
1/ 2
 B2 
= CW  0 
 4 πω 
−iΨ M
1
(1 + RW ) ,
δ1/W2
(66)
1/ 2
 πω 
Cout = 4  2n 
 B0 
iΨM
1 − RW
e
− RM e
.
= δW
−iΨ M
iΨ M
1 + RW
e
+ RM e
(67)
Sout =
2
ρ c
ω2 − k z2 cAM
δ W = M AM
.
ρ W cSW ω ( ω2 − k 2 c 2 )1/ 2
t AM
1
Cout
4
2
=
4πωnν Wn τ gn
B02
2
CM .
(71)
Сопоставляя формулы (52), (62) и (71), с учетом
обозначений (70) приходим к естественным с физической точки зрения соотношениям
По порядку величины, согласно соотношению
(1), δ W ~ cSW / cAM << 1.
Малая величина параметра δW имеет ключевое значение для развиваемой в настоящей работе теории.
Sout = 2ν WnWR , QA = 2ν AnWR .
Из них следует, что
∂WR
= − ( Sout + QA ) = −2ν nWR .
∂t
Собственные моды магнитосферного резонатора
Для собственных мод магнитосферного резонатора граничное условие при x→∞ есть условие убегания, т. е. отсутствие падающей на магнитопаузу
волны. Из общей формулы (65) такое решение можно получить с помощью предельного перехода
CW→0 при условии CW RW ≡ Cout = const. Здесь Cout
– амплитуда убегающей волны. В этом пределе
RW → ∞ и соотношение (67) приводится к виду
Энергия собственных мод резонатора теряется
по двум каналам – убегает в солнечный ветер через
не полностью отражающую магнитопаузу и диссипирует в окрестности альфвеновского резонанса.
Возбуждение магнитосферного резонатора падающей волной
Для полного решения задачи о падении на магнитосферу монохроматической волны из солнечного
ветра нам осталось определить коэффициент отражения RW и выразить амплитуду волны в магнитосфере
CM через амплитуду волны в солнечном ветре CW.
Коэффициент отражения определяется из уравнения (67). В первом порядке по малым параметрам
δA и δW имеем
% ( ω) + i ( δ − δ )
tg Ψ
M
A
W
(72)
RW = −
.
% ( ω) + i ( δ + δ )
tgΨ
A
W
(68)
Здесь
% (ω) = Ψ (ω) + φ M = Ψ (ω) − π + α M .
Ψ
M
M
M
2
4 2
δ1/W2 CM .
Поток ее энергии
Здесь обозначено
% (ω) + i (δ + δ )  = 0.
tg  Ψ
M
W
A 
(70)
(69)
В нулевом порядке по малым параметрам δW, δA,
αM уравнение (68) сводится к соотношению (54). Это
означает, что описанные выше собственные моды
действительно в главном порядке по малым параметрам представляют собой собственные колебания
магнитосферы. Тот факт, что в этом главном порядке сшивка с решением в солнечном ветре дает условие (53), или, что то же самое, ξ x |xM = 0, является
В соответствии с аргументами предыдущего раздела здесь и далее считаем, что
% (ω) = Ψ (ω) − iπ / 4.
Ψ
M
M
Из (72) следует:
вполне естественным. В пределе ρW/ρM→∞ можно
считать, что плазма солнечного ветра обладает
очень большой (бесконечной) инерцией и поэтому
граница с ней неподвижна.
В следующем порядке по указанным малым параметрам уравнение (68) дает
RW
2
=
% (ω) + ( δ − δ )2
tg 2 Ψ
M
A
W
.
2
2 %
tg Ψ M (ω) + ( δ A + δ W )
(73)
Величина
(
QW = S W 1 − RW
α
( ω − ωn ) τ gn + M = −i ( δ W + δA ) .
2
2
)
представляет собой поток энергии из солнечного
ветра в магнитосферу. Ее можно трактовать как
мощность накачки магнитосферного резонатора. Из
(73) имеем
Малая величина αM приводит к малой поправке к
реальной части частоты, не представляющей особо60
МГД-резонатор в лобовой части магнитосферы и его возбуждение…
QW = S W
4δ A δ W
% (ω) + ( δ + δ )
tg Ψ
M
A
W
2
2
полной энергии в окрестности альфвеновского резонанса:
2δ A δ W / γ A
,
(79)
WA = S W 2
2
%
tg Ψ (ω) + ( δ + δ )
(74)
.
Амплитуда волны в резонаторе CM определяется
уравнением (66). Подставляя в него выражение (72),
имеем
1/ 2
 B2 
CM = iCW  0 
 4πω 
M
δ1/W2
.
% (ω) + i ( δ + δ )
tgΨ
M
A
W
QA = 2γ AWA = S W
Отсюда
CM
2
= CW
B02
δW
.
% (ω) + ( δ + δ )2
4πω tg 2 Ψ
M
A
W
2δ W
1
,
2
%
v
(
tg Ψ M (ω) + ( δ A + δ W ) g x, ω)
W ( x) = S W
2
W ( x) = S W
2δ A δ W / γ A
% (ω) + ( δ + δ )2
tg Ψ
M
A
W
2
ε lA
π ( x − xA ) 2 + ε 2 lA2 
% (ω) + ( δ + δ )2
tg Ψ
M
A
W
.
(80)
(75)
×
ν AWR = γ AWA ,
Q ( x) = 2γ AW ( x) =
4δ A δ W
ε P lA
% (ω) + ( δ + δ )2 π ( x − xA ) 2 + ε 2P lA2 
tg 2 Ψ
M
A
W


фоне. Это означает, что все они зависят от частоты
примерно одинаково. Рассмотрим для примера одну
из них – полную энергию резонатора WR . Точнее,
мы рассмотрим отношение
. (77)
Проинтегрировав выражения (75)–(77) по координате x, получаем выражения для полной энергии
магнитосферного резонатора:
2δ W τ g
% (ω) + ( δ + δ )
tg 2 Ψ
M
A
W
2
,
(82)
устанавливающую связь между полными энергиями,
заключенными в резонаторе и в окрестности альфвеновского резонанса.
Все функции, фигурирующие в (81), имеют одинаковые резонансные знаменатели, которые и определяют их главную зависимость от частоты ω. Более
слабая зависимость от ω остальных множителей –
δ A , δ W , τ g , γ A – не играет существенной роли на этом
(76)
.
(81)
Они означают, что мощность, поступающая из
солнечного ветра в магнитосферный резонатор, равна мощности, теряемой резонатором в результате
просачивания энергии в область альфвеновского
резонанса, которая, в свою очередь, равна мощности,
диссипируемой в окрестности альфвеновского резонанса. Таким образом, вся энергия, проникающая в
резонатор из солнечного ветра, в конечном счете
диссипирует в окрестности альфвеновского резонанса. Особо отметим часть соотношения (95)
Схематический график зависимости W от x приведен на рис. 4. Плотность мощности, диссипируемой в окрестности альфвеновского резонанса, дается выражением
WR = S W
4δ A δ W
2
QW = 2ν AWR = 2γ AWA = QA .
и в окрестности альфвеновского резонанса:
= SW
W
Из соотношений (74) и (78)–(80) следуют равенства
2
Используя это выражение в формулах (51) и (60),
получаем выражения для плотности энергии в области прозрачности:
×
A
и полной мощности, диссипируемой в этой окрестности:
tR = WR / S W ,
(83)
которое не зависит от амплитуды падающей волны
и является характеристикой резонатора как такового.
Величина tR имеет размерность времени. Ее
можно трактовать как время, необходимое для того,
чтобы падающий поток, если бы он полностью проникал в резонатор, накачал его до необходимого
равновесного уровня. Чем больше tR, тем выше уровень накачки резонатора. Если бы резонатор в магнитосфере отсутствовал (т. е. не было скачка плотности на магнитопаузе, что соответствует значению
δW~1), то характерным значением для величины tR
было бы время пробега волны по магнитосфере, т. е.
tR~τg. Наличие резонатора существенно меняет ситуацию, приводя к сильной зависимости величины tR от
частоты ω.
Главное в этой зависимости – наличие острых
максимумов на собственных частотах ω=ωn. При
ω − ω n << τ −gn1 имеем
(78)
tR (ω) =
2ν Wn
.
(ω − ωn ) 2 + (ν An + ν Wn ) 2
(84)
Ширина резонансного пика равна полному декременту затухания собственной моды ν n = ν An + ν Wn .
Рис. 4. График плотности энергии колебаний W ( x) в магнитосфере.
61
В.А. Мазур
Высота пика
tR (ωn ) =
и узкие минимумы в точках ω=ωn, где
2ν Wn
.
(ν An + ν Wn ) 2
2
RW (ωn ) =
(ν An −ν Wn ) 2
.
(ν An + ν Wn ) 2
Если ν An ≤ ν Wn , то tR (ωn ) ~ 1/ ν n >> τ gn . Таким
Отметим любопытный случай ν An = ν Wn . Тогда
образом, на резонансных частотах наличие резонатора приводит к значительному увеличению уровня
возмущения в магнитосфере.
Функция tR(ω) имеет минимумы в точках ω = ωn ,
% (ω) = ∞ . Легко увидеть, что часв которых tg 2 Ψ
RW (ωn ) = 0 , т. е. вся падающая на магнитопаузу
2
энергия проникает в магнитосферу, – имеет место
своего рода «просветление» магнитопаузы. Впрочем,
это утверждение является лишь приближенным.
При учете следующих членов разложения по малому параметру δA=δW имеем в этом случае
M
тоты ωn чередуются с частотами ωn. Обозначим
τgn = τ g (ωn ), ν Wn = δ W / ωn , по порядку величины τgn ~ τ gn , ν Wn ~ ν Wn .
Вблизи
минимума,
2
RW (ωn ) = δ A2 / 4 .
при
Заключение
Основные результаты работы можно сформулировать следующим образом.
1. Представлена одномерно-неоднородная модель лобовой части магнитосферы и прилегающей
застойной области солнечного ветра (magnetosheath),
в рамках которой решается задача о падении на магнитосферу гидромагнитной волны из солнечного
ветра. Одномерная неоднородность модели по оси x
отражает неоднородность магнитосферы поперек
магнитных оболочек, в том числе и скачок параметров на магнитопаузе. По координате z плазма в модели ограничена проводящими плоскостями, изображающими магнитосопряженные области ионосферы. Ось y моделирует азимутальную координату
магнитосферы, ее цикличность учитывается периодическими граничными условиями. Магнитное поле
в модели направлено по оси z. Зависимости скорости Альфвена cA и скорости звука cS от координаты
x передают их изменение поперек магнитных оболочек в экваториальной плоскости.
2. В рамках описанной модели и в приближении
идеальной МГД установлено существование в магнитосфере резонатора для быстрого магнитного звука. По координате x волна запирается между магнитопаузой и точкой поворота во внутренней магнитосфере, по координате z – между проводящими ионосферными плоскостями, по координате y условие
периодичности приводит к квантованию волнового
вектора ky. В результате собственные моды резонатора определяются тремя волновыми числами –
волновым числом n по координате x (оно равно числу узлов собственной моды), азимутальным волновым
числом m и продольным (вдоль магнитного поля) волновым числом N. Таким образом, собственные частоты
ω=ωnmN. В рамках приближения ВКБ найдены выражения для собственных частот и собственных мод резонатора. Для основной моды n=m=N=1 и приведенных в тексте характерных значений параметров среды
имеем f=ω/2π≈ktcAM/2π≈5 мГц. Затухание собственных
мод обусловлено двумя причинами: слабым убеганием волны через не полностью отражающую магнитопаузу в солнечный ветер и поглощением в точке
альфвеновского резонанса. Малость последнего обусловлена широким барьером непрозрачности между
резонатором и точкой альфвеновского резонанса.
3. Проникновение в магнитосферу гидромагнитной волны из солнечного ветра и уровень накач-
ω − ω n << τgn , имеем
2
tR (ω) = 2ν Wn τgn2 1 + ( ω − ωn ) τgn2  .


(85)
Значение в минимуме tR (ωn ) = 2ν Wn τgn2 много
меньше значения в максимуме:
2
tR (ωn )
~ (ν n τ gn ) << 1 .
tR (ωn )
Более того, оно много меньше значения, характерного в отсутствие резонатора: tR (ωn ) << τ gn , т. е.
существование резонатора в областях между его
собственными частотами уменьшает уровень возбуждения магнитосферы.
Из (85) видно, что характерная ширина минимума τgn−1 по порядку величины совпадает с расстоянием между собственными частотами. Таким образом,
функция tR(ω) представляет собой набор острых
максимумов на собственных частотах с широкими
минимумами между ними. Ее схематический график изображен на рис. 5.
2
Величина RW как функция ω, наоборот, имеет
широкие максимумы в точках ω = ωn , где ее значение близко к единице:
2
RW (ωn ) = 1 − 4δ A δ W ,
Рис. 5. График уровня накачки магнитосферного резонатора tR(ω).
62
МГД-резонатор в лобовой части магнитосферы и его возбуждение…
ки магнитосферного резонатора носят ярко выраженный резонансный характер. Энергия колебаний в резонаторе, рассматриваемая как функция частоты падающей волны, имеет острые и высокие пики на собственных частотах резонатора. Энергия, локализованная в
окрестности альфвеновского резонанса, пропорциональна энергии резонатора: WA = (ν An / γ A ) WR . Декремент затухания альфвеновских волн γA довольно
мал, γA/ω~10–1÷10–3, так что отношение vAn/γA обычно можно считать большим. Это означает, что энергия колебаний внутри магнитосферы, в основном,
сосредоточена вблизи точки альфвеновского резонанса на малом пространственном масштабе
∆x~(γA/ω)lA. Но роль магнитосферного резонатора и
в этом случае велика – его существование определяет
резонансную зависимость уровня колебаний от частоты, а его собственные частоты – места расположения
точек альфвеновского резонанса: xAn=xA(ωn).
После изложения результатов работы обсудим
вопрос, какие из этих результатов можно считать
свойствами колебаний в реальной магнитосфере, а
какие есть свойства слишком простой модели среды.
Во-первых, обсудим вывод о существовании резонатора. Представляется, что одномерная неоднородность модели правильно описывает эффект запирания волны поперек магнитных оболочек. В направлении вдоль магнитного поля колебания в нашей модели ограничены проводящими ионосферами.
В реальной магнитосфере значительное возрастание
скорости Альфвена вдоль силовой линии от экватора
к ионосфере приводит к тому, что магнитозвуковые
волны не достигают ионосферы, а отражаются от точек (поверхности) поворота и в этом направлении. Но
этот факт не отменяет запирания колебания вдоль силовой линии и в рамках нашей модели может быть
учтен эффективным уменьшением параметра L.
Более неопределенная ситуация с координатой y.
Если азимутальная неоднородность магнитосферы
такова, что при движении от лобовой точки к флангам скорость Альфвена нарастает, то возможно запирание волны и по азимуту. В нашей модели этот
факт можно учесть ограничением волнового числа
m ≥ 2π / ∆φ , где ∆φ – азимутальный размер резонатора. При ∆φ~60÷120° имеем m ≥ 3 ÷ 6. Если же
азимутальная неоднородность носит противоположный характер или просто недостаточна для запирания волны, то резонатор представляет собой тороидальную область, охватывающую Землю, что, собственно, и подразумевается в нашей модели. В таком случае возможны меньшие значения m, но реально имеющаяся азимутальная неоднородность
делает самые малые значения m=0, 1 вряд ли приемлемыми. Таким образом, вывод о существовании
МГД-резонатора во внешней магнитосфере представляется вполне правдоподобным, что находит
свое подтверждение в расчетах, выполненных в более реалистичных моделях среды [Lee, Lysak, 1991;
Leonovich, Mazur, 2001].
Во-вторых, обсудим реалистичность рассмотренных механизмов затухания собственных мод.
Плазма солнечного ветра и магнитосферы является
в высшей степени бесстолкновительной, поэтому
63
пренебрежение столкновительной диссипацией вполне
оправдано. Бесстолкновительное затухание Ландау
быстрого магнитного звука и альфвеновской волны
в холодной плазме также чрезвычайно мало [Ахиезер и др., 1974]. Остаются два учтенных нами механизма. Поглощение в точке альфвеновского резонанса имеет неустранимый характер. Сколь бы малы
ни были эффекты диссипации, это поглощение
обеспечивается одним (наибольшим) из них. В окрестности альфвеновского резонанса колебание
представляет собой альфвеновскую волну, которая,
в отличие от быстрого магнитного звука, даже в
продольно неоднородной плазме достигает ионосферы [Dungey, 1954]. Поэтому наиболее важным
диссипативным эффектом оказывается омическое
поглощение в ионосфере. Затухание собственных
мод, обусловленное их убеганием в солнечный ветер, имеет бездиссипативную природу. Его существование не вызывает сомнения, но возникает вопрос:
не является ли оно слишком большим? Представляется, что наша одномерная модель вполне разумно
описывает явление частичного отражения волны от
магнитопаузы и дает правильную оценку параметра
δW, характеризующего «проницаемость» магнитопаузы для магнитозвуковых волн.
Третий вопрос, который следует обсудить, – вывод о резонансном характере возбуждения колебаний в магнитосфере. Представляется, что этот вывод
является очень общим и имеет место в любой модели среды при условии существования резонатора.
Аналогичными свойствами обладают резонаторы
самой разной природы – от электротехнических
[Каценеленбаум, 1966] до атомных и ядерных объектов [Ландау, Лифшиц, 2004]. Во всех этих случаях
зависимость уровня возбуждения резонатора от частоты падающей волны характеризуется наличием
максимумов на собственных частотах резонатора.
Высота этих максимумов определяется потоком падающей энергии и «проницаемостью» границы резонатора, т. е. тем же параметром δW, что определяет
излучательную часть затухания собственных мод.
Эти общие свойства целиком находят свое отражение в результатах нашей работы.
Последний вопрос, который мы обсудим, связан
с азимутальной асимметрией падения гидромагнитной волны на магнитопаузу. Застойная область солнечного ветра примыкает только к лобовой части
магнитосферы. На ее флангах скорость течения солнечного ветра велика, и ее наличие нельзя игнорировать. В этих областях важную роль играет неустойчивость Кельвина–Гельмгольца. С ночной стороны солнечный ветер вообще отсутствует. Столь
же асимметрично и падение волны на магнитопаузу.
Если это падение приблизительно нормально в ее
лобовой части, то, конечно, далеко от нормального в
других ее частях. Поэтому мы ограничили рассмотрение лобовой частью магнитосферы. Если резонатор в магнитосфере также ограничен по азимуту
лобовой ее частью, то наше рассмотрение самосогласованно (с учетом ограничения снизу волнового
числа m). Если же резонатор имеет тороидальную
форму, охватывая все значения азимутальной координаты, то проникновение в него энергии будет
В.А. Мазур
происходить только через своего рода «окно» в лобовой части. Соответственно, ниже будет уровень
накачки резонатора, что можно учесть пропорциональным уменьшением параметра δW (но только в
отношении проникновения волны извне, для описания убегания волны следует сохранить прежнее значение этого параметра).
Проведенный критический анализ применимости
использованной в работе одномерной модели среды
дает основания утверждать, что она адекватно отражает основные свойства рассматриваемого явления,
дает правильные качественные выводы и правильные оценки важнейших параметров, описывающих
это явление.
в ионосферных слоях с учетом того, что Ez=0 и
εˆ ⊥ = i (4π / ω) σˆ ⊥ , принимает вид
r
r
4πω r
∆ E⊥ − ∇ ⊥ (∇ ⊥ E⊥ ) = −i 2 σˆ ⊥ E⊥ .
c
(П2)
Проинтегрируем его по z в интервале (–∆, 0).
Пренебрегая изменением электрических полей
внутри этого интервала, имеем
r
r
r
∂E⊥
∂E⊥
|z = 0 −
|z =−∆ = −iKˆ 0 E⊥ |z = 0 .
(П3)
∂z
∂z
)
Тензор K 0 имеет вид
 KP KH 
Kˆ 0 = 
,
 −KH KP 
Приложение. Учет конечной проводимости
ионосферы
Конечная проводимость ионосферы приводит к
изменению граничных условий на ионосферных
плоскостях z=0 и z=L. Для вывода этих граничных
условий конкретизируем модель среды вне магнитосферного слоя, т. е. при z<0 и z>L. Будем считать,
что ионосферные слои имеют малую, но конечную
толщину ∆ и занимают области (–∆, 0) и (L, L+∆).
Поперечные компоненты тензора проводимости в
этих областях
где
0
K P, H =
4πω
Σ P, H , Σ P, H = ∫ σ P, H dz ,
c2
−∆
Σ P, H – интегральные проводимости ионосферы. Ин-
декс «0» у тензора K̂ 0 означает, что он относится к
ионосфере z=0. Величины KP и KH имеют размерность волнового вектора. Будем предполагать, что
проводимость ионосферы столь велика, что
K P, H >> k⊥ . Тогда, учитывая (П1), приводим (П3) к
 σP σH 
σ̂ ⊥ = 
,
 −σ H σ P 
виду
где σP – педерсеновская, σH – холловская проводимости. Как и ранее, будем полагать, что продольная
проводимость σ|| = ∞ . Еще далее от магнитосферно-
r
r
∂E⊥
|z = 0 = −iKˆ 0 E⊥ |z = 0 .
∂z
r
Разрешая это соотношение относительно E⊥ |z = 0 ,
имеем
r
r
∂
E
1
−
⊥
|z = 0 .
(П4)
E⊥ |z = 0 = i Kˆ 0
∂z
го слоя, при z<–∆ и z>L+∆, находятся области, моделирующие нейтральную атмосферу, в которых
σ̂ ⊥ = 0, σ|| = 0 .
Мы дадим вывод граничных условий в предположении, что проводимость ионосферы достаточно
велика, чтобы в главном порядке ее можно было
считать идеально проводящей, но ее толщина настолько мала, что возмущенное электрическое поле
в пределах ионосферного слоя мало меняется по
координате z (однако производная от поля может
измениться сильно).
Атмосфера является однородной средой. Для
возмущений вида плоских волн дисперсионное
уравнение в ней имеет вид k 2 ≡ k x2 + k y2 + k z2 = 0 . Это
Очевидно,
Kˆ 0−1 =
1
K + K H2
2
P
 KP − KH 

.
 KH KP 
Совершенно аналогично получаем граничное условие для сопряженной ионосферы:
r
r
∂E
E⊥ |z = L = −i Kˆ L−1 ⊥ |z = L .
(П5)
∂z
В пределе идеально проводящей ионосферы, т. е.
при KP, H→∞, правые части в (П4) и (П5) исчезают.
Поэтому далее мы будем их рассматривать как малые поправки к приближению идеальной ионосферы.
Выясним, как эти поправки изменяют основное
уравнение колебаний в магнитосфере (26). Действуя
в духе теории возмущений, положим
r
r
r
E⊥ ( x, z ) = E% ⊥ ( x) U N ( z ) + e⊥ ( x, z ) ,
(П6)
r
где e⊥ ( x, z ) – поправка, обусловленная конечной
проводимостью ионосферы. Без умаления общности
можно считать, что она как функция z ортогональна
решению нулевого приближения:
есть дисперсионное уравнение для электромагнитной волны в вакууме, в котором пренебрегается током смещения, т. е. рассматриваются столь низкочастотные колебания, что можно полагать ω/c→0.
Таким образом, k z2 = − k⊥2 , где k⊥2 = k x2 + k y2 . Возмущение, ограниченное при z→–∞, имеет зависимость
от z вида exp(k⊥ z ) . В таком случае на границе ионосферы и атмосферы
r
r
∂E⊥
|z =−∆ = k⊥ E⊥ |z =−∆ .
(П1)
∂z
Общее уравнение монохроматических колебаний
r ω2 r
rot rot E = 2 εˆ E
c
L
∫U
0
64
N
r
( z ) e⊥ ( x, z )dz = 0 .
(П7)
МГД-резонатор в лобовой части магнитосферы и его возбуждение…
Здесь
Линеаризация граничных условий по этой поправке дает
r
r
e⊥ |z = 0, L = ±iU N′ |z = 0, L ⋅Kˆ 0,−1L E% ⊥ ( x) .
(П8)
2
q% = q 2 + ik z2 ε P .
Параметр εP по порядку величины
Отметим, что U N′ |z = 0 = k z и U N′ |z = L = (−1) N k z .
r
Уравнения для E⊥ можно также получить из
εP ~
r
r
ω2 r
∆E⊥ + 2 E⊥ − ∇ ⊥ (∇ ⊥ E⊥ ) = 0 .
cA
ω = ω0 − iγ, ω0 = k cA , γ =
Подставим в это уравнение разложение (П6), умножим на UN(z) и проинтегрируем по z от 0 до L. С
учетом условия ортогональности (П7) получим:
r
r
r 
L  d 2 E% ⊥
 2 + ( q 2 − k y2 ) E% ⊥ − ∇ ⊥ (∇ ⊥ E% ⊥ )  +
2  dx


r
L
∂ 2 e ( x, z )
+∫U N ( z)
dz = 0.
∂z 2
0
Работа выполнена при поддержке РФФИ, гранты
№ 07-05-00185 и № 09-02-00082.
СПИСОК ЛИТЕРАТУРЫ
Арцимович Л.А., Сагдеев Р.З. Физика плазмы для физиков. М.: Атомиздат, 1979. 320 с.
Ахиезер А.И., Ахиезер И.А., Половин Р.В. Электродинамика плазмы. М.: Наука, 1974. 720 с.
Гульельми А.В., Троицкая В.А. Геомагнитные пульсации и диагностика магнитосферы. М.: Наука, 1973. 208 с.
Гульельми А.В., Потапов А.С., Д’Коста А. К теории
возбуждения геомагнитных пульсаций типа Рс3 // Исследования по геомагнетизму, аэрономии и физике Солнца.
1976. Вып. 39. С. 27–32.
Каценеленбаум Б.З. Высокочастотная электродинамика. М.: Наука, 1966. 244 с.
Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. М.: Наука, 1982. 620 с.
Ландау Л.Д. О колебаниях электронной плазмы //
ЖЭТФ. 1946. Т. 16. С. 574–593.
Линь Цзя-Цзяо. Теория гидродинамической устойчивости. М.: ИЛ, 1958.194 с.
Ландау Л.Д., Лифшиц Е.М. Квантовая механика (нерелятивистская теория). М.: Физматлит, 2004. 800 с.
Справочник по специальным функциям / Под ред.
М. Абрамовица и И. Стигана М.: Наука, 1979. 830 с.
Baumjohann W., Junginger H., Haerendel G., Bauer O.H.
Resonant Alfven waves excited by sudden impulse // J. Geophys. Res. 1984. V. 89. P. 2765–2769.
Duhau S., Gratton J. Kelvin-Helmholtz instability of anisotropic plasma in a magnetic field // J. Plasma Phys. 1975. V. 13.
P. 451–479.
Dungey J.W. Electrodynamics of the outer atmosphere //
Penn. State Univ. Ionos. Res. Sci. Rept. 1954. 69 p.
Fujita S., Tanaka T. Magnetospheric plasma processes
during a sudden commencement revealed from a global MHD
simulation // Geophysical Monograph. 2006. V. 169. P. 31–50.
Guglielmi A.V. Diagnostics of the magnetosphere and interplanetary medium by means of pulsations // Space Sci. Rev.
1974. V. 16. P. 331–345.
Kivelson M.G., Southwood D.J. Coupling of global magnetospheric MHD eigenmodes to field line resonances // J.
Geophys. Res. 1986. V. 91. P. 4345.
Kivelson M.G., Southwood D.J. Resonant ULF waves: A
new interpretation // Geophys. Res. Lett. 1985. V. 12. P. 49.
Lysak R.L., Lee D.H. The response of the dipole magnetosphere to pressure pulses // Geophys. Res. Lett. 1992.
V. 19. P. 937–940.
)
y
ε̂ ⊥ =
2 ˆ −1 ˆ −1  ε P − ε H 
K0 + K L = 
,
L
 εH εP 
(
)
где
εP = ∑
2KP
(K
2
P
+K
2
H
)L
, εH = ∑
2KH
(K
2
P
+ K H2 ) L
,
а знак суммы означает суммирование по сопряженным ионосферам. С помощью этих обозначений
уравнение (П9) можно переписать в виде
r
r
(П10)
( q 2 + ik 2ε) ) E% = −∇% (∇% E% ) .
z
⊥
⊥
⊥
⊥
⊥
Из соотношений (5), (7), (8) и (25) нетрудно получить, что
r
% E% = i 4 πω ς% .
∇
(П11)
⊥ ⊥
cB0
Это означает, что правая часть уравнения (П10)
выражается через функцию ς% . Но тогда, рассматривая его как алгебраическое уравнение на компоненr
ты вектора E% , можно выразить их через функцию
⊥
ς% . Подставляя найденные таким образом значения в
соотношение (П11), получим уравнение на эту
функцию:
d 2 ς% (q% 2 )′ dς%
− 2
+ (q% 2 − k y2 )ς% = 0 .
dx 2
q% dx
k z2
ε P ω0 .
2k 2
В области альфвеновского резонанса ω0=kzсА,
γ ≡ γ A = (1 / 2)ε P ω0 .
Интегрируя в последнем члене по частям и используя соотношения (П8), приводим уравнение к виду
r
r
r
d 2 E% ⊥
2
2
% − ∇ (∇ E% ) +
q
k
E
+
−
(
)
y
⊥
⊥
⊥ ⊥
dx 2
r
2k 2
(П9)
+i z Kˆ 0−1 + Kˆ L−1 E% ⊥ = 0.
L
Введем обозначения
% = (−ik , d / dx ) ,
∇
⊥
1
c2
=
.
K P L 4π ω L ΣP
Поскольку K P >> k⊥ и k⊥ L ~ k⊥ l >> 1 , то тем более K P L >> 1 и, следовательно, параметр ε P << 1 .
В однородной плазме из уравнения (П12) следует дисперсионное уравнение
(П2), если учесть, что в магнитосфере тензор σ̂ ⊥
имеет только диагональные компоненты, равные
−i (ω c 2 / 4 π cA2 ) :
(
(П13)
(П12)
65
В.А. Мазур
Leonovich A.S., Mazur V.A. ULF hydromagnetic oscillations with the discrete spectrum as eigenmodes of MHDresonator in the near – Earth part of the plasma sheet // Annales Geophysicae. 2006. V. 24. P. 1639–1648.
Lee D.H., Lysak R.L. ULF wave excitation in the dipole
magnetosphere // J. Geophys. Res. 1991. V. 96. P. 5811–5823.
Leonovich A.S., Mazur V.A. On the spectrum of magnetosonic eigenoscillations of an axisymmetric magnetosphere //
J. Geophys. Res. 2001. V. 106. P. 3919–3928.
Mann I.R., Wright A.N., Cally P.S. Coupling of magnetosperic cavity modes to field line resonances: A study of resonance width // J. Geophys. Res. 1995. V. 104. P. 441.
Newton R.S., Southwood D.J., Hughes W.J. Damping of
geomagnetic pulsation by the ionosphere // Planet. Space Sci.
1978. V. 26. P. 201.
Potapov A.S., Mazur V.A. Pc3 pulsations: From the
source in the upstream region to Alfven resonances in the
magnetosphere. Theory and observations // Geophysical
Monograph. 1994. V. 81. P. 135–145.
Southwood D.J., Hughes W.J. Theory of hydromagnetic waves in the magnetosphere // Space Sci. Rev. 1983.
V. 35. P. 301.
Walker A.D.M. Excitation of field line resonances by
MHD waves originating in the solar wind // J. Geophys. Res.
2002. V 107. doi: 10.1029/2001JA009188
Walker A.D.M. Excitation of magnetohydrodynamic cavities in the magnetosphere // J. Atmos. Solar-Terr. Phys. 1998.
V. 60. P. 1279.
Институт солнечно-земной физики, Иркутск
66
Download