We have done nth order perturbation analysis of the Navier-Stokes equation with the presence of turbulent viscosity in context of thin accretion flow around a black hole. In order to find the stability criteria, we used Green's function to solve the nth order velocity perturbation equation. Finally, we have proposed a velocity stability criterion which helps to find whether the base flow to be stable for higher order perturbation.

Entangled magnetic fields entering into an accretion flow would very soon be stretched into a dominant toroidal component due to strong differentially rotating motion inside the accretion disc. This is particularly true for weakly viscous, low angular momentum transonic or advective discs. We study the trajectories of toroidal flux tubes inside a geometrically thick flow which undergoes a centrifugal force supported shock. We also study effects of these flux tubes on the dynamics of the inflow and the outflow. We use a finite difference method (Total Variation Diminishing) for this purpose and specifically focussed on whether these flux tubes significantly affect the properties of the outflows such as its collimation and the rate. It is seen that depending upon the cross-sectional radius of the flux tubes which control the drag force, these field lines may move towards the central object or oscillate vertically before eventually escaping out of the funnel wall (pressure zero surfaces) along the vertical direction. A comparison of results obtained with and without flux tubes show these flux tubes could play pivotal role in collimation and acceleration of jets and outflows.

We study time evolution of rotating, axisymmetric, two dimensional inviscid accretion flows around black holes using a grid based finite difference method. We do not use reflection symmetry on the equatorial plane in order to inspect if the disk along with the centrifugal barrier oscillated vertically. In the inviscid limit, we find that the CENtrifugal pressure supported BOundary Layer (CENBOL) is oscillating vertically, more so, when the specific angular momentum is higher. As a result, the rate of outflow produced from the CENBOL, also oscillates. Indeed, the outflow rates in the upper half and the lower half are found to be anti-correlated. We repeat the exercise for a series of specific angular momentum {\lambda} of the flow in order to demonstrate effects of the centrifugal force on this interesting behaviour. We find that, as predicted in theoretical models of disks in vertical equilibrium, the CENBOL is produced only when the centrifugal force is significant and more specifically, when {\lambda} > 1.5. Outflow rate itself is found to increase with {\lambda} as well and so is the oscillation amplitude. The cause of oscillation appears to be due to the interaction among the back flow from the centrifugal barrier, the outflowing winds and the inflow. For low angular momentum, the back flow as well as the oscillation are missing. To our knowledge, this is the first time that such an oscillating solution is found with an well-tested grid based finite difference code and such a solution could be yet another reason of why Quasi-Periodic Oscillations should be observed in black hole candidates which are accreting low angular momentum transonic flows.

In this paper, we present the time evolution of a rotationally axisymmetric gas ring around a non rotating black hole using two dimensional grid-based hydrodynamic simulation. We show the way in which angular momentum transport is included in simulations of non-self-gravitating accretion of matter towards a black hole. We use the Shakura-Sunyaev {\alpha} viscosity prescription to estimate the turbulent viscosity for all major viscous stress tensors. We investigate how a gas ring which is initially assumed to rotate with Keplerian angular velocity is accreted on to a black hole and hence forms accretion disc in the presence of turbulent viscosity. We show that the centrifugal pressure supported sub-Keplerian flow with shocks forms when the ring starts to disperse with inclusion of relatively small amount of viscosity. But, if the viscosity is above the critical value, the shock disappears altogether and the whole disc becomes Kepleiran which is subsonic everywhere except in a region close to the horizon, where it supersonically enters to the black hole. We discovered a multiple valued Mach number solution and the corresponding density distributions that connects matter (a) from the initial Keplerian gas ring to a sub-Keplerian disc with shocks in presence of small amount of viscosity and (b) from the sub-Keplerian flow to a Keplerian disc in presence of huge amount of viscosity. We calculate the temporal variations of the magnitude of various time scales which ensure us about the stability of the flow.

A black hole accretion is necessarily transonic. In presence of sufficiently high viscosity and cooling effects, a low-angular momentum transonic flow can become a standard Keplerian disc except close to the where hole where it must pass through the inner sonic point. However, if the viscosity is not high everywhere and cooling is not efficient everywhere, the flow cannot completely become a Keplerian disc. In this paper, we show results of rigorous numerical simulations of a transonic flow having vertically varying viscosity parameter (being highest on the equatorial plane) and optical depth dependent cooling processes to show that the flow indeed segregates into two distinct components as it approaches a black hole. The component on the equatorial plane has properties of a standard Keplerian disc, though the flow is not truncated at the inner- most stable circular orbit. This component extends till the horizon as a sub-Keplerian flow. This standard disc is found to be surrounded by a hot, low angular momentum component forming a centrifugal barrier dominated oscillating shock wave, consistent with the Chakrabarti-Titarchuk two component advective flow configuration.

Two component advective flows are the most physical accretion disks which arise from theoretical consideration. Since viscosity is the determining factor, we investigate the effects of viscous stresses on accretion flows around a nonrotating black hole. As a consequence of angular momentum transfer by viscosity in an accretion flow, the angular momentum distribution is modified. We include cooling effects and found that a Keplerian disk is produced on the equatorial plane and the flow above and below remains sub Keplerian. This gives a complete picture to date, of the formation of a Two component advective flow around a black hole.

We carry out numerical simulation of accretion flows around a black hole in presence of both viscous and cooling effects. Instead of using a constant $\alpha$ parameter throughout the simulation grid, we assume that $\alpha$ is maximum on the equatorial plane and gradually goes down in the perpendicular direction We show that when the injected sub-Keplerian flow angular momentum is high enough and/or the viscosity and also cooling is high enough, Two Component {\bf Accretion} Flow (TCAF) would be formed, otherwise the sub-Keplerian flow would remain sub-Keplerian. We see that a Keplerian disk is produced on the equatorial plane. Time variations of the total, Keplerian and Sub-Keplerian matter are studied with respect to various flow parameters.

We carry out a series of numerical simulations of viscous accretion flows having a reasonable spatial distribution of the viscosity parameter. We add the power-law cooling throughout the flow. We show that, in agreement with the theoretical solutions of viscous transonic flows, matter having the viscosity parameter above a critical value becomes a Keplerian disk while matter having lesser viscosity remains a low angular momentum, sub-Keplerian flow. The latter component produces centrifugal pressure supported shock waves. Thus, for instance, a flow having sufficiently high viscosity on the equatorial plane and low viscosity above and below, would produce a Two Component Advective Flow (TCAF) where a Keplerian disk is surrounded by a rapidly infalling sub-Keplerian halo. We find that the post- shock region of the relatively cooler Keplerian disk is evaporated and the overall configuration is quite stable. This agrees with the theoretical model with two components which attempt to explain the spectral and timing properties of black hole candidates.

We study the time evolution of a rotating, axisymmetric, viscous accretion flow around black holes using a grid based finite difference method. We use the Shakura-Sunyaev viscosity prescription. However, we compare with the results obtained when all the three independent components of the viscous stress are kept. We show that the centrifugal pressure supported shocks became weaker with the inclusion of viscosity. The shock is formed farther out when the viscosity is increased. When the viscosity is above a critical value, the shock disappears altogether and the flow becomes subsonic and Keplerian everywhere except in a region close to the horizon, where it remains supersonic. We also find that as the viscosity is increased, the amount of outflowing matter in the wind is decreased to less than a percentage of the inflow matter. Since the post-shock region could act as a reservoir of hot electrons or the so-called 'Compton cloud', the size of which changes with viscosity, the spectral properties are expected to depend on viscosity strongly: the harder states are dominated by low-angular momentum and the low viscosity flow with a significant outflows while the softer states are dominated by the high viscosity Keplerian flow having very little outflows.

We carry out a time dependent numerical simulation where both the hydrodynamics and the radiative transfer are coupled together. We consider a two-component accretion flow in which the Keplerian disk is immersed inside an accreting low angular momentum flow (halo) around a black hole. The injected soft photons from the Keplerian disk are reprocessed by the electrons in the halo. We show that in presence of an axisymmetric soft-photon source, the spherically symmetric Bondi flow losses its symmetry and becomes axisymmetric. The low angular momentum flow was observed to slow down close to the axis and formed a centrifugal barrier which added new features into the spectrum. Using the Monte Carlo method, we generated the radiated spectra as functions of the accretion rates. We find that the transitions from a hard state to a soft state is determined by the mass accretion rates of the disk and the halo. We separate out the signature of the bulk motion Comptonization and discuss its significance. We study how the net spectrum is contributed by photons suffering different number of scatterings and spending different amounts of time inside the Compton cloud. We study the directional dependence of the emitted spectrum as well.

We study the accretion processes on a black hole by numerical simulation. We use a grid based finite difference code for this purpose. We scan the parameter space spanned by the specific energy and the angular momentum and compare the time-dependent solutions with those obtained from theoretical considerations. We found several important results (a) The time dependent flow behaves close to a constant height model flow in the pre-shock region and a flow with vertical equilibrium in the post-shock region. (c) The infall time scale in the post-shock region is several times higher than the free-fall time scale. (b) There are two discontinuities in the flow, one being just outside of the inner sonic point. Turbulence plays a major role in determining the locations of these discontinuities. (d) The two discontinuities oscillate with two different frequencies and behave as a coupled harmonic oscillator. A Fourier analysis of the variation of the outer shock location indicates higher power at the lower frequency and lower power at the higher frequency. The opposite is true when the analysis of the inner shock is made. These behaviours will have implications in the spectral and timing properties of black hole candidates.