CC BY
7
Кинетическая модель двумерного цилиндрического токового слоя
И.Ю. Васькоа, В.Ю. Попов6
Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра математики. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
Институт космических исследований РАН.
Россия, 117997, Москва, ул. Профсоюзная, д. 84/32.
E-mail: а [email protected], ь [email protected] Статья поступила 28.05.2011, подписана в печать 13.09.2011
Работа посвящена построению класса стационарных решений системы уравнений Власова-Максвелла в виде двумерных цилиндрических токовых слоев с током вдоль азимутального направления / = jg(p,z)eg. Магнитное поле такой системы обладает двумя компонентами В = Bz{p,z)ez + Bp(p,z)ep. Математически задача сводится к нахождению решений нелинейного уравнения в частных производных для функции u{p,z): д2и/дх2 + x^ld2u/dz2 = , где х = р2. Методами теории групп найдено трехпараметрическое семейство точных решений этого уравнения и рассмотрены асимптотики решений при больших значениях р и вблизи оси цилиндрической системы координат (р~0). Обсуждается применение построенного класса решений для описания токовых слоев, наблюдаемых в магнитосферах планет Солнечной системы, магнитные диполи которых лежат в плоскости эклиптики.
Ключевые слова: система уравнений Власова-Максвелла, нелинейные уравнения в частных производных, цилиндрический токовый слой, магнитосферы планет.
УДК: 533.95. PACS: 52.25.Dg.
Введение
Ключевым вопросом физики космической плазмы является вопрос о структуре и динамике токовых слоев, которые наблюдаются в гелиосфере, в частности в магнитосферах планет Солнечной системы и в солнечном ветре. Эти плазменные объекты обладают уникальным свойством накапливать и затем взрывным образом высвобождать энергию магнитного поля. В результате именно токовые слои играют ключевую роль в крупномасштабных магнитных возмущениях (суббурях) в магнитосферах планет, во вспышках на Солнце и ускорении частиц в солнечном ветре (см., например, [1-3]).
Наиболее простая конфигурация токового слоя, рассматриваемая в большинстве моделей токовых слоев, представляет собой плоский токовый слой, характерный для магнитосфер планет, чей магнитный диполь направлен почти перпендикулярно к плоскости эклиптики. Примером является магнитосфера Земли, чей магнитный момент образует небольшой угол (~ 10°) с осью вращения, которая в свою очередь почти перпендикулярна к плоскости эклиптики. Силовые линии земного диполя при взаимодействии с потоком солнечного ветра вытягиваются на ночной стороне магнитосферы в направлении от Земли и образуют так называемый магнитосферный хвост. В центре его находится практически плоский токовый слой, разграничивающий разнонаправленные магнитные поля. Магнитосферу, схожую с земной, имеют также планеты Меркурий, Юпитер и Сатурн. Среди моделей плоских токовых слоев можно отметить широко известную модель Харриса [4] и ее двумерное обобщение [5, 6]. Кроме того, известен целый класс одномерных моделей плоских токовых слоев, построенных с учетом адиабатических инвариантов частиц (см. [7] и цит. лит.).
Более сложной структурой обладают токовые слои,
наблюдаемые в магнитосферах планет, чьи магнитные диполи лежат почти в плоскости эклиптики (примером являются магнитосферы Урана и Нептуна [8], а также астероида Плутона). Исследования спутника «Voya-ger-2» показали, что ось вращения Урана лежит близко к плоскости эклиптики, а магнитный диполь составляет с осью вращения угол ~60°. В результате во время вращения планеты один из полюсов магнитного диполя периодически принимает положение вдоль направления скорости солнечного ветра. В этом случае обтекание солнечным ветром приводит к симметричному вытягиванию силовых линий, входящих в магнитный полюс, на ночную сторону от планеты. А это в свою очередь приводит к формированию в хвосте магнитосферы цилиндрического токового слоя [8], у которого максимум плотности тока сосредоточен в окрестности гигантской цилиндрической поверхности (или поверхности вращения) на ночной стороне магнитосферы, а вдоль оси хвоста формируется магнитное поле с аксиальной симметрией.
Одномерный цилиндрический токовый слой в рамках кинетической теории был рассмотрен в работе [9]. В этой модели плотность тока направлена вдоль азимутального направления / = je(p)eg (здесь и в дальнейшем будем использовать цилиндрические координаты г = (р,в,г), где ось г направлена вдоль оси хвоста), а магнитное поле В = Вг(р)ег, причем на некоторой поверхности р = const магнитная компонента Вг обращается в ноль. Данная модель не учитывает изменения плазменных характеристик в направлении 2, а также радиальную компоненту Вр магнитного поля.
В настоящей работе построен класс двумерных цилиндрических токовых слоев, в которых плотность тока имеет азимутальное направление j = jg(p,z)eg, а магнитное поле имеет две компоненты
В = В,{р,г)е, + Вр(р,г)ер. С математической точки зрения задача заключается в построении решений одного нелинейного уравнения в частных производных. Групповые методы позволяют построить трехпараметри-ческое семейство решений и тем самым получить класс двумерных токовых слоев.
1. Постановка задачи
Рассматриваются двумерные стационарные решения системы уравнений Власова-Максвелла для электрон-ионной бесстолкновительной (сильно разреженной) плазмы. Пусть система обладает аксиальной симметрией, т.е. все величины зависят от (р,г). Электрическое и магнитное поля определяются скалярным потенциалом ф{р,г) и компонентой векторного потенциала Ад(р,г):
У^ + ^-Щ+кхфхАее, ]])^=0,
плотности тока
дг
А А№ =
47Г х -
(ve0)fa(r,v, t) d3v, fa(r,V, t) d3V.
Здесь a = i,e — сорт частиц и заряды q-L = Условие квазинейтральности плазмы
fa(r,V, t) d3V = О,
(1)
q.
(2)
позволяет определить связь между скалярным и векторным потенциалами. Решением первого уравнения системы (1) является произвольная функция интегралов движения частицы в полях —Чф и [V х А$е$]. Из условий стационарности и двумерности задачи следует сохранение двух интегралов движения — полной энергии На и обобщенного момента Рда:
Н„ =
1
2 тп
Рр,с + (Ре'а - у РЛб) + + ЧпФ'
(3)
Ре,а = Ре,а + ~~рАе = тарщ + ^рА0.
Следовательно, стационарную функцию распределения плазмы можно задавать как функцию интегралов движения: 1а = ¡а(На, Ре,а) ■ КоНКрвТНЫЙ ВИД ЭТОЙ фуНКЦИИ определяется видом источников плазмы и различными магнитосферными процессами. Выберем функцию распределения частиц сорта а аналогично работе [7] вида
( ГПа \Ш ( На-ПаРв,п\ = п0 —— ехр <----} , (4)
fa
\2жТа
ехр
Та
где Та — температура частиц, щ — характерное значение концентрации частиц плазмы, Оа — постоянная угловая скорость их вращения. Условие квазинейтральности (2) приводит к уравнению, аналогичному полученному в [7] для одномерной системы:
Ф = Р2
тп
2q(T + Те) \ Т
т-Д]
mett2e Т,
а
Т,
'Т,
. (5)
Введем обозначения О = О; — Ов, Т = Т + Те. Функция распределения (4) задает единственную компоненту
je = qtiopil ехр \ —рА0 +
(qO, О*<jmT? + meT}) л I сТ 2Т3
Р2 ■ (6)
Далее введем безразмерные переменные р» = р/ро, 2, =z/p0 и А* = Ав/А0, где
Ро
с2Т сТ
го, А0 =
X =
_ П2р2(ггцТ2 + теТ2)
щ2щ$12' qQpo, " 2Т3 ■
Второе уравнение системы (1), представляющее собой уравнение для векторного потенциала, принимает вид
д 1 д л» д2А* , г л» 21
ф: * + ж = ехр + Ар*/ ■(7)
Производя повторную замену переменных х = р\, ( = 22* и и = - (р«Л* + Ар2), получим следующее уравнение для функции и:
д2и 1 д2и _ _и
+ % д? "е ■ 1 }
Безразмерные компоненты магнитного поля обозначим через В* и В*. Единицей измерения поля является величина Ао/ро ■ Компоненты магнитного поля выражаются через функцию и следующим образом:
дА* _ 2 ди
D* _
(9)
Безразмерные концентрацию и плотность тока обозначим п* = и jg = у/х соответственно. Единицей измерения концентрации является величина п0, а плотности тока — qrioilpo.
Уравнение силовых линий, лежащих в плоскости (р., 2»), имеет вид р*Л* = const и может быть записано через функцию и: и + Ах = const.
Уравнение (8) является нелинейным уравнением относительно функции и. Далее рассмотрим построение его частных решений.
2. Построение семейства частных решений уравнения (8)
2.1. Сведение уравнения (8) к обыкновенному дифференциальному уравнению второго порядка
Уравнение (8) допускает две однопараметрические группы, инфинитезимальные операторы которых имеют вид [10]
Решение, инвариантное относительно оператора (10), находится из условия F(x, и) = 0, где Х2F = 0, и имеет вид и = 2 In х + tp(0 и £ = (2/х, £ е R+ . Подставляя эту функцию в уравнение (8), получим обыкновенное дифференциальное уравнение для функции (р
m + W + = 2 + (11)
Заменой г/ = + 4) (г/ е (0,1)) уравнение (11)
приводится к виду
(1^2)%эд = 4(2 + е^). (12)
Явное интегрирование уравнения (12) представляется затруднительным. Однако несложно вычислить асимптотики решения данного уравнения в двух предельных случаях г/ и 0 и г/ га 1, что отвечает рассмотрению системы соответственно вдали от оси цилиндрической системы координат и вблизи нее.
Применяя к функции и группу переносов вдоль оси ( (группа с оператором ), получим однопара-метрическое семейство решений
иа = 2\пх + р((( + 2а)2/х), аеШ1.
(13)
2.2. Асимптотика функции <р при г/ % О
Уравнение (12) перепишем в виде <-рт = 4(2 +
+ е-¥!)/(1 -г]2)2. В окрестности точки г] = 0 правая часть уравнения голоморфна, поэтому асимптотика с точностью 0(г]3) имеет вид
íp = ao + alri + 4[ 1 +
r]2 + 0{rf). (14)
Здесь ао,1 " произвольные константы. Значениям г] к, 0 соответствуют значения £гаО, причем г1= \л/1+0(ф12), поэтому (р = а0+щ ,
где ао i — новые произвольные константы и ц = 1 + + е^а°'/2 Ф 0.
2.3. Асимптотика функции <р при г/ % 1
Из уравнения (12) видно, что <рщ > 8/(1 — rj2)2,
поэтому при r¡ -¥ 1 очевидно, что ip-¥ оо. Тогда при
íj га 1 имеем е-ч> <С 1 и главный член асимптотики
функции ip (обозначаемый далее £¿?as) находится из уравнения
d2ipas 8
dr¡2
"
(15)
(l-rfy
Явное интегрирование (15) дает ipas = —2(1 — - (1 -rj))ln(l-т]) + 2rjln(l+rj) + Ci(1 -т]) + С0, где C0ji — произвольные константы. Тогда с точностью до 0((1 — г])2) асимптотика функции ip примет вид <р = -21п(1 -г)) + 2(1 -rj)ln(l -ri) + С0 + Ci(l - rj) + +0((1 —г])2). Значениям г)х, 1 соответствуют значения £ 1, причем г/ = 1 — 2/£ + О ( 1 /Ç2 ), следовательно
4 /1
ip = 2 In £ + С0 + - In
f+o(iin(a), de)
где Coj — новые произвольные константы.
2.4. Асимптотики основных величин при р» и О и р, » 1
Асимптотики основных величин системы (магнитные поля, плотность плазмы и тока) приводятся для функции и0 семейства (13). Асимптотики этих же величин для функций иа могут быть получены подстановкой +2а (г, ^ г, + а).
Рассмотрим решение вдали от оси цилиндрической системы, когда jc> 1 и, следовательно, ^гаО. Для функции и0 вдали от оси, согласно (14), можно записать
ICI С2
u0 = 2\nx + a0 + al^L+fi— + 0(лГ3/2). (17)
у X X
Заметим, что эта же асимптотика будет справедлива и в другом предельном случае: х = const, Чтобы
избежать разрыва производной дщ/дС,, будем рассматривать решение при ( > 0.
Используя формулы (9) для характеристик плазмы и компонент поля, получаем асимптотики
п* = е^р:4 ■
fí = e^pf + ...,
4 _„ 2at . С (18)
В* = -2X - -г + ..., В* = —i- + 4// s, + ... .
Pi ' Pi Pi
Как показывают асимптотики (18), концентрация и плотность тока вдали от оси достаточно быстро убывают. Компонента В* выходит на отрицательную константу, а В* стремится к нулю.
Рассмотрим теперь систему в противоположном предельном случае: лсгаО, ( = const ф 0 и, следовательно, ■цк, 1. Для функции щ вблизи оси, согласно (16), имеем
u0 = 21nC2 + C0+^ln^^ +С1^ + 0(х21пх). (19)
Заметим, что эта же асимптотика будет справедлива в другом предельном случае: х = const, 1 (с заменой в (19) 0(х2\пх) ^0(С2\п()).
Вновь используя формулы (9), получим при х га О
-С„
tí* =
В* =
16 -2À ■
41пр»
le
В* =
16 4
-p*z» +.
(20)
Асимптотики (18) показывают, что плотность плазмы вблизи оси симметрии падает с увеличением г», а плотность тока Ц стремится к нулю. При этом компоненты магнитного поля имеют сингулярности. Существование этих сингулярностей следует из баланса давления, которое в направлении оси 2 имеет вид const = —jgB*. Из данного уравнения видно, что соотношения дп*/дг*Ф 0 и Ц ~ р» приводят к сингулярности В*~р~1. Наличие сингулярности другой компоненты магнитного поля В* объясняется следующим образом. Из формул для асимптотик плотности тока (18), (20) следует, что она всюду регулярна. С другой стороны, jg ~ дВ*/др* - dB*/dz* и второе слагаемое этого выражения вблизи оси ~р«1- Отсюда заключаем, что первое слагаемое этого выражения обязано содержать слагаемое ~ In р».
Отметим, что полученные асимптотики (18), (20) указывают на наличие поверхности, на которой В* = 0. Это следует из того, что и 5*| ^<0.
Кроме того, асимптотика для концентрации п* (20) показывает, что плотность плазмы убывает при 2» -¥ оо при любом фиксированном значении р».
Далее приведены численное решение уравнения (12) и построение двумерных профилей концентрации и плотности тока.
2.5. Профили концентрации и плотности тока
Задача Коши для уравнения (12) с начальными условиями <р(0) = (ро и <р'(0) = (р'0 решалась численно
на множестве г] е [0, 1) методом Рунге-Кутта 4-го порядка в системе Ма1ХаЬ. Из численного решения (12) находились функции иа. Далее при помощи формулы (9) строились графики различных физических величин в плоскости (р*,2*) (на рис. 1-3 оси обозначены соответственно как (р,г)).
Заметим, что построенное семейство решений (13), которое описывает двумерные цилиндрические токовые слои, является трехпараметрическим. Параметрами служит тройка величин {щ,{р'0,а). Выбирая эти параметры различным образом, можно получать различные токовые слои. Приведем графики физических величин для конкретных значений величин (сро,(р'0,а) и укажем на характер изменения графиков при изменении параметров ((^0,(^0, а).
Все физические характеристики будем изображать при «г* > 0, считая, что центр планеты расположен в отрицательной области оси «г*. Для построения силовых линий необходимо задать численное значение параметра Л. Однако форма силовых линий качественно не зависит от Л и далее используется Л = 0.1.
Примем следующие значения параметров ((р0, (р'0, а) = (-3,0, 3). График концентрации п*/п^ах
Рис. 2. Силовые линии магнитного поля для (<Ро,<Ро> а) = (-3,0, 3)
-40
-40
Рис. 1. Концентрация плазмы п*/п^ах (а) и плотность тока ¡е/Цтш (б) для (<р0, а) = (-3,0, 3)
Рис. 3. Плотность тока для (щ^^а)
= (-3,0,5) (а) и (^.^.¿) = (-5,0,3) (б)
плазмы показан на рис. 1 ,а, плотности тока /|//|тах — на рис. 1,6. Концентрация частиц плазмы и плотность тока нормированы на максимальные значения п*тш и /|тах соответственно. Вид силовых линий изображен на рис. 2.
Чтобы показать качественные изменения в профилях физических величин при изменении параметров {(р0,(р'0,а), на рис. 3,а приведен профиль плотности тока для параметров (ср0, (р'0, а) = (^3,0,5).
Видно, что с увеличением параметра а область, в которой плотность тока отлична от нуля, расширяется в направлениях р* иг», причем расширение в направлении р* гораздо сильнее. Кроме того максимум тока смещается в сторону больших значений р*.
Изменение параметра <р0 приводит в основном лишь к смещению максимума плотности тока в сторону больших значений р*. Это демонстрирует рис. 3, б, на котором приведена плотность тока в слое с (<р0,<р'0,а) = (-5,0,3).
Численный анализ показывает, что профили физических величин слабо зависят от параметра (р'0, так что в действительности имеется два параметра, которые определяют структуру токового слоя.
Наличие поверхности, на которой Вг = 0, подтверждается рис. 3. Кроме того, видно, что вдали от оси
в;^ о.
Заключение
В работе получен класс двумерных решений уравнения Власова-Максвелла, описывающих цилиндрические токовые слои, и приведены самосогласованные профили магнитного поля, плотностей плазмы и тока. Показано, что рассматриваемая система обладает тремя свободными параметрами {(ро,(р'0,а), выбирая которые можно получать различные токовые слои. Приведены
примеры таких слоев для некоторых конкретных значений параметров.
Основной особенностью построенного класса кинетических равновесий является сингулярность магнитного поля вблизи оси симметрии системы. Это указывает, что выбор функции распределения (4) не вполне удовлетворителен с физической точки зрения. Тем не менее построенная модель хорошо описывает токовый слой всюду, кроме области, непосредственно примыкающей к оси системы.
Авторы выражают благодарность академику РАН JI. М. Зеленому, д. ф.-м.н. Х.В. Маловой и к. ф.-м.н. А. В. Артемьеву за полезные обсуждения и внимание к работе.
Работа выполнена при финансовой поддержке РФФИ (грант 11-02-01166-а), грантами НОЦ 14.740.11.0086 и НШ-3200.2010.2.
Список литературы
1. Плазменная гелиогеофизика / Под ред. J1.M. Зелёного, И. С. Веселовского. М„ 2008.
2. Вайнштейн С.И., Быков A.M., Топтыгин И.Н. Турбулентность, токовые слои и ударные волны в космической плазме. М., 1989.
3. Зелёный U.M., Малова Х.В., Артемьев A.B. и др. // Физика плазмы. 2011. 37, № 2. С. 137.
4. Harris E.G. // Nuovo Cimento. 1962. 23. P. 115.
5. Kan J.R. // Geophys. Res. 1973. 78. P. 3773.
6. Lembege В., Pellat R. // Phys. Fluids. 1982. 25. P. 1995.
7. Zelenyi L.M., Malova H.V., Popov V.Y. et al. // Nonlinear Proc. Geophys. 11. P. 579.
8. Bagenal F. // Institute of Physics Publishing 2001. Dirac House, Temple Back, Bristol, BS16BE, UK, 2001.
9. Laval G., Pellat R., Vuillemin M. // Plasma Phys. and Controlled Fusion Res. 1966. 2. P. 259.
10. Ибрагимов H.X. // Матем. и кибернетика. 1989. № 8.
The kinetic model of the two dimensional cylindrical current sheet I. Y. Vasko", V.Y. Popov
Department of Mathematics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
Space Research Institute, Russian Academy of Sciences, Moscow 117997, Russia. E-mail: a [email protected], b [email protected].
In this paper we develop a one class of solutions of the steady Vlasov-Maxwell equations, which describes two dimensional cylindrical current sheets with current directed azimuthally / = jg(p,z)eg. Magnetic field of these sheets has two components В = Bz(p,z)ez + Bp(p, z)ep. From mathematical point of view, we find solutions of the nonlinear equation in partial derivatives for some function u(p,z): д2и/дх2 d2u/dz2 = e^u, where x = p2. We apply methods of group theory to develop three-parameter class of solutions. We also derive asymptotic behavior of these solutions for large values of p and for p«0. We discuss applications of these solutions for description of current sheets in magnetospheres of planets with magnetic dipoles located near the ecliptic plane.
Keywords: system of Vlasov-Maxwell equations, nonlinear equations in partial derivatives, cylindrical current sheet, magnetospheres of planets. PACS: 52.25.Dg. Received 28 May 2011.
English version: Moscow University Physics Bulletin 1(2012).
Сведения об авторах
1. Васько Иван Юрьевич — аспирант; e-mail: [email protected].
2. Попов Виктор Юрьевич — докт. физ.-мат. наук, доцент; e-mail: [email protected].
