close

Вход

Забыли?

вход по аккаунту

?

Инверсии в модели геодинамо управляемой 6-ячейковой конвекцией.

код для вставкиСкачать
Вестник КРАУНЦ. Физ.-мат. науки. 2015. № 2(11). C. 45-54. ISSN 2079-6641
DOI: 10.18454/2079-6641-2015-11-2-45-54
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
УДК 517.9
ИНВЕРСИИ В МОДЕЛИ ГЕОДИНАМО, УПРАВЛЯЕМОЙ 6-ЯЧЕЙКОВОЙ
КОНВЕКЦИЕЙ
Г.М. Водинчар1, 2 , Л.К. Фещенко1
1
Институт космофизических исследований и распространения радиоволн ДВО РАН,
684034, Камчатский край, п. Паратунка, ул. Мирная, 7
2 Камчатский государственный университет имени Витуса Беринга, 683032,
г. Петропавловск-Камчатский, ул. Пограничная, 4
E-mail: gvodinchar@gmail.com
Описана крупномасштабная модель геодинамо, основанная на косвенных данных о неоднородности плотности в жидком ядре Земли. Конвективная структура соотнесена со
сферической гармоникой Y42 , определяющей основную полоидальную составляющую скорости. Кориолисов снос этой моды определяет тороидальную составляющую скорости.
Это формирует 6-ячейковую конвективную структуру. В модели учитывается эффект
обратного влияния магнитного поля на конвекцию. Установлено, что в модели реализуется устойчивый режим генерации поля. Скорость конвекции и величина поля в этом
режиме согласуются с реальными.
Ключевые слова: геодинамо, инверсии геомагнитного поля, конвекция
© Г.М. Водинчар, Л.К. Фещенко, 2015
MSC 97A80
MATHEMATICAL MODELING
REVERSALS IN THE 6-CELLS CONVECTION DRIVEN
G.M. Vodinchar1, 2 , L.K. Feshenko2
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: gvodinchar@gmail.com
We describe the large-scale model geodynamo, which based on indirect data of inhomogeneities
in the density of the Earth’s core. Convection structure is associated with spherical harmonic
Y42 , which defines the basic poloidal component of velocity. Coriolis drift of this mode
determines the toroidal component of velocity. Thus, 6 convective cells are formed. The
model takes into account the feedback effect of the magnetic field on convection.It was
ascertained that the model contains stable regimes of field generation. The velocity of
convection and the dipole component of the magnetic field are close to the observed ones.
Key words: geodynamo, geomagnetic field reversals, convection
© Vodinchar G.M., Feshenko L.K., 2015
45
ISSN 2079-6641
Водинчар Г.М., Фещенко Л.К.
Введение
Процесс формирования магнитных полей планет и звезд успешно объясняется
теорией гидродмагнитного динамо. Существующие модели конвекции в жидких ядрах планет и конвективных зонах звезд допускают структуры течений, способные
генерировать магнитные поля, подобные реальным [1, 2, 3].
Современные вычислительные системы не позволяют вести прямое численной моделирование на разностных сетках нестационарных трехмерных задач планетарного
динамо на временных масштабах порядка времени существования поля (∼ 109 лет).
Поэтому существующие численные модели либо воспроизводят МГД-течения с хорошим разрешением по пространству на относительно небольших (по геологическим
масштабам) промежутках времени (∼ 104 лет), либо воспроизводят длительную эволюцию наиболее крупномасштабных пространственных структур. В моделях первого
типа пространственная структура течений и поля просчитывается, а для моделей второго типа эта структура должны изначально задаваться. При этом возникает вопрос
о том, какова реальная крупномасштабная структура конвекции.
Косвенные данные об этой структуре для ядра Земли можно получить из данных
о пространственной неоднородности плотности жидкого ядра. В работе [4] проанализированы результаты исследований функций расщепления собственных колебаний
Земли после сильнейшего Чилийского землетрясения. Распределение этих функций
характеризует неоднородность ядра Земли на различных глубинах. Проанализировав
эти функции, автор [4] пришел к выводу, что в жидком ядре существует 12-зонная
шахматного типа неоднородность, которая естественным образом описывается в первом приближении сферической гармоникой Y42 . Им же была выдвинута гипотеза о
том, что в ядре существует конвекция с 6 областями подъема вещества и 6 областями
опускания, т.е. о 6-ячейковой конвекции. В такой конвективной структуре в каждой
из полярных областей расположены по 2 ячейки и еще 2 лежат в экваториальной
области.
Отметим также, что 6-ячейковая структура возникает в результате прямого численного моделирования конвекции в ядре Земли при некоторых значениях параметров ядра [5]. В этих расчетах в каждом полушарии расположены по 3 ячейки. Такая
структура описывается сферической гармоникой Y43 .
Известным свойством реальных космических динамо-систем является наличие
инверсий – смены знака дипольной составляющей поля. При этом большой интерес
представляют инверсии не связанные с перестройкой конвекции.
В настоящей работе исследуется вопрос о генерации поля 6-ячейковой структурой, описываемой гармоникой Y42 с учетом турбулентного α-эффекта, и об инверсиях
в такой системе.
Основные уравнения
Рассматриваем сферическую оболочку вязкой несжимаемой жидкости (жидкое
ядре), вращающуюся вокруг z-оси с постоянной угловой скоростью . В сферической
системе координат (r, θ , ϕ) внутрення граница r = ri и внешняя граница r = ro . Температура на внутренней и внешней границах постоянна и равна, соответственно, Ti
и To .
46
Инверсии в модели геодинамо, управляемой 6-ячейковой . . .
ISSN 2079-6641
Уравненя геодинамо, учитывающие α-эффект в приближении Буссинеска имеют
вид:
E
∂
+ v · ∇ v = E∇2 v − ∇p − 2ez × v + RaPmT r + rotB × B,
Pm ∂t
Pm 2
∂T
+ (v · ∇) (T + Ts ) =
∇ T,
∂t
Pr
∂B
= rot (v × B) + Rα rot (αB) + ∇2 B,
∂t
∇ · v = 0,
(1)
∇ · B = 0.
Основными безразмерными параметрами (параметрами подобия) являются число
Экмана E = ν/Ωro2 , число Релея Ra = β g0 δ Tro /νΩ, число Прандтля Pr = ν/κ, магнитное число Прандтля Pm = ν/η и амплитуда α-эффекта Rα . В этих выражения ν
обозначает кинематическую вязкость, β – коэффициент объемного расширения, g0 –
ускорение свободного падения на внешней границе, δ T = Ti − To – разность значений
температуры на внутренней и внешней границах, κ – коэффициент температуропроводности и η – коэффициент омической диссипации. Такой набор параметров
√
возникает, если ro , ro2 /η, δ T и µ0 ρωη, где ρ – плотность и µ – магнитная проницаемость, являются единицей длины, единицей времени, единицей температуры и
единицей магнитного поля, соответственно. Радиус внутреннего ядра, выраженный
в такой единице длины, составляет ri = 0.35.
В уравнениях (1) переменная T обозначает отклонение температуры от равновес1/r − 1
+ Ti − 1.
ного гиперболического r-профиля Ts =
1/ri − 1
Мы считаем, что турбулентность в ядре изотропная и используем скалярную
параметризацию α-эффекта в виде α = a (r) cos θ , где max |a(r)| = 1.
Для скорости v принимаются граничные условия прилипания. Также мы считаем,
что магнитная проницаемость внутреннего и внешнего ядра одинаковы и среда вне
ядра (r > ro ) непроводящая. Таким образом, для магнитного поля принимаются вакуумные граничные условия на внешней границе и условия ограниченности в центре
Земли. Наконец, отклонения температуры T на границах внешнего ядра считаются
нулевыми.
Далее, мы будем использовать специальные представления для скорости, температуры и магнитного поля, основанные на некоторых спектральных разложениях.
Используя эти представления, методом Галеркина получим уравнения нашей модели.
Конвективная часть модели
Исходим из следующего основного требования – гидродинамическая часть нашей
модели должна совпасть с классической системой Лоренца для конвекции в плоском
слое. Это совпадение должно быть не только формальным, но и соответствующим
самой идее построения системе Лоренца.
Напомним, что при выводе системы Лоренца используют одномодовое приближение для скорости и двухмодовое приближение для температуры в задаче свободной
47
ISSN 2079-6641
Водинчар Г.М., Фещенко Л.К.
конвекции в плоском слое [6]. В качестве моды скорости используют одну из собственных мод свободных колебаний жидкости в слое. Вертикальная компонента этой
моды не имеет нулей между границами слоя, т.е. обеспечивает перенос жидкости от
нижней границы к верхней. В качестве температурных мод используются две собственные моды свободных затухающих колебаний отклонений температуры от стационарного профиля. Одна из них пространственно согласована с модой скорости,
другая изменяется только в вертикальном направлении.
В нашей модели мы будем использовать подобные моды. Важно, что при любом
одномодовом приближении скорости галеркинская процедура удаляет из уравнения
Навье-Стокса кориолисов член. Поэтому необходимо использовать моды, сама пространственная структура которых содержит информацию о вращении слоя. Одним
из возможных вариантов является использование собственных мод свободных колебаний вязкой вращающейся жидкости, т.е. решения спектральной задачи
E
µv + E∇2 v − 2ez × v − ∇p = 0,
Pm
∇ · v = 0,
(2)
v(r = ri ) = v(r = ro ) = 0.
Для удобства мы использовали те же безразмерные параметры, что и в уравнениях (1).
Невязкий аналог этой задачи известен как задача Пуанкаре и является косоэрмитовым [7]. Далее мы будем называть решения задачи (2) модами Пуанкаре. В вязком
случае оператор задачи не обладает эрмитовой или косоэрмитовой симметрией, что
создает значительные трудности. Точные решения в настоящее время не известны.
Есть только некоторые оценки спектра, доказана полнота системы собственных и
присоединенных полей этой задачи [8]
Мы будем использовать простую аппроксимацию решений задачи (2), основанную на хорошо известных модах свободных колебаний невращающейся оболочки.
Формально они определяются как решения спектральной задачи
E
µv + ∇2 v − ∇p = 0, ∇ · v = 0,
Pm
v(r = ri ) = v(r = ro ) = 0.
(3)
Эта задача является эрмитовой, ее решения образуют полную ортогональную
систему и распадаются на классы тороидальных и полоидальных мод.
Тороидальные и полоидальные моды, соответственно, имеют следующий вид:
1
∂
∂
T
T
vk,n,m = Rk,n (r)
eθ
− eϕ
Ynm (θ , ϕ),
sin θ ∂ ϕ
∂θ
P
Rk,n (r) m
vPk,n,m = n(n + 1)
Yn (θ , ϕ)er +
(4)
r
d
1
∂
1
∂
+
+
RPk,n (r) eθ
+
eϕ
Y m (θ , ϕ),
dr r
∂ θ sin θ ∂ ϕ n
k = 0, 1, . . . , n = 1, 2, . . . , m = −n, . . . , n,
где RTk,n (r) и RPk,n (r) некоторые комбинации степенных функций и сферических функций Бесселя, а Ynm (θ , ϕ) – сферические гармоники.
48
Инверсии в модели геодинамо, управляемой 6-ячейковой . . .
ISSN 2079-6641
Следующие множества соленоидальных полей, равных нулю на границах, образуют разложение пространства соленоидальных полей в прямую сумму ортогональных
подпространств:
H0T = vTk1 ,1,0 , vPk2 ,2,0 , vTk3 ,3,0 , vPk4 ,4,0 , . . . ,
H0P = vPk1 ,1,0 , vTk2 ,2,0 , vPk3 ,3,0 , vTk4 ,4,0 , . . . ,
HmT = vTk0 ,m,±m , vPk1 ,m+1,±m , vTk2 ,m+2,±m , vPk3 ,m+3,±m , . . . ,
HmP = vPk0 ,m,±m , vTk1 ,m+1,±m , vPk2 ,m+2,±m , vTk3 ,m+3,±m , . . . ,
(5)
ki = 0, 1, 2, . . . , m = 1, 2, 3, . . .
Каждое из этих подпространств инвариантно относительно оператора задачи (2).
Это легко проверить, если вычислить матрицу оператора в базисе (4). Поэтому каждая из мод Пуанкаре полностью лежит в одном из подпространств (5).
Для представления 6-ячейковой конвекции важно подпространство H2P , содержащее полоидальную моду vP0,4,±2 . Из уравнений (4) хорошо видно, что радиальная
компонента этой моды в (θ , ϕ)-направлении определяется сферической гармоникой
Y4±2 .
В нашей модели мы используем только одну моду скорости v0 , являющуюся простейшей аппроксимацией мод Пуанкаре из подпространства H2P . Мы определим v0 в
виде
v0 = β1 vT1,3,−2 + β2 vT1,3,2 + β3 vP0,4,−2 + β4 vP0,4,2 + β5 vT0,5,−2 + β6 vT0,5,2 ,
(6)
где коэффициенты βi находятся из задачи (2) методом Галеркина.
Ясно, что радиальная компонента v0 имеет в (θ , ϕ)-направлении структуру гармоники Y4±2 .
Теперь рассмотрим собственные моды тепловой диффузии во внешнем ядре. Они
задаются в виде Tk,n,m = Xk,n (r)Ynm , где Xk,n (r) являются линейными комбинациями
сферических функций Бесселя первого и второго родов.
В нашей модели главная температурная мода T0 должна быть согласована с радиальной компонентой скорости, Поэтом определим ее в выражением
T0 = X0,4 (r) β3Y4−2 + β3Y42 = β3 T0,4,−2 + β4 T0,4,2 ,
(7)
где β3 и β4 взяты из формулы (6).
Следующая температурная мода T1 выбирается в виде T1 = T1,0,0 = X1,0Y00 подобном
моде из системы Лоренца.
Окончательно, мы получаем следующее представление для скорости и температуры в нашей модели
v(r,t) = u(t)v0 (r),
T (r,t) = θ0 (t)T0 (r) + θ1 (t)T1 (r).
49
(8)
ISSN 2079-6641
Водинчар Г.М., Фещенко Л.К.
Магнитная часть модели
Для представления магнитного поля мы используем некоторые моды омической
диссипации, т.е. решения спектральной задачи
ηB + ∇2 B = 0,
B(r = 0) 6= ∞,
∇ · B = 0,
BT (r ≥ 1) = 0,
rotBP (r ≥ 1) = 0,
(9)
где BT и BP тороидальная и полоидальная составляющие магнитного поля, соответственно.
Тороидальные и полидальные решения этой задачи
q
m
T
m
T
ηkn r Yn (θ , ϕ)r ,
Tkn = akn rot jn
q
(10)
m
P
m
P
ηkn r Yn (θ , ϕ)r ,
Pkn = akn rotrot jn
T and
где jn (·) сферические функции Бесселя первого рода. Собственные значения ηkn
P определяются из граничных условий, а aT и aP являются нормирующими коηkn
kn
kn
эффициентами. В дальнейшем считаем эти коэффициенты выбранными так, чтобы
среднеквадратическое значение мод было равно единице. Все эти моды в совокупности образуют ортогональную систему.
Для отбора магнитных мод мы использовали следующую схему, предложенную в
работе [9].
Прежде всего, мы упорядочим моды по возрастанию собственных значений, т.е.
по скорости диссипации. Получим следующую последовательность, где прямоугольниками выделены дублеты – пары мод с одинаковыми собственными значениями:
−2..2
−3..3
−3..3 −4..4
−2..2
−5..5
−1..1
, T−2..2
, P−1..1
, T−1..1
, T−4..4
,
, T−1..1
P01
01 , P02
02 , P03
11 , T03 , P04
11 , P12
04 , P05
....
Затем, выберем несколько мод с наименьшими собственными значениями, обозначив их Bk (r).
Теперь мы запишем магнитное поле в виде B(r,t) = ∑ gk (t)Bk (r) и подставим
k
это выражение в уравнение индукции (трете из уравнений (1)). Скорость задаем в
виде (8), где временно считаем, что u(t) = u0 = const > 0. Моду скорости v0 считаем
нормированной так, чтобы ее среднеквадратическое значение равнялось единице.
Тогда u0 можно интерпретировать как магнитное число Рейнольдса Rem .
На следующем шаге применим к уравнению индукции метод Галеркина и получим
систему
dgk
= Rem ∑ Wki gi + Rα ∑ Aki gi − ηi gi ,
dt
i
i
Z
Wki =
rot (v0 × Bi ) Bk dV,
(11)
Z
Aki =
rot (a(r) cos θ Bi ) Bk dV,
где ηi – это собственное значение для Bi .
Фактически, мы получим маломодовую аппроксимацию для задачи кинематического динамо с шестью конвективными ячейками. Пусть λi обозначает собственное значение матрицы системы (11). Кинематическое динамо генерирует поле, если
50
Инверсии в модели геодинамо, управляемой 6-ячейковой . . .
ISSN 2079-6641
max Reλi > 0. Мода с таким собственным значением является лидирующей, растущей быстрее других. Если собственное значение лидирующей моды комплексное, то
магнитное поле будет осциллировать.
Далее, мы будем постепенно наращивать число магнитных мод, пока не получим
осциллирующее динамо.
Такова схема отбора магнитных мод.
Разумеется собственные значения зависят от параметров Rem и Rα , а также от
формы a(r) радиального профиля α-эффекта. Мы провели серию расчетов, используя
два выражения для профиля a(r) = 1 и a(r) = r. Оказалось, что результаты отличаются незначительно.В этих расчетах параметры Rem и Rα равномерно варьировались
в логарифмической шкале в диапазоне [10−1 ; 103 ]. После получения осциллирующего
динамо мы отбрасывали некоторые моды, если это не меняло качественно результата
или не меняло существенно порог генерации.
±2
±2
0
В результате мы остановились на магнитных модах P001 , P±2
03 , P11 , T04 , P05 .
Полученная схема генерации магнитных мод крупномасштабной конвекцией и
α-эффектом показана на рис. 1, а области осциллирующего и неосциллирующего
динамо на плоскости параметров (Rem , Rα ) показаны на рис. 2.
Рис. 1. Схема динамо. Сплошные стрелки – крупномасштабный генератор; пунктир-
ные стрелки – α-генератор.
Рис. 2. Области генерации поля. Красные точки – неосциллирующее динамо; зеле-
ные точки – осциллирующее динамо. a(r) = 1 (вверху) и a(r) = r (внизу).
51
ISSN 2079-6641
Водинчар Г.М., Фещенко Л.К.
На этих рисунках видно, что для достаточно больших Rα возможна генерация
только за счет α-эффекта, но при этом не происходит генерации дипольной компоненты. В свою очередь, при достаточно больших Rem возможна генерация поля и без
α-эффекта, но при этом не будет инверсий.
Нелинейная модель МГД-конвекции
В вышеописанном кинематическом динамо возможен только неограниченный рост
поля. Для получения устойчивой генерации необходимо водить механизм подавления. В связи с этим обратим внимание на области на рис. 2 выделенные эллипсами.
Зафиксируем Rα ∼ 15. Если точка (Rem , Rα ) изменяется вдоль большой оси эллипса,
то значение Rem также варьируется. При этом в крайних положениях динамо не
работает. Так можно получить генерацию поля ограниченной величины.
Напомним, что в нашей конструкции Rem является фиксированным уровнем амплитуды u(t) моды скорости v0 . Поэтому мы можем составить теперь комбинированную модель МГД-конвекции с переменной амплитудой скорости. Для ее получения
подставим выражения для скорости, температуры и магнитного поля в первые три
уравнения (1) и применим метод Галеркина.
Получим нелинейную динамическую систему, описывающую 6-ячейковую МГДконвекция в ядре Земли:
8
E du(t)
= −Eµu(t) + RaPmSθ0 (t) + ∑ Li j gi (t)g j (t),
Pm dt
i, j=1
dθ0 (t)
Pm
= −Fu(t)θ1 (t) + Hu(t) −
ζ0 θ0 (t),
dt
Pr
Pm
dθ1 (t)
= Fu(t)θ0 (t) −
ζ1 θ1 (t),
dt
Pr
(12)
8
8
dgk (t)
= u(t) ∑ Wki gi (t) + Rα ∑ Aki gi (t) − ηi gi (t),
dt
i=1
i=1
k = 1, . . . , 8,
где ζi собственные значеия температурных мод Ti , а µ > 0, S, Li j , F, H являются
коэффициентами Галеркина.
Результаты вычисленных экспериментов с моделью
Опишем теперь некоторые результаты вычислительных экспериментов с моделью (12). Мы брали турбулентные значения для диссипативных коэффициентов
ν = 102 м2 /с, κ = 10−2 м2 /с, η = 20 м2 /с [10]. Внешний радиус ядра ro = 3480
км а угловая скорость Ω = 7.29 × 10−5 рад/с. Соответствующие значения параметров
модели: E = 10−7 , Pr = 104 , Pm = 5. Также мы использовали Rα ∼ 15, следуя оценкам
из работы [11] и различные значения числа Релея из диапазона Ra = 102 ÷ 104 .
Система уравнений (12) численно решалась с применением стандартных процедур
пакета SciLab.
Были найдены устойчивые режимы динамо с инверсиями в магнитном поле для
Ra ∼ 103 . На рис. 3 показаны графики амплитуд скорости и дипольной части магнитного поля для Ra = 5000.
52
Инверсии в модели геодинамо, управляемой 6-ячейковой . . .
ISSN 2079-6641
5
4
axial dipol, 10-4 T
3
2
1
0
-1
-2
-3
-4
-5
0
0.5
0
0.5
1
1.5
2
1
1.5
2
t, in units of L2/η
0
-2
u(t), 10-4 m/s
-4
-6
-8
-10
-12
-14
t, in units of L2/η
Рис. 3. Амплитуды мод скорости и вертикального диполя для Ra = 5000
На этом рисунке видно, что скорость и магнитное полей осциллируют. Но, если
магнитное поле периодически меняет знак, то смены знака скорости не происходит.
Таким образом, в модели реализуется режим инверсий без перестройки структуры
конвекции.
Разумеется, к конкретным численным значениям, полученным в такой простой
модели следует относиться с осторожностью. Тем не менее отметим, что характерное значение скорости на рис. 3 составляет 5 × 10−4 м/с. Это хорошо согласуется
с известными оценками реальной скорости конвекции ∼ 10−4 м/м [1]. Характерное
значение поля в модели 5 × 10−4 Тл также хорошо согласуется с экстраполяцией
реальной величины геомагнитного поля на границу ядро–мантия [12]. Временной
интервал между инверсиями в модели 3200 лет, что коррелирует со средней величиной интервала полярности для палеомагнитных шкал инверсий [12].
Выводы
Разработана маломодовая модель геодинамо, управляемая 6-ячейковой конвекцией в ядре Земли и учитывающая турбулентный механизм генерации поля. Модель
основана на косвенных данных о крупномасштабной структуре конвекции.
Для описания пространственной структуры течений предложена одномодовая аппроксимация скорости. Эта мода является простейшей аппроксимацией одной из
собственных мод свободных колебаний жидкого ядра.
53
ISSN 2079-6641
Водинчар Г.М., Фещенко Л.К.
Численное моделирование показало, что в модели можно реализовать устойчивый
режим генерации поля, сопровождающийся инверсиями поля. Важно, что при этом
не происходит перестройки структуры конвекции.
Характерные значения скорости и поля в модели совпадают по порядку величины
с известными оценками реальных скорости и поля.
Инверсии поля в модели носят регулярный характер, тогда как реальная последовательность инверсий хаотическая. Однако авторам представляется, что получение
хаотических инверсий в подобной простой модели в принципе невозможно, без введения шумовых возмущений модели при моделировании.
Библиографический список
1. Jones C.A. Convection-driven geodynamo models // Phil. Trans. R. Soc. Lond. A. 2000. Vol. 358.
pp. 873-897.
2. Kono M., Roberts P.H. Recent geodynamo simulations and observations of the field // Reviews of
Geophysics. 2002. Vol. 40. pp. B1-B41.
3. Соколов Д.Д., Степанов Р.А., Фрик П.Г. // Динамо: на пути от астрофизических моделей к
лабораторному эксперименту // Успехи. физ. наук. 2014. Т. 184. № 3. С. 313-335.
4. Kuznetsov V.V. The anisotropy of properties of the Earth’s inner core // Physics-Uspekhi. 1997.
Vol. 40. pp. 951–961.
5. Гореликов А.В., Ряховский А.В., Фокин А.С. Численное исследование некоторых нестационарных
режимов естественной конвекции во вращающемся сферическом слое // Вычислительная механика
сплошной среды. 2012. Т. 5. № 2. С. 184-192.
6. Lorenz E.N. Deterministic nonperiodic flow // J. Atmos. Sci. 1963. Vol. 20. pp. 130–141.
7. Greenspan H.P. The Theory of rotating Fluids. New York: Cambridge University Press, 1968.
8. Резников Е.Л., Розенкноп Л.М. О собственных колебаниях вращающейся вязкой жидкости во
внешнем ядре Земли // Вопросы геодинамики и сейсмологии (Вычислительная сейсмология. Вып.
30). М.: Геос, 1998. С. 121-132.
9. Соколов Д.Д., Нефедов С.Н. Маломодовое приближение в задаче звездного динамо // Вычислительные методы и программирование. 2007. Т. 8. С. 195-204.
10. Решетняк М.Ю. Моделирование процессов динамо в геофизике / Дисс. на соиск. уч. степ. д-ра
физ.-мат. наук, Москва: ОИФЗ РАН, 2003.
11. Ануфриев А.П., Решетняка М.Ю., Соколов Д.Д. Оценка динамо-числа в модели турбулентного
α-эффекта для жидкого ядра Земли // Геомагнетизм и аэрономия. 1997. Т. 37. С. 141-146.
12. 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.
Поступила в редакцию / Original article submitted: 17.11.2015
54
Документ
Категория
Без категории
Просмотров
3
Размер файла
501 Кб
Теги
конвекция, управляемое, ячейковой, геодинамо, модель, инверсия
1/--страниц
Пожаловаться на содержимое документа