We perform two-dimensional hydrodynamic simulations to study the non-rotating gas flow in the region from 1-1000 parsec. This region connects the galaxy scale and the accretion disk scale. The gas is irradiated by the photons generated by the accretion disk system near the central black hole with $10^8$ solar mass. We assume that the luminosity of the central accretion system is constant and equals to $0.3$ Eddington luminosity. Gas with density higher than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally unstable. Two phase, cold and hot, gas flow will form due to thermal instability. We calculated the ratio of cold gas accretion rate to hot gas accretion rate. This ratio is a function of gas density and generally increases with the increase of gas density. The maximum value of this ratio is $16$, when gas density is $10^{-21} {\rm g \ cm^{-3}} $. Gas with density lower than $10^{-24} {\rm g \ cm^{-3}}$ is found to be thermally stable and the gas flow is one-phase. The applications of the results are briefly discussed.

Numerical simulations of radiative two-temperature hot accretion flows (HAFs) around Neutron stars (NSs) are performed. We assume that all of the energy carried by the HAF around a NS will be thermalized and radiated out at the surface of the NS. The thermal photons will propagate outwards radially and cool the HAF vis Comptonization. We define $\dot m$ as the mass accretion rate at the surface of the central object in unit of Eddington accretion rate ($\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ being Eddington luminosity and speed of light, respectively). When $\dot m$ is lower than $\sim 10^{-4}$, the cooling of the HAF is not important and outflows are very strong. When $\dot m > \sim 10^{-3}$, cooling becomes important and outflows are significantly weak. In the range $10^{-4} < \dot m < 10^{-3}$, the HAFs transients from a strong outflow phase to a very weak outflow phase with increase of $\dot m$. The properties of the HAF around a NS are also compared to those of the HAF around a BH. We find that with a similar $\dot m$, the dynamical properties of the HAF around a NS are quite similar as those of the HAF around a BH. However, the emitted spectrum of a HAF around a NS can be quite different from that of a HAF around a BH due to the presence of a thermal soft X-ray component coming from the surface of the NS.

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.

A hot corona is suggested to be above the standard thin disk. The anisotropy of hard X-ray emission in radio-quiet active galactic nuclei implies that the corona is not static and probably moves outwards like winds. We perform two-dimensional magnetohydrodynamical simulations to study the outflowing corona driven by magnetic field and radiation force. In our simulations, as the initial state and the boundary condition at the disk surface, the corona temperature is set to $10^9$ K inside 10 Schwarzschild radius ($r_{\rm s}$) while the corona temperature is set to $10^7$ K a 10 $r_{\rm s}$. We employ a weak poloidal magnetic field as the initial magnetic field. A collimated outflow and a wide-angle ordered outflow are observed in our simulations. The collimated outflow is around the rotational axis and has a bulk velocity of $\sim$0.03--0.3 $c$ ( $c$ is speed of light) at 90 $r_{\rm s}$, while their mass outflow rate is very low. The collimated outflow is a weak jet. The wide-angle ordered outflow is distributed at the middle and high latitudes and moves outwards with a velocity of $10^2$--$10^4$ km s$^{-1}$. The outflow velocity depends on the disk luminosity. The gas around the disk surface is turbulent, especially outside of 10 $r_{\rm s}$. The other properties of outflows are discussed in detail.

We perform simulations to study the effects of active galactic nuclei (AGNs) radiation and wind feedback on the properties of slowly rotating accretion flow at parsec scale. We find that when only radiative feedback is considered, outflows can be produced by the radiation pressure due to Thomson scattering. The mass flux of outflow is comparable to that of inflow. Although strong outflow is present, the luminosity of AGN can be easily super-Eddington. When wind feedback is also taken into account, the mass flux of outflow does not change much. Consequently, the luminosity of the central AGN can still be super-Eddington. However, observations show that the luminosity of most AGNs is sub-Eddington. Some other mechanisms are needed to reduce the AGNs luminosity. Although the mass outflow rate does not change much by wind feedback, other properties of outflow (density, temperature, velocity and kinetic power) can be significantly changed by wind feedback. In the presence of wind feedback, the density of outflow becomes significantly lower, temperature of outflow becomes significantly higher, velocity of outflow is increased by one order of magnitude, the kinetic power of outflow is increased by a factor of $40-100$.

We perform as the first time hydrodynamic simulations to study the properties of hot accretion flow (HAF) around a neutron star (NS). The energy carried by the HAF will eventually be radiated out at the surface of the NS. The emitted photons can propagate inside the HAF and cool the HAF via Comptonization. We find that the Compton cooling can affect the properties of HAF around a NS significantly. We define the Eddington accretion rate as $\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ being the Eddington luminosity and the speed of light, respectively. We define $\dot m$ as the mass accretion rate at the NS surface in unit of $\dot M_{\rm Edd}$. When $\dot m > 10^{-4}$, Compton cooling can effectively cool the HAF and suppress wind. Therefore, the mass accretion rate is almost a constant with radius. The density profile is $\rho \propto r^{-1.4}$. When $\dot m < 10^{-4}$, the Compton cooling effects become weaker, wind becomes stronger, accretion rate is proportional to $r^{0.3-0.5}$. Consequently, the density profile becomes flatter, $\rho \propto r^{-1 \sim -0.8}$. When $\dot m < 10^{-6}$, the Compton cooling effects can be neglected. We find that with a same accretion rate, the temperature of HAF around a NS is significantly lower than that of HAF around a black hole (BH). Also, the Compton $y-$parameter of HAF around a NS is significantly smaller than that of HAF around a BH. This result predicts that HAF around a NS will produce a softer spectrum compared to HAF around a BH, which is consistent with observations.

Mass accretion rate is a key parameter in accretion disk theory. It determines black hole accretion mode. In large scale cosmological simulations studying galaxy formation and evolution, Bondi radius can at most be marginally resolved. In those simulations, Bondi accretion formula is always used to estimate black hole accretion rate. Bondi solution is too simple, which cannot represent the real accretion process. We perform simulations of hot accretion flow at parsec scale irradiated by a low luminosity active galactic nucleus (AGN). We perform 77 simulations with varying density and temperature at outer boundary (10 parsec). Our purpose is to find a formula to calculate real black hole accretion rate based on the gas density and temperature at parsec scale. We define Eddington accretion rate to be $\dot M_{\rm Edd}=10L_{\rm Edd}/c^2$, with $L_{\rm Edd}$ and $c$ been Eddington luminosity and speed of light respectively. We set black hole mass to be $10^8M_{\odot}$, $M_{\odot}$ is solar mass. We find that black hole accretion rate can be expressed as $\frac{\dot M_{\rm BH} }{\dot M_{\rm Edd}}=10^{-3.11} (\rho_0/10^{-22}{\rm g cm}^{-3})^{1.36} (T_0/10^7 {\rm K})^{-1.9}$, with $\rho_0$ and $T_0$ being density and temperature at parsec scale, respectively. We find the formula can accurately predict the luminosity of observed low-luminosity AGNs (with black hole mass $\sim 10^8M_{\odot}$). This formula can be used in the sub-grid models in large scale cosmological simulations with a black hole mass of $\sim 10^8M_{\odot}$.

Observations of many dim galactic nuclei in local universe give good estimation of gas density and temperature at Bondi radius. If we assume the black hole accretes at Bondi accretion rate and radiates at efficiency of low-luminosity hot accretion flow, the predicted nuclei luminosity can be significantly higher than observations. Therefore, the real black hole mass accretion rate in these sources may be significantly smaller than the Bondi value. Active galactic nuclei (AGN) feedback may be responsible for decreasing black hole accretion rate to values much smaller than Bondi rate. We perform two-dimensional simulations of low angular momentum accretion flow at parsec and sub-parsec scales around low-luminosity active galactic nuclei (LLAGNs). We take into account the radiation and wind feedbacks of the LLAGN. The cross section of particle-particle interaction can be several orders of magnitude larger than that of photon-particle interaction. Therefore, we find that for the LLAGNs, effects of radiation feedback in decreasing black hole accretion rate are small. However, wind feedback can effectively decrease the black hole mass accretion rate. Due to the decrease of accretion rate, the black hole luminosity can be decreased by a factor of $\sim 33-400$. The results may be useful to explain why many galactic nuclei in the local universe are so dim.

Observations show that winds can be driven from the innermost region (inside a 50 Schwarschild radius) of a thin disk. It is interesting to study the winds launched from the innermost region. A hot corona above the black hole (BH) thin disk is irradiated by the disk. We perform two-dimensional hydrodynamical simulations to study the winds driven by radiation force from the corona in the innermost regions. The hard X-ray spectrum from active galactic nuclei (AGNs) suggests that the corona temperature is about $10^9$ K, so that we mainly analyze the properties of winds (or outflows) from the $10^9$ K corona. The disk luminosity plays an important role in driving the outflows. The more luminous the disk, the stronger the outflows. Mass outflow rate ($\dot{M}_{\rm out}$) at a 90 Schwarschild radius depends on disk luminosity, which can be described as $\dot{M}_{\rm out}\propto 10^{3.3 \Gamma}$ ($\Gamma$ is the ratio of the disk luminosity to the Eddington luminosity). In the case of high luminosity (e.g. $\Gamma=0.75$), the supersonic outflows with maximum speed $1.0 \times 10^4$ Km s$^{-1}$ are launched at $\sim17^{o}$ --$30^{o}$ and $\sim50^{o}$ --$80^{o}$ away from the pole axis. The Bernoulli parameter keeps increasing with the outward propagation of outflows. The radiation force keeps accelerating the outflows when outflows move outward. Therefore, we can expect the outflows to escape from the BH gravity and go to the galactic scale. The interaction between outflows and interstellar medium may be an important AGN feedback process.

It is often suggested that a super massive black hole is embedded in a nuclear bulge of size of a few $10^2$ $\text{ parsec}$ . The nuclear stars gravity is not negligible near $\sim 10 \text{ parsec}$. In order to study the effect of nuclear stars gravity on quasar radiation feedback on the parsec scale, we have simulated the parsec scale flows irradiated by a quasar by taking into account the gravitational potential of both the black hole and the nuclear star cluster. We find that the effect of nuclear stars gravity on the parsec-scale flows is related to the fraction of X-ray photons in quasar radiation. For the models in which the fraction of X-ray photons is not small (e.g. the X-ray photons contribute to $20\%$ of the quasar radiation), the nuclear stars gravity is very helpful to collimate the outflows driven by UV photons, significantly weakens the outflow power at the outer boundary and significantly enhances the net accretion rate onto the black hole. For the models in which X-ray photons are significantly decreased (e.g. the X-ray photons contribute to $5\%$ of the quasar radiation), the nuclear stars gravity can just slightly change properties of outflow and slightly enhance the net accretion rate onto the black hole.

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 study the effects of accretion environment (gas density, temperature and angular momentum) at large radii ($\sim 10$pc) on luminosity of hot accretion flows. The radiative feedback effects from the accretion flow on the accretion environment are also self-consistently taken into account. We find that the slowly rotating flows at large radii can significantly deviate from Bondi accretion when radiation heating and cooling are considered. We further find that when the temperature of environment gas is low (e.g. $T=2\times 10^7$K), the luminosity of hot accretion flows is high. When the temperature of gas is high (e.g. $T\geq4\times 10^7$K), the luminosity of hot accretion flow significantly deceases. The environment gas density can also significantly influence the luminosity of accretion flows. When density is higher than $\sim 4\times 10^{-22}\text{g} \text{cm}^{-3}$ and temperature is lower than $2\times 10^7$K, hot accretion flow with luminosity lower than $2\%L_{\text{Edd}}$ is not present. Therefore, the pc-scale environment density and temperature are two important parameters to determine the luminosity. The results are also useful for the sub-grid models adopted by the cosmological simulations.

We study slowly rotating accretion flow at parsec and sub-parsec scale irradiated by a low luminosity active galactic nuclei. We take into account the Compton heating, photoionization heating by the central X-rays. The bremsstrahlung cooling, recombination and line cooling are also included. We find that due to the Compton heating, wind can be thermally driven. The power of wind is in the range $(10^{-6}-10^{-3}) L_{\rm Edd}$, with $L_{\rm Edd}$ being the Eddington luminosity. The mass flux of wind is in the range $(0.01-1) \dot M_{\rm Edd}$ ($\dot M_{\rm Edd}= L_{\rm Edd}/0.1c^2$ is the Eddington accretion rate, $c$ is speed of light). We define the wind generation efficiency as $\epsilon = P_W/\dot {M}_{\rm BH}c^2$, with $P_W$ being wind power, $\dot M_{\rm BH}$ being the mass accretion rate onto the black hole. $\epsilon$ lies in the rage $10^{-4}-1.18$. Wind production efficiency decreases with increasing mass accretion rate. The possible role of the thermally driven wind in the active galactic feedback is briefly discussed.

We perform two-dimensional simulations to study how the wind strength changes with accretion rate. We take into account bremsstrahlung, synchrotron radiation and the Comptonization. We find that when the accretion rate is low, radiative cooling is not important, the accretion flow is hot accretion flow. For the hot accretion flow, wind is very strong. The mass flux of wind can be $\sim 50\%$ of the mass inflow rate. When the accretion rate increases to a value at which radiative cooling rate is roughly equal to or slightly larger than viscous heating rate, cold clumps can form around the equatorial plane. In this case, the gas pressure gradient force is small and wind is very weak. Our results may be useful for the sub-grid model of active galactic nuclear feedback study.

Numerical simulations of black hole hot accretion flows have shown the existence of strong wind. Those works focus only on the region close to black hole thus it is unknown whether or where the wind production stops at large radii. To address this question, Bu et al. (2016) have performed hydrodynamic (HD) simulations by taking into account the gravitational potential of both the black hole and the nuclear star clusters. The latter is assumed to be $\propto \sigma^2 \ln(r)$, with $\sigma$ being the velocity dispersion of stars and $r$ be the distance from the center of the galaxy. It was found that when the gravity is dominated by nuclear stars, i.e., outside of radius $R_A\equiv GM_{\rm BH}/\sigma^2$, winds can no longer be produced. That work, however, neglects the magnetic field, which is believed to play a crucial dynamical role in the accretion and thus must be taken into account. In this paper, we revisit this problem by performing magneto-hydrodynamical (MHD) simulations. We confirm the result of Bu et al. (2016), namely wind can't be produced at the region of $R>R_A$. Our result, combined with the results of Yuan et al. (2015), indicates that the formula describing the mass flux of wind $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$ can only be applied to the region where the black hole potential is dominant. Here $\dot{M}_{\rm BH}$ is the mass accretion rate at the black hole horizon and the value of $R_A$ is similar to the Bondi radius.

Previous works have clearly shown the existence of winds from black hole hot accretion flow and investigated their detailed properties. In extremely low accretion rate systems, the collisional mean-free path of electrons is large compared with the length-scale of the system, thus thermal conduction is dynamically important. When the magnetic field is present, the thermal conduction is anisotropic and energy transport is along magnetic field lines. In this paper, we study the effects of anisotropic thermal conduction on the wind production in hot accretion flows by performing two-dimensional magnetohydrodynamic simulations. We find that thermal conduction has only moderate effects on the mass flux of wind. But the energy flux of wind can be increased by a factor of $\sim 10$ due to the increase of wind velocity when thermal conduction is included. The increase of wind velocity is because of the increase of driving forces (e.g. gas pressure gradient force and centrifugal force) when thermal conduction is included. This result demonstrates that thermal conduction plays an important role in determining the properties of wind.

Previous works show strong winds exist in hot accretion flows around black holes. Those works focus only on the region close to the black hole thus it is unknown whether or where the wind production stops at large radii. In this paper, we investigate this problem by hydrodynamical simulations. We take into account gravities of both the black hole and the nuclear star clusters. For the latter, we assume that the velocity dispersion of stars is a constant and its gravitational potential $\propto \sigma^2 \ln (r)$, where $\sigma$ is the velocity dispersion of stars and $r$ is the distance from the center of the galaxy. We focus on the region where the gravitational potential is dominated by the star cluster. We find, same as the accretion flow at small radii, the mass inflow rate decreases inward and the flow is convectively unstable. However, trajectory analysis shows that there is very few wind launched from the flow. Our result, combined with the results of Yuan et al. (2015), indicates that the mass flux of wind launched from hot accretion flow $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$, with $r\la R_A\equiv GM_{\rm BH}/\sigma^2$. Here $\dot{M}_{\rm BH}$ is accretion rate at black hole horizon. $R_A$ is similar to Bondi radius. We argue that the inward decrease of inflow rate is not because of mass loss via wind, but because of convective motion. The disappearance of wind outside of $R_A$ must be because of the change of the gravitational potential, but the exact reason remains to be probed.

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.

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind, such as mass flux and poloidal speed, and the mechanism of wind production. For this aim, we make use of a three dimensional GRMHD simulation of hot accretion flows around a Schwarzschild black hole. The simulation is designed so that the magnetic flux is not accumulated significantly around the black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Largrangian particles from simulation data. We find two types of real outflows, i.e., a quasi-relativistic jet close to the axis and a sub-relativistic wind subtending a much larger solid angle. Most of the wind originates from the surface layer of the accretion flow. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows $v_{\rm p, wind}(r)\approx 0.25 v_k(r)$. The mass flux of jet is much lower but the speed is much higher, $v_{\rm p,jet}\sim (0.3-0.4) c$. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. We find that the wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by magnetic pressure gradient. Finally, we find that the wind production efficiency $\epsilon_{\rm wind}\equiv\dot{E}_{\rm wind}/\dot{M}_{\rm BH}c^2\sim 1/1000$, in good agreement with the value required from large-scale galaxy simulations with AGN feedback.

A pair of giant gamma-ray bubbles have been revealed by the {\it Fermi} LAT. In this paper we investigate their formation mechanism. Observations have indicated that the activity of the supermassive black hole located at the Galactic center, Sgr A*, was much stronger than the present time. Specifically, one possibility is that while Sgr A* was also in the hot accretion regime, the accretion rate should be $10^3-10^4$ times higher during the past $\sim 10^7$ yr. On the other hand, recent MHD numerical simulations of hot accretion flows have unambiguously shown the existence of strong winds and obtained their properties. Based on these knowledge, by performing three-dimensional hydrodynamical simulations, we show in this paper that the Fermi bubbles could be inflated by winds launched from the ``past' hot accretion flow in Sgr A*. In our model, the active phase of Sgr A* is required to last for about 10 million years and it was quenched no more than 0.2 million years ago. The Central Molecular Zone (CMZ) is included and it collimates the wind orientation towards the Galactic poles. Viscosity suppresses the Rayleigh-Taylor and Kelvin-Helmholtz instabilities and results in the smoothness of the bubble edge. The main observational features of the bubbles can be well explained. Specifically, the {\it ROSAT} X-ray features are interpreted by the shocked interstellar medium and the interaction region between winds and CMZ gas. The thermal pressure and temperature obtained in our model are in good consistency with the recent {\it Suzaku} observations.

The regular satellites found around Neptune ($\approx 17~M_{\Earth}$) and Uranus ($\approx 14.5~M_{\Earth}$) suggest that past gaseous circumplanetary disks may have co-existed with solids around sub-Neptune-mass protoplanets ($< 17~M_{\Earth}$). These disks have been shown to be cool, optically thin, quiescent, with low surface density and low viscosity. Numerical studies of the formation are difficult and technically challenging. As an introductory attempt, three-dimensional global simulations are performed to explore the formation of circumplanetary disks around sub-Neptune-mass protoplanets embedded within an isothermal protoplanetary disk at the inviscid limit of the fluid in the absence of self-gravity. Under such conditions, a sub-Neptune-mass protoplanet can reasonably have a rotationally supported circumplanetary disk. The size of the circumplanetary disk is found to be roughly one-tenth of the corresponding Hill radius, which is consistent with the orbital radii of irregular satellites found for Uranus. The protoplanetary gas accretes onto the circumplanetary disk vertically from high altitude and returns to the protoplanetary disk again near the midplane. This implies an open system in which the circumplanetary disk constantly exchanges angular momentum and material with its surrounding prenatal protoplanetary gas.

If the specific angular momentum of accretion gas at large radius is small compared to the local Keplerian value, one usually believes that there exists a "circularization radius" beyond which the angular momentum of accretion flow is almost a constant while within which a disk is formed and the angular momentum roughly follows the Keplerian distribution. In this paper, we perform numerical simulations to study whether the picture above is correct in the context of hot accretion flow. We find that for a steady accretion flow, the "circularization radius" does not exist and the angular momentum profile will be smooth throughout the flow. However, for transient accretion systems, such as the tidal disruption of a star by a black hole, a "turning point" should exist in the radial profile of the angular momentum, which is conceptually similar to the "circularization radius". At this radius, the viscous timescale equals the life time of the accretion event. The specific angular momentum is close to Keplerian within this radius, while beyond this radius the angular momentum is roughly constant.

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{\theta\phi}$ component of the viscous stress and consider various values of the viscous parameter $\alpha$. We find that when $T_{\theta\phi}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $\alpha$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $\rho(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $\alpha$ ranging from $\alpha=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $\alpha$, the flow is marginally stable (when $\alpha$ is small) or unstable (when $\alpha$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.

We study the effects of initial and boundary conditions, taking two-dimensional hydrodynamical numerical simulations of hot accretion flow as an example. The initial conditions considered include a rotating torus, a solution expanded from the one-dimensional global solution of hot accretion flows, injected gas with various angular momentum distributions, and the gas from a large-scale numerical simulation. Special attention is paid to the radial profiles of the mass accretion rate and density. Both can be described by a power-law function, $\dot{M}\propto r^s$ and $\rho\propto r^{-p}$. We find that if the angular momentum is not very low, the value of $s$ is not sensitive to the initial condition and lies within a narrow range, $0.47\la s \la 0.55$. However, the value of $p$ is more sensitive to the initial condition and lies in the range $0.48\la p \la 0.8$. The diversity of the density profile is because different initial conditions give different radial profiles of radial velocity due to the different angular momentum of the initial conditions. When the angular momentum of the accretion flow is very low, the inflow rate is constant with radius. Taking the torus model as an example, we have also investigated the effects of inner and outer boundary conditions by considering the widely adopted "outflow" boundary condition and the "mass flux conservation" condition. We find that the results are not sensitive to these two boundary conditions.

Numerical simulations of hot accretion flows have shown that the mass accretion rate decreases with decreasing radius. Two models have been proposed to explain this result. In the adiabatic inflow-outflow solution (ADIOS), it is thought to be due to the loss of gas in outflows. In the convection-dominated accretion flow (CDAF) model, it is explained as because that the gas is locked in convective eddies. In this paper we use hydrodynamical (HD) and magnetohydrodynamical (MHD) simulations to investigate which one is physical. We calculate and compare various properties of inflow (gas with an inward velocity) and outflow (gas with an outward velocity). Systematic and significant differences are found. For example, for HD flows, the temperature of outflow is higher than inflow; while for MHD flows, the specific angular momentum of outflow is much higher than inflow. We have also analyzed the convective stability of MHD accretion flow and found that they are stable. These results suggest that systematic inward and outward motion must exist, i.e., the ADIOS model is favored. The different properties of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD flows are buoyancy associated with the convection and the centrifugal force associated with the angular momentum transport mediated by the magnetic field, respectively. The latter mechanism is similar to the Blandford & Payne mechanism but no large-scale open magnetic field is required here. Possible observational applications are briefly discussed.

Previous hydrodynamical (HD) and magnetohydrodynamical (MHD) numerical simulations of hot accretion flows have shown that the mass accretion rate decreases with decreasing radius. Two models have been proposed to explain this result. In the ADIOS model, the inward decrease of accretion rate is because of the loss of gas in the outflow. In the CDAF model, the gas is assumed to be locked in convective eddies, which results in the inward decrease of the accretion rate. We investigate the nature of inward decrease of accretion rate using HD and MHD simulations. We calculate various properties of inflow and outflow, including the mass flux, radial and rotational velocities, temperature, and the Bernoulli parameter ($Be$). Systematic and significant differences between inflow and outflow are found. These results suggest that the inflow and outflow are not dominated by convective turbulence, but are systematic inward and outward motion. We have also analyzed the convective stability of MHD accretion flow and found that they are convectively stable. These results indicate that the ADIOS scenario is favored. The different properties of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD flows are buoyancy and centrifugal force associated with the magnetic field, respectively. The latter mechanism is similar to the Blandford & Payne mechanism. We also study the effect of initial conditions in the simulations. We find that the value of $Be$ is mainly determined by the value of $Be$ of the initial condition. We discuss some possible observational applications of our outflow model. These observations include the Fermi bubble observed in the Galaxy center, and winds widely observed in AGNs and black hole X-ray binaries.

Numerical simulations of hot accretion flow have shown that the mass accretion rate decreases with decreasing radius; consequently the density profile of accretion flow becomes flatter compared to the case of a constant accretion rate. This result has important theoretical and observational implications. However, because of technical difficulties, the radial dynamic range in almost all previous simulations usually spans at most two orders of magnitude. This small dynamical range, combined with the effects of boundary conditions, makes the simulation results suspectable. Especially, the radial profiles of density and accretion rate may not be precise enough to be used to compare with observations. In this paper we present a "two-zone" approach to expand the radial dynamical range from two to four orders of magnitude. We confirm previous results and find that from $r_s$ to $ 10^4r_s$ the radial profiles of accretion rate and density can be well described by $\dot{M}(r)\propto r^s$ and $\rho\propto r^{-p}$. The values of (s, p) are (0.48, 0.65) and (0.4, 0.85), for viscous parameter $\alpha=0.001$ and 0.01, respectively. We have looked up numerical simulation works in the literature and found that the values of $s$ and $p$ are all similar, no matter a magnetic field is included or not and what kind of initial conditions are adopted. The density profile we obtain is in good quantitative agreement with that obtained from the detailed observations and modeling to Sgr A* and NGC 3115. The origin of such a accretion rate profile will be investigated in a subsequent paper.

The effects of viscosity on the circumplanetary disks residing in the vicinity of protoplanets are investigated through two-dimensional hydrodynamical simulations with the shearing sheet model. We find that viscosity can affect properties of the circumplanetary disk considerably when the mass of the protoplanet is $M_p \lesssim 33M_\oplus$, where $M_\oplus$ is the Earth mass. However, effects of viscosity on the circumplanetary disk are negligibly small when the mass of the protoplanet $M_p \gtrsim 33M_\oplus$. We find that when $M_p \lesssim 33M_\oplus$, viscosity can disrupt the spiral structure of the gas around the planet considerably and make the gas smoothly distributed, which makes the torques exerted on the protoplanet weaker. Thus, viscosity can make the migration speed of a protoplanet lower. After including viscosity, size of the circumplanetary disk can be decreased by a factor of $\gtrsim 20%$. Viscosity helps to transport gas into the circumplanetary disk from the differentially rotating circumstellar disk. The mass of the circumplanetary disk can be increased by a factor of 50% after viscosity is taken into account when $M_p \lesssim 33M_\oplus$. Effects of viscosity on the formation of planets and satellites are briefly discussed.

In a hot, dilute, magnetized accretion flow, the electron mean-free path can be much greater than the Larmor radius, thus thermal conduction is anisotropic and along magnetic field lines. In this case, if the temperature decreases outward, the flow may be subject to a buoyancy instability (the magnetothermal instability, or MTI). The MTI amplifies the magnetic field, and aligns field lines with the radial direction. If the accretion flow is differentially rotating, the magnetorotational instability (MRI) may also be present. Using two-dimensional, time-dependent magnetohydrodynamic simulations, we investigate the interaction between these two instabilities. We use global simulations that span over two orders of magnitude in radius, centered on the region around the Bondi radius where the infall time of gas is longer than the growth time of both the MTI and MRI. Significant amplification of the magnetic field is produced by both instabilities, although we find that the MTI primarily amplifies the radial component, and the MRI primarily the toroidal component, of the field, respectively. Most importantly, we find that if the MTI can amplify the magnetic energy by a factor $F_t$, and the MRI by a factor $F_r$, then when the MTI and MRI are both present, the magnetic energy can be amplified by a factor of $F_t \cdot F_r$. We therefore conclude that amplification of the magnetic energy by the MTI and MRI operates independently. We also find that the MTI contributes to the transport of angular momentum, because radial motions induced by the MTI increase the Maxwell (by amplifying the magnetic field) and Reynolds stresses. Finally, we find that thermal conduction decreases the slope of the radial temperature profile. The increased temperature near the Bondi radius decreases the mass accretion rate.

How many fraction of gas available at the outer boundary can finally fall onto the black hole is an important question. It determines the observational appearance of accretion flows, and is also related with the evolution of black hole mass and spin. Previous two-dimensional hydrodynamical simulations of hot accretion flows find that the flow is convectively unstable because of its inward increase of entropy. As a result, the mass accretion rate decreases inward, i.e., only a small fraction of accretion gas can fall onto the black hole, while the rest circulates in the convective eddies or lost in convective outflows. Radiation is usually neglected in these simulations. In many cases, however, radiative cooling is important. In the regime of the luminous hot accretion flow (LHAF), radiative cooling is even stronger than the viscous dissipation. In the one dimensional case, this implies that the inward increase of entropy will become slower or the entropy even decreases inward in the case of an LHAF. We therefore expect that convective instability becomes weaker or completely disappears when radiative cooling is important. To examine the validity of this expectation, in this paper we perform two-dimensional hydrodynamical simulations of hot accretion flows with strong radiative cooling. We find that compared to the case of negligible radiation, convection only becomes slightly weaker. Even an LHAF is still strongly convectively unstable, its radial profile of accretion rate correspondingly changes little. We find the reason is that the entropy still increases inward in the two-dimensional case.

How many fraction of gas available at the outer boundary can finally fall onto the black hole is an important question. It determines the observational appearance of accretion flows, and is also related with the evolution of black hole mass and spin. Previous two-dimensional hydrodynamical simulations of hot accretion flows find that the flow is convectively unstable because of its inward increase of entropy. As a result, the mass accretion rate decreases inward, i.e., only a small fraction of accretion gas can fall onto the black hole, while the rest circulates in the convective eddies or lost in convective outflows. Radiation is usually neglected in these simulations. In many cases, however, radiative cooling is important. In the regime of the luminous hot accretion flow (LHAF), radiative cooling is even stronger than the viscous dissipation. In the one dimensional case, this implies that the inward increase of entropy will become slower or the entropy even decreases inward in the case of an LHAF. We therefore expect that convective instability becomes weaker or completely disappears when radiative cooling is important. To examine the validity of this expectation, in this paper we perform two-dimensional hydrodynamical simulations of hot accretion flows with strong radiative cooling. We find that compared to the case of negligible radiation, convection only becomes slightly weaker. Even an LHAF is still strongly convectively unstable, its radial profile of accretion rate correspondingly changes little. We find the reason is that the entropy still increases inward in the two-dimensional case.

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.