Математическое моделирование комплексных процессов / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1982.

" / * № - * - !&И £ )* ( 2 1 ) 3 3 ( 2 3 ) п=»1 где у. —коэффициент вязкости ионов, Xj и \ е - коэффициенты теплопроводности ионов и электронов, I - угол наклона магнитно го поля, М = 3 / 2 X N , 5 / 3 ; остальные величины определены выше. При записи ( 2 1 ) - ( 2 2 ) использованы выражения тензора вяз ких напряжений ионов и векторов тепловых потоков ионов и электро нов из / 1 5 / с учетом сильной замагниченности заряженных частиц; при выводе коэффициентов вязкости и теплопроводности ионов у. и A,j нами учтены только упругие столкновения ионов с ионами. При выведении коэффициента теплопроводности электронов Хе учиты вались упругие столкновения электронов с электронами и нейтраль- Шыми частицами. Определяемые системой уравнений ( 2 0 ) - ( 2 3 ) искомые величи ны N s V„ Т | , Те являются функциями двух переменных: времени t и расстояния от земной поверхности вдоль силовой линии h . Время t с использованием скорости U c может быть связано с расстоя нием, проходимым началом подвижной системы координат вдоль ли нии конвекции. Поэтому отыскание нестационарных решений системы { 2 0 ) - ( 2 3 ) равносильно отысканию величин N , V, Т] , 1 е как функций двух пространственных переменных: расстояния от земной поверхности Ь и расстояния вдоль линии конвекции. Поскольку можно взять большое число линий конвекции, чтобы они покрывали собой земную поверхность достаточно густо и отстоял*, друг от дру га не слишком далеко, тс это позволяет фактически построить про странственно трехмерное распределение искомых величин N , V, Т- , Те над земной поверхностью. Таким образом, физические условия позволяют задачу об отыскании пространственно трехмерных решений системы уравнений переноса ионосферной плазмы расщепить на две значительно менее сложные задачи: решение зависящих от двух пе ременных уравнений переноса; отыскание линий конвекции и скорос ти движения плазмы вдоль них. Нам удалось разработать методику численного решения систе мы уравнений ( 2 0 ) - ( 2 3 ). Эти уравнения решаются относительно искомых величин N 9 V, Ті , Те численно методом конечных разнос тей в пределах расстояний h от земной поверхности от 1 0 0 до 12