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

(<3to - полноесечениеупругогорассеяния), чтосоответствуетрассеяниюна малыеуглы. Хотяэтопредположениефизическиоправданолишьдлявысокоэнер гичныхэлектронов, установлено/12-14/, чтоприрешениикинетическогоурав нениядлявысыпающихсяиубегающихэлектроновзаменаизотропноготипаупру гогорассеяниянакулоновскийнедаетбольшихотличийвопределенииоснов ных характеристиквобластиэнергииЕ> I эВ. Конечно-разностнаяаппроксимацияпроизводныхпо s и Е сводитуравне ние (I) квиду: - ф - + А Ц - СФ = D, (5) гдекоэффициентыА,С, D являютсяфункциямипеременных jx., s,E. Приэтом, од нако, необходимообеспечитьвыполнениеусловияустойчивостиС^о. Последнее достигаетсяаппроксимациейтипа: (Ф--Ф” 8 ) / д а , >1>0 8 Ф Те ( 6 ) (ф+а-ф)/дв, О Ш - • <’ > где<р- различныефункции, присутствувдиеввыражениях(I), (2), (3), ипод боромтакихшагов дай дЕ, прикоторыхвыполняетсянеравенство as саЕ/|е|. Граничныеусловияпо а иЕимеютпростойвид: Ф| = 0 ; Ф"а = 0 ; Ф+в к > в Г max 3“ 3В1ІП J*> о ( 8 ) Д<0 гдеiff определяетсяэнергетическимипитч-угловымспектромвторгающихсяча стицнаверхнейгранице. Областьпитч-углов д<0 соответствуетнаправлению движениявниз. Граничныеусловияпоугловойпеременнойнаходятсяприпод становкезначенийju.= +1 вуравнение (5) сучетомконкретноговидакоэффи циентовA,C, d : = ЭФ = 0 , (9) Ѳ =0 30 Ѳ = УГ ЭФ аѳ гдеѲ = arcсоа Ji. Дополнительныеусловия, необходимостькоторыхсвязанасраздельнымрешением уравнения ( 5 ) для ,#?0 иЛ <0 (выражение ( 6 ), можнополучить, еслифункция ие производнаяпри >=0 будутнепрерывными. Алгоритмчисленногорешенияуравнения(5) основаннапримененииконеч- но-разностнойсхемы, сводящей(5) ксистемелинейных алгебраическихурав ненийснеизвестнымиФ^, гдеиндексксоответствуетнекоторомузначению питч-углаѲ. Даннаясистемарешаетсянакаждомe-слоеметодомпрогонки. В областии> 0 вычисленияпроводятсяприпоследовательномувеличениикоордина тыя, авобласти jti< 0 - припоследовательноме уменьшении. Длякаждого значения Евычисленияповторяютсявитерационномцикледотехпор, покана всехуровняхневыполнитсяусловиесовпадениярешениянаграницеобластей ( > = 0). Полученноерешениезапоминаетсяииспользуетсянапоследующем (ни жележащем) энергетическомуровне. Исходныеданныеирезультатымодельныхрасчетов. Сечениянеунругих столкновенийэлектроновсатомамиимолекуламирассчитывалисьнаосновера бот / 1 5 , 1 6 /. Дифференциальныесеченияударнойионизацииаппроксимировались 20