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

налсначальногоидоконечногомо мента т . Приэтомзадачасводитсяк вычислениюискомойфункции N»N(z,t) навременномслое t приусловии, что известныеезначениянавсехпреды дущихслояхt- At, t -2 it,... At, 0 . Предположимтакже, чтопопространст веннойпеременной z областьотыскания решенияG разбитанарядслоевzB , Z jj — az , .., az +zH , Предлагаемый здесьспособзаданияверхнегогранич ногоусловиядляуравнения (II) под разумевает, чтооновычисляетсяв процессепереходакновымвременным слоям, иоснованнапримененииявной схемыповременидлявычисленияконцентрациизаряженныхчастицнаверхней границе, независимооттого, какаясхемаповременииспользуетсянавнутрен нихвысотныхслоях. Заметим, чтоалгоритмы, основанныенаприменениинеявных схемвовнутреннихузлахрасчетныхобластейвсочетаниисявнымсчетомгра ниц, успешноприменялисьприрешениизадачгазовойдинамики/4/. Такимобразом, значением наверхнейграницевычисляемвпроцессере шениязадачисиспользованиемявнойсхемыпоуравнению (5): РисЛ. Областьотысканиярешений G (z,t ). Выделеныучасткиграницы, на которыхставятсяначальныеигранич ныеусловия. H U B ,t+ At) . N U B ,t)+ At[-N*? 2 - V < f 3 +R0]Zi!it, (13) где 3Uz <f 3 (t). (14) Будемсчитать<f 2 (t) заданнойвеличинойтакже, каки^ ( t ) всоотношении (12). Величину<p 3 (t) вправойсчасти (13) такжебудемсчитатьизвестнойибудемза даватье следующимобразом. Рассмотримуравнениенеразрывности (5) встаци онарномслучае, когда эв ТЕ (15) иподставимвнегои из ( 10 ), тогданавысоте z„ получимсоотношение: ь В + bif 3 + (H<f 2 -R0 ) » 0 . (16) Используя(16) определимср3: Ь -b +j/b^-4a(if 2 - R 0 /H) (17) Очевидно, чтозасчетвыборазначения можнодобитьсяналичияу (16) хотя быодногоотрицательноговещественногокорня, таккакименноэтоготребуют физическиеусловиявовнешнейионосфере, гдеконцентрациязаряженныхчастиц убываетсвысотой. Подставляя найденноетакимспособомзначение<f 3 вмомент t вправуючасть (13), определимзначениеконцентрациизаряженныхчастицна верхнейграниценаследующемвременномслое t + At. Послеэтого, используя найденноезначениеи(г£ , t+ At), можемопределитьзначенияН навнутренних высотныхслояхвмомент t+ At. 108