Complete integrability of geodesic motion, the well known feature of the fields of isolated stationary black holes, can easily be "spoilt" by the presence of some additional source (even if highly symmetric). In previous papers, we used various methods to show how free time-like motion becomes chaotic if the gravitational field of the Schwarzschild black hole is perturbed by that of a circular disc or ring, considering specifically the inverted first disc of the Morgan-Morgan counter-rotating family and the Bach-Weyl ring as the additional sources. The present paper focuses on two new points. First, since the Bach-Weyl thin ring is physically quite unsatisfactory, we now repeat some of the analysis for a different, Majumdar-Papapetrou--type (extremally charged) ring around an extreme Reissner-Nordstr\"om black hole, and compare the results with those obtained before. We also argue that such a system is in fact more relevant astrophysically than it may seem. Second, we check numerically, for the latter system as well as for the Schwarzschild black hole encircled by the inverted Morgan-Morgan disc, how indicative is the geometric (curvature) criterion for chaos suggested by Sota, Suzuki & Maeda (1996). We also add a review of the literature where the relevance of geometric criteria in general relativity (as well as elsewhere) has been discussed for decades.

We perform 1D/2D/3D relativistic hydrodynamical simulations of accretion flows with low angular momentum, filling the gap between spherically symmetric Bondi accretion and disc-like accretion flows. Scenarios with different directional distributions of angular momentum of falling matter and varying values of key parameters such as spin of central black hole, energy and angular momentum of matter are considered. In some of the scenarios the shock front is formed. We identify ranges of parameters for which the shock after formation moves towards or outwards the central black hole or the long lasting oscillating shock is observed. The frequencies of oscillations of shock positions which can cause flaring in mass accretion rate are extracted. The results are scalable with mass of central black hole and can be compared to the quasi-periodic oscillations of selected microquasars (such as GRS 1915+105, XTE J1550-564 or IGR J17091-3624), as well as to the supermassive black holes in the centres of weakly active galaxies, such as Sgr $A^{*}$.

We study the X-ray emission of the microquasar XTE J1550-564 and analyze the properties of its light curves using the recurrence analysis method. The indicators for non-linear dynamics of the accretion flow are found in the very high state and soft state of this source. The significance of deterministic variability depends on the energy band. We discuss the non-linear dynamics of the accretion flow in the context of the disc-corona geometry and propagating oscillations in the accretion flow.

The high energy radiation emitted by black hole X-ray binaries originates in an accretion disk, hence the variability of the lightcurves mirrors the dynamics of the disc. We study the time evolution of the emitted flux in order to find evidences, that low dimensional non-linear equations govern the accretion flow. Here we test the capabilities of our novel method to find chaotic behaviour on the two numerical time series describing the motion of a test particle around a black hole surrounded by a thin massive disc, one being regular and the other one chaotic.

We study the low angular momentum accretion of matter onto Schwarzschild or Kerr black hole using fully relativistic numerical simulations. Our aim is to include the rotation of both the gas and the black hole and the presence of the magnetic field into the picture. We compare resulting structure of transonic flow with results of 1D pseudo-Newtonian semi-analytical computations of non-magnetized flow.

The black hole candidates exhibit fast variability of their X-ray emission on a wide range of timescales. The short, coherent variations, with frequencies above 1 Hz, are referred to as quasi-periodic oscillations, and are generally explained by resonant effects in the black hole accretion flow. The purely stochastic variability that occurs due to turbulent conditions in the plasma, is quantified by the power density spectra and appears practically in all types of sources and their spectral states. The specific kind of quasi-periodic flares is expected, when the global structure of the accretion flow, governed by the nonlinear hydrodynamics, induces fluctuations around a fixed point solution. These conditions, which occur at high accretion rates, lead to the variability in the sense of deterministic chaos. We study the nonlinear behaviour of X-ray sources using the recurrence analysis method. We estimate quantitatively the indications for deterministic chaos, such as the Renyi's entropy, for the observed time series, and we compare them with the surrogate data. Using the data collected by RXTE satellite, we reveal the oscillations pattern and the observable properties of six black hole systems. We confirm that both the well known microquasar GRS 1915+105, as well as its recently discovered analogue, IGR J17091-3624, exhibit variability characteristic to deterministic chaotic system. We conclude therefore that the underlying nature of the process must be intrinsically connected in these sources with the accretion flow instability, that leads to the limit cycle oscillations around a fixed point. Furthermore, we studied four other black hole candidates with high accretion rates. For GX 339-4, GRO J1655-40 and XTE J1550-564 we obtained significant results for some observations, whereas for XTE J1650-500 no such data set was found.

"Almost all" seems to be known about isolated stationary black holes in asymptotically flat space-times and about the behaviour of {\em test} matter and fields in their backgrounds. The black holes likely present in galactic nuclei and in some X-ray binaries are commonly being represented by the Kerr metric, but actually they are not isolated (they are detected only thanks to a strong interaction with the surroundings), they are not stationary (black-hole sources are rather strongly variable) and they also probably do not live in an asymptotically flat universe. Such "perturbations" may query the classical black-hole theorems (how robust are the latter against them?) and certainly affect particles and fields around, which can have observational consequences. In the present contribution we examine how the geodesic structure of the static and axially symmetric black-hole space-time responds to the presence of an additional matter in the form of a thin disc or ring. We use several different methods to show that geodesic motion may become chaotic, to reveal the strength and type of this irregularity and its dependence on parameters. The relevance of such an analysis for galactic nuclei is briefly commented on.

We derive the conditions for shock formation in a quasi-spherical, slightly rotating flows. We verify the results of semi-analytical, stationary calculations with the time evolution studied by numerical hydro-simulations, and we study the oscillations of the shock position. We also study the behaviour of flows with varying specific angular momentum, where the 'hysteresis' type of loop is found when passing through the multiple sonic points region. Our results are in agreement with the timescales and shapes of the luminosity flares observed in Sgr A*. These models may also be applicable for the Galactic stellar mass black holes, like GX 339-4 or GRS 1915+105, where periodic oscillations of X-ray luminosity are detected.

We continue the study of time-like geodesic dynamics in exact static, axially and reflection symmetric space-times describing the fields of a Schwarzschild black hole surrounded by thin discs or rings. In the first paper of this series, the rise (and decline) of geodesic chaos with ring/disc mass and position and with test particle energy was revealed on Poincar\'e sections and on time series of position or velocity and their power spectra. In the second paper we compared these results with those obtained by two recurrence methods, focusing on "sticky" orbits whose different parts show different degrees of chaoticity. Here we complement the analysis by using several Lyapunov-type coefficients which quantify the rate of orbital divergence. After comparing the results with those obtained by the previous methods, we specifically consider a system involving a black hole surrounded by a small thin disc or a large ring, having in mind the configuration which probably occurs in galactic nuclei. Within the range of parameters which roughly corresponds to our Galactic center, we found that the black-hole accretion disc does not have a significant gravitational effect on the dynamics of free motion at larger radii, while the inner circumnuclear molecular ring (concentrated above 1 parsec radius) can only induce some irregularity in motion of stars ("particles") on smaller radii if its mass reaches 10 to 30% of the central black hole (which is the upper estimate given in the literature), if it is sufficiently compact (which does not hold but maybe for its inner rim) and if the stars can get to its close vicinity. The outer dust ring between 60 and 100 parsecs appears to be less important for the geodesic dynamics in its interior.