Two-particle kinetic model of the solar wind in approximation of frozen-in magnetic field .pdf Современные кинетические модели солнечного ветра основываются на уравнениях для одно- и двухчастичных функций распределения электронов и протонов плазмы по скоростям и формулируются в приближении стационарного потока полностью ионизованной квазинейтральной плазмы (например, [1, 2]). При этом замкнутые кинетические модели течения солнечной плазмы строятся на основе уравнений Власова, позволяющих найти как функции распределения электронов и ионов по скоростям, так и векторы напряженности электрического и магнитного полей. Решение и анализ задачи о распространении солнечного ветра в такой самосогласованной постановке затруднены, вследствие чего вводятся упрощающие предположения относительно магнитного поля. В одних случаях его влиянием пренебрегается [1, 2], а в других - магнитное поле определяется на основе ряда допущений [3 - 5], сформулированных Паркером [6]. Первый из этих двух подходов позволяет вычислить параметры потока плазмы однозначно - на основе функции распределения частиц по скоростям (как статистические моменты). Во втором случае постановка задачи получается переопределенной, так что скорость потока плазмы может быть найдена как путем вычисления первого статистического момента, что является общепринятым (например, [3, 4]), так и на основе уравнений Максвелла в рамках принятых допущений о поле скоростей и магнитном поле [5, 7]. Последний из названных подходов представлен в настоящей работе в применении к двухчастичной кинетической модели солнечного ветра.Рассмотрим солнечный ветер в приближении стационарного осесимметрично-го бесстолкновиетльного потока квазинейтральной плазмы, состоящей из электронов и протонов. Плазма испускается Солнцем, вращающимся с постоянной угловой скоростью Q, и распространяется в гравитационном, магнитном полях и в электрическом поле поляризации плазмы. Тогда уравнение для двухчастичнойфункции распределения электронов и протонов по скоростям f (r, ue, up) запишется в сферической системе координат, связанной с вращающимся Солнцем, в следующем виде:(1)где Fe = -meVp- eE - (e/c) ue x B + AFe, Fp = -mpV(p + eE + (e/c) up x B + AFpобозначают силы, действующие соответственно на протон и электрон в гравитационном поле с потенциалом tp(r), в поляризационном и магнитном полях с ос-редненными векторами напряженности Е и В (предполагается, что плазма распространяется в вакууме, вследствие чего векторы магнитной напряженности и магнитной индукции принимаются равными H = B ). Фиктивные (дополнительные) силы AF имеют следующие компоненты во вращающейся сферической системе координат r, 9, ф (j = e,p):2 +2AFjr =Я Uij + 2uЛ Q sin 9 + rQ2 sin2 9 ;(2)AF,q =-+ ctg9 + 2uQ cos 9 + rQ2 sin 9 cos 9 ;(3)rrд/г.ф = -J^JL -J^JLctge-2u]rQsin9-2ujeQcos9 .(4)Здесь me, ue, mp, up - массы и скорости электронов и протонов соответственно;±e - заряды протона и электрона. Индексы e, p относятся к электронам и протонам соответственно.Предполагается, что солнечный ветер истекает со сферической поверхности радиуса r0, расположенной в нижних слоях короны Солнца (r0~Rs, Rs - радиус Солнца), где плазму можно рассматривать как равновесную:Л = f (Г0 > Ue > Up ) = N0 (me m pi2nkBTQ ) eXP (-( me Ue + Mp Up ^Ik^ ) .(5)Здесь N0 , T0 - числовая плотность температура плазмы при r = r0. При этом предполагается, что испускаемые электроны и протоны имеют одинаковые температуры (Te0=Tp0=T0). Принятое допущение о квазинейтральности плазмы означает, что числовые плотности электронов и протонов при r = r0 также одинаковы (Ne0 = Np0 = 0,5 N0).Решение уравнения (1) - (4) методом характеристик позволяет найти двухчастичную функцию распределения как произвольную функцию шести первых интегралов:f ( r, ue > Up ) = f (Ee > Me > He > Ep > Mp > Ир ) ,(6)выражающих соответственно законы сохранения энергии E, момента количества движения M (при осевой симметрии потока) и магнитного момента р, (первый адиабатический инвариант) вдоль траекторий электронов и протонов:2 /2 22 /Ej = mjuj 2 + nijф-qjy-mj r sin 9-Q /2 = const ;(7) Mj = r sin 9 (+ r sin 9 Q + cm j ) = const ;(8)p.j = (cmjjqj)(v2/В) s constгде My - модуль вектора скорости частицы j-го сорта; у - потенциал электрического поля (E = V\|/), А - векторный потенциал магнитного поля (B = rot A ), Аф - его азимутальная компонента; vy - линейная скорость ларморовского вращения частицы j-го сорта. Как известно, первый адиабатический потенциал (9) выводится впредположении, что магнитное поле слабо меняется вдоль траектории частицы за период ее ларморовского вращения.Из граничных условий (5) определяется вид функции (6):f = No2 ((mm:/2nkBT0 / exp(-(meM2/2 + mpu2 /2 + (П -П0 ))/kBro) ,(10)где П = П(г) = mtp(r) - m(r sin 9-Q)2/2 обозначает суммарный потенциал электрона и протона, П0 = П(г0), m = me + mp. Границы области определения двухчастичной функции распределения (10) в фазовом пространстве скоростей электронов и протонов зависят от поляризационного потенциала у и азимутальной компоненты векторного потенциала магнитного поля Аф. Следовательно, статистические моменты (средние параметры потока - плотность, скорость и др.), вычисляемые путем интегрирования функции распределения по ее области определения, также зависят от у и Аф. Однако поляризационный потенциал у можно полностью исключить из решения задачи, если рассмотреть модель, описывающую статистику нейтральных динамических пар электрон - протон [2]. В этом случае как двухчастичная функция распределения f так и область ее определения зависят только от суммы первых интегралов (7):f = f (Ee + Ep , Me , Pe , Mp , Ир ) .(П)Этот подход используется ниже для получения верхней оценки стационарного потока частиц плазмы Nu(r) (17).Магнитное поле задается на основе ключевых предположений Паркера о межпланетном магнитном поле, как вмороженном в стационарный осесиммет-ричный поток плазмы, истекающий с нулевой начальной скоростью с поверхности вращающегося Солнца и имеющий нулевую меридиональную компоненту скорости. Отметим, что в отличие от Паркера угол истечения плазмы не фиксируется и скорость потока не считается постоянной. Указанные допущения позволяют вывести из уравнений Максвелла (div B = 0 и rot B = (4п / c)J) не только зависимости для компонент вектора магнитной индукции B [6]:Br = Bre(9)(r0/r)2, B6 = 0, Вф = Br (r - r0)Qsin9/Ur, (12) но и следующие выражения для компонент скорости солнечного ветра U:Ur (r, 9) = \Br0 (9)| sin2 9-g (r), U6 = 0, иф = (r - r0 )Q sin 9 .(13)При выводе зависимости Ur(r, 9) (13) учитывалось, что при стационарном радиальном расширении квазинейтральной плазмы радиальный ток отсутствует (Jr = 0). Если меридиональный ток также равен нулю (Jr = 0, J& = 0), что соответствует гипотезе о вмороженном магнитном поле, когда поперечным дрейфом заряженных частиц можно пренебречь, то из уравнений Максвелла находится зависимость радиальной скорости от гелиоцентрического расстояния:Ur (r, 9) = а | Br о (9)| sin2 9 (1 - r0f r). (14)Здесь J - вектор тока, gi(r), a - соответственно произвольные функция и константа интегрирования уравнений Максвелла.Заметим, что если рассмотреть полностью бестоковое приближение (Jr = 0, Je = 0, J,], = 0, т.е. азимутальный ток также принять равным нулю), то в выражении для радиальной скорости солнечного ветра (14) коэффициент Д-о(9) = Bro = const и радиальная компонента магнитной индукции (12) остается функцией только гелиоцентрического расстояния r. Согласно данным измерений миссии Ulysses (сплошная линия на рис. 1 [5]), профиль скорости солнечного ветра U(oo, 9) имеет форму бабочки («butterfly» diagram) - уменьшение скорости потока солнечной плазмы наблюдается вблизи экваториальной плоскости Солнца (9~л/2) и при приближении к его полюсам (9~0, 9~л). Множитель sin29 обеспечивает качественное согласование зависимости (14) с указанными наблюдениями для всех гелиоширот, кроме приэкваториальной области. Наблюдаемое уменьшение скорости U в этой зоне может быть достигнуто, если радиальная компонента магнитной индукции на сфере истечения плазмы r = r0 Br0 = Br0(9) имеет минимум вблизи экватора. Такой характер указанной зависимости подтверждается данными наблюдений для магнитного поля фотосферы в период активного Солнца [8], которые могут быть аппроксимированы следующей формулой [5]:а| Br0(9)| = а0(1 + 16|sin79|). (15) Если предположить, что зависимость (15) можно спроецировать (экстраполировать) в нижние слои солнечной короны, где задаются граничные условия для истекающей плазмы (r = r0), то полученный профиль предельной скорости солнечного ветра (14) с учетом (15)Ur(oo, 9) = а0(1 + 16 |sin7 9|) sin2 9 (16)качественно согласуется с данными наблюдений миссии Ulysses (рис. 1) [5]. На рис. 1 теоретическая зависимость (16) и данные Ulysses представлены соответственно пунктирной и сплошной линиями. Константа a0 = 26,5 рассчитана на основе эмпирического значения предельной скорости солнечного ветра (450 км/с) в плоскости эклиптики [9]. Теоретические результаты (16) не рассматриваются вблизи экваториальной плоскости (9 » п/2) и полюсов Солнца (9 » 0, 9 » п), поскольку приближение вмороженного магнитного поля в этих областях неприменимо. Анализ уравнения для поперечного дрейфа частиц плазмы в заданном магнитном поле (12) показывает, что в указанных областях скорости движения быстрых частиц поперек силовых магнитных линий превышают скорости их продольного движения. Дрейф частиц также значителен при r ~ r0, что приводит к рассогласованию теоретических результатов с данными наблюдений у поверхности истечения плазмы.Уравнение для дрейфа частиц приводит также к соотношению, связывающему поперечные температуры испускаемых электронов Teio и протонов Tpl0 с гравитационным потенциалом ср0, линейной скоростью вращения Солнца ua0 = r0Q и углом выхода магнитных силовых линийВфо/во при r = r0:Te±0 + Tp±0 = (mK|/6kB)(1 -«дo/h-0,5Бфо/Бо2) .(17)Для Солнца Uqq/|ср0| 1,5r0 в пределах 35%. Кроме того, радиальный профиль скорости Ur(r, 9i) лежит в пределах разброса данных измерений для высокоширотного солнечного ветра (9 * 30°) [11].Зависимость числовой плотности от меридионального угла (20) на уровне орбиты Земли (r = 1 а.е.) качественно согласуется с данными наблюдений Ulysses, что иллюстрирует рис. 2, на котором теоретическая кривая и наблюдательные данные представлены пунктирной и сплошной линиями соответственно [5]. Сравнение проводилось при следующих значениях параметров: r0=3Rs, температура T0 рассчитывалась по формуле (18), плотность N0 вычислялась из согласования теоретического (20) и наблюдательного значений плотности солнечного ветра на уровне орбиты Земли.Таким образом, представленная кинетическая модель стационарного осе-симметричного солнечного ветра в заданном магнитном поле паркеровского типа позволяет получить с использованием уравнений Максвелла пространственную зависимость для скорости, потока частиц и числовой плотности плазмы, которые согласуются с данными наблюдений солнечного ветра как в плоскости эклиптики, так и на высоких гелиоширотах. Анализ расхождений теоретических результатов и наблюдательных данных с использованием уравнений Максвелла и уравнения дрейфа частиц в заданном магнитном поле позволяет выявить роль допущений об отсутствии токов в плазме и предложить соответствующую корректировку модели.

Скачать электронную версию публикации