Задачи с малым параметром и распространение фронтов в теории галактического динамо
Е. А. Михайлов
Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра математики. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
E-mail: [email protected]
Статья поступила 07.11.2014, подписана в печать 28.11.2014.
Исследован вопрос о распространении магнитного поля на периферию галактического диска и возможности его переноса с помощью динамо-волны. Для этого рассмотрена система из двух уравнений с малым параметром при старших производных. Ее существенной особенностью является то, что при определенном значении независимой переменной происходит смена устойчивости корней вырожденного уравнения. Для данной задачи численно получены скорости распространения переходных слоев, исследован вопрос о возможности их проникновения за пределы точки смены устойчивости. Кроме того, изучено близкое по смыслу скалярное уравнение, для которого могут быть получены асимптотические оценки скорости распространения фронта, хорошо согласующиеся с результатами численного моделирования.
Ключевые слова: магнитные поля галактик, динамо-волна, задачи с малым параметром.
УДК: 517.958. PACS: 02.30.Jr.
Введение
Вопрос о существовании магнитных полей на больших расстояниях от центра галактики является весьма важным с астрофизической точки зрения. Твердо установлено [1], что на расстояниях до 10-15 кпк от центра галактики присутствуют поля напряженностью порядка 1 мкГс. Тем не менее было бы неправильным считать, что галактический диск имеет жесткую границу: плотность межзвездного газа, движения которого и определяют генерацию магнитных полей [2], достаточно плавно спадает по мере удаления от центра, поэтому логично предположить, что поле тоже проникает достаточно далеко и превышает значение межгалактического фона, составляющего порядка 10~9 мкГс [3]. Генерация поля определяется механизмом динамо, работа которого связана с совокупным действием дифференциального вращения и альфа-эффекта. Этот механизм имеет пороговый характер и работает лишь при определенных значениях параметров, описывающих межзвездную среду [2]. Как правило, в центральных областях галактик этот порог превышен, а на периферии все обстоит несколько хуже: действие альфа-эффекта серьезно ослаблено. Однако следует учитывать другую возможность, когда поле генерируется в центральных частях галактики и там возникает волна, выносящая его во внешние области. Кроме того, представляет интерес ситуация, когда в центральных частях галактики образуются два концентрических кольца с противоположными направлениями магнитного поля [4]. Подобные явления качественно исследовались в [5], некоторые оценки получены также в [6]. Существование таких волн для конкретных параметров межзвездной среды Млечного Пути подтверждается результатами моделирования [7]. Тем не менее многие детали данного процесса в случае больших расстояний от центра галактики, где значения параметров находятся на границе или даже ниже порога генерации магнитного поля, остаются неясными.
С другой стороны, в математике хорошо известны процессы нелинейного распространения волн [8-11]. Часто они возникают в нелинейных параболических уравнениях, содержащих малый параметр при лапласиане. Такие уравнения встречаются в различных задачах физики, химии и биологии. То же самое можно сказать и про уравнения магнитной гидродинамики, описывающие эволюцию космических магнитных полей. В роли малого параметра там выступает коэффициент турбулентной диффузии, характеризующий диссипацию крупномасштабного магнитного поля. В ряде случаев для подобных уравнений могут быть получены асимптотические разложения, помогающие понять природу таких явлений, как распространение фронтов в галактиках. Для более сложных задач результаты могут быть уточнены численно.
Мы рассматриваем модельную задачу об эволюции магнитного поля в галактике, численно решая систему из двух уравнений с малым параметром (каждая из неизвестных функций характеризует одну из компонент магнитного поля). Данная задача характерна тем, что в определенной точке происходит смена устойчивости: устойчивые корни становятся неустойчивыми и уходят в комплексную область, а неустойчивый корень становится устойчивым. Моделирование показывает, что фронт может проникнуть на небольшое расстояние за пределы точки смены устойчивости.
Далее рассматривается упрощенная скалярная задача, для которой могут быть получены асимптотические скорости распространения волны. Результаты сравниваются с численным моделированием, проведенным для разных значений малого параметра, входящего в задачу. Показано, что результаты, полученные асимптотически и численно, совпадают с очень большой точностью.
1. Уравнения галактического динамо
Если учесть соленоидальность магнитного поля и ввести безразмерные переменные (время измерять в расстояния в радиусах галактики R), то основ-
ные уравнения динамо для компонент Вг и Вф в рамках планарного приближения примут вид [13, 14]
Решения (ф, ф) выглядят следующим образом:
dB = -Ra (1 - Jü) Вф - B + а2ab;
dB, ~dt
an
dr
(1)
v = Ry^Br - П-Вф + Х2АгВф.
Параметр Ra характеризует действие альфа-эффекта, Ru — дифференциальное вращение, А = h/R — полутолщину диска, В* — поле насыщения, определяемое плотностью энергии турбулентных движений. Как правило, полагают, что Ra = 0(1), Ru = 0(10), А = 0(10~2). Таким образом, в ряде случаев А может выступать как малый параметр, описывающий возникновение контрастных структур и движение переходных слоев в галактике.
Возможность роста магнитного поля определяется в первую очередь соотношением между коэффициентами Ra и Ru .В том случае если RaRu > 7, в линейной области В с В* происходит устойчивая генерация поля за времена, сопоставимые с временем жизни галактики. В противном случае поле затухает либо его рост крайне медленный. Эти коэффициенты меняются в зависимости от расстояния до центра галактики: в центральных областях, как правило, их произведение превышает порог генерации, за счет чего возникает магнитное поле, наличие которого подтверждается наблюдениями [1]. Однако во внешних областях галактики коэффициенты имеют намного меньшие значения. Таким образом, рост магнитного поля в этих областях навряд ли возможен. Тем не менее результаты численного моделирования [7] показывают, что в центральной области может сформироваться фронт, который движется в сторону возрастания r и переносит магнитное поле, сгенерированное в центральных частях галактики, в ее внешние части. Вопрос о формировании и распространении таких фронтов является основным для настоящего исследования.
2. Математическая задача и анализ ее численных решений
Путем преобразований, приводящих задачу (1) к более удобному виду, и небольших упрощений, связанных с тем, что мы исследуем случай больших расстояний от центра галактики, систему уравнений (1) можно привести к виду
ut = -A(x)v(1-
u2 + v2 U 2(x)
)-U+S2Uxx, Vt = -B(x)u-V+£2Vx
(2)
с условиями
u(0, t) = u0, v(0, t) = v0, u(x, t) ^ 0, v(x, t) ^ 0 при x ^ TO, u(x, 0) = f (x), v(x, 0) = g(x).
(3)
Из физических соображений следует [7], что функции А(х), В(х) и и(х) монотонно затухают с ростом х. Положим е = 0 и рассмотрим стационарные решения задачи (2) (| = 0):
0 = -A(x)v(1 -
u2 + v2 U 2(x)
) - u, 0 = —B(x)u - v.
ф1 (x) = U (x)
■I
A(x)B(x) - 1 A(x)B(x)(1 + B2(x))'
фl(x) = -U (x)
ф2(x) = 0, ф2(x)= 0,
'f3(x) = -U (x)
/
A(x)(1 + B2(x)) '
■I
A(x)B(x) - 1 A(x)B(x)(1 + B2(x))'
ФзМ = U (x)
/
B(x)(A(x)B(x) - 1) A(x)(1 + B2(x)) ■
При АВ > 1 решения 1 и 3, подставленные в (2), устойчивы, нулевое решение 2 неустойчиво. Физически это означает, что любое небольшое магнитное поле будет усиливаться механизмом динамо вплоть до достижения определенного предельного значения. В тех областях, где АВ < 1, ситуация меняется: устойчивым становится только нулевое решение 2, а решения 1 и 3 уходят в комплексную область и не имеют физического смысла.
Результаты численного моделирования показывают (рис. 1), что при определенных начальных условиях в области, где АВ > 1, может образоваться переходный слой, соединяющий два решения (ф1, ф1) и (ф3, ф3). Данный слой может достаточно медленно двигаться в сторону возрастания х, что с астрофизической точки зрения соответствует распространению динамо-волны во внешние области галактики. Возникает вопрос о том, сможет ли данный фронт распространиться за пределы точки смены устойчивости хс^. С другой стороны, отдельный интерес представляет вычисление скорости распространения данного фронта. Эта задача решалась численно. Рассматривались следующие значения входящих в задачу функций:
А(х) = ехр(-0.3466х), В(х) = 2 ехр(-0.3466х), и (х) = ехр(-х).
и, V
Рис. 1. Типичное поведение решения задачи (2): е = 0.02, t = 50. Сплошная линия показывает и, пунктирная — V
Значение малого параметра е варьировалось, принимая в разных случаях значения порядка 10~2. При таких начальных условиях точка смены устойчивости xch ~ 1. Для решения задачи использовались следующие начальные и граничные условия:
и(0, ^ = v(0, 0, u(x, 0) = 0,
-0.05 при 0.1 < x <0.2, v(x, 0) = { 0.05 при 0.3 < x <0.4, ,0 при остальных x.
(4)
На начальном этапе происходит рост функций u и v в областях (0.1,0.2) и (0.3,0.4), их значения там становятся близкими к стационарным. Затем формируется переходный слой между двумя стационарными решениями, локализованный на начальном этапе в окрестности точки Xfr = 0.25. После этого фронт начинает смещаться влево. Мы вычислили скорости распространения фронта для различных е, полученные значения проиллюстрированы в таблице. С большой точностью зависимость скорости распространения фронта w от малого параметра е описывается квадратичным законом: w(е) « 31е2. Был изучен также вопрос о возможности распространения фронта за пределы точки смены устойчивости x > х^. Численный расчет показывает, что переход фронта в данную область возможен (рис. 2), а глубина проникновения растет с увеличением е.
Скорость распространения переходного слоя между устойчивыми корнями т и переходного слоя между устойчивым и неустойчивым корнем ш
Задача Двумерная система (2) Скалярное уравнение (6)
е да ш w ш
0.01 2.9 • 10~4 1.40 • 10~2 5.0 • 10~5 1.46 • 10~2
0.02 1.16 • 10~3 2.7 • 10~2 2.0 • 10~4 2.9 • 10~2
0.03 2.7 • 10~3 4.0 • 10~2 4.5 • 10~4 4.3 • 10~2
0.04 4.9 • 10~3 5.3 • 10~2 7.8 • 10~4 5.7 • 10~2
0.05 8.0 • 10~3 6.7 • 10~2 1.18 • 10~3 7.1 • 10~2
-0.6
Рис. 2. Распространение фронта, соединяющего устойчивые корни, за пределы точки смены устойчивости при е = 0.03. Сплошная линия показывает момент времени t = 160, пунктирная — t = 240, штриховая — t = 320, штрихпунктирная — t = 400
Другая важная для приложений задача связана с переходным слоем, соединяющим устойчивое при x < х^ решение (ф3(х), ф3(х)) и неустойчивое нулевое решение. В данном случае интерес представляют те же вопросы: скорость ш распространения переходного слоя вправо и возможность его проникновения в область x > х^. Для этого численно решалась задача (2) с условиями вида
и(0, ^ = v(0, 0, u(x, 0) = 0,
{0 при x <0.1,
0.05 при 0.1 < x <0.2,
0 при x >0.2.
(5)
На рис. 3 показан процесс распространения фронта для разных моментов времени. Отметим, что переходный слой проникает в правую область на небольшое расстояние, где через определенное время останавливается. Глубина проникновения оказывается тем большей, чем больше е, стабилизировавшиеся решения показаны на рис. 4. Скорости ш распространения фронта в окрестности точки x = 0.5 проиллюстрированы в таблице. С большой точностью зависимость можно аппроксимировать законом ш = 1.34е.
Рис. 3. Распространение фронта, соединяющего устойчивый и неустойчивый корни, при е = 0.03. Сплошная линия показывает момент времени t = 10, пунктирная — t = 20, штриховая — t = 50
Рис. 4. Конечная конфигурация поля при начальных условиях вида (5). Сплошная линия — е = 0.01, пунктирная — е = 0.03, штриховая — е = 0.05
3. Скалярная задача и построение асимптотики
Построение асимптотики и вычисление скорости распространения фронта в двумерной задаче (2) представляет определенные сложности. В то же время качественно поведение магнитного поля может быть описано и более простым скалярным уравнением
щ = ^(х){ 1 — и2) и + е2
(6)
С помощью перехода к другому масштабу времени ^ = t/е2 задачу можно свести к стандартному виду
е2Щ' = е2ихх + !(х, и),
I(х, и) = ^(х)(1 — и2)и
(7)
с граничными условиями
и(0, ?) = 0, и(х, ^) ^ 0 при х ^ то,
где ^(х) монотонно убывает с ростом х (возможно, принимая также отрицательные значения). Введем растянутую переменную т = , где х* — точка локализации погранслоя. Решения будем искать в виде следующих комбинаций [15]:
и(-) = Щ(-)(х, Г, е) + Q(-)(т, е)+ Я(-) при х < х*,
1(+) = Щ+Цх, ^ е) + <3(+)(т, е)+ Я(+)
(8)
при х > х ,
где Q(± — погранслойные функции, описывающие поведение решения в окрестности точки х*, описывают поведение решения вблизи границ области, а стационарные решения имеют вид
Щ(-) = —1, Щ(+) = 1.
При х < х* :
! = —21(хП(-\ 1 + Q(-У) (2 + Q(-У), а производная решения по времени выглядит так [15]:
и, = ± Q(-Чт ) = Q'(-\т) % = — 1 wQ'(-). Тогда в нулевом порядке получаем уравнение
Q"(-)(т) — 27(х0^(-) (т) (1 + Q(-))(2 + Q(-)) = 0.
Аналогичное уравнение можно составить и для х > х*. Решения для пограничных функций в нулевом порядке выглядят следующим образом:
Q0-) = 1 +№
Ш-
Введем функцию
Q(+^)= — 1 +Ш
1
ш-
Ф(т) = ^ = _, 1. Ф(т) дт V 2 сь2 (туДУ
Неоднородность преобразуется к виду
+Ь (т Д") [Щ(х) + Qo(т), х) = 1(х) 2\ .
сЬ (тУ1)
Скорость движения фронта может быть найдена при помощи интеграла [16]:
/ 1х(щ(х) + Qo(т), х^тФ(т) йт
да0 = —
(9)
] Ф2 (т) йт
При вычислении интеграла (9) получается результат w0 = — 27'(х0)/7(х0). При возврате к «физическому» времени t получим в нулевом порядке скорость
2
распространения фронта, равную да0 = — ^'(х0)/ч(х0). Проверим полученные результаты численно, рассмотрев задачу (6) с коэффициентом х) = ехр(—х) и начальными условиями вида
Г 0.05 при 0.1 < х <0.2, и(х,0) = Г — 0.05 при 0.3 < х <0.4, (10) 0 при остальных х.
Значение малого параметра е принимало различные величины порядка 10-2. Численные результаты показаны в таблице и соответствуют аналитической формуле,
2
которая при таком выборе ^(х) даст да0 = у. Таким образом, характер зависимости (скорость квадратична по е) повторяет то, что было получено в предыдущем разделе для двумерной задачи.
Другой вопрос связан с распространением фронта, соединяющего устойчивое решение и = 1 и неустойчивое решение и = 0. Скорость распространения фронта для таких уравнений была получена, например, в [5], где утверждается, что ш « 2е^/^(х), т.е. скорость распространения переходного слоя зависит от х. Мы решили также численную задачу с начальными условиями
и(х, 0) =
0.05 при 0.1 < х <0.2, 0
(11)
при остальных х.
Скорость распространения переходного слоя в окрестности точки х = 0.5 согласно аналитической формуле должна составлять величину ш « 1.2е. Результаты численного моделирования (см. таблицу) вполне близки к этому результату.
Заключение
Нами исследована задача распространения фронта магнитного поля во внешние области галактики. Показано, что динамо-волна может переносить магнитное поле в том числе и во внешние области галактики, где невозможна его самостоятельная генерация. Также изучен процесс движения во внешние части галактики переходного слоя между областями, в которых направление поля противоположно. Для простейшего скалярного случая получены асимптотические скорости распространения фронтов, характер зависимости которых от малого параметра е очень близок к тому, что получено в результате компьютерного моделирования как для скалярного, так и для более сложного двумерного случая. Продемонстрировано, что скорость распространения переходного слоя между двумя устойчивыми решениями пропорциональна е2, между устойчивым и неустойчивым — е.
Автор выражает благодарность профессору Д. Д. Соколову и профессору Н. Н. Нефёдову за ценные рекомендации при подготовке работы.
Работа выполнена при частичной финансовой поддержке РФФИ (грант 13-01-00200) и Программы поддержки аспирантов и молодых ученых без степени Фонда некоммерческих программ «Династия».
и
Список литературы
1. Beck R., Brandenburg A., Moss D. et al. // Ann. Rev. Astron. Astrophys. 1996. 34. P. 155.
2. Arshakian T., Beck R., Krause M., Sokoloff D. // Astron. Astrophys. 2009. 494. P. 21.
3. Neronov A, Semikoz D. // Phys. Rev. D. 2009. 80. 123012.
4. Moss D., Stepanov R., Arshakian T. et al. // Astron. Astrophys. 2012. 537. P. 68.
5. Moss D., Shukurov A., Sokoloff D. // Geophys. Astrophys. Fluid Dynam. 1998. 89. P. 285.
6. Sokoloff D., Petrov A, Moss D. // Proc. of the Int. Conf. «Plasma Turbulence and Energetic Particles in Astrophysics». Cracow, Poland, 1999. P. 92.
7. Mikhailov E., Kasparova A., Beck R. et al. // Astron. Astrophys. 2014. 568. A66.
8. Колмогоров А.Н., Петровский И.Г., Пискунов Н.С. // Бюлл. МГУ. Матем. и мех. 1937. 1, № 6. C. 1.
9. Nefedov N., Sakamoto K. // Hiroshima Mathematical J. 2003. 33, N 3. P. 391.
10. Нефёдов Н.Н., Никитин А.Г., Петрова М.А., Рекке Л. // Дифференц. уравн. 2011. 47. С. 1305.
11. Nefedov N.N., Recke L, Schnieder K.R. // J. of Mathem. Analysis and Applications. 2013. 405. P. 90.
12. Phillips A. // Geophys. Astrophys. Fluid. Dynam. 2001. 94. P. 135.
13. Moss D. // Mon. Not. R. Astr. Soc. 1995. 275. P. 191.
14. Михайлов Е.А. // Письма в Астрон. журн. 2013. 39. С. 474.
15. Божевольнов Ю.В., Нефёдов Н.Н. // Журн. вычисл. матем. и матем. физ. 2010. 50, № 2. С. 276.
16. Volkov V.T., Nefedov N.N. Lecture Notes in Computer Science. B., 2013. 8236. P. 524.
Problems with a small parameter and propagation of fronts in the galactic dynamo theory E. A. Mikhailov
Department of Mathematics, Faculty of Physics, M. V. Lomonosov Moscow State University,
Moscow 119991, Russia.
E-mail: [email protected].
Propagation of the magnetic field to the periphery of the galactic disc and the possibility of its transfer using a dynamo wave are studied. For this purpose, a system of two equations with a small parameter of high-order derivatives is considered. Its important feature is that for a certain value of the independent variable the stability of the roots of the degenerate equation changes. For this problem the numerical velocity values of transition layers are obtained. The problem of their possible penetration beyond the bounds of the point of stability change is studied. A similar scalar equation is studied, for which asymptotic estimates of the front propagation velocity can be obtained. These estimates are in good agreement with the numerical simulation results.
Keywords: magnetic fields of galaxies, dynamo wave, problems with a small parameter.
PACS: 02.30.Jr.
Received 7 November 2014.
English version: Moscow University Physics Bulletin 2(2015).
Сведения об авторе
Михайлов Евгений Александрович — аспирант; тел.: (495) 939-10-33, e-mail: [email protected].