Влияние магнитной и перекрестной спиральности на свободное вырождение МГД-турбулентности Текст научной статьи по специальности «Физика»

2510

Фундаментальные проблемы теоретической и прикладной механики Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (5), с. 2510-2512

УДК 537.84

ВЛИЯНИЕ МАГНИТНОЙ И ПЕРЕКРЕСТНОЙ СПИРАЛЬНОСТИ НА СВОБОДНОЕ ВЫРОЖДЕНИЕ МГД-ТУРБУЛЕНТНОСТИ

© 2011 г. Р.А. Степанов

Институт механики сплошных сред УрО РАН, Пермь

rodion@icmm.ru

Поступила в редакцию 24.08.2011

Изучается долговременное свободное вырождение магнитогидродинамической (МГД) турбулентности при больших числах Рейнольдса с использованием каскадных моделей. Статистические свойства численных решений получены на основе моделирования 128 реализаций процесса вырождения продолжительностью, превышающей время оборота крупномасштабного вихря на 5 порядков. Результаты расчетов демонстрируют возможность реализации двух сценариев эволюции системы при слабой вариации начального состояния. Ранее известный первый сценарий, характерный для большинства реализаций, соответствует быстрому накоплению перекрестной спиральности, что приводит к состоянию полной коррелированности магнитного поля и поля скорости. Следуя второму, более редкому, сценарию вырождения, система остается при слабом уровне перекрестной спиральности и достигает существенно более низкого значения кинетической энергии. Впервые удалось установить, что второй сценарий вырождения обусловлен влиянием магнитной спиральности, которая блокирует часть магнитной энергии в крупномасштабной части спектра.

Ключевые слова: магнитная гидродинамика, теория турбулентности, свободное вырождение, каскадные модели, законы сохранения.

Введение

Свободное вырождение МГД-турбулентнос-ти представляет интерес по двум основным причинам. Во-первых, оно имеет отношение к физике межзвездного пространства и космологии в контексте проблем эволюции начального магнитного поля и его вклада в наблюдаемые магнитные поля во Вселенной. Во-вторых, МГД-турбу-лентность отличается от обычной турбулентности своими законами сохранения, что может приводить к принципиально иному характеру эволюции турбулентности в электропроводящей среде. В пределе идеальной МГД-системы имеется три закона: сохранения полной энергии E = Eи + Eb, перекрестной спиральности Hc = ££[ v-b и магнитной спиральности Hb =£j^ a-b, где v — поле скорости, b — магнитное поле, a — векторный потенциал (rot a = b).

Первая попытка моделирования поведения МГД-турбулентности на поздних стадиях вырождения с помощью каскадных моделей была предпринята в работе [1], где и было показано, что, несмотря на общую тенденцию к росту со временем степени корреляции поля скорости и магнитного поля, существуют отдельные реализации с сильно отличающейся скоростью диссипации энергии пульсаций. В настоящей работе пробле-

ма исследуется с упором на выявление роли спи-ральностей в формировании конкретного сценария эволюции системы.

Результаты численного моделирования

Эффективным инструментом исследования статистических свойств развитой мелкомасштабной турбулентности зарекомендовали себя каскадные модели (см., например, [2]), описывающие процессы спектрального переноса в развитой турбулентности с помощью небольшого числа переменных, каждая из которых является коллективной характеристикой амплитуд пульсаций поля скорости и магнитного поля в некотором интервале волновых чисел. Уравнения для коллективных переменных записываются так, чтобы они воспроизводили «базовые» свойства исходных уравнений движения, а именно имели те же интегралы движения и тот же вид нелинейности. В настоящем исследовании использовалась каскадная модель, предложенная в [3], отличающаяся более адекватным способом описания магнитной спиральности.

Уравнения каскадной модели численно интегрировались при гидродинамическом и магнитном числах Рейнольдса, равных 105. В начальный момент времени энергия была сконцентрирована

на самом крупном масштабе, а затем, в процессе каскада по спектру, достигала диссипативного масштаба. Во всех реализациях начальные состояния отличались на некоторую случайную величину, так что Eu ~ Eb ~ 1, \Hc\ < 10-4, \Hb\ < 10-4.

Эволюция системы отслеживалась по изменению полной энергии и нормированной перекрестной спиральности C = Hc/E, которые показаны соответственно на рис. 1а, б. Выделяются реализации, где C не достигает предельных значений.

E

10-1

10-2

10-3

10-4

10-5

1

101 102 103 104 t

а)

С

-0.5

-1.0 L

чае система стремится к предельному значению нормированной перекрестной спиральности C = = ±1, а во втором - к предельному значению нормированной магнитной спиральности Cb = ±1. При этом наступление первого сценария даже при начально нулевой перекрестной спиральности наиболее вероятно за счет ее генерации вблизи диссипативного масштаба и последующего обратного каскада. Реализация второго сценария связана с законом сохранения магнитной спиральности, которой не свойственен прямой каскад. Магнитная спиральность, которая имеется изначально или возникает за счет диссипации, приводит к задержке части магнитной энергии.

102 103

б)

Рис. 1

Распределение состояний по значениям нормированной перекрестной C и магнитной спи-ральности Cb = Hb/Eb , а также по уровням магнитной и кинетической энергии показано на рис. 2. Выясняется, что второе предельное состояние характеризуется максимальным уровнем магнитной спиральности. На рисунках обозначено: крестики - t = 100, кружочки - t = 1000, точки - t = 104. Также следует отметить, что кинетическая энергия по-прежнему способна переноситься по спектру и диссипировать.

Можно выделить два сценария свободного вырождения МГД-турбулентности: в первом слу-

Рис. 2

В настоящем исследовании также принимал участие П. Г. Фрик.

Список литературы

1. Антонов Т.Ю., Фрик П.Г, Соколов Д.Д. // Вычислительные методы и программирование. 2000. Т. 1, №1. С. 14-18.

2. Фрик П.Г. Турбулентность: подходы и модели. М.-Ижевск: Регулярная и хаотическая динамика, 2010. 332 с.

3. Мизева И.А., Степанов Р.А., Фрик П.Г. // Докл. РАН. 2009. Т 424, №4. С. 479-483.

THE EFFECT OF MAGNETIC AND CROSS HELICITIES ON THE FREE DECAYING MHD TURBULENCE

R.A. Stepanov

Free decay of MHD turbulence for large Reynolds numbers is studied numerically using a shell model. The statistical properties are studied based on a representative sample of realizations (128 realizations for each type of initial conditions) over the period of 5 times the order of magnitude of large-scale turbulence turnover times. The performed simulations show that the free-decaying non-helical MHD turbulence can demonstrate two different scenarios of evolution in spite of similar initial conditions. According to the first scenario, the cross-helicity accumulation is so fast that the energy cascade vanishes before

any significant magnetic energy dissipates. Then the system approaches the state of maximal cross-helicity. In the second scenario, the cascade process remains active until the time 10000 in the units of the large-scale turnover time. Then the magnetic field becomes vastly helical due to magnetic helicity conservation. Thus, the magnetic energy does not dissipate with kinetic energy.

Keywords: magnetohydrodynamics, theory of turbulence, free decay, shell models, conservation laws.