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

+ (23) гдеѵ- продольнаякомпонентаскорости; ju,- коэффициент е я з к о с т и заряженных частиц. Вуравнениинеразрывности (21) первыйчленвправойчастиявляетсячленом сжатия (илиразрежения), второй- конвективнымчленом. Предположим, чтонадоотыскатьрешениесистемы (2І)-(23) вобласти T(b,t)»[h H4 h 6 h ^ x j p ^ t 4 t]. Таккакучитывающеесилыинерцииивязкости уравнениедвижения (22) иуравнениетеплопроводности (23) являютсядифферен циальнымиуравнениямивчастныхпроизводныхвторогопорядкапараболического типа, длянихкорректнойявляетсяпостановкасмешаннойкраевойзадачисна чальнымиусловиямипри t»0 играничнымиусловияминанижнейиверхнейграни цахhH иhg. Будемсчитать, чтонанижнейграницеизвестныскоростиv(hH ,t) итемператураT(hH ,t) какфункциивремени, анаверхнейгранице- производ ные: Уравнениенеразрывности (21) являетсяквазилинейнымдифференциальным уравнениемвчастныхпроизводныхпервого порядка, интегрированиекоторого эквивалентнорешениюхарактеристической системыобыкновенныхдифференциаль ныхуравнений: Каждое решениеэтойсистемызадаетхарактеристикувпространствепеременных t. h, N /8/. Признакопеременнойпродольнойскоростиѵ задачадляуравнения неразрывности (21) вобластиY(h,t) являетсяпоставленнойкорректно, если задаются значениярешения на тех, итолькотехучасткахграницыобласти X(h,t), скоторыххарактеристикиидутвнутрьэтойобласти/8/ (рис.2). Мыжеобсудимтолькоконечно-разностныеметодырешенияоистемы(21)— (23 когдаэтирешенияотыскиваютсякакфункцииh нанекоторойпоследовательно стивременныхслоев, начинаяснулевогоикончаяконечныммоментомт. Поэто муалгоритмотысканияреиенияуравнениянеразрывностинанекоторомвременном слоеt долженпредусматриватьпроверкузнакаѵ науровняхhj и^ ив зави симостиотзнака(рис.2) использоватьобаграничныхусловиянаhH и t^c) либотольконаhgCtg6 t ^т^), либотолькона 1 ^ ( 0 4 либонеиспользоватьсовсем(т^ t 6 i 2) дляотысканияконцентрациивовну треннихузлахобластиy(h,t). Независимооттого, какаясхемаповремени применяетсядлявычислениязначенийК вовнутреннихузлахииспользуютсяли граничныезначенияннаданномвременномслое, будемвычислятьзначение нанижнейграницеизусловия фотохимическогоравновесия, анавер хнейграницепоуравнению(21) появнойсхеме: (24) 110