We integrate numerically axially symmetric stationary Einstein equations describing self-gravitating disks around spinless black holes. The numerical scheme is based on a method developed by Shibata, but contains important new ingredients. We derive a new general-relativistic Keplerian rotation law for self-gravitating disks around spinning black holes. Former results concerning rotation around spin-less black holes emerge in the limit of a vanishing spin parameter. These rotation curves might be used for the description of rotating stars, after appropriate modification around the symmetry axis. They can be applied to the description of compact torus--black hole configurations, including active galactic nuclei or products of coalescences of two neutron stars.

In this paper we investigate a steady accretion within the Einstein-Straus vacuole, in the presence of the cosmological constant. The dark energy damps the mass accretion rate and --- above certain limit --- completely stops the steady accretion onto black holes, which in particular is prohibited in the inflation era and after (roughly) $10^{12}$ years from Big Bang (assuming the presently known value of the cosmological constant). Steady accretion would not exist in the late phases of the Penrose's scenario - known as the Weyl curvature hypothesis - of the evolution of the Universe.

We analyze stationary accretion of selfgravitating gas onto a compact center within general-relativistic radiation hydrodynamics. Spherical symmetry and thin gas approximation are assumed. Numerical investigation shows that transonic flows exist for small redshifts and they cease to exist for high redshifts and high luminosities. There exist two branches of flows (subsonic or supersonic) that originate at a bifurcation point and that embrace the set of subsonic solutions. The morphology of the set of subsonic solutions is essentially independent of redshifts and flows that belong to their boundary provide estimates of the gas abundance of subsonic solutions. It appears that prescribed boundary data guarantee uniqueness only of the bifurcation point, and that the latter has maximal luminosity.

Quasi-stationary flows of gas accreting onto a compact center are analyzed in the framework of general-relativistic radiation hydrodynamics, under assumptions of spherical symmetry and thin gas approximation. Numerical investigation shows that luminosity, redshift and gas abundance are correlated. The gas can constitute up to one third of the total mass of brightest low-redshift sources, but its abundance goes down to 1/30 for sources with luminosities close to the Eddington limit.

We investigate newtonian description of accreting compact bodies with hard surfaces, including luminosity and selfgravitation of polytropic perfect fluids. This nonlinear integro-differential problem reduces, under appropriate boundary conditions, to an algebraic relation between luminosity and the gas abundance in stationary spherically symmetric flows. There exist, for a given luminosity, asymptotic mass and the asymptotic temperature, two sub-critical solutions that bifurcate from a critical point. They differ by the fluid content and the mass of the compact centre.

Corrected misprints and omissions.