Двухчастичные и многочастичные статистические модели потока солнечной плазмы
На правах рукописи
Минькова Наталья Романовна ДВУХЧАСТИЧНЫЕ И МНОГОЧАСТИЧНЫЕ СТАТИСТИЧЕСКИЕ МОДЕЛИ ПОТОКА СОЛНЕЧНОЙ ПЛАЗМЫ Специальность – 01.02.05 – механика жидкости, газа и плазмы
Автореферат диссертации на соискание ученой степени кандидата физико-математических наук
Томск - 2005 2
Работа выполнена на кафедре математической физики физико технического факультета Томского государственного университета.
Научный консультант: доктор физико-математических наук, профессор Васенин Игорь Михайлович
Официальные оппоненты: доктор физико-математических наук, профессор Глазунов Анатолий Алексеевич доктор физико-математических наук, профессор Рудяк Валерий Яковлевич
Ведущая организация: Институт прикладной механики, УрО РАН.
Защита состоится 23 декабря 2005 года в 14.30 часов на заседании диссертационного совета Д 212.267.13 при Томском государственном университете по адресу: 634050, г. Томск, пр. Ленина, 36, ауд 503, корпус №10 (НИИ ПММ).
С диссертацией можно ознакомиться в Научной библиотеке Томского государственного университета по адресу: г. Томск, пр. Ленина, 34а.
Автореферат разослан 17 ноября дата
Ученый секретарь диссертационного совета д.т.н. _ /Ю.Ф. Христенко/
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
.
Актуальность темы. В настоящее время активно развиваются астрофизические приложения механики газа и плазмы. В частности, многие исследования, проводимые в этом направлении, посвящены моделированию расширяющегося потока плазмы, испускаемого Солнцем и получившего название солнечного ветра. Солнечный ветер характеризуется резким изменением параметров с увеличением гелиоцентрического расстояния. Его числовая плотность и температура изменяются в периоды низкой солнечной активности соответственно от значений порядка 1014 м-3 и 106K у основания короны до 106 м-3 и 105 K на уровне орбиты Земли. При этом скорость потока плазмы возрастает от почти нулевых значений до величин порядка 400- км/с. Вспышки на Солнце сопровождаются увеличением скорости солнечного ветра в 2-3 раза и соответствующим уменьшением его плотности.
Такой характер потока солнечной плазмы ограничивает возможности приближения сплошной среды при его описании и приводит к необходимости применения методов статистической механики.
Актуальность изучения солнечного ветра обусловлена его воздействием на магнитосферу Земли. Возникающие в результате взаимодействия быстрого потока солнечной плазмы с геомагнитным полем возникают магнитные бури, оказывающие влияние на работу навигационных систем (в частности, известны случаи потери связи и управления околоземными спутниками), электрических и магнитных приборов, линий электропередачи, здоровье людей, урожайность сельскохозяйственных культур и т.д. Необходимо учитывать также влияние солнечного ветра на бортовые системы космических аппаратов в межпланетном пространстве, на проводимые на их борту измерения. Все это обуславливает важность и актуальность изучения природы солнечного ветра, основных его закономерностей и прогнозирования его параметров на основе адекватных моделей солнечной плазмы.
Кроме того, изучение Солнца, наиболее близкой к Земле звезды и наиболее доступной для наблюдения, имеет большое фундаментальное значение, поскольку способствует развитию теории звезд, звездного и межзвездного ветра.
Современная теория солнечного ветра за последние 40 лет достигла значительных успехов в моделировании процессов, происходящих в гелиосфере. Однако несмотря на многочисленность исследований факторов, влияющих на солнечный ветер, остается актуальным вопрос о природе ускорения солнечного ветра и его прогнозировании.
Цель и задачи исследования. Теоретическое исследование распространения плазмы солнечного ветра в широком диапазоне значений гелиоцентрического расстояния. Построение двухчастичных кинетических моделей солнечного ветра в приближении сферически симметричного течения квазинейтральной бесстолкновительной водородной плазмы без учета влияния магнитного поля. Получение в рамках этих моделей аналитических зависимостей средних параметров потока от гелиоцентрического расстояния и анализ возможности ускорения солнечного ветра за счет теплового движения частиц высокотемпературной плазмы у основания короны Солнца.
Разработка многочастичной статистической модели потока газа и плазмы, основанной на теореме Лиувилля, и методологии вывода в рамках этой модели среднестатистических характеристик такого потока. Получение с использованием данной модели аналитических выражений для числовой плотности и скорости плазмы солнечного ветра. Интерпретация двухчастичных моделей с точки зрения многочастичного подхода.
Объект и предмет исследования. Объектом исследования диссертационной работы является конкретное природное явление - поток бесстолкновительной плазмы. Предмет исследования – анализ возможностей двух и более частичных статистических подходов к описанию бесстолкновительной плазмы, в частности к моделированию ускорения ее сферически расширяющегося потока.
Методы исследования и фактический материал. Проведенное исследование основывается на методах статистической механики и методах решения дифференциальных уравнений в частных производных первого порядка.
Фактические материалы – данные астрофизических наблюдений, опубликованные в литературе.
Достоверность. Полученные аналитические решения для функции распределения частиц медленного и быстрого солнечного ветра по скоростям, для числовой плотности и скорости солнечного ветра сравниваются с известными теоретическими и наблюдательными данными, что демонстрирует их количественное согласование с последними.
Для проверки методики и результатов, представленных для многочастичной статистической модели, проведены решения тестовых классических задач, известных в литературе, а также сравнение с двухчастичной кинетической моделью.
Научная новизна работы заключается в следующем:
1. Построены двухчастичные кинетические модели бесстолкновительной полностью ионизованной водородной плазмы, которые применены к описанию солнечного ветра. В рамках этих моделей получены точные аналитические решения соответствующих кинетических уравнений для функций распределения вероятностей и аналитических пространственных зависимостей скорости и плотности солнечного ветра для случая стационарного сферически симметричного распространения квазинейтральной бесстолкновительной водородной плазмы в пренебрежении влиянием магнитного поля.
2. Показано в рамках рассматриваемых моделей, что ускорение потока солнечной плазмы в гравитационном поле Солнца обусловлено уменьшением доли частиц атмосферы (кинетическая энергия которых меньше гравитационного потенциального барьера) по сравнению с долей убегающих частиц (кинетической энергии которых достаточно для преодоления гравитационного потенциального барьера).
3. Построена кинетическая полуэмпирическая модель плазмы быстрого солнечного ветра, дающая результаты по плотности и скорости потока частиц, согласующиеся с данными наблюдений.
4. Построена многочастичная статистическая модель стационарного потока газа и плазмы, основанная на теореме Лиувилля, обобщенной на случай открытой системы с постоянным числом частиц, находящейся в динамическом равновесии с окружающей средой.
Разработана методика вычисления статистических моментов в рамках этой модели. Показано, что для стационарного сферически симметричного потока двухкомпонентной плазмы, истекающей со сферы заданного радиуса, многочастичный статистический подход, основанный на уравнении Лиувилля, сводится к двухчастичной кинетической модели, если движение частиц статистически независимо. Указаны частные случаи, когда понижение размерности функции распределения, описывающей течение плазмы (газа) возможно лишь до значений, больших числа компонент плазмы.
Теоретическая значимость полученных результатов. Впервые представлены двухчастичные кинетические модели сферически симметричного стационарного бесстолкновительного потока полностью ионизованной водородной плазмы. Их применение к описанию солнечного ветра позволило вывести аналитические пространственные зависимости для плотности и скорости солнечного ветра, не содержащие параметров согласования, которые в целом лучше согласуются с данными наблюдений по сравнению с одночастичными кинетическими моделями при тех же допущениях. Получены верхние и нижние оценки для профилей скорости и плотности на множестве реализаций поляризационного потенциала плазмы, обеспечивающих существование этих стационарных зависимостей в рамках рассматриваемых приближений. На основе указанных моделей показан механизм ускорения медленного солнечного ветра на фоне торможения испускаемых Солнцем частиц плазмы в его гравитационном поле. На основе представленной полуэмпирической двухчастичной кинетической модели потока плазмы в расширяющейся трубке тока исследовано влияние различных факторов на ускорение быстрого солнечного ветра.
Разработан многочастичный статистический подход к описанию потока газа и плазмы, основанный на уравнении Лиувилля, и предложена в его рамках методика вычисления статистических моментов для газа и полностью ионизованной многокомпонентной плазмы. Показано, что зависимости для плотности и скорости двухкомпонентного (электрон-протонного) солнечного ветра, полученные на основе данного подхода, совпадают с результатами двухчастичной кинетической модели для стационарного сферически симметричного бесстолкновительного потока плазмы.
Значимость работы для практики. Полученные аналитические результаты по радиальным профилям числовой плотности и скорости солнечного ветра позволяют в рамках принятых допущений прогнозировать изменение этих параметров при различных значениях гелиоцентрического расстояния в зависимости от наблюдаемых значений температуры и плотности плазмы у основания короны, а также восстанавливать последние по измерениям на других расстояниях от Солнца (например, вблизи земной орбиты, что более доступно для различных измерительных комплексов).
Результаты настоящей работы включены в отчеты по фундаментальным научным исследованиям Томского государственного университета «Исследования по математике и моделям естественно протекающих процессов» и по гранту РФФИ 01-01-00983, а также могут быть использованы в учебном процессе.
Апробация работы. Материалы, представленные в диссертации, докладывались и получили признание на 23 конференциях, в том числе:
1. Международная конференция «Математические модели и методы их исследования» (Красноярск, 2001).
2. Пятая Сибирская школа-семинар «Математические проблемы механики сплошных сред» (Новосибирск, 2001).
3. IV Всероссийский семинар «Моделирование неравновесных систем 2001», (Красноярск, 2001).
4. Конференция «Проблемы и перспективы технологий атомной промышленности». (Томск, 2001).
5. 11-ый Международный конгресс по физике плазмы. (Сидней, Австралия, 2002 г.) 6. Международная конференция Strongly Coupled Coulomb Systems (Santa Fe, New Mexico, USA 2002 г.) 7. IV Международная научно-практическая конференция «Фундаментальные и прикладные проблемы приборостроения, информатики, экономики и права» (Сочи, 2002) 8. V всероссийский семинар «Моделирование неравновесных систем 2002» (Красноярск, 2002 г.) 9. Конференции стран СНГ и Прибалтики «Актуальные проблемы физики солнечной и звездной активности» (Нижний Новгород, 2003 года).
10. International Union of Geodesy and Geophysics (IUGG) General Assembly (Sapporo, Japan, 2003).
11. VI Всероссийский семинар “Моделирование неравновесных систем” (Красноярск, 2003 г.) 12. XXXI Звенигородская конференция по физике плазмы и управляемому термоядерному синтезу (г. Звенигород, Московской обл., 2004 г.).
13. Всероссийская Астрономическая Конференция ВАК-2004 «Горизонты Вселенной» (Москва, 2004 г.).
14. XXXII Звенигородская конференция по физике плазмы и управляемому термоядерному синтезу (г. Звенигород, Московской обл., 2005 г.).
15. Solar Wind 11 - SOHO 16 «Connecting Sun and Heliosphere». (Whistler, Canada, 2005).
16. International Conference on Strongly Coupled Coulomb Systems (SCCS 2005) (Moscow, 2005).
Основные результаты диссертации опубликованы в 17 тезисах и статьях и докладах, список которых приведен в конце автореферата.
На защиту выносятся следующие положения:
1. Новые двухчастичные кинетические модели солнечного ветра, построенные в приближении стационарного сферически симметричного бесстолкновительного потока квзаинейтральной водородной плазмы (без учета влияния магнитного поля).
2. Новые аналитические результаты по зависимости числовой плотности и скорости медленного и быстрого солнечного ветра от гелиоцентрического расстояния.
3. Объяснение, в рамках представленных моделей, ускорения потока солнечной плазмы в гравитационном поле Солнца уменьшением доли частиц атмосферы по сравнению с долей убегающих частиц.
4. Новый многочастичный статистический подход к моделированию потока частиц газа и плазмы, основанный на теореме Лиувилля, обощенной на случай системы с постоянным числом частиц, находящейся в динамическом равновесии с окружающей средой;
новая методология расчета (в рамках этой модели) осредненных характеристик потока частиц в предположении о неразличимости частиц в объеме, масштаб которого обусловлен разрешением проводимых измерений или принятыми модельными допущениями.
5. Новые аналитические результаты по зависимости числовой плотности и скорости потока газа и плазмы, полученные на основе указанного выше мнгочастичного подхода.
6. Эквивалентность многочастичной модели, основанной на уравнении Лиувилля, и двухчастичной кинетической модели для стационарного сферически симметричного потока двухкомпонентной квазинейтральной бесстолкновительной плазмы, истекающей со сферы заданного радиуса.
7. Понижение размерности многочастичной статистической модели описывающей стационарное течение k-компонентной плазмы (газа) до значения, равного числу ее компонент k, в случае задания одинакового граничного условия для всех частиц потока, движение которых статистически независимо.
Структура и объем диссертации. Диссертация состоит из введения, четырех глав и заключения, одного приложения и списка использованной литературы. Работа изложена на 131 странице, содержит 13 рисунков, одну таблицу, список литературы включает 145 наименований.
СОДЕРЖАНИЕ РАБОТЫ
Во введении сформулированы актуальность работы, цель и задачи исследования, их научная новизна, практическая и теоретическая значимость результатов, положения, выносимые на защиту, а также приведен список конференций, на которых были представлены результаты работы.
В главе 1 проведен обзор исследований по теме диссертации, затрагивающий избранные вопросы современного состояния теории солнечного ветра. Настоящий обзор отражает в основном развитие кинетической концепции солнечного ветра в соответствии с темой диссертации и не претендует на исчерпывающий характер (по причине обширности литературы в данной области исследований). Дана краткая характеристика развития кинетических теорий в приближении статической короны Солнца в 1920-ые гг., когда было получено поле поляризации гидростатически равновесной плазмы Паннекока-Росселанда, разработаны модели Чепменом, Гунном, Коулингом для нижних слоев солнечной атмосферы. Ряд исследователей с конца XIX века развивали также динамические модели солнечной короны как совокупности потоков частиц, основываясь на многочисленных визуальных наблюдениях солнечных затмений, демонстрирующих существование таких потоков (Шэберле, Бредихин в XIX в., численные расчеты Штёрмера движения частиц короны с учетом воздействия дипольного магнитного поля в 1910-1030-ых гг., теория Альвена ускорения заряженных частиц в неоднородном магнитном поле в 1940 г., модель Кипенхойера солнечной короны как совокупности газовых шаров).
В 1950-ые гг. на основе наблюдений, моделей статической равновесной атмосферы Солнца и механической теории короны формируются представления о солнечном ветре (непрерывно расширяющейся короне).
Характеристика современной теории солнечного ветра начинается с работ, обсуждавших вопрос о роли электростатического поля, которое порождается разделением зарядов в плазме, возникающим вследствие значительно бльших скоростей легких электронов по сравнению с тяжелыми протонами (Северный, Альвен, Ван де Хулст, Пикельнер). Современный этап в развитии теории солнечного ветра открыли в конце 1950-ых гг. работы Паркера, базирующиеся на приближении сплошной среды и рассматривающие солнечную атмосферу, расширяющуюся непрерывно до сверхзвуковых скоростей, в противоположность гидростатически равновесной модели Чепмена. Для достижения наблюдаемых значений скорости и плотности солнечного ветра на орбите Земли Паркер фактически ввел дополнительный источник энергии, предположив изотермическое течение плазмы вблизи Солнца. Следует подчеркнуть, что указанная проблема (дополнительных энергетических источников), вообще, характерна для гидродинамических моделей (основанных на допущении сплошности среды). Это объясняется тем, что плотность солнечной плазмы падает экспоненциально с характерным пространственным масштабом ~0.1R от значений ~1014 м-3 у основания короны (число Кнудсена Kn1) до значений ~103 м-3 на орбите Земли, удаленной от Солнца на расстояние в 1 а.е.215R (Kn1). Для описания таких течений (с нарушением сплошности среды) используются кинетические и смешанные подходы. Развитие кинетических теорий медленного солнечного ветра представлено работами Чемберлена 1960 г., Оепика и Зингера 1959-1961 гг., Йенсена 1963 г., Джокерса 1970 г., Лемера и Шерера 1970 г., Пиеррард и др. 2003 г. и некоторыми другими работами. К началу 1970-ых гг. проведенные исследования показали, что как гидродинамический, так и одночастичный кинетический подходы позволяют моделировать ускоряющийся поток солнечной плазмы. Однако скорость медленного солнечного ветра, рассчитанная по соответствующим моделям при задании наблюдаемых значения температуры и плотности потока плазмы у основания короны, оказывается заниженной по сравнению с данными измерений на орбите Земли, если, конечно, не вводятся явно или неявно дополнительные источники энергии.
Исследованию возможных механизмов ускорения солнечного ветра посвящены многочисленные публикации (новый импульс работе в этом направлении дали данные наблюдений, полученные в 1990-ые гг. на борту космического аппарата Ulysses для быстрого солнечного ветра, а также исследования в рамках научной программы №16 отделения физических наук РАН «Плазменные процессы в солнечной системе»). В частности, рассматривается подпитка среднего движения плазмы за счет энергии МГД турбулентности, альвеновских волн, ионно-циклотронного резонанса, процессов пересоединения магнитных силовых линий, ударных волн.
Специфические условия истечения быстрого солнечного ветра из корональных дыр - а именно, расширяющиеся магнитные трубки (и соответственно трубки тока), наличие фракций тяжелых ионов с высокой температурой, повышенная концентрация высокоскоростных частиц и другие - учитываются во многих гидродинамических и кинетических моделях, опубликованных в последнее десятилетие. Полученные в их рамках результаты воспроизводят более резкое и более значительное изменение параметров быстрого солнечного ветра по сравнению с медленным солнечным ветром.
Несмотря на многочисленность исследований, учитывающих различные факторы, влияющие на солнечный ветер, вопрос о природе источников энергии, обеспечивающих повышение скорости потока солнечной плазмы до наблюдаемых величин, до сих пор остается открытым.
В главе 2 формулируется двухчастичная кинетическая модель бесстолкновительного стационарного сферически симметричного расширения квазинейтральной полностью ионизованной водородной плазмы, истекающей со сферической поверхности заданного радиуса и движущейся в поле гравитационных сил. Влиянием магнитного поля пренебрегается. В рамках этой модели находятся радиальные зависимости плотности и скорости потока частиц солнечного ветра, содержащие потенциал поляризации плазмы как неизвестную функцию. Получена область допустимых значений этого потенциала, при котором существуют указанные стационарные зависимости. Теоретические результаты сравниваются с данными наблюдений по медленному солнечному ветру.
На основе одночастичных кинетических уравнений для функций распределения вероятностей электронов fe и протонов fp по скоростям выводится двухчастичное кинетическое уравнение для функции распределения протонов и электронов f, которая в предположении статистической независимости движения электронов и протонов определяется как f= fe fp. В рамках указанных выше допущений это уравнение принимает следующий вид в сферической системе координат:
u p f f f f f u e + re + rp = 0, (1) r e p r u e r u p где r – гелиоцентрическое расстояние, ue, u p, u re, u rp - тангенциальные (модуль результирующего вектора меридиальной и азимутальной компонент) и радиальные компоненты скорости электронов и протонов, 2 2 2 e = m e (u re + u e ) / 2, p = m p (u rp + u p ) / 2 - кинетические энергии электрона и протона соответственно, me, mp,- массы электрона и протона, d d d d, re = e rp = e me mp -силы, действующие на протоны и dre dre dr p dr p электроны в поле поляризации плазмы с неизвестным потенциалом (r) и в гравитационном поле Солнца с потенциалом (r) = M r (M -масса Солнца, - гравитационная постоянная). Предположим, что сечение истечения находится у основания солнечной короны (r0~R, R – радиус Солнца), где состояние истекающей плазмы принимается равновесным:
2 m j ( u rj 0 + u j 0 ) 3/ mj 2 kT j f = f e 0 f p 0, f j 0 = f j (r0, u rj 0, u j 0 ) = 2 N j 0 e. (2) 2kT j Здесь urj, 0 u j, индексом 0 помечены значения переменных на границе r=r0. Начальные температуры электронов и протонов предполагаются одинаковыми Te0=Tp0=T0. Здесь индекс j=e,p обозначает соответственно параметры электронов и протонов. При этом fe0, fp нормированы таким образом, что их интегрирование по фазовому пространству начальных скоростей дает соответственно число электронов Ne и протонов Np0 в единице объема в сечении r=r0.
В модели не рассматриваются частицы, приходящие с бесконечности, и частицы, движущиеся по замкнутым траекториям, не выходящие на границы рассматриваемой области.
Решением задачи (1), (2), полученным методом характеристик, является функция распределения f = f (r, ure, ue, urp, u p ) :
m ( 0 ) (me m p ) 3 / 2 kT f = N0 e, (3) (2kT0 ) определенная в области:
D = {0 u e u e, ure ure ;
0 u p u p, urp urp }, (4) ) ( ) ( 2 rj j ( j - j ), ) ( 2 r u j = urj = r= j = m j ± e, где,, mj r 2 1 mj r j = j 0 j, j=e,p, m=me+mp,. Решение (3), (4) представлено на рис.1.
Рис.1.
Эволюция двухчастичной функции распределения f (3)-(4) вдоль радиального направления (каждый из трех графиков функции f представлен в своем относительном масштабе).
В рамках предположения о квазинейтральности плазмы (означающего равенство локальных значений числовой плотности электронов и протонов Ne=Np) и отсутствии тока (то есть при равенстве удельных потоков электронов и протонов Nue=Nup) следующие интегралы от двухчастичной функции распределения f представляют собой квадраты плотности плазмы N и ее удельного расхода Nu:
(Nu (r ) )2 = u re u rp f (r,U )dU, N 2 (r ) = (5) f ( r, U ) dU, D D где U обозначает совокупность переменных u кe, u e, u rp, u p, dU = 2u e du e du re 2u p du p du rp, D - область определения решения (4). На основе (5), (3), (4) выведены аналитические зависимости числовой плотности N(r), удельного потока частиц Nu(r) и скорости u (r ) = Nu (r ) N (r ) солнечного ветра от гелиоцентрического расстояния r, содержащие в качестве параметра потенциал поляризации плазмы (r), для определения которого необходимо введение дополнительной гипотезы. Получена область значений потенциала, при которых эти стационарные зависимости существуют:
e 0 1 e e 0 1 e + (1 ), 0 1, (6) m 0 r m m 0 r m r где 0=(r0), 0=(r0). В предположении гидростатического равновесия плазмы при r=r0 (приближение Паннекока Росселанда) e 0 = m 0 1 e неравенство (6) принимает вид:. Причем в области 2r m 2r допустимых значений потенциала поляризации (6) любое решение рассматриваемой задачи удовлетворяет неравенству:
N (r ) N (r ), N (r ) N (r ), Nu (r ) Nu (r ), u (r ) u (r ), (7) где N (r ) - числовая плотность убегающих частиц (кинетической энергии которых достаточно для преодоления потенциального барьера), звездочкой помечены зависимости, рассчитанные при условии e (r ) = m(r ) 2, соответствующем приближению Паннекока-Росселанда.
Предельное выражение для числовой плотности плазмы при rr N = 0.5 N 0 exp(0.50 ) (0.5 e 0 )(0.5 p 0 ) r 2 с учетом зависимости для Nu (r ) = 0.5 N 0 exp(0.50 ) u e 0 u p 0 (1 e 0 )(1 p 0 ) r - удельного потока приводит к следующему соотношению для предельного значения скорости солнечного ветра:
(1 e 0 )(1 p 0 ) u = u e0 u p 0, (8) (0.5 e 0 )(0.5 p 0 ) 0 = m 0 kT0, j 0 = m j 0 ± e 0, u j 0 = 2kT0 m j ( j = e, p ). В Здесь приближении Паннекока-Росселанда из формулы (8) получаем соотношение u = u e 0 u p 0 (2 0 ) (1 0 ), которое для наблюдаемых значений температуры T0=1.010 K при r0=1.5·R дает u 1.031 u e 0 u p 0 540 км / с.
На уровне орбиты Земли (ra.e215R) согласно полученным зависимостям скорость солнечного ветра ua.e равна примерно 460 км/c. Обе приведенные оценки для скорости солнечного ветра u и u согласуются с данными наблюдений. Заметим, что возрастание скорости u (r ) обусловлено уменьшением доли частиц атмосферы в потоке солнечной плазмы с удалением от Солнца. Частицы атмосферы играют роль балласта, поскольку их суммарный стационарный поток равен нулю, и, таким образом, они не дают вклад в общий поток плазмы Nu. Как следствие, u (r ) стремится к значению средней скорости убегающих частиц при r.
При r0=1.5R и T0=1.210 K, согласующихся с данными наблюдений, расхождение между теоретической зависимостью N (r ) и эмпирической зависимостью Рубцова N emp (r ) не превышает 8%. На рис.2 представлено сравнение теоретических результатов по скорости солнечного ветра U(r/R) при тех же значениях параметров r0 и T. (нижняя сплошная линия) с данными наблюдений, опубликованными в монографии Якубова, и другими моделями.
Представленные результаты демонстрируют адекватность предлагаемой модели для описания таких средних параметров солнечного ветра, как его числовая плотность, поток частиц и скорость. Они также позволяют заключить, что на больших гелиоцентрических расстояниях (включая орбиту Земли r/R 215) солнечный ветер состоит в основном из убегающих частиц (~64% при r=1а.е.) и, следовательно, может моделироваться бесстолкновительным потоком плазмы.
Рис. Зависимости скорости солнечного ветра от безразмерного гелио центрического расстояния u2 (r ) (штриховая линия) и u 2 (r ) (15) (верхняя сплошная линия) при T0=1.2 106K и r0=1.5R, в сравнении с зависимостью u (r ) (нижняя сплош ная линия) и с данными наблюде ний, опубликованными Якубовым.
Область, соответствующая нера венству (13) u (r ) u (r ) u 2 (r ), за крашена серым цветом. Пунктирная и штрих-пунктирная линии соот ветствуют скоростям солнечного ветра, рассчитанным по двухжид костной гидродинамической модели Хартля-Барнеса и по одночастич ной кинетической модели Лемера.
В главе 3 предлагается другая модель стационарного сферически симметричного потока электронно-протонной плазмы, основанная на бесстолкновительном двухчастичном кинетическом уравнении. В ней вводится дополнительное замыкающее допущение о том, что решение задачи зависит только от суммарной кинетической энергии динамической пары электрон-протон и, таким образом, не содержит потенциал поляризации плазмы. Показывается, что эта модель дает верхнюю оценку (мажоранту) для радиальной зависимости числовой плотности и нижнюю оценку для радиальной зависимости скорости потока частиц солнечного ветра, полученных в рамках «квазинейтральной» модели в предыдущей главе.
Теоретические результаты сравниваются с данными наблюдений по медленному солнечному ветру. Исследуется возможность применения модификации рассматриваемой модели к быстрому солнечному ветру.
Предположение о квазинейтральности плазмы означает, как видно из главы 2, что движение электронов и протонов описывается двухчастичной функцией распределения, зависящей от суммы энергий электронов и протонов f = f (r, r, ), если начальная функция распределения f0 и граничные условия для нее обладают таким же свойством. Равновесное начальное состояние плазмы (максвелловское распределение частиц по скоростям (2)) соответствует этому требованию. Решение кинетического уравнения определено в фазовом пространстве скоростей в области, ограниченной характеристиками этого уравнения. Предположим, что в отличие от модели, рассмотренной в главе 2, первые интегралы вдоль граничных характеристик также зависят только от суммы энергий электронов и протонов. Тогда двухчастичное кинетическое уравнение с равновесными граничным условием при r = r0 запишется в следующем виде:
f f f + ( re + rp ) = 0, (9) r r (me m p )3 exp 0, r = r0 : f 0 = f (r0, u rp 0, u p 0 ) = 4 N 0 (10) kT (2kT0 )3 2 2 2 r = (meure + m purp ) 2,. = (meu e + m pu p ) 2, = e + p = r +, где re + rp = (me + m p ) d dr. Построенная таким образом модель замкнута, так как не зависит от потенциала поляризации плазмы.
Сравнивая задачи (1), (2) и (9), (10), заметим, что в рамках первой из них процесс описывается статистическими характеристиками электронов и протонов в отдельности в соответствующем четырехмерном фазовом пространстве скоростей - uкe, u e, urp, u p. Зависимость модели (1), (2) от поляризационного потенциала обусловлена граничными характеристиками (4) и условием для потенциального барьера, разделяющим частицы атмосферы и убегающие частицы: Модель (9), (10) не включает в себя характеристики кулоновского взаимодействия частиц и описывает статистику нейтральных псевдочастиц, каждая из которых, не имеет постоянного состава, а является динамическим образованием, состоящим из пары противоположно заряженных частиц. С физической точки зрения такое преобразование (re, e, rp, p) (r, ) можно интерпретировать как перенос условий квазинейтральности с уровня электростатического взаимодействия коллективов электронов и протонов (через среднее поляризационное поле) на уровень их индивидуального, попарного взаимодействия, неявно предполагает введение дополнительных источников энергии.
Решение задачи (9), (10) имеет также вид (3), однако область ее определения:
) ( D 2 = {0 u u 2, u r 2 u r }, (11) расширяется в четырехмерном фазовом пространстве скоростей U = {uкe, u e, urp, u p } по сравнению с областью определения D (4) первой модели. Иначе говоря, область определения решения «квазинейтральной» модели D является подобластью области определения решения «нейтральной» модели D2 для любых реализаций поляризационного потенциала (r ), допустимых в рамках рассматриваемой стационарной задачи (6):
D D2 (r ) D, где D, D2 U. (12) 2 r m(0 ) ( ) ( m ), u 2 =, ur 2 = Здесь u r = ± r u =,, r 2 m m r = r r0, D - область (6) значений потенциала. Поскольку по областям D и D2 проводится интегрирование положительно определенной функции f при вычислении плотности числа частиц (5) в рамках соответственно первой N и второй N2 моделей, то, следовательно, N2(r) является верхней оценкой (мажорантой) для N(r) при любых допустимых реализациях потенциала :
Проведенный в главе анализ показывается с учетом (7), что имеют место следующие неравенства, связывающие аналитические зависимости, полученные в рамках двух моделей (r ) D при r r0 :
N (r ) N (r ) N 2 (r ), N atm (r ) N atm (r ) N 2 atm (r ), N ( r ) N ( r ) N 2 ( r ), Nu (r ) Nu (r ) Nu 2 (r ), u (r ) u (r ) u 2 (r ). (13) Здесь индексом 2 помечены средние параметры потока плазмы, полученные в рамках второй («нейтральной») модели, а индексом atm - числовые плотности частиц атмосферы (то есть частиц, не обладающих достаточной кинетической энергией для преодоления потенциального барьера).
Неравенство (13) для скоростей солнечного ветра иллюстрирует рис.2 в виде области, закрашенной серым цветом и заключенной между графиками функций u (r ) (нижняя сплошная линия) и u 2 (r ) (штриховая линия).
Отметим, что зависимость Nu 2 (r ) не удовлетворяет условию 4r 2 Nu 2 (r ) = const, которое должно выполняться для расхода частиц через сферическое сечение в случае стационарного течения. Это объясняется тем, что данная модель не обеспечена какими-либо физическими механизмами реализации упомянутой выше «попарной» квазинейтральности. Однако указанное условие выполняется для предельной зависимости Nu 2 (r ) :
( 0 2 ) + 2 N Nu 2 (r ) 0 u ep 0 exp 0, r 1. (14) 2 2 r На рис.2 верхняя сплошная линия соответствует зависимости u 2 (r ), вычисленной по формуле (14):
2 + (2 0 ) 2 Nu 2 (r ) u 2 (r ) = = u ep 0, (15) 2(r ) r N 2 (r ) 0 0 1 1 + (1 0 ) (r ) = 2 exp( )1 exp 2 1 r2 2r r 1 r Предельные соотношения (14), (15), как показывается в работе, удовлетворяют следующим неравенствам:
Nu (r ) Nu (r ) Nu 2 (r ) Nu 2 (r ), u (r ) u (r ) u 2 (r ) u 2 (r ), (16) последнее из которых иллюстрирует рис.2. Все зависимости на рис. построены при температуре T0=1.2 106K, при которой получено согласование теоретических и наблюдательных данных в предыдущей главе (r0=1.5R).
Однако мажорантные зависимости, полученные в рамках рассматриваемой «нейтральной» модели, также согласуются с данными наблюдений по медленному солнечному ветру как по числовой плотности, так и по скорости при температуре, несколько меньшей приведенного выше значения, а, именно, при T0=1.1510 К, r0=1.5R максимальное отклонение теоретической зависимости N 2 (r ) от эмпирической зависимости Рубцова N emp (r ) не превышает 13%, а теоретические значения скорости u2 (r ) лежат в пределах разброса данных измерений, опубликованных в монографии Якубова.
Сравнение теоретических зависимостей с наблюдательными данными проведено выше для граничного гелиоцентрического расстояния r0=1.5R.
Согласно имеющимся наблюдениям можно ожидать, что приближение равновесной плазмы (максвелловское распределение частиц по скоростям f0), заложенное в рассматриваемые модели, является удовлетворительным на этом расстоянии. Однако частота столкновений в этой области солнечной короны настолько высока, что предположение о бесстолкновительном характере течения на столь малых гелиоцентрических расстояниях представляется необоснованным. В то же время, зависимости для плотности и скорости солнечного ветра, полученные выше для бесстолкновительного приближения, согласуются количественно с данными наблюдений. Для потока частиц Nu это можно объяснить тем, что он формируется убегающими частицами, относящимися к высокоскоростной части функции распределения, а, следовательно, согласно проведенным в главе 2 оценкам относительные энергетические потери таких частиц при столкновениях на расстояниях полутора солнечных радиусов составляют порядка 20-25% и быстро убывают с расстоянием. Что касается числовой плотности частиц, то влияние столкновений на эту среднестатистическую величину обусловлено изменением числа частиц в области интегрирования D при учете этого фактора. Такое изменение числа частиц происходит за счет их перераспределения через фазовую поверхность, определяющую потенциальный барьер, и оно относительно мало. Однако влияние столкновений на вторые статистические моменты, например, температуру, существенно, и им нельзя пренебречь в отличие от рассмотренных выше первых моментов (плотности и скорости).
Учитывая, что медленный солнечный ветер можно рассматривать как бесстолкновительный поток, начиная с расстояний примерно в 6 радиусов Солнца (согласно оценкам Лемера и Шерера), и, с другой стороны, максвелловское распределение частиц по скоростям приемлемо в области, где столкновения играют значительную роль (и обеспечивают равновесное состояние), сдвинем задание граничных условий до r0=5R и сравним теоретические зависимости с эмпирическими, полученными В.Кенляйном в результате обработки результатов различных наблюдений медленного солнечного ветра. При r0=5R температура и плотность принимают в соответствии с эмпирическими зависимостями следующие значения:
T08.9105K и N(r0)4.51010м-3. При заданных таким образом входных параметрах различие между теоретическими зависимостями N2(r), u2(r) и соответствующими эмпирическими зависимостями при rr0 составляет менее 30% для плотности и менее 16% для скорости потока. Интересно отметить, что в интервале 30Rr370R выполняются неравенства u*(r)uemp(r)u2(r) и Nemp(r)N*(r)N2(r), совпадающие с неравенствами (13) и (16), если эмпирические значения плотности Nemp(r) и скорости uemp(r) отождествить с плотностью N(r) и скоростью u(r), вычисленными по первой модели при значении электростатического потенциала из области его допустимых значений.
В работе рассмотрено приложение «нейтральной» модели к так называемому быстрому солнечному ветру, который, представляет собой поток плазмы, истекающей из корональных дыр на высоких гелиоширотах.
Он характеризуется большим градиентом скорости потока и меньшей плотностью на малых гелиоцентрических расстояниях по сравнению с медленным солнечным ветром, который наблюдается вблизи экваториальной плоскости Солнца. Теоретическая зависимость N2(r) согласуется с данными наблюдений для северных и южных корональных дыр при задании наблюдаемых значений начальной температуры T0 и электронной плотности Ne(r0)=0.5N2(r0) в качестве входных параметров модели. (r0=1.6R, T0=1.10106K, Ne=6.911011 м-3). В то же время зависимость u2(r) (15) возрастает значительно медленнее, чем показывают данные измерений скорости солнечного ветра, однако на больших гелиоцентрических расстояний она достигает наблюдаемые значения (что свидетельствует о ведущей роли убегающих высокоэнергетических электронов и протонов в формировании быстрого солнечного ветра вне короны Солнца).
Для адаптации модели к описанию быстрого солнечного ветра в ней учитываются с помощью феноменологических формул процессы, приводящие к дополнительному подводу энергии на малых гелиоцентрических расстояниях, и (или) влияние геометрической трансформации потока. Проведенный анализ показывает, что согласование теоретических и наблюдательных данных достигается как в короне, так и на уровне орбиты Земли при включении в «нейтральную» модель дополнительного источника в уравнение энергии. При этом согласование достигается при большей мощности этого источника, вне того диапазона значений, который обычно рассматривается в МГД моделях (однако в этих моделях учитываются одновременно и другие процессы, служащие источниками энергии, как, например, альфвеновские волны). Введение в «нейтральную» модель фактора расширения потока, вытекающего из корональной дыры (то есть соответствующая модификация уравнения сохранения моментов количества движения) ухудшает согласование теоретических результатов по плотности с данными измерений на уровне орбиты Земли. Внесение изменений в уравнение сохранения момента количества движения в предположении, что течение плазмы, испускаемой Солнцем, формируется в виде сопла, также позволяет получить согласование теоретических и доступных наблюдательных данных для скорости и плотности быстрого солнечного ветра в короне и на уровне земной орбиты.
Следует отметить, что размеры (параметры) упомянутого выше «сопла», позволяющего согласовать теоретические и наблюдательные результаты в отсутствии дополнительного подвода энергии, не соответствуют имеющимся наблюдениям по структуре магнитного поля и потока плазмы. Можно предположить, что адекватная модель быстрого солнечного ветра может быть построена на основе сочетания двух факторов: 1) источников энергии известной физической природы (например, альвеновские волны, обычно включаемые в МГД модели) и 2) ускорения потока в магнитном «сопле», которое может быть сформировано, например, в течении с пересоединением силовых линий магнитного поля. В этом случае, вполне возможно не потребуется введение дополнительных источников энергии неизвестной природы, как это практикуется в современных МГД моделях быстрого солнечного ветра.
Таким образом, для применения двухчастичной кинетической модели к описанию быстрого солнечного ветра необходимо учесть процессы, приводящие к дополнительному подводу энергии на малых гелиоцентрических расстояниях, или геометрическую трансформациию потока в корональных дырах. Из анализа этих двух факторов на основе аналитических зависимостей для плотности и скорости потока солнечной плазмы и соответствующих феноменологических формул, широко используемых в МГД моделях, следует, что они оказывают близкое по характеру влияние на параметры солнечного ветра. Причем совместный учет этих двух факторов в указанной форме приводит к рассогласованию теоретических и наблюдательных данных при увеличении степени суперрадиального расширения трубок тока.
В главе 4 предложен многочастичный статистический подход к описанию потока газа и плазмы, основанный на теореме Лиувилля, и метод вычисления статистических моментов в рамках этого подхода при предположении о неразличимости частиц в объеме, масштаб которого обусловлен разрешением измерений и(или) принятыми модельными допущениями. Указанный подход протестирован на классических задачах.
Рассмотрено его приложение к описанию стационарного сферически симметричного бесстолкновиетльного потока однокомпонентного газа и квазинейтральной полностью ионизованной многокомпонентной плазмы.
Распределение частиц по скоростям в сечении истечения принято равновесным. На основе полученной многочастичной функции распределения частиц газа и плазмы по скоростям, находятся соответствующие средние параметры (плотность и скорость) потока газа и солнечной плазмы. Полученные радиальные зависимости для плотности и скорости в случае двухкомпонентной (электрон-протонной) плазмы согласуются с результатами двухчастичной модели, предложенной в главе 2, и с данными наблюдений.
Существующие кинетические модели, описывающие солнечный ветер в широком диапазоне изменения гелиоцентрического расстояния, строятся в приближении стационарного сферически симметричного потока квазинейтральной бестоковой полностью ионизованной плазмы и основываются на уравнениях для одно- и двухчастичных функций распределения.
Рассмотрим статистическое описание потока частиц на основе совместной многочастичной функции распределения вероятности f N (r1, p1,..., rN, p N ). Здесь ri, pi - пространственные координаты и импульсы выделенных N частиц (N принимается меньшим полного числа частиц системы N). Заметим, что при классическом подходе принятое определение статистически средней локальной числовой плотности системы N частиц (размеры частиц принимаются нулевыми) базируется на определении микроскопической (истинной) локальной плотности как суммы вероятностей каждой частицы иметь точечные отличимые друг от друга координаты N n(r, t ) = (ri (t ) r ). Это предположение принимается независимо от того, i = каковы масштабы различимости положения отдельной частицы в физическом пространстве, связанные, например, с разрешением проводимых измерений, или с воздействием частиц на какой-либо физический объект, размер которого превышает среднее расстояние между частицами. В уравнении Леонтовича для функции распределения, сглаженной по физически малому объему, осреднение по этому объему проводится также на основе суммирования вероятностей. Родственный подход предлагает В.Я.Рудяк, проводя осреднение (сглаживание) функции распределения по эффективному радиусу межмолекулярного взаимодействия, который неразличим с точки зрения статистического описания разреженного газа как системы точечных частиц.. При этом все частицы по-прежнему принимаются различимыми.
Более обоснованным представляется рассматривать совместную вероятность нахождения частиц в объеме, соответствующем масштабу разрешения измерений и/или масштабу взаимодействия потока частиц с макрообъектом, и на основе совместной функции распределения вероятностей рассчитывать среднестатистические локальные характеристики системы частиц (к указанным масштабам «неразличимости» можно добавить и масштаб «неразличимости» принятого приближения, модели).
Итак, используем N-частичную функцию распределения f N (r1, u1,..., rN, u N, t ) для расчета числовой плотности частиц потока, исходя из предложенного выше подхода об определении совместной вероятности нахождения частиц в объеме (r ), внутри которого координаты частиц неразличимы (в рамках условий измерений и/или принятых приближений):
(r )( )( ) F (X ) ~ N ~ ~ ~ ~ ~ f N (X N, t ) = ri p ji pi F X N, t dX N = P, ( N, N ), t dX N N ji N 1 ji N i =1 N j1... j N D DN (r )( ) ~ N ~ N (r1,..., rN ) = ri p ji pi где - (сингулярная) динамическая ji 1 j i N i = j1... j N N-частичная фазовая плотность, X N = (r1, p1,..., rN, p N ), функция, или ~ ~~~ ~ ~~ ~ ~ ~ dX N = dr1...drN dp1...dp N, dX N N = drN +1...drN dp N +1...dp N, - число P, ( N, N ) возможных перестановок из полного числа частиц всей системы N по N местам в объеме (r ), дающие эквивалентные макросостояния. Функция распределения вероятностей F для системы из N частиц в ограниченном объеме, находящейся в динамическом равновесии с окружающей средой, удовлетворяет уравнению Лиувилля и нормирована по ее области определения {ri, pi D,1 i N }. Интегрирование в (1) проводится по всей области определения функции распределения F за исключением объема конфигурационного пространства и области фазового пространства для N частиц из этого объема. Ниже принимается, что рассматриваемые системы частиц стационарны, и на поверхности задания граничных условий они имеют максвелловское распределение по скоростям). Пространственными корреляциями пренебрегается: F (X N, t ) = f1 (r1, p1 )... f1 (rN, p N ). Например, для ~ тестовой задачи о флуктуациях плотности в ограниченном объеме, занятом однородным идеальным газом, приходим к биноминальному распределению вероятности актуальных значений числа частиц в выделенном объеме (которое, как известно, в пределе при больших значениях числа частиц переходит в распределение Пуассона, а при относительно малых флуктуациях в распределение Максвелла). Кроме того, для задачи об идеальном газе в гравитационном поле получаем также биноминальное распределение со средним значением числовой плотности, равным классическому распределению Больцмана.
Для стационарного сферически симметричного бесстолкновительного потока частиц одного сорта в гравитационном поле получено биноминальное распределение вероятности нахождения N частиц в пробном объеме.
Средняя числовая плотность для этого распределения совпадает с выражением, выведенным на основе одночастичной функции распределения.Многочастичный подход применен к описанию солнечного ветра в приближении бесстолкновительного потока квзинейтральной водородной плазмы. Из предположения о неразличимости частиц одного сорта по их свойствам, следует, что эквивалентные макросостояния реализуются при перестановках частиц одного сорта независимо от частиц других сортов. Исходя из этого, вероятность нахождения N частиц плазмы в объеме определяется как произведение вероятностей для каждого сорта частиц и выражается в виде полиноминального распределения. Выражение для средней плотности двухкомпонентной плазмы совпадает с аналитической зависимостью, полученной на основе двухчастичной кинетической модели, которая удовлетворительно согласуется с данными наблюдений.
Совпадение пространственных зависимостей для плотности и скорости солнечного ветра, полученных на основе многочастичной и двухчастичной статистических моделей, объясняется тем, что для стационарного сферически симметричного бесстолкновительного потока плазмы движение всех частиц статистически независимо и одинаково, а, следовательно, полночастичная функция распределения представляется как произведение одинаковых в рамках каждой фракции одночастичных функций распределения.
Аналогично для такого же потока k-компонентной плазмы результаты многочастичного подхода совпадают с результатами k-частичной модели. В случае стационарного сферически симметричного столкновительного течения однокомпонентного идеального газа, рассмотренного в настоящей главе, совпадают результаты много- и одночастичных моделей, так как движение частиц в этом случае также статистически независимо. Однако ситуация меняется, когда, например, задаются различные граничные условия для частиц одной фракции (в частности, если они поступают в рассматриваемый объем через сферу радиуса r0 и с бесконечности, тогда фактически их можно рассматривать принадлежащими двум разным фракциям, и многочастичный подход сведется к k+1-частичной модели).
Другим примером может служить k-компонентная система статистически зависимых частиц (взаимодействие между которыми нельзя ограничить мгновенными актами столкновений). Обратим внимание на то, что в случае разного числа частиц во фракциях, возникает проблема определения средних параметров потока в целом из произведения их значений для каждой фракции, которое получается на основе мультипликативного (многочастичного) подхода.
Заметим, что полученные теоретические результаты (на базе многочастичной и двухчастичной статистических моделей) для числовой плотности солнечного ветра демонстрируют лучшее количественное согласование с данными наблюдений по сравнению с одночастичными кинетическими моделями, также задающими распределение Максвелла на поверхности истечения солнечного ветра. Последнее может свидетельствовать в пользу многочастичного подхода, основанного на гипотезе неразличимости частиц в пробном объеме, размер которого определяется масштабом разрешения проводимых измерений.
ОБЩИЕ ВЫВОДЫ ДИССЕРТАЦИИ Представленная диссертационная работа отражает результаты исследования автора в области решения следующих задач: 1) построение двухчастичных кинетических моделей стационарного потока двухкомпонентной солнечной плазмы;
2) разработка многочастичного подхода к описанию потока газа и плазмы;
3) получение на основе вышеуказанных моделей аналитических точных и приближенных формул для первых статистических моментов потока частиц при ряде упрощающих предположений;
4) исследование поведения характеристик потока плазмы при изменении условий ее течения на основе выведенных зависимостей.
Исследования проводились в Томском государственном университете в течение 2000-2005 гг.
Основные результаты настоящей работы сводятся к следующему:
1. Представлена новая двухчастичная кинетическая модель стационарного потока бесстолкновительной квазинейтральной двухкомпонентной (водородной) плазмы, не учитывающая влияние магнитного поля и зависящая от потенциала поляризации плазмы как неизвестной функции, нахождение которой требует введения дополнительных гипотез.
Указанная модель применена к солнечному ветру, получены точные аналитические решения соответствующего кинетического уравнения для функции распределения вероятностей и на его основе выведены аналитические пространственные зависимости скорости и плотности потока солнечной плазмы, которые количественно согласуются с данными наблюдений по медленному солнечному ветру в приближении Паннекока Росселанда для поляризационного потенциала.
Исследовано поведение полученных пространственных зависимостей от потенциала поляризации плазмы и показано, что приближение Паннекока Росселанда для этого потенциала, полученное в предположении термодинамически и гидростатически равновесной плазмы, доставляет максимум для плотности и удельного потока частиц и минимум для его скорости.
2. Представлена замкнутая двухчастичная кинетическая модель стационарного потока бесстолкновительной двухкомпонентной (водородной) плазмы как статистической системы нейтральных динамических электрон протонных пар, не зависящая от поляризационного потенциала плазмы (без учета влияния магнитного поля).
Предлагаемая модель применена к солнечному ветру, получены точные аналитические решения соответствующего кинетического уравнения для функции распределения вероятностей и на его основе выведены аналитические пространственные зависимости скорости и плотности потока солнечной плазмы, которые не содержат параметров согласования и полностью определяются температурой и плотностью втекающего через границу потока плазмы.
Показано, что зависимости для плотности и скорости потока частиц, полученные на основе данной модели служат верхними оценками для соответствующих зависимостей, выведенных в рамках модели из п.1, и количественно согласуются с данными наблюдений по медленному солнечному ветру при меньших значениях температуры плазмы, втекающей в рассматриваемую область (у основания солнечной короны).
3. На основе модели, указанной в п.2, предложена полуэмпирическая модель течения солнечной плазмы в расширяющихся трубках тока, исходящих из корональных дыр, с подводом энергии на малых гелиоцентрических расстояниях от источника неизвестной природы.
Получены в рамках указанной модели точные аналитические решения кинетического уравнения и аналитические зависимости для плотности и скорости потока плазмы. В рамках данной модели проведен анализ влияния указанных факторов на характер течения плазмы и согласование теоретических результатов с данными наблюдений.
4. Разработан многочастичный статистический подход к моделированию стационарного неоднородного потока газа и плазмы, который основан на теореме Лиувилля, и методология вывода в рамках этого подхода статистических моментов потока газа и плазмы в предположении о неразличимости частиц в объеме, масштаб которого обусловлен разрешением измерений или принятыми модельными допущениями. Данный подход протестирован на решении классических задач статистической физики для столкновительного идеального газа (о флуктуациях в однородном идеальном газе и об атмосфере в гидростатическом равновесии).
5. В рамках модели, указанной в п.4, получены аналитические выражения для числовой плотности стационарного сферически симметричного бесстолкновительного потока однокомпонентного газа и для флуктуационной функции распределения его частиц по скоростям. Выведены аналитические выражения для числовой плотности и скорости стационарного сферически симметричного потока многокомпонентной плазмы, в том числе двухкомпонентного солнечного ветра (в пренебрежении столкновениями).
Показано, что для стационарного сферически симметричного потока k компонентной квазинейтральной плазмы, истекающей со сферы заданного радиуса, полночастичный статистический подход, основанный на уравнении Лиувилля, сводится к k-частичной кинетической модели (при статистической независимости движения частиц). Указаны частные случаи, когда понижение размерности функции распределения, описывающей течение плазмы (газа), возможно лишь до значений, больших числа ее компонент.
6. Предлагаемые модели потока солнечной плазмы количественно согласуются с данными наблюдений по числовой плотности и скорости медленного солнечного ветра в широком диапазоне изменений гелиоцентрического расстояния. При этом для достижения наблюдаемых значений указанных параметров потока не требуется введения дополнительных источников энергии в отличие от гидродинамических и одночастичных кинетических моделей.
Показано в рамках рассматриваемых моделей для случая медленного солнечного ветра, что ускорение потока солнечной плазмы в гравитационном поле Солнца обусловлено уменьшением доли частиц атмосферы (кинетическая энергия которых меньше гравитационного потенциального барьера) по сравнению с долей убегающих частиц (кинетической энергии которых достаточно для преодоления гравитационного потенциального барьера). При этом скорости всех частиц убывают с удалением от Солнца.
7. Полученные теоретические результаты (на базе многочастичной и двухчастичной статистических моделей) для числовой плотности и скорости солнечного ветра демонстрируют лучшее количественное согласование с данными наблюдений по сравнению с одночастичными кинетическими моделями, также задающими распределение Максвелла на поверхности истечения солнечного ветра (у основания короны). Последнее может свидетельствовать в пользу многочастичного подхода, основанного на гипотезе неразличимости частиц в пробном объеме, размер которого определяется масштабом разрешения проводимых измерений Результаты диссертационной работы могут быть использованы при решении следующих задач: 1) прогнозирование параметров солнечного ветра у орбиты Земли по данным наблюдений в солнечной короне или для восстановления параметров солнечного ветра на заданном гелиоцентрическом расстоянии по измерениям, проведенным вблизи Земли;
2) построение кинетических моделей потоков газа и плазмы;
3) вычисление коэффициентов переноса для таких систем при описании их в рамках приближений сплошной среды.
СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ
1. Васенин И.М., Минькова Н.Р. Точное решение бесстолкновительных кинетических уравнений в применении к солнечному ветру.
Математические модели и методы их исследования: // Труды международной конференции. - Красноярск: 2001. т.1, с. 137-141.
2. Васенин И.М., Минькова Н.Р. Оценка параметров стационарного солнечного ветра путем решения двухчастичного кинетического уравнения. // Материалы Пятой Сибирской школы-семинара «Математические проблемы механики сплошных сред» - Новосибирск, 2001. C.34.
3. Васенин И.М., Минькова Н.Р. Оценка параметров стационарного солнечного ветра на основе кинетического уравнения для нейтральной плазмы. // Материалы IV Всероссийского семинара «Моделирование неравновесных систем-2001». - Красноярск, 2001. C.48-49.
4. Васенин И.М., Минькова Н.Р. Кинетическая модель солнечного ветра. // Конференция «Проблемы и перспективы технологий атомной промышленности». Сборник тезисов - Томск, 2001. C.27.
5. Васенин И.М., Минькова Н.Р. Кинетическая модель солнечного ветра. // Известия Томского политехнического университета. – 2002. - Т.305, вып.3. - С.331-335.
6. Васенин И.М., Минькова Н.Р. Кинетическая модель солнечного ветра с учетом влияния магнитного поля. // Материалы V Всероссийского семинара «Моделирование неравновесных систем-2002». - Красноярск, 2002. C.26-27.
7. Васенин И.М., Минькова Н.Р. Некоторые проблемы оценки скорости солнечного ветра. // Материалы IV Международной научно-практической конференции «Фундаментальные и прикладные проблемы приборостроения, информатики, экономики и права». Сочи, 2002. – M., 2002. C.70-74.
8. Vasenin, Y.M., Minkova, N.R., Shamin A.V. A kinetic model of solar wind.
2002. // 11th International Congress on Plasma.Physics. ICPP-2002. Abstract book. P.153.
9. Vasenin, Y.M., Minkova, N.R. Two-particle quasi-neutral kinetic model of collisionless solar wind. // Strongly Coupled Coulomb Systems. Santa Fe, New Mexico, USA, 2-6 September 2002.
Abstract
book, P.52.
10. Васенин И.М., Минькова Н.Р. Оценка параметров стационарного солнечного ветра путем решения двухчастичного кинетического уравнения. // Динамика сплошной среды: Сб. науч. тр. / РАН. Сиб. отд ние. Ин-т гидродинамики, Новосибирск, 2002. Вып. 120. С. 84-88.
11. Vasenin, Y.M., Minkova, N.R., Shamin A.V. A kinetic model of solar wind. // AIP Conference Proceedings, 2003, v. 669, series PLASMA PHYSICS: 11th International Congress on Plasma Physics: ICPP2002. P.516-519.
12. Vasenin, Y.M., Minkova, N.R. Two-particle quasi-neutral kinetic model of collisionless solar wind. // Journal of Physics A. Mathematical and General. – 2003. V.36, Issue 22. P.6215-6220.
13. Минькова Н.Р. К моделированию солнечного ветра двухчастичными кинетическими уравнениями. // Конференция стран СНГ и Прибалтики "Актуальные проблемы физики солнечной и звездной активности", 2- июня 2003 года. Тезисы докладов - Нижний Новгород, 2003. C. 14. Vasenin, Y.M., Minkova, N.R. Two-particle kinetic model of the fast solar wind. // IUGG General Assembly, Sapporo, Japan, 30 June-11 July 2003.
Abstracts Book. – Sapporo, 2003. P. А.328.
15. Vasenin, Y.M., Minkova, N.R. Shamin A.V. Two-particle kinetic model in application to the fast solar wind. // Материалы VI Всероссийского семинара “Моделирование неравновесных систем-2003”.- Красноярск, 2003. C.53- 16. Минькова Н.Р. К моделированию солнечного ветра двухчастичными кинетическими уравнениями. // Актуальные проблемы физики солнечной и звездной активности, / Конференция стран СНГ и Прибалтики (Нижний Новгород, 2-7 июня 2003 г.), Сборник докладов в двух томах. Том 1. Нижний Новгород: ИПФ РАН, 2003, с 209-212.
17. Васенин И.М., Минькова Н.Р. Многочастичная квазинейтральная кинетическая модель полностью ионизованной плазмы в приложении к медленному солнечному ветру. // Тезисы докладов XXXI Звенигородской конференции по физике плазмы и управляемому термоядерному синтезу.
г.Звенигород, Московской обл., 16 - 20 февраля 2004 года.– M.,2004. С. 18. Васенин И.М., Минькова Н.Р. Многочастичная квази-нейтральная кинетическая модель полностью ионизованной плазмы.// Сборник тезисов докладов Всероссийской Астрономической Конференции ВАК- «Горизонты Вселенной» (Москва, 2-10 июня 2004 года). – M., 2004.С. 19. Minkova N.R. On Multiparticle Kinetic Models of the Solar Wind. // 10th International conference and School on Plasma Physics and Controlled Fusion.
Alushta (Cream), Ukraine, 2004. Book of Abstracts. P. 20. Васенин И.М., Минькова Н.Р., Шамин А.В. Двухчастичная кинетическая модель в применении к быстрому солнечному ветру. // Журнал «Проблемы эволюции открытых систем». – 2004. - Т.1, вып. 6. - С.73-78.
21. Минькова Н.Р. Многочастичный статистический подход к описанию потока солнечной плазмы.// Материалы VII Всероссийского семинара “Моделирование неравновесных систем-2004”. - Красноярск, 2004.C.107 108.
22. Минькова Н.Р. Многочастичный статистический подход к моделированию солнечной бесстолкновительной плазмы. //. Известия вузов. Физика. 2004 - Т.47 - №10. Приложение. (Прикладные проблемы сплошных сред.
Тематический выпуск). - С.73-80.
23. Минькова Н.Р. Многочастичная статистическая модель солнечного ветра.
// Тезисы докладов XXXII Звенигородской конференции по физике плазмы и управляемому термоядерному синтезу. г.Звенигород, Московской обл., февраль, 2005. – M., 2005. C.256.
24. Minkova N.R. Statistical Model of the Solar Wind based on the Liouville equation. // Solar Wind 11 - SOHO 16 «Connecting Sun and Heliosphere».
Whistler, Canada, 12-17 June 2005. Abstracts. P.206.
25. Minkova N.R. Multiparticle statistical approach to solar wind modeling. // International Conference on Strongly Coupled Coulomb Systems (SCCS-2005).
Moscow, June 20-25, 2005. Abstract book. P.51.