CC BY
22
DOI: 10.18454/2079-6641-2015-11-2-55-60
УДК 517.958:537.84
инверсии магнитного поля в модели крупномасштабного
аО -динамо
Г.М. Водинчар1, 2, А.Н. Годомскяя1, О.В. Шереметьева1, 2
1 Институт космоФизичєских исследований и распространения радиоволн ДВО РАН,
684034, Камчатский край, с. Паратунка, ул. Мирная, 7
2 Камчатский государственный университет имени Витуса Беринга,
683032, г. Петропавловск-Камчатский край, ул.Пограничная, 4
E-mail: [email protected], [email protected], [email protected]
В работе исследуется вопрос о возможности возникновения инверсий в рамках маломодовой модели а О -динамо. Определены параметры МГД-системы, при которых возможны инверсии магнитного поля при относительном постоянстве поля скоростей вязкой проводящей замагниченной жидкости. Обсуждаются результаты численного решения системы в предположении различного вида зависимостей амплитуды а-эффекта от радиуса.
Ключевые слова: крупномасштабная модель динамо, альфа-омега динамо, магнитное поле, инверсии.
(с) Водинчар Г.М., Годомская А.Н., Шереметьева О.В., 2015
MSC 65N80
reversal of magnetic field in the large-scale аО -DINAMO
G.M. Vodinehar1, 2, A.N. Godomskaya1, O.V. Sheremetyeva1, 2
1 Institute of Cosmophysical Researches and Radio Wave Propagation Far-Eastern Branch, Russian Academy of Sciences, 684034, Kamchatskiy Kray, Paratunka, Mirnaya st., 7, Russia
2 Vitus Bering Kamchatka State University, 683031, Petropavlovsk-Kamchatsky, Pogranichnaya st., 4, Russia
E-mail: [email protected], [email protected], [email protected]
In the paper investigates the question of the possibility of a reversal in the framework of low-mode model, а О -dinamo. The parameters of the MHD system in which the possible reversal of the magnetic field in the relative constancy of the velocity field are defined. There are results of numerical solution of the assumption of various type of а -effect amplitude dependence from the radius.
Key words: а О-dynamo, magnetic field, reversal, large-scale dynamo.
(c) Vodinchar G.M., Godomskaya A.N., Sheremetyeva O.V., 2015
Введение
Для астрофизических объектов (планет, звезд, галактических дисков) типичен механизм аО-динэмо [1]. Данный вид динамо предполагает дифференциальное вращение объекта и турбулентный характер движения проводящей среды в этом объекте [2, 3, 4, 5]. В реальных космических динамо-системах наблюдаются инверсии, в том числе и на Земле. Информацию об инверсиях геомагнитного поля получают из палеомагнитных записей, на основе которых строится шкала геомагнитной полярности. Интервалы времени между инверсиями (интервалы полярности) отличаются на несколько порядков, существуют длительные интервалы без инверсий [6]. Кроме того, следует отметить, что периоды инверсий в магнитном поле значительно меньше нежели в поле скорости хорошо проводящей вязкой жидкости. Построение модели динамо, полностью воспроизводящей реальную палеомагнитную шкалу, является трудновыполнимой задачей. Различные модели динамо позволяют получить случайные последовательности инверсий, свойства которых сильно отличаются в зависимости от выбора параметров.
Основной целью работы в рамках принятой модели динамо является получение решений системы магнитогидродинамического типа, при которых могут возникать инверсии в магнитном поле при относительном постоянстве поля скорости.
Постановка задачи
Будем рассматривать модель аО-динэмо в предположении аксиальной симметричности поля скорости вязкой жидкости u и магнитного поля B. Тогда пространственная структура среднего поля U простая и можно ограничиться одномодовыми приближениями для тороидальной и полоидэльной компонент этих полей, которые можно описать скалярными функциями:
v = UT vi(r) + UP V2(r), (1)
B = BT Bi(r) + BPB2(r), (1)
где компоненты поля скорости и магнитного поля считаются независимыми от времени и составляющая BP(r) является дипольной
vi(r) = VT (r) V2(r) = vP(r), Bi(r) = BT (r), B2(r) = BP(r).
(2)
Полагаем, что среднее течение U носит характер дифференциального вращения [2, 3, 4, 5]. С учетом сделанных предположений нижеприведенную систему МГД-уравнений крупномасштабного динамо
dV + (vV) v + 20. х v + О x (О х r)
Д* = V х (v х B) + VmAB,
Vv = VB = 0,
VAv
1 VP + — (V x B) x B + f, P PM
(3)
где v - поле скорости вязкой жидкости, B - магнитное поле, f - внешняя сила (источник полоидэльной скорости), P - давление, r - радиус-вектор, р - плотность, v
- кинематическая вязкость, vm - магнитная вязкость, д - магнитная проницаемость, можно записать в маломодовом приближении, используя метод Галёркина [3], в следующем виде
duT тр Rem T , Rem r, P . t r>T r>P
—— = RemB112u U----— Діи + р12u + L\\2B b ,
dt Re ERe
du
P
dt
dBT
dt
dBP
dt
= RemB211(uT)2---B2uP + F m P21uT + L211 (bT)2 + L222(bF)2 + f,
Re E ■ Re
= Rem(W112UTBP + W121UPBT)+ RaW02BP - V1BT,
= RemW222UPBP + Ra Wa1BT - V2BP,
(4)
где uT - тороидальная мода скорости, ир - полоидальная мода скорости, BT - тороидальная мода магнитного поля, Bp - полоидальная мода магнитного поля, Re -число Рейнольдса, Rem - магнитное число Рейнольдса, Ra - амплитуда а-эффекта, E - число Экмана, f - внешняя сила, д - коэффициент вязкой диссипации, Vi -коэффициент магнитной диссипации. Случайным образом задаются значения Re = 10 ^ 104, Rem = 104 ^ 1010, E = 10-10 ^ 10-5 в границах изменения их возможных значений. Коэффициенты Bijk,Wijk,Wij,P12 - это объемные интегралы от рассматриваемых полей, значения которых для принятой модели принимаются равными:
B112 = 0.365, B211 = -0.365, B222 = -0.170, P12 = 0.193, Lm = 0.487, L211 = 0.211, W112 = -0.487, W121 = -0.211, W222 = 0.170, v1 = 33.218, V2 = 9.870, д1 = 28.159, Д2 = 86.573.
В зависимости от определения алгебраической функции для задания амплитуды а-эффекта коэффициенты Wa принимают следующие значения:
1) если зависимость постоянная a(r) = 1 : W^ = 1.874,W2°1 = 1.037,
2) если зависимость линейная a(r) = r: W/2 = 1.073, W2°1 = 0.788,
3) если зависимость от радиуса синусоидальная a(r) = sinr: W/2 = -0.506, W2a = 0.117.
Результаты моделирования
В предположении относительного постоянства поля скорости примем значения мод скорости постоянными. Тогда наличие инверсий в магнитном поле будет возможным лишь при условии наличия осциллирующего решения для двух последних дифференциальных уравнений системы (4):
' dB т_ dt
' dBP dt
Rem(W112uTBP + W121uPBT) + RaW1a2BF - V1BT,
RemW222uPBP + RaW2°1BT - V2BP,
(5)
с неизвестными параметрами Ra, RemuP, RemuT, области возможного изменения которых в логарифмическом масштабе приняты следующими
lgRa Є [-1,3], lg(RemuP) є [3,7], lg(RemuT) Є [2,8].
Осциллирующее решение системы (5) может быть получено лишь при условии отрицательности дискриминанта ее характеристического уравнения
RemW121 up - Vl RemWn2UT + RaW1a2
RaW21 RemW222UP - V2
0
и получаемое решение будет устойчивым по Ляпунову [8], если положительны коэффициенты характеристического уравнения, записанного в общем виде Я2+рЯ + q = 0:
p = 43.088 + 0.041RemUP,
q = -0.035870Rem(uP)2 - 3.564490RemUP + 0.487W2iRa RemUT + 327.861660 - W^R^'W12,
Исходя из принятых ограничений получена область в пространстве рассматриваемых параметров Ra, RemUP, RemUT, в которой магнитное поле имеет колебательный характер (рис. 1).
Рис. 1. Проекции области осциллирующего решения системы (5) в пространстве параметров Ra, RemUP, RemUT
Численные расчеты показали, что инверсии в магнитном поле возникают в большинстве точек полученной области. Однако, как правило, и поле скорости, и магнитное поле быстро затухают, причем периоды осцилляций в рассматриваемых полях практически одинаковы (рис. 2 - 3).
а)
б)
Рис. 2. Изменения тороидальной и полоидальной мод поля скорости и магнитного поля при задании амплитуды a-эффекта линейной зависимостью a(r) = 1
а)
б)
Рис. 3. Изменения тороидальной и полоидальной мод поля скорости и магнитного поля при задании амплитуды а-эффекта линейной зависимостью а(r) = r: а) зависимость компонент рассматриваемых полей от времени; б) инверсии компонент рассматриваемых полей
Инверсии в магнитном поле, возникающие на фоне более медленных изменений в поле скорости, характерны для случая синусоидальной зависимости амплитуды а-эффекта от радиуса (рис. 4). Но и в этом случае магнитное поле быстро затухает и его колебания происходят вблизи нулевого значения.
а)
б)
Рис. 4. Изменения тороидальной и полоидальной мод поля скорости и магнитного поля при задании амплитуды а-эффекта нелинейной зависимостью a (r) = sin r
Заключение
Исходя из принятой маломодовой модели аП-динамо в предположении аксиальной симметрии полей были получены инверсии в магнитном поле при реалистичных областях изменения параметров МГД-системы (4). Однако данная модель не позво-
лила в пределах принятых ограничений параметров получить инверсии на фоне практически постоянного поля скорости. Во всех случаях, за исключением нелинейной зависимости амплитуды а-эффекта от радиуса, периоды инверсии поля скоростей сравнимы с периодами инверсий магнитного поля.
Модели данного класса предполагают значительные упрощения и отражают наиболее общие свойства, поэтому для описания механизма инверсий необходим корректный выбор мод рассматриваемых полей. Полученные результаты позволяют определить направление уточнения выбранной мэломодовой модели аП-динэмо для описания рассматриваемого механизма инверсий.
Библиографический список
1. Parker E.N. Hydromagnetie dynamo models // Astrophys. J. 1955. № 122. p. 293-314.
2. Krause, F. and RAdler, K.-H. Mean-filed magnetohydrodynamies and dynamo theory. Pergamon Press. 1980.
3. Steenbek M., Krause F. and Radler K.-H. Bereehnung der mittlerer Lorentz-Field Starke v x B fUr ein elektriseh leitendes Medium in turbulenter, dureh Coriolis—KrAfte beenflusster Bewegung // Z. Naturforseh. 1996. № 21. p. 369-376.
4. Steenbek M. and Krause F. Zur Dynamotheorie stellarer and planetarer Magnetfelder. I. Bereehnunug sonnenahnlieher Weehselfeldgeneratoren // Astron. Naehr. 1969. № 291. p. 49-84.
5. Zeldovieh, Ya.B., Rusmaikin, A.A., and Sokoloff, D.D.: Magnetic fields in astrophysics. The Fluid Mechanics of Astrophysics and Geophysics, Gordon and Breach, New York, 1983.
6. Merril R.T., MeElhinny M.W. and MeFadden P.L. The Magnetic Field of the Earth: Paleomagnetism, the Core, and the Deep Mantle. Aeademie Press. London. 1996.
7. Гледзер Е.Б., Должэнский Ф.В., Обухов А.М. Системы гидродинамического типа и их применение. М.: Наука, 1981. 368 с.
8. Эльсгольц Л.Э. Дифференциальные уравнения и вариационное исчисление. М: Наука, 1969. 424 с.
Поступила в редакцию / Original article submitted: 16.11.2015
