Курс теоретической астрофизики

Таким образом, полный поток излучения (т.е. поток излучения, проинтегрированный по всему спектру) в сферически-симметричной фотосфере обратно пропорционален квадрату расстояния от центра звезды. Соотношение (1.23), как и уравнение (1.17), является следствием отсутствия источников и стоков энергии в фотосфере.

Как уже говорилось, почти все звёзды обладают фотосферами, толщина которых очень мала по сравнению с радиусом звезды. Для этих звёзд уравнения (1.20) и (1.23) могут быть сильно упрощены. Этого нельзя сделать лишь для звёзд особых типов (например, для звёзд типа Вольфа — Райе).

Рис. 2

Если толщина фотосферы гораздо меньше радиуса звезды, то фотосферные слои могут считаться не сферическими, а плоскопараллельными (рис. 2). В этом случае угол не меняется вдоль луча и вместо уравнения (1.20) получаем

cos

dI

dr

=-

I

+

.

(1.24)

Так как расстояние r от центра звезды меняется в фотосфере в очень небольших пределах, то вместо уравнения (1.23) имеем

0

H

d

=

const.

(1.25)

Таким образом, при рассмотрении поля излучения в фотосферах «обычных» звёзд следует пользоваться уравнениями (1.24) и (1.17) или уравнениями (1.24) и (1.25).

§ 2. Теория фотосфер при коэффициенте поглощения, не зависящем от частоты

1. Основные уравнения.

Первоначально в теории фотосфер делалось предположение о независимости коэффициента поглощения от частоты, ведущее к существенному упрощению теории. В дальнейшем, однако, было установлено, что это предположение является весьма грубым. Тем не менее теория фотосфер при коэффициенте поглощения, не зависящем от частоты, продолжает сохранять своё значение, так как она может рассматриваться как первое приближение к более строгой теории.

Считая, что коэффициент поглощения не зависит от частоты (т.е. =), вместо уравнения переноса излучения (1.24) и уравнения лучистого равновесия (1.17) получаем

cos

dI

dr

=-

I

+

,

(2.1)

4

0

d

=

d

0

I

d

.

(2.2)

Введём обозначения

0

I

d

=

I

,

0

d

=

.

(2.3)

Величину I можно назвать полной интенсивностью излучения, а величину — полным коэффициентом излучения.

Проинтегрировав уравнение (2.1) по всем частотам, находим

cos

dI

dr

=-

I

+

,

(2.4)

а уравнение (2.2) переписывается в виде

4

=

I

d

.

(2.5)

При исследовании переноса излучения в любой среде целесообразно переходить от геометрических расстояний к оптическим расстояниям. В данном случае удобно ввести оптическую глубину , определяемую формулой

=

r

dr

(2.6)

Положим также

=

S

.

(2.7)

Тогда уравнения (2.4) и (2.5) принимают вид

cos

dI

d

=

I-S

,

S

=

I

d

4

.

(2.8)

Таким образом, мы получили два уравнения для определения двух неизвестных функций I и S.

В системе уравнений (2.8) величина I является функцией от и , а величина S — функцией от . Учитывая, что d=sin d d, и производя интегрирование по в пределах от 0 до 2, вместо (2.8) получаем

cos

dI(,)

d

=

I(,)

-

S

,

S

=

1/2

0

I(,)

sin

d

.

(2.9)

К системе уравнений (2.9) необходимо добавить ещё граничное условие. Оно выражает тот факт, что нет излучения, падающего на звезду извне, т.е.

I(0,)

=

0

при

>

2

.

(2.10)

Кроме того, для получения вполне определённого решения системы уравнений (2.9) при граничном условии (2.10) следует задать ещё полный поток излучения в фотосфере, равный

H

=

L

4R^2

,

(2.11)

где L — светимость звезды (т.е. полное количество энергии, излучаемое звездой за 1 с) и R — радиус звезды.

Системы уравнений типа (2.9) весьма часто встречаются в астрофизике. С такими же уравнениями приходится иметь дело и в геофизике (при изучении рассеяния света в земной атмосфере и в водных бассейнах). К аналогичным уравнениям приводят и некоторые проблемы физики (например, проблема диффузии нейтронов). Поэтому системы уравнений типа (2.9) были предметом многочисленных исследований и для их решения предложен ряд методов (см. [4] и [5]).

Ниже излагаются некоторые из этих методов, представляющих наибольший интерес для астрофизики.

2. Приближённое решение уравнений.

Для решения системы уравнений (2.9) были предложены приближённые методы, основанные на усреднении интенсивности излучения по направлениям. Первый из этих методов принадлежит Шварцшильду и Шустеру, второй — Эддингтону. Мы сейчас решим систему уравнений (2.9) при помощи каждого из указанных методов.

Метод Шварцшильда — Шустера. Обозначим через I среднюю интенсивность излучения, идущего снизу вверх, и через I — среднюю интенсивность излучения, идущего сверху вниз. Эти величины равны

I

=

/2

0

I(,)

sin

d

,

I

=

/2

I(,)

sin

d

.

(2.12)

Умножая первое из уравнений (2.9) на sin d и интегрируя в пределах от 0 до /2, получаем

d

d

/2

0

I(,)

cos

sin

d

=

I

-

S

.

(2.13)

Интеграл в левой части этого уравнения приближённо представим в виде

/2

0

I(,)

cos

sin

d

=

1/2

I

,

(2.14)

т.е. вынесем за знак интеграла среднее значение cos в верхней полусфере, равное 1/2 . Тогда вместо (2.13) будем иметь

1

2

dI

d

=

I

-

S

.

(2.15)

Умножая первое из уравнений (2.9) на sin d и интегрируя в пределах от /2 до , аналогично находим

-

1

2

dI

d

=

I

-

S

.

(2.16)

Второе из уравнений (2.9) при помощи величин I и I переписывается так:

S

=

1/2 [

I

+

I

]

(2.17)