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.

We show that the solutions obtained in the paper `An exact solution for arbitrarily rotating gaseous polytropes with index unity' by Kong, Zhang, and Schubert represent only approximate solutions of the free-boundary Euler-Poisson system of equations describing uniformly rotating, self-gravitating polytropes with index unity. We discuss the quality of such solutions as approximations to the rigidly rotating equilibrium polytropic configurations.

We analyze the properties of a polytropic fluid which is radially accreted into a Schwarzschild black hole. The case where the adiabatic index gamma lies in the range 1 < gamma <= 5/3 has been treated in previous work. In this article we analyze the complementary range 5/3 < gamma <= 2. To this purpose, the problem is cast into an appropriate Hamiltonian dynamical system whose phase flow is analyzed. While for 1 < gamma <= 5/3 the solutions are always characterized by the presence of a unique critical saddle point, we show that when 5/3 < gamma <= 2, an additional critical point might appear which is a center point. For the parametrization used in this paper we prove that whenever this additional critical point appears, there is a homoclinic orbit.

We formulate new general-relativistic extensions of Newtonian rotation laws for self-gravitating stationary fluids. They have been used to re-derive, in the first post-Newtonian approximation, the well known geometric dragging of frames. We derive two other general-relativistic weak-field effects within rotating tori: the recently discovered dynamic anti-dragging and a new effect that measures the deviation from the Keplerian motion and/or the contribution of the fluids selfgravity. One can use the rotation laws to study the uniqueness and the convergence of the post-Newtonian approximations, and the existence of the post-Newtonian limits.

We investigate accreting disk systems with polytropic gas in Keplerian motion. Numerical data and partial analytic results show that the self-gravitation of the disk speeds up its rotation -- its rotational frequency is larger than that given by the well known strictly Keplerian formula that takes into account the central mass only. Thus determination of central mass in systems with massive disks requires great care -- the strictly Keplerian formula yields only an upper bound. The effect of self-gravity depends on geometric aspects of disk configurations. Disk systems with a small (circa $10^{-4}$) ratio of the innermost radius to the outermost disk radius have the central mass close to the upper limit, but if this ratio is of the order of unity then the central mass can be smaller by many orders of magnitude from this bound.

The Keplerian motion of accretion disks in Active Galactic Nuclei (AGN) is usually believed to be generated by a heavy central mass. We investigate accreting disk systems with polytropic gas in Keplerian rotation and obtain a phenomenological formula that relates the Keplerian angular frequency to the ratio of disk and central masses. Central mass approaches the Keplerian value, if the inner boundary of a disk is close to the minimal stable orbit of a black hole. These results are applied to NGC 4258, the unique AGN with a finely measured Keplerian rotation curve of the central disk, with the conclusion that its rotation curve is, in fact, determined by the central black hole. The mass of the accretion disk exceeds 100 solar masses.

All real solutions of the Lane-Emden equation for n = 5 are obtained in terms of Jacobian and Weierstrass elliptic functions. A new family of solutions is found. It is expressed by remarkably simple formulae involving Jacobian elliptic functions only. The general properties and discrete scaling symmetries of these new solutions are discussed. We also comment on their possible applications.

For stationary, barotropic fluids in Newtonian gravity we give simple criteria on the equation of state and the "law of motion" which guarantee finite or infinite extent of the fluid region (providing a priori estimates for the corresponding stationary Newton-Euler system). Under more restrictive conditions, we can also exclude the presence of "hollow" configurations. Our main result, which does not assume axial symmetry, uses the virial theorem as the key ingredient and generalises a known result in the static case. In the axially symmetric case stronger results are obtained and examples are discussed.

A continuum version of the virial theorem is derived for a radiating self-gravitating accretion disc around a compact object. The central object is point-like, but we can avoid the regularization of its gravitational potential. This is achieved by applying a modified Pohozaev-Rellich identity to the gravitational potential of the disk only. The theorem holds for general stationary configurations, including discontinuous flows (shock waves, contact discontinuities). It is used to test numerical solutions of a model of self-gravitating radiative accretion discs. The presented virial theorem should be useful in the analysis of those (possibly radiating) hydrodynamical systems in astrophysics where the central mass and the mass of the fluid are comparable and none of them can be neglected.

We studied a steadily accreting, geometrically thick disk model that selfconsistently takes into account selfgravitation of the polytropic gas, its interaction with the radiation and the mass accretion rate. The accreting mass is injected inward in the vicinity of the central $z=0$ plane, where also radiation is assumed to be created. The rest of the disk remains approximately stationary. Only conservation laws are employed and the gas-radiation interaction in the bulk of the disk is described in the thin-gas approximation. We demonstrate that this scheme is numerically viable and yields a structure of the bulk that is influenced by the radiation and (indirectly) by the prescribed mass accretion rate. The obtained disk configurations are typical for environments in Active Galactic Nuclei (AGN), with the central mass of the order of $10^7 M_{\astrosun}$ to $10^8 M_{\astrosun}$, quasi-Keplerian rotation curves, disk masses ranging from about $10^6 M_{\astrosun}$ to $10^7 M_{\astrosun}$, and the luminosity ranging from $10^6 L_{\astrosun}$ to $10^9 L_{\astrosun}$. These luminosities are much lower than the corresponding Eddington limit.

Corrugation instabilities occurring for solutions of the Riemann problem in relativistic hydrodynamics in which the fluid moves with a non-zero velocity tangent to the initial discontinuity are studied numerically. We perform simulations both for ultrarelativistic and perfect gas equations of state. We focus on a set of problems with moderately relativistic velocities but exhausting all possible wave patterns of solutions. Perturbations are applied to the shape of the initial discontinuity. Instabilities that develop are only restricted to a region around a contact discontinuity. Both shock and rarefaction waves appear to be stable.

We estimate the minimal mass of selfgravitating polytropic disks using the famous Sobolev inequality. This bound resembles the well known mass formula for Lane-Emden stars. For ideal gas with the polytropic index n = 3 the minimal mass is not smaller than the Jeans mass. The accuracy of the estimate is verified in a number of numerical examples. The bound works well for heavy selfgravitating disks and is less useful for light disks.

An analysis of the propagation of non stationary waves in the adiabatic region of stellar interior is presented. An equation of motion with an effective potential is derived, similar to the Zerilli equation known in the propagation of gravitational waves. The Huyghens principle is violated in this case and the energy diffusion outward null cones is expected. Numerical calculations demonstrate that the diffusion is weak for the case of standard Solar model; thus no significant effect corresponding to quasinormal modes can be expected. The likely reason for the absence of stronger features is the restriction of our analysis to adiabatic regions only, where the breakdown of the Huyghens principle is insignificant.