Two-dimensional radiation-hydrodynamic model for limit-cycle oscillations of luminous accretion disks

We investigate the time evolution of luminous accretion disks around black holes, conducting the two-dimensional radiation-hydrodynamic simulations. We adopt the alpha prescription for the viscosity. The radial-azimuthal component of viscous stress tensor is assumed to be proportional to the total pressure in the optically thick region, while the gas pressure in the optically thin regime. The viscosity parameter, alpha, is taken to be 0.1. We find the limit-cycle variation in luminosity between high and low states. When we set the mass input rate from the outer disk boundary to be 100 L_E/c^2, the luminosity suddenly rises from 0.3L_E to 2L_E, where L_E is the Eddington luminosity. It decays after retaining high value for about 40 s. Our numerical results can explain the variation amplitude and duration of the recurrent outbursts observed in microquasar, GRS 1915+105. We show that the multi-dimensional effects play an important role in the high-luminosity state. In this state, the outflow is driven by the strong radiation force, and some part of radiation energy dissipated inside the disk is swallowed by the black hole due to the photon-trapping effects. This trapped luminosity is comparable to the disk luminosity. We also calculate two more cases: one with a much larger accretion rate than the critical value for the instability and the other with the viscous stress tensor being proportional to the gas pressure only even when the radiation pressure is dominant. We find no quasi-periodic light variations in these cases. This confirms that the limit-cycle behavior found in the simulations is caused by the disk instability