CC BY
74
УДК 550.388.2
Н. М. Кащенко, С. А. Ишанов, С. В. Мациевский
РАЗВИТИЕ НЕУСТОЙЧИВОСТИ РЭЛЕЯ - ТЕЙЛОРА В ЭКВАТОРИАЛЬНОЙ ИОНОСФЕРЕ ИЗ СТОХАСТИЧЕСКИХ ИОНОСФЕРНЫХ НЕОДНОРОДНОСТЕЙ
Рассмотрены процессы развития неустойчивости Рэлея — Тейлора в экваториальной области ионосферы из случайно заданных малых возмущений нижней части F-слоя ионосферы.
Processes of development of instability of Rayleigh-Taylor in the equatorial region of the ionosphere from accidentally given small perturbations of lower part of F-layer of the ionosphere.
Ключевые слова: численная модель, неустойчивость Рэлея-Тейлора, плазменные пузыри.
Key words: numeral model, the Rayleigh-Taylor instability, plasma bubbles.
Введение
Известно, что в ионосферной плазме на фоне регулярных распределений электронной концентрации имеются случайные возмущения с характерными масштабами от ~10 м до ~100 км. Особенности таких возмущений исследованы в ряде работ, например в [1; 2], на основе данных вертикального зондирования ионосферы при измерениях флуктуаций уровней радиоизлучения звезд, а также при использовании спутниковых измерений. В экваториальной ионосфере одной из причин появления среднемасштабных неоднородностей в вечернее время стало усиление малых неоднородностей механизмом неустойчивости Рэлея — Тейлора, приводящее к появлению ионосферных пузырей [3].
1. Численная модель экваториального F-слоя ионосферы
В данной работе приведены результаты численных экспериментов по развитию плазменных пузырей из малых начальных возмущений со случайными параметрами. Для проведения численных экспериментов использована пятиионная модель экваториальной F-области ионосферы в двумерном диффузионном приближении, уравнения которой имеют вид [4]
dn, ~dt
( дТ. _
dt
—n.k 2 j
_L + V±(n,V,) = Q, - Lj, (1)
- + V,V±)T + P,VV, +V±q, = G, -P,, (2)
13
V±(aV±U) = V± A, (3)
© Н. М. Кащенко, С. А. Ишанов, С. В. Мациевский, 2015
Вестник Балтийского федерального университета им. И. Канта. 2015. Вып. 10. С. 13 — 17.
14
В1=-УАи, (4)
Е + У] хВ = 0, (5)
где щ, V -, Qj, Ь, р, Т, у - С, Р,) — соответственно концентрация, дрейфовая
скорость, скорость образования и потерь, давление, температура, плотность теплового потока, скорость нагрева и скорость охлаждения заряженных частиц сорта р; к — постоянная Больцмана; Е — напряженность электрического поля; ст — тензор проводимостей; А — источники тока; и — потенциал электрического поля; В — вектор индукции магнитного поля Земли. В модели использованы ионы О2, ,N2,, N0+ ,0+ ,Н+. Для вычисления внешних параметров, входящих в коэффициенты уравнений, использовалась глобальная термосферная модель МБ1Б [5—6].
Начальные значения при решении этих нестационарных уравнений задавались решением уравнений, полученным в результате численных расчетов на установление при заданном электрическом поле и при использовании глобальной модели термосферы МБ1Б для задания параметров нейтральной атмосферы.
Уравнения модели (1) — (5) решались в экваториальной плоскости конечно-разностными методами второго порядка точности на квазиравномерных регулярных сетках, шаги сеток в центральной части области были взяты равными 1 км. Для уравнений непрерывности ионов (1) и уравнений теплопроводности (2) использовались разностные схемы ТУО типа [7]. При этом для функции /(г), ограничивающей потоки [7], выбран следующий вид:
' гТ2 Л
Дг) = тах
0,
VI + г2
В соответствии с [7] такой ограничитель находится в области нелинейных монотонных схем с ТУО свойством.
Уравнение потенциала электрического поля (3) имеет эллиптический тип с несамосопряженным оператором и переменными коэффициентами и решается на каждом временном шаге стандартным многосеточным методом. Для получения относительной погрешности порядка 10-12 при решении этим методом в среднем требовалось примерно 30 итераций.
Для улучшения качества совместного решения этих сильно связанных уравнений использовались итерации по нелинейности отдельно в каждом уравнении теплопроводности, по связям между уравнениями теплопроводности для разных сортов заряженных частиц и по связям между уравнениями непрерывности концентраций и уравнениями теплопроводности.
2. Результаты численных экспериментов
В работе [8] было проведено исследование конкуренции небольшого количества среднемасштабных неоднородностей на основе двумерной модели. Было показано, что при достаточно малых расстояниях экваториальные ионосферные пузыри взаимодействуют между собой,
что приводит к изменению параметров развития и даже подавлению некоторой части пузырей. Более того, начальные малые неоднородности концентрации плазмы могут иметь различные причины возникновения и различные параметры.
Для моделирования этой ситуации в данной работе рассмотрены процессы развития системы ионосферных пузырей в условиях, когда в некоторый начальный момент времени на фоновую концентрацию нижней части Б-области экваториальной ионосферы наложено поле случайных малых возмущений. Пусть у — горизонтальная, а 2 — вертикальная координата области интегрирования. Тогда начальные малые возмущения будем задавать по формуле
п, (У, 2) = п, о( У, 2)
1 - М У )ехр
( ( 2 * (У)А2 ^
2 - /2( У) /э( У)
//
где п,о(У, 2) — невозмущенные ионные концентрации; /\(у), /2(у), /з(у) — функции, описывающие структуру возмущения. Значения этих функций задавались псевдослучайными с помощью встроенного в язык программирования С++, на котором реализована модель генератора псевдослучайных чисел, использующего как известно линейный конгруэнтный метод. Кроме этого в отдельных сравнительных расчетах был использован метод Фибоначчи с запаздываниями. Плотности вероятности для /2 и /3 показаны на рисунке 1, а значения /1 распределены равномерно на отрезке [0; 0,2].
Рис. 1. Плотности вероятностей для значений функций /2 (слева) и /3 (справа)
15
Численные эксперименты проводились для условий, соответствующих равноденствию, уровню солнечной активности Р10.7 = 150 и уровню геомагнитной активности кр = 3, при этом начало процесса соответствовало 18.30 местного времени. Была проведена серия численных экспериментов с разными последовательностями псевдослучайных чисел и на конечноразностных сетках с различными пространственными шагами.
На рисунке 2 показано распределение электронной концентрации для одного из численных экспериментов в моменты времени 2200 с, 2600 с, 3000 с, 3400 с, 3800 с и 4000 с, считая от начала процесса. На этом
рисунке по горизонтали отложена долгота относительно центра области, а по вертикали отложена высота. Для лучшего отображения картинки вытянуты по горизонтали.
16
Рис. 2. Распределение десятичного логарифма электронной концентрации в различные моменты времени от начала процесса
Анализ результатов проведенных экспериментов показывает, что в описанных условиях процесс развития идет по следующему сценарию. На начальной, линейной, стадии процесса, которая в условиях экспериментов имеет длительность 1400 с, в области интегрирования (протяженностью по долготе 800 км) формируется примерно 80 зон пониженной концентрации. На следующей нелинейной стадии развития неустойчивости Рэлея — Тейлора из этих зон формируется примерно 40 пузырей. Наконец на последней стадии эксперимента, соответствующей выходу наиболее развитых пузырей во внешнюю ионосферу, остается около 20 конкурирующих между собой развитых пузырей.
Таким образом, на первой стадии плотность неоднородностей составляет примерно одна неоднородность на 10 км, на второй стадии плотность пузырей составляет примерно один пузырь на 20 км, и наконец на последней стадии плотность развитых пузырей составляет примерно один пузырь на 40 км. Расчеты, проведенные с другими псевдослучайными последовательностями и с другими разрешениями разностных сеток, дают для указанных гелио-геофизических условий примерно такие же плотности ионосферных пузырей.
Кроме этих количественных особенностей процесса следует выделить следующие качественные особенности:
— сильная долготная асимметрия пузырей, существенно большая, чем при развитии одиночных пузырей;
— возникновение изолированных пузырей существенно раньше, чем при развитии одиночных пузырей, что приводит в дальнейшем к быстрому развитию отдельных менее развитых пузырей.
Работа выполнена при финансовой поддержке РФФИ по проекту 14-01-00020.
Список литературы
1. Гершман Б. Н, Ерухимов Л. М., Яшин Ю. Я. Волновые явления в ионосфере и космической плазме. М., 1984.
2. Гершман Б. Н., Казимировский Э. С., Кокоуров В. Д., Чернобровкина H.A. Явление F-рассеяния в ионосфере. М., 1984.
3. Ossakow S. L., Zalesak S. Т., McDonald B. E., Chaturvedi P. K. Nonlinear equatorial spread-F: Dependence of altitude of the F-peak and bottomside background electron density gradient scale length // J. Geophys. Res. 1979. Vol. A84, № 1. P. 17-39.
4. Кащенко Н. М., Мациевский С. В. Неустойчивость экваториального F-слоя ионосферы в условиях переменного электрического поля // Вестник Балтийского федерального университета им. И. Канта. 2014. Вып. 10. С. 30 — 35.
5. 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. Vol. 82, № A1. P. 2139-2147.
6. 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 // ibid.
7. Сафронов А. В. Оценка точности и сравнительный анализ разностных схем сквозного счета повышенного порядка // Вычислительные методы и программирование. 2010. Т. 11. С. 137—143.
8. Кащенко Н. М., Мациевский С. В., Никитин М. А. Динамика системы множественных рэлей-тейлоровских ионосферных пузырей // Геомагнетизм и аэрономия. 1990. № 2. С. 281—286.
17
Об авторах
Николай Михайлович Кащенко — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград. E-mail: [email protected]
Сергей Александрович Ишанов — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Калининград. E-mail: [email protected]
Сергей Валентинович Мациевский — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, Калининград. E-mail: [email protected]
About the authors
Dr Nikolay Kashchenko — ass. prof., I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
Prof. Sergey Ishanov — I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
Dr Sergey Matsievsky — ass. prof., I. Kant Baltic Federal University, Kaliningrad. E-mail: [email protected]
