We present a novel, relativistic accretion model for accretion onto a Schwarzschild black hole. This consists of a purely hydrodynamical mechanism in which, by breaking spherical symmetry, a radially accreting flow transitions into an inflow-outflow configuration. The spherical symmetry is broken by considering that the accreted material is more concentrated on an equatorial belt, leaving the polar regions relatively under-dense. What we have found is a flux-limited accretion regime in which, for a sufficiently large accretion rate, the incoming material chokes at a gravitational bottleneck and the excess flux is redirected by the density gradient as a bipolar outflow. The threshold value at which the accreting material chokes is of the order of the mass accretion rate found in the spherically symmetric case studied by Bondi and Michel. We describe the choked accretion mechanism first in terms of a general relativistic, analytic toy model based on the assumption of an ultrarelativistic stiff fluid. We then relax this approximation and, by means of numerical simulations, show that this mechanism can operate also for general polytropic fluids. Interestingly, the qualitative inflow-outflow morphology obtained appears as a generic result of the proposed symmetry break, across analytic and numeric results covering both the Newtonian and relativistic regimes. The qualitative change in the resulting steady state flow configuration appears even for a very small equatorial to polar density contrast ($\sim 0.1\,\%$) in the accretion profile. Finally, we discuss the applicability of this model as a jet-launching mechanism in different astrophysical settings.

Steady state, spherically symmetric accretion flows are well understood in terms of the Bondi solution. Spherical symmetry however, is necessarily an idealized approximation to reality. Here we explore the consequences of deviations away from spherical symmetry, first through a simple analytic model to motivate the physical processes involved, and then through hydrodynamical, numerical simulations of an ideal fluid accreting onto a Newtonian gravitating object. Specifically, we consider axisymmetric, large-scale, small amplitude deviations in the density field such that the equatorial plane is over dense as compared to the polar regions. We find that the resulting polar density gradient dramatically alters the Bondi result and gives rise to steady state solutions presenting bipolar outflows. As the density contrast increases, more and more material is ejected from the system, attaining speeds larger than the local escape velocities for even modest density contrasts. Interestingly, interior to the outflow region, the flow tends locally towards the Bondi solution, with a resulting total mass accretion rate through the inner boundary $choking$ at a value very close to the corresponding Bondi one. Thus, the numerical experiments performed suggest the appearance of a maximum achievable accretion rate, with any extra material being ejected, even for very small departures from spherical symmetry.

We present a novel analytic model of relativistic wind accretion on to a Schwarzschild black hole. This model constitutes a general relativistic extension of the classical model of wind accretion by Bondi, Hoyle, and Lyttleton (BHL). As in BHL, this model is based on the assumptions of steady state, axisymmetry, and ballistic motion. Analytic expressions are provided for the wind streamlines while simple numerical schemes are presented for calculating the corresponding accretion rate and density field. The resulting accretion rate is greater in the relativistic model when compared to the Newtonian BHL one. Indeed, it is two times greater for asymptotic wind speeds $v_\infty \ge 0.4\,\mathrm{c}$ and more than an order of magnitude greater for $v_\infty \ge 0.8\,\mathrm{c}$. We have compared this full relativistic model versus numerical simulations performed with the free GNU Public Licensed hydrodynamics code AZTEKAS and found a good agreement for the streamlines in the upstream region of the flow and also, to within 10 per cent, for the corresponding accretion rates.

We propose an approximate approach for studying the relativistic regime of stellar tidal disruptions by rotating massive black holes. It combines an exact relativistic description of the hydrodynamical evolution of a test fluid in a fixed curved spacetime with a Newtonian treatment of the fluid's self-gravity. Explicit expressions for the equations of motion are derived for Kerr spacetime using two different coordinate systems. We implement the new methodology within an existing Newtonian Smoothed Particle Hydrodynamics code and show that including the additional physics involves very little extra computational cost. We carefully explore the validity of the novel approach by first testing its ability to recover geodesic motion, and then by comparing the outcome of tidal disruption simulations against previous relativistic studies. We further compare simulations in Boyer--Lindquist and Kerr--Schild coordinates and conclude that our approach allows accurate simulation even of tidal disruption events where the star penetrates deeply inside the tidal radius of a rotating black hole. Finally, we use the new method to study the effect of the black hole spin on the morphology and fallback rate of the debris streams resulting from tidal disruptions, finding that while the spin has little effect on the fallback rate, it does imprint heavily on the stream morphology, and can even be a determining factor in the survival or disruption of the star itself. Our methodology is discussed in detail as a reference for future astrophysical applications.

We present a simple, analytic model for the accretion flow of an incompressible wind onto a gravitating object. This solution corresponds to the Newtonian limit of a previously known relativistic model for a fluid obeying a stiff equation of state for which the sound speed is constant everywhere and equal to the speed of light. The new solution should be useful as a benchmark test for numerical hydrodynamics codes and, moreover, it can be used as an illustrative example in a gas dynamics course.

We perform hydro- and magnetohydrodynamical general relativistic simulations of a tidal disruption of a $0.1\,M_\odot$ red dwarf approaching a $10^5\,M_\odot$ non-rotating massive black hole on a close (impact parameter $\beta=10$) elliptical (eccentricity $e=0.97$) orbit. We track the debris self-interaction, circularization, and the accompanying accretion through the black hole horizon. We find that the relativistic precession leads to the formation of a self-crossing shock. The dissipated kinetic energy heats up the incoming debris and efficiently generates a quasi-spherical outflow. The self-interaction is modulated because of the feedback exerted by the flow on itself. The debris quickly forms a thick, almost marginally bound disc that remains turbulent for many orbital periods. Initially, the accretion through the black hole horizon results from the self-interaction, while in the later stages it is dominated by the debris originally ejected in the shocked region, as it gradually falls back towards the hole. The effective viscosity in the debris disc stems from the original hydrodynamical turbulence, which dominates over the magnetic component. The radiative efficiency is very low because of low energetics of the gas crossing the horizon and large optical depth that results in photon trapping.

Solitary stars that wander too close to their galactic centres can become tidally disrupted, if the tidal forces due to the supermassive black hole (SMBH) residing there overcome the self-gravity of the star. If the star is only partially disrupted, so that a fraction survives as a self-bound object, this remaining core will experience a net gain in specific orbital energy, which translates into a velocity "kick" of up to $\sim 10^3$ km/s. In this paper, we present the result of smoothed particle hydrodynamics (SPH) simulations of such partial disruptions, and analyse the velocity kick imparted on the surviving core. We compare $\gamma$ = 5/3 and $\gamma$ = 4/3 polytropes disrupted in both a Newtonian potential, and a generalized potential that reproduces most relativistic effects around a Schwarzschild black hole either exactly or to excellent precision. For the Newtonian case, we confirm the results of previous studies that the kick velocity of the surviving core is virtually independent of the ratio of the black hole to stellar mass, and is a function of the impact parameter $\beta$ alone, reaching at most the escape velocity of the original star. For a given $\beta$, relativistic effects become increasingly important for larger black hole masses. In particular, we find that the kick velocity increases with the black hole mass, making larger kicks more common than in the Newtonian case, as low-$\beta$ encounters are statistically more likely than high-$\beta$ encounters. The analysis of the tidal tensor for the generalized potential shows that our results are robust lower limits on the true relativistic kick velocities, and are generally in very good agreement with the exact results.

We present a generalized Newtonian description of particle dynamics valid for any spherically symmetric, static black hole spacetime. This approach is derived from the geodesic motion of test particles in the low-energy limit. It reproduces exactly the location of the marginally stable, marginally bound, and photon circular orbits; the radial dependence of the energy and angular momentum of circular orbits; parabolic motion; pericentre shift; and the spatial projection of general trajectories. As explicit examples of the new prescription, we apply it to the Schwarzschild, Schwarzschild-de Sitter, Reissner-Nordstr\"om, Ay\'on Beato-Garc\'ia, and Kehagias-Sfetsos spacetimes. In all of these examples, the orbital and epicyclic frequencies are reproduced to better than 10%. The resulting equations of motion can be implemented easily and efficiently within existing Newtonian frameworks.

A generalized Newtonian potential is derived from the geodesic motion of test particles in Schwarzschild spacetime. This potential reproduces several relativistic features with higher accuracy than commonly used pseudo-Newtonian approaches. The new potential reproduces the exact location of the marginally stable, marginally bound, and photon circular orbits, as well as the exact radial dependence of the binding energy and the angular momentum of these orbits. Moreover, it reproduces the orbital and epicyclic angular frequencies to better than 6%. In addition, the spatial projections of general trajectories coincide with their relativistic counterparts, while the time evolution of parabolic-like trajectories and the pericentre advance of elliptical-like trajectories are both reproduced exactly. We apply this approach to a standard thin accretion disc and find that the efficiency of energy extraction agrees to within 3% with the exact relativistic value, while the energy flux per unit area as a function of radius is reproduced everywhere to better than 7%. As a further astrophysical application we implement the new approach within a smoothed particle hydrodynamics code and study the tidal disruption of a main-sequence star by a supermassive black hole. The results obtained are in very good agreement with previous relativistic simulations of tidal disruptions in Schwarzschild spacetime. The equations of motion derived from this potential can be implemented easily within existing Newtonian hydrodynamics codes with hardly any additional computational effort.

We present a relativistic model for the stationary axisymmetric accretion flow of a rotating cloud of non-interacting particles falling onto a Kerr black hole. Based on a ballistic approximation, streamlines are described analytically in terms of timelike geodesics, while a simple numerical scheme is introduced for calculating the density field. A novel approach is presented for describing all of the possible types of orbit by means of a single analytic expression. This model is a useful tool for highlighting purely relativistic signatures in the accretion flow dynamics coming from a strong gravitational field with frame-dragging. In particular, we explore the coupling due to this between the spin of the black hole and the angular momentum of the infalling matter. Moreover, we demonstrate how this analytic solution may be used for benchmarking general relativistic numerical hydrodynamics codes by comparing it against results of smoothed particle hydrodynamics simulations for a collapsar-like setup. These simulations are performed first for a ballistic flow (with zero pressure) and then for a hydrodynamical one where we measure the effects of pressure gradients on the infall, thus exploring the extent of applicability of the ballistic approximation.

We construct a general relativistic model for the accretion flow of a rotating finite cloud of non-interacting particles infalling onto a Schwarzschild black hole. The streamlines start at a spherical shell, where boundary conditions are fixed, and are followed down to the point at which they either cross the black hole horizon or become incorporated into an equatorial thin disc. Analytic expressions for the streamlines and the velocity field are given, in terms of Jacobi elliptic functions, under the assumptions of stationarity and ballistic motion. A novel approach allows us to describe all of the possible types of orbit with a single formula. A simple numerical scheme is presented for calculating the density field. This model is the relativistic generalisation of the Newtonian one developed by Mendoza, Tejeda, Nagel, 2009 and, due to its analytic nature, it can be useful in providing a benchmark for general relativistic hydrodynamical codes and for exploring the parameter space in applications involving accretion onto black holes when the approximations of steady state and ballistic motion are reasonable ones.

We construct a steady analytic accretion flow model for a finite rotating gas cloud that accretes material to a central gravitational object. The pressure gradients of the flow are considered to be negligible and so, the flow is ballistic. We also assume a steady flow and consider the particles at the boundary of the spherical cloud to be rotating as a rigid body, with a fixed amount of inwards radial velocity. This represents a generalisation to the traditional infinite gas cloud model described by Ulrich (1976). We show that the streamlines and density profiles obtained deviate largely from the ones calculated by Ulrich. The extra freedom in the choice of the parameters on the model can naturally account for the study of protostars formed in dense clusters by triggered mechanisms, where a wide variety of external physical mechanisms determine the boundary conditions. Also, as expected, the model predicts the formation of an equatorial accretion disc about the central object with a radius different from the one calculated by Ulrich (1976).