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.

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 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.