ABSTRACT

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor g = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are < c/\sqrt 3 in the weakly magnetized simulations and < 0.3c in the strongly magnetized simulations. The Alfven wave speed is < 0.07c in the weakly magnetized simulations and < 0.56c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

ABSTRACT

We have performed several simulations of black hole systems (non-rotating, black hole spin parameter a=0.0 and rapidly rotating, a=0.95) with a geometrically thin Keplerian disk using the newly developed RAISHIN code. The simulation results show the formation of jets driven by the Lorentz force and the gas pressure gradient. The jets have mildly relativistic speed (~0.4c). The matter is continuously supplied from the accretion disk and the jet propagates outward until each applicable terminal simulation time (non-rotating: t/tau_S = 275 and rotating: t/tau_S=200, $tau_S \equiv r_S/c$). It appears that a rotating black hole creates an additional, faster, and more collimated matter-dominated inner outflow (~0.5c) formed and accelerated by the twisted magnetic field resulting from frame-dragging in the black hole ergosphere. This is the first known simulation confirming the formation of an inner magnetically-driven, matter-dominated jet by the frame-dragging effect from a black hole co-rotating with a thin Keplerian disk threaded by a vertical magnetic field. This result indicates that jet kinematic structure depends on black hole rotation and on the initial magnetic field configuration and strength.

RAISHIN: A High-Resolution Three-Dimensional General Relativistic Magnetohydrodynamics Code

Yosuke Mizuno, Ken-Ichi Nishikawa, Shinji Koide, Philip Hardee, & Gerald, J. Fishman

Submitted to ApJS

ABSTRACT

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.