Хрущинский А. А. Распространение гидромагнитных волн в неоднородной магнитосфере / Хрущинский А. А., Пудовкин А. И., Мальков М. В. ; АН СССР, Кол. фил. им. С. М. Кирова, Поляр. геофиз. ин-т. - Препринт ПГИ-88-01-61. - Апатиты, 1988.

SdhHm .nJ- 'JU- l m,/i)+[-al(SsU.m+1n)-jy (i,,n-i, n ))/2 & + *Q}(b„(i,iri,n+i)-l)g(i,m,n~l))/2.8y+уJг -ZSt Vl (i*J i m ,n )-V 2(L-i, m, n )* [ o f’fp j& i n ) ) 2 S y + (6r) +py (h (L ,m ,n + l)- 62(L ,m ,ii--t))/2 fy -CLz -2 Ti<(pI,&t (L,tv/ n )+ (6ц) *fyij (i,m,n) (pxjj -Pyjx)]2$t. Исследование устойчивости вычислительной схемы с уравнениями ( 6 ) приведено в приложении. Следует отметить, что приведенный алгоритм требует достаточной гладкости функций, описывающих источник: в противном случае быстро растет вычислительная неустойчивость. Для обеспечения требуемой глад кости приходится соответственно увеличить размеры области источника (по числу узлов сетки), что, в свою очередь, влечет увеличение сеточ ной области, в которой ищется решение. Так, минимальный приемлемый размер области определения решения как показали оценочные расчеты, составляет 200 х 200 узлов сетки. Если учесть, что каждый узел описы вается пятью комплексными амплитудами, то весь объем данных составля ет 3,2 мегабайта, что превышает реальные возможности доступных вычис лительных средств. Поэтому алгоритм был построен таким образом, что бы сохранить основную часть информации в файле прямого доступа на диске, а не в оперативной памяти. Каждая запись содержит все данные одной строки сеточного квадра та для моментов времени и £ . В оперативной памяти хранятся данные только трех последовательных записей для t и одна запись для t -Si , что позволяет по проведенным выше формулам определить новые значения полей для средней из трех строк для t *SH . Эти новые значе ния записываются в файл, обеспечивая таким образом его постепенное построчное обновление. Конкретный вид функций, описывающих источник, следующий. J ct ~Jcto ■f(по) ■ ) ■ ехр(скг) где: при ^ ' при 6