Научная статья на тему 'Неустойчивость экваториального F-слоя ионосферы в условиях переменного электрического поля'

Неустойчивость экваториального F-слоя ионосферы в условиях переменного электрического поля Текст научной статьи по специальности «Физика»

CC BY
184
40
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИОНОСФЕРА / F-СЛОЙ / НЕУСТОЙЧИВОСТЬ РЭЛЕЯ-ТЕЙЛОРА / ЭЛЕКТРИЧЕСКОЕ ПОЛЕ / МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ / IONOSPHERE / F-LAYER / RAYLEIGH-TAYLOR INSTABILITY / ELECTRIC FIELD / MATHEMATICAL MODEL / NUMERICAL SIMULATION

Аннотация научной статьи по физике, автор научной работы — Кащенко Николай Михайлович, Мациевский Сергей Валентинович

Приведены результаты численного моделирования неустойчивости Релея-Тэйлора в условиях переменного электрического поля на основе электродинамически согласованной математической модели экваториального F-слоя ионосферы.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Кащенко Николай Михайлович, Мациевский Сергей Валентинович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

An instability of the equatorial F-layer of the ionosphere under conditions of alternating electric fields

The results of numerical simulation of Rayleigh-Taylor instability under conditions of alternating electric fields on the basis of electrodynamically consistent mathematical model of the Equatorial F-layer of the ionosphere.

Текст научной работы на тему «Неустойчивость экваториального F-слоя ионосферы в условиях переменного электрического поля»

УДК 519.6

Н. М. Кащенко, С. В. Мациевский

НЕУСТОЙЧИВОСТЬ ЭКВАТОРИАЛЬНОГО F-СЛОЯ ИОНОСФЕРЫ В УСЛОВИЯХ ПЕРЕМЕННОГО ЭЛЕКТРИЧЕСКОГО ПОЛЯ

Приведены результаты численного моделирования неустойчивости Релея-Тэйлора в условиях переменного электрического поля на основе электродинамически согласованной математической модели экваториального F-слоя ионосферы.

The results of numerical simulation of Rayleigh-Taylor instability under conditions of alternating electric fields on the basis of electrodynamically consistent mathematical model of the Equatorial F-layer of the ionosphere.

Ключевые слова: ионосфера, F-слой, неустойчивость Рэлея-Тейлора, электрическое поле, математическая модель, численное моделирование.

Key words: ionosphere, F-layer, Rayleigh-Taylor instability, electric field, mathematical model, numerical simulation.

1. Математическая модель экваториального F-слоя ионосферы

Наиболее подробное описание движения плазмы основано на использовании функции распределения молекул, удовлетворяющей уравнению Больцмана. Однако такое описание для ряда задач моделирования

© Кащенко Н. М., Мациевский С. В., 2014

Вестник Балтийского федерального университета им. И. Канта. 2014. Вып. 10. С. 30-35.

ионосферы слишком подробно, поэтому будем использовать уравнения, основанные на приближении Грэда [1]. Плазма на ионосферных высотах является тепловой с длинами свободного пробега частиц существенно меньше характерных пространственных масштабов изменения ее параметров. Согласно Гершману [2], общепринятую систему уравнений Максвелла и гидродинамических уравнений с учетом электромагнитных сил можно записать в виде, включающем п уравнений непрерывности, 3п уравнений движения ионов и электронов, п уравнений теплопроводности ионов и электронов, 3 уравнения непрерывности электрического тока и 3 уравнения потенциальности электрического поля [3]:

дп ■ -

—'-+у(пу1) = а -Ц,

дЬ - - ' '

^+(те=-пт+т{Ё+х ё)^^- *)^- *)+1,

3 JdTt

-n k 2 j

Л

-+(V V)T

dt K j ' '

+ PW + Vq^ = Gj - P,

_J_

v 6

V x E = 0,

V = e n V =0,

где n j, V, Qj, L j, m, e,, p,, Vjn, v p T , q , G , P j — соответственно концентрация, дрейфовая скорость, скорости образования и потерь, масса, заряд, давление, частоты соударений с нейтралами, частоты столкновений между заряженными частицами, температура, плотность теплового потока, скорость нагрева и скорость охлаждения частиц сорта j; k — постоянная Больцмана; j — плотность тока; E — напряженность электрического поля.

Вследствие замагниченности ионосферной плазмы F-области процессы переноса вдоль магнитного поля будут определяться столкновениями, а поперек поля — дрейфовым движением. Из-за сильной анизотропии, обусловленной магнитным полем Земли, процессы диффузии и теплопроводности в области F и во внешней ионосфере происходят в основном вдоль силовых линий геомагнитного поля. Выберем ди-польную систему координат (а, ф, р), так как магнитное поле Земли аппроксимируется дипольным приближением [4]. Координатные составляющие этой системы определяются через сферические координаты по формулам:

r _ cos 8 а = —^^ р = ——, Ф = Ф, sin 8 r

где r — расстояние до центра Земли; 8 — коширота; ф — долгота. Кроме того, при моделировании используем условия диффузионного приближения

dV -¡- = 0 dt

31

Н.М. Кащенко, С. В. Мациевский

32

и условие нейтральности

X е,п, =0

Помимо этого, будем считать магнитное поле постоянным по времени и дипольным [5].

Наконец, благодаря условию электростатики

Ух Е = 0

электрическое поле потенциально:

Е = -УФ,

где Ф — электрический потенциал. Плазма вдоль его геомагнитного поля на высотах экваториальной области Б является высокопроводящей средой, поэтому будем считать, что силовые линии геомагнитного поля эквипотенциальны, то есть

Ф = Ф(а, ф).

Понизим размерность трехмерного уравнения для потенциала, проинтегрировав его вдоль силовых линий. При этом область интегрирования учитывает параметры Е-области: ее силовые линии должны начинаться под областью Е в южном полушарии, проходить через геомагнитный экватор и заканчиваться под областью Е в северном полушарии. Тогда на концах силовых линий можно задать естественные граничные условия, заключающиеся в непроницаемости для электрического тока нижней границы Е-области.

Большинство достаточно развитых среднемасштабных экваториальных неоднородностей сильно вытянуты вдоль силовых линий и остаются такими при своем движении. Будем полагать эти неоднородности двумерными и описывать их динамику в плоскости геомагнитного экватора двумерными уравнениями. Тогда можно считать, что все переменные зависят только от координат а и ф, и в результате получаем двумерную модель развитых неоднородностей:

дп- -

+у—) = д,-е, ,

3 ,ГдТ.

—п к 2 '

Л

+ уЦ j + р..у±у. + уц. = с, -Р.,

У± (ст Ухф) = У± А.

Для вычисления параметров нейтральных частиц использовалась глобальная термосферная модель МБК [6 — 7].

2. Решение нелинейных систем уравнений

Уравнения вышеприведенной модели в двумерном приближении решались численно конечно-разностными методами на квазиравномерных сетках, сгущающихся к центру области решения. Начальные

значения задавались путем решения низкоширотной модели ионосферы до получения периодического решения. Эти начальные значения дополнялись искусственно введенными начальными возмущениями. Решение конечно-разностных аналогов исходных уравнений выполнялось по следующей схеме:

1. Сначала решался конечно-разностный аналог двумерного эллиптического уравнения, при этом использовался многосеточный метод до получения относительной погрешности 10-12.

2. Совместно решались уравнения непрерывности концентраций и теплопроводности.

Для улучшения качества совместного решения этих сильно связанных уравнений использовались итерации по нелинейности отдельно в каждом уравнении теплопроводности, по связям между уравнениями теплопроводности для разных сортов заряженных частиц и по связям между уравнениями непрерывности концентраций и уравнениями теплопроводности.

33

3. Результаты численного моделирования

Численное моделирование проводилось в двумерной области, ограниченной снизу высотой 100 км, а сверху — 1300 км, по горизонтали протяженность области интегрирования составляла 400 км. Сетка выбрана так, что в центральной области шаги равны 2 км.

Условия расчетов соответствовали среднему уровню солнечной активности с F10.7 = 150 и уровню геомагнитной активности кр = 3.

Для потенциала граничные условия заданы через фоновое электрическое поле. Для концентраций и температур на нижней границе заданы соответственно условия химического равновесия и равенства нейтральной температуре, вверху и на боковых границах задано условие равенства нулю потоков. Восточная компонента фонового электрического поля задавалась модельно: в начале процесса положительным значением (в представленных расчетах 1 мв/ м), а начиная с некоторого момента времени эта компонента меняла знак.

Была проведена серия расчетов для моментов перемены знака восточной компоненты электрического поля в моменты времени из диапазона 1800 — 2500 с, считая от начала процесса. Результаты вычислений показали, что неустойчивость Релея-Тэйлора развивается по классическому сценарию с выходом образовавшегося пузыря во внешнюю ионосферу для интервала времени перемены знака вертикального дрейфа, начиная от 1850 с. Для меньших интервалов знакоположительности фонового вертикального дрейфа развитие неустойчивости прекращается. На рисунке 1 показаны результаты моделирования процесса для времени изменения знака вертикального дрейфа 1850 с и 1900 с. Цифрами указаны значения десятичного логарифма концентрации (см-3).

Н. М. Кащенко, С. В. Мациевский

ь

км 1300

ь

34

горизонтальное расстояние, км

расстояние по горизонтали, кы

а б

Рис. 1. Стадия развитого пузыря для момента времени 3000 с при перемене знака вертикального дрейфа в момент времени 1900 с (а) и 1850 с (б)

Анализ результатов показывает, что изменение сценария развития неустойчивости Рэлея-Тейлора соответствует моменту времени, когда передний фронт ионосферного пузыря находится вблизи максимума Б-слоя ионосферы. Этот вывод остается также справедливым и при других условиях развития одиночных ионосферных пузырей в экваториальной Б-области.

Адекватность использованной математической модели и применяемых численных расчетов подтверждается многочисленными аналогичными расчетами авторов, а также зарубежными источниками, избранный список которых можно найти в [8].

Список литературы

1. Грэд Г. О кинетической теории разреженных газов // Механика. 1952. Вып. 4. С. 71-79 ; Вып. 5. С. 61-96.

2. Гершман Б. Н. Динамика ионосферной плазмы. М., 1974.

3. Кащенко Н. М., Мациевский С. В. Математическое моделирование неус-тойчивостей экваториального Б-слоя ионосферы // Вестник КГУ. 2003. Вып. 3. С. 59-68.

4. Рыбин В. В., Поляков В. М. Об амбиполярности движений ионосферной плазмы // Ионосферные исследования. 1983. № 33. С. 5—44.

Неустойчивость экваториального F-слоя ионосферы

5. Фаткуллин М. Н., Ситнов Ю. С. Диполярная система координат и ее некоторые особенности // Геомагнетизм и аэрономия. 1972. Т. 12, № 2. С. 333 — 335.

6. Hedin A. E., Salah J. E., Evans J. V. et al. A global thermospheric model based on mass spectrometer and incoherent scatter data MSIS 1. N2 density and temperature / / J. Geophys. Res. 1977. V. 82, N. A1. P. 2139—2147.

7. Hedin A. E., Reber C. A., Newton G. P. et al. A global thermospheric model based on mass spectrometer and incoherent scatter data MSIS 2. Composition // J. Geophys. Res. 1977. V. 82, N. A1. P. 2148—2156.

8. Мациевский С. В., Кащенко Н. М., Ишанов С. А. и др. 3Б-моделирование экваториального F-рассеяния: сравнение моделей MI3 и SAMI3 / / Вестник Балтийского федерального университета им. И. Канта. 2013. № 4. С. 102 — 105.

Об авторах

Николай Михайлович Кащенко — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград.

E-mail: kaschtschenko@mail.ru

35

Сергей Валентинович Мациевский — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград. E-mail: matsievsky@newmail.ru

About the authors

Dr Nikolay Kashchenko — Ass. Prof., I. Kant Baltic Federal University, Kaliningrad. E-mail: kaschtschenko@mail.ru

Dr Sergey Matsievsky — Ass. Prof., I. Kant Baltic Federal University, Kaliningrad. E-mail: matsievsky@newmail.ru

i Надоели баннеры? Вы всегда можете отключить рекламу.