Распределение электронов и физические процессы в полярной ионосфере / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Поляр. геофиз. ин-т. – Апатиты : [б. и.], 1981. – 144 с.
СкоростьпотерьприхимическихреакцияхL вноситнелиней ностьвуравнениенеразрывности(3). Приинтегрировани уравне нийгазовойдинамикиотрешениекабсолютноустойчивым, диверген тнымиточнымсхемамприводиткнеобходимостиитерацийпонели нейности. Существуетмнение, чтометодпредшстор-корректорпозво ляетполучатьабсолютноустойчивыедивергентныесхемывторого порядкаточности, неприбегая китерациям. Вприменени курав нениямвчастныхпроизводных онрассматривалсявработах/2-4/ идр. Сформулируемэтотметодприменительнокуравнениюнеразрыв ности(3). Напервомэтапе(предиктор) определяетсярешениеи*, отнесеноекнекоторомумоментувремениt*: tn < t* < t n+1 с помощьюнекоторойабсолютноустойчивойаппроксимации, недивер гентной s скоэффициентами, определенымидлямоментавремени t=tn. Послеэтогоприменяетсякорректор, т.е, схема, основаная наприменени законасохранения. Приреализаци этогометодавоз можныразличныевариантыпредиктораикорректора. Дляконкретныхфизическихусловий, осуществляющихсявF-об ластиполярнойионосферы, которыеподробноописанывработе/I/, былипроведены расчетысразличнымивариантамидивергентных схемсцельювыборанаиболе простыхвреализациииэкономичных всмыслезатратмашиноговремени. Рассмотримрезультатырасчетов когдапродольнаяскоростьнаправленавверх. Всответстви сфизи ческойпостановкойзадачи, подробноописанойв/I/, мыотыскива емстационарноерешениеуравнениянеразрывности(3) методомуста новления. Стационарноерешениесчитаетсянайденымивычисления прекращаются, когдамаксимальнаяотносительнаяневязкастановит сяменьшенапередзаданойвеличины, какэто продробноописано вработе/5/. ВрасчетахпродольнаяскоростьV являетсязаданойфункцией; е высотный профиль представленкривой4 нарис. I работы /I/. ВкачественачальногоприближенияKQ выбранагладкаяфункция, имеющаямаксимумнавысоте300 хш. Напервомэтаперасчетыпрово дилисьподвухшаговойдивергентнойсхеме: вспомогательныйшаг _ t _ l _ + -----Ш _----U + . Qol , ( 4 ) навторомшагеприменяетсякорректор: н*+ 1 -ы* (HV)* 1 -(м)! 1 * t+1 + tt V t .. f c l + . ,о1. (5) где f* вit. 29
Made with FlippingBook
RkJQdWJsaXNoZXIy MTUzNzYz