АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 2015, том 70, № 2, с. 224-230
УДК 524.352:524.1- 782
ВЗРЫВНЫЕ ВОЛНЫ С КОСМИЧЕСКИМИ ЛУЧАМИ
2015 Б. Арбутина*
Белградский университет, Белград, 11000 Сербия
Поступила в редакцию 8 декабря 2014 года; принята в печать 1 апреля 2015 года
Взрывные волны присутствуют во многих астрофизических явлениях, таких как сверхновые. В данной статье мы обсуждаем взрывные волны с космическими лучами, т.е. с компонентой со степенной функцией распределения плотности N(р) гс р-Г, которая может быть особенно важной при описании эволюции остатков сверхновых. Мы подтверждаем некоторые предыдущие результаты о том, что существенное количество энергии космических лучей выделяется в направлении центра остатка.
Ключевые слова: ударные волны — сверхновые
1. ВВЕДЕНИЕ
Взрывные волны появляются при многих астрофизических явлениях, таких как сверхновые. В данной статье мы обсуждаем взрывные волны с космическими лучами (КЛ), т.е. с компонентой со степенной функцией распределения плотности N(р) гс р-Г, где р — импульс частицы, а Г — так называемый индекс энергии. Такие объяснения ударных волн могут быть особенно важны при описании эволюции остатков вспышек сверхновых, которые считаются основным источником галактических КЛ вплоть до энергий 1015 эВ.
Ударные волны с давлением космических лучей уже обсуждались в работах Шевалье [1,2], Топтыгина [3], Петухова и др. [4] и, недавно, Белла [5]. Как и в работах Шевалье, мы будем полагать, что КЛ ускоряются на фронте ударной волны и адиабатически остывают в области за фронтом. Тем не менее, мы расширили анализ путем использования немного модифицированных условий скачка и проанализировали решения с разными индексами энергии.
2. АНАЛИЗ И РЕЗУЛЬТАТЫ
2.1. Анализ подобия
Хорошо известно, что уравнения Эйлера позволяют автомодельные решения с введенной переменной £ = crt-x, откуда получаем распространение радиуса как R гс tx. Чтобы оценить эволюцию ударной волны со временем, можно воспользоваться теорией подобия. В предположении, что остаток сверхновой «забыл» точные
E-mail: [email protected]
начальные условия взрыва сверхновой, только параметры Eo и po имеют значение. По размерности, Eo = [эрг] = [гсм2 с-2], а po = [гсм-3]. Исключая [г], получаем Eo х [см-2с2] = po х [см3], т.е. [см5] = (Eo/po) х [с2] ,поэтому R гс (Eo/po)1/5t2/5 и Л = 2/5.
Чтобы найти полное решение, перепишем переменную подобия в виде:
£ = r/R(t), R(t) = B(Eo/po)1512/5, (1)
где для нахождения B нужно решить уравнения Эйлера. Если мы просто имеем дело с обычным газом с Yg = 5/3, получим известное решение Седова с R = 1.15(Eo/po)1/512/5 [6]. Если подобие сохраняется, ввод новой компоненты, связанной с космическими лучами, не изменит зависимость R = R(t), а повлияет только на значение B.
2.2. Уравнение состояния космических лучей
Предположим, что КЛ имеют степенное распределение плотности N(р) гс р-Г, т.е., что функция распределения в импульсном пространстве имеет вид f = kp-r-2, где k — некая константа. По определению, давление и внутренняя энергия на единицу объема есть
Рж
р = f 4жр>у;Пр, (2)
Р0
рж
u = j 4np2(E — mc2)f dp, (3)
Р0
224
ВЗРЫВНЫЕ ВОЛНЫ С КОСМИЧЕСКИМИ ЛУЧАМИ
225
где энергия задается релятивистским выражением E2 = p2c2 + m2c4, m — масса частицы, c — скорость света, а р0 и pж — некоторые нижний и верхний пределы импульса соответственно.
Тогда для давления имеем
P
4пк
3
2—Г 2 -р С
Р0
\/р2с2 + т2с4
dp. (4)
p
1
Интегрируя по частям, находим плотность внутренней энергии
и = [ 4тгкр~г (\/р2с2 + m2c4 - тс2) dp
Р0
3
Г - 1
P +
Атгк Г - 1
рО г
\JpIc2 + т2с4 — mcY
—Рсо Г (л/рХ2 + w-2c4 — тс2)
(5)
Если po ^ mc и p^ ^ то, для 2 < Г < 3 имеем
3
(6)
Т.к. и = —-—Р,
Y — 1
Г + 2
7 = 7cr = —g—• (7)
Видно, что КЛ с Г 2 ведут себя как ультрарелятивистский газ (7 = 4/3), тогда как КЛ с Г 3 — как обычный газ с 7 = 5/3.
2.3. Условия скачка
Если разделить ударную волну на три области:
(0) дальнюю прифронтовую,
(1) ближнюю прифронтовую (предшественницу),
(2) зафронтовую,
то сохранение массы, импульса и энергии дают соотношения Ренкина—Гюгонио для сильных ударных волн
Мы не предполагали сохранение энергии в субволне (между областями (2) и (1)); положим Pcr2 = Pcri [7], т.е., будем считать, что давление КЛ непрерывно — высокоэнергетические частицы не видят ударной волны, т.е. субволны.
Закон сохранения энергии также можно написать в виде
1
Х2 +
7 Р2 7 - 1 р2
До,
1
(11)
где мы определили эффективный 7-фактор таким образом, что
7
7 - 1
7д
79 - 1
в+
7с R
Ycr - 1
(1 - в),
(12)
где в = Pg2/P2, P2 = Pg2 + Pcr2. Для полного скачка давления тогда имеем
Xtot
Р2
po
7 + 1 7 - 1
(13)
и
Рд2 = ^уРоХ (14)
PcR 2 = ^ ^РоХ (15)
Y + 1
V2 = —y-rVs, (16)
Y + 1
где vs = R = —v0 — скорость ударной волны. Мы также знаем из теории ускорения рассеянных волн [8—10], что (см. уравнение (7))
X _ Р2 _ Г + 2 _ 7CR Pi Г — 1 7cr - 1 ’
(17)
а из закона сохранения импульса мы можем получить соотношение
X — 1 = e(Xtot — 1). (18)
Совмещая (6), (7), (11) и (12), мы окончательно имеем (см. также рис. 1)
в
1
Д _ р7й ~7cr / Yg — 1
О7д — 7cr 7д - 1
(19)
Y = 1 + 2e(YcR — 1). (20)
p2v2 = Plv1 = p0v0, p2v2 + Pg2 + PCR2 = Plv2 + PCR1
1 . Yg Pg2 . YCR PCR2
Kv2 4 7-----1 7-----
2 Yg — 1 p2 Ycr — 1 p2
(8)
po v0, (9)
1 2 2V°- (10)
2.4. Автомодельный вид уравнений Эйлера
Используя уравнения (13)—(16), мы можем определить безразмерные плотность F, давление газа G, давление КЛ H и скорость U следующим образом:
15 АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 2 2015
226
АРБУТИНА
Рис. 1. Зависимость в и y от ycr.
Р = ^PoFiO
Pg
PCR =
Y - 1
_w_
Y + 1 2(1 - в) Y2
Y + 1 2
poR 2G(£) po R 2h (0
Y+1
R U (£).
На фронте ударной волны (£ = 1) граничные условия будут F(1) = G(1) = U(1) = 1. С этими последними выражениями и частными производными
<L = i>(d_-iJL
т \дп Пд£
д 1 д dr Rd£’
мы можем перейти к сферически-симметричным уравнениям Эйлера
dp 1 д , 2 ч
gf + Ya;(r'"’)=0-
Dv
dv
dv
PDt Pdt+pVdr dr D(Pg p~la ) d
dPg dPcR
dr
Dt
d
+vd~r (Pgp 79) = °>
D{Pckp-^)_ = 8 {PcRp-ICR)
Dt
d
+v— (PcrP“7cr) = 0,
(30)
D
(21)
(22)
(23)
(24)
где — — производная Лагранжа или субстанциональная производная.
Уравнение непрерывности (27) дает нам
!-^)4+г/'иг/=°- <з1>
а с помощью (27), уравнения энтропии (29) и (30) можно проинтегрировать и получить
^ (Г _ Y + i) = 1_7_ (32|
HFl-lcn(lL_l±l\ =
U 2 ) 2С3 •
(33)
(25)
(26)
£ 2 ) 2£3
Во всех выражениях мы будем использовать сле дующую систему обозначений:
/ = |[, f = {F,G,H,U}.
Наконец, напишем автомодельный вид уравнения импульса (28):
(27) £
(28)
7 + U ) 2{
[@G' + (1 - в)Н'] = 0.
(34)
(29)
Если использовать уравнения (32) и (33), чтобы найти G и H', и подставить их в (34), то с небольшими преобразованиями уравнения (34) и(31)становятся
v
АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 2 2015
ВЗРЫВНЫЕ ВОЛНЫ С КОСМИЧЕСКИМИ ЛУЧАМИ
227
<Ш_ _ 1
(7 - 1)2
^~* 2 ( -(7 + 1)С - 2Ъи ) + (1 - Р)ГГ°*
2 ' ^(7 + 1)£ -27CRи
(7 - 1)2 * [в7дFYg 2 + (1 - e)7cRFycr-2} + 4£2
7 + 1
и-1—с 2 s
т 3
3(7 + 1)c3u
+-
7+1
и-1—с
2 s
2
(7 - 1)2 [в7дFYg 2 + (1 - e)7cRFycr-2} + 4£2
7+1
и-Ц-с
3
(35)
d F
d£
F
-(7 + 1)(7 - I)2 [l3F^~2 + (1 - 13)Г'1™~2]
u-l±±t 2 s
(7 - 1)2 [в7дFYg-2 + (1 - e)7cRFYcr-2} + 4£2
CU [8U - (7 + 1)]
+-
7+1
и-1—с
2 s
2
(7 - 1)2 [в7дFYg 2 + (1 - e)7cRFYCR-2] + 4£2
7 +1
U- 1—с 2 s
(36)
3
3
Эти два последних уравнения нужно решать численно; для этого, мы воспользовались интегратором Рунге—Кутты шестого порядка [11]. Результаты показаны на рис. 2. По мере приближения к центру C 0, если в = 1, то ско-
7 +1
рость становится приблизительно U ~-----£, а
27cr
F х C3/(YCR-1), G х С3(ъ-YCR)/(ycr-1), H (х const. Если в = 1,получаем результаты решения Седова:
U зз F ос и G ос const.
27д
2.5. Интеграл энергии
Рассмотрим теперь движущуюся поверхность с C = const. Энергия в пределах этой поверхности
равна
E (C)
1 2 1 1
~PV Н тРд Н т
2 7д - 1 7cr - 1
PcR^ 4лг2 dr
8л
7 + 1
poR 2R
2 ^ I I —-—FU2 + —rG + ——-Я] £2 d£
7 - 1 7д - 1 7cr - 1
32?r B5Eq
25(7 + 1) J V7-1 0
1 FU2 + —^—G + 1 13 h] i2 d£,
7д - 1 7cr - 1 J
0
0
(37)
АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №2 2015
15
АРБУТИНА
228
Рис. 2. Полные решения для разных индексов энергии Г. Обозначения одинаковы на всех рисунках.
где мы воспользовались тем фактом, что
R2R3 = 5
25ро
Если для движущейся поверхности мыпримем
£ = 1 (фронт ударной волны), уравнение (37) даст нам полную энергию взрыва Ео в виде суммы кинетической энергии, внутренней энергии газа и энергии КЛ:
АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 2 2015
ВЗРЫВНЫЕ ВОЛНЫ С КОСМИЧЕСКИМИ ЛУЧАМИ
229
Параметры решений для разных индексов энергии Г. Eo — энергия взрыва, To — полная кинетическая энергия, U0,g — полная тепловая энергия, U0,CR — полная энергия КЛ. В последних двух колонках представлены соотношения давлений в центре и на фронте ударной волны для газа и КЛ — Pg0/Pg2 и PCR0/PCR2. Другие параметры см. в тексте.
Г Ys Ycr /3 Y X Mot В To/Eo Uo,g/Ea Yo,cr/E0 Pgo/Pg2 -Pcro/ PcR2
2 5/3 4/3 1 1.67 4 4 1.152 0.283 0.717 0 0.306 0
2.1 5/3 1.37 0.76 1.56 3.73 4.59 1.090 0.245 0.356 0.399 0 1.401
2.2 5/3 7/5 0.69 1.55 3.5 4.62 1.092 0.247 0.324 0.429 0 1.088
2.3 5/3 1.43 0.65 1.56 3.31 4.57 1.099 0.252 0.315 0.433 0 0.941
2.4 5/3 1.47 0.61 1.57 3.14 4.5 1.107 0.257 0.315 0.428 0 0.849
2.5 5/3 3/2 0.59 1.59 3 4.41 1.115 0.262 0.319 0.418 0 0.784
2.6 5/3 1.53 0.56 1.60 2.88 4.33 1.123 0.267 0.325 0.407 0 0.733
2.7 5/3 1.57 0.54 1.62 2.76 4.24 1.131 0.272 0.333 0.395 0 0.687
2.8 5/3 8/5 0.53 1.63 2.67 4.16 1.138 0.276 0.341 0.383 0 0.623
2.9 5/3 1.63 0.51 1.65 2.58 4.08 1.145 0.279 0.350 0.371 0 0.501
3 5/3 5/3 0.5 1.67 2.5 4 1.152 0.283 0.359 0.359 0.306 0.306
E0 — T0 + Uo,g + Uq, cr r
0
\pv2 4---l~:Pg H------—tPcr ) 4тгr2 dr
2 Yg - 1 Ycr - 1
32л
25(y + 1)
что позволит нам найти B:
B5Eq
1 -FU2 + -^—rG + ——Lh] C2d{,
B-5 —
32л
25(y + 1) J \y - 1 0
Y - 1 Yg - 1 Ycr - 1'
1 FU2 + —^G + ——(2d^,
Yg - 1 Ycr - 1
i
0
1
(38)
(39)
а также доли Tq/Eq, Uo,g/Eq и U0,cr/Eq (см. таблицу).
3. ЗАКЛЮЧЕНИЕ
В данной работе мы обсудили ударные волны с космическими лучами, т.е. с компонентой со степенной функцией распределения плотности N(р) « р-Г, которая приводит к возникновению дополнительного давления Pcr, связанного с КЛ.
Мы проанализировали решения с разными индексами энергии Г и вычислили параметры этих решений.
Мы подтверждаем некоторые более ранние сведения о том, что значительная доля энергии космических лучей выделяется в сторону центра остатка [2]. В большинстве наших моделей давление КЛ в центре сравнимо с давлением на фронте ударной волны. Если плотности энергии КЛ и магнитного поля распределены
1 Г + 1
равномерно [12, 13], ив = —В2 зз-----uCR, т.к.
8л 4
АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №2 2015
230
АРБУТИНА
синхротронная объемная излучательная способность ev х ugrB 1+аv а, где v — частота, а а = (Г — 1)/2 — спектральный индекс, то в дополнение к излучению от фронта ударной волны возможно также наличие сравнимого излучения из центра. Поэтому такое выделение энергии КЛ может служить объяснением для радиоизлучения некоторых остатков сверхновых с заполненным центром.
БЛАГОДАРНОСТИ
Во время работы над данной статьей автору была оказана финансовая поддержка Министерства образования, науки и технологического развития Республики Сербия (проекты 176004 «Звездная физика» и 176005 «Эмиссионные туманности: структура и эволюция»). Автор также хотел бы воспользоваться возможностью поблагодарить As kin Ankay и Deniz Go ker за гостеприимство и возможность представить данную тему на семинаре физического факультета Босфорского университета Стамбуле в феврале 2015 г..
СПИСОК ЛИТЕРАТУРЫ
1. R. A. Chevalier, in Proc. 18th Int. Cosmic Ray Conf. (Tata Institute of Fundamental Research, 1983), Vol. 2, p. 314.
2. R. A. Chevalier, Astrophys. J. 272, 765 (1983).
3. I. N. Toptygin, Astronomy Letters 26, 356 (2000).
4. Yu. V. Petukhov, A. V. Razin, and V. A. Razin, Astronomy Letters 32, 747 (2006).
5. A. R. Bell, Monthly Notices Royal Astron. Soc. 447, 2224 (2015).
6. L. I. Sedov, Similarity and Dimensional Methods in Mechanics (Academic Press, New York, 1959).
7. L. O. C. Drury, Rep. Prog. Phys. 46, 973 (1983).
8. A. R. Bell, Monthly Notices Royal Astron. Soc. 182, 147(1978).
9. A. R. Bell, Monthly Notices Royal Astron. Soc. 182, 443 (1978).
10. R. D. Blandford and J. P. Ostriker, Astrophys. J. 221, L29 (1978).
11. H. A. Luther, Math. Comp. 22, 434 (1968).
12. B. Arbutina, D. Urosevic, M. M. Andjelic, et al., Astrophys. J. 746,79(2012).
13. B. Arbutina, D. Uros;evic, M. M. Vucetic, et al., Astrophys. J. 777, 31 (2013).
Перевод Е. Чмырёвой
Blast Waves with Cosmic Rays B. Arbutina
Blast waves appear in many astrophysical phenomena, such as supernovae. In this paper we discuss blast waves with cosmic rays, i.e., with a component with a power-law number density distribution function N(p) х p-r that may be particulary important in describing the evolution of supernova remnants. We confirm some previous findings that a significant amount of cosmic ray energy is deposited towards the center of a remnant.
Keywords: shock waves—supernovae: general
АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70
№2
2015