We prescribe a method to study the effects of self-gravity of accretion disk around a black hole associated with long Gamma Ray Bursts (GRBs) in an evolving background Kerr metric. This is an extension to our previous work where we presented possible constraints for the final masses and spins of these astrophysical black holes. Incorporating the self-force of the accreting cloud around the black hole is a very important aspect due to the transient nature of the event, in which a huge amount of mass is accreted and changes the fundamental black hole parameters i.e. its mass and spin, during the process. Understanding of the GRBs engine is important because they are possible sources of high-energy particles and gravitational waves as most of the energy released from the dynamical evolution is in the form of gravitational radiation. Here, we describe the analytical framework we developed to employ in our numerical model. The numerical studies are planned for the future work.

We present the relativistic 2D and 3D GRMHD simulation of axisymmetric, inviscid, hydrodynamic accretion flows in a fixed Kerr black hole gravitational field. The flow is having low angular momentum with respect to Keplerian one. A relativistic fluid where its bulk velocity is comparable to the speed of light, flowing in the accretion disk very close to the horizon should be described by adiabtic index: 4=3 < g < 5=3 .The time dependent evolution of shock position and respective effect on mass accretion rate and oscillation frequency with varying adiabatic index has been studied. Here we present some of the results for adiabatic index = 1.4 in a 2D and 3D model.

We compute the evolution of a quasi-spherical, slowly rotating accretion flow around a black hole, whose mass and spin evolve adequately to the mass-energy transfer through the horizon. Our model is relevant for the central engine driving a long gamma ray burst, that originates from the collapse of a massive star. Our results show how much mass and spin a newly formed black hole should possess during collapsar to launch long GRB.

We study the role of adiabatic index in determining the critical points in the transonic low angular momentum accretion flow onto a black hole. We present the general relativistic 2D hydrodynamic simulations of axisymmetric, inviscid accretion flows in a fixed Kerr black hole gravitational field. A relativistic fluid where its bulk velocity is comparable to the speed of light, flowing in the accretion disk very close to the horizon can be described by an adiabatic index of 4/3 < {\gamma} < 5/3. The time dependent evolution of the shock position and respective effects on mass accretion rate and oscillation frequency with varying adiabatic index is discussed in the context of the observed microquasars.

We compute the evolution of a quasi-spherical, slowly rotating accretion flow around a black hole, whose mass and spin evolve adequately to the mass-energy transfer through the horizon. Our model is relevant for the central engine driving a long gamma ray burst, that originates from the collapse of a massive star. The computations of a GRB engine in a dynamically evolving spacetime metric are important specifically due to the transient nature of the event, in which a huge amount of mass is accreted and changes the fundamental black hole parameters, its mass and spin, during the process. We discuss the results in the context of angular momentum magnitude of the collapsing star. We also study the possible formation and evolution of shocks in the envelope, which may temporarily affect accretion. Our results are important for the limitations on the mass and spin range of black holes detected independently by electromagnetic observations of GRBs and gravitational waves. We speculate on the possible constraints for the final masses and spins of these astrophysical black holes. It is shown that the most massive BHs were rahter not formed in a powerful GRB explosion, if the cores of their progenitors were only weakly rotating.