Вестник КРАУНЦ. Физ.-мат. науки. 2022. Т. 41. №4. C. 107-119. ISSN 2079-6641
УДК 519.254, 517.958:537.84
Научная статья
Динамика изменения режимов генерации магнитного поля в зависимости от частоты осцилляций процесса подавления а-эффекта энергией поля в модели аО-динамо
О. В. Шереметьева
Институт космофизических исследований и распространения радиоволн ДВО РАН, 684034, Камчатский край, с. Паратунка, ул. Мирная, 7, Россия E-mail: sheremeteva@ikir.ru
В статье рассматривается крупномасштабная модель аО-динамо в маломодо-вом приближении. Интенсивность а-эффекта регулируется процессом с эреди-тарными свойствами («памятью»), зависящим от энергии магнитного поля. Характер воздействия процесса определяется знакопеременным ядром с варьируемой частотой затухания и постоянным коэффициентом затухания равным 0.1. По результатам численного моделирования режимов генерации магнитного поля определены ограничения для значений параметров модели, при которых нет осцилляций в поле скорости вязкой жидкости или они незначительны. Результаты численного моделирования режимов генерации магнитного поля при различных значениях частоты затухания отражены на фазовой плоскости управляющих параметров, в которых заложена информация о крупномасштабном и турбулентном генераторах. В работе исследуется вопрос о динамике изменения картины на фазовой плоскости в зависимости от значения частоты затухания и времени ожидания. Проводится сравнение с результатами, полученными ранее при постоянной интенсивности а-эффекта и при её изменении, которое определяется показательным ядром.
Ключевые слова: маломодовая модель динамо, аО-динамо, а-эффект, магнитное поле, инверсии.
Для цитирования. Шереметьева О. В. Динамика изменения режимов генерации магнитного поля в зависимости от частоты осцилляций процесса подавления а-эффекта энергией поля в модели аО-динамо // Вестник КРАУНЦ. Физ.-мат. науки. 2022. Т. 41. № 4. C. 107-119. d DOI: 10.26117/2079-6641-2022-41-4-107-119
Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.Org/licenses/by/4.0/deed.ru)
© Шереметьева О. В., 2022
Финансирование. Работа выполнялась в рамках государственного задания по теме «Физические процессы в системе ближнего космоса и геосфер при солнечных и литосферных воздействиях» (No АААА-А21-121011290003-0).
d DOI: 10.26117/2079-6641-2022-41-4-107-119
Поступила в редакцию: 29.11.2022
В окончательном варианте: 06.12.2022
Введение
Данная статья продолжает цикл работ ИКИР ДВО РАН по исследованию мало-модовой модели динамо земного типа [1]-[8]. Задача исследования различных приближений на предмет генерации геомагнитного поля, близкого к наблюдаемому, несмотря на успехи, достигнутые в данной области, остаётся актуальной. Земное динамо хорошо согласуется с моделью аП-динамо, характерными свойствами которого является сильное дифференциальное вращение объекта и турбулентный характер движения проводящей среды [4],[9]-[13]. Характерным свойством динамо-системы земного типа является наличие хаотических инверсий без существенной перестройки движения проводящей среды [14]. В модели аП-динамо, которая описывается магнитогидродинамическими (далее МГД) уравнениями [2, 4], включающими уравнение Навье-Стокса, уравнение индукции магнитного поля В с учётом турбулентного а-эффекта, условия неразрывности поля скорости V и соленоидаль-ности магнитного поля, а также граничные условия прилипания для поля скорости V, генерация магнитного поля и поля скорости определяется работой П-гене-ратора и а-генератора и интенсивность воздействия последнего изменяется во времени под влиянием внешних сил. В рамках данной работы эта детерминированная динамическая модель рассматривается с аддитивным регулятором интенсивности а-эффекта в виде процесса с конечной памятью, которая определяется знакопеременной функцией |("Ь) = еЪ1соза± с постоянным коэффициентом затухания Ь = 0.1 и различными значениями частоты а затухающих колебаний [5]-[8].
Для исследования используется перемасштабированная и обезразмеренная МГД-система [6]. В качестве единицы времени принято время затухания магнитного поля. Использование такого временного масштаба позволяет смоделировать основные динамические характеристики поля с помощью маломодовых пространственных приближений. В работе приняты одномодовое приближение для поля скорости вязкой жидкости и трёхмодовое - для магнитного поля, что обеспечивает работу динамо с инверсиями в магнитном поле.
Целью данного исследования маломодовой модели аП-динамо является определение параметров модели, при которых могут быть получены режимы генерации магнитного поля как регулярного, так и нерегулярного характера, в том числе с инверсиями, на фоне слабо изменяющегося поля скорости. Сравнительный анализ результатов численного эксперимента для допустимых параметров модели с ранее полученными результатами при постоянной интенсивности а-эффекта и интенсивности, изменяющейся по показательному закону с заданным постоянным значением коэффициента затухания.
Математическая модель
Подробное описание перемасштабирования и обезразмеривания МГД-системы проведено в статье [6], поэтому в рамках данной работы приведём уже преобразованную систему МГД-уравнений с включённой аддитивной поправкой 7("Ь) интен-
сивности турбулентного эффекта
3v 1
— = PmAv - VP - E-1 Pm(ez xv) + fout + ( V x В) x В,
ot
ЯВ
— = Rem[V x (v x В)] + (Ra - Z(t))[V x a(r, 0)B)] + ДВ,
9t (1)
W = 0, W
VB = 0,
< rH B ? ) -.
где безразмерные величины: v - амплитуда поля скорости, В - амплитуда индукции магнитного поля, fout - массовая плотность внешних сил, P - давление, Pm = v/vm - магнитное число Прандтля, E = v/(OL2) - число Экмана, a(r,0) = aoa(r)cos0 - скалярная параметризация a-эффекта, ao - интенсивность a-эффекта, управляющие параметры: Rem = uoL/vm - магнитное число Рейнольд-са, которое отвечает за воздействие крупномасштабного генератора, Ra = aoL/vm -амплитуда a-эффекта, в которой заложена информация о турбулентном генераторе, выражены через величины: L [м] - радиус внешнего ядра, v [м2/с] - кинематическая вязкость, po = 7 ■ 103 кг/м3 - плотность, vm [м2/с] - магнитная вязкость, ri и г2 - радиусы-векторы внутренней и внешней границ сферической оболочки жидкого ядра, uo [м/с] - характерная величина скорости, П. [1/с] - угловая скорость вращения, единица времени L2/vm [с] - время диссипации магнитного поля с периодом 104 лет.
Процесс Z(t) регулирует степень влияния турбулентного эффекта, силу воздействия процесса определяет квадратичная форма Q(B(т),В(т)) - энергия магнитного поля, а характер воздействия - безразмерное ядро J(t-т) - знакопеременная функция, задающая затухающие колебания с коэффициентом затухания b = 0.1 и частотой a:
Z(t) =
J(t - t)Q(B (t),B (T))dT = [e-b(t-T)cos(a(t - т))]в 2(T)dT. (2)
Воздействие энергии поля максимально в момент времени а влияние в предыдущие моменты времени постепенно уменьшается с удалением от 1 к начальному моменту времени (нулю).
Для решения МГД-системы (1) используем маломодовое приближение галёр-кинского типа [3, 4, 15]. Известно, что источником для тороидальной (поло-идальной) компоненты магнитного поля является нелинейное взаимодействие по-лоидальной (тороидальной) компоненты этого поля с жидким потоком, поэтому для численного решения используются четыре моды - одна гидродинамическая и три магнитных: одна тороидальная и две полоидальных [4]:
V = иЩ У0 = и^На! Ут,1)0 + + а3^0,3,0 + + > (3)
В = ВО (1)в0До (г) + ВР «В ^(г) + В3Р «В ^(г),
где уо - мода Пуанкаре [3], |уо! = 1, и(Ъ) - амплитуда скорости, Вр1 о (г) дипольная составляющая магнитного поля, которая под воздействием дифференциального вращения порождает тороидальную В0 2 о (г) и полоидальную Вр 3 0(г) составляющие. Компоненты поля скорости и магнитного поля в данной модели независимы.
Подставим в систему (1) разложения (3) и функционал (2) и применим метод Галёркина. Результатом этих преобразований является система, которая используется для численной реализации с помощью неявного метода Эйлера,
9и ,2
_=_рти(1;)£ аЛк + ^иъ + Ц агЦкВ^Вк , ог к у ,к
9В •
^ = Кгти(Ъ^ о^кВк - Мг + (И« - г) £ ШОкВк , (4)
з,к к
ддг = ГВк-Ы-аг8, дг = аг-ъг8, м,к = 1,2,3. к
где Лг - собственное значение моды Пуанкаре [16], ц-г - коэффициент вязкой диссипации [16], коэффициенты Ь^к, М^к, М- - объемные интегралы от рассматриваемых полей [17], значения коэффициентов а вычислялись из задачи Пуанкаре о свободных малых колебаниях вязкой вращающейся жидкости [3].
Принятое значение шага Н на два порядка меньше наименьшего из значений (1/Лг) и (1/Цг), т.е. времён затухания полей скорости и магнитного соответственно. Выбор начальных условий обеспечивает выход системы (4) из точки покоя в начальный момент времени Ъ = 0 и дальнейшую взаимную генерацию компонент магнитного поля и поля скорости
и(о) = 1, В2(0) = о, Вр(0) = 1, В3р(0) = о, г(0) = о, г*(о) = о.
Затухающий эффект в поле скорости компенсируется воздействием внешних сил, поэтому в численной реализации значение массовой плотности внешних сил принято равным коэффициенту при и(Ъ) [6]. Численная модель рассматривается для управляющих параметров, изменяющихся в диапазонах Иет € (0, 1000] и € (0, 100], с шагом равным пяти.
Характеристики функции ядра
Регулирование воздействия турбулентного генератора в рассматриваемой модели задаётся процессом г(Ъ) с конечным временем ожидания Ъо («памяти»), которое определяется выбором знакопеременной затухающей функции в качестве ядра |(Ъ) (рис.1). При постоянной интенсивности а-эффекта в результате действия турбулентного генератора были получены три режима генерации магнитного поля в зависимости от параметров, определяющих влияние крупномасштабного и турбулентного генераторов. Наибольшая из трёх областей - область неограниченного возрастания магнитного поля (рис.2а), поэтому регулирование воздействия турбулентного генератора сводится к подавлению интенсивности а-эффекта.
Выбор процесса г(Ъ) в виде (2) определяет его суммарное воздействие на промежутке времени [0, Ъ], основной вклад в которое вносится в течение времени
0.1 t _ 0.1 t _ 0.1 t _ 0.1 t e e cos(0.1 t) e cos(0.5 t) e cos(t)
li01
llllta
(ffpw
1"................................!t
0 10 20 30 40 50 60 1
o.i t
e cos (5 f)
o.i t
e cos(lOt)
Рис. 1. Графики функции ядра J(t) = e0-1tcosat функционала Z(t). [Figure 1. The graphs of the kernel function J(t) = e0-ltcosat of the functional Z(t).]
ожидания to. Воздействие процесса максимально в момент времени t и постепенно ослабевает при удалении по временной оси от этой точки к нулю. Так как функция ядра J(t) - знакопеременная, то эффект от воздействия - подавление или подкачка - определяется знаком полученной интегральной характеристики процесса Z(t).
Для вычисления основной характеристики процесса, времени ожидания to, используется нормированная функция J(t) от модуля ядра J(t): J(t) = k|J(t)|, где k -нормирующий множитель, определялся из условия нормировки
J(t)dt =
k|J(t)|dt = k
-bt
|co s(at)|dt = 1.
Тогда математическое ожидание случайной величины T, заданной знакоположительной функцией J(t), это время ожидания to:
to = M(T) =
tJ(t) dt =
tke—bt|co s (at)|dt = k
te—bt|co s(at)|dt.
(5)
e
—oo
—oo
Из (5) время ожидания Ъо можно выразить следующим образом с помощью интеграла
t _ t> _
t0 _ ю _
te-bt|cos(at)|dt. (6)
Несобственный интеграл в правой части равенства (6) вычисляется аналитически. Раскрытие модуля преобразует исходный интеграл в сумму трёх интегралов
^ я/2а 5я/2а+2як/а
]Чв-м|соз(аъ) = $ ъв-мсоз(аъ) аъ + $ ъв-мсоз(аъ) аъ-
0 0 3я/2а+2як/а
3я/2а+2як/а ( )
- $ Ъе-ъъсс>э(аЪ)аъ, к = 0,1,2,...
я/2а+2як/а
К каждому из интегралов в правой части равенства (7) применяется интегрирование по частям
я/2а
bt , (b2 - а2) (п(а2 + b2)+4ab)e-nb/2a
te соs(at)dt +1-iiaW-, (8)
0
5я/2а+2як
J te-btcos(at)dt _ ((3п(а2 + b2) + 4ab + 4п(а2 + b2)k)e-3nb/2a+
3я/2а+2як (9)
л e-2nbk/a V '
+ (5п( а2 + b2) + 4ab + 4п( а2 + b2) к) e-5nb/2a) -—,,
v у /2(а2 + b2)2
3я/2а+2як
J te-bt соs(at)dt _ - ((3я(а2 + b2) + 4ab + 4п(a2 + b2)k) e-3nb/2a+
я/2а+2як (10)
+ (п( а2 + b2) + 4ab + 4п( а2 + b2) k) e-nb/2a) ■
e -2nbk/a
7 2(а2 + Ъ2)2'
Подставляем решения (8)-(10), зависящие от параметров а и Ъ функционала г(Ъ), в выражение (7) для вычисления несобственного интеграла
1 /
_ _^__(2(Ъ2- а2^
а 1 , £ te-bt|cos(at)|dt _ 2(а2 +b2)2 (2(b2- а2) +
+(п( а2 + b2) + 4ab) e-nb/2a + + (2п( а2 + b2) (e-3nb/2a + 2e-5nb/2a) +
а2 + Ъ2)+ 4аЬ)е пЪ/2а + + | 2П(а2
^ 7 (11) +(п(а2 + Ъ2) + 4аЪ) (е-пЪ/2а + 2е-3пЪ/2а + г-5пЪ/2а)) £ е-2пЪк/а+
' к=0
ОО \
+4п( а2 + Ъ2) (е-пЪ/2а + 2е-3пЪ/2а + е-5пЪ/2а) £ ке-2пЪк/а).
к=0 '
Решение (11) содержит сходящиеся ряды. Первый ряд является бесконечно убывающей геометрической прогрессией с первым членом равным единице и частным
q = e 2nb/a, а его сумма имеет вид
k=0
-2nbk/a_
1
1 — e-2nb/a'
(12)
Второй ряд может быть записан в следующей форме
L ke—
2nbk/a _
LL
e
2nbp/a
к=0 к=1 р=к
и сведён к вычислению суммы бесконечно убывающих геометрических прогрессий с первым членом равным е-2пЬр/а и частным q = е-2пЬ/а
L ke
k=0
—2nbk/a _
LL
k=1 p=k
e
2nbp/a _
-2nb/a
(1 — e—2nb/a)2'
(13)
Замена рядов на их суммы (12) и (13) в интеграле (11) даёт аналитическое решение для вычисления времени ожидания to, зависящее от параметров a и b ядра J(t),
J 1
to = J" te—bt|co s(at)|dt = 2 2 2 0 2(a2 + b2)2
2(b2 — a2)+
+(n(a2 + b2) + 4ab)e—nb/2a + (2n(a2 + b2) (e—3nb/2a + 2e—5nb/2a) + +(n( a2 + b2) + 4ab) (e—nb/2a + 2e—3nb/2a + e—5nb/2a))
+4n( a2 + b2) (e—nb/2a + 2e—3nb/2a + e—5nb/2a)
1
7 1 — e—2nb/a
e—2nb/a
(14)
(1 - e-2nb/a)2y
В Таблице 1 представлены вычисленные значения времени ожидания to с использованием решения (14) при выбранных значениях частоты затухания a и постоянном значении коэффициента затухания b = o.1, а также интегральные характеристики функции J(t) и её модуля.
Таблица 1
Характеристики функции ядра J(t) = e-o1^ras(at) [Characteristics of the kernel
function J(t) = e-o1tcos(at)]
a to J°|J(t)|dt o J J|J(t)|dt 0 J J(t)dt 0 J jJ(t)dt 0
1 2 3 4 5 6
0 100 9.999 10 9.999 10
0.1 57.5 7.147 7.173 4.978 5
0.5 56.6 6.385 6.406 0.3855 0.3846
1 59.4 6.359 6.376 0.1 0.099
5 62.7 6.3552 6.367 0.0038 0.0039
10 63.2 6.3548 6.366 0.00096 0.00099
e
Сравнение значений интегралов от функции |(Ь) и её модуля показывает, что суммарное воздействие |(Ь) (таб.1, пятый столбец) в течение времени ожидания Ьо меньше суммарного воздействия ||(Ь)| (таб.1, третий столбец) и при увеличении значений частоты затухания а это различие увеличивается на несколько порядков. Следовательно, несмотря на увеличение времени ожидания Ьо, при увеличении значений частоты а, эффект от воздействия процесса на турбулентный генератор будет уменьшаться.
Результаты моделирования и обсуждение
Исследуем влияние частоты затухающих колебаний а функционала подавления Z(t) на характер генерации магнитного поля при постоянном коэффициенте затухания b = 0.1. Результаты численного моделирования при заданных параметрах модели представлены на фазовой плоскости управляющих параметров Rem, Ra (рис.2). На рисунке 2a изображены области трёх режимов генерации магнитного поля (затухающий, регулярный, бесконечное возрастание магнитного поля), полученные при постоянной интенсивности а-эффекта [6]. В сравнении рассмотрим остальные результаты моделирования, когда введено подавление турбулентного генератора процессом Z(t). Отметим, что использование в качестве ядра функций J(t) = e0-1t [6] и J(t) = e^^cosO.lt привело к уменьшению области неограниченного роста (белая область) магнитного поля и увеличению количества режимов генерации магнитного поля (рис.2 б, в). Схожие результаты (есть небольшое отличие конфигураций областей с режимами васцилляции (жёлтая) и динамо-всплеска (фиолетовая)) определяются близкими значениями числовых характеристик (таб.1) выбранных функций ядра (рис.1), а влияние «хвостов» -незначительно.
Таблица 2
Наличие осцилляций в поле скорости [Presence of oscillations in the
velocity field]
Режимы генерации магнитного поля а
0.1 0.5 1 5 10
затухающий н н н н н
стационарный - н н н б
динамо-всплеск н м м б б
васцилляция н м м б б
м б б
регулярные н м (б: [300,470]) (м: [180,300]) (м: [165,210]
[350,370])
хаотические - б б б б
Обозначения: н - нет; м - есть малой амплитуды (менее 0.05), б - есть большой амплитуды.
В скобках указаны интервалы по оси Rem. [Notation: n - no; m - there is a small amplitude (less than 0.05), b - there is a large amplitude. The intervals along the axis Rem are indicated in brackets]
70 60 50 40 30
lot
R
70
д) 60; I ( е)
_ _ m П 1ПП ООП ^ПП 4ПП ^ПП £ПП 7ПП ЙПП ООП 1П|
R
0 100 200 300 400 500 600 700 800 900 1000
Ren,
0 100 200 300 400 500 600 700 800 900 1000 "
Рис. 2. Режимы генерации магнитного поля на плоскости управляющих параметров Rem и Ra: белая область - неограниченно возрастающий, красная - затухающий, зелёная - регулярный, синяя - стационарный, жёлтая - васцил-ляция, фиолетовая - динамо-всплеск, серая -хаотический. Интенсивность а-эффекта: а) постоянна ао [6] или определяется функцией Z(t) с ядром J(t): б) e_0,1t [6], в) e"0-ucos0.1t, г) e"0Jtcos0.5t, д) e-0-1tcost, е) e"0Jtcos5t, ж) e"0Jtcos10t.
[Figure 2. The modes of magnetic field generation on plane of the parameters Rem and Ra: the white region is the increases infinitely magnetic field, the red one is the damped, the green one is the steady, the blue one is the steady-state, the yellow one is the vacillation, the lilac one is the dynamo-burst, the grey one is the chaotic. The а-effect intensity is: a) the а0 [6] or defined by the function Z(t) with the kernel J(t): б) e-0^ [6], в) e-°.1tcos0.1t, г) e-°.1tcos0.5t, д) e"0.ucost, е) e"0.ucos5t, ж) e"0.ucos 10t.]
При увеличении частоты затухающих колебаний (а > 0.5) области всех режимов генерации магнитного поля на плоскости управляющих параметров модели сокращаются и уменьшается разнообразие генерируемых режимов. На границе с областью неограниченного роста магнитного поля появляются области хаотического режима, всегда сопровождающиеся осцилляциями в поле скорости. Начиная со значения а = 1 при генерации регулярного режима магнитного поля, появляются области с сильными осцилляциями в поле скорости (таб.2).
Таким образом, можно сделать вывод, что включение в МГД-систему процесса Z(t), регулирующего интенсивность a-эффекта, приводит к сокращению области неограниченного роста магнитного поля при 0<а<1 (рис.2в, г). В этих границах изменения значений частоты а затухающих колебаний получено большее разнообразие режимов генерации магнитного поля на фоне слабо изменяющегося поля скорости по сравнению со случаем постоянной интенсивности a-эффекта (рис.2 а) [6]. Если сравнивать с результатами для функции ядра J(t) = e0-1t (рис.2 б) [6], то при а = 0.1 наблюдается увеличение областей регулярного режима и васцилляции, размеры которых сокращаются при увеличении а.
На основании полученных результатов численного моделирования и анализа характеристик процесса Z(t) применимость модели при коэффициенте затухания b = 0.1 ограничена значениями частоты а < 1, а также значениями управляющих параметров модели Rem, Ra, из которых должны быть исключены промежутки с большими амплитудами осцилляций в поле скорости.
Заключение
Использование в крупномасштабной модели aO-динамо в качестве регулятора интенсивности a-эффекта процесса Z(t), который задаётся функционалом со знакопеременным ядром J(t) = e^^cosat, позволило получить такие виды генерации магнитного поля, как стационарный, регулярный, васцилляция, динамо-всплеск, на фоне слабо изменяющегося поля скорости для значений частоты затухающих осцилляций а < 1 при условии исключения определённых значений управляющих параметров Rem, Ra.
Область генерации бесконечно возрастающего магнитного поля уменьшилась в сравнении с картиной на фазовой плоскости, полученной при постоянной интенсивности a-эффекта [6]. При увеличении значений частоты a до значения единица области ограниченных незатухающих режимов генерации магнитного поля постепенно уменьшаются по сравнению с картиной на фазовой плоскости для функции ядра J(t) = ea1t [6].
Увеличение значений частоты a приводит к увеличению времени ожидания, но уменьшению суммарного накопленного эффекта и, как следствие, ослаблению подавляющего воздействия процесса Z(t) на турбулентный генератор.
Конкурирующие интересы. Конфликтов интересов в отношении авторства и публикации нет.
Авторский вклад и ответственность. Автор несёт полную ответственность за предоставление окончательной версии статьи в печатном виде. Окончательный вариант рукописи был одобрен автором.
Список литературы
1. Водинчар Г. М., Паровик Р. И., Пережогин А. С., Шереметьева О. В. Работы по моделированию физических процессов и систем в институте космофизических исследований и распространения радиоволн ДВО РАН, История науки и техники, 2017. №8, С. 100-112.
2. Водинчар Г.М. Использование собственных мод колебаний вязкой вращающейся жидкости в задаче крупномасштабного динамо, Вестник КРАУНЦ. Физ.-мат. науки, 2013. №2(7), С. 33-42 DOI: 10.18454/2079-6641-2013-7-2-33-42.
3. Feschenko L. K., Vodinchar G. M. Reversals in the large-scale aQ-dynamo with memory, Nonlinear Processes in Geophysics, 2015. vol.22, no. 4, pp. 361-369 DOI: 10.5194/npg-22-361-2015.
4. Vodinchar G. M., Feshchenko L. K. Model of geodynamo dryven by six-jet convection in the Earth's core, Magnetohydrodynamics, 2016. vol.52, no. 1, pp. 287-300.
5. Godomskaya A. N., Sheremetyeva O. V. Reversals in the low-mode model dynamo with aQ-generators, E3S Web of Conferences, 2018. vol.62, pp. 02016 DOI: 10.1051/ e3sconf/ 20186202016.
6. Шереметьева О. В., Годомская А. Н. Моделирование режимов генерации магнитного поля в маломодовой модели aQ-динамо с изменяющейся интенсивностью a-эффекта, Вестн. ЮУр-ГУ. Сер. Матем. моделирование и программирование, 2021. Т. 14, №2, С. 27-38 DOI: 10.14529/mmp210203.
7. Шереметьева О. В. Режимы генерации магнитного поля в маломодовой модели aQ-динамо с динамическим подавлением a-эффекта энергией поля, Вестник КРАУНЦ. Физ.-мат. науки, 2021. №4(37), С. 92-103 DOI: 10.26117/2079-6641-2021-37-4-92-103.
8. Godomskaya A. N., Sheremetyeva O. V.The modes of magnetic field generation in a low-mode model of aQ-dynamo with a-generator varying intensity regulated by a function with an alternating kernel, EPJ Web of Conferences, 2021. vol.254, pp. 02015 DOI: 10.1051/epjconf/202125402015.
9. Parker, E.N. Hydromagnetic dynamo models, Astrophys. J., 1955. no. 122, pp. 293-314.
10. Steenbek M., Krause F. Zur Dynamotheorie stellarer and planetarer Magnetfelder. I. Berechnunug sonnenahnlicher Wechselfeldgeneratoren, Astron. Nachr., 1969. no. 291, pp. 49-84.
11. Zeldovich Ya. B., Rusmaikin A. A., Sokoloff D. D. Magnetic fields in astrophysics. The Fluid Mechanics of Astrophysics and Geophysics. New York: Gordon and Breach, 1983. 382 pp.
12. Krause F., Radler K.H. Mean-filed magnetohydrodynamics and dynamo theory. Oxford: Pergamon Press, 1980. 271 pp.
13. Соколов Д. Д., Нефедов С. Н. Маломодовое приближение в задаче звездного динамо, Выч. мет. программирование, 2007. Т. 8, №2, С. 195-204.
14. Merril R. T., McElhinny M. W., McFadden P. L. The Magnetic Field of the Earth: Paleomag-netism, the Core, and the Deep Mantle. London: Academic Press, 1996. 531 pp.
15. Гледзер Е. Б., Должанский Ф. В., Обухов А. М. Системы гидродинамического типа и их применение. М.: Наука, 1981.368 с.
16. Водинчар Г. М. База данных «Параметры собственных мод свободных колебаний МГД полей в ядре Земли». гос. рег. 2019620054, 2019.
17. Водинчар Г. М., Фещенко Л.К. Библиотека программ для исследования «Маломодовой модели геодинамо» «LowModedGeodinamoModel». гос. рег. 50201100092, 2011.
Шереметьева Ольга ВладимировнаА - кандидат технических наук, научный сотрудник лаборатории моделирования физических процессов, Паратунка, Елизовский район, Камчатский край, Россия, (ЖСГО 0000-0001-9417-9731.
Sheremetyeva Olga Vladimirovna& - PhD (Tech.), Research Scientist, Laboratory of Physical Process Modeling, Institute of Cosmophysical Research and Radio Wave Propagation FEB RAS, Paratunka, Kamchatka, Russia, ORCID 0000-0001-9417-9731.
Vestnik KRAUNC. Fiz.-Mat. Nauki. 2022. vol. 41. no. 4. pp. 107-119. ISSN 2079-6641
MSC 65N80 Research Article
Dynamics of generation modes changes in magnetic field depending on the oscillation frequency of the a-effect suppression process by field energy in the aO-dynamo model
O. V. Sheremetyeva
Institute of Cosmophysical Research and Radio Wave Propagation, FEB RAS, 684034 Kamchatka region, Elizovskiy district, Paratunka, Mirnaya str., 7, Russia E-mail: sheremeteva@ikir.ru
This article considers a aO-dynamo large-scale model in a low-mode approximation. The intensity of the a-effect is regulated by a process with hereditarity properties, depending on the energy of the magnetic field. The nature of the process impact is determined by an alternating kernel function with a variable damping frequency and a constant damping coefficient equal to 0.1. Based on the results of numerical simulation of the magnetic field generation modes, the boundaries for the values of the model parameters are determined for which there are no oscillations in the velocity field of a viscous liquid or they are insignificant. The results of magnetic field generation modes simulation at various damping frequency values are displayed on the phase plane of the control parameters. These parameters contain information about large-scale and turbulent generators. The paper investigates the question about the dynamics of changes in the pattern on the phase plane depending on the value of the damping frequency and the waiting time. A comparison is made with the results obtained earlier for the constant intensity of the a-effect and for the intensity, the change of which is determined by the process with exponential kernel.
Key words: low-mode dynamo model, aO-dynamo, a-effect, magnetic field, reversals.
d DOI: 10.26117/2079-6641-2022-41-4-107-119
Original article submitted: 29.11.2022 Revision submitted: 06.12.2022
For citation. Sheremetyeva O. V. Dynamics of generation modes changes in magnetic field depending on the oscillation frequency of the a-effect suppression process by field energy in the aO-dynamo model. Vestnik KRAUNC. Fiz.-mat. nauki. 2022,41: 4,107-119. d DOI: 10.26117/2079-6641-2022-41-4-107-119
Competing interests. The author declare that there is no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. The author is solely responsible for providing the final version of the article in print. The final version of the manuscript was approved by the author. The content is published under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/deed.ru)
© Sheremetyeva O.V., 2022
Funding. The work was carried out within the framework of the state assignment on the topic «Physical processes in the system of near space and geospheres under solar and lithospheric influences» (No. AAAA-A21-121011290003-0).
References
[1
[2
[3 [4
[5]
[6]
[7 [8
[9 10
11
12
13
14
15
16 17
Vodinchar G. M., Parovik R. I., Perezhogin A. S., Sheremetyeva O. V. Simulation of physical processes and systems in institute of cosmophysical research and radio wave propagation FEB RAS [Raboty po modelirovaniyu fizicheskih processov i sistem v institute kosmofizicheskih issledovanij i rasprostraneniya radiovoln DVO RAN], History of Science and Engineering [Istoriya nauki i tehniki], 2017, 8, 100-112. (In Russian) Vodinchar G. M. Using Modes of Free Oscillation of a Rotating Viscous Fluid in the Large-Scale Dinamo, Bulletin KRASEC. Phys. & Math. Sci., 2013, 2(7), 33-42. (In Russian) DOI: 10.18454/2079-6641-2013-7-2-33-42
Feschenko L. K., Vodinchar G. M. Reversals in the large-scale aQ-dynamo with memory, Nonlinear Processes in Geophysics, 2015, 22: 4, 361-369 DOI: 10.5194/npg-22-361-2015 Vodinchar G. M., Feshchenko L. K. Model of geodynamo dryven by six-jet convection in the Earth's core, Magnetohydrodynamics, 2016, 52:1, 287-300.
Godomskaya A. N., Sheremetyeva O. V. Reversals in the low-mode model dynamo with aQ-generators, E3S Web of Conferences, 2018, 62, 02016. DOI: 10.1051/ e3sconf/ 20186202016 Sheremetyeva O. V., Godomskaya A. N. Modelling the magnetic field generation modes in the low-mode model of the aQ-dynamo with varying intensity of the a-effect, Bulletin of the South Ural State University, Series «Mathematical Modelling, Programming & Computer Software», 2021, vol. 14, no. 2, pp. 27-38 (In Russian DOI: 10.14529/mmp210203). Sheremetyeva O.V. Modes of magnetic field generation in the low-mode aQ-dynamo model with dynamic regulation of the a-effect by the field energy. Vestnik KRAUNC. Fiz.-mat. nauki. 2021, 37:4, 92-103 (In Russian DOI: 10.26117/2079-6641-2021-37-4-92-103). Godomskaya A. N., Sheremetyeva O. V. The modes of magnetic field generation in a low-mode model of aQ-dynamo with a-generator varying intensity regulated by a function with an alternating kernel, EPJ Web of Conferences, 2021, 254, 02015. DOI: 10.1051/epjconf/202125402015
Parker, E.N. Hydromagnetic dynamo models, Astrophys. J., 1955, 122, 293-314. Steenbek M., Krause F. Zur Dynamotheorie stellarer und planetarer Magnetfelder I. Berechnung sonnenahnlicher Wechselfeldgeneratoren, Astron.Nachr., 1969, 291, 49-84. Zeldovich Ya. B., Rusmaikin A. A., Sokoloff D. D. Magnetic fields in astrophysics. The Fluid Mechanics of Astrophysic and Geophysics. Gordon and Breach: New York, 1983, 382. Krause F., Raadler K.H. Mean-filed magnetohydrodynamics and dynamo theory. Pergamon Press: Oxford, 1980, 271.
Sokoloff D. D., Nefedov S. N. A small-mode approximation in the stellar dynamo problem, Num. Meth. Prog., 2007, 8:2, 195-204 (In Russian http://mi.mathnet.ru/vmp485). Merril R. T., McElhinny M. W., McFadden P. L. The Magnetic Field of the Earth: Paleomagnetism, the Core, and the Deep Mantle. Academic Press: London, 1996, 531. Gledzer E. B., Dolzhanskiy F. V., Obukhov A. M. Sistemy gidrodinamicheskogo tipa i ikh primenenie [Hydrodynamic Type Systems and Their Application]. Nauka: Moscow, 1981, 368 (In Russian).
Vodinchar G. M. Baza dannyh «Parametry sobstvennyh mod svobodnyh kolebanij MGD polej v yadre Zemli» [Database «Parameters of natural modes of normal oscillations of MHD fields in the Earth's core»], 2019620054, 2019 (In Russian).
Vodinchar G. M., Feshenko L. K. Biblioteka programm dlya issledovaniya «Malomodovoj modeli geodinamo» [Library of programs for researching «Low-mode geodynamo model» «LowModedGeodinamoModel»], 50201100092, 2011 (In Russian).