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

N2+ + ° 2 - N2 +0 j , k4 = 10-16 О 2 + n 2 - Ncf- N0, I to 10-21 N0++ е —* N + 0 ; о* = 1 4.1 10"13 О 2 + е —♦- 0 + 0 ; *2 * 2.2 10~13. Скорость ионов определялась из системы тр ех уравнений; ! г * i f r ! r ~ * Ni“i 9 - ( eEi-beC5i"'B] z)Ni * M|NiS»ji,(Vnr V i - е Е , - е [ ѵ . * в ] х . М і Г ) і п ( Ѵ п Г Ѵ і х ) і - = Е ѵ - е [ ? ; * Т ] у - M , S > i „ ( V - Ѵ Д (2 ) где Pj = N j kT j ; Pe = NekTe; Pj , — парциальные давления. ny P i- N e = £ N j — концентрация электронов, Mj — масса ионов; g — ускорение силы тяж ести ; Е — ^іектрическое попе; В — магнитное поле. В = ( - В cos I , О, - В s in І ) ; В = о. 16775-10-4 У і + 3 c o s 2 © ' , гд е I — наклонение магнитного поля; 6 ’ — дополнение к геомагнитной широте; Мп, Ѵп — масса и скорость нейтральных частиц; л 0 j п— частота столкновений ионов с нейтралами. Вертикальную составляющую скорости нейтралов полагаем равной нулю. Полученные из системы ( 2) значения составляющих скоростей ионов поп- являем в уравнение ( 1 ) . Получаем систему тр ех дифференциальных уравнений в частных про­ изводных Еторого порядка по координате z (по высоте) и первого порядка по време­ ни. Входящие в систему производные заменяю тся конечными разностями, и получен­ ная алгебраическая система решается методом прогонки На нижней границе концентрации ионов N 0 + и вычислялись из условия ; а- венства скоростей ионообразования и рекомбинации, а концентрация ионов 0 + полага­ лась равной нулю. На верхней границе потоки ионов N0+ и считались равными 46