ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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). Вековая составляющая возникает при этом и в выражении для конвективного по-