CC BY
55
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
УДК 519.6
С. А. Ишанов, Л. В. Зинин, Е. В. Зубков, С. В. Мациевский
КОЛЕБАТЕЛЬНО-ВОЗБУЖДЕННЫЙ МОЛЕКУЛЯРНЫЙ АЗОТ В ИОНОСФЕРНОЙ ПЛАЗМЕ
Представлена математическая модель для исследования нестационарных процессов, протекающих в ионосфере и плазмосфере Земли. Она позволяет исследовать динамику взаимодействия заряженной и нейтральной компонент ионосферной плазмы, а также выявить вклад различных процессов в формирование ионосферы. Результаты вычислительного эксперимента указывают на сильную зависимость электронной концентрации от высотного распределения азота. Учет колебательного возбуждения молекулярного азота приводит к заметным вариациям электронной концентрации и температуры на высотах слоя F2 ионосферы и улучшает согласие результатов модельных расчетов вариаций электронной концентрации с экспериментальными данными. Показано, что вклад малых компонент в аэрономические процессы существенен.
The mathematical model for research of the non-stationary processes proceeding in ionosphere and plasmasphere of Earth is presented. This model allows to investigate dynamics of interaction charged and neutral components of ionospheric plasma, and also to reveal a contribution of various processes to formation of ionosphere. Results of computing experiment indicate strong dependence of electronic concentration on high-rise distribution of nitrogen. The accounting of oscillatory excitement of molecular nitrogen leads to noticeable variations of electronic concentration and temperature at ionosphere F2 layer heights and improves a consent of results of model calculations of variations of electronic concentration with experimental data. It is shown that the contribution of small component in aeronomiche processes is essential.
Ключевые слова: математическая модель, нестационарный процесс, ионосфера, плазмосфера, вычислительный эксперимент, электронная концентрация, азот, молекулярный азот, аэродинамика, экспериментальные данные.
Key words: mathematical model, non-stationary process, ionosphere, plasma-sphere, computing experiment, electronic concentration, nitrogen, molecular nitrogen, aerodynamics, experimental data.
Введение
Непрерывное возрастание техногенного (химического, электромагнитного) воздействия на ионосферу Земли требует глубокого понимания физико-химических процессов в околоземной плазме, умения прогнозировать возможное поведение этой многокомпонентной среды в различных геофизических условиях.
7
© Ишанов С. А., Зинин Л. В., Зубков Е. В., Мациевский С. В., 2014
Вестник Балтийского федерального университета им. И. Канта. 2014. Вып. 10. С. 7—15.
Возмущения естественного и антропогенного происхождения в атмосфере обусловливают образование возбужденных молекул на различных высотах.
Благодаря запасу внутренней энергии и большим сечениям взаимодействия эти компоненты играют заметную роль в кинетике микропроцессов ионосферной плазмы. Влияние колебательно-возбужденного азота N2^ на динамику и структуру ионосферы исследовалось в работах [1; 2; 4]. Колебательное возбуждение молекул азота влияет на скорость ионно-молекулярной реакции
O+ + N2 ^ N0+ + N (М)
которая вместе с реакцией
O+ + O2 ^ O+ + O
определяет рекомбинацию заряженных частиц в F2-области.
От распределения молекул азота по колебательным уровням зависит скорость охлаждения электронного газа, а, следовательно, и электронная температура [2].
Цель данной работы состоит в получении высотно-временного распределения концентраций колебательно-возбужденного азота N2^ (V = 1, 10) для среднеширотной ионосферы и в изучении его влияния на концентрацию электронов в F-области.
Корректный расчет концентраций возможен на основе совместного решения системы уравнений, описывающей самосогласованную динамику нейтральной, заряженной и колебательно возбужденной компоненты ионосферной плазмы.
1. Нестационарная математическая модель тепловой плазмы в магнитных силовых трубках
Математическая модель [4—7] построена для исследования нестационарных процессов, происходящих в ионосфере и плазмосфере Земли. Модель основана на квазигидродинамическом описании плазмы и позволяет вычислять концентрации, потоковые скорости и температуры основных заряженных составляющих. Модель описывает частично ионизированную плазму, состоящую из свободных электронов, ионов
0+, Н+, О^, N0+, N2, Не+ и нейтральных компонентов
0(1Б), 0(1Б), ^Б), N(^0), N 2^, 02(Х), N0, О, 02,
Все параметры среды рассчитываются вдоль силовой линии геомагнитного поля, взятого в дипольном приближении, что позволяет свести уравнения модели к системе нестационарных одномерных уравнений. Силовая линия магнитного поля рассматривается как замкнутая, коротирующая вместе с Землей (модель среднеширотной ионосферы и плазмосферы).
8
Математическая модель включает уравнения непрерывности, движения и теплового баланса заряженных частиц. Уравнения непрерывности используются в виде
дпа 1 дАпаиа _ _ 1 ...
-я-_ Ьа ' (1)
дt А дэ
где па — концентрация заряженных частиц сорта а (индекс а принимает значения 0+, Н+, ...); иа — переносная скорость частиц; Ьа — скорость гибели и Qа — скорость рождения частиц; t — время; э — пространственная координата, отсчитываемая вдоль силовой линии магнитного поля; А(э) — так называемая расходимость силовых линий. В правых частях уравнений учитывается ионизация солнечным ультрафиолетовым излучением, реакции перезарядки и рекомбинации.
Уравнения движения имеют вид
диа диа 1 дРа . .
-дг + иа~дТ +--7Т_ Га , (2)
дt дэ тапа дэ
здесь та — масса частицы сорта а; Ра — давление, равное кпаТа, где к — постоянная Больцмана; Та — температура; Ра — внешняя сила, включающая силу трения, силу тяжести и электрическую силу. Индекс а в уравнении (2) означает сорт положительно заряженного иона, как и в уравнении (1).
В силу предположения о квазинейтральности плазмы электронная концентрация и потоковая скорость электронов вычисляются по формулам
№е _Упа , ие .
е ^ а е №
Уравнения теплового баланса используются в модели в виде
3 кпа(дГ, + иа ^ + 1_д_ГАХа дГа] _ Эа, (3)
2 Ч дt а дэ ) А дэ А дэ У а дэ ) а
где Ха — коэффициент теплопроводности; Эа — разность скоростей нагрева и охлаждения частиц сорта а. В уравнении (3) индекс а принимает значения е , О+, Н+, ... Правая часть уравнения теплового баланса электронов включает в себя и скорость нелокального нагрева посредством сверхтепловых электронов.
Высотно-временное распределение малых компонентов в предположении, что в рассматриваемой области высот силовая линия совпадает с местной вертикалью, описывается уравнением непрерывности вида (1).
Полное описание коэффициентов уравнений модели ионосферы и плазмосферы, схемы фотохимических реакций, способа построения разностной сетки и других деталей математической модели приведено в [4 — 6].
В [4 — 7] описаны различные вычислительные алгоритмы решения уравнений модели. В частности, в [3] показаны границы применимости так называемого диффузионного приближения. В качестве базовых в модели используются неявные разностные схемы. Линеаризованные
9
10
уравнения непрерывности и движения для каждого сорта заряженных частиц решаются совместно, при этом пространственные производные аппроксимируются несимметричными конечными разностями. В модели предусмотрен двухступенчатый итерационный процесс, состоящий из итераций по нелинейности, выполняемых отдельно при численном решении уравнений теплового баланса и при решении систем уравнений непрерывности и движения для каждого сорта заряженных частиц, и внешнего итерационного цикла, включающего решение всех уравнений модели на текущем временном шаге.
Построенные таким образом вычислительные алгоритмы устойчивы при достаточно больших шагах интегрирования по времени, что позволяет применять их для исследования ионосферных процессов с характерными временами от нескольких минут до десятков суток.
Нестационарность модели и отсутствие необходимости задания верхних граничных условий позволили использовать ее для изучения различных процессов в ионосферной плазме [4 — 7].
2. Уравнение непрерывности для колебательно-возбужденного молекулярного азота
Уравнение, описывающее высотно-временное распределение колебательно-возбужденного молекулярного азота N2"' для каждого колебательного уровня (у), представим в виде
дпу д ( \ ^
—- + —(«V Vпг ) = -«упу + Оу , (4)
дt дг
где — концентрация N2"' на уровнях V = 1, ..., 9; упг — вертикальная составляющая скорости молекулярной диффузии; а — члены, учитывающие образования и потери N2"' в фотохимических реакциях (табл.).
Таблица химических реакций колебательно-возбужденного N
Реакция Константа скорости (см3 с-1)
0(1Б) + N. —— N2'' + 0(3Р) к1 = 5 10-11
N(4S) + N0 ———• N2'' + 0(3Р) к2 = 3 10-11
N2Vl) + о(3р' ———• о(3р'+ы2,1-1> к3 = 1,2 10-13ехр(-23/Т1/3)
N2Vl)+^2'2) ^ N2Vl-1) + ^2+1) к4 = 3 10-12
N2^ + е ^ + е к5 = ехр(Р1 + Р2 + Р3/ Те + Р4/ Те)
О + + N2Vl) —— N0 + + N % кб = Яу + ЬуТ« + с уТ«2 + йу-16
Примечание. Р1, Р2 Р3, Р4 Яу, Ьу, Су, йу — коэффициенты [5].
Скорость молекулярной диффузии определяется как
( 1 — 1 дТп 1 ^ --- +--- + —
дг Тп дг Н
кТ
где Н = —— — шкала высот №; т — масса молекулы N2; Тп — темпе-mg
ратура нейтрального газа.
Коэффициент диффузии О не зависит от номера колебательного уровня и рассчитывался по формуле:
5М. ( ^ I"5 см2. с-,,
^ М ^ 300)
где М — суммарная концентрация основных нейтральных составляющих в 02, N и О.
После подстановки (5) в (4) получим
дпг='дг п^+Qv, (6)
( 1 дТ 1Л здесь А = О; р = О п
--— + —
V Тп дг Н,
Система уравнений (6) дополняется нижними и верхними граничными условиями. Нижние граничные условия определяются из предположения фотохимического равновесия, а верхние — из условия равенства нулю потока частиц. При решении уравнений (6) использовался потоковый вариант метода прогонки.
Термодинамическая характеристика ионосферной плазмы — колебательная температура К/*1 рассчитывается по формуле
й
Т = —
1п
(V) Л '
1 +
I 8=0 N 2 I 8=1N2^
где 0 = 3390К.
Колебательная температура учитывается при определении скорости охлаждения электронного газа при возбуждении колебательных уровней молекул К2.
Фотохимические источники и стоки Qv и аv для колебательно-возбужденного азота рассчитываются согласно работам [2; 4; 5] с учетом таблицы.
В качестве основных источников N2'' приняты: экзотермическая реакция атомарного азота с окисью азота, гашение метастабильных атомов кислорода в состоянии ^ на молекулярном азоте, взаимодействие с атомарным кислородом 0(3Р), с тепловыми, быстрыми электронами.
Необходимо отметить, что на высотах порядка 200 км начинает играть роль молекулярная диффузия для N2
11
12
3. Результаты вычислительных экспериментов
По модели МКЬМБКЕ-ОО [8] вычислялись концентрации нейтральных компонентов и температура нейтральной атмосферы.
Высотно-временные распределения ионосферных параметров рассчитывались в дипольной системе координат, для средних широт (параметр Мак-Илвейна Ь = 3) для зимних условий, средней солнечной активности (Р107 = 150) и высокой солнечной активности (Р10.7 = 234). Другие входные данные брались согласно работам [5 — 7].
В работе возмущенные условия моделировались путем задания в начальных условиях завышенных значений температур и концентраций относительно фона.
Такие начальные условия могут соответствовать сильной солнечной вспышке или нагреву электронного газа мощным КВ-радиоизлу-чением.
Уравнения модели интегрировались методом конечных разностей с шагом по времени т = 0,01 с, затем с временным шагом т = 0,1 с в диффузионном приближении.
Нижние граничные условия — фотохимическое равновесие.
В качестве возмущенных начальных условий в полдень задавалось высотное распределение концентраций ионов, электронов, малых компонентов и электронной температуры с Те тах = 1,3 эв на высоте 300 км.
На рисунке 1 приводятся полученные на основе модельных расчетов высотные профили электронной Те и колебательной Ту температур.
Рис. 1. Высотные профили колебательной Т и электронной Те температур для дневных условий в невозмущенной среднеширотной ионосфере
В случае возмущения температура ТV может достигать 3000 К, вклад реакции (И1) в потери ионов 0+ в Р2-области заметно увеличивается, что приводит к существенному понижению электронной концентрации.
Важную роль в энергетическом балансе ионосферы играют электронно-возбужденные атомы кислорода 0(1Э) и колебательно-возбужденные молекулы азота. На рисунке 2 представлены высотно-временные профили ^(к). В расчетах электронная концентрация определялась из условия квазинейтральности. Через 30 мин от начала возмущения типа солнечной вспышки электронная концентрация релаксирует к фоновым значениям. На рисунке 3 показаны высотные распределения 0(1Э) в зависимости от момента времени.
105 2 4 10 4 10' Ке, см3
Рис. 2. Высотные профили электронной концентрации Nе в различные моменты времени
Рис. 3. Высотные профили концентрации 0(1Б) в различные моменты времени
В качестве основных источников для 0(1Э) взяты: диссоциативная рекомбинация 0+, соударения фотоэлектронов с атомным кислородом 0(3Р), реакция взаимодействия К(2Э) с 02. Переход атомарного кислорода 0(1Э) в основное состояние дает красную линию X = 6300 А и X = 6364 А.
Время установления концентрации 0(1Э) от начальных условий до фоновых значений составляет 30—40 минут.
На рисунке 4 показано высотно-временное распределение колебательно-возбужденного азота на первом колебательном уровне.
13
h, км
ivf
400
300
200
100
14
10б 2 4 6 8 106 2
4 6 8 106 2 4 6 8
N, см3
Рис. 4. Высотные профили распределения концентрации колебательно-возбужденного азота на первом уровне
В начальные моменты времени и до 30 мин сохраняются достаточ-
Представленная модель позволяет исследовать динамику взаимодействия заряженной и нейтральной компонент ионосферной плазмы, а также выявить вклад различных процессов в формирование ионосферы.
Результаты решения рассмотренной модельной задачи могут быть использованы для интерпретации существующих экспериментальных данных как по концентрации, так и по интенсивности свечения мета-стабильных компонент.
Результаты вычислительного эксперимента указывают на сильную зависимость электронной концентрации от высотного распределения Ы(2У).
Учет колебательного возбуждения молекулярного азота приводит к достаточно заметным вариациям электронной концентрации и температуры на высотах слоя Б2 ионосферы.
Например, при потоке солнечного радиоизлучения Бю.7 = 234 летом в полуденное время концентрация Б2-области уменьшается в 1,7 раза.
Показано, что вклад малых компонент в аэрономические процессы существенен. Учет зависимости скорости ионно-молекулярной реакции (И1) от колебательных состояний молекулярного азота улучшает согласие результатов модельных расчетов вариаций электронной концентрации с экспериментальными данными.
Работа выполнена при финансовой поддержке РФФИ по проекту № 14-01-00020а.
1. Павлов А. В., Намгаладзе А. А. Колебательно-возбужденный молекулярный азот в верхней атмосфере // Геомагнетизм и аэрономия. 1988. Т. 28, № 5. С. 705 — 721.
2. Pavlov А. V. The role of vibrationally excited nitrogen in the formation of the mid-latitude ionization through // Ann. Geophys. 1993. Vol. 11. P. 479 — 484.
но высокие концентрации n(21) относительно фоновых значений.
Заключение
Список литературы
3. Власов М. Н, Григорьев С. А., Ишанов С. А. и др. Сравнительный анализ различных гидродинамических приближений для описания ионосферно-магнитосфер-ной плазмы // Космические исследования. 1991. Т. 29, № 3. С. 404 —413.
4. Medvedev V. V., Ishanov S. A., Zenkin V. I. Effect of vibrationally excited nitrogen on recombination in ionosphere plasma // Geomagnetism and Aeronomy. 2003. Vol. 43, № 2. С. 248—255.
5. Ишанов С. А., Медведев В. В. Математическое моделирование метастабиль-ных компонентов в ионосфере Земли // Инженерно-физический журнал. Национальная Академия наук Беларуси. 2005. Т. 78, № 6. С. 26 — 33.
6. Власов М. Н., Ишанов С. А., Григорьев С. А. Моделирование эффектов динамических возмущений плотности и температуры ионосферно-магнитосферной плазмы на плазмосферных высотах // Космические исследования. 1997. Т. 35, № 3. С. 248 — 252.
7. Ишанов С. А., Леванов Е. И., Медведев В. В. и др. Магнитосферно-ионосфер-ные изменения, вызванные полетами космических аппаратов // Инженерно-физический журнал. 2006. Т. 79, № 6. С. 11—15.
8. Picone J. M., Hedin A. E., Drob D. P. et al. NRLMSISE-00 empirical model of the atmosphere: statistical comparisons and scientific issues // J. Geophys. Res. 2002. Vol. 107, NA 12. P. 1468-1483.
Об авторах
Сергей Александрович Ишанов — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: [email protected]
15
Леонид Викторович Зинин — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград. E-mail: [email protected]
Евгений Вячеславович Зубков — асп., Балтийский федеральный университет им. И. Канта, Калининград. E-mail: [email protected]
Сергей Валентинович Мациевский — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград. E-mail: [email protected]
About the authors
Prof. Sergey Ishanov — I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
Prof. Leonid Zinin — I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
Evgeny Zubkov — PhD student, I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
Dr Sergey Matsievsky — Ass. Prof., I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
