We present fully time-dependent solutions of accretion processes on black holes and compare them with analytical solutions. We use Smoothed Particle Hydrodynamics and Total Variation Diminishing mathods as the numerical tenchniques. Apart from steady state solutions, we also obtain time dependent results where the shock oscillates with the characteristic period of quasi-periodic oscillations (QPOs) observed in galactic and extragalactic black hole candidates. This oscillation produces 10-15% variation in the hard X-rays in galactic candidates and in soft X-rays in Extragalactic candidates.

We extend our previous numerical simulation of accretion disks with shock waves when cooling effects are also included. We consider bremsstrahlung and other power law processes: $\Lambda \propto T^{\alpha} \rho^2$ to mimic cooling in our simulation. We employ {\it Smoothed Particle Hydrodynamics} technique as in the past. We observe that for a given angular momentum of the flow, the shock wave undergoes a steady, radial oscillation with the period is roughly equal to the cooling time. Oscillations seem to take place when the disk and cooling parameters (i.e., accretion rate, cooling process) are such that the infall time from shock is of the same order as the post-shock cooling time. The amplitude of oscillation could be up to ten percent of the distance of the shock wave from the black hole when the black hole is accreting. When the accretion is impossible due to the centrifugal barrier, the amplitude variation could be much larger. Due to the oscillation, the energy output from the disk is also seen to vary quasi-periodically. We believe that these oscillations might be responsible for the quasi periodic oscillation (QPO) behaviors seen in several black hole candidates, in neutron star systems as well as dwarf novae outbursts such as SS Cygni and VW Hyi.

Results of numerical simulations of shock solutions in a geometrical thin accretion disk around a Kerr black hole (BH) are presented. Using the smoothed particle hydrodynamics (SPH) technique, the influence of the central object is included by means of an effective potential, We first present the theory of standing shock formation in accretion disks around a Kerr black hole, and show that the results of our numerical simulation agree very well with the theoretical results. We find that the shocks in an inviscid flow are very stable. We also remove the ambiguity prevalent regarding the location and stability of shocks in adiabatic flows. Finally we sketch some of the astrophysical consequences of our findings in relation to accretion disks in Active Galactic Nuclei (AGN) and Quasars.