Файл: Горбацкий, В. Г. Новоподобные и новые звезды.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 24.10.2024

Просмотров: 66

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

144

ГЛ. VI. ВСПЫШКИ ЗВЕЗД

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

&ть

0,0

 

з января 7082

 

0.5

 

 

 

0,6

 

 

 

0,0

.

6

января

0,7

 

 

 

1.3Иянваря

1,0

 

 

 

 

 

 

 

 

1,7

 

 

 

 

 

 

!2января

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,0

0,2

Ofi

0,0

0,8

0,0

0.2

0,0

0,6 Р

Рис. 39. Кривая блеска системы U Gem во время вспышки [13].

[95 — 97]. Тем самым указанное «преимущество» гипоте­ зы о вспышках горячей звезды в системах U Gem отпа­ дает. Вместе с тем, становятся понятными многие из спект­ ральных данных, трудно объяснимые в рамках гипотезы, связывающей вспышки с горячей звездой, и выясняется причина радикальных отличий вспышек U Gem от вспы­ шек новых. Поэтому первоначальный вывод [13] о локали­ зации вспышек типа U Gem в холодном спутнике системы представляется более предпочтительным. В следующих разделах исследуется тепловая неустойчивость, возни­ кающая в такой звезде.

§ 2. Нестационарная конвекция в звезде—компоненте тесной двойной системы

Во внешних областях звезды карлика позднего спект­ рального класса (G — М), которая, в соответствии со ска­ занным выше, считается источником вспышечной актив­ ности типа U Gem, перенос энергии; осуществляется путем


2.

НЕСТАЦИОНАРНАЯ КОНВЕКЦИЯ

145

конвекции.

Условия

конвективного переноса

энергии

в звезде — компоненте

тесной двойной системы — отли­

чаются от условий в одиночной звезде. Одиночная звезда является самобалансирующейся системой, и конвектив­ ный перенос в ней приспосабливается к производительно­ сти внутренних источников энергии. В том случае, когда звезда входит в состав тесной двойной системы, конвек­ тивный поток энергии, помимо внутренних условий, зави­ сит еще и от внешних по отношению к звезде обстоятельств, и приспособления к имеющимся источникам энергии мо­ жет и не быть. В частности, оказывается возможной веко­ вая неустойчивость конвективного потока, приводящая к накоплению энергии во внутренних слоях конвективной зоны. Вопрос о том, в какой мере этим эффектом объяс­ няются вспышки типа U Gem, рассматривается в следую­ щем параграфе. Здесь же излагаются результаты теорий нестационарной (зависящей от времени) конвекции в пе­ риодически меняющемся гравитационном поле. Эта теория была развита в работах Л. Н. Иванова [96 — 98].

При асинхронности вращения вокруг оси и орбитально­ го обращения звезды, входящей в состав тесной двойной системы, ускорение силы тяжести на ее поверхности g периодически изменяется. Для простоты в дальнейшем внешние области звезды моделируются плоским слоем газа, в котором величина g не меняется с глубиной. Так считать допустимо, если протяженность конвективной зоны мала по сравнению с радиусом звезды. Изменение g со временем задается выражением

 

g (t) = g0 (1. + е cos 2 соt),

(1 .6)

где со =

со0бр — швр. Для тесных двойных

систем звезд

карликов

е ^ 0,4.

 

Конвективный поток энергии Fc существенно зависит от величины ускорения силы тяжести и, значит, в данном случае он должен меняться со временем. Задача заклю­ чается в определении характера изменений при условии (1.6). Как предполагается в дальнейшем, поток энергии, подводимой к нижней границе конвективного слоя, по­ стоянен. Это означает, что к конвективной зоне примы­ кает зона лучистого равновесия, в которой можно прене­ бречь влиянием переменности g на величину потока лучи­ стой энергии.

' 3/зб В. Г. Горбацкий


146

ГЛ . V I. ВСПЫ Ш КИ ЗВ Е ЗД

Конвекция во внешних слоях звезд имеет турбулент­ ный характер. Поскольку строгой теории турбулентной конвекции до сих пор нет, в работах [96 — 98] использо­ ван метод расчета конвективного потока, основанный на понятии длины перемешивания при учете модификации метода, предложенной Унно [99]. Скорости конвективных движений считаются дозвуковыми, и поэтому применяется приближение Буссинеска, при котором сжимаемость сре­ ды учитывается только при вычислении силы плавучести, действующей на конвективный элемент. Общее выражение зависимости потока энергии Fc от постоянного ускорения силы тяжести предложено в работе [1 0 0 ].

Вычисление зависимости потока F с от времени, выпол­ ненное при указанных предположениях в [96], дало сле­ дующий результат:

Fc(0 = Fс (go)

где

я гГ

1 + а г ____sin (at + а) ,

( 2.6)

У 1 + 7

 

и а — сдвиг по фазе,

(3.6)

а величина г)-1 определяет характерное время релаксации конвективных элементов и, следовательно, также и кон­ вективного потока. Выражение для ц имеет вид:

■4 =

- 1 + 1

2&г„^ V/Jl

(4.6)

TY* j J ’

 

где R c — размер конвективных элементов, %— коэффи­ циент теплопроводности и [J—сверхадиабатический темпе­ ратурный градиент

'

d T

(5.6)

Р = -

dz

При р не зависящем от времени для сжимаемой среды ве­ личина а, входящая в (2.6), равна 3/2.

В том случае, когда время релаксации конвективного потока очень велико по сравнению с периодом колебаний, из (2 .6) получается

F c ( t ) ^ F c (g0),

т. е. конвективный поток не успевает реагировать на изме­ нения ускорения силы тяжести. Если же q<^. 1, то поток


§ 2. НЕСТАЦИОНАРНАЯ КОНВЕКЦИЯ

147

сразу следует в своем изменении за g (2) и поэтому зависит от времени по синусоидальному закону с амплитудой порядка е.

При выводе соотношения (2.6) изменения конвектив­ ного потока считались адиабатическими, т. е. не прини­ малась во внимание их взаимосвязь с состоянием среды, в которой переносится энергия. Если же учитывать неадиабадичность элементов, то структура конвективной зоны очень существенно зависит от изменений g (2) и опре­ деление величины конвективного потока оказывается го­ раздо более сложным. Задача о зависимости конвективно­ го потока от времени в этом случае рассмотрена в работах [97, 98] при тех же предположениях, т. е. для плоского конвективного слоя с постоянным потоком энергии Н, приходящей к его внутренней границе, и величине g (2) в слое, меняющейся согласно (1 .6).

Распределения температуры, плотности и давления в плоском конвективном слое в зависимости от глубины и времени можно представить в следующем виде:

Т (т\ t) =

Т0(то) ^1 +

е sin со2 + - ~

j ,

Р {m; 2) =

ро (то) (l +

- у - 8 sin at + -jy) ,

(6.6)

p (to.; 2) =

p0(m) (1 +

e sin at).

 

где то — масса, содержащаяся в столбе единичного сече­ ния от внутренней поверхности до данной точки, прини­

маемая за лагранжеву координату,

и Т0 (то),

р0 (то) и

Ро im) — значения

соответствующих

величин

в стацио­

нарном состоянии,

у — показатель

адиабаты.

Через

и pj обозначены изменения температуры и плотности вслед­ ствие неадиабатичности. Членами, содержащими множи­ тель е sin со2, учитывается действие динамических при­ ливов.

В тех случаях, когда период изменения g намного пре­ восходит время механической релаксации слоя, вместо уравнения движения можно использовать уравнение гид­ ростатического равновесия:

dp (/гс; /)

- g i t ) .

(7.6)

dm

s*


-148

ГЛ . V I. ВСПЫ Ш КИ З В Е ЗД

В карликовой компоненте тесной двойной системы со-1 ^ 5h -f- 10й, а время механической релаксации кон­ вективной зоны порядка времени распространения звука сквозь нее и составляет всего несколько минут. Поэтому использование (7.6) при исследовании конвекции в этих условиях вполне допустимо.

В предположении об адиабатичности приливных воз­ действий из уравнения теплопроводности получается уравнение для определения величины Тх\

д_(Т_i_\

___ 1__ д£_

( 8. 6)

dt \ То )

СрТо

dm

 

где F — общий поток энергии,

как

конвективный F с,

так и переносимой излучением Fr\

 

F =

Fr 4 F c.

(9.6)

При помощи уравнения состояния идеального газа находится также уравнение, описывающее изменение со временем рх:

д_

= (l + 2 ^

е sin at

1

dF

( 10. 6)

dt

СрТо

dm

В (8 .6) и (9.6) dldt означает локальную (лагранжеву) про­ изводную.

Для потока лучистой энергии в [98] принято следующее выражение:

т

T s

d T

 

( 11.6)

Fr

р

dm

где L — некоторый коэффициент.

Величина F с опреде­

ляется по (2 .6).

 

 

rp1

 

 

 

 

и —

При аналитическом решении задачи величины

удобно представить в виде

 

 

1 о

Ро

 

 

 

 

- j t = бх + 0;

 

=

°i + о,

(12.6)

где 0Хи ах — периодические составляющие, а 0 и о — мед­ ленно меняющиеся (вековые) составляющие, присутствие которых непосредственно устанавливается при численном решении системы (8 .6) — (11.6). Вековая составляющая возникает при этом и в выражении для конвективного по-