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.

In this article we develop a Primitive Variable Recovery Scheme (PVRS) to solve any system of coupled differential conservative equations. This method obtains directly the primitive variables applying the chain rule to the time term of the conservative equations. With this, a traditional finite volume method for the flux is applied in order avoid violation of both, the entropy and "Rankine-Hugoniot" jump conditions. The time evolution is then computed using a forward finite difference scheme. This numerical technique evades the recovery of the primitive vector by solving an algebraic system of equations as it is often used and so, it generalises standard techniques to solve these kind of coupled systems. The article is presented bearing in mind special relativistic hydrodynamic numerical schemes with an added pedagogical view in the appendix section in order to easily comprehend the PVRS. We present the convergence of the method for standard shock-tube problems of special relativistic hydrodynamics and a graphical visualisation of the errors using the fluctuations of the numerical values with respect to exact analytic solutions. The PVRS circumvents the sometimes arduous computation that arises from standard numerical methods techniques, which obtain the desired primitive vector solution through an algebraic polynomial of the charges.