CC BY
248
МАТЕМАТИКА
УДК 550.338
С. А. Ишанов, С. В. Клевцур
ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ ПРИ МОДЕЛИРОВАНИИ ВЫСОТНО-ШИРОТНЫХ РАСПРЕДЕЛЕНИЙ
Рассмотрено влияние смешанных производных при трехмерном моделировании ионосферы разностными методами.
A mixed derivatives influence in during three-dimensional simulation of ionosphere by difference methods is considered.
Ключевые слова: трехмерное моделирование, разностные методы, ионосфера, смешанные производные.
Key words: three-dimension simulation, difference methods, ionosphere, mixed derivatives.
В задачах моделирования ионосферной плазмы возникает необходимость отыскания периодических решений уравнения диффузии для различных компонент. Анизотропия среды и зависимость коэффициента диффузии от направления приводят к необходимости учета смешанных производных. Зафиксируем в трехмерном уравнении диффузии, записываемом в дивергентной форме в сферической системе координат в шаровом слое [1], географическую долготу, тогда получим двумерную нестационарную модель вида:
где N (і = 1, 2, 3, 4, 5) — концентрации ионов 0+, Н+, 02, К+, N0+ соответственно; и и V — меридиональные и зональные компоненты нейтрального ветра; Qi — функция ионообразования; аі — скорость потерь ионов в результате химических реакций.
Коэффициенты, входящие в уравнения, имеют следующий вид:
106
ИОНОСФЕРНЫХ ПАРАМЕТРОВ
(1)
(2)
(3)
Вестник Балтийского государственного университета им. И. Канта. 2011. Вып. 5. С. 106-110.
Вычислительный эксперимент при моделировании ионосферных параметров
Wu - RXX +
Wv -
т > rv
Pn r
pn rpn sin © 5X
Здесь pn — плотность нейтрального газа; Pn — давление; ^ — коэф-
возникающие за счет соударений нейтральных частиц с ионами.
Решение (1 — 3) ищем в области G - {r е [R0, R1], © е [©0, ©1], X - const} при t е [0, T], в данном случае t имеет смысл местного (локального) времени (LT), Г — граница области G. Дополним систему уравнений (1 — 3) начальными (t - 0) и граничными (r, ©) е Г условиями. Поскольку решение ищется в секторе, ограниченном коширотами ©0 и ©1, то краевые условия на северной (© - ©0 - 20°) и южной (© - ©1 - 70°) границах получим как решения упрощенной системы уравнений (1 — 3), в которой дополнительно зафиксирована географическая широта ©:
Система уравнений (1 — 3) совместно с соответствующими начальными и граничными условиями образует теоретическую нестационарную высотно-широтную модель ионосферы.
Разностные уравнения, аппроксимирующие систему (1 — 3), решались «а— р« итерационным методом [3], а разностные уравнения, аппроксимирующие систему (4 — 6), решались методом прогонки. Как следует из анализа уравнений (1), член при старшей производной в горизонтальном направлении имеет второй порядок малости по г, и можно предположить, что он не дает существенного вклада в распределение заряженных частиц.
Однако члены при смешанных производных имеют первый порядок малости по г и, таким образом, горизонтальная неоднородность ионосферы проявляется в уравнениях непрерывности для ионов через смешанные производные. Для оценки их влияния численно решались системы (1 — 3) с учетом и без учета этих производных. Корректировка неточности задания начальных условий проводилась методом установления до выхода решения на периодический суточный режим.
На нижней границе по г (г = Я0) задавались условия фотохимического равновесия ^ = аЫ). На верхней границе (Яг = 600 км) уравнение
фициент вязкости; Q — угловая скорость вращения Земли; Rni — силы,
107
(4)
(5)
(6)
С. А. Ишанов, С. В. Клевцур
108
непрерывности для ионов дополнялись заданием через нее потоков частиц, которые, в свою очередь, получались из расчетов по модели геомагнитной силовой трубки [2]. Потоки ионов на различных широтах получены путем интерполирования результатов расчетов для геомагнитных силовых трубок с параметрами Мак-Илвейна от Ь = 2,42 до Ь = 4,95 в узлы разностной сетки (гі, ©¡). Некоторые результаты численных расчетов представлены на рисунках 1 — 3.
Рис. 1. Суточная вариация логарифма электронной концентрации на кошироте 45°. Сплошные кривые соответствуют расчету без учета смешанных производных, пунктирные — с их учетом.
Зимний сезон, ?м,7 = 150, Кр = 3
На рисунке 1 показаны высотно-временные вариации логарифма электронной концентрации по кошироте © = 45°. Сплошные кривые относятся к расчетам без учета в уравнениях непрерывности смешанных производных, а пунктирные — к расчетам с учетом этих производных. Сравнительный анализ полученных результатов показывает, что влияние смешанных производных наиболее существенно проявляется в восходно-заходный периоды суток на высотах свыше 250 км, то есть в тех пространственных областях, где доминируют процессы переноса. Смешанные производные в уравнениях непрерывности через ионное торможение оказывают влияние и на нейтральный ветер, однако это влияние малосущественно (менее 10 %).
На рисунке 2 показаны высотно-широтные вариации электронной концентрации в 00 ч ЬТ (местного времени), то есть в то время, когда влияние смешанных производных наиболее значительно. Видно, что
Вычислительный эксперимент при моделировании ионосферных параметров
влияние производных практически не проявляется ниже 250 км. В то же время расхождение значений концентрации ионов, полученных в двух вариантах расчетов, уменьшается к северной (©0 = 20°) и южной (©і = 70°) границам области решения. Подобный эффект является результатом влияния краевых условий, которые, как уже отмечалось, задаются без учета горизонтальной неоднородности ионосферы. Данный факт служит свидетельством того, что кроме верхних граничных условий на электронную концентрацию в Б-области могут влиять и другие краевые условия при многомерном моделировании ионосферы.
Рис. 2. Высотно-широтные изолинии логарифма электронной концентрации с учетом (пунктирные кривые) и без учета (сплошные кривые) смешанных производных в 00 ч ЬТ. Зимний сезон, ґ10,7 = 150, Кр = 0
На рисунке 3 показано процентное относительное отклонение концентраций положительных ионов в максимуме Б2-слоя, полученных в двух вариантах расчетов для некоторых коширот. Видно, что смешанные производные максимально проявляются на кошироте © = 45°, то есть в средней части области решения, в околополуночные часы. Это позволяет сделать вывод о том, что при многомерном моделировании ионосферы необходимо сохранить смешанные производные по пространственным переменным в уравнениях диффузии ионов. Показано, что неучет может привести к 15%-ному отклонению в электронной концентрации Б2-слоя ионосферы на средних широтах в околопо-луночные часы. Влияние смешанных производных будет еще более значительным при моделировании ионосферных возмущений.
12 18 00 оь UlT
Рис. 3. Относительное процентное отклонение электронной концентрации в максимуме F2-слоя полученных в двух вариантах расчетов с учетом и без учета смешанных производных на некоторых коширотах (цифры у кривых)
Работа выполнена при финансовой поддержке РФФИ по проектам № 11-01- 00098-а и № 11-01- 00558-а.
Список литературы
1. Ишанов С. А., Клевцур С. В. Математическое моделирование ионосферы с учетом ее трехмерной неоднородности // Вестник Российского государственного университета им. И. Канта. 2010. Вып. 10. С. 152—158.
2. Латышев К. С., Зинин Л. В., Ишанов С. А. Математическое моделирование околоземной космической плазмы / / Энциклопедия низкотемпературной плазмы. 2008. Т. 7-1, ч. 3. С. 337—349.
3. Ишанов С. А., Клевцур С. В., Латышев К. С. Алгоритм а — Р итераций в задачах моделирования ионосферной плазмы // Математическое моделирование. 2009. Т. 21, № 1. С. 33—45.
Об авторах
Сергей Александрович Ишанов — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта, e-mail: math@dekan. albertina.ru.
Сергей Владимирович Клевцур — канд. физ.-мат. наук, доц., Балтийский федеральный университет им. И. Канта.
Authors
Dr. Sergey Ishanov — assistant professor, I. Kant Baltic Federal University, e-mail: sergey.ishanov@ yandex.ru.
Dr. Sergey Klevtsur — assistant professor, I. Kant Baltic Federal University.
