A magnetohydrodynamical model for the formation of episodic jets

Episodic ejection of plasma blobs has been observed in many black hole systems. While steady, continuous jets are believed to be associated with large-scale open magnetic fields, what causes the episodic ejection of blobs remains unclear. Here by analogy with the coronal mass ejection on the Sun, we propose a magnetohydrodynamical model for episodic ejections from black holes associated with the closed magnetic fields in an accretion flow. Shear and turbulence of the accretion flow deform the field and result in the formation of a flux rope in the disc corona. Energy and helicity are accumulated and stored until a threshold is reached. The system then loses its equilibrium and the flux rope is thrust outward by the magnetic compression force in a catastrophic way. Our calculations show that for parameters appropriate for the black hole in our Galactic centre, the plasmoid can attain relativistic speeds in about 35 min.     

Magnetohydrodynamical non-radiative accretion flows in two dimensions

We present the results of axisymmetric, time-dependent magnetohydrodynamic simulations of accretion flows around black holes. The calculations begin from a rotationally supported thick torus which contains a weak poloidal field. Accretion is produced by growth and saturation of the magnetorotational instability (MRI) provided that the wavelength of the fastest growing mode is less than the thickness of the torus. Using a computational grid that spans more than two decades in radius, we compare the time-averaged properties of the flow with previous hydrodynamical simulations. The net mass accretion rate is small compared with the mass inflow and outflow rates at large radii associated with turbulent eddies. Turbulence is driven by the MRI rather than convection. The two-dimensional structure of the time-averaged flow is significantly different compared with the hydrodynamical case. We discuss the limitations imposed on our results by the assumption of axisymmetry and the relatively small radial domain.     

A semi-analytic solution for the radial and vertical structure of accretion discs with a magnetic wind

A semi-analytic method is presented for solving for the radial and vertical structures of an accretion disc, with a magnetically channelled wind flowing from its surfaces. Both magnetic and turbulent viscous effects are taken into account, and the essential wind properties are related to the disc structure. The angular momentum removed by the wind plays a major part in driving the inflow through the disc, with photospheric temperatures being sufficient to generate the required wind mass flux. The magnetic field is generated by an αω-dynamo, but the method of solution should have application with other magnetic field sources. Self-consistent disc-wind solutions result, with rms turbulent Mach numbers which are in good agreement with those found in simulations of turbulence generated from magnetic shearing instabilities.     

Studies of accretion flows around rotating black holes – II. Standing shocks in the pseudo-Kerr geometry

Standing, propagating or oscillating shock waves are common in accretion and winds around compact objects. We study the topology of all possible solutions using the pseudo-Kerr geometry. We present the parameter space spanned by the specific energy and angular momentum and compare it with that obtained from the full general relativity to show that the potential can work satisfactorily in fluid dynamics also, provided the polytropic index is suitably modified. We then divide the parameter space depending on the nature of the solution topology. We specifically study the nature of the standing Rankine–Hugoniot shocks. We also show that as the Kerr parameter is increased, the shock location generally moves closer to the black hole. In future, these solutions can be used as guidelines to test numerical simulations around compact objects.     

Mass loss from the inner regions of accretion discs due to centrifugally driven magnetic wind flows

Wind flows and collimated jets are believed to be a feature of a range of disc accreting systems. These include active galactic nuclei, T Tauri stars, X-ray binaries and cataclysmic variables. The observed collimation implies large-scale magnetic fields and it is known that dipole-symmetry fields of sufficient strength can channel wind flows emanating from the surfaces of a disc. The disc inflow leads to the bending of the poloidal magnetic field lines, and centrifugally driven magnetic winds can be launched when the bending exceeds a critical value. Such winds can result in angular momentum transport at least as effective as turbulent viscosity, and hence they can play a major part in driving the disc inflow.
It is shown here that if the standard boundary condition of vanishing viscous stress close to the stellar surface is applied, together with the standard connection between viscosity and magnetic diffusivity, then poloidal magnetic field bending increases as the star is approached with a corresponding increase in the wind mass loss rate. A significant amount of material can be lost from the system via the enhanced wind from a narrow region close to the stellar surface. This occurs for a Keplerian angular velocity distribution and for a modified form of angular velocity, which allows for matching of the disc and stellar rotation rates through a boundary layer above the stellar surface. The enhanced mass loss is significantly affected by the behaviour of the disc angular velocity as the stellar surface is approached, and hence by the stellar rotation rate. Such a mechanism may be related to the production of jets from the inner regions of disc accreting systems.     

Stability of accretion discs threaded by a strong magnetic field

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The three-dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou . In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outwardly and allows mass to accrete inwardly. The instability takes the form of a m =1 rotating 'crescent', reminiscent of the purely hydrodynamic non-linear instability previously found in pressure-supported discs. A model where the initial surface mass density Σ( r ) and B _z ( r ) decrease with radius as power laws shows transient mass accretion during about six orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.     

An accretion disc model with a magnetic wind and turbulent viscosity

A model is presented for an accretion disc with turbulent viscosity and a magnetically influenced wind. The magnetic field is generated by a dynamo in the disc, involving the turbulence and radial shear. Disc-wind solutions are found for which the wind mass flux is sufficient to play a major part in driving an imposed steady inflow, but small enough for most material to be accreted on to the central object. Constraints arise for the magnetic Reynolds and Prandtl numbers in terms of the turbulent Mach number and vertical length-scale of the disc's horizontal magnetic field. It is shown that the imposition of a stellar boundary condition enhances the wind mass flux in the very inner region of the disc and may result in jet formation.     

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.     

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.     

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.     

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.