This article reports on three-dimensional (3-D) MHD simulations of non-rotating and rapidly rotating black holes and the adjacent black hole accretion disk magnetospheres. A particular emphasis is placed on the vertical magnetic flux that is advected inward from large radii and threads the equatorial plane near the event horizon. In both cases of non-rotating and rotating black holes, the existence of a significant vertical magnetic field in this region is like a switch that creates powerful jets. There are many similarities in the vertical flux dynamics in these two cases in spite of the tremendous enhancement of azimuthal twisting of the field lines and enhancement of the jet power because of an "ergospheric disk" in the Kerr metric. A 3-D approach is essential because two-dimensional axisymmetric flows are incapable of revealing the nature of vertical flux near a black hole. Poloidal field lines from the ergospheric accretion region have been visualized in 3-D and much of the article is devoted to a formal classification of the different manifestations of vertical flux in the Kerr case.

A mechanism for the transition between low/hard, high/soft, and steep power law (SPL) spectral states in black hole X-ray binaries is proposed. The low/hard state is explained by the development of a magnetically arrested accretion disk attributable to the accumulation of a vertical magnetic field in a central bundle. This disk forms powerful jets and consists of thin spiral accretion streams of a dense optically thick plasma surrounded by hot, magnetized, optically thin corona, which emits most of the energy in hard X-rays. State transition occurs because of the quasi-periodic or random inversion of poloidal magnetic fields in the accretion flow supplied by the secondary star. The inward advection of the inverted field results in a temporal disappearance of the central bundle caused by the annihilation of the opposed fields and restoration of the optically thick disk in the innermost region. This disk represents the high/soft state. The SPL state develops at the period of intensive field annihilation and precedes the high/soft state. The continuous supply of the inverted field leads to a new low/hard state because of the formation of another magnetically arrested disk.

The dynamics and structure of accretion disks, which accumulate the vertical magnetic field in the centers, are investigated using two- and three-dimensional MHD simulations. The central field can be built up to the equipartition level and disrupts a nearly axisymmetric outer accretion disk inside a magnetospheric radius, forming a magnetically arrested disk (MAD). In the MAD, the mass accretes in a form of irregular dense spiral streams and the vertical field, split into separate bundles, penetrates through the disk plane in low-density magnetic islands. The accreting mass, when spiraling inward, drags the field and twists it around the axis of rotation, resulting in collimated Poynting jets in the polar directions. These jets are powered by the accretion flow with the efficiency up to ~1.5% (in units \dot{M}c^2). The spiral flow pattern in the MAD is dominated by modes with low azimuthal wavenumbers m~1-5 and can be a source of quasi-periodic oscillations in the outgoing radiation. The formation of MAD and Poynting jets can naturally explain the observed changes of spectral states in Galactic black hole binaries. Our study is focused on black hole accretion flows; however, the results can also be applicable to accretion disks around nonrelativistic objects, such as young stellar objects and stars in binary systems.

Spherical (nonrotating) accretion flows with small-scale magnetic fields have been investigated using three-dimensional, time-dependent MHD simulations. These simulations have been designed to model high-resolution (quasi) steady accretion flows in a wedge computational domain that represents a small fraction of the full spherical domain. Subsonic and supersonic (super-fast-magnetosonic) accretion flows have been considered. Two accretion regimes have been studied: conservative, or radiatively inefficient, and nonconservative, in which the heat released in magnetic reconnections is completely lost. The flows in both regimes are turbulent. They show the flattened radial density profiles and reduction of the accretion velocities and mass accretion rates in comparison with hydrodynamic Bondi flows. In the conservative regime, the turbulence is more intensive and supported mostly by thermal convection. In the nonconservative regime, the turbulence is less intensive and supported by magnetic buoyancy and various magnetic interactions. We have concluded that steady, supersonic spherical accretion cannot be developed in the presence of small-scale magnetic fields.

We calculate synchrotron radiation in three-dimensional pseudo-Newtonian magnetohydrodynamic simulations of radiatively inefficient accretion flows. We show that the emission is highly variable at optically thin frequencies, with order of magnitude variability on time-scales as short as the orbital period near the last stable orbit; this emission is linearly polarized at the 20-50 % level due to the coherent toroidal magnetic field in the flow. At optically thick frequencies, both the variability amplitude and polarization fraction decrease significantly with decreasing photon frequency. We argue that these results are broadly consistent with the observed properties of Sgr A* at the Galactic Center, including the rapid infrared flaring.

We consider the effects of advection and radial gradients of pressure and radial drift velocity on the structure of accretion disks around black holes with proper description of optically thick/thin transitions. We concentrated our efforts on the models with large accretion rate. Contrary to disk models neglecting advection, we find that continuous solutions extending from the outer disk regions to the inner edge exist for all accretion rates we have considered. We show that the sonic point moves outward with increasing accretion rate, and that in the innermost disk region advection acts as a heating process that may even dominate over dissipative heating. Despite the importance of advection on it's structure, the disk remains geometrically thin. Global solutions of advective accretion disks, which describe continuously the transition between optically thick outer region and optically thin inner region are constructed and analyzed.

We consider an accretion flow model originally proposed by Bisnovatyi-Kogan & Ruzmaikin (1974), which has been confirmed in recent 3D MHD simulations. In the model, the accreting gas drags in a strong poloidal magnetic field to the center such that the accumulated field disrupts the axisymmetric accretion flow at a relatively large radius. Inside the disruption radius, the gas accretes as discrete blobs or streams with a velocity much less than the free-fall velocity. Almost the entire rest mass energy of the gas is released as heat, radiation and mechanical/magnetic energy. Even for a non-rotating black hole, the efficiency of converting mass to energy is of order 50% or higher. The model is thus a practical analog of an idealized engine proposed by Geroch and Bekenstein.

We present three-dimensional MHD simulations of rotating radiatively inefficient accretion flows onto black holes. In the simulations, we continuously inject magnetized matter into the computational domain near the outer boundary, and we run the calculations long enough for the resulting accretion flow to reach a quasi-steady state. We have studied two limiting cases for the geometry of the injected magnetic field: pure toroidal field and pure poloidal field. In the case of toroidal field injection, the accreting matter forms a nearly axisymmetric, geometrically-thick, turbulent accretion disk. The disk resembles in many respects the convection-dominated accretion flows found in previous numerical and analytical investigations of viscous hydrodynamic flows. Models with poloidal field injection evolve through two distinct phases. In an initial transient phase, the flow forms a relatively flattened, quasi-Keplerian disk with a hot corona and a bipolar outflow. However, when the flow later achieves steady state, it changes in character completely. The magnetized accreting gas becomes two-phase, with most of the volume being dominated by a strong dipolar magnetic field from which a thermal low-density wind flows out. Accretion occurs mainly via narrow slowly-rotating radial streams which `diffuse' through the magnetic field with the help of magnetic reconnection events.

Convection-dominated accretion flow (CDAF) is a promising model to explain underluminous accreting black holes in X-ray binaries and galactic nuclei. I discuss effects of angular momentum transport in viscous hydrodynamical and MHD CDAFs. In hydrodynamical CDAFs, convection transports angular momentum inward, and this together with outward convection transport of thermal energy determine the radial structure of the flow. In MHD CDAFs, convection can transport angular momentum either inward or outward, depending on properties of turbulence in rotating magnetized plasma, which are not fully understood yet. Direction of convection angular momentum transport can affect the law of rotation of MHD CDAFs.

Radiatively inefficient accretion flows onto black holes are unstable due to both an outwardly decreasing entropy (`convection') and an outwardly decreasing rotation rate (the `magnetorotational instability'; MRI). Using a linear magnetohydrodynamic stability analysis, we show that long-wavelength modes are primarily destabilized by the entropy gradient and that such `convective' modes transport angular momentum inwards. Moreover, the stability criteria for the convective modes are the standard Hoiland criteria of hydrodynamics. By contrast, shorter wavelength modes are primarily destabilized by magnetic tension and differential rotation. These `MRI' modes transport angular momentum outwards. The convection-dominated accretion flow (CDAF) model, which has been proposed for radiatively inefficient accretion onto a black hole, posits that inward angular momentum transport and outward energy transport by long-wavelength convective fluctuations are crucial for determining the structure of the accretion flow. Our analysis suggests that the CDAF model is applicable to a magnetohydrodynamic accretion flow provided the magnetic field saturates at a sufficiently sub-equipartition value (plasma beta >> 1), so that long-wavelength convective fluctuations can fit inside the accretion disk. Numerical magnetohydrodynamic simulations are required to determine whether such a sub-equipartition field is in fact obtained.

We consider the radial structure of radiatively inefficient hydrodynamic accretion flows around black holes. We show that low-viscosity flows consist of two zones: an outer convection-dominated zone and an inner advection-dominated zone. The transition between these two zones occurs at ~50 Schwarzschild radii.

We present three-dimensional numerical magnetohydrodynamic simulations of radiatively inefficient spherical accretion onto a black hole. The simulations are initialized with a Bondi flow, and with a weak, dynamically unimportant, large-scale magnetic field. The magnetic field is amplified as the gas flows in. When the magnetic pressure approaches equipartition with the gas pressure, the field begins to reconnect and the gas is heated up. The heated gas is buoyant and moves outward, causing line stretching of the frozen-in magnetic field. This leads to further reconnection, and more heating and buoyancy-induced motions, so that the flow makes a transition to a state of self-sustained convection. The radial structure of the flow changes dramatically from its initial Bondi profile, and the mass accretion rate onto the black hole decreases significantly. Motivated by the numerical results, we develop a simplified analytical model of a radiatively inefficient spherical flow in which convective transport of energy to large radii plays an important role. In this "convection-dominated Bondi flow" the accretion velocity is highly subsonic and the density varies with radius as ~R^{-1/2} rather than the standard Bondi scaling ~R^{-3/2}. We estimate that the mass accretion rate onto the black hole is significantly less than the Bondi accretion rate. Convection-dominated Bondi flows may be relevant for understanding many astrophysical phenomena, e.g. post-supernova fallback and radiatively inefficient accretion onto supermassive black holes, stellar-mass black holes and neutron stars.

Recent hydrodynamical simulations of radiatively inefficient black hole accretion flows with low viscosity have demonstrated that these flows differ significantly from those described by an advection-dominated model. The black hole flows are advection-dominated only in their inner parts, but convectively dominated at radii R>100R_g. In such flows, the radiative output comes mostly from the convection part, and the radiative efficiency is independent of accretion rate and equals ~0.001. This value gives a limit for how dim an accreting black hole could be. It agrees with recent Chandra observations which indicate that accreting black holes in low-mass X-ray binaries are by factor about 100 dimmer that neutron stars accreting with the same accretion rates.

Recent numerical simulations of radiatively inefficient accretion flows onto compact objects have shown that convection is a general feature in such flows. Dissipation of rotational and gravitational energies in the accretion flows results in inward increase of entropy and development of efficient convective motions. Convection-dominated accretion flows (CDAFs) have a structure that is modified significantly in comparison with the canonical advection-dominated and Bondi-like accretion flows. The flows are characterized by the flattened radial density profiles, ~R^{-1/2}, and have reduced mass accretion rates. Convection transports outward a significant amount of the released binding energy of the accretion flow. We discuss basic dynamical and observational properties of ADAFs using numerical models and self-similar analytical solutions.

We argue that the proper way to treat thin-thick accretion-disk transitions should take into account the 2-D nature of the problem. We illustrate the physical inconsistency of the 1-D vertically integrated approach by discussing a particular example of the convective transport of energy.

Advection-dominated accretion flows (ADAFs) are known to be convectively unstable for low values of the viscosity parameter alpha. Two-dimensional axisymmetric hydrodynamic simulations of such flows reveal a radial density profile which is significantly flatter than the rho~r^{-3/2} expected for a canonical ADAF. The modified density profile is the result of inward transport of angular momentum by axisymmetric convective eddies. We present three-dimensional hydrodynamic simulations of convective ADAFs which are free of the assumption of axisymmetry. We find that the results are qualitatively and quantitatively similar to those obtained with two-dimensional simulations. In particular, the convective eddies are nearly axisymmetric and transport angular momentum inward.

We present a systematic numerical study of two-dimensional axisymmetric accretion flows around black holes. The flows have no radiative cooling and are treated in the framework of the hydrodynamical approximation. The models calculated in this study cover the large range of the relevant parameter space. There are four types of flows, determined by the values of the viscosity parameter $\alpha$ and the adiabatic index $\gamma$: convective flows, large-scale circulations, pure inflows and bipolar outflows. Thermal conduction introduces significant changes to the solutions, but does not create a new flow type. Convective accretion flows and flows with large-scale circulations have significant outward-directed energy fluxes, which have important implications for the spectra and luminosities of accreting black holes.

Global solutions of optically thick advective accretion disks around black holes are constructed. The solutions are obtained by solving numerically a set of ordinary differential equations corresponding to a steady axisymmetric geometrically thin disk. We pay special attention to consistently satisfy the regularity conditions at singular points of the equations. For this reason we analytically expand a solution at the singular point, and use coefficients of the expansion in our iterative numerical procedure. We obtain consistent transonic solutions in a wide range of values of the viscosity parameter alpha and mass acretion rate. We compare two different form of viscosity: one takes the shear stress to be proportional to the pressure, while the other uses the angular velocity gradient-dependent stress. We find that there are two singular points in solutions corresponding to the pressure-proportional shear stress. The inner singular point locates close to the last stable orbit around black hole. This point changes its type from a saddle to node depending on values of alpha and accretion rate. The outer singular point locates at larger radius and is always of a saddle-type. We argue that, contrary to the previous investigations, a nodal-type inner singular point does not introduce multiple solutions. Only one integral curve, which corresponds to the unique global solution, passes simultaneously the inner and outer singular points independently of the type of inner singular point. Solutions with the angular velocity gradient-dependent shear stress have one singular point which is always of a saddle-type and corresponds to the unique global solution. The structure of accretion disks corresponding to both viscosities are similar.

We consider height-integrated equations of an advection-dominated accretion flow (ADAF), assuming that there is no mass outflow. We include convection through a mixing length formalism. We seek self-similar solutions in which the rotational velocity and sound speed scale as R^{-1/2}, where R is the radius, and consider two limiting prescriptions for the transport of angular momentum by convection. In one limit, the transport occurs down the angular velocity gradient, so convection moves angular momentum outward. In the other, the transport is down the specific angular momentum gradient, so convection moves angular momentum inward. We also consider general prescriptions which lie in between the two limits. When convection moves angular momentum outward, we recover the usual self-similar solution for ADAFs in which the mass density scales as rho ~ R^{-3/2}. When convection moves angular momentum inward, the result depends on the viscosity coefficient alpha. If alpha>alpha_{crit1} ~ 0.05, we once again find the standard ADAF solution. For alpha<alpha_{crit}, however, we find a non-accreting solution in which rho ~ R^{-1/2}. We refer to this as a "convective envelope" solution or a "convection-dominated accretion flow". Two-dimensional numerical simulations of ADAFs with values of alpha<0.03 have been reported by several authors. The simulated ADAFs exhibit convection. By virtue of their axisymmetry, convection in these simulations moves angular momentum inward, as we confirm by computing the Reynolds stress. The simulations give rho ~ R^{-1/2}, in good agreement with the convective envelope solution. The R^{-1/2} density profile is not a consequence of mass outflow.

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 onto 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 onto black holes and show that such, necessarily transonic, flows may have either positive or negative Bernoulli function depending on the flow viscosity. Flows with low viscosities (alpha<0.1 in the alpha-viscosity model) have a negative Bernoulli function. Without exception, all 2-D and 1-D 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, 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 2-D 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.

Using two-dimensional simulations of non-radiative viscous rotating black hole accretion flows, we show that the flows with alpha~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 onto black hole. Observational implications are discussed.

Numerical, two-dimensional, time-dependent hydrodynamical models of geometrically thick accretion discs around black holes are presented. Accretion flows with non-effective radiation cooling (ADAFs) can be both convectively stable or unstable depending on the value of the viscosity parameter \alpha. The high viscosity flows (\alpha~1) are stable and have a strong equatorial inflow and bipolar outflows. The low viscosity flows (\alpha<0.1) are convectively unstable and this induces quasi-periodic variability.

Observed spectra of Active Galactic Nuclei (AGN) and luminous X-ray binaries in our Galaxy suggest that both hot (~10^9 K) and cold (~10^6 K) plasma components exist close to the central accreting black hole. Hard X-ray component of the spectra is usually explained by Compton upscattering of optical/UV photons from optically thick cold plasma by hot electrons. Observations also indicate that some of these objects are quite efficient in converting gravitational energy of accretion matter into radiation. Existing theoretical models have difficulties in explaining the two plasma components and high intensity of hard X-rays. Most of the models assume that the hot component emerges from the cold one due to some kind of instability, but no one offers a satisfactory physical explanation for this. Here we propose a solution to these difficulties that reverses what was imagined previously: in our model the hot component forms first and afterward it cools down to form the cold component. In our model, accretion flow has initially a small angular momentum, and thus it has a quasi-spherical geometry at large radii. Close to the black hole, the accreting matter is heated up in shocks that form due to the action of the centrifugal force. The hot post-shock matter is very efficiently cooled down by Comptonization of low energy photons and condensates into a thin and cold accretion disk. The thin disk emits the low energy photons which cool the hot component.

We numerically construct slim, global, vertically integrated models of optically thin, transonic accretion discs around black holes, assuming a regularity condition at the sonic radius and boundary conditions at the outer radius of the disc and near the black hole. In agreement with several previous studies, we find two branches of shock-free solutions, in which the cooling is dominated either by advection, or by local radiation. We also confirm that the part of the accretion flow where advection dominates is in some circumstances limited in size: it does not extend beyond a certain outer limiting radius. New results found in our paper concern the location of the limiting radius and properties of the flow near to it. In particular, we find that beyond the limiting radius, the advective dominated solutions match on to Shapiro, Lightman & Eardley (SLE) discs through a smooth transition region. Therefore, the full global solutions are shock-free and unlimited in size. There is no need for postulating an extra physical effect (e.g. evaporation) for triggering the ADAF-SLE transition. It occurs due to standard accretion processes described by the classic slim disc equations.

A new statistical approach is presented to study the thermal instability process of optically thin unmagnetized plasma. In this approach the time evolution of mass distribution function over temperature is calculated. This function characterizes the statistical properties of the multiphase medium of arbitrary spaced three-dimensional structure of arbitrary temperature perturbations. We construct our theory under the isobarical condition (P=const over space), which is satisfied in the short wavelength limit. The developed theory is illustrated in the case of thermal instability of a slowly expanding interstellar cloud. Numerical solutions of equations of the statistical theory are constucted and compared with hydrodynamical solutions. The results of both approaches are identical in the short wavelength range when the isobarity condition is satisfied. Also the limits of applicability of the statistical theory are estimated. The possible evolution of initial spectrum of perturbations is discussed. The proposed theory and numerical models can be relevant to the formation of the two-phases medium in the ~1pc region around quasars. Then small warm (T~10000K) clouds are formed as the result of thermal instability in an expanded gas fragment, which is a product of either a star-star or star-accretion disk collision.

The morphology of adiabatic accretion flows with small non-zero specific angular momentum has been investigated in the axisymmetric and nonviscous limit. For an initial state characterized by a Bondi flow with the specific angular momentum distributed with respect to polar angle, a travelling shock wave forms which propagates more rapidly in the equatorial plane than in the plane perpendicular to it resulting in the formation of a hot torus. In cases where the incoming flow is restricted to lie near the equatorial plane, a strong wind forms directed away from this plane with the tendency for the formation of a nonsteady shock structure. As the height of the incoming flow is increased it is found that the resultant wind weakens. The parameter regime which delineates accretion flows characterized by a travelling shock and a nearly standing shock structure is presented.

A simple re-scaling of velocities calculated from the Paczy{\'n}ski \& Wiita (1980) pseudo-Newtonian potential makes them consistent with special relativity and greatly improves the agreement with the exact relativistic calculations. The improvement is relevant in calculations of the Doppler effect, spectra and radiation transfer which all involve effects of special relativity.

Two-dimensional accretion flows near black holes have been investigated by time-dependent hydrodynamical calculations. We assume that the flow is axisymmetric and that radiative losses of internal energy are negligible, so that the disc is geometrically thick and hot. Accretion occurs due to the overflow of the effective potential barrier near the black hole, similar to the case of the Roche lobe overflowing star in a binary system. We make no pre-assumptions on the properties of the flow, instead our models evolve self-consistently from an initially non-accreting state. The viscosity is due to the the small-scale turbulence and it is described by the $\alpha$-viscosity prescription. We confirm earlier suggestions that viscous accretion flows are convectively unstable. We found that the instability produces transient eddies of various length-scales. The eddies contribute to the strength of the viscosity in the flow by redistributing the angular momentum. They also introduce low amplitude oscillatory variations which have a typical frequency about $100\; \msun/M$ Hz for a system of mass $M$. This may be relevant to the high frequency ($\sim 4-10$ Hz) quasi-periodic oscillations observed in the Galactic black hole candidate X-ray sources.