NUMERICAL TREATMENT OF THIN ACCRETION DISK DYNAMICS AROUND ROTATING BLACK HOLES

In the present study, we perform the numerical simulation of a relativistic thin accretion disk around the nonrotating and rapidly rotating black holes using the general relativistic hydrodynamic code with Kerr in Kerr–Schild coordinate that describes the central rotating black hole. Since the high energy X-rays are produced close to the event horizon resulting the black hole–disk interaction, this interaction should be modeled in the relativistic region. We have set up two different initial conditions depending on the values of thermodynamical variables around the black hole. In the first setup, the computational domain is filled with constant parameters without injecting gas from the outer boundary. In the second, the computational domain is filled with the matter which is then injected from the outer boundary. The matter is assumed to be at rest far from the black hole. Both cases are modeled over a wide range of initial parameters such as the black hole angular momentum, adiabatic index, Mach number and asymptotic velocity of the fluid. It has been found that initial values and setups play an important role in determining the types of the shock cone and in designating the events on the accretion disk. The continuing injection from the outer boundary presents a tail shock to the steady state accretion disk. The opening angle of shock cone grows as long as the rotation parameter becomes larger. A more compressible fluid (bigger adiabatic index) also presents a bigger opening angle, a spherical shock around the rotating black hole, and less accumulated gas in the computational domain. While results from [J. A. Font, J. M. A. Ibanez and P. Papadopoulos, Mon. Not. R. Astron. Soc.305 (1999) 920] indicate that the tail shock is warped around for the rotating hole, our study shows that it is the case not only for the warped tail shock but also for the spherical and elliptical shocks around the rotating black hole. The warping around the rotating black hole in our case is much smaller than the one by [J. A. Font, J. M. A. Ibanez and P. Papadopoulos, Mon. Not. R. Astron. Soc.305 (1999) 920], due to the representation of results at the different coordinates. Contrary to the nonrotating black hole, the tail shock is slightly warped around the rotating black hole. The filled computational domain without any injection leads to an unstable accretion disk. However much of it reaches a steady state for a short period of time and presents quasi-periodic oscillation (QPO). Furthermore, the disk tends to loose mass during the whole dynamical evolution. The time-variability of these types of accretion flowing close to the black hole may clarify the light curves in Sgr A*.