Кинематическая генерация длинномасштабного магнитного поля течениями с зеркальной антисимметрией Текст научной статьи по специальности «Физика»

www.volsu.ru

ФИЗИКА И АСТРОНОМИЯ

DOI: https://doi.Org/10.15688/mpcm.jvolsu.2018.2.7

УДК 537.84+52-337 ББК В253.350

кинематическая генерация длинномасштабного магнитного поля

________о _______ о л

течениями с зеркальной антисимметриеи 1

Александр Андреевич Андриевский

Системный архитектор,

Общество с ограниченной ответственностью «ЦПР РТСофт» [email protected]

ул. Никитинская, 3, стр. 1, 105037 г. Москва, Российская Федерация

оо Владислав Александрович Желиговский

о Доктор физико-математических наук, главный научный сотрудник,

, Институт теории прогноза землетрясений и математической геофизики РАН

^ [email protected]

х ул. Профсоюзная, 84/32, 117997 г. Москва, Российская Федерация

«

о m

о

& Роман Александрович Чертовских

^ Кандидат физико-математических наук, доктор философии (PhD), исследователь,

< Исследовательский центр систем и технологий (SYSTEC), Университет Порто,

Португалия

m

s [email protected]

«

Rua Dr. Roberto Frias, s/n 4200-465 Porto, Portugal; ^ старший научный сотрудник, Самарский национальный исследовательский

5 университет им. акад. С.П. Королева

^ Московское шоссе, 34, 443086 г. Самара, Российская Федерация

<

3S

g Аннотация. Работа посвящена численному моделированию кинематиче-

о

g скои генерации магнитного поля, в котором присутствуют длинные простран-

S

ственные масштабы, короткомасштабным течением несжимаемой электропроводной жидкости с определенной антисимметрией. Рассмотренный механизм @ генерации работает в широком диапазоне магнитных чисел Прандтля, что

важно для астрофизических приложений. Численное моделирование проведено для течения, поле скорости которого имеет поточечно нулевую кинетическую спиральность; показано, что это не препятствует генерации магнитного поля.

Ключевые слова: магнитогидродинамика, магнитное поле, магнитный а-эффект, кинематическое динамо, спиральность.

Введение

В настоящее время общепризнано, что источник магнитного поля многих астрофизических объектов (например, Солнца) — течение электропроводной жидкости [3]. Магнитогидродинамическая генерация описывается уравнениями Навье — Стокса и магнитной индукции (в совокупности, например, с уравнением теплопередачи), которые, в принципе, можно использовать для моделирования. Однако такой подход упирается в необходимость проводить вычисления с очень тонким пространственным разрешением, что требует чрезмерных вычислительных ресурсов. Желательно ограничить рассмотрение задачи на достаточно большие масштабы. Для этого теория электродинамики средних полей [2; 8] предлагает различные оценки воздействия малых масштабов, не участвующих в моделировании, на большие (аналогично замыканиям метода крупных вихрей в вычислительной гидродинамике [10]). Недостаток этого подхода состоит в необходимости принимать дополнительные гипотезы для вывода таких оценок. Однако при сильном разделении масштабов их вывод можно провести математически строго асимптотическими методами осреднения эллиптических операторов (см. [12] и ссылки в этой монографии). Это позволяет при изучении эволюции длинномасштабных линейных и слабо-нелинейных возмущений в теории магнитогидродинамической устойчивости независимо вычислять решения так называемых вспомогательных задач, описывающих короткомасштабные процессы и их воздействие на длинномасштабные процессы, и амплитудных уравнений, определяющих поведение длинномасштабных средних полей возмущений [5]. В данной работе мы следуем этому подходу к численному моделированию генерации магнитного поля стационарным течением несжимаемой электропроводной жидкости в рамках задачи о кинематическом динамо [3].

При рассмотрении кинематической генерации магнитного поля при наличии существенного разделения пространственных масштабов в рамках теории магнитогидродина-мической устойчивости к длинномасштабным возмущениям сделан вывод (см., например, [12]), что два наиболее важных для генерации явления — а-эффект и вихревая диффузия — действуют на разных временных масштабах. При наличии магнитного а-эффекта эволюция среднего магнитного поля Ь(Х) подчиняется в главном порядке уравнению

д Ь

— = Ух х (АЬ) (1)

(отсутствие в нем диссипативных операторов второго порядка, отвечающих молекулярной или турбулентной диффузии — важное обстоятельство). Здесь Та = — так называемое медленное время; малый параметр задачи е — отношение пространственных масштабов, описываемых быстрой пространственной переменной х порядка характерного размера структур генерирующего потока и медленной переменной X = ех; £ — быстрое время порядка периода оборота вихря в жидкости; А — тензор магнитного а-эффекта.

(Поля, зависящие только от быстрой пространственной переменной х, будем называть короткомасштабными, а зависящие также и от медленной переменной X = ех — длин-номасштабными.) Если же а-эффект отсутствует (например, когда течение имеет центр симметрии), то эволюция среднего магнитного поля в главном порядке удовлетворяет уравнению

Здесь Dmk — тензор магнитной вихревой («турбулентной») диффузии.

Рассмотрим (следуя [5]) эволюцию средних магнитных полей, периодических по пространству, которые в начальный момент имеют вид длинномасштабных гармоник Фурье Не1фХ (где H — вектор-константа) и являются собственными функциями операторов а-эффекта или вихревой диффузии (задаваемых правыми частями уравнений (1) и (2) соответственно). Спектр оператора а-эффекта симметричен относительно мнимой оси [1], так как если гармонике Не1чХ отвечает собственное значение А, то гармонике Не1тч-Х, где т — целое (при этом условии гармоника с растянутым в т раз волновым вектором удовлетворяет условию периодичности для первой), отвечает собственное значение тА. Таким образом, если действительная часть А не равна 0, то новая гармоника Heimq-X при соответствующем выборе знака т вырастает в е раз за время (e|mReА|)-1 (в масштабе быстрого времени), а начальное поле, имеющее вид бесконечной суммы таких гармоник (то есть бесконечного ряда Фурье), будет расти сверхэкспоненциально. Оператор магнитной вихревой диффузии — второго порядка, поэтому если гармонике Не1чХ отвечает собственное значение А, то гармонике Не1тфХ (где т — также произвольное целое число) отвечает собственное значение т2А. В отличие от случая оператора а-эффекта, растяжение в т раз волнового вектора q не меняет знак собственного значения; однако если собственное значение А c положительной действительной частью существует (то есть имеет место явление отрицательной магнитной вихревой диффузии), то после такого растяжения новая гармоника вырастает в е раз за (быстрое) время (гт2КеА)-1, а начальное поле, имеющее вид бесконечной суммы таких собственных функций (бесконечного ряда Фурье), растет сверхэкспоненциально, как и в случае оператора а-эффекта.

С математической точки зрения задача кинематического динамо — это часть задачи о линейной устойчивости немагнитного магнитогидродинамического состояния: поскольку магнитное поле воздействует на течение электропроводной жидкости посредством силы Лоренца, квадратичной по магнитному полю, задача о линейной устойчивости чисто гидродинамического состояния распадается на независимые подзадачи об устойчивости к гидродинамическим и магнитным возмущениям; вторая из них и есть задача кинематического динамо. Когда амплитуда начального возмущения достаточно вырастает, из-за действия силы Лоренца течение подвергается значительному возмущению, и оно, в частности, может утерять способность генерировать магнитное поле. Из приведенного выше анализа вытекает, что при действии а-эффекта самоуничтожение динамо вследствие неустойчивости может произойти на временах (в масштабе быстрого времени) порядка O(e-1), а при его отсутствии и действии отрицательной вихревой диффузии — на временах порядка O(e-2). Таким образом, в рамках гипотезы о магнито-гидродинамическом характере источника космических магнитных полей, по-видимому, астрофизическим приложениям в наибольшей степени отвечает такое устройство динамо, при котором все собственные значения оператора а-эффекта чисто мнимы, а рост

(2)

магнитного поля на больших временах обеспечивает отрицательная магнитная вихревая диффузия.

Работа посвящена изучению генерации магнитного поля течениями с полем скорости, имеющим зеркальную антисимметрию относительно некоторой плоскости (в дальнейшем считаем, что это координатная плоскость х1 = 0). Поле ¥ = (/1, /2, /з) называют зеркально-симметричным относительно плоскости х1 = 0 (или симметричным относительно декартовой координаты ж1), если при всех г выполнено

/¿(-Ж1,Ж2,Жз) = (-1)61 /¿(х),

и зеркально-антисимметричным относительно этой плоскости (или антисимметричным относительно ж1), если при всех г выполнено

/¿(-Ж1,Ж2,Жз) = (-1)1-81 /г(х);

здесь — символ Кронекера. Мы показываем, что течения из этого класса удовлетворяют сформулированному выше условию «максимально долгой работы динамо». Под действием а-эффекта в таких течениях среднее магнитное поле совершает гармонические колебания с конечным периодом в медленном времени Т1, и возникает анизотропная магнитная вихревая диффузия, способная при любой величине магнитной молекулярной диффузии п генерировать длинномасштабное поле в масштабе медленного времени Т2. Астрофизические динамо функционируют при магнитных числах Прандтля Рт из очень широкого интервала (см., например, [3;6]), от весьма малых (например, в межзвездной среде) до весьма больших (например, геодинамо), при этом динамо при Рт ^ 1 более сложны для анализа (см. [4]). Поскольку рассматриваемый здесь механизм генерации работоспособен при произвольной молекулярной диффузии, он может представлять определенную ценность для астрофизических приложений.

В теории электродинамики средних полей [2] показано, что если время корреляции турбулентного течения V мало (по сравнению с характерным временем диффузии магнитного поля), то след тензора а-эффекта пропорционален средней кинетической спи-ральности V • (Ух V) (являющейся спиральностью завихренности Ух V). Вместе с тем известно [9], что способность гладкого (ламинарного) течения генерировать магнитное поле не связана с наличием у него ненулевой кинетической спиральности. Чтобы проиллюстрировать это, мы выбрали для численного исследования течение из так называемого семейства Ь (определенного в [9]), которое ортогонально завихренности в каждой точке пространства, то есть кинетическая спиральность которого поточечно равна 0.

1. Формализм теории усреднения

Задача о генерации магнитного поля стационарным потоком v(x) эквивалентна задаче на собственные значения

£Ь = ЛЬ, (3)

где

£Ь = пУ2Ь + Ух (V х Ь) — (4)

оператор магнитной индукции. Короткомасштабную скорость течения v(x) и длинно-масштабные магнитные моды Ь(х, X) считаем соленоидальными и 2п-периодическими

по быстрым переменным кроме того, v(x) имеет нулевое среднее по ячейке периодичности Т3 быстрых переменных. Средняя по мелким масштабам и флуктуирующая части поля ¥(х, X) определяются равенствами

<f (x, X)) = / f (x, X) = E (f)k , if} = f - (f)

J T3 (2п)3 f-1

к=1

где ек — орты декартовой системы координат.

В дальнейшем будем использовать следующие обозначения: дифференциальные операторы по медленным переменным маркируем нижним индексом X, а в операторе £, определенном равенством (4), предполагаем дифференцирование только по быстрым переменным.

Разложим решение уравнения (3) в степенные ряды по малому параметру задачи:

те

-П

Ь(х, Х) = £ Ьга(х, Х)еп, (5.1)

п=0 те

Л = ^ лга£га. (5.2)

п=0

Подстановка этих рядов в (3) дает

те / п \

Е I ^Ьп + п(2(У ■ Ух)Ьп-1 + У^Ьп-2) + Vx х (V х Ьп-1) — Е Л п—т11т еп = 0. п=0 \ т=0 /

(6)

Аналогично из условия соленоидальности магнитной моды находим, что для всех п > 0 выполнены равенства

Vx -(Ьп) = 0, (7.1)

V ■ Ьп + Vx ■ {Ьп-1} = 0. (7.2)

Приравняв к нулю коэффициенты ряда в левой части (6), получим иерархию уравнений, из которой с учетом (7.1) можно последовательно определить все коэффициенты разложений (5) [12].

1.1. Уравнение порядка е0

При п = 0 получаем из (6) первое уравнение в иерархии:

£Ь0 = Л0Ь0, (8)

усреднением которого находим 0 = Л0 (Ь0). При малых е собственное значение (5.2) мало отличается от Л0, поэтому, если действительная часть Л0 не равна нулю, то при малых е длинномасштабная мода (5.1) растет или затухает одновременно с короткомасштабной модой Ь0, отвечающей собственному значению Л0. Следовательно, при таком Л0 длинно-масштабность моды существенно не влияет на физический механизм работы изучаемого динамо. Случай чисто мнимого Л0 возможен, но не является случаем общего положения, и потому также не представляет особенного интереса. Соответственно, мы рассмотрим

длинномасштабное динамо, для которого Ло = 0. Поскольку векторные поля-константы входят в ядро оператора, сопряженного к £, размерность ядра £ не менее 3. Ограничимся рассмотрением случая общего положения, когда ядро короткомасштабного оператора £ трехмерно.

При этих предположениях решение (8) имеет вид

з

Ьо = ^ <Ьо)Л , (9)

к=1

где соленоидальные короткомасштабные нейтральные моды (х) — решения вспомогательных задач типа I:

£Бк = 0, ) = ек, У • = 0

(доказательство их существования приведено в [12]).

1.2. Уравнение порядка е1

При п =1 получаем из (6) следующее уравнение из иерархии:

£Ь1 + 2п(У • Ух)Ьо + Ух х (V х Ьо) = Л^. (10)

Подстановкой (9) и усреднением находим условие его разрешимости (совпадающее с условием ортогональности в пространстве Лебега Ь2(Т3) неоднородности ядру сопряженного оператора, то есть векторным полям-константам):

Ух х (А <Ьо)) = Л1 <Ьо) , (11)

где тензор магнитного а-эффекта А состоит из столбцов

Ак = <v х ).

Для поля скорости, антисимметричного относительно декартовой координаты ж1, он имеет следующую структуру:

0 A

A = 0 A2 . (12)

A -A3 0

Следуя [7;9; 12], рассмотрим пространственно-периодическое поле

(ho) = He^X (13.1)

Волновой вектор q удобно задать в сферических координатах с осью, направленной вдоль декартовой оси х\:

q = (cos 0, sin 0 cos ф, sin 0 sin ф). (13.2)

Соленоидальность (Ь0) (см. (7.1) для п = 0) эквивалентна ортогональности постоянных векторов Н и q. С учетом этой ортогональности, подставив (12) и (13) в (11), находим

А± = i ± (Afo - A^s))

H+ H-

(-(Al cos ф - Al sin ф) sin 0, Al cos 0, -A2 cos 0), (0, - sin ф, cos ф).

(14.1)

(14.2)

Таким образом, собственные значения оператора а-эффекта для течения, антисимметричного относительно декартовой координаты, чисто мнимы.

После подстановки (13.1) флуктуирующая часть уравнения (10) теперь имеет вид

Lhi + eiqX Отсюда

t

\ к=1

Нк(q ■ V)sfc + iq X < v X

hi

fc=i \

< v X Y^ HkSfc > - Ai Y Hk(Sfc} I

l k=1 J k=1 )

(hi)fc sk + eiqXHk Aiyfc + i ^ Qmgmk

(15)

где ук(х) и gmk(х) — короткомасштабные пространственно-периодические решения с нулевым средним вспомогательных задач типов II и II

д б

Lgmk = -2п~--ет х {V х }, £ук = {Б*,}.

ОХ™

1.3. Уравнение порядка е2

При п = 2 получаем из (6) уравнение

£Ь2 + 2л(У ■ Ух)Ь + пУХЬ + Ух х (V х Ь^ = Л^1 + Л2Ьо. Его усреднение после подстановки (9) и (15) дает

Vx х (A (hi))+ieiq^Xq х

£

k=i \

)

к | AiDDк + i ^ qmDmk) = Ai (hi) + (А2 + лУфХН, (16)

m=1

где тензоры вихревой диффузии состоят из столбцов

© к = (V х ук) , ©тк = (V х %тк) .

(17)

Для поля скорости, антисимметричного относительно декартовой координаты ж1, их структура такова:

D

D

D 2 DD 1

2 о Di 32

i -D 2 0

0 Di2 Dii

Di D12 0 D?s

D3 Dii D2 Dis 0

(18.1)

D1 —D1 —D1

D21 D22 d2S

D22 D2S

d2

d22

D2

d2S

d2

d2S

Ds

D2S

D1 —D1 —d1

DS1 DS2 Dss

D1 D32

D1 D33

D2 D32

D2 D33

D2 D33

D3 D33

(18.2)

В (18.2) каждая матрица размера 3x3 состоит из столбцов Dmk (17) при фиксированном т, возрастающем слева направо от 1 до 3.

Решение (16) имеет вид гармоники Фурье: (hi) = eiq^XH'. Постоянный вектор H' ортогонален q в силу соленоидальности (7.1) для п = 1 и удовлетворяет уравнению

iq x ^ нк(AiDк + i J] qmDmA = -iq x (AH') + AiH' + (Л2 + n)H.

k=i V m=1 /

Отсюда с учетом (12), (14) и (18) находим, что при Л+ = Л- (случай общего положения) Л± = Qi + Q2 + (Qi - Q2) cos 20 ± Q3 sin 20 - n, (19.1)

где

Qi = - 1(D 3a2 + D2s), (19.2)

Q2 = 1 ((Dia + + AiD2 - A3D3) cos 2Ф + (D^ - ^2 + A?Di + A2Di) sin 2Ф+

+ Dia - D32 - AiDi - A?D^ , (19.3)

^(D2a + D22 - Dia - A?D2 - A3D?) - ^(D3a + D32 - D^ - A3Di - Ai>D3)-

2AiD3a + 2A3D22 - (a3(2D32 + D23 - D22) + Ai(2D23 + D32 - D33)) sin 2Ф

4(A2 sin ф - A3 cos ф) '

(19.4)

Соответственно, когда Л± мнимы (в том числе когда поле скорости антисимметрично по декартовой координате — см. (14.1)), Л± действительны.

2. Численные результаты

В [9] рассмотрены шесть семейств течений несжимаемой жидкости с нулевой кинетической спиральностью. Течения так называемого семейства L имеют следующий вид:

v(x) = AVB - BVA. (20)

Если потенциалы Монжа А и В — скалярные собственные функции лапласиана, отвечающие одному и тому же собственному значению, поток (20) соленоидален. Для численного исследования выбран поток данного вида с 2п-периодическими по всем декартовым координатам потенциалами, отвечающими собственному значению -18. Они состоят из гармоник Фурье с волновыми числами (±3, ±3, 0) и (±4, ±1, ±1) и их перестановками (всего 36 волновых векторов). Когда А нечетно, а В четно по Х\, поток (20) антисимметричен по этой координате. Соответственно, выбранные потенциалы — линейные комбинации из 16 (в случае А) или 20 (в случае В) произведений тригонометрических функций

П {1}

i=1 ^ J

с псевдослучайными коэффициентами, которые равномерно распределены на отрезке [-1,1], где А не содержит cos п\Х\, а В — sin п\Х\. Поле нормализовано так, что его среднеквадратичная величина равна 1. Рисунки изоповерхностей плотностей кинетической энергии i|v|2 (рис. 1а) и завихренности V x v (рис. 1б) рассматриваемого течения показывают, что оно имеет достаточно нетривиальную структуру, что, по-видимому, положительно сказывается на способности этого течения генерировать магнитное поле.

Будучи доминирующими (то есть имеющими максимальный инкремент роста), ко-роткомасштабные нейтральные магнитные моды sk (x) вычислены кодом [11] с применением псевдоспектральных методов. График коэффициента магнитной вихревой диффузии

neddy(Ф) = - max Л±

(рис. 2) иллюстрирует сингулярное поведение инкрементов роста магнитного поля Л± (19) как функции азимутального направления ф в некоторой малой окрестности критического азимутального направления волнового вектора q, при котором знаменатель дроби, входящей в выражение (19.4) для Q3, становится равен нулю. Такое сингулярное поведение имеет место при любом коэффициенте магнитной молекулярной диффузии п. При этом, в любой малой окрестности критического азимутального направления знаменатель дроби, входящей в (19.4), принимает как положительные, так и отрицательные значения. Соответственно, в некоторой малой окрестности критического азимутального направления один из инкрементов Л± положителен, каково бы ни было п. Отметим также жесткую фокусировку рассматриваемого динамо: вне некоторого малого интервала ф вокруг критического направления инкремент роста (убывания) магнитного поля Л± (19) очень слабо зависит от азимутального направления, то есть генерируемое магнитное поле имеет вид волны, распространяющейся в пространстве вдоль четко выраженного конуса направлений.

а б

Рис. 1. Изоповерхности плотности кинетической энергии поля скорости течения (а) и его завихренности (б) на уровнях 15 % и 25 % максимумов соответственно

Выводы

Нами описано двумасштабное динамо, работа которого основана на совместном функционировании двух важнейших механизмов генерации длинномасштабного магнитного поля: а-эффекта и анизотропной отрицательной вихревой диффузии, ответственных за гармонические колебания магнитного поля и его рост соответственно.

Это динамо работает при любых коэффициентах магнитной молекулярной диффузии, когда разделение масштабов настолько велико, что первые члены асимптотических разложений (5) достаточно точно описывают эволюцию доминирующей магнитной моды.

Обращение в ноль кинетической спиральности течения не препятствует функционированию ни магнитного а-эффекта, ни вихревой диффузии. Представляется, что никакие функционалы вида спиральности некоторого поля не могут описать способность гладкого течения генерировать магнитное поле.

-0.5 -

Рис. 2. Зависимость коэффициента вихревой диффузии пеаау (вертикальная ось) от азимутального направления ф (горизонтальная ось) волнового вектора (13.2) при п = 0, 3

ПРИМЕЧАНИЕ

1 P.A. Чертовских был поддержан проектом P0CI-01-0145-FEDER-006933/SYSTEC, финансируемым ERDF (European Regional Development Fund) в рамках программы COMPETE 2020 (Programa Operacional Competitividade e Internacionalizacao) и фондом FCT (Fundacao para a Ciencia e a Tecnologia, Португалия). Основной объем вычислений выполнен на кластере «Сергей Королев» Самарского университета.

СПИСОК ЛИТЕРАТУРЫ

1. Вишик, М. М. Периодическое динамо II. Численное моделирование и анализ геофизических процессов / М. М. Вишик // Вычислительная сейсмология. — 1987. — Т. 20. — С. 12-22.

2. Краузе, Ф. Магнитная гидродинамика средних полей и теория динамо / Ф. Краузе, К.-Х. Рэдлер. — М. : Мир, 1984. — 320 с.

3. Моффат, Г. Возбуждение магнитного поля в проводящей среде / Г. Моффат. — М. : Мир, 1980. — 339 с.

4. Brandenburg, A. Dissipation in dynamos at low and high magnetic Prandtl numbers / A. Brandenburg // Astron. Nachr. - 2011. - Iss. 332. - P. 725-731.

5. Chertovskih, R. Large-scale weakly nonlinear perturbations of convective magnetic dynamos in a rotating layer / R. Chertovskih, V. Zheligovsky // Physica D. — 2015. — Iss. 313. — P. 99-116.

6. Moffatt, H. K. Helicity and celestial magnetism / H. K. Moffatt // Proc. R. Soc. A. — 2015. — Iss. 472. — Article ID: 20160183. — DOI: 10.1098/rspa.2016.0183.

7. Negative magnetic eddy diffusivities from test-field method and multiscale stability theory / A. Andrievsky, A. Brandenburg, A. Noullez, V. Zheligovsky // Astrophysical J. — 2015. — Iss. 811. — P. 135-164.

8. Radler, K.-H. Mean-field dynamo theory: early ideas and today's problems. / K.-H. Radler // Magnetohydrodynamics. — Dordrecht : Springer, 2007. — P. 55-72. — DOI: 10.1007/978-1-4020-4833-3_4.

9. Rasskazov, A. Magnetic field generation by pointwise zero-helicity three-dimensional steady flow of an incompressible electrically conducting fluid / A. Rasskazov, R. Chertovskih, V. Zheligovsky // Phys. Rev. E. — 2018. — Iss. 97. — Article ID: 043201. — DOI: 10.1103/PhysRevE.97.043201.

10. Sagaut, P. Large eddy simulation for incomressible flows / P. Sagaut. — Berlin : Springer-Verlag, 2006. — 558 p.

11. Zheligovsky, V. Numerical solution of the kinematic dynamo problem for Beltrami flows in a sphere / V. Zheligovsky // J. Scientific Computing. — 1993. — Iss. 8. — P. 41-68.

12. Zheligovsky, V. Large-scale perturbations of magnetohydrodynamic regimes: linear and weakly nonlinear stability theory / V. Zheligovsky. — Heidelberg : Springer-Verlag, 2011. — 330 p. — DOI: 10.1007/978-3-642-18170-2.

REFERENCES

1. Vishik M.M. Periodicheskoe dinamo II. Chislennoe modelirovanie i analiz geofizicheskikh protsessov [Periodic Dynamo II. Numerical Modelling and Analysis of Geophysical Processes]. Vychislitelnaya seysmologiya, 1987, vol. 20, pp. 12-22.

2. Krause F., Radler K.-Kh. Magnitnaya gidrodinamika srednikh poley i teoriya dinamo [Mean-Field Magnetohydrodynamics and Dynamo Theory]. Moscow, Mir Publ., 1984. 320 p.

3. Moffatt H.K. Vozbuzhdenie magnitnogo polya v provodyashchey srede [Magnetic Field Generation in Electrically Conducting Fluids]. Moscow, Mir Publ., 1980. 339 p.

4. Brandenburg A. Dissipation in Dynamos at Low and High Magnetic Prandtl Numbers. Astron. Nachr, 2011, iss. 332, pp. 725-731.

5. Chertovskih R., Zheligovsky V. Large-Scale Weakly Nonlinear Perturbations of Convective Magnetic Dynamos in a Rotating Layer. Physica D, 2015, iss. 313, pp. 99-116.

6. Moffatt H.K. Helicity and Celestial Magnetism. Proc. R. Soc. A, 2015, iss. 472, article ID: 20160183. DOI: 10.1098/rspa.2016.0183.

7. Andrievsky A., Brandenburg A., Noullez A., Zheligovsky V. Negative Magnetic Eddy Diffusivities From Test-Field Method and Multiscale Stability Theory. Astrophysical J., 2015, iss. 811, pp. 135-164.

8. Radler K.-H. Mean-field dynamo theory: early ideas and today's problems. Magnetohydrodynamics. Dordrecht, Springer, 2007, pp. 55-72. DOI: 10.1007/978-1-4020-4833-3_4.

9. Rasskazov A., Chertovskih R., Zheligovsky V. Magnetic Field Generation by Pointwise Zero-Helicity Three-Dimensional Steady Flow of an Incompressible Electrically Conducting Fluid. Phys. Rev. E, 2018, iss. 97, article ID: 043201. DOI: 10.1103/PhysRevE.97.043201.

10. Sagaut P. Large eddy simulation for incomressible flows. Berlin, Springer-Verlag, 2006. 558 p.

11. Zheligovsky V. Numerical Solution of the Kinematic Dynamo Problem for Beltrami Flows in a Sphere. J. Scientific Computing, 1993, iss. 8, pp. 41-68.

12. Zheligovsky V. Large-scale perturbations of magnetohydrodynamic regimes: linear and weakly nonlinear stability theory. Heidelberg, Springer-Verlag, 2011. 330 p. DOI: 10.1007/9783-642-18170-2.

KINEMATIC DYNAMO BY PARITY-ANTISYMMETRIC FLOWS

Aleksandr Andreyevich Andrievsky

Software Architect,

SDC RTSoft Ltd

[email protected]

Nikitinskaya St., 3 bld. 1, 105037 Moscow, Russian Federation

Vladislav Aleksandrovich Zheligovsky

Doctor of Physical and Mathematical Sciences, Chief Researcher,

Institute of earthquake prediction theory and mathematical geophysics, Russian

Academy of Sciences

[email protected]

Profsoyuznaya St., 84/32, 117997 Moscow, Russian Federation

Roman Aleksandrovich Chertovskikh

Candidate of Physical and Mathematical Sciences, PhD, Researcher, Research Center for Systems and Technologies (SYSTEC), University of Porto [email protected]

Rua Dr. Roberto Frias, s/n 4200-465 Porto, Portugal;

Senior Researcher, Samara State Aerospace University (National Research University) Moskovskoe shosse, 34, 443086 Samara, Russian Federation

Abstract. The paper is devoted to mathematical and numerical modelling of kinematic generation of magnetic field, involving large spatial scales, by a small-scale flow of incompressible electrically conducting fluid featuring a mirror antisymmetry. Direct numerical simulation demonstrates that such flows can support a magnetic field generation in presence of two main generation mechanisms, namely, the magnetic alpha-effect and negative eddy diffusivity. The magnetic field generation can be described as follows:

• a-effect creates large-scale field (h0) of amplitude O(1), oscillating on a time scale O(e-1).

• Fluctuations {h0} of this field have an amplitude O(1).

• Small-scale flow creates {10} with amplitude O(e).

• Interaction of this field with small-scale flow creates an electromotive force (v x {h1}) of amplitude O(e).

• This electromotive force gives rise to an eddy diffusivity that supports a growth of a mean field (h0) on a time scale O(e-1).

Here е is the characteristic spatial scale ratio.

It may be important for applications that the mechanism for generation considered here does work in a wide range of magnetic Prandtl numbers. Numerical simulation for a flow, which velocity has a zero kinetic helicity everywhere in space, shows that the absence of helicity does not affect magnetic field generation.

Key words: magnetohydrodynamics, magnetic field, magnetic a-effect, kinematic dynamo, helicity.