Hydrodynamical non-radiative accretion flows in two dimensions

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows that cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function ν . We study the flows generated as the amplitude and form of ν are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by and a large fraction of the mass participating in the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long-time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. The two-dimensional structure of the time-averaged flow is marginally stable to the Høiland criterion, indicating that convection is efficient. Near the equatorial plane, the radial profiles of the time-averaged variables are power laws with an index that depends on the radial scaling of the shear stress. A stress in which ν ∝ r ^1/2 recovers the widely studied self-similar solution corresponding to an ' α -disc'. We find that, regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.     

On the absence of winds in advection-dominated accretion flows

We show that recently published assertions that advection-dominated accretion flows (ADAFs) require the presence of strong winds are unfounded because they assume that low radiative efficiency in flows accreting at low rates on to black holes implies vanishing radial energy and angular momentum fluxes through the flow (which is also formulated in terms of the 'Bernoulli function' being positive). This, however, is a property only of self-similar solutions which are an inadequate representation of global accretion flows. We recall general properties of accretion flows on to black holes and show that such, necessarily transonic, flows may have either d positive or negative Bernoulli function depending on the flow viscosity. Flows with low viscosities ( α ≲0.1 in the α -viscosity model) have a negative Bernoulli function. Without exception, all 2D and 1D numerical models of low-viscosity flows constructed to date experience no significant outflows. At high viscosities the presence of outflows depends on the assumed viscosity, on the equation of state and on the outer boundary condition. The positive sign of the Bernoulli function invoked in this context is irrelevant to the presence of outflows. As an illustration, we recall 2D numerical models with moderate viscosity that have positive values of the Bernoulli function and experience no outflows. ADAFs, therefore, do not differ from this point of view from thin Keplerian discs: they may have, but they do not have to have, strong winds.     

Behaviour of dissipative accretion flows around black holes

We investigate the behaviour of dissipative accreting matter close to a black hole, as this provides important observational features of galactic and extragalactic black hole candidates. We find a complete set of global solutions in the presence of viscosity and synchrotron cooling. We show that advective accretion flow can have a standing shock wave and the dynamics of the shock is controlled by the dissipation parameters (both viscosity and cooling). We study the effective region of the parameter space for standing as well as oscillating shock. We find that the shock front always moves towards the black hole as the dissipation parameters are increased. However, viscosity and cooling have opposite effects in deciding the solution topologies. We obtain two critical cooling parameters that separate the nature of the accretion solution.     

Self-similar solution of hot accretion flows with ordered magnetic field and outflow

Observations and numerical magnetohydrodynamic (MHD) simulations indicate the existence of outflows and ordered large-scale magnetic fields in the inner region of hot accretion flows. In this paper, we present the self-similar solutions for advection-dominated accretion flows (ADAFs) with outflows and ordered magnetic fields. Stimulated by numerical simulations, we assume that the magnetic field has a strong toroidal component and a vertical component in addition to a stochastic component. We obtain the self-similar solutions to the equations describing the magnetized ADAFs, taking into account the dynamical effects of the outflow. We compare the results with the canonical ADAFs and find that the dynamical properties of ADAFs such as radial velocity, angular velocity and temperature can be significantly changed in the presence of ordered magnetic fields and outflows. The stronger the magnetic field is, the lower the temperature of the accretion flow will be and the faster the flow rotates. The relevance to observations is briefly discussed.     

Unipolar outflows and global meridional circulations in rotating accretion flows

Using two-dimensional simulations of non-radiative viscous rotating black hole accretion flows, we show that the flows with α ∼0.1–0.3 self-organize to form stationary unipolar or bipolar outflows accompanied by global meridional circulations. The required energy comes, with efficiency ∼0.001–0.01, from the matter directly accreted on to the black hole. Observational implications are discussed.     

Iron Kα line profiles and the inner boundary condition of accretion flows

Recent X-ray observations have shown evidence for exceptionally broad and skewed iron Kα emission lines from several accreting black hole systems. The lines are assumed to be due to fluorescence of the accretion disc illuminated by a surrounding corona and require a steep emissivity profile increasing into the innermost radius. This appears to question both standard accretion disc theory and the zero-torque assumption for the inner boundary condition, both of which predict a much less extreme profile. Instead it argues that a torque may be present due to magnetic coupling with matter in the plunging region or even to the spinning black hole itself. Discussion so far has centred on the torque acting on the disc. However, the crucial determinant of the iron line profile is the radial variation of the power radiated in the corona. Here we study the effects of different inner boundary conditions on the coronal emissivity and on the profiles of the observable Fe Kα lines. We argue that in the extreme case where a prominent highly redshifted component of the iron line is detected, requiring a steep emissivity profile in the innermost part and a flatter one outside, energy from the gas plunging into the black hole is being fed directly to the corona.     

Non-axisymmetric instabilities in shocked adiabatic accretion flows

We investigate the linear stability of a shocked accretion flow on to a black hole in the adiabatic limit. Our linear analyses and numerical calculations show that, despite the post-shock deceleration, the shock is generally unstable to non-axisymmetric perturbations. The simulation results of Molteni, Tóth & Kuznetsov can be well explained by our linear eigenmodes. The mechanism of this instability is confirmed to be based on the cycle of acoustic waves between the corotation radius and the shock. We obtain an analytical formula to calculate the oscillation period from the physical parameters of the flow. We argue that the quasi-periodic oscillation should be a common phenomenon in accretion flows with angular momentum.     

Fractal features in accretion discs

Fractal concepts have been introduced in the accretion disc as a new feature. Due to the fractal nature of the flow, its continuity condition undergoes modifications. The conserved stationary fractal flow admits only saddle points and centre-type points in its phase portrait. Completely analytical solutions of the equilibrium point conditions indicate that the fractal properties enable the flow to behave like an effective continuum of lesser density, and facilitate the generation of transonicity. However, strongly fractal flows inhibit multitransonicity from developing. The mass accretion rate exhibits a fractal scaling behaviour, and the entire fractal accretion disc is stable under linearized dynamic perturbations.     

Magnetized accretion flows: effects of gas pressure

We study how axisymmetric magnetohydrodynamic (MHD) accretion flows depend on γ adiabatic index in the polytropic equation of state. This work is an extension of Mościbrodzka & Proga, where we investigated the γ dependence of two-dimensional Bondi-like accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if simulations for various γ can give us insights into the problem of various modes of accretion observed in several types of accretion systems, such as black hole binaries (BHBs), active galactic nuclei (AGN) and gamma-ray bursts. We find that for  γ≳ 4/3  , the fast-rotating flow forms a thick torus that is supported by rotation and gas pressure. As shown before for  γ= 5/3  , such a torus produces a strong, persistent bipolar outflow that can significantly reduce the polar funnel accretion of a slowly rotating flow. For low γ, close to 1, the torus is thin and is supported by rotation. The thin torus produces an unsteady outflow which is too weak to propagate throughout the polar funnel inflow. Compared to their HD counterparts, the MHD simulations show that the magnetized torus can produce an outflow and does not exhibit regular oscillations. Generally, our simulations demonstrate how the torus thickness affects the outflow production. They also support the notion that the geometrical thickness of the torus correlates with the power of the torus outflow. Our results, applied to observations, suggest that the torus ability to radiatively cool and become thin can correspond to a suppression of a jet as observed in the BHBs during a transition from a hard/low to soft/high spectral state and a transition from a quiescent to hard/low state in AGN.     

Energy-dependent variability from accretion flows

We develop a formalism to calculate energy-dependent fractional variability (rms) in accretion flows. We consider rms spectra resulting from radial dependences of the level of local variability (as expected from the propagation of disturbances in accretion flows) assuming the constant shape of the spectrum emitted at a given radius. We consider the cases when the variability of the flow is either coherent or incoherent between different radial zones. As an example of local emission, we consider blackbody, Wien and thermal Comptonization spectra. In addition to numerical results, we present a number of analytical formulae for the resulting rms. We also find an analytical formula for the disc Wien spectrum, which we find to be a very good approximation to the disc blackbody. We compare our results to the rms spectrum observed in an ultrasoft state of GRS 1915+105.     

Model dependence of transonic properties of accretion flows around black holes

We analytically study how the behaviour of accretion flows changes when the flow model is varied. We study the transonic properties of the conical flow, a flow of constant height and a flow in vertical equilibrium, and show that all these models are basically identical, provided that the polytropic constant is suitably changed from one model to another. We show that this behaviour is extendible even when standing shocks are produced in the flow. The parameter space where shocks are produced remains roughly identical in all these models when the same transformation among the polytropic indices is used. We present applications of these findings.