Курс теоретической астрофизики
Шрифт:
(2.25)
Вынесем за знак интеграла среднее значение cos^2 на сфере, равное ^1/ т.е. приближённо положим
0
I(,)
cos^2
sin
d
=
1
3
0
I(,)
sin
d
.
(2.26)
Тогда вместо (2.25) при учёте второго из уравнений (2.9) находим
4
3
dS
d
=
H
.
(2.27)
Так
S
=
3
4
H
+
C
,
(2.28)
где C — произвольная постоянная.
Для нахождения C напишем выражение для величин S и H при =0. Принимая во внимание граничное условие (2.10), находим
S(0)
=
1
2
0
I(0,)
sin
d
,
(2.29)
а также приближённо
H
=
0
I(0,)
sin
d
.
(2.30)
Поэтому имеем
S(0)
=
H
2
.
(2.31)
При условии (2.31) для постоянной C получаем
C
=
H
2
.
(2.32)
Подстановка (2.32) в (2.28) даёт
S
=
F
3
4
+
1
2
,
(2.33)
где, как и раньше, использовано обозначение (2.23).
Мы видим, что выражение (2.33) для функции S не сильно отличается от выражения (2.24), полученного предыдущим методом.
3. Применение квадратурных формул.
Изложенные выше приближённые методы нашли довольно широкое применение в астрофизике. Однако точность результатов, получаемых этими методами, сравнительно невелика. Поэтому получил распространение другой приближённый метод, основанный на замене интегрального члена уравнения лучистого равновесия суммой Гаусса для численных квадратур. Уравнение переноса излучения пишется при этом для тех значений cos, которые являются точками деления интервала в квадратурной формуле. Это позволяет свести задачу к системе линейных дифференциальных уравнений с постоянными коэффициентами.
Преимущество этого метода состоит в том, что можно повышать точность результатов, увеличивая число членов квадратурной формулы. Однако и при небольшом числе членов этой формулы получаются удовлетворительные результаты благодаря высокой точности замены интеграла суммой Гаусса.
Указанный метод был подробно разработан Чандрасекаром [4]. Мы сейчас применим этот метод к решению системы уравнений (2.9).
Предварительно перепишем эту систему в виде одного уравнения:
dI(,)
d
=
I(,)
–
1
2
+1
– 1
I(,')
d'
,
(2.34)
где обозначено =cos.
Представим интегральный член уравнения (2.34) в виде суммы согласно квадратурной формуле Гаусса:
+1
– 1
I(,)
d
=
n
j=-n
a
j
I(,
j
)
.
(2.35)
Здесь – n,…,– 1,1,…n суть корни полинома Лежандра P2n и aj — некоторые весовые множители (a– j=aj). Представление (2.35) тем точнее, чем больше n
В n-м приближении уравнение (2.34) заменяется системой линейных дифференциальных уравнений порядка 2n:
i
dIi
d
=
I
i
–
1
2
j
a
j
I
j
(i
=
±1,
±2,
…,
±n
),
(2.36)
где для краткости I(,) обозначено через Ii.
Произвольные постоянные, входящие в общее решение этой системы, определяются из следующих условий: 1) отсутствует излучение, падающее на фотосферу извне, т.е. I– i=0 при =0 (i=1,2, ,n); 2) не может быть членов, экспоненциально возрастающих с , 3) задан поток излучения H=F.
После нахождения величин Ii из уравнений (2.36) основная искомая функция S определяется по формуле
S
=
1
2
a
j
I
j
.
(2.37)
Найдём в виде примера функцию S в первом приближении. В данном случае 1=-– 1=1/3, a1=a– 1=1. Поэтому вместо (2.36) получаем
1
3
dI1
d
=
I
1
–
1
2
(
I
1
+
I
– 1
),
–
1
3
dI– 1
d
=
I
– 1
–
1
2
(
I
1
+
I
– 1
).
(2.38)
Система уравнений (2.38) должна быть решена при условиях, что I– 1=0 при =0 и
2
3
(
I
1
+
I
– 1
)=
F
.
(2.39)