Кинетическое моделирование распределения межзвездных атомов водорода в гелиосфере и анализ экспериментальных данных SOHO/SWAN по рассеянному солнечному Лайман-альфа излучению тема автореферата и диссертации по астрономии, 01.03.03 ВАК РФ

Катушкина, Ольга Александровна АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2013 ГОД ЗАЩИТЫ
   
01.03.03 КОД ВАК РФ
Диссертация по астрономии на тему «Кинетическое моделирование распределения межзвездных атомов водорода в гелиосфере и анализ экспериментальных данных SOHO/SWAN по рассеянному солнечному Лайман-альфа излучению»
 
Автореферат диссертации на тему "Кинетическое моделирование распределения межзвездных атомов водорода в гелиосфере и анализ экспериментальных данных SOHO/SWAN по рассеянному солнечному Лайман-альфа излучению"

На правах рукописи

КАТУШКИНА ОЛЬГА АЛЕКСАНДРОВНА

КИНЕТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРЕДЕЛЕНИЯ МЕЖЗВЕЗДНЫХ АТОМОВ ВОДОРОДА В ГЕЛИОСФЕРЕ И АНАЛИЗ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ SOHO/SWAN ПО РАССЕЯННОМУ СОЛНЕЧНОМУ ЛАЙМАН-АЛЬФА ИЗЛУЧЕНИЮ

Специальность - 01.03.03 - Физика Солнца 01.02.05 - Механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата физико-математических наук

005540833

МОСКВА

- 2013

005540833

Работа выполнена на кафедре аэромеханики и газовой динамики механико-математического факультета федерального государственного бюджетного образовательного учреждения высшего профессионального образования "Московский государственный университет имени М.В.Ломоносова".

Научный руководитель: доктор физико-математических наук,

Измоденов Владислав Валерьевич

Официальные оппоненты: Курт Владимир Гдалевич, д.ф.-м.н., заведующий

отделом «Квантовой астрофизики» Астрокосмического центра Физического института им. П.Н. Лебедева РАН

Угольников Олег Станиславович,

к.ф.-м.н., старший научный сотрудник

ФГБУН Института космических исследований РАН

Ведущая организация: Федеральное государственное бюджетное

учреждение науки Институт земного магнетизма, ионосферы и распространения радиоволн им. Н.В. Пушкова РАН

Защита состоится 20 ноября 2013 года в 11 часов на заседании диссертационного совета Д 003.113.03 при ИКИ РАН по адресу: 117997, г. Москва, ул. Профсоюзная 84/32, Институт Космических исследований РАН, конференц-зал.

С диссертацией можно ознакомиться в библиотеке Института космических исследований РАН.

Автореферат разослан « /I » ¿¿/^^ 2013 г.

Ученый секретарь диссертационного совета кандидат физико-математических наук

Буринская Т.М

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы

В диссертации проводится теоретический анализ эволюции функции

распределения межзвездных атомов водорода по скоростям внутри гелиосферы -области космического пространства, занятой солнечным ветром. Разработана трехмерная нестационарная кинетическая модель распределения атомов водорода, в которой учитываются эффекты, связанные 1) с 11-летним циклом солнечной активности; 2) с зависимостью параметров солнечного ветра от гелиошироты; 3) с изменением функции распределения межзвездных атомов при их прохождении через область взаимодействия солнечного ветра с межзвездной средой. Построенная численная модель применяется для анализа экспериментальных данных по рассеянному солнечному Лайман-альфа излучению, полученных на космическом аппарате (КА) SOHO (прибором SWAN). Актуальность данной работы обусловлена, во-первых, интересом к исследованию границы гелиосферы, возросшим в последнее время благодаря появлению новых уникальных данных с КА Вояджер-1 и 2, а также IBEX - Interstellar Boundary Explorer, и, во-вторых, большим количеством экспериментальной информации о рассеянном солнечном Лайман-альфа излучении, накопленной за последние 25 лет.

Исторически, именно данные о рассеянном Лайман-альфа излучении привели к идее о существовании межзвездного ветра, а также к развитию теоретических моделей проникновения межзвездных нейтральных атомов внутрь гелиосферы. В самом начале космической эры (Shklovsky, 1959) при первых ракетных запусках было обнаружено излучение в ультрафиолетовом (УФ) диапазоне, которое невозможно было объяснить излучением от ближайших к Солнцу звезд или рассеянием солнечного излучения на нейтральных атомах в атмосфере Земли. В работах Курта & Гермогеновой (1967) и Курта & Сюняева (1967) была высказана гипотеза о галактическом происхождении измеряемого излучения. В дальнейшем более детальные измерения на аппарате OGO-5 на основании обнаруженного параллакс-эффекта показали, что единственным возможным источником

рассеянного Лайман-альфа излучения могут быть межзвездные атомы водорода, которые проникли внутрь Солнечной системы за счет относительного движения Солнца и окружающей его межзвездной среды. Такое явление получило название межзвездного ветра.

Таким образом, межзвездные атомы водорода проникают из Локальной Межзвездной Среды (ЛМС) внутрь гелиосферы, где они могут быть измерены с помощью прямых или косвенных методов. Внутри гелиосферы распределение атомов водорода является источником информации как о параметрах ЛМС, так и о свойствах переходной области между солнечным и межзвездным ветрами (в литературе эту область называют гелиосферным ударным слоем). Кроме того, солнечное Лайман-альфа излучение рассеивается на межзвездных атомах водорода. Спектральные свойства рассеянного излучения определяются функцией распределения атомов водорода по скоростям. Следовательно, измерения рассеянного солнечного Лайман-альфа излучения открывают возможности для удаленной диагностики распределения межзвездных атомов водорода в гелиосфере.

Начиная с 1980-х годов прошлого века основной теоретической моделью распределения атомов водорода в гелиосфере являлась, так называемая, классическая горячая модель (Fahr, 1971), которая использовалась в течение многих лет для анализа и интерпретации данных по рассеянному солнечному Лайман-альфа излучению (Lallement et al., 1984, 1985 и др.). Тем не менее, эта модель обладает рядом недостатков. А именно, в ней не учитываются нестационарные и трехмерные эффекты, связанные с 11-летнем циклом солнечной активности и с зависимостью параметров солнечного ветра от гелиошироты, а также не учитывается существование области гелиосферного ударного слоя.

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

Разработанная модель применяется для исследования спектральных характеристик рассеянного солнечного Лайман-альфа излучения в гелиосфере (глава 2), а также для анализа и интерпретации экспериментальных данных, полученных на космическом аппарате SOHO (прибором SWAN) (глава 3).

Цели и задачи диссертационной работы состоят в следующем:

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

2. Исследовать, как немаксвелловские свойства функции распределения атомов водорода во внешней гелиосфере (на 90 а.е. от Солнца) влияют на распределение (по пространству и скоростям) атомов водорода вблизи Солнца. Немаксвелловский характер функции распределения во внешней гелиосфере связан с процессом перезарядки в области гелиосферного ударного слоя.

3. Проанализировать спектральные характеристики рассеянного солнечного Лайман-альфа излучения в гелиосфере. Определить, как немаксвелловские свойства в распределении атомов водорода во внешней гелиосфере влияют на спектральные характеристики рассеянного Лайман-альфа излучения.

4. Провести теоретическое моделирование интенсивности рассеянного солнечного Лайман-альфа излучения и сравнить результаты теоретических расчетов с экспериментальными данными, полученными на КА SOHO (прибором SWAN).

5. Провести исследование зависимости потока массы солнечного ветра от гелиошироты на основании полных карт неба в интенсивности рассеянного излучения, полученных из данных SOHO/SWAN.

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

эффекты, а также кинетические эффекты, связанные с перезарядкой на границе

гелиосферы. На основе этой модели были проанализированы данные ЗОНО/ЗАУАЫ.

Автором впервые получены и выносятся на защиту следующие основные

положения:

1. Впервые классифицированы немаксвелловские свойства функции распределения атомов водорода на расстоянии 90 а.е. от Солнца, возникающие за счет процесса перезарядки межзвездных атомов на протонах ЛМС и солнечного ветра в области гелиосферного ударного слоя.

2. Разработана трехмерная нестационарная кинетическая модель распределения межзвездных атомов водорода внутри гелиосферы с учетом эффектов, связанных с изменением функции распределения атомов в области гелиосферного ударного слоя. Данная модель является эффективным инструментом для анализа многочисленных экспериментальных данных.

3. Исследовано влияние гелиосферного ударного слоя на моменты функции распределения атомов водорода вблизи Солнца, а также на характеристики рассеянного солнечного Лайман-альфа излучения. Показано, что наиболее чувствительными к немаксвелловским свойствам функции распределения атомов во внешней гелиосфере являются кинетическая температура атомов и ширина линии рассеянного Лайман-альфа излучения.

4. Проведен детальный теоретический анализ зависимости ширины линии рассеянного Лайман-альфа излучения от направления луча зрения наблюдателя. Показано, что немонотонное поведение ширины линии, которое ранее было обнаружено в данных 80Н0/8\УАМ, может быть частично объяснено в рамках текущих моделей только с учетом многокомпонентного характера плазмы межзвездной среды.

5. На основании разработанной в настоящей диссертации трехмерной и нестационарной модели распределения атомов водорода в гелиосфере был проведен анализ полных карт неба в интенсивности рассеянного Лайман-альфа

излучения, полученных на КА ЗОНО/БХУАЫ в 1996-2009 гг. При этом впервые в модели использовалась реальная зависимость параметров солнечного ветра от гелиошироты, основанная на различных экспериментальных данных.

6. Установлено, что для объяснения данных ЗОНО/БХУАН по интенсивности рассеянного Лайман-альфа излучения в период максимума солнечной активности (с 2001 по 2005 гг.) необходимо предположить, что график зависимости потока массы солнечного ветра от гелиошироты имеет два максимума на средних широтах (± 50-60°). Информация о гелиоширотных вариациях потока массы солнечного ветра имеет принципиальное значение для глобального моделирования гелиосферного ударного слоя и исследования процессов, происходящих на границе гелиосферы.

Достоверность результатов, полученных в диссертации, основана на следующем:

1. Для моделирования распределения атомов водорода использовался наиболее адекватный в данном случае, кинетический подход. Необходимость использования кинетического подхода связана с большой длиной свободного пробега атомов водорода, сравнимой с характерным размером гелиосферы. Также для задания в модели граничного условия на расстоянии 90 а.е. от Солнца использовались результаты глобальных самосогласованных кинетико-газодинамических моделей взаимодействия солнечного ветра с межзвездной средой (Вагапоу & Ма1ата, 1993; 1ипос1епоу й а1., 2009 и др.), которые отлично зарекомендовали себя ранее тем, что с их помощью удается объяснить множество экспериментальных фактов.

2. Разработанный численный алгоритм тестировался на классических аналитических решениях в стационарном осесимметричном случае. Также проводилось дополнительное тестирование программы в нестационарном трехмерном случае с целью выявления наиболее оптимальных параметров расчетных сеток.

3. Теоретические результаты, полученные на основании разработанной в

диссертации модели, хорошо согласуются с экспериментальными данными

SOHO/SWAN по рассеянному Лайман-альфа излучению.

Практическая ценность работы заключается в том, что разработанную в диссертации модель можно эффективно применять для анализа экспериментальных данных по измерениям потоков межзвездных нейтральных атомов (водорода, гелия и др.), а также параметров рассеянного на атомах солнечного излучения с уже существующих и будущих космических аппаратов. В частности, эта модель представляет собой удобный инструмент для анализа данных по потокам межзвездных атомов водорода и гелия, измеряемых с 2009 г. на КА IBEX прибором IBEX-Lo. Кроме того, модель можно использовать для определения научных задач будущих российских и зарубежных космических миссий. В частности, в настоящее время результаты, полученные в диссертации, используются для составления научной программы прибора "Лаймус", разрабатываемого в ИКИ РАН. Задачей прибора является измерение интенсивности и спектра рассеянного солнечного Лайман-альфа излучения. В случае успешного запуска прибора, разработанная в диссертации модель будет также использоваться и для анализа полученных на "Лаймусе" данных.

Апробация

Работы, вошедшие в диссертацию, докладывались и обсуждались на научно-исследовательских семинарах по аэромеханике и газовой динамике механико-математического факультета МГУ им. М.В. Ломоносова (рук. - акад. Черный Г.Г.), семинарах Института космических исследований РАН (рук. — академик Зеленый Л.М.), семинарах лаборатории физической газовой динамики Института проблем механики им. А.Ю. Ишлинского РАН (рук.- проф. Баранов В.Б.).

Основные положения и результаты, вошедшие в диссертацию, докладывались на российских и международных конференциях, в том числе:

• на У-ІХ конференциях молодых ученых «Фундаментальные и прикладные космические исследования» (ИКИ РАН, Москва, 2009-2012 гг.);

• на Международном научном форуме студентов, аспирантов и молодых ученых «Ломоносов» (МГУ, Москва, 2009-2011 гг.);

• на конференции Ломоносовские чтения (МГУ, Москва, 2010 г.);

• на конференции «Физика плазмы в солнечной системе» (ИКИ РАН, Москва, 2009-2011 гг.);

• на Всероссийской школе-семинаре «Аэрофизика и физическая механика классических и квантовых систем» (ИПМех РАН, Москва, 2010-2011 гг.);

• на X Всероссийском съезде по фундаментальным проблемам теоретической и прикладной механики (Нижний Новгород, 2011 г.);

• на 18-й конференции молодых ученых по астрономии и космической физике (Национальный университет им. Тараса Шевченко, Киев, Украина, 2011 г.);

• на Всероссийской ежегодной конференции по физике Солнца "Солнечная и солнечно-земная физика - 2008" (Главная астрономическая обсерватория РАН, Санкт-Петербург, 2008 г.)

• на ежегодных конференциях Европейского геофизического общества (Ейи) (г. Вена, Австрия, 2010, 2012 гг.);

• на 12-й международной конференции "Солнечный ветер" (г. Сан-Мало, Франция, 2009);

• на 11-й Научной Ассамблее Ассоциации Геомагнетизма и Аэрономии (ІАОА) (г. Сопрон, Венгрия, 2009 г.);

• на 37-й, 38-й и 39-й научных ассамблеях С08РАЯ (г. Монреаль, Канада, 2008; г. Бремен, Германия, 2010; г. Май сор, Индия, 2012);

• на XXVIII Генеральной Ассамблее Международного астрономического общества (1Аи) (г. Пекин, Китай, 2012 г.).

Публикации и личный вклад автора

Представленные в диссертации результаты опубликованы в четырех статьях в рецензируемых российских и международных журналах из перечня ВАК, в главе из коллективной монографии, напечатанной в издательстве Springer, а также в трех статьях из сборников трудов конференций.

Вклад автора во все работы по теме диссертации является основным. Все результаты, выносимые на защиту, были получены лично автором диссертации. Автором осуществлялись: разработка и тестирование всех численных программ; проведение всех расчетов; анализ полученных результатов; сравнение результатов с экспериментальными данными; подготовка и написание текстов публикаций; а также переписка с редакциями журналов и рецензентами. Структура и объем работы

Диссертация состоит из введения, трех глав и заключения. Общий объем диссертации составляет 147 страниц. В диссертации содержится 32 иллюстрации. Список используемой литературы состоит из 127 библиографических ссылок.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

Глава 1 диссертации посвящена описанию кинетической модели распределения межзвездных атомов водорода в гелиосфере. Рассматриваются как упрощенная постановка задачи — осесимметричная стационарная, так и полная - трехмерная нестационарная.

В разделе 1.1 представлен исторический обзор работ, связанных с моделированием распределения межзвездных атомов водорода в гелиосфере. Впервые задача о движении облака частиц в гравитационном поле звезды с учетом конечной температуры была решена аналитически в работе Danby & Camm (1957). В

дальнейшем (в работах Meier, 1977 и Wu & Judge, 1979) это решение было обобщено и адаптировано для описания распределения межзвездных атомов водорода в гелиосфере. Такая модель стала называться классической горячей моделью. В подразделе 1.1.2 представлено аналитическое решение для функции распределения атомов водорода по скоростям в случае классической горячей модели.

В подразделе 1.1.3 описываются современные представления о структуре области взаимодействия солнечного ветра с межзвездной средой. В подразделе 1.1.4 приводится обзор работ, посвященных современному развитию классической горячей модели.

В разделе 1.2 представлена разработанная автором диссертации кинетическая модель распределения межзвездных атомов водорода в гелиосфере в осесимметричном стационарном случае. Эта модель совмещает в себе сравнительную простоту классической горячей модели с возможностью учитывать немаксвелловские свойства в функции распределения атомов водорода на границе гелиосферы, связанные с существованием области гелиосферного ударного слоя. В подразделе 1.2.1 приводятся основные уравнения для описания распределения атомов водорода. Кинетическое уравнение для функции распределения атомов по скоростям выглядит следующим образом:

н, .fM+F(r) ШьцИ* -т • /(г, „).

Здесь f(r,w) - функция распределения атомов водорода по скоростям w в зависимости от пространственной координаты г. Результирующая сила F(r), действующая на атом вблизи Солнца, является суммой силы гравитационного притяжения к Солнцу (Fg) и силы радиационного отталкивания (Frad). Результирующая сила определяется безразмерным параметром И =| Frad | /1 Fg \, который в общем случае зависит от радиальной компоненты скорости атома, времени и гелиошироты, но в рассматриваемом в разделе 1.2 стационарном осесимметричном случае полагалось, что jJ=const. Параметр Р(г) представляет собой

суммарную частоту ионизации, складывающуюся из частоты фотоионизации и частоты перезарядки. Общие выражения для вычисления частоты перезарядки и фотоионизации приведены в подразделе 1.2.2 диссертации.

В подразделе 1.2.3 обсуждаются немаксвелловские эффекты в функции распределения межзвездных атомов водорода на расстоянии 90 а.е. от Солнца, вызванные процессом перезарядки первичных межзвездных атомов водорода на протонах ЛМС и солнечного ветра в области гелиосферного ударного слоя. В результате перезарядки в этой области образуются, так называемые, вторичные межзвездные атомы, функция распределения которых зависит от локальных свойств плазмы в области их рождения.

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

В подразделе 1.2.5 описан метод характеристик для решения шестимерного кинетического уравнения.

В подразделе 1.2.6 представлены результаты, полученные в рамках разработанной модели для различных способов задания граничного условия. Также проводится сравнение полученных результатов с результатами самосогласованной кинетико-газодинамической модели взаимодействия солнечного ветра с межзвездной средой (Вагапоу & Ма1ата, 1993). В целом результаты проведенного сравнения свидетельствуют о том, что изменения функции распределения атомов водорода в гелиосферном ударном слое (которые вообще не учитываются в классической горячей модели) приводят к значительным изменениям в распределении (как в пространстве, так и по скоростям) атомов водорода внутри гелиосферы. Показано, что наиболее чувствительной характеристикой к эффектам границы гелиосферы является кинетическая температура атомов

ад: |/Хг,и>У(Г,(г)-у<г,)гс1».

В разделе 1.3 диссертации описана трехмерная нестационарная кинетическая модель. В рамках этой модели решается нестационарное кинетическое уравнение, в котором учитывается зависимость частоты ионизации атомов Ре и параметра Ц от времени I и гелиошироты X. Кроме того, в модели учитывается зависимость параметра ц от радиальной компоненты скорости атома \\'г. В подразделе 1.3.1 представлена математическая постановка задачи (уравнения и граничные условия) в трехмерном нестационарном случае. В подразделе 1.3.2 описаны способы задания входных параметров задачи: балансового параметра ц(1,А.,\уг), частот фотоионизации РрКеОЛ) и перезарядки ре)(Е(1:Д) на 1 а.е. Для задания этих зависимостей используются различные экспериментальные данные о параметрах солнечного ветра и излучения. В частности, используются базы данных ОММ-2, 80ЬАЯ2000, ЬАЭР. Гелиоширотная зависимость частоты ионизации атомов задается на основании результатов анализа данных ЗОНО/БХУАЫ по рассеянному Лайман-альфа излучению ((ЗиетегаБ й а1.2006а).

В подразделе 1.3.3 приводятся результаты численного моделирования распределения атомов водорода в гелиосфере в трехмерном нестационарном случае. Показано, что на расстоянии 1 а.е. от Солнца нестационарные эффекты приводят к значительным (до 30%) колебаниям концентрации, средней скорости и средней температуры атомов в течение 11-летнего цикла солнечной активности, а на расстоянии 10 а.е. от Солнца параметры атомов остаются практически постоянными в течение цикла. На рис. 1 показаны концентрация, средняя скорость Уг и

кинетическая температура Тг(г): ^(г,ю)-(Уг(г)~м>г)1(1}*> первичных и

вторичных межзвездных атомов водорода на расстоянии 1 а.е. от Солнца в зависимости от гелиошироты и гелиодолготы. Видно, что зависимость частоты перезарядки от гелиошироты приводит к раздвоению максимумов и минимумов в распределении концентрации и температуры атомов водорода.

В разделе 1.4 диссертации приводится заключение к первой главе.

вторичные атомы

первичные атомы

■ 7350 й||7250 Г™17000 6500 6250 ' '16000

......15500

5000

■ 4500 4250 4000

Рис.1. Концентрация (А, О), радиальная скорость (В, Е) и радиальная кинетическая температура Тг (С, Р) вторичных (графики слева) и первичных (графики справа) межзвездных атомов на расстоянии 1 а.е. от Солнца в зависимости от направления луча зрения наблюдателя. По осям отложены гелиоширота и гелиодолгота в гелиографических координатах. Эти результаты получены в условиях солнечного минимума (1996-го года), когда гелиоширотная зависимость частоты перезарядки проявляется особенно сильно.

Глава 2 диссертации посвящена теоретическому исследованию спектральных характеристик рассеянного солнечного Лайман-альфа излучения. Описана математическая модель переноса излучения, которая позволяет вычислять спектры рассеянного Лайман-альфа излучения в произвольном направлении луча зрения, если известно распределение атомов водорода вдоль луча зрения. Исследуется вопрос о том, как немаксвелловские эффекты (связанные с процессом перезарядки в гелиосферном ударном слое) в функции распределения атомов водорода проявляются в спектрах рассеянного солнечного Лайман-альфа излучения на 1 а.е.

В разделе 2.1 приводится исторический обзор работ, связанных с измерением и исследованием рассеянного солнечного Лайман-альфа излучения.

В разделе 2.2 описана модель переноса излучения. В подразделе 2.2.1 приводится способ получения формального аналитического решения уравнения переноса излучения. В данной работе рассматривается только излучение за счет однократно рассеянных фотонов, при этом учитывается поглощение рассеянного излучения между точкой рассеяния и наблюдателем. В этом случае решение уравнения переноса записывается аналитически в виде интеграла вдоль луча зрения.

В разделе 2.3 представлены результаты расчетов спектральных характеристик рассеянного солнечного Лайман-альфа излучения на основании распределений межзвездных атомов водорода. Указанные распределения были получены в рамках различных кинетических моделей из первой главы диссертации.

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

В подразделе 2.3.2 проводится детальное теоретическое исследование зависимости ширины линии (Tios) рассеянного Лайман-альфа излучения от направления луча зрения наблюдателя. Изучается влияние на эту зависимость многокомпонентного характера плазмы ЛМС, а также трехмерных и нестационарных эффектов. В работах Costa et al. (1999) и Quemerais et al. (2006b) на основании анализа экспериментальных данных SOHO/SWAN было показано, что в

период 1996-1997 гг. наблюдался минимум ширины линии рассеянного излучения при 0 = 5О'-6О° ^ тогда как классическая однокомпонентная горячая модель предсказывает монотонный рост ширины линии при увеличении угла 9. Результаты исследования, проведенного в диссертации, показывают, что ни одна из имеющихся на данный момент моделей распределения атомов водорода в гелиосфере не приводит к образованию минимума в T|OS, сравнимого по размеру с наблюдаемым в данных SOHO/SWAN. Тем не менее, учет многокомпонентного характера плазмы JIMC в модели позволяет добиться качественного согласия с экспериментальными данными (см. рис. 2). Многокомпонентный характер плазмы учитывается в модели с помощью задания граничного условия, основанного на результатах многокомпонентной модели гелиосферного ударного слоя (Malama et al., 2006).

Рис.2. Графики зависимости ширины линии Tios рассеянного Лайман-апьфа излучения от угла 9. Угол 0 отсчитывается от направления навстречу вектору скорости ЛМС. Данные результаты получены в рамках осесимметричных стационарных моделей распределения атомов водорода в гелиосфере с граничными условиями в виде трехмерного нормального распределения. Пунктирная кривая соответствует модели без учета многокомпонентного характера плазмы ЛМС; сплошная кривая - модели с учетом многокомпонентного характера плазмы.

В подразделе 2.3.3 проводится исследование зависимости моментов спектра

рассеянного излучения от значений определяющих параметров задачи ц и рЕ в осесимметричном стационарном случае. В разделе 2.4 приводится заключение к главе 2.

Глава 3 посвящена исследованию зависимости потока массы солнечного ветра от гелиошироты на основании анализа карт неба в интенсивности рассеянного солнечного Лайман-альфа излучения.

В разделе 3.1 приведено введение к данной проблеме.

В Ре,Е.2(и) (Э1) 8 Ц.Е-07

1998 2002 2006 2010

время

Рис.3. Частота перезарядки на 1 а.е. в зависимости от времени и гелиошироты. График А соответствует результатам первого способа, основанного на анализе данных по рассеянному

Лайман-альфа излучению ( Рех E l (t, Я) ). График В соответствует результатам второго способа, основанного на данных по межпланетным мерцаниям, а также измерениях Ulysses ( Рех,Е,г

В разделе 3.2 описаны два имеющихся на сегодняшний день способа задания зависимости частоты перезарядки (или потока массы) солнечного ветра от времени и гелиошироты. Первый способ основан на методе, изложенном в (ЗиетегаІБ еі аі. (2006а), и связан с анализом данных 80Н0/8\УАЫ. Частоту перезарядки, полученную данным способом, мы будем обозначать Д). Второй способ

основан на результатах работы Бокої еі аі. (2012). В этой работе были найдены гелиоширотные и временные вариации скорости солнечного ветра (исходя из

1998 2002 2006 2010

время

измерений межпланетных мерцаний), а также плотности солнечного ветра (исходя из измерений на КА Ulysses и некоторых предположений). Частоту перезарядки, полученную этим способом, мы будем обозначать Дх.едО Д). Отметим, что частота перезарядки напрямую определяется потоком массы солнечного ветра. На рис. 3 показаны графики зависимости полученных частот перезарядки от времени и гелиошироты. Видно, что эти функции качественно различны. В частности, функция , имеет характерные максимумы на средних широтах в 2000-2005 гг. Эта особенность отсутствует во второй функции Д-^едОД). Как будет видно в дальнейшем, это различие в частотах перезарядки приводит к качественным различиям в картах интенсивности рассеянного Лайман-альфа излучения.

В разделе 3.3 приводятся результаты расчетов карт неба в интенсивности рассеянного солнечного Лайман-альфа излучения. Чтобы выяснить, насколько

различия в частоте перезарядки (т.е. между функциями А.,,, 1 и Рех,Е 2 ) сказываются на распределении рассеянного солнечного Лайман-альфа излучения, в диссертации были проведены расчеты распределения атомов водорода в гелиосфере в рамках разработанной трехмерной нестационарной кинетической модели, описанной в главе 1. При этом зависимость частоты перезарядки от времени и гелиошироты задавалась либо в виде Рех,к,\, либо в виде Ра,Е,2 ■ После этого были вычислены карты интенсивности рассеянного солнечного Лайман-альфа излучения в рамках описанной в главе 2 модели переноса излучения. Расчеты были проведены отдельно для периодов минимума (1997 и 2009 гг.) и максимума (2002-2005 гг.) солнечной активности. Показано, что в период минимума солнечной активности обе модели с различными способами задания частоты перезарядки приводят к качественно похожим результатам и к хорошему согласию с данными SWAN. Однако, в период максимума солнечной активности результаты качественно различаются.

01/06/2003

модель с (3,

модель с (3,

01/06/2002

модель С Р,

модель С р,

экл. долгота

Рис.4. Результаты расчетов в период максимума солнечной активности (2002-2003 гг.). Графики A-F представляют полные карты неба в интенсивности рассеянного солнечного Лайман-альфа излучения (по осям отложены эклиптическая широта и долгота). Все интенсивности нормированы таким образом, что для каждой карты нормированная интенсивность в одном выбранном направлении (характеризуемым 80.5° по широте и 120.5° по долготе) равна 1. Графики А и D показывают экспериментальные данные SOHO/SWAN. Отметим, что полоса с нулевой интенсивностью, проходящая через всю карту в данных SWAN, не имеет физического смысла и связана с тенью аппарата. Графики В и Е представляют теоретические результаты, полученные в рамках нашей модели с частотой перезарядки Рех,Е,\ ■ Графики С и F представляют теоретические результаты, полученные в рамках той же модели, но с частотой перезарядки PeXiE,i •

На рис. 4 показаны карты неба в интенсивности рассеянного Лайман-альфа излучения, полученные для 2002-2003 гг. Верхняя строка соответствует данным SWAN, средняя строка - результатам модели с частотой перезарядки РСх,е,\ > а нижняя строка - результатам модели с частотой перезарядки Рех,е,2 • Видно, что в данных SWAN, как и в результатах модели с максимум интенсивности

располагается в носовой части гелиосферы при небольших гелиоширотах. В то же время, в результатах модели с Рех,Е,г максимум интенсивности располагается ближе к северному эклиптическому полюсу. Такие различия связаны с зависимостью частоты ионизации атомов от гелиошироты в период максимума солнечной активности. А именно, в данных Рех.Е.i в это время присутствует два максимума на

средних широтах, а в данных Pa,E,i - один максимум в районе солнечного экватора (т.е. при нулевой гелиошироте). Были проведены специальные тестовые расчеты, которые подтвердили, что изменение других параметров модели (граничных условий и параметра ц) не позволяет добиться качественного совпадения в картах интенсивности между результатами модели с и данными SWAN. При этом

результаты модели с Рех,Е,\ приводят к отличному совпадению с экспериментальными данными. Это говорит о том, что характерные максимумы частоты перезарядки fia-E,i на средних широтах в период максимума солнечной активности действительно должны присутствовать, чтобы было возможно объяснить данные SOHO/SWAN по интенсивности рассеянного Лайман-альфа излучения в рамках существующих моделей.

В разделе 3.4 приводится заключение к главе 3.

Заключение содержит основные результаты и выводы диссертационной работы.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Разработана новая трехмерная нестационарная кинетическая модель распределения межзвездных атомов водорода по пространству и скоростям в

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

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

3. Проанализировано влияние нестационарных и трехмерных эффектов на распределение атомов водорода в гелиосфере. Впервые в модели были учтены реальные данные о гелиоширотной зависимости параметров солнечного ветра, полученные на основании экспериментальных данных. Показано, что зависимость частоты ионизации атомов от гелиошироты приводит к качественным изменениям в концентрации и температуре атомов на 1 а.е.

4. Проведено исследование влияния кинетических эффектов гелиосферного ударного слоя на измеряемые с орбиты Земли характеристики рассеянного солнечного Лайман-альфа излучения. Показано, что использование ранее общепринятой "классической горячей модели" приводит к ошибкам в интенсивности Лайман-альфа излучения (от 2% в носовой части гелиосферы до 25% в хвостовой части), а также к ошибкам в ширине линии рассеянного излучения (до 18%).

5. Проведено детальное исследование зависимости ширины линии рассеянного Лайман-альфа излучения от направления луча зрения наблюдателя. Показано, что в рамках имеющихся на сегодняшний день моделей образование минимума, обнаруженного в данных SOHO/SWAN, возможно только при учете многокомпонентного характера плазмы межзвездной среды.

6. Проведено моделирование полных карт неба в интенсивности рассеянного Лайман-альфа излучения и сравнение с данными SOHO/SWAN. Моделирование проводилось для двух имеющихся на сегодняшний день способов задания временных и гелиоширотных вариаций потока массы солнечного ветра. Показано, что пространственное распределение интенсивности рассеянного солнечного Лайман-альфа излучения на 1 а.е. определяется зависимостью частоты ионизации атомов водорода от гелиошироты. Это означает, что экспериментальные данные по Лайман-альфа излучению можно использовать для диагностики зависимости параметров солнечного ветра от широты.

7. Показано, что разработанная трехмерная нестационарная модель с учетом гелиоширотной зависимости параметров солнечного ветра, основанной на результатах Quemerais et al. (2006а), приводит к хорошему совпадению с данными SOHO/SWAN по пространственному распределению интенсивности рассеянного Лайман-альфа излучения. В то же время, результаты аналогичной модели с учетом гелиоширотной зависимости параметров солнечного ветра, основанной на результатах Sokol et al. (2012), противоречат данным SOHO/SWAN в период максимума солнечной активности (2002-2005 гг.).

ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ

1. Katushkina O.A., Izmodenov V.V., Quemerais Е., Sokol J.M., Heliolatitudinal and time variations of the solar wind mass flux: inferences from the solar Lyman-alpha intensity maps // Journal of Geophysical Research. -2013. -V. 118. -P. 2800-2808.

2. Катушкина O.A., Измоденов B.B., Влияние эффектов гелиосферного ударного слоя на параметры рассеянного солнечного Лайман-альфа излучения,

измеряемого на орбите Земли // Космические исследования. -2012. -Т.50. -№ 2. -С.147-157.

3. Katushkina О.A., Izmodenov V.V., Spectral properties of backscattered solar Ly-alpha radiation in the heliosphere: a theoretical search of the heliospheric boundaries effects // Advances in Space Research. -2011. -V. 48. -Is. 12. -C. 1967-1979.

4. Катушкина O.A., Измоденов B.B., Влияние эффектов гелиосферного интерфейса на распределение параметров атомов межзвездного водорода внутри гелиосферы // Письма в Астрон. Ж. -2010. -Т. 36. -№ 4. -С.310-319.

5. Izmodenov V.V., Katushkina О.А., Quemerais Е., Bzowski М., Distribution of Interstellar Hydrogen Atoms in the Heliosphere and Backscattered Solar Lyman-alpha // chapter 2 in "Cross-Calibration of Past and Present Far UV Spectra of Solar System Objects and the Heliosphere", ISSI Scientific Report Series, ed. R.M. Bonnet, E. Quemerais, M. Snow. -2013. -V.13. -P.7-65. Springer.

6. Катушкина O.A., Кинетические эффекты в распределении межзвездных атомов водорода в межпланетном пространстве: анализ данных по рассеянному лайман-альфа излучению // Вестник Нижегородского государственного университета им. Н.И. Лобачевского. №4. -часть 5 (сборник трудов X Всероссийского съезда по фундаментальным проблемам теоретической и прикладной механики). -2011. -С.313-315.

7. Катушкина О.А., О спектральных свойствах рассеянного солнечного лайман-альфа излучения: теоретическое исследование эффектов, связанных с границей гелиосферы // Сборник трудов 4-ой Всероссийской школы-семинара "Аэрофизика и физическая механика классических и квантовых систем", ИПМех РАН.-2011.-С.71-77.

8. Катушкина О.А., Измоденов В.В., Теоретическое моделирование движения атомов водорода из межзвездной среды внутрь гелиосферы // Сборник трудов XIII Школы молодых ученых "Актуальные проблемы физики", ФИАН. -2010. -С.128-129.

ЦИТИРУЕМАЯ ЛИТЕРАТУРА

1. Shklovsky I.S., On hydrogen emission in the night glow // Planet. Space Sci. -1959. -V.l. -P.63-65.

2. Курт В.Г., Гермагенова T.A., Рассеяние солнечного Ly-a излучения на галактическом водороде // Астрон. Ж. -1967. -Т.44. -№2. -С. 352-357.

3. Курт В.Г., Сюняев Р.А., Наблюдения и интерпретация ультрафиолетового излучения галактики // Астрон. Ж. -1967. -Т.44. -№6. -С. 1157-1161.

4. Fahr H.J., The Interplanetary Hydrogen Cone and its Solar Cycle Variations // Astronomy and Astrophysics. -1971.-V.14. -P.263.

5. Lallement, R., Bertaux, J.L., Kurt, V.G., and Mironova, E.N. Observed perturbations of the velocity distribution of interstellar H atoms in the solar system with Prognoz Lyman a measurements // Astron. Astrophys. -1984. -V.140. -P.243-250.

6. Lallement, R., Bertaux, J.L., and Kurt, V.G. Solar wind decrease at high heliographic latitudes detected from Prognoz interplanetary Lyman-alpha mapping // J. Geophys. Res.-1985.-V.90.-P.1413-1423.

7. Baranov V.B., Malama Yu.G., Model of the solar wind interaction with the local interstellar medium: numerical solution of self-consistent problem // J. Geophys. Res. -1993. -V. 98. -No. A9. -P.15,157-15,163.

8. Izmodenov V.V., Malama Yu.G., Ruderman M.S., Chalov S.V., Alexashov D.B., Katushkina O.A., Provornikova E.A., Kinetic-gasdynamic modeling of the heliospheric interface // Space Science Reviews. -2009. -V. 146. -Is. 1-4. -P. 329-351.

9. Danby J.M.A., and Camm G.L., Statistical dynamics and accretion // Monthly Notices of the Royal Astron. Soc. -1957. -V.l 17. -P.50-71.

10.Meier R.R., Some optical and kinetic properties of the nearby interstellar gas // Astronomy and Astrophysics.-1977. -V. 55. -No. 2.-P. 211-219.

11 .Wu F.M., and Judge D.L., Temperature and flow velocity of the interplanetary gases along solar radii // Astrophys. J. -1979. -V.231. -P. 594-605.

12.Quemerais, E., et al., Interplanetary hydrogen absolute ionization rates: Retrieving the solar wind mass flux latitude and cycle dependence with SWAN/SOHO maps // J. Geophys. Res. -2006a. -V.l 11. -P.A09114.

13.Costa J., Lallement R., Quemerais E. Heliospheric interstellar H temperature from SOHO/SWAN H cell data // Astron. Astrophys. -1999. -V. 349. -P. 660-672.

14.Quemerais E., Lallement R., Bertaux J.-L., et al. Interplanetary Lyman-alpha line profiles: variations with solar activity cycle // Astronomy and Astrophysics. -2006b. -V.455.-ÎS.3.-P.1135-1142.

15.Malama Y.G., Izmodenov V.V., Chalov S.V., New model of the heliospheric interface: three-fluid plasma approximation // Astron. Astrophys. -2006. -V.445. -P.693-701.

lô.Sokol J. M., et al., Heliolatitude and Time Variations of Solar Wind Structure from in situ Measurements and Interplanetary Scintillation Observations // Solar Physics. -2012.-doi:10.1007/sl 1207-012-9993-9.

Подписано в печать: 01.10.2013 Объем: 1,0 п. л. Тираж: 100 экз. Заказ № 165 Отпечатано в типографии «Реглет» 119526, г. Москва, пр-т Вернадского, д. 39 (495) 363-78-90; www.reglet.ru

 
Текст научной работы диссертации и автореферата по астрономии, кандидата физико-математических наук, Катушкина, Ольга Александровна, Москва

Московский государственный университет имени М.В.Ломоносова Механико-математический факультет

На правах рукописи

ПХ ЧП4 -lJ.lt и Л иТ4.и I I

Катушкина Ольга Александровна

Кинетическое моделирование распределения межзвездных атомов водорода в гелиосфере и анализ экспериментальных данных 8(ЭНО/8ЛОД1Ч по рассеянному солнечному Лайман-альфа излучению

01.03.03 — физика Солнца 01.02.05 — механика жидкости, газа и плазмы

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель — д.ф.-м.н., профессор Измоденов Владислав Валерьевич

Москва 2013

Оглавление

Введение: общая характеристика работы 4

Определения, сокращения и обозначения 11

1 Кинетическая модель распределения межзвездных атомов

водорода в гелиосфере 13

1.1 Введение к главе 1........................ 13

1.1.1 Первые идеи о проникновении межзвездных атомов водорода в окрестность Солнца............. 13

1.1.2 Аналитическое решение для функции распределения атомов по скоростям в случае классической горячей модели........................... 17

1.1.3 Современные представления о структуре гелиосфер-ного ударного слоя.................... 19

1.1.4 Современное развитие горячей модели......... 25

1.2 Распределение атомов водорода по пространству и скоростям

в гелиосфере в осесимметричном стационарном случае ... 28

1.2.1 Постановка задачи: уравнения............. 28

1.2.2 Упрощение правой части уравнения Больцмана ... 30

1.2.3 Общие выражения для частоты перезарядки и фотоионизации ......................... 31

1.2.4 Кинетические эффекты в распределении межзвездных атомов на расстоянии 90 а.е. от Солнца..... 32

1.2.5 Постановка задачи: граничные условия........ 37

1.2.6 Метод решения...................... 43

1.2.7 Результаты: влияние эффектов гелиосферного ударного слоя на распределение атомов водорода внутри гелиосферы........................ 44

1.3 Распределение атомов водорода по пространству и скоростям

в трехмерном нестационарном случае............. 48

1.3.1 Постановка задачи: уравнения и граничные условия . 48

1.3.2 Задание определяющих параметров Рех,Е^,Х), РрИ.Е^, и \ щ) в трехмерном нестационарном случае........................... 51

1.3.3 Результаты численного моделирования: параметры атомов водорода в гелиосфере ............. 54

1.4 Выводы к главе 1......................... 63

2 Теоретическое исследование спектров рассеянного солнечного Лайман-альфа излучения в гелиосфере 65

2.1 Введение к главе 2........................ 65

2.2 Постановка задачи о нахождении спектра рассеянного солнечного Лайман-альфа излучения............... 72

2.2.1 Получение формального аналитического решения

уравнения переноса излучения............. 80

2.3 Результаты............................ 82

2.3.1 Влияние эффектов гелиосферного ударного слоя на параметры рассеянного солнечного Лайман-альфа излучения на орбите Земли................ 82

2.3.2 Исследование зависимости ширины линии рассеянного Лайман-альфа излучения от направления луча зрения наблюдателя..................... 92

2.3.3 Исследование зависимости моментов спектра рассеянного излучения от значений параметров ц, и Рм,е • • НО

2.4 Выводы к главе 2......................... 112

3 Рассеянное солнечное Лайман-альфа излучение как способ

диагностики гелиоширотных и временных вариаций потока

массы солнечного ветра 114

3.1 Введение.............................. 114

3.2 Два способа задания частоты перезарядки Рех^, А) .....116

3.3 Результаты: карты интенсивности рассеянного солнечного Лайман-альфа излучения, сравнение с данными ЗОНО/8\¥А]М120

3.4 Заключение к главе 3 ...................... 130

Заключение 132

Литература

135

Введение: общая характеристика работы

Актуальность работы

В диссертации проводится теоретический анализ эволюции функции распределения межзвездных атомов водорода по скоростям внутри гелио-сферы - области космического пространства, занятой солнечным ветром. Разработана трехмерная нестационарная кинетическая модель распределения атомов водорода, в которой учитываются эффекты, связанные 1) с 11-летним циклом солнечной активности; 2) с зависимостью параметров солнечного ветра от гелиошироты; 3) с изменением функции распределения межзвездных атомов при их прохождении через область взаимодействия солнечного ветра с межзвездной средой. Построенная численная модель применяется для анализа экспериментальных данных по рассеянному солнечному Лайман-альфа излучению, полученных на космическом аппарате (КА) SOHO (прибором SWAN). Актуальность данной работы обусловлена, во-первых, интересом к исследованию границы гелиосферы, возросшим в последнее время благодаря появлению новых уникальных данных с КА Вояджер-1 и 2, а также IBEX - Interstellar Boundary Explorer, и, во-вторых, большим количеством экспериментальной информации о рассеянном солнечном Лайман-альфа излучении, накопленной за последние 25 лет. Исторически, именно данные о рассеянном Лайман-альфа излучении привели к идее о существовании межзвездного ветра, а также к развитию теоретических моделей проникновения межзвездных нейтральных атомов внутрь гелиосферы. В самом начале космической эры (Shklovsky, 1959) при первых ракетных запусках было обнаружено излучение в ультрафиолетовом (УФ) диапазоне, которое невозможно было объяснить излучением от ближайших к Солнцу звезд или рассеянием солнечного излучения на нейтральных атомах в атмосфере Земли. В работах Курта & Гермогеновой (1967) и Курта & Сюняева (1967) была высказана гипотеза о галактическом происхождении измеряемого излучения. В дальнейшем более детальные измерения на аппарате OGO-5 на основании обнаруженного параллакс-эффекта показали, что единственным возможным источником рассеянного Лайман-альфа излучения могут быть межзвездные атомы водорода, которые проникли

внутрь Солнечной системы за счет относительного движения Солнца и окружающей его межзвездной среды. Такое явление получило название межзвездного ветра. Таким образом, межзвездные атомы водорода проникают из Локальной Межзвездной Среды (ЛМС) внутрь гелиосферы, где они могут быть измерены с помощью прямых или косвенных методов. Внутри гелиосферы распределение атомов водорода является источником информации как о параметрах ЛМС, так и о свойствах переходной области между солнечным и межзвездным ветрами (в литературе эту область называют гелиосферным ударным слоем). Кроме того, солнечное Лайман-альфа излучение рассеивается на межзвездных атомах водорода, и поэтому спектральные свойства рассеянного излучения определяются функцией распределения атомов водорода по скоростям. Следовательно, измерения рассеянного солнечного Лайман-альфа излучения открывают возможности для удаленной диагностики распределения межзвездных атомов водорода в гелиосфере. Начиная с 1980-х годов прошлого века основной теоретической моделью распределения атомов водорода в гелиосфере являлась, так называемая, классическая горячая модель (Fahr, 1971), которая использовалась в течение многих лет для анализа и интерпретации данных по рассеянному солнечному Лайман-альфа излучению (Lallement et. al., 1984, 1985 и др.). Тем не менее, эта модель обладает рядом недостатков. А именно, в ней не учитываются нестационарные и трехмерные эффекты, связанные с 11-летнем циклом солнечной активности и с зависимостью параметров солнечного ветра от гелиошироты, а также не учитывается существование области гелиосферного ударного слоя. В диссертационной работе разработана (глава 1) новая кинетическая модель распределения атомов водорода в гелиосфере, которая, с одной стороны, позволяет избавиться от перечисленных недостатков классической горячей модели, а с другой стороны, является достаточно простой и эффективной с вычислительной точки зрения. Разработанная модель применяется для исследования спектральных характеристик рассеянного солнечного Лайман-альфа излучения в гелиосфере (глава 2), а также для анализа и интерпретации экспериментальных данных, полученных на космическом аппарате SOHO (прибор SWAN) (глава 3).

Цели и задачи работы состоят в следующем: 1. Разработать трехмерную нестационарную кинетическую модель рас-

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

2. Исследовать, как немаксвелловские свойства функции распределения атомов водорода во внешней гелиосфере (на 90 а.е. от Солнца) влияют на распределение (по пространству и скоростям) атомов водорода вблизи Солнца. Немаксвелловский характер функции распределения во внешней гелиосфере связан с процессом перезарядки в области гелиосферного ударного слоя.

3. Проанализировать спектральные характеристики рассеянного солнечного Лайман-альфа излучения в гелиосфере. Определить, как немаксвелловские свойства в распределении атомов водорода во внешней гелиосфере влияют на спектральные характеристики рассеянного Лайман-альфа излучения.

4. Провести теоретическое моделирование интенсивности рассеянного солнечного Лайман-альфа излучения и сравнить результаты теоретических расчетов с экспериментальными данными, полученными на КА SOHO (прибором SWAN).

5. Провести исследование зависимости потока массы солнечного ветра от гелиошироты на основании полных карт неба в интенсивности рассеянного излучения, полученных из данных SOHO/SWAN.

Научная новизна работы

Впервые была разработана детальная кинетическая модель распределения атомов водорода в гелиосфере, которая учитывает одновременно локальные трехмерные и нестационарные эффекты, а также кинетические эффекты, связанные с перезарядкой на границе гелиосферы. На основе этой модели были проанализированы данные SOHO/SWAN.

Автором впервые получены и выносятся на защиту следующие основные положения:

I. Впервые классифицированы кинетические, немаксвелловские свойства функции распределения атомов водорода на расстоянии 90 а.е. от Солнца, возникающие за счет процесса перезарядки межзвездных атомов на протонах ЛМС и солнечного ветра в области гелиосферного ударного слоя.

2. Разработана трехмерная нестационарная кинетическая модель распределения межзвездных атомов водорода внутри гелиосферы с учетом эффектов, связанных с изменением функции распределения атомов в области гелиосферного ударного слоя. Данная модель является эффективным инструментом для анализа многочисленных экспериментальных данных.

3. Исследовано влияние гелиосферного ударного слоя на моменты функции распределения атомов водорода вблизи Солнца, а также на характеристики рассеянного солнечного Лайман-альфа излучения. Показано, что наиболее чувствительными к немаксвелловским свойствам функции распределения атомов во внешней гелиосфере являются кинетическая температура атомов и ширина линии рассеянного Лайман-альфа излучения.

4. Проведен детальный теоретический анализ зависимости ширины линии рассеянного Лайман-альфа излучения от направления луча зрения наблюдателя. Показано, что немонотонное поведение ширины линии, которое ранее было обнаружено в данных ЗОНО/ЗШАИ, может быть частично объяснено в рамках текущих моделей только с учетом многокомпонентного характера плазмы межзвездной среды.

5. На основании разработанной в настоящей диссертации трехмерной и нестационарной модели распределения атомов водорода в гелиосфере впервые был проведен анализ полных карт неба в интенсивности рассеянного Лайман-альфа излучения, полученных на КА SOHO/SWAN в 1996-2009 гг. При этом впервые в модели использовалась реальная зависимость параметров солнечного ветра от гелиошироты, основанная на различных экспериментальных данных.

6. Установлено, что для объяснения данных SOHO/SWAN по интенсивности рассеянного Лайман-альфа излучения в период максимума солнечной активности (с 2001 по 2005 гг.) необходимо предположить, что график зависимости потока массы солнечного ветра от гелиошироты имеет два максимума на средних широтах (±50 — 60°). Информация о гелиоширотных вариациях потока массы солнечного ветра имеет принципиальное значение для глобального моделирования гелиосферного ударного слоя и исследования процессов, происходящих

на границе гелиосферы. Достоверность результатов

Достоверность представленных в диссертации результатов, основана на следующем:

1. Для моделирования распределения атомов водорода использовался наиболее адекватный в данном случае кинетический подход. Необходимость использования кинетического подхода связана с большой длиной свободного пробега атомов водорода, сравнимой с характерным размером гелиосферы. Также для задания в модели граничного условия на расстоянии 90 а.е. от Солнца использовались результаты глобальных самосогласованных кинетико-газодинамических моделей взаимодействия солнечного ветра с межзвездной средой (Baranov & Malama, 1993; Izmodenov et al., 2009), которые отлично зарекомендовали себя ранее тем, что с их помощью удается объяснить множество экспериментальных фактов.

2. Разработанный численный алгоритм тестировался на классических аналитических решениях в стационарном осесимметричном случае. Также проводилось дополнительное тестирование программы в нестационарном трехмерном случае с целью выявления наиболее оптимальных параметров расчетных сеток.

3. Теоретические результаты, полученные на основании разработанной в диссертации модели, хорошо согласуются с экспериментальными данными SOHO/SWAN по рассеянному Лайман-альфа излучению.

Практическая ценность

Практическая ценность работы заключается в том, что разработанную в диссертации модель можно эффективно применять для анализа экспериментальных данных по измерениям потоков межзвездных нейтральных атомов (водорода, гелия и др.), а также параметров рассеянного на атомах солнечного излучения с уже существующих и будущих космических аппаратов. В частности, эта модель представляет собой удобный инструмент для анализа данных по потокам межзвездных атомов водорода и гелия, измеряемых с 2009 г. на КА IBEX прибором IBEX-Lo. Кроме того, модель можно использовать для определения научных задач будущих российских и зарубежных космических миссий. В частности, в настоящее

время результаты, полученные в диссертации, используются для составления научной программы прибора "Лаймус", разрабатываемого в ИКИ РАН. Задачей прибора является измерение интенсивности и спектра рассеянного солнечного Лайман-альфа излучения. В случае успешного запуска прибора, разработанная в диссертации модель будет также использоваться и для анализа полученных на "Лаймусе" данных.

Апробация работы

Работы, вошедшие в диссертацию, докладывались и обсуждались на научно-исследовательских семинарах по аэромеханике и газовой динамике механико-математического факультета МГУ им. М.В. Ломоносова (рук. -акад. Черный Г.Г.), семинарах Института космических исследований РАН (рук. - академик Зеленый Л.М.), семинарах лаборатории физической газовой динамики Института проблем механики им. А.Ю. Ишлинского РАН (рук.- проф. Баранов В.В.). Основные положения и результаты, вошедшие в диссертацию, докладывались на российских и международных конференциях, в том числе:

• на V-IX конференциях молодых ученых "Фундаментальные и прикладные космические исследования "(ИКИ РАН, Москва, 2009-2012 гг.);

• на Международном научном форуме студентов, аспирантов и молодых ученых "Ломоносов"(МГУ, Москва, 2009-2011 гг.);

• на конференции Ломоносовские чтения (МГУ, Москва, 2010 г.);

• на конференции "Физика плазмы в солнечной системе"(ИКИ РАН, Москва, 2009-2011 гг.);

• на Всероссийской школе-семинаре "Аэрофизика и физическая механика классических и квантовых систем"(ИПМех РАН, Москва, 20102011 гг.);

• на X Всероссийском съезде по фундаментальным проблемам теоретической и прикладной механики (Нижний Новгород, 2011 г.);

• на 18-й конференции молодых ученых по астрономии и космической физике (Национальный университет им. Тараса Шевченко, Киев, Украина, 2011 г.);

• на Всероссийской ежегодной конференции по физике Солнца "Солнечная и солнечно-земная физика - 2008" (Главная астрономическая обсерватория РАН, Санкт-Петербург, 2008 г.)

• на ежегодных конференциях Европейского геофизического общества (EGU) (г. Вена, Австрия, 2010, 2012 гг.);

• на 12-й международной конференции "Солнечный ветер"(г. Сан-Мало, Франция, 2009);

• на 11-й Научной Ассамблее Ассоциации Геомагнетизма и Аэрономии (IAGA) (г. Сопрон, Венгрия, 2009 г.);

• на 37-й, 38-й и 39-й научных ассамблеях COSPAR (г. Монреаль, Канада, 2008; г. Бремен, Германия, 2010; г. Майсор, Индия, 2012);

• на XXVIII Генеральной Ассамблее Международного астрономического общества (IAU) (г. Пекин, Китай, 2012 г.).

Публикации и личный вклад автора

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