24 / 32
M
≈
3
M
J
, окруженной околозвездным диском радиусом
R
d
≈
40
а.е.
Диск пронизан полоидальным магнитным полем, напряженность кото-
рого составляет около 0,06 Э. На систему аккрецирует сверхзвуковой
поток вещества с темпом аккреции около
5
∙
10
−
5
M
J
/
год. Перпен-
дикулярно экваториальной плоскости диска формируется коллимиро-
ванное истечение вещества мощностью
0
,
65
∙
10
−
5
M
J
/
год, скорость
потока составляет 500 км/с, угол раствора джета — около
10
◦
.
Модель ускорения плазмы в канале над горячим аккреционным
диском.
В рамках РМГД-задачи схема остается неизменной, в рассмо-
трение включается давление излучения тонкого диска.
Полная РМГД-система в квазистационарном моноэнергетическом
приближении имеет вид (1), (3)–(7).
Рассмотрим томпсоновское рассеяние излучения на свободных
электронах. Индикатриса такого рассеяния равна рэлеевской
Γ (
t, x, ω, ω
0
) =
3
4
1 + (
ω
∙
ω
0
)
2
, а сечение рассеяния равно
σ
T
=
= 6
,
65210
−
29
см
2
(коэффициент рассеяния
β
=
nσ
T
, где
n
— концен-
трация вещества).
Модельные условия для газовой и магнитной части РМГД-системы
в целом повторяют условия для МГД-модели. Предполагается, что ме-
ханизм образования ветра над тонким диском функционирует посто-
янно, вне зависимости от угла наклона магнитных силовых линий.
Подобное условие обеспечивается выталкивающим действием излу-
чения внутри тонкого диска.
Кроме того, предполагается, что окрестности центрального гра-
витирующего объекта (часть тонкого диска) излучают с интенсивно-
стью излучения абсолютно черного тела с характерной температурой
7
∙
10
4
K. При этом излучение сфокусировано внутрь канала: излучаю-
щей является граница
z
= 0
,
0
≤
r
≤
0
,
2
, излучение распространяется
вдоль направлений, для которых
cos
θ >
0
,
9
, где
θ
— полярный угол
луча.
Примем, что ослабление интенсивности излучения в процессах
рассеяния существенно выше интенсивности в процессах поглощения
фотонов веществом. Другими словами, вещество достаточно разреже-
но и достаточно слабо поглощает, поэтому можно принять коэффици-
ент поглощения
α
(
t, x, ω
) = 0
.
Численный метод и его программная реализация.
Для решения
РМГД-системы применен уже упомянутый метод дробных шагов, до-
полненный этапом решения УПИ. С позиции вычислений наиболее
требовательной к ресурсам является процедура решения УПИ. Для
интегрирования УПИ выбран метод дискретных направлений.
Для решения уравнения (7) в соответствии с методом дискретных
направлений (МДН) необходимо провести интегрирование УПИ вдоль
86
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2015. № 2