Numerical simulation of rotating accretion disk around the Schwarzschild black hole using GRH code

The 2D time dependent solution of thin accretion disk in a close binary system have been presented on the equatorial plane around the Schwarzschild black hole. To do that, the special part of the general relativistic hydrodynamical (GRH) equations are solved using high resolution shock capturing (HRSC) schemes. The spiral shock waves on the accretion disk are modeled using perfect fluid equation of state with adiabatic indices γ = 1.05, 1.2 and 5/3. The results show that the spiral shock waves are created for gammas except the case γ = 5/3. These results are consistent with the results from Newtonian hydrodynamic code except those which are close to black hole. Newtonian approximation does not give good solution when the matter is close to the black hole. Our simulations illustrate that the spiral shock waves are created close to black hole and the location of inner radius of spiral shock wave is around 6M and it depends on the specific heat rates. We also find that the smaller γ is the more tightly packed in the spiral winds.