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

Пі ~ІЕ(піѵі )+ І (піті )'(іі * Li (5) напередзаданноймалойвеличины £. Расчетыначинаютсясмалымишагами повремени, прикоторыхискомыевеличиныустанавливаютсявсамомнизурас сматриваемоговысотногоинтервалаипочтинеменяютсявверхнейчасти. По достижениинанижнемслоеусловия полученноезначение наэтом слоесчитаетсяустановившимсярешениемификсируется, аслойисключаетсяиз дальнейшихвычислений. Вычислениянавышележащихслояхпроводятсяужесболь шимшагомповремени, поканаочередномслоенебудетдостигнутоусловие Послеэтогопроцедураповторяется. Этотспособопробованврасчетахпоопределениюконцентрациипятиупоми навшихсявышесортовионов. Максимальныешагиповременисоставилинесколько десятковсекунд, замораживаниераспространилосьвверхдовысоты310 км. Ди намикаизмененияотносительныхневязокиконцентрации иона0+ внекоторых узлахпространственнойсеткивходерасчетовпоказанавтаблице. Окончатель ныерезультатырасчетовпрификсированнойгоризонтальнойкоординатеX приве денынарис.2-4. Предлагаемыйспособпозволилсократитьчислошаговс73200, которыепришлосьбысде латьбеззамораживания слоевдляобеспечениянуж нойточностирешения, до 529 шагов. ДляЭВМ ЕС-1022 этипоказателисо ответствуютсокращению машинноговременис2200 до 16 ч. Помимозадач отыс каниястационарныхреше ний, используемыхнамиобычновкаче ственачальныхусловийдлямоделиру ющихповедениеионосфернойплазмы уравнений, приходитсятакжерешать эволюционныезадачипоотысканиюне стационарныхрешенийуравненийпри изменяющихсявовременикоэффициен тахиграничныхусловиях. Таккак характерныевременаизмененияуправ- новившиеся’(сплошная6лигая^Твысот^шеТа- лякщихпараметровионосферныхмоделей, — *--------------- — -----^ какправило, составляютнесколькоми нут/6/, топредотавляетинтересре шениенаотрезкахвременидесятки иболееминут. Примаксимальнодопус тимыхшагахповремени, которые, на пример, длясистемыуравненийнераз рывности (I) дляпятисортовионов навысоте100-400 кмсоставляютЮ-^с, С необходимоделатьболее10 шаговпри _ . „ / , использованииявныхразностныхсхем. Рис.4. Начальные (пунктир) иуста- * новившиеся(сплошнаялиния) высотные Намивыработанспособсокращения профиликонцентрацииионов N£ и ND+. количествашаговприрешениинестаци- Рис.2. Начальный (I) иустановившийся(2) высотныепрофиликонцентрации0+. профышконцентрацииионов0+ иN 52