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

3. АЛГОРИТМЫ ЧИСЛЕННОГО РЕШЕНИЯ. Раснишвм систему (2) в координатах: Vx= ~а- Р</ ~ > (За) Чу-сёрх (Чх&у &x)~<fPx jz' * ( 36 ) Уг*-а-%(L*&* " Ъ 4 УО- %(Ъ&г ~L (P*]d -PJJ (зв) £Х~ Ц ,( И, ру ~VyPx)~i-KЪ рх ’ (Зг) ЧгР* -V(Vxfr~ V*p,)> (ЗД) / г = V* П 4 Р х )" Vy ( Иг P J > t 3 e) где: о , -Ща.1/с , а у компонент плотности стороннего тока jeT опущен индекс "ст ". Полученная система уравнений (3) является основой для построения разностной схемы, приближавдей дифференциальную задачу. Для создания более экономичного алгоритма произведем некоторые преобразования системы (3 ). Умножим уравнение (За) на /э и вычтем одно из другого. Получим уравнение: р_ I. •, (36) на W = —< 2 ^ ( Ч ц Sly1х) + J 2 (4а) где W * \/хру -VyPx. Уравнения (Зв) - (Зе) перепишем, используя последнее обозначение Йг = О'ЧРх ч &г *ру Vy t>i)-CLzi« (Рх£* *R ,b (46) b = Vy W 'i * V z P x (4в) Ду z~y,\J~iK Vt Ру (4г) Ьг = ^ ( ^ P x ) r ^y(PyV2)- (4д) В полученных уравнениях отражена ортогональность вектора скорос ти V к невозмущенному магнитному полю В0, так что три неизвестные компоненты скорости V„ , V, , Уг заменены на две: W и V2 . Уравнения системы (4) позволяют непосредственно приступить к построению разностной схемы. Разностный аналог системы (4) был выбран явным, трехслойным по времени / 4 / . 4