# Jin Matsumoto - nagataki-lab.riken.jpnagataki-lab.riken.jp/workshop/SNGRB2014/ Jin Matsumoto...

date post

29-Jul-2018Category

## Documents

view

221download

1

Embed Size (px)

### Transcript of Jin Matsumoto - nagataki-lab.riken.jpnagataki-lab.riken.jp/workshop/SNGRB2014/ Jin Matsumoto...

Jin MatsumotoRIKEN

Astrophysical Big Bang Laboratory

Numerical Experiments of GRB Jets

Collaborators: Nagataki, Ito, Mizuta, Barkov, Dainotti, Teraki (RIKEN)

Schematic Picture of the GRB JetMeszaros 01

Central engine MHD process

thermal process

Extraction of the rotational energy of the BH or accretion disk through B-field

-annihilation -annihilation fireball

relativistic jet is launched from the central engine associated the death of the massive

star

The dynamics and stability of the relativistic jet is important in order to understand the GRB emissions.

First Simulation of GRB Jet PropagationL120 RELATIVISTIC JETS FROM COLLAPSARS Vol. 531

Fig. 1.Contour maps of the logarithm of the rest-mass density after (a) 3.87 s and (b) 5.24 s and (c) of the Lorentz factor after 5.24 s. The X and Y axesmeasure distance in centimeters. The dashed and solid arcs mark the stellar surface and the outer edge of the exponential atmosphere, respectively. The othersolid line encloses the matter whose radial velocity is greater than 0.3c and whose specific internal energy density is greater than ergs g .19 !15# 10

Fig. 2.Rest-mass density (top) and Lorentz factor (bottom) vs. radius alongthe symmetry axis for model C50 at s (long dashed line), st = 0 t = 1.44(dotted line), s (dashed line), and s (solid line).t = 3.87 t = 5.24

3. RESULTS

3.1. Constant Small-Energy Deposition Rate (Model C50)

For a constant ergs s , a relativistic jet forms50 !1E = 10within a fraction of a second and starts to propagate along therotation axis with a mean speed of cm s (Fig. 1).9 !17.8# 10The jet exhibits all the morphological elements of the Blandford& Rees (1974) jet model in the context of classical doubleradio sources: a terminal bow shock, a narrow cocoon, a contactdiscontinuity separating stellar and jet matter, and a hot spot.Figure 1 shows that the density structure of the star does not

change noticeably during the whole evolution. This, a poste-riori, justifies our neglect of the self-gravity of the star.The propagation of the jet is unsteady because of density

inhomogeneities in the star. The Lorentz factor of the jet, G,increases nonmonotonically with time, while the density dropsto 10!6 g cm (Fig. 2). The density profile shows large var-!3iations (up to a factor of 100) due to internal shock waves.The mean density in the jet is 10!1 g cm . Some of the!3internal biconical shocks, which develop during the jets prop-agation, recollimate the beam. They may provide the internalshocks that are proposed to explain the observed gamma-rayemission. A particularly strong recollimation shock wave(hardly evident at low resolution) forms early in the evolution.A very strong rarefaction wave behind this recollimation shockcauses the largest local acceleration of the beam material, giv-ing rise to a maximum in the Lorentz factor. When the jetencounters a region along the axis where the density gradientis positive (at and ), the jets head islog r 8.1 log r 8.6decelerated, while a central channel in the beam is cleaned byoutflow into the cocoon through the head, which acceleratesthe beam. The combination of both effects (deceleration of thehead and beam acceleration) increases the strength of the in-ternal shocks. Within the jet, the mean value of the specificinternal energy is 10201021 ergs g , or . The mean!1 2O(c )temperature is K (well below the pair creation thresh-85# 10old), implying that the pressure is dominated by radiation inaccordance with our simplified EOS.The relativistic treatment of the hydrodynamics leads to a

qualitatively similar (i.e., formation of a jet) but a quantitativelyvery different evolution than is found in MW99. According totheir Figure 27, the jet propagates 7000 km within the first0.82 s. Furthermore, MW99 infer an asymptotic G of 10 andfind a half-opening angle Q for their jet of 10!. In our sim-ulation, at the same time for the same angular resolution (2!)and , the head reaches a radius of 30,000 km, but the max-Eimum Lorentz factor ( ) is only 4.62 at 12,200 km. SuchGmax

Aloy et al. 2000

inside star: collimation by cocoon outside star: drastic acceleration

adiabatic expansion: energy conversion from to Bernoulli equation:

3D calculation: Zhang+ 2004

0.2 0.4 0.6 0.8 1.0 1.2

-0.2

0.0

0.2Model 2T, 9s

0.2 0.4 0.6 0.8 1.0 1.2

-0.2

0.0

0.2

y (1

0^11

cm

) Model 3A, 9s 0.2 0.4 0.6 0.8 1.0 1.2

-0.2

0.0

0.2Model 3BS, 9s

0.2 0.4 0.6 0.8 1.0 1.2z (10^11 cm)

-0.2

0.0

0.2Model 3BL, 8s

-8.0 -6.5 -5.0 -3.5 -2.0 -0.5 1.0log rho (g/cm^3)

The dynamics of the jet does not drastically change in between 2D and 3D

2D

3D

3D

3D

the initial conditions have cylindrical symmetry

1% asymmetry

10% asymmetry

3D calculation: Zhang+ 2004

0.2 0.4 0.6 0.8 1.0 1.2

-0.2

0.0

0.2

Model 3P3, 8s

0.2 0.4 0.6 0.8 1.0 1.2

-0.2

0.0

0.2

y (1

0^11

cm

)

Model 3P5, 8s

0.2 0.4 0.6 0.8 1.0 1.2z (10^11 cm)

-0.2

0.0

0.2

Model 3P10, 10s

0.2 0.4 0.6 0.8 1.0 1.2

Model 3P3, 8s

0.2 0.4 0.6 0.8 1.0 1.2

Model 3P5, 8s

0.2 0.4 0.6 0.8 1.0 1.2

Model 3P10, 10s

log rho

-8.0 -6.5 -5.0 -3.5 -2.0 -0.5 1.0

gamma

1.0 4.2 7.3 10.5 13.7 16.8 20.0with precession inclination

angle

3

5

10

The jet is decelerated due to the material mixing between the jet and surrounding medium.

3D calculation: Lopez-Camara+ 20132D 3D

The relativistic jet can propagate and break out the progenitor star while remaining relativistic without the dependence of the resolution.

The amount of turbulence and variability observed in the simulations is greater at higher resolutions.

Jet properties are only marginally affected by the dimensionality.

Oscillation Induced RTI and RMI

growth of the Rayleigh-Taylor and Richtmyer-Meshkov instabilities radial oscillation motion of the jet

jet cross section

JM & Masada 2013

These instabilities grow in GRB jet?

3D simulation: propagation of the relativistic jet

focusing on the impact of the oscillation-induced Rayleigh-Taylor and Richtmyer-Meshkov instabilities (JM & Masada 2013) on the 3D jet propagation.

progenitor: 16TI model (Woosley & Heger 06)

specific enthalpy ratio of specific heats Lorentz factor

energyconservation

massconservation

momentum conservation

Basic Equations

3D HD GRB Jet Propagation

Finger-like structures inside star

Finger like structures are excited at the interface between the jet and cocoon.

Finger like structures are excited at the interface between the jet and cocoon.

Finger-like structures inside star

3D simulations: evolution of the cross section of the relativistic jet

(a) Side View

cocoon

shocked ambient medium

bow shock

un-shocked ambient medium

Jet

Zz1 z2 z3

contact discontinuity (CD)

P3

P1P2

reconfinementregion

(b) Top View(b2) Contraction Phase (I)

expanding shock

contracting CD(b3) Contraction Phase (II)

contracting CD

[P1 < P2 < P3][P1 < P2 < P3] [P2 < P3 < P1]

(b1) Expansion Phase

expanding shock

expanding CD

P1P2

P3

[z=z1] [z=z2] [z=z3]

contracting shock

many numerical works in order to investigate the propagation dynamics of the relativistic jet (e.g., Marti+ 97, Aloy+ 00, Zhang+ 03,04, Mizuta+ 06, Perucho+ 08, Morsony+07, Lazzati+ 09, Lopez- Camara+ 13) reconfinement shock in the

collimated jet (Norman et al. 1982; Sanders 1983)

JM& Masada 13

adiabatic cooling

radial oscillating motion and repeated excitation of the reconfinement region (e.g., Gomez+ 97, JM+ 12, Mizuta+ 14)

Growth of RTI in SN Explosion

It is well known that a contact discontinuity, formed where a heavy fluid is supported above a light fluid against gravity, is unstable to perturbations to the boundary between the two fluids.

heavy fluid

light fluid

Although the gravity is negligible in SN explosions, the acceleration of the gas works as a gravity.

, (Chevalier 76)

The condition for the RTI

Growth of RTI in SN Explosion

, (Chevalier 76)

Ono+13

Ebisuzaki+89

The condition for the RTI for a compressible fluid in the absence of the gravity is given by

He/H interface

C+O/He interface

The mechanism of the growth of RTI

In addition to the growth of the Rayleigh-Taylor instability, the growth of the Richtmyer-Meshkov instability is also contributed to the finger like structures.

relativistically hot plasma:

effective inertia:

The effective inertia is important.

(a) Side View

cocoon

shocked ambient medium

bow shock

un-shocked ambient medium

Jet

Zz1 z2 z3

contact discontinuity (CD)

P3

P1P2

reconfinementregion

(b) Top View(b2) Contraction Phase (I)

expanding shock

contracting CD(b3) Contraction Phase (II)

contracting CD

[P1 < P2 < P3][P1 < P2 < P3] [P2 < P3 < P1]

(b1) Expansion Phase

expanding shock

expanding CD

P1P2

P3

[z=z1] [z=z2] [z=z3]

contracting shock

many numerical works in order