We investigate the behaviour of low angular momentum viscous accretion flows around black holes using Smooth Particle Hydrodynamics (SPH) method. Earlier, it has been observed that in a significant part of the energy and angular momentum parameter space, rotating transonic accretion flow undergoes shock transition before entering in to the black hole and a part of the post-shock matter is ejected as bipolar outflows, which are supposed to be the precursor of relativistic jets. In this work, we simulate accretion flows having injection parameters from the inviscid shock parameter space, and study the response of viscosity on them. With the increase of viscosity, shock becomes time dependent and starts to oscillate when the viscosity parameter crosses its critical value. As a result, the in falling matter inside the post-shock region exhibits quasi-periodic variations and causes periodic ejection of matter from the inner disc as outflows. In addition, the same hot and dense post-shock matter emits high energy radiation and the emanating photon flux also modulates quasi-periodically. Assuming a ten solar mass black hole, the corresponding power density spectrum peaks at the fundamental frequency of few Hz followed by multiple harmonics. This feature is very common in several outbursting black hole candidates. We discuss the implications of such periodic variations.

We examine the low angular momentum flow model for Sgr A* using two-dimensional hydrodynamical calculations based on the parameters of the specific angular momentum and total energy estimated in the recent analysis of stellar wind of nearby stars around Sgr A*. The accretion flow with the plausible parameters is non-stationary and an irregularly oscillating shock is formed in the inner region of a few tens to a hundred and sixty Schwarzschild radii. Due to the oscillating shock, the luminosity and the mass-outflow rate are modulated by several per cent to a factor of 5 and a factor of 2-7, respectively, on time-scales of an hour to ten days. The flows are highly advected and the radiative efficiency of the accreting matter into radiation is very low, 10^{-5}--$10^{-3}, and the input accretion rate of 4.0* 10^{-6} solar mass/yr results in the observed luminosities -- 10^{36} erg/s of Sgr A* if a two-temperature model and the synchrotron emission are taken into account. The mass-outflow rate of the gas originating in the post-shock region increases with the increasing input specific angular momentum and ranges from a few to 99 per cent of the input accreting matter, depending on the input angular momentum. The oscillating shock is necessarily triggered if the specific angular momentum and the specific energy belong to or are located just nearby in the range of parameters responsible for a stationary shock in rotating inviscid and adiabatic accretion flow. The time variability may be relevant to the flare activity of Sgr A*.

Current theories on planetary formation establish that giant planet formation should be contextual to their quick migration towards the central star due to the protoplanets-disc interactions on a timescale of the order of $10^5$ years, for objects of nearly 10 terrestrial masses. Such a timescale should be smaller by an order of magnitude than that of gas accretion onto the protoplanet during the hierarchical growing-up of protoplanets by collisions with other minor objects. These arguments have recently been analysed using N-body and/or fluid-dynamics codes or a mixing of them. In this work, inviscid 2D simulations are performed, using the SPH method, to study the migration of one protoplanet, to evaluate the effectiveness of the accretion disc in the protoplanet dragging towards the central star, as a function of the mass of the planet itself, of disc tangential kinematics. To this purpose, the SPH scheme is considered suitable to study the roles of turbulence, kinematic and boundary conditions, due to its intrinsic advective turbulence, especially in 2D and in 3D codes. Simulations are performed both in disc sub-Keplerian and in Keplerian kinematic conditions as a parameter study of protoplanetary migration if moderate and consistent deviations from Keplerian Kinematics occur. Our results show migration times of a few orbital periods for Earth-like planets in sub-Keplerian conditions, while for Jupiter-like planets estimates give that about $10^4$ orbital periods are needed to half the orbital size. Timescales of planet migration are strongly dependent on the relative position of the planet with respect to the shock region near the centrifugal barrier of the disc flow.

We examine the jets and the disc of SS 433 at super-Eddington luminosities with 600 times Eddington critical accretion rate by time-dependent two-dimensional radiation hydrodynamical calculations, assuming alpha-model for the viscosity. One-dimensional supercritical accretion disc models with mass loss or advection are used as the initial configurations of the disc. As a result, from the initial advective disc models with alpha =0.001 and 0.1, we obtain the total luminosities 2.5x10^{40} and 2.0x10^{40} erg/s. The total mass-outflow rates are 4x10^{-5} and 10^{-4} solar-mass/yr and the rates of the relativistic axial outflows in a small half opening angle of 1 degree are about 10^{-6} solar-mass/yr: the values generally consistent with the corresponding observed rates of the wind and the jets, respectively. From the initial models with mass loss but without advection, we obtain the total mass-outflow and axial outflow rates smaller than or comparable to the observed rates of the wind and the jets respectively, depending on alpha. In the advective disc model with alpha=0.1, the initially radiation-pressure dominant, optically thick disc evolves to the gas-pressure dominated, optically thin state in the inner region of the disc, and the inner disc is unstable. Consequently, we find remarkable modulations of the disc luminosity and the accretion rate through the inner edge. These modulations manifest themselves as the recurrent hot blobs with high temperatures and low densities at the disc plane, which develop outward and upward and produce the QPOs-like variability of the total luminosity with an amplitude of a factor of 2 and quasi-periods of 10 -- 25 s. This may explain the massive jet ejection and the QPOs phenomena observed in SS 433.

The heating of the ion-neutral (or ambipolar) diffusion may affect the thermal phases of the molecular clouds. We present an investigation on the effect of this heating mechanism in the thermal instability of the molecular clouds. A weakly ionized one dimensional slab geometry, which is allowed for self-gravity and ambipolar diffusion, is chosen to study its thermal phases. We use the thermodynamic evolution of the slab to obtain the regions where slab cloud becomes thermally unstable. We investigate this evolution using the model of ambipolar diffusion with two-fluid smoothed particle hydrodynamics, as outlined by Hosking & Whitworth. Firstly, some parts of the technique are improved to test the pioneer works on behavior of the ambipolar diffusion in an isothermal self-gravitating slab. Afterwards, the improved two-fluid technique is used for thermal evolution of the slab. The results show that the thermal instability may persist inhomogeneities with a large density contrast at the intermediate parts of the cloud. We suggest that this feature may be responsible for the planet formation in the intermediate regions of a collapsing molecular cloud and/or may also be relevant to the formation of star forming dense cores in the clumps.

We present an investigation on effect of the ion-neutral (or ambipolar) diffusion heating rate on thermal phases of a molecular cloud. We use the modeling of ambipolar diffusion with two-fluid smoothed particle hydrodynamics, as discussed by Nejad-Asghar & Molteni. We take into account the ambipolar drift heating rate on the net cooling function of the molecular clouds, and we investigate the thermal phases in a self-gravitating magnetized one dimensional slab. The results show that the isobaric thermal instability criterion is satisfied in the outer parts of the cloud, thus, these regions are thermally unstable while the inner part is stable. This feature may be responsible for the planet formation in the outer parts of a collapsing molecular cloud and/or may also be relevant for the formation of star forming dense cores in the clumps.

We revisit the modeling of ion-neutral (or ambipolar) diffusion with two fluid smoothed particle hydrodynamics, as discussed by Hosking & Whitworth. Some parts of the technique are optimized to testify the pioneer works on behavior of the ambipolar diffusion in an isothermal self-gravitating layer. The frictional heating by ambipolar diffusion is examined, and its effect on fragmentation of the layer is studied. The results are compared to the thermal phases of instability as obtained by Nejad-Asghar.

We numerically examine centrifugally supported shock waves in 2D rotating accretion flows around a stellar-mass (10M_sun) and a supermassive (10^6M_sun) black holes over a wide range of input accretion rates of 10^7 >\dot M/\dot M_E>10^{-4}. The resultant 2D-shocks are unstable with time and the luminosities show quasi-periodic oscillations (QPOs) with modulations of a factor of 2-3 andwith periods of a tenth seconds to several hours, depending on the black hole masses. The shock oscillation model may explain the intermediate frequency QPOs with 1-10 Hz observed in the stellar-mass black hole candidates and also suggest the existence of QPOs with the period of hours in AGNs. When the accretion rate is low, the luminosity increases in proportion to the accretion rate. However, when the accretion rate exceeds greatly the Eddington critical rate \dot M_E, the luminosity is insensitive to the accretion rate and is kept constantly around 3 L_E. On the other hand, the mass-outflow rate increases in proportion to the accretion rate and it amounts to about a few percent of theinput mass-flow rate.

We present a simple physical model of the central source emission in the M87 galaxy. It is well known that the observed X-ray luminosity from this galactic nucleus is much lower than the predicted one, if a standard radiative efficiency is assumed. Up to now the main model invoked to explain such a luminosity is the ADAF (Advection-Dominated-Accretion-Flow) model. Our approach supposes only a simple axis-symmetric adiabatic accretion with a low angular momentum together with the bremsstrahlung emission process in the accreting gas. With no other special hypothesis on the dynamics of the system, this model agrees well enough with the luminosity value measured by Chandra.

We discuss a special case of formation of axisymmetric shocks in the accretion flow of ideal gas onto a Schwarzschild black hole: when the total energy of the flow is negative. The result of our analysis enlarges the parameter space for which these steady shocks are exhibited in the accretion of gas rotating around relativistic stellar objects. Since keplerian disks have negative total energy, we guess that, in this energy range, the production of the shock phenomenon might be easier than in the case of positive energy. So our outcome reinforces the view that sub-keplerian flows of matter may significantly affect the physics of the high energy radiation emission from black hole candidates. We give a simple procedure to obtain analytically the position of the shocks. The comparison of the analytical results with the data of 1D and 2D axisymmetric numerical simulations confirms that the shocks form and are stable.

We present the simulation of 3D time dependent flow of rotating ideal gas falling into a Schwarzschild black hole. It is shown that also in the 3D case steady shocks are formed in a wide range of parameters (initial angular momentum and thermal energy). We therefore highlight the stability of the phenomenon of shock formation in sub keplerian flows onto black holes, and reenforce the role of the shocks in the high luminosity emission from black hole candidates. The simulations have been performed using a parallelized code based on the Smoothed Particles Hydrodynamics method (SPH). We also discuss some properties of the shock problem that allow its use as a quantitative test of the accuracy of the used numerical method. This shows that the accuracy of SPH is acceptable although not excellent.

We present results of several numerical simulations of two dimensional advective flows which include cooling processes. We show that the computed light curve is similar to the $\chi$ state in GRS 1915+105. The power density spectrum (PDS) also shows presence of QPOs near the break frequency.

Super-Eddington accretion discs with 3 and 15 dot M_E around black holes with mass 10 M_sun are examined by two-dimensional radiation hydrodynamical calculations extending from the inner disc edge to 5*10^4 r_g and lasting up to \sim 10^6 r_g/c. The dominant radiation-pressure force in the inner region of the disc accelerates the gas vertically to the disc plane, and jets with 0.2 -- 0.4$c$ are formed along the rotational axis. In the case of the lower accretion rate, the initially anisotropic high-velocity jet expands outward and becomes gradually isotropic flow in the distant region. The mass-outflow rate from the outer boundary is as large as \sim 10^{19} -- 10^{23} g s^{-1}, but it is variable and intermittent with time; that is, the outflow switches occasionally to inflow in the distant region. The luminosity also varies as \sim 10^{40} -- 10^{42} erg s^{-1} on a long time-scale. On the other hand, the jet in the case of the higher accretion rate maintains its initial anisotropic shape even after it goes far away. The mass-outflow rate and the luminosity attain to steady values of 3*10^{19} g s^{-1} and 1.3*10^{40} erg s^{-1}, respectively. In accordance with the local analysis of the slim accretion disc model, the disc is thermally unstable in the case of 3 \dot M_E} but stable in the case of 15 \dot M_E. The super-Eddington model with 15 \dot M_E is promising to explain a small collimation degree of the jet and a large mass-outflow rate observed in the X-ray source SS 433.

It is well known that the rotating inviscid accretion flows with adequate injection parameters around black holes could form shock waves close to the black holes, after the flow passes through the outer sonic point and can be virtually stopped by the centrifugal force. We examine numerically such shock waves in 1D and 2D accretion flows, taking account of cooling and heating of the gas and radiation transport. The numerical results show that the shock location shifts outward compared with that in the adiabatic solutions and that the more rarefied ambient density leads to the more outward shock location. In the 2D-flow, we find an intermediate frequency QPO behavior of the shock location as is observed in the black hole candidate GRS 1915+105.

We show that accretion disks, both in the subcritical and supercritical accretion rate regime, may exhibit significant amplitude luminosity oscillations. The luminosity time behavior has been obtained by performing a set of time-dependent 2D SPH simulations of accretion disks with different values of alpha and accretion rate. In this study, to avoid any influence of the initial disk configuration, we produced the disks injecting matter from an outer edge far from the central object. The period of oscillations is 2 - 50 s respectively for the two cases, and the variation amplitude of the disc luminosity is 10^38 - 10^39 erg/s. An explanation of this luminosity behavior is proposed in terms of limit cycle instability: the disk oscillates between a radiation pressure dominated configuration (with a high luminosity value) and a gas pressure dominated one (with a low luminosity value). The origin of this instability is the difference between the heat produced by viscosity and the energy emitted as radiation from the disk surface (the well-known thermal instability mechanism). We support this hypothesis showing that the limit cycle behavior produces a sequence of collapsing and refilling states of the innermost disk region.

We present the results of several numerical simulations of two dimensional axi-symmetric accretion flows around black holes using Smoothed Particle Hydrodynamics (SPH) in the presence of cooling effects. We consider both stellar black holes and super-massive black holes. We observe that due to {\it both radial and vertical oscillation} of shock waves in the accretion flow, the luminosity and average thermal energy content of the inner disk exhibit very interesting behaviour. When power density spectra are taken, quasi-periodic variabilities are seen at a few Hz and also occasionally at hundreds of Hz for stellar black holes. For super-massive black holes, the time scale of the oscillations ranges from hours to weeks. The power density spectra have a flat top behavior with average {\it rms} amplitude of a few percent and a broken power-law behavior. The break frequency is generally found to be close to the frequency where the shock oscillates.

We show that a standard Shakura-Sunyaev accretion disc around a black hole with an accretion rate lower than the critical Eddington limit does show the instability in the radiation pressure dominated zone. We obtain this result performing time-dependent simulations of accretion disks for a set of values of the viscosity parameter and accretion rate. In particular we always find the occurrence of the collapse of the disc: the instability develops always towards a collapsed gas pressure dominated disc and not towards the expansion. This result is valid for all initial configurations we tested. We find significant convective heat flux that increases the instability development time, but is not strong enough to inhibit the disc collapse. A physical explanation of the lack of the expansion phase is proposed considering the role of the radial heat advection. Our finding is relevant since it excludes the formation of the hot comptonizing corona -often suggested to be present- around the central object by the mechanism of the Shakura-Sunyaev instability. We also show that, in the parameters range we simulated, accretion disks are crossed by significant amplitude acoustic waves.

We present results of several numerical simulations of two dimensional axi-symmetric accretion flows around black holes using the Smoothed Particle Hydrodynamics (SPH). We consider both stellar black holes and as well as super-massive black holes. We assume bremsstrahlung to be the only source of cooling as it is simpler to implement in numerical simulations. We observe that due to both radial and vertical oscillation of shock waves in the accretion flow, the luminosity and average thermal energy content of the inner disk exhibit very interesting behaviors. When power density spectra are taken, quasi-periodic oscillations are seen at a few Hz and also occasionally at hundreds of Hz for stellar black holes. For super-massive black holes, the time-scale of the oscillations ranges from hours to weeks. The power density spectra have usual flat top behavior with average {\it rms} amplitude a few percent and a broken power-law behavior. The break frequency is generally found to be close to the QPO frequency as seen in the observed power spectra of black holes.

We demonstrate that generically an accretion disk around a compact object could have a new type of instability in that the accretion flow need not be symmetric with respect to the equatorial plane even when matter is supplied symmetrically farther out. We find that this behaviour is mainly due to interaction of outgoing matter bounced off the centrifugal barrier and the incoming accretion. We believe that X-ray variability could be developed by this instability.

Using simple arguments we show that the canonical thin keplerian accretion disks cannot smoothly match any plain advection dominated flow (ADAF) model. By 'plain' ADAF model we mean the ones with zero cooling. The existence of sonic points in exact solutions is critical and imposes constraints that cannot be surpassed adopting 'reasonable' physical conditions at the hypothetical match point. Only the occurrence of new critical physical phenomena may produce a transition. We propose that exact advection models are a class of solutions which don't necessarily involve the standard thin cool disks and suggest a different scenario in which good ADAF solutions could eventually occur.

We present the results of systematic comparison of 2D and 3D numerical models of mass transfer in semidetached binaries. It is shown that only for case of \gamma \sim 1 (near-isothermal case) the obtained 2D and 3D solutions are qualitatively similar. For higher value of \gamma the 3D flow structure is drastically changed and for \gamma=1.2 the accretion disk is not formed in the system. Numerical results show that for case of \gamma=1.2 the 2D and 3D solutions are different. We discuss the spiral-shaped shock waves obtained in numerical models. It is shown that these shocks are not intrinsic spirals and are caused by collisions of the gaseous flows in the system.

Analytical studies and numerical simulations of time dependent axially symmetric flows onto black holes have shown that it is possible to produce stationary shock waves with a stable position both for ideal inviscid and for moderately viscous accretion disks. We perform several two dimensional numerical simulations of accretion flows in the equatorial plane to study shock stability against non-axisymmetric azimuthal perturbations. We find a peculiar new result. A very small perturbation seems to produce an instability as it crosses the shock, but after some small oscillations, the shock wave suddenly transforms into an asymmetric closed pattern, and it stabilizes with a finite radial extent, despite the inflow and outflow boundary conditions are perfectly symmetric. The main characteristics of the final flow are: 1) The deformed shock rotates steadily without any damping. It is a permanent feature and the thermal energy content and the emitted energy vary periodically with time. 2) This behavior is also stable against further perturbations. 3) The average shock is still very strong and well defined, and its average radial distance is somewhat larger than that of the original axially symmetric circular shock. 4) Shocks obtained with larger angular momentum exhibit more frequencies and beating phenomena. 5) The oscillations occur in a wide range of parameters, so this new effect may have relevant observational consequences, like (quasi) periodic oscillations, for the accretion of matter onto black holes. Typical time scales for the periods are 0.01 and 1000 seconds for black holes with 10 and 1 million solar mass, respectively.

The results of 3D hydrodynamic simulation of mass transfer in semidetached binaries of different types (cataclysmic variables and low-mass X-ray binaries) are presented. We find that taking into account of a circumbinary envelope leads to significant changes in the stream-disc morphology. In particular, the obtained steady-state self-consistent solutions show an absence of impact between gas stream from the inner Lagrangian point L1 and forming accretion disc. The stream deviates under the action of gas of circumbinary envelope, and does not cause the shock perturbation of the disc boundary (traditional `hotspot'). At the same time, the gas of circumbinary envelope interacts with the stream and causes the formation of an extended shock wave, located on the stream edge. We discuss the implication of this model without `hotspot' (but with a shock wave located outside the disc) for interpretation of observations. The comparison of synthetic light curves with observations proves the validity of the discussed hydrodynamic model without `hotspot'. We also consider the influence of a circumbinary envelope on the mass transfer rate in semidetached binaries. The obtained features of flow structure in the vicinity of L1 show that the gas of circumbinary envelope plays an important role in the flow dynamics, and that it leads to significant (in order of magnitude) increasing of the mass transfer rate. The comparison of gaseous flows structure obtained in 2D and 3D approaches is presented. We discuss the common features of the flow structures and the possible reasons of revealed differences.

Viscous Keplerian discs become sub-Keplerian close to a black hole since they pass through sonic points before entering into it. We study the time evolution of polytropic viscous accretion discs (both in one and two dimensional flows) using Smoothed Particle Hydrodynamics. We discover that for a large region of the parameter space, when the flow viscosity parameter is less than a critical value, standing shock waves are formed. If the viscosity is very high then the shock disappears. In the intermediate viscosity the disc oscillates very significantly in viscous time-scale. Our simulations indicate that these centrifugally supported high density region close to a black hole plays an active role in the flow dynamics, and consequently, the radiation dynamics.

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 characterize the nature of thin, axisymmetric, inviscid, accretion flows of cold adiabatic gas with zero specific energy in the vicinity of a black hole by the specific angular momentum. Using two-dimensional hydrodynamic simulations in cylindrical geometry, we present various regimes in which the accretion flows behave distinctly differently. When the flow has a small angular momentum $(\lambda\lsim\lambda_b)$, most of the material is accreted into the black hole forming a quasi-spherical flow or a simple disk-like structure around it. When the flow has a large angular momentum (typically, larger than the marginally bound value, $\lambda\gsim\lambda_{mb}$), almost no accretion into the black hole occurs. Instead, the flow produces a stable standing shock with one or more vortices behind it and is deflected away at the shock as a conical outgoing wind of higher entropy. If the flow has an angular momentum somewhat smaller than $\lambda_{mb}$ $(\lambda_{u}\lsim\lambda\lsim\lambda_{mb})$, a fraction (typically, $5-10$\%) of the incoming material is accreted into the black hole, but the the flow structure formed is similar to that as for $\lambda\gsim\lambda_{mb}$. Some of the deflected material is accreted back into the black hole, while the rest is blown away as an outgoing wind. These two cases with $\lambda\gsim\lambda_u$ correspond those studied in the previous works by Molteni, Lanzafame, \& Chakrabarti (1994) and Ryu \etal (1995). However, the flow with an angular momentum close to the marginally stable value $(\lambda_{ms})$ is found to be unstable. More specifically, if $\lambda_b\lsim\lambda\sim\lambda_{ms}\lsim\lambda_u$, the flow displays a distinct periodicity in the sense that the inner part of the disk is built and

We compare the results of numerical simulations of thin and quasi-spherical (thick) accretion flows with existing analytical solutions. We use a Lagrangian code based on the Smooth Particle Hydrodynamics (SPH) scheme and an Eulerian finite difference code based on the Total Variation Diminishing (TVD) scheme. In one-dimensional thin flows, the results of the simulations, with or without shocks, agree very well with each other and with analytical solutions. In two-dimensional thick flows, the general features, namely the locations and strengths of centrifugal and turbulent pressure supported shocks, centrifugal barriers, and the funnel walls which are expected from analytical models, agree very well, though the details vary. Generally speaking, the locations of the shocks may be better obtained by SPH since the angular momentum is strictly preserved in SPH, but the shocks themselves are better resolved by TVD. The agreement of these code test results with analytical solutions provides us with confidence to apply these codes to more complex problems which we will discuss elsewhere.

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.

We study the evolution of viscous, isothermal, rotating, thin, axisymmetric accretion disks around a compact object using Smoothed Particle Hydrodynamics. We emphasize the effects of different choices of viscosity prescriptions on the evolution of angular momentum as well as other physical quantities of the disk. We show that a flow with the Shakura-Sunyaev viscosity prescription may produce only {\it shear} shock, where the angular momentum changes significantly across the shock. We present and study the effects of other prescriptions of viscosity which render the angular momentum continuous across the shock waves. In general, it is observed that for flows with a small viscosity, the shocks are weaker, form farther away and are wider as compared to the shocks in inviscid flows. If the viscosity is high, shocks do not form at all. The flow remains subsonic and Keplerian throughout the disk and becomes supersonic only very close to the horizon.

We present results of numerical simulation of inviscid thick accretion disks and wind flows around black holes. We use Smoothed Particle Hydrodynamics (SPH) technique for this purpose. Formation of thick disks are found to be preceded by shock waves travelling away from the centrifugal barrier. For a large range of the parameter space, the travelling shock settles at a distance close to the location obtained by a one-and-a-half dimensional model of inviscid accretion disks. Occasionally, it is observed that accretion processes are aided by the formation of oblique shock waves, particularly in the initial transient phase. The post-shock region (where infall velocity suddenly becomes very small) resembles that of the usual model of thick accretion disk discussed in the literature, though they have considerable turbulence. The flow subsequently becomes supersonic before falling into the black hole. In a large number of cases which we simulate, we find the formation of strong winds which are hot and subsonic when originated from the disk surface very close to the black hole but become supersonic within a few tens of the Schwarzschild radius of the blackhole. In the case of accretion of high angular momentum flow, very little amount of matter is accreted directly onto the black hole. Most of the matter is, however, first squeezed to a small volume close to the black hole, and subsequently expands and is expelled as a strong wind. It is quite possible that this expulsion of matter and the formation of cosmic radio jets is aided by the shock heating in the inner parts of the accretion disks.

We present results of numerical simulation of thin accretion disks and winds. We use the Smoothed Particle Hydrodynamics (SPH) technique for this purpose. We show that the simulation agrees very well with the recent theoretical work on the shock formation. The most significant conclusion is that shocks in an inviscid flow are extremely stable. For the first time, our work also removes the ambiguity in terms of the location and stability of shocks in adiabatic flows.