Распределение электронов и физические процессы в полярной ионосфере / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Поляр. геофиз. ин-т. – Апатиты : [б. и.], 1981. – 144 с.

Длякаждого'значенияi (±= I, 5) числоуравненийв(13) равночислуузловсетки, которойпокрытаобластьи, т.е. числу неизвестныхпараметров , илиразмерностиискомоговекторанА- из(14). Однаковвидутого, чтонаграницеs расматриваемой областииискомыеконцентрациип. мызадаемкакфункци времени, чтоозначаетзадание величин (t) втехузлах( j,к), кото­ рыележатнагранице S, системауравнений(13) или, чтотоже самое,(14) являетсяпереопределенной. Поэтомумыпонижаемпоря­ доксистемы(13), илиразмерностьискомоговектора N. из(14), начислоузлов, лежащихнагранице s, припомощиприема, подроб­ ноописаногов/25/. Приэтомвидсистемыуравненийостается прежним. Длярешениясистемыобыкновеныхдиференциальныхуравне­ ний(14) применяемсхемупредиктор-корректор: N * +1/2- n ! Mi ---= V tlX “**«N5)-Ki(tl)Nb 1+1 1 М. N1 - V tl+1/2* ) H $ врезультатечегополучаем значенияN±(t) впоследовательныемо­ ментывременисзаданнымшагом«с. Возникающиеприэтомсистемы линейныхалгебраическихуравненийрешаемметодомматричнойпро­ гонки/30/. Приведемтеперьнекоторыерезультатырасчетов, полученые прирешенииописаннойвышезадачи. Будемсчитатьначалокорди­ нат0 двигающимсяпо поверхностиЗемливсеверномполушари та­ кимобразом, чтовкоординатахгеомагнитнаяширота- местноегео­ магнитноевремяточка0 неподвижна, причеме кординатыследую­ щие: геомагнитнаяширота85°, местноегеомагнитноевремя12 ч. СчитаемосьОхнаправленнойкэкваторуилежащейвплоскости геомагнитногомеридиана. Возьмемхл= 0, хп= 1000 км, = 100км, z = 420 км. Такимобразом, расматриваемаяобластьилежитв плоскости полуденогогеомагнитногомеридиана так, чтосечение дневногокаспадолжнопопадатьввыбраныепределыизменения координатых. Величинышаговсеткизадаемтакими: h1 = 50 км, ь2 = Ю км. Концентрацииотдельныхсортовнейтральныхчастиц задаемоднороднымипох, температурыотдельныхсортовнейтралов задаемодинаковымииоднороднымипох, аэкзосфернуютемпературу т'-1 принимаемравной1200°К. Скоростиотдельныхсортовнейтра­ ловзадаемодинаковымииоднороднымипох, причемвертикальную проекцию считаемравнойнулю. Скоростиобразованияионовв реакцияхфотоионизации(4) вычисляемдляусловийлетавсеверном 117