Моделирование физических процессов в полярной ионосфере / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Поляр. геофиз. ин-т. – Апатиты : Кольский филиал АН СССР, 1979. – 148 с.

нулю, а поток ионов 0 + полагался постоянным, направленным вверх. Верхняя гра­ ница 420 хм выбрана так низко для то го , чтобы получать решения при разумной з а ­ трате машинного времени. Уравнения баланса энергии для определения Т и Т- не включены в решение. Мы полагали, что Те = Tj =. Т п , хотя в программе преду­ смотрен расчет зависимости температуры ионов от величины электрического поля и соответственно зависимости констант ионно-молекулярных реакций 0 + ♦ N 2 — - N 0 ++ N, 0 + + 02 ----О 2 + О от температуры ионов аналогично описанным в работе / 13 / . При расчетах исполь­ зовали модель нейтральной атмосферы Cl R A -7 2 /1 4 / для условий равноденствия и вы. сокой солнечной активности. Потоки солнечного излучения в диапазоне 40 — 1026 А и соответствующие сечения ионизации и поглощения для О, O j . N j взяты из работы / 15 / . Зависимость скорос­ ти ионообразования от зенитного угла Солнца рассчитывалась с помощью аппроксима­ ции функции Чепмена из работы /1 6 /. В отличие от работы / 8 / мы не полагали ско­ рость нейтрального ветра равной нулю. Система нейтральных ветров была взята по данным работы /1 7 / и схематично на уровне 300 км показана на рис. 2. На этом же рисунке показаны зоны высыпания авроральных частиц (электронов) на дневной и ночной стороне. При расчете были взяты следующие параметры потоков; на дневной стороне интенсивность потока І 0 = 5 . 4 ‘ 1о8 с м ~ 2 с ~ \ энергия электрона E Q= 150 эВ , на ночной стороне l Q = 1.2-10® с м ~ 2 с - 1 ( е Q = q . 5 к эВ . Скорость ионо­ образования потоком электронов рассчитывалась по формуле Лазарева /1 8 /. Одним из основных вопросов при постановке задачи является вопрос о выборе модели электрического поля, или схемы конвекции. Как уже отмечалось ранее, по­ явившиеся в последних работах модели электрических полей хотя в общих чертах и не противоречат друг другу, м о гут существенно отличаться как величинами электри­ ческих полей, так и направлением. Последнее особенно важно, так как при учете склонения магнитного поля зональная компонента электрического поля очень сильно влияет на профиль электронной концентрации. Сравнение модельных расчетов с экс­ периментальными данными в такой постановке задачи практически невозможно, так как получить мгновенную картину электронной концентрации в большом диапазоне ши­ рот и местного времени современные экспериментальные средства не позволяют. Для расчетов нами была выбрана модель электрического поля K a wa s a k i /1 0 /, которая позволяет по аналитическим формулам вычислять величину электричес­ кого поля а зависимости от геомагнитной широты (г _ расстояние от полюса) и местного времени X — д о л го та ). Эти формулы имеют следующий вид; r=A[a(^) 2 cosx 4 KFj) 2 + 1 } sinx]T - 4 a(^)2sinX+{b® 2+1} cosX]x для г< г 0;