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.