We present the two-dimensional inflow-outflow solutions of radiation hydrodynamic (RHD) equations of supercritical accretion flows. Compared with prior studies, we include all components of the viscous stress tensor. We assume steady state flow and use self-similar solutions in the radial direction to solve the equations in $ r-\theta $ domain of the spherical coordinates. The set of differential equations have been integrated from the rotation axis to the equatorial plane. We find that the self-similarity assumption requires that the radial profile of density is described by $ \rho(r) \propto r^{-0.5} $. Correspondingly, the radial profile of the mass inflow rate decreases with decreasing radii as $ \dot{M}_\mathrm{in} \propto r $. Inflow-outflow structure has been found in our solution. In the region $ \theta > 65^{\circ} $ there exist inflow while above that flow moves outward and outflow could launch. The driving forces of the outflow are analyzed and found that the radiation force is dominant and push the gas particles outwards with poloidal velocity $ \sim 0.25 c $. The properties of outflow are also studied. The results show that the mass flux weighted angular momentum of the inflow is lower than that of outflow, thus the angular momentum of the flow can be transported by the outflow. We also analyze the convective stability of the supercritical disk and find that in the absence of the magnetic field, the flow is convectively unstable. Our analytical results are fully consistent with the previous numerical simulations of the supercritical accretion flow.

We perform two-dimensional hydrodynamical simulations of slowly rotating accretion flows in the region of $ 0.01-7\, \mathrm{pc} $ around a supermassive black holes with $ M_\mathrm{BH} = 10^{8} M_{\odot} $. The accretion flow is irradiated by the photons from the central active galactic nucleus (AGN). In addition to the direct radiation from the AGN, we have also included the "re-radiation", i.e., the locally produced radiation by Thomson scattering, line and bremsstrahlung radiation. Compare to our previous work, we have improved the calculation of radiation force due to the Thomson scattering of X-ray photons from the central AGN. We find that this improvement can significantly increase the mass flux and velocity of outflow. We have compared the properties of outflow --- including mass outflow rate, velocity, and kinetic luminosity of outflow --- in our simulation with the observed properties of outflow in AGNs and found that they are in good consistency. This implies that the combination of line and re-radiation forces is the possible origin of observed outflow in luminous AGNs.

Observations indicate that wind can be generated in hot accretion flow. By performing numerical simulations, Yuan et al. studied the detailed properties of wind generated from weakly magnetized accretion flow. However, properties of wind generated from strongly magnetized hot accretion flow have not been studied. In this paper, we study the properties of wind generated from both weakly and strongly magnetized accretion flow. We focus on how the magnetic field strength affects the wind properties. We solve time-steady two-dimensional magnetohydrodynamic (MHD) equations of black hole accretion in the presence of large-scale magnetic field. We assume self-similarity in radial direction. The magnetic field is assumed to be evenly symmetric with the equatorial plane. We find that wind exists in both weakly and strongly magnetized accretion flow. When magnetic field is weak (magnetic pressure is more than 2 orders of magnitude smaller than gas pressure), wind is driven by gas pressure gradient and centrifugal forces. When magnetic field is strong (magnetic pressure is slightly smaller than gas pressure), wind is driven by gas pressure gradient and magnetic pressure gradient forces. The power of wind in strongly magnetized case is just slightly larger than that in weakly magnetized case. The power of wind lies in a range $ P_W \sim 10^{-4}-10^{-3} \dot M_{\rm in} c^2 $, with $ \dot M_{\rm in} $ and $ c $ being mass inflow rate and speed of light, respectively. The possible role of wind in active galactic nuclei feedback is briefly discussed.

We solved the set of two-dimensional magnetohydrodynamic (MHD) equations for optically thin black hole accretion flows incorporating toroidal component of magnetic field. Following global and local MHD simulations of black hole accretion disks, the magnetic field inside the disk is decomposed into a large scale field and a fluctuating field. The effects of the fluctuating magnetic field in transferring the angular momentum and dissipating the energy are described through the usual $ \alpha $ description. We solved the MHD equations by assuming steady state and radially self-similar approximation in $ r-\theta $ plane of spherical coordinate system. We found that as the amount of magnetic field at the equatorial plane increases, the heating by the viscosity decreases. In addition, the maximum amount of the heating by the viscous dissipation is produced at the mid-plane of the disk, while that of the heating by the magnetic field dissipation is produced at the surface of the disk. Our main conclusion is that in terms of the no-outflow solution, thermal equilibrium still exists for the strong magnetic filed at the equatorial plane of the disk.

We solve the radiation-hydrodynamic (RHD) equations of supercritical accretion flows in the presence of radiation force and outflow by using self similar solutions. Compare with the pioneer works, in this paper we consider power-law function for mass inflow rate as $ \dot{M} \propto r^{s} $. We found that $ s = 1 $ when the radiative cooling term is included in the energy equation. Correspondingly, the effective temperature profile with respect to the radius was obtained as $ T_{\text{eff}} \varpropto r^{-1/2} $. In addition, we investigated the influence of the outflow on the dynamics of the accretion flow. We also calculated the continuum spectrum emitted from the disk surface as well as the bolometric luminosity of the accretion flow. Furthermore, our results show that the advection parameter, $ f $, strongly depends on mass inflow rate.

We solve the two-dimensional magnetohydrodynamic (MHD) equations of black hole accretion with the presence of magnetic field. The field includes a turbulent component, whose role is represented by the viscosity, and a large-scale ordered component. The latter is further assumed to be evenly symmetric with the equatorial plane. The equations are solved in the $r-\theta$ plane of a spherical coordinate by assuming time-steady and radially self-similar. An inflow-wind solution is found. Around the equatorial plane, the gas is inflowing; while above and below the equatorial plane at a certain critical $\theta$ angle, $\theta\sim 47^{\circ}$, the inflow changes its direction of radial motion and becomes wind. The driving forces are analyzed and found to be the centrifugal force and the gradient of gas and magnetic pressure. The properties of wind are also calculated. The specific angular momentum of wind is found to be significantly larger than that of inflow, thus wind can transfer angular momentum outward. These analytical results are compared to those obtained by the trajectory analysis based on MHD numerical simulation data and good agreements are found.

The main aim of this paper is studying the effect of toroidal magnetic field on the structure of Advection-Dominated Accretion Flows (ADAF) in the presence of the turbulence viscosity and diffusivity due to viscosity and magnetic field respectively. We use self-similar assumption in radial direction to solve the magnetohydrodynamic (MHD) equations for hot accretion disk. We use spherical coordinate $ (r, \theta, \varphi) $ to solve our equation. The toroidal component of magnetic field is considered and all three components of the velocity field $ \mathbf{v}\equiv (v_{r}, v_{\theta}, v_{\varphi}) $ are present in our work. We reduce the equations to a set of differential equations about $ \theta $ and apply the symmetric boundary condition at the equatorial plane of the disk. Our results indicate that the outflow region, where the redial velocity becomes positive in a certain inclination angle $ \theta_{0} $, always exist. The results represent that the stronger the magnetic field, the smaller the inclination angle, $ \theta_{0} $ becomes. It means that a magnetized disk is thinner compared to a non-magnetized disk. According to the work by \citealt{jw}, we can define three regions. The first one is called inflow region, which starts from the disk midplane to a certain inclination $ \theta_0 $ where $ v_{r}(\theta_{0}) = 0 $. In this region, the velocity has a negative value and the accretion material moves toward the central object. The outflow region, where $ v_{r}(\theta) > 0 $, is placed between $ \theta_{0} $ and surface of the disk, $ \theta_{0} < \theta < \theta_{s} $. In this area, the accretion flow moves away from the central object. The third region, which is located between the surface of the disk and the polar axis, is called wind region.

We have explored the structure of hot flow bathed in a general large-scale magnetic field. The importance of outflow and thermal conduction on the self-similar structure of a hot accretion flows has been investigated. We consider the additional magnetic parameters $ \beta_{r,\varphi,z}\big[= c^2_{r,\varphi,z}/(2 c^2_{s}) \big] $, where $ c^2_{r,\varphi,z} $ are the Alfv$\acute{e}$n sound speeds in three direction of cylindrical coordinate. In comparison to the accretion disk without winds, our results show that the radial and rotational velocities of the disk become faster however it become cooler because of the angular momentum and energy flux which are taking away by the winds. but thermal conduction opposes the effect of winds not only decrease the rotational velocity but also increase the radial velocity as well as the sound speed of the disk. In addition we study the effect of global magnetic field on the structure of the disk. Our numerical results show that all components of magnetic field can be important and they have a considerable effect on velocities and vertical structure of the disk.

The existence of outflow and magnetic field in the inner region of hot accretion flows have been confirmed by observations and numerical magnetohydrodynamic simulations (MHD). We present self-similar solutions for radiation inefficiently accretion flows (RIAF) around black holes in the presence of outflow and global magnetic field. The influence of outflow is taken into account by adopting a radius dependent of mass accretion rate $ \dot{M} = \dot{M}_{0}(r/r_{0})^{s} $ with $ s > 0 $. Also we consider convection through a mixing length formalism to calculate convection parameter $ \alpha_{con} $. Moreover we consider the additional magnetic field parameters $ \beta_{r,\varphi,z}\big[= c^2_{r,\varphi,z}/(2 c^2_{s}) \big] $, where $ c^2_{r,\varphi,z} $ are the Alfv$\acute{e}$n sound speeds in three direction of cylindrical coordinate. Our numerical results show that by increasing all components of magnetic field, the surface density and rotational velocity increase although the sound speed and radial infall velocity of the disc decrease. Also we have found out that the existence of wind will lead to reduction of surface density as well as rotational velocity. Moreover the radial velocity, sound speed, advection parameter and the vertical thickness of the disc will increase when outflow becomes important in the radiation inefficiently accretion flow.

We have studied the structure of hot accretion flow bathed in a general large-scale magnetic field. We have considered magnetic parameters $ \beta_{r,\varphi,z}[=c^2_{r,\varphi,z}/(2c^2_{s})] $, where $ c^2_{r, \varphi, z} $ are the Alfv\'{e}n sound speeds in three direction of cylindrical coordinate $ (r,\varphi,z) $. The dominant mechanism of energy dissipation is assumed to be the magnetic diffusivity due to turbulence and viscosity in the accretion flow. Also, we adopt a more realistic model for kinematic viscosity $ (\nu=\alpha c_{s} H) $, with both $ c_{s} $ and $ H $ as a function of magnetic field. As a result in our model, the kinematic viscosity and magnetic diffusivity $ (\eta=\eta_{0}c_{s} H) $ are not constant. In order to solve the integrated equations that govern the behavior of the accretion flow, a self-similar method is used. It is found that the existence of magnetic resistivity will increase the radial infall velocity as well as sound speed and vertical thickness of the disk. However the rotational velocity of the disk decreases by the increase of magnetic resistivity. Moreover, we study the effect of three components of global magnetic field on the structure of the disk. We found out that the radial velocity and sound speed are Sub-Keplerian for all values of magnetic field parameters, but the rotational velocity can be Super-Keplerian by the increase of toroidal magnetic field. Also, Our numerical results show that all components of magnetic field can be important and have a considerable effect on velocities and vertical thickness of the disk.

Force due to the self-gravity of the disc in the vertical direction is considered to study its possible effects on the structure of a magnetized advection-dominated accretion disc. We present steady-sate self similar solutions for the dynamical structure of such a type of the accretion flows. Our solutions imply reduced thickness of the disc because of the self-gravity. It also imply that the thickness of the disc will increase by adding the magnetic field strength.