Научная статья на тему 'Взрывные волны с космическими лучами'

Взрывные волны с космическими лучами Текст научной статьи по специальности «Физика»

CC BY
186
38
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Астрофизический бюллетень
WOS
Scopus
ВАК
Область наук
Ключевые слова
ударные волны—сверхновые / shock waves—supernovae: general

Аннотация научной статьи по физике, автор научной работы — Б. Арбутина

Взрывные волны присутствуют во многих астрофизических явлениях, таких как сверхновые. В данной статье мы обсуждаем взрывные волны с космическими лучами, т.е. с компонентой со степенной функцией распределения плотности N(p) ∝ p−Γ , которая может быть особенно важной при описании эволюции остатков сверхновых. Мы подтверждаем некоторые предыдущие результаты о том, что существенное количество энергии космических лучей выделяется в направлении центра остатка.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Blast Waves with Cosmic Rays

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−Γ 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.

Текст научной работы на тему «Взрывные волны с космическими лучами»

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 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: arbo@matf.bg.ac.rs

начальные условия взрыва сверхновой, только параметры 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

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)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Совмещая (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)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

d F

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

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

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.