УДК 517.958:537.84
БСТ: 10.14529/ mmp210203
МОДЕЛИРОВАНИЕ РЕЖИМОВ ГЕНЕРАЦИИ МАГНИТНОГО ПОЛЯ В МАЛОМОДОВОЙ МОДЕЛИ аП-ДИНАМО С ИЗМЕНЯЮЩЕЙСЯ ИНТЕНСИВНОСТЬЮ а-ЭФФЕКТА
О.В. Шереметьева1, А.Н. Годомская2
1 Институт космофизических исследований и распространения радиоволн ДВО РАН, с. Паратунка, Российская Федерация
2Муниципальное бюджетное учреждение дополнительного образования «Центр «Луч>, г. Елизово, Российская Федерация
В рамках модели аО-динамо рассматривается магнитогидродинамическая система (далее МГД-система) с введенной аддитивной поправкой интенсивности а-эффекта в виде функции Z(£). Изменение интенсивности а-эффекта со временем определяется показательным ядром J(£) функции Z(£). Проведен предварительный анализ влияния изменения значений массовой плотности внешних сил и интенсивности а-эффекта на значения магнитного поля и поля скорости. Для упрощения численной модели МГД-система перемасштабирована и в качестве единицы времени принято время диссипации магнитного поля, которое составило порядка 104 лет. Выбор временной единицы определяется необходимостью исследования поведения магнитного поля на больших временных масштабах. С целью сокращения количества варьируемых параметров МГД-система обезразмерена таким образом, что управляющими параметрами выступают число Рейнольдса, которое несет информацию о крупномасштабном генераторе, и амплитуда а-эффекта, характеризующая турбулентный генератор. Результаты численного моделирования режимов генерации магнитного поля отражены на фазовой плоскости управляющих параметров, и исследуется вопрос о динамике изменения картины на фазовой плоскости в зависимости от времени ожидания, определяемого масштабным коэффициентом показательного ядра J(£) функции Z(£).
Ключевые слова: аО-динамо; магнитогидродинамическая система; 'режимы генерации магнитного поля; инверсия.
Введение
Магнитное поле является одним из важных аспектов наличия жизни на планете. Исследование природы возникновения магнитного поля и моделирование вариаций поля является интенсивно разрабатываемым разделом теории динамо. Изучение проводится как с использованием прямого численного моделирования, так и на основе упрощенных моделей. Магнитные поля реальных космических объектов (планет, звезд, галактик) с течением времени претерпевают изменения, которые носят как хаотический, так и периодический характер. Известным свойством магнитного поля является резкая смена полярности дипольной составляющей поля - инверсия. Большой интерес представляют инверсии без существенной перестройки движения проводящей среды. Протяженность интервалов времени между инверсиями (интервалы полярности) может отличаться на несколько порядков, существуют длительные интервалы без инверсий [1]. Информацию об инверсиях магнитного поля Земли получают из па-леомагнитных записей, на основе которых строится шкала геомагнитной полярности. Построение модели динамо, полностью воспроизводящей реальную палеомагнитную
шкалу, является трудновыполнимой задачей. Прямое численное моделирование маг-нитогидродинамических уравнений позволяет воспроизвести многочисленные режимы инверсий, но не отвечает на вопрос об их причине. Кроме этого, полные уравнения содержат очень много параметров, оценки которых часто расходятся на много порядков, либо вообще не существуют. Поэтому исследователи с помощью простых динамических систем малой размерности пытаются объяснить физическую причину, признаки, важнейшие свойства этого явления [1-5].
Для систем космического динамо (планет, звезд, галактик) обычно рассматривают три механизма аП, а2 и а2П, из которых третий является наиболее общим [3,6,7]. Различные модели динамо позволяют получить случайные последовательности инверсий, свойства которых сильно отличаются в зависимости от выбора параметров. Первые два используются, когда эффективности а- и П-генераторов резко отличаются. В частности, для Земли характерны сильное дифференциальное вращение объекта и турбулентный характер движения проводящей среды, при которых генерацию магнитного поля описывают с помощью аП-динамо [2,3,6-8]. В своей работе мы используем маломодовую модель аП-динамо [9] для исследования смены полярности магнитного поля, не связанной с перестройкой конвекции. В рамках данной модели такое исследование удобнее проводить во временных масштабах, сравнимых со временем затухания магнитного поля. Использование новой временной шкалы приведет к перемасштабированию магнитогидродинамической системы, описывающей динамо. Кроме того, целесообразным представляется обезразмеривание системы в целом, что снизит возможность возникновения расчетных погрешностей на длительных временных интервалах, перегрузки численных решений большими значениями и ухода на «машинную» бесконечность. Данная статья посвящена решению задач перемаштаби-рования и обезразмеривания МГД-системы в рамках принятых ограничений модели аП-динамо, численному моделированию МГД-системы с учетом введенной аддитивной поправки интенсивности а-эффекта в виде интегральной функции с показательным ядром и анализу результатов моделирования.
1. Постановка задачи
Опишем основные положения и ограничения модели. Исследование проводится в рамках модели аП-динамо, где жидкое ядро Земли полагается ограниченным сферической оболочкой с центром в начале координат. Сферическая оболочка вращается вокруг оси Oz с постоянной угловой скоростью Q. Поле скорости v вязкой несжимаемой жидкости внешнего ядра и порождаемое магнитное поле считаем аксиально симметричными. Внутренняя граница оболочки жидкого ядра имеет радиус r = r 1, а внешняя - радиус r = r2 (r1 < r2), на которых поле скорости вязкой жидкости v нулевое (граничные условия прилипания). Принимаем магнитную проницаемость внутреннего и внешнего ядер одинаковыми, а среду за границей внешнего ядра (r > r2) - непроводящей. Таким образом, для магнитного поля принимаются вакуумное граничное условие на внешней границе и условие ограниченности в центре Земли. Физические параметры оболочки считаем неизменными и турбулентность в ядре изотропной. Полагаем также, что турбулентный а-эффект антисимметричен относительно экватора, поэтому используем его скалярную параметризацию в виде функции а^,9) = аой^) cos 9 = ао^^,9), где положительный коэффициент ао
определяет интенсивность а-эффекта, а1 (г, в) - безразмерная величина, радиальная составляющая а(г) удовлетворяет условию тах | а (г) | ~ 1. В соответствии с принятой моделью динамо считаем, что среднее течение V носит характер дифференциального вращения в поле внешних сил с массовой плотностью ^оиг [2,3,7,8]. Тогда маломодо-вая модель будет описываться МГД-уравнениями, включающими уравнение Навье -Стокса, уравнение индукции магнитного поля В с учетом турбулентного а-эффекта, условие неразрывности поля скорости V, условие соленоидальности магнитного поля и граничные условия. В приближении Буссинеска МГД-система принимает следующий вид:
дv 1
— + (уУ)У + £ = !/АУ - -УР - Ъ + ^ + Ъ,
дь Р0
дВ
- = Ух(ухВ) + Ух(а(г)р) + г/тДВ,
= 0, УВ = 0,
v(гl) = v(г2) = 0,
где ускорение центробежной силы
Ъ = О х (О х г), (2)
значение угловой скорости О изменяется в пределах 10-1 ^ 10 1/с, V - кинематическая вязкость (пределы изменения 10-6 ^ 102 м2/с), р0 = 7 • 103 кг/м3 - плотность, Р -давление, vm - магнитная вязкость (изменяется в пределах 1 ^ 20 м2/с), массовая плотность силы Кориолиса
= 2О х v, (3)
- массовая плотность поля внешних сил (источник полоидальной скорости), массовая плотность силы Лоренца
= —-—(V х В) х В, (4)
Ро^о ^
= 4п • 10-7 Г/м - магнитная постоянная, ^ =1 - магнитная проницаемость, г1 и г2 - радиусы-векторы внутренней и внешней границ сферической оболочки жидкого ядра.
При постоянной интенсивности а-эффекта а0 в рамках используемой модели система (1) определяет три режима магнитного поля: затухающий, регулярный и неограниченное увеличение (расходимость поля) (рисунок а). Интенсивность а0 может изменяться со временем под действием внешних сил. Такое изменение определим путем ввода в динамическую систему (1) аддитивной поправки для интенсивности а-эффекта в виде функции
г
2(Ь) = ] ^ - т)3(В(т), В(т))^т. (5)
0
Функция (5) имеет размерность скорости [м/с], как и функция а(г, в), и используется с целью регулирования степени влияния турбулентного эффекта за счет задержки
О 100 200 300 400 500 600 700 800 900 1000 0 100 200 300 400 500 600 700 800 900 1000
0 100 200 300 400 500 600 700 800 900 1000 0 100 200 300 400 500 600 700 800 900 1000
Характер генерации магнитного поля в зависимости от параметров Явт и Яа. Интенсивность а-эффекта: а) постоянная - ао; в остальных случаях определяется функцией Z(г) с ядром: ь) з(г) = в-0'1£, с) з(г) = в-0'5£, а) з(г) = е-, е) з(г) = в-5£, о з(г) = в-10£. Белая область - магнитное поле бесконечно возрастает, разреженная наклонная штриховка - магнитное поле затухает, частая наклонная штриховка - регулярный режим генерации магнитного поля, вертикальная штриховка - стационарный режим, черная область -васцилляция, серая область - динамо-всплеск
или поддержания его воздействия. Характер такого воздействия определяется безразмерным ядром J(t — т). Квадратичная форма QВ(т), В(т)) в подынтегральном выражении имеет размерность [м/с2], а ее коэффициенты qj [м/(сТл)2]. С учетом аддитивной поправки для интенсивности а-эффекта МГД-система (1) запишется следующим образом
d v 1
— + (vV)v + fc = uAv - -VP -ÎK + îout + fL,
dt Po
dB
— = V x (v x В) + V x (a0 - Z{t)) a\r, 9)В + umAB,
Vv = 0, VB = 0,
v(ri) = v(r2) = 0.
2. Свойства МГД-системы
Прежде, чем переходить к обезразмериванию, а в дальнейшем к численному решению системы (6), преобразуем ее и исследуем влияние на результат численного моделирования изменения основных параметров системы уравнений в А раз.
Поскольку при постоянной угловой скорости Q ускорение центробежной силы (2) в (1) представимо в следующем виде
12
fc = Г2 х (Г2 х г) =--У(Г2 х г)2, (7)
то его можно объединить вместе с давлением в одно потенциальное слагаемое —VP*,
Po
где редуцированное давление будет задаваться выражением
12
Р* = Р + -р0(^ хг)2. (8)
Тогда система (6) в соответствии с принятыми обозначениями (3), (4), (8) и учетом того факта, что (vV)v = 0, т. к. в модели рассматриваются собственные колебания вязкой вращающейся жидкости малой амплитуды [9], может быть записана в следующем виде
dv 1 1
— = vAv--VP* -2ilxv + îout +-(V x В) x В,
dt po po VoV
dB
— = V x (v x B) + (a0 - Z{t)) (V x a'(r, 9) B) + umAB, Vdvt = 0,
VB = 0,
v(ri) = v(r2) = 0.
Исследуем, сохранится ли инвариантность формы записи МГД-системы (9) при изменении массовой плотности внешних сил в А раз
iout
А îout. (10)
Выразим из (10) подставим в систему (9) и, учитывая вид функции Z(¿) (5), получим систему (9) в преобразованном виде
дV 1 1 ~ — = 1/Ду--УР - 2П х V + ^ +-(V х В) х В,
д£ Ро Ро^о ^
® = ^х(ухВ) + Ц- (V х с/(г, В) + итАВ, (п)
Уу- = 0, УВ = 0,
У(Г1) = У(Г2) = 0.
Общий вид всех уравнений, кроме уравнения индукции, в результате перемасштабирования не изменится. Таким образом, увеличение или уменьшение воздействия внешних сил в Л раз приводит к изменению масштаба магнитного поля соответственно в \/Л раз, а полей скорости и давления в Л раз. При увеличении значений Л увеличиваются амплитуды этих полей, однако коэффициент 1/А при этом уменьшается, что приводит к уменьшению влияния функции Z (¿) на интенсивность а-эффекта, а значит к уменьшению возможности изменения режимов генерации магнитного поля на фоне незначительных изменений поля скорости по сравнению со случаем, когда интенсивность а-эффекта постоянна.
К аналогичным выводам приходим, если проводить исследование системы (9) при изменении поля скорости или магнитного поля в А раз. Неинвариантность формы уравнения магнитной индукции при таком преобразовании может повлиять на режим генерации магнитного поля и изменить картину на фазовой плоскости. Отметим также, что появление нормировочного коэффициента А в ядре функции Z(¿)
г
Z(í) = y, А7(* - т)д(В(т), В(т))^т (12)
о
равносильно перемасштабированию магнитного поля г г
гц) = J J(t — т)Q(\^B(т),VXB(т))dт = ! Л1-т)С}{Щт)}Щт))(1т. (13) оо
Поэтому произойдет такое же перемасштабирование остальных параметров, но в данном случае в левой части уравнения индукции появится коэффициент 1/А только перед первым слагаемым, что не окажет влияния на функцию Z(¿) и не изменит степени ее воздействия на амплитуду а-эффекта. Таким образом, введение нормировочного коэффициента в функцию Z (¿) приведет к перемасштабированию поля скорости и магнитного поля, а также повлияет на режимы генерации магнитного поля.
3. Обезразмеривание уравнений МГД-системы
Введем характерные величины скорости и0 [м/с] и линейного размера области Ь [м] (радиус внешнего ядра), а в качестве составных будем рассматривать единицы магнитной индукции В0 [Тл], давления Р0 [Па], массовой плотности внешних сил /о [Н/кг] и угловой скорости П0 [1/с]. Исходя из принятых предположений модели, определим поле скорости как линейную комбинацию V = и0 V', где V' - безразмерное поле, индукцию магнитного поля как В = В0 В', где В' - безразмерное поле, угловую
скорость Q = ez, где ez - единичный вектор, направленный вдоль оси вращения Oz, давление P* = P0 P', где P' - безразмерная величина, и массовую плотность внешних сил fout = f0 fout, где f'out - безразмерное поле. В качестве единицы времени используем время затухания магнитного поля при отсутствии внешних воздействий (подкачки) - время диссипации магнитного поля L2/vm [с], порядок которого составляет 104 лет. Тогда переход к новой временной шкале определяется соотношением L2
t = —t'. Проведем перемасштабирование системы (9), принимая во внимание вве-
Vm
денные характерные величины, ограничения и зависимости, и обезразмерим ее 7^77 = — Av--VP--x v ) +-W+
dt Vm P0 Uo Vm Vm Uo Vm
+-^-(VxB')xB',
^ ^0 Po Vm Uo
^ = ^[V x (v' x B')] + (— - —Z(t')\ [V x a'(r,0)B'] + AB', (14)
dt Vm ^ Vm Vm '
Vv' = 0,
VB' = 0,
v'g)=v'(e2) = 0.
В качестве управляющих парамтеров системы примем безразмерные параметры:
—---это магнитное число Рейнольдса Rem и Ra = —--безразмерная амплитуда
Vm V Vm V
ск-эффекта. Учитывая, что — = Рт - магнитное число Прандтля и-- = Е - число
Vm ^L2
Экмана, коэффициент в кориолисовом члене можно выразить следующим образом
Hob2 n0L2 и
-=--= Е Рт. Для упрощения системы и дальнейшего ее использования
Vm V Vm
для численного моделирования удобно принять равными единицам коэффициенты
P0L L2 f0 B2L
1, -= 1, -= 1, выбрав соответствующие значения магнитной
индукции В0, давления Р0 и массовой плотности внешних сил /0. В результате перехода к новой системе единиц в обезразмеренном виде МГД-система (9) запишется следующим образом
д V' 1 , ,
— = РтАлг' - УР' - Е~ Рт(ег х V) + С + (V х В') х В',
дБ' , , , ,
— = Дет[У х (у' х В')] + (Па - г'Ц'Ш х а'{г, 0)В')] + ДВ',
' = 0, ( ) УБ' = 0,
4. Численное моделирование
Для численной реализации системы (15) используем маломодовое приближение, включающее минимальное количество мод скорости и магнитной индукции, позволяющее получить инверсии в модели аП-динамо с переменной интенсивностью а-
генератора [10-12]. В данном случае ограничиваемся одномодовым приближением поля скорости
V = и(г)уо = ! , 0 + a2vp2 , 0 + «3^, 3 , 0 + ап ^^^, 0 + а^Р 3 , 0), (16)
и трехмодовым приближением магнитного поля
В = В^В^М + вр (¿)В0Р1;0(г) + вр (¿)В0Рз;0(г), (17)
где v0(г) - мода Пуанкаре, представленная в виде разложения в гильбертовом подпространстве НТ = ^т10, v0:'2O, vT30, ... } [4] таким образом, чтобы ^0(г)| = 1, и(£) -амплитуда скорости, Вр1 0(г) - дипольная составляющая магнитного поля, которая под воздействием дифференциального вращения порождает тороидальную ВТ2 0(г) и полоидальную В03 0(г) составляющие. Компоненты поля скорости и магнитного поля считаются независимыми. Функция Z(£) (5) задается как
Z (t) = e-b(t-T) В2(т )dr, (18)
где характер воздействия определяется показательным ядром 3(£) = е-64. Выбор в качестве ядра 3 (£) показательной функции определяет процесс с конечным временем ожидания или памяти, подавляющее влияние которого на изменение интенсивности а-эффекта включается в начальный момент времени и постепенно накапливается в течение времени ожидания. Для различных значений масштабного коэффициента Ь время ожидания £0 вычислялось, исходя из соотношения
4о те
[ е-6^ = 0.95 / е-64(19)
Подставим в систему (15) разложения (16), (17) и функцию (18). После применения метода Галеркина получаем систему
дм ~dt дВ.
г,, - -Pmu{t)Y^OL2k\k + fout + J] aiLijkBjBk, dt k i, j, k
Remu(t) £ aj WjkBk - ^Вг + (Ra - Z) £ WgBk, (20)
dt j , k
где Ai - собственное значение моды Пуанкаре, коэффициенты , Wijk, - объемные интегралы от рассматриваемых полей, - коэффициент вязкой диссипации.
Для численной реализации системы (20) использованы неявный метод Эйлера и явный четырехэтапный метод Рунге - Кутты. Шаг дискретизации к выбран равным 10-5. Приняты начальные условия, обеспечивающие выход системы из точки покоя в момент времени £ = 0 и генерацию остальных компонент магнитного поля,
и(о) = 1, вТ(0) = о, вр(0) = 1, вр(0) = 0, Z(0) = 0. (21)
t
0
0
0
Значение массовой плотности внешних сил /ои принято равным коэффициенту при м(£) для обеспечения компенсации затухающего эффекта в поле скорости. Управляющие параметры модели - магнитное число Рейнольдса Дет и амплитуда а-эф-фекта задавались в диапазонах (0, 1000] и (0, 100] соответственно. Вычисления проводились для значений масштабного коэффициента Ь, представленных в таблице. Оба метода при заданных параметрах модели дали идентичные результаты, которые представлены на фазовой плоскости параметров Дет, Ка (рисунок Ь - :£). Полученные результаты показывают, что с увеличением значения масштабного коэффициента Ь ядра функции Z (¿) на фазовой плоскости уменьшается область расходимости магнитного поля, по сравнению со случаем, когда интенсивность а-эффекта постоянна (рисунок а). Появляются области, где генерируются новые режимы: стационарный, васцилляция, динамо-всплеск. С ростом значений масштабного коэффициента Ь значительно увеличиваются области регулярного и стационарного режимов генерации магнитного поля, в то время как области васцилляций и динамо-всплеска постепенно уменьшаются. Такие изменения на фазовой плоскости связаны с поведением ядра функции Z(¿).
Время ожидания ¿о при увеличении значений масштабного коэффициента - уменьшается, как и интеграл от функции 3(¿) на промежутке [0, то) (см. таблицу). Таким образом, при одинаковых в остальном условиях численного эксперимента разнообразие режимов генерации магнитного поля в условиях слабо изменяющегося поля скорости увеличивается при уменьшении значений масштабного коэффициента Ь, когда а-эффект подавляется процессом, описываемым функцией Z(¿), в течение длительного времени ожидания. Воздействие растянуто во времени, что позволяет на интервале 0 < Дет < 500 подавить расходимость поля при Ка < 50. Когда значения масштабного коэффициента увеличиваются (Ь > 1), на первых отсчетах времени сосредоточивается большее по значению тормозящее воздействие процесса Z (¿) и незначительно изменяется в течение остального времени. Это воздействие уменьшает влияние а-эффекта и значительно увеличивает область незатухающих ограниченных осцилля-ций магнитного поля, но приводит к уменьшению разнообразия режимов генерации.
Таким образом, введение аддитивной поправки интенсивности а-эффекта с ядром в виде показательной функции приводит к увеличению области генерации незатухающих осцилляций магнитного поля, в том числе с инверсиями, на фоне слабо изменяющегося поля скорости и сокращению области расходимости магнитного поля. Дополнительное подавляющее а-эффект воздействие процесса, определяемого функцией Z(¿), естественным образом не приводит к уменьшению области затухания магнитного поля, а лишь незначительно увеличивает ее при Ь < 1 .
Заключение
В рамках принятой маломодовой модели аП-динамо и предположении аксиальной симметрии поля скорости и магнитного поля проведено перемасштабирование и обез-
Таблица
Характеристики функции 3 (¿)
ь ¿0 интеграл от
0,1 29,957 10
0,5 5,992 2
1 2,996 1
5 0,599 0,2
10 0,3 0,1
размеривание МГД-системы. Временной масштаб выбран равным времени диссипации магнитного поля под воздействием внешней силы и при отсутствии в дальнейшем дополнительных воздействий. Система приведена к виду, упрощающему дальнейшую численную реализацию математической модели.
Исследование свойств МГД-системы показало, что увеличение воздействия внешних сил приведет к увеличению амплитуд поля скорости и индукции магнитного поля. Причем рост значений поля скорости будет интенсивнее в \/Л раз, что может приводить к значительным изменениям его амплитуды. В то же время влияние функции Z(t) будет сводиться к минимуму и не приведет к значительным изменениям режимов генерации магнитного поля по сравнению со случаем постоянной интенсивности а-эффекта. Следовательно, при реализации численной схемы обезразмеренной МГД-системы (15) для обеспечения компенсации затухающего эффекта в поле скорости достаточно принять значение массовой плотности внешних сил равным коэффициенту при u(t) в уравнении Навье - Стокса, а для уменьшения вариаций амплитуды поля скорости ввести нормировочный коэффициент в ядро функции Z (t).
Результаты расчетов показали, что торможение интенсивности а-эффекта путем введения процесса, описываемого функцией Z(t) с показательным ядром, приводит по сравнению со случаем постоянной интенсивности ао к появлению новых режимов генерации магнитного поля, в том числе с инверсиями. Воздействие процесса, определяемое ядром с б о льшим временем ожидания, приводит к большему разнообразию режимов генерации магнитного поля, а с уменьшением времени ожидания -значительно уменьшается область расходимости магнитного поля.
Литература
1. Merril, R.T. The Magnetic Field of the Earth: Paleomagnetism, the Core, and the Deep Mantle / R.T. Merril, M.W. McElhinny, P.L. McFadden. - London: Academic Press, 1996.
2. Steenbek, M. Zur Dynamotheorie stellarer und planetarer magnetfelder. Berechnunug sonnenahnlicher wechselfeld generatoren / M. Steenbek, F. Krause // Astronomische Nachrichten. - 1969. - № 291. - P. 49-84.
3. Zeldovich, Ya.B. Magnetic Fields in Astrophysics. The Fluid Mechanics of Astrophysics and Geophysics / Ya.B. Zeldovich, A.A. Rusmaikin, D.D. Sokoloff. - New York: Gordon and Breach, 1983.
4. Водинчар, Г.М. Инверсии в модели геодинамо, управляемой 6-ячейковой конвекцией / Г.М. Водинчар, Л.К. Фещенко // Вестник КРАУНЦ. Физико-математические науки. -2015. - Т. 11, № 2. - С. 45-54.
5. Гледзер, Е.Б. Системы гидродинамического типа и их применение / Е.Б. Гледзер, Ф.В. Должанский, А.М. Обухов. - М.: Наука, 1981.
6. Parker, E.N. Hydromagnetic Dynamo Models / E.N. Parker // The Astrophysical Journal. -1955. - № 122. - P. 293-314.
7. Krause, F. Mean-Filed Magnetohydrodynamics and Dynamo Theory / F. Krause, K.-H. Raedler. - Oxford: Pergamon Press, 1980.
8. Steenbek, M. Berechnung der mittlerer Lorentz-Field starke о x B für ein elektrisch leitendes medium in Turbulenter, Durch Coriolis-Krafte beenflusster bewegung / M. Steenbek, F. Krause, K.-H. Raedler // Zeitschrift fur naturforschung. - 1996. - № 21. - P. 369-376.
9. Водинчар, Г.М. Использование собственных мод колебаний вязкой вращающейся жидкости в задаче крупномасштабного динамо / Г.М. Водинчар // Вестник КРАУНЦ. Физико-математические науки. - 2013. - Т. 7, № 2. - С. 33-42.
10. Godomskaya, A.N. Reversals in the Low-Mode Model Dynamo with aQ-generators / A.N. Godomskaya, O.V. Sheremetyeva // E3S Web of Conferences. - 2018. - V. 62. -Article ID: 02016. - 6 p.
11. Godomskaya, A.N. Modes of Magnetic Field Generation in Models of a aQ-dynamo with a Power Type a-generator / A.N. Godomskaya, O.V. Sheremetyeva // E3S Web of Conferences. - 2019. - V. 127. - Article ID: 02023. - 7 p.
12. Годомская, А.Н. Режимы генерации магнитного поля в модели aQ-динамо с a-генератором степенного типа / А.Н. Годомская, О.В. Шереметьева // Вестник КРАУНЦ. Физико-математические науки. - 2019. - Т. 29, № 4. - C. 58-66.
Ольга Владимировна Шереметьева, кандидат технических наук, научный сотрудник, лаборатория моделирования физических процессов, Институт космофизических исследований и распространения радиоволн ДВО РАН (Камчатский край, с. Пара-тунка, Российская Федерация), sheremeteva@ikir.ru.
Анна Николаевна Годомская, педагог дополнительного образования, естественнонаучный отдел, Муниципальное бюджетное учреждение дополнительного образования «Центр «Луч> (г. Елизово, Российская Федерация), anna_antonenko@mail.ru.
Поступила в редакцию 17 сентября 2020 г.
MSC 65N80 DOI: 10.14529/mmp210203
MODELLING THE MAGNETIC FIELD GENERATION MODES IN THE LOW-MODE MODEL OF THE aQ-DYNAMO WITH VARYING INTENSITY OF THE a-EFFECT
O.V. Sheremetyeva1, A.N. Godomskaya2
1 Institute of Cosmophysical Research and Radio Wave Propagation FEB RAS, Paratunka, Russian Federation
2Municipal Budgetary Institution of Supplementary Education "Center "Luch", Yelizovo, Russian Federation
E-mails: sheremeteva@ikir.ru, anna_antonenko@mail.ru
In the framework of the a^-dynamo model a magnetohydrodynamic system (MHD-system) with an included additive correction of the a-effect intensity in the form of a function Z(t) is considered. The variation of the a-effect intensity with time is determined by the exponential kernel J(t) of the function Z(t). A preliminary analysis of the changes influence in the values of the mass density of external forces and the a-effect intensity on the values of the magnetic field and the velocity field is carried out. To simplify the numerical model, the MHD-system was rescaled and the time of the magnetic field dissipation, which
is of the order of 104 years, was taken as a unit of time. The time unit choice is determined by the need to study the magnetic field modes on a large time scales. In order to reduce the number of variable parameters, the MHD-system was made dimensionless so that the control parameters are the Reynolds number, which carries information about the large-scale generator, and the amplitude of the a-effect, which characterizes the turbulent generator. The results of numerical simulation of the modes of magnetic field generation are displayed on the phase plane of the control parameters. The question of the dynamics of the change in the pattern on the phase plane depending on the waiting time determined by the scale factor of the exponential kernel J(t) of the function Z(t) is investigated.
Keywords: aQ-dynamo; magnetohydrodynamic system; modes of magnetic field generation; reversal.
References
1. Merril R.T., McElhinny M.W., McFadden P.L. The Magnetic Field of the Earth: Paleomagnetism, the Core, and the Deep Mantle. London, Academic Press, 1996.
2. Steenbek M., Krause F. Zur dynamotheorie stellarer and planetarer magnetfelder. Berechnunug sonnenahnlicher wechselfeldgeneratoren. Astronomische Nachrichten, 1969, no. 291, pp. 49-84. (in German)
3. 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.
4. Vodinchar G.M., Feshenko L.K. [Reversals in the 6-Cells Convection Driven]. Bulletin KRASEC. Physical and Mathematical Sciences, 2015, vol. 11, no. 2, pp. 41-50. (in Russian) DOI: 10.18454/2313-0156-2015-11-2-41-50
5. Gledzer E.B., Dolzhanskiy F.V., Obukhov A.M. Sistemy gidrodinamicheskogo tipa i ikh primenenie [Hydrodynamic Type Systems and Their Application]. Moscow, Nauka, 1981. (in Russian)
6. Parker E.N. Hydromagnetic Dynamo Models. The Astrophysical Journal, 1955, no. 122, pp. 293-314.
7. Krause F., Raedler K.-H. Mean-Filed Magnetohydrodynamics and Dynamo Theory. Oxford, Pergamon Press, 1980.
8. Steenbek M., Krause F., Raedler K.-H. Berechnung der mittlerer lorentz-field starke ü x B für ein elektrisch leitendes medium in turbulenter, durch Coriolis-Krafte beenflusster bewegung. Zeitschrift fur Naturforschung, 1996, no. 21, pp. 369-376. (in German)
9. Vodinchar G.M [Using Modes of Free Oscillation of a Rotating Viscous Fluid in the Large-Scale Dinamo]. Bulletin KRASEC. Physical and Mathematical Sciences, 2013, vol. 7, no. 2, pp. 33-42. (in Russian) DOI: 10.18454/2079-6641-2013-7-2-33-42
10. Godomskaya A.N., Sheremetyeva O.V. Reversals in the Low-Mode Model Dynamo with aQ-Generator. E3S Web of Conferences, 2018, vol. 62, 6 p. DOI: 10.1051/e3sconf/20186202016.
11. Godomskaya A.N., Sheremetyeva O.V. Modes of Magnetic Field Generation in Models of a aQ-dynamo with a Power Type a-generator. E3S Web of Conferences, 2019, vol. 127, article ID: 02023, 7 p. DOI: 10.1051/e3sconf/201912702016
12. Godomskaya A.N., Sheremetyeva O.V. [Dependence of the Magnetic Field Generation in the Model of a aQ-Dynamo on the Generator Intensity of a Power Type]. Bulletin KRASEC. Physical and Mathematical Sciences, 2019, vol. 29, no 4, pp. 58-66. (in Russian) DOI: 9.26107/2079-6641-2019-29-4-58-66
Received September 17, 2020