Self-Similar Solution of Hot Accretion Flow with Anisotropic Pressure

For the accretion flow in extremely low-luminosity active galactic nuclei, such as our Galactic center (Sgr A*) and M 87, the collisional mean-free path of ions may be much larger than its gyroradius. In this case, the pressure parallel to the magnetic field is different from that perpendicular to the field; therefore, the pressure is anisotropic. We study the effects of anisotropic pressure on the dynamics of accretion flow by assuming the flow is radially self-similar. We find that in the case where the outflow is present, the radial and rotational velocities, the sound speed, and the Bernoulli parameter of the accretion flow are all increased when the anisotropic pressure is taken into account. This result suggests that it becomes easier for the accretion flow to generate outflow in the presence of anisotropic pressure.


Introduction
According to the temperature of the gas, black hole accretion disk/flow can be divided into two groups, namely cold accretion disk and hot accretion flow.There are many numerical simulations studying thin disks.One of the hot topics being studied by simulations is jet launching from a thin disk (e.g., Casse and Keppens [1,2]; Zanni et al. [3,4]; Sheikhnezami et al. [5]; Sheikhnezami and Fendt [6,7]).The acceleration of jets from a thin disk is also being studied by analytical works (e.g., Takahashi et al. [8]).Simulations also studied relativistic thin disks (e.g., De Villiers and Hawley [9,10]; Gammie et al. [11]; Noble et al. [12]; Sadowski et al. [13]).Hawley et al. [14] gave a very extensive review for accretion disks.
Observations show that most of the galactic nuclei in the local universe are very dim (Ho [15] and the references therein).Radiatively inefficient hot accretion flows operate in these sources (Narayan and Yi [16]; see the reviews by Ho [15] and Yuan and Narayan [17]).Hot accretion flows also operate in the quiescent and hard states of black hole X-ray binaries (Done et al. [18]; Qiao and Liu [19]; Gilfanov [20]; Zhang et al. [21]; Yan and Yu [22]; see the review by Yuan and Narayan [17]).The temperature of the gas in the hot accretion flow almost equals the virial temperature.In the innermost region, the gas can have a temperature of approximately 10 11 K (Yuan and Narayan [17]).The density is very low and depends on the mass accretion rate.In the accretion flow in the Galactic center, in the innermost region, gas density is approximately 10 −18 g/cm 3 .
There are some extremely low accretion rate accretion flows.For example, the accretion rate of the accretion flow in the Galactic center (Sgr A*) is 2 × 10 −3 L Edd /c 2 = 10 21 g/s, with L Edd and c being the Eddington luminosity and the speed of light, respectively.The other example is the accretion flow in the M87 galaxy.The accretion rate of the accretion flow in M87 is 1.8 × 10 25 g/s (Yuan and Narayan [17]).In these extremely low accretion rate hot accretion flows, such as the accretion flow in the Galactic center (Sgr A*) and M87, the Coulomb collision mean free paths of both electrons and ions are much larger than the typical length-scale of the accretion flows, ∼GM/c 2 , where G and M are the gravitational constant and black hole mass, respectively (Tanaka and Menou [54]; Johnson and Quataert [55]; Foucart et al. [56]).At first glance, magneto-hydrodynamic (MHD) approximations cannot be applied to these plasmas, and we should solve the distribution function of both ions and electrons.However, particle-in-cell simulations show that the wave-particle interactions can effectively increase the collision rate between particles (Kunz, Schekochihin, and Stone [57]; Riquelme, Quataert, and Verscharen [58]; Sironi and Narayan [59]).The plasmas in these extreme low accretion rate accretion flows are just weakly collisional.For the weakly-collisional plasmas, the non-ideal effects can be treated as perturbations relative to an ideal fluid, by the inclusion of anisotropic thermal conduction and anisotropic pressure (Foucart et al. [56]).
In weakly-collisional accretion flows, anisotropic thermal conduction is important and along the magnetic field lines (Parrish and Stone [60]; Sharma, Quataert, and Stone [61]; Bu, Yuan and Stone [62]).Thermal conduction can transport thermal energy from the inner (hotter) to the outer (cooler) regions.Previous works showed that at large radii, the gas temperature can be increased to be above the virial temperature.Strong outflows can be launched, and the black hole mass accretion rate can be significantly decreased (Bu et al. [63]).
As introduced above, in weakly-collisional accretion flows, one should include an anisotropic pressure as perturbations to an ideal fluid.In weakly-collisional accretion flows, the ions collision mean free path can be much larger than its gyroradius.For example, the ions' mean free path in the accretion flow in the Galactic center (Sgr A*) is ∼10 17 cm (Tanaka and Menou [54]).If we assume that magnetic pressure is 0.1-times the gas pressure, we will have that ions gyro-radius is ∼4 ×10 6 cm (Bu et al. [63]).In these systems, pressure parallel to the magnetic field is different from that perpendicular to the magnetic field (Quataert et al. [64]; Sharma et al. [65]; Chandra et al. [66]).The pressure anisotropy can be modeled by an anisotropic viscosity with respect to magnetic field lines (Balbus [67]; Islam and Balbus [68]).In Wu et al. [69], the effects of anisotropic pressure on properties of hot accretion flows were studied.We note that in that paper, only a very weak magnetic field was considered, and the magnetic pressure was at least four orders of magnitude smaller than the gas pressure.They found that when magnetic fields are extremely weak, Maxwell stress is negligibly small.Anisotropic pressure can transfer angular momentum, and a hot accretion flow can be formed.Furthermore, they found very weak outflow.
Currently, it is quite difficult to measure or even to constrain the strength of the magnetic field.Among the existing methods, the Faraday rotation measure (RM) can provide the integration of electron density and magnetic field along the line of sight.RM combined with spectrum can roughly give magnetic field strength.For example, the RM of Sgr A*, together with the radio-up to-millimeter spectrum, constrains the magnetic field strength in the accretion flow to be 20 Gauss at 10 Schwarzschild radii (Yuan et al. [70]).
The magnetic stress for angular momentum transfer is induced by turbulence driven by magnetorotational instability (Balbus and Hawley [71,72]).The stress is B r B φ /(4π) ∼ B 2 , with B r and B φ being the radial and toroidal components of the magnetic field, respectively.The viscous coefficient (α 1 ) is stress divided by gas pressure.Therefore, we have α 1 ∼ B 2 /p, with p and B being the gas pressure and magnetic field, respectively.Observations of some accretion systems indicate that the viscosity coefficient α 1 ∼ 0.1 (see King et al. [73] and the references in that paper).Since α 1 ∼ B 2 /p, such a value of α implies that the magnetic pressure is approximately 1/10-times the gas pressure.The magnetic field is not significantly weak compared to the gas pressure.The effects of anisotropic pressure in hot accretion flows with a relatively strong magnetic field have not been studied previously.
In this work, we study the effects of pressure anisotropy (or anisotropic viscosity) on the properties of hot accretion flows in this case.We set the magnetic pressure to be 0.1-times the gas pressure.
In Section 2, we will describe the basic equations and assumptions.In Section 3, we will present the results.We discuss and summarize our results in Section 4.

Basic Equations and Assumptions
We study time-steady and axisymmetric (∂/∂t = ∂/∂φ = 0) accretion flows in cylindrical coordinates (r, φ, z).The basic equations describing accretion are as follows, In the equations, ρ is the gas density, v is velocity, e is gas specific internal energy, Φ is the gravitational potential, q vis is the heating rate of the gas, and q rad represents the radiative cooling rate.For hot accretion flows, the radiative cooling is not important.For simplicity, we assume that q rad = 0 and f = 1.In this paper, we use the Newtonian potential.d/dt = ∂/∂t + v • ∇.
Numerical simulations found that the magnetic field in an accretion flow/disk can be decomposed into a large-scale ordered component and a turbulent component (e.g., Machida et al. [74]; Hirose et al. [75]; Bai and Stone [76]).In Equation ( 2), we use B to represent the large-scale ordered magnetic field component.Simulations showed that the main body of the accretion flow is governed by a toroidal magnetic field (see Figure 6 in Hirose et al. [75]).In this paper, we assume that the large-scale ordered magnetic field has only the toroidal component (B φ ).Following Lovelace et al. [77], we assume that B φ is an odd function of z.
H is the half-thickness of the accretion flow.The large-scale magnetic field only has a toroidal component, and we assume axisymmetry; therefore, the divergence-free condition of the magnetic field is satisfied.We note that for the strongly-magnetized accretion flow (MAD disk, Narayan et al. [25]), a substantial poloidal field is present.We plan to study accretion flow with the poloidal magnetic field in the future.
In Equation ( 2), the term ∇ • T represents the angular momentum transfer by the turbulent magnetic field component.For the stress T, we use the usual α description, In this equation, ν 1 = α 1 c 2 s /Ω K .c s is sound speed; Ω K is the Keplerian angular velocity; Ω is the angular velocity of the gas.Numerical simulations showed that the dissipation heating mainly comes from the magnetic reconnection associated with the turbulent magnetic field (Hirose et al. [78]).In Equation ( 4), we use q vis,1 = −TdΩ/dr to represent the heating due to the magnetic reconnection due to the turbulent component of the magnetic field.
Following Chandra et al. [66], we use anisotropic viscous stress tensor Π to model the effects of ion pressure anisotropy (see also Braginskii [79]), where b is the unit vector in magnetic field direction, and I is the unit tensor.bb is a dyadic tensor defined as, and bb : ∇v In Equation ( 8), We use α 2 to denote the strength of anisotropic pressure.In this paper, the magnetic field only has a toroidal component.Correspondingly, the anisotropic pressure only has three non-zero components, In Equation ( 4), the term q vis,2 = −Π : ∇v represents the heating by the work done by anisotropic pressure.

Vertical Integrated Equations
In this paper, we use several parameters (ξ 1 , ξ 2 , ξ 3 ; see the definitions below) to characterize the properties of outflows.We assume that the z component of gas velocity is zero.We also assume that the gas is vertically isothermal.The temperature or sound speed is only a function of r.At a fixed radius, the gas temperature or sound speed does not change with height.We vertically integrate Equations ( 1)-( 3) to obtain one-dimensional equations.In the presence of outflow, the vertically-integrated mass conservation Equation (1) will be as follows, ṁ is constant.The parameter s is in the range −1/2 ≤ s ≤ 1/2 (see the explanation below).
In the presence of outflow, the equations describing inflow are given in Xie and Yuan [80].The radial momentum equation is, In this equation, Σ is the surface density and defined as Σ = 2ρH.In Equation ( 13), v r,w is the radial velocity of outflow.When the radial velocity of outflow is different from that of inflow, radial momentum will be taken away from or deposited into inflow by the outflow.Correspondingly, the inflow velocity will be changed.The first, second, and third terms on the right-hand side of Equation ( 13) are centrifugal, gravitational, and gas pressure gradient forces, respectively.The fourth term is the Lorentz force due to the large-scale ordered magnetic field.The last two terms are the forces of divergence of anisotropic pressure (∇ • Π in Equation ( 2)).
The azimuthal momentum equation is, In this equation, v K = GM r is the Keplerian rotational velocity.v φ,w is the rotational velocity of outflow.If the rotational velocity of outflow is different from that of inflow, angular momentum will be taken away from or deposited into the inflow by the outflow.Correspondingly, the angular momentum of inflow will be changed.The term on the right-hand side is the angular momentum transfer by the turbulent magnetic field (∇ • T in Equation ( 2)).
We define a parameter β to describe the strength of the large-scale ordered magnetic field, The gas in the vertical direction achieves static equilibrium, and the z−component of gravitational force is balanced by gas pressure gradient and magnetic pressure gradient forces, We note that the force of divergence of anisotropic pressure (∇ • Π in Equation ( 2)) only exists in the radial momentum equation.This is because the anisotropic pressure only has the Π rr , Π zz , Π φφ components.Anisotropic pressure cannot transfer angular momentum, because its Π rφ component is zero.
The energy equation is, γ is the ratio of the specific heats, which varies from the relativistic value 4/3 to the typical monatomic gas value 5/3.In this paper, we set γ = 4/3.We find that the results are not sensitive to the value of γ. e w is the specific internal energy of outflow.If the specific internal energy of outflow is different from that of inflow, internal energy will be taken away from or deposited into the inflow by the outflow.
Correspondingly, the internal energy of inflow will be changed.The first term on the right-hand side is the heating due to the reconnection of turbulent magnetic field (q vis,1 in Equation ( 4)).The second term on the right-hand side is the heating due to the work done by anisotropic pressure (q vis,2 in Equation ( 4)).
We parameterize the outflow properties.We assume that the velocities of outflow are proportional to those of inflow, with v r,w = ξ 1 v r and v φ,w = ξ 2 v φ .The specific internal energy of outflow is proportional to that of inflow, with e w = ξ 3 e.

Self-Similar Assumptions
Following Narayan and Yi [16], we set the following radially self-similar assumptions, The radial velocity v r ∝ r −1/2 .The mass accretion rate Ṁ ∝ r 2 ρv r ∝ ρr 3/2 .If outflow is absent, the mass accretion rate will be a constant with radius, and we will have ρ ∝ r −3/2 .When outflow is present, the mass accretion rate will decrease with decreasing radius.Numerical simulations found that when outflow is present, the density profile is very shallow ρ ∝ r −Γ , with 1/2 ≤ Γ ≤ 3/2 (Stone et al. [81]).The surface density Σ = 2ρH.For hot accretion flows, the scale height of the flow is roughly equal to the radius H ∼ r.Therefore, we have Σ = Σ 0 r −Γ+1 = Σ 0 r s .Σ 0 is a constant.The parameter s = −Γ + 1 is in the range −1/2 ≤ s ≤ 1/2.s = −1/2 corresponds to the case where outflow is absent.

Results
In this paper, we set β = 0.1 and α 1 = 0.1.We focused on studying the effects of anisotropic pressure.For a hot accretion flow, the gas temperature was too high such that gas was fully ionized.Observationally, it is very hard to detect outflow directly through the absorption line.There is only some indirect evidence that outflow should be present for a hot accretion flow (Wang et al. [82]; Ma et al. [44]; Park et al. [45]).The properties (velocity, temperature) of outflows cannot be given by observations.Numerical simulations give us the properties of outflows.The values of ξ were mainly from numerical simulations' results (Yuan et al. [23]).

Accretion Flow without Outflow
In this subsection, we focus on the case where outflows are absent.In this case, the parameter s = −1/2.We show the results in Figure 1.The top-left panel plots radial infall velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The top-right panel plots rotational velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The lower panel plots sound speed (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The viscous coefficient α 1 was set to be 0.1.The ratio of magnetic to gas pressure β was set to be 0.1.For the self-similar solution, at any radii, the ratio of velocities to Keplerian velocity was constant.In the figure, we plot the ratio of velocities to Keplerian velocity; therefore, the results (velocities plotted in the top-left panel of Figure 1) can be applied at any radii. .The lower panel plots sound speed (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The viscous coefficient α 1 was set to be 0.1.The ratio of magnetic to gas pressure β was set to be 0.1.
The last term on the right-hand side of the radial momentum Equation ( 19) is the force of divergence of anisotropic pressure.It is clear, when s = −1/2, that the force is negative; this force points towards the central black hole.Therefore, from the top-left panel of Figure 1, we see that with the increase of strength of anisotropic pressure, the radial infall velocity of the accretion flow increases.If outflow is absent, viscosity due to turbulent magnetic field is the only mechanism for angular momentum transfer (see Equations ( 14) and ( 20)).The inward advection of angular momentum by inflowing gas is balanced by the outward angular momentum transportation by viscosity.With the increase of radial infall velocity, the inward advection of angular momentum flux increases (see Equation ( 20)).Correspondingly, the outward angular momentum flux also needs to increase.The increase of outward angular momentum flux can be realized by increasing the sound speed (see Equation ( 20)), because the viscosity T is proportional to c 2 s .Therefore, with the increase of strength of anisotropic pressure, the sound speed (or gas temperature) becomes bigger.From the energy Equation ( 22), we see that the viscous heating (first term on the right-hand side) and heating due to work done by anisotropic pressure (second term on the right-hand side) are balanced by the advection of energy by the inflow (first term on the left-hand side).With the increase of c 1 , heating due to work done by anisotropic pressure increases faster than the increase of energy advection.Therefore, in order to keep the accretion flow stable, the viscous heating term needs to decrease with an increase of c 1 .Therefore, with an increase of α 2 , the rotational velocity decreases.
From Figure 1, we see that with the increase of α 2 , both the radial infall velocity and sound speed increase.Comparing the top-left and lower panels, we see that the infall velocity is always smaller than the sound speed by one order of magnitude at any given value of α 2 .Therefore, the infall velocity is sub-sonic.We assume that magnetic pressure is 10-times smaller than gas pressure.Correspondingly, the Alfvenic velocity is smaller than gas sound speed by a factor of 3.3.Therefore, the gas infall velocity is sub-Alfvenic.

Accretion Flow with Outflow
We first studied the case where the specific velocity and internal energy of outflow are the same as those of inflow (ξ 1 = ξ 2 = ξ 3 = 1).The results are shown by the blue lines in Figure 2.With the increase of strength of anisotropic pressure (α 2 ), the infall velocity (c 1 ) increases.The reason is the same as explained for the case where the outflow is absent.The force of the anisotropic pressure gradient is negative.Therefore, with an increase of α 2 , gas infall velocity increases.With the increase of infall velocity (c 1 ), the sound speed increases.The reason is also the same as in the case where the outflow is absent, as introduced in the last subsection.We do not explain it again here.We now explain why the rotational velocity increases with an increase of α 2 .As introduced above, the viscous heating and heating due to work done by anisotropic pressure are balanced by the advection of energy by the inflow (see Equations ( 17) and ( 22)).With the increase of c 1 , the energy advection term increases faster than the increase of work done by anisotropic pressure.Therefore, in order to keep the accretion flow stable, the viscous heating term needs to increase with an increase of c 1 .Therefore, with an increase of α 2 , the rotational velocity increases.Now, we study the cases where the specific velocity and internal energy of outflow are different from those of inflow.We first study the effects of the discrepancy of the radial velocity between inflow and outflow.We find that for a reasonable change of values of ξ 1 , the effects are of minor importance.This result is consistent with Xie and Yuan [80].Therefore, for the models described below, we set MHD numerical simulation by Yuan et al. [23] has found that the specific angular momentum of outflow is higher than that of inflow.We study the case where the specific angular momentum of outflow is higher than that of inflow by setting ξ 2 = 1.2.The results are shown in Figure 2 by the yellow lines.We first focus on the angular momentum Equation (20).In this equation, as introduced above, the first term on the left-hand side is the angular momentum advection by inflow gas.The second term on the left-hand side is the angular momentum flux taken away by outflow.The term on the right-hand side corresponds to the angular momentum outward transportation by viscosity.The inward advection of angular momentum by gas inflow is balanced by the outward transportation of angular momentum by both outflow and viscosity.In the model with ξ 2 = 1.2, outflow can take away angular momentum, so the outward flux of angular momentum transported by viscosity (denoted by sound speed c 3 in Equation ( 20)) becomes smaller.Therefore, with the same value of α 2 , the sound speed in model ξ 2 = 1.2 is smaller than that in the model ξ 2 = 1.With the help of outflow (ξ 2 = 1.2), the total outward transportation rate of angular momentum by both outflow and viscosity is larger than that in the case ξ 2 = 1.Therefore, the inward advection term of angular momentum (denoted by radial velocity c 1 in Equation ( 20)) in the model with ξ 2 = 1.2 should be larger than that in the model with ξ 2 = 1.Therefore, with the same value of α 2 , the radial velocity in the model with ξ 2 = 1.2 is larger than that in the model with ξ 2 = 1.From the energy Equation ( 22), with an increase of c 1 , the energy advection term increases faster than the increase of work done by anisotropic pressure.Therefore, in order to keep the accretion flow stable, the viscous heating term (denoted by 9/4c 2  2 ) needs to increase with an increase of c 1 .Therefore, with the same value of α 2 , c 2 in the model with ξ 2 = 1.2 is larger than that in the model with ξ 2 = 1. .The lower panel plots sound speed (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The viscous coefficient α 1 was set to be 0.1.The ratio of magnetic to gas pressure β was set to be 0.1.ξ 1 is the ratio of outflow radial velocity to inflow radial velocity.ξ 2 is the ratio of outflow rotational velocity to inflow rotational velocity.ξ 3 is the ratio of outflow sound speed to inflow sound speed.We set ξ 1 = ξ 3 = 1.The blue and yellow lines correspond to ξ 2 = 1 and ξ 2 = 1.2, respectively.
The specific internal energy of outflow was also found to be higher than that of inflow by simulations (e.g., Yuan et al. [23]).In Figure 3, we show the results for the flow in which outflow specific internal energy is higher than that of inflow by the yellow lines.When we calculated the results shown by yellow lines, we set ξ 3 = 1.2.Note that the results shown by blue lines in this figure were calculated with ξ 1 = ξ 2 = ξ 3 = 1.They are the same as the results shown by blue lines in Figure 2. Now, we see the energy Equation ( 22).The first term on the left-hand side is the advection of energy by gas inflow.The second term on the left-hand side is the energy taken away by outflow.The two terms on the right-hand side are heating by viscosity and work done by anisotropic pressure.One portion of the heating energy is advected to the black hole by gas inflow; the other portion is taken away by outflow.In the model with ξ 3 = 1.2, the energy can be taken away by outflow; therefore, the energy advection term (denoted by c 1 [1/(γ − 1) + (s − 1)]) becomes smaller.Therefore, with the same value of α 2 , the radial velocity of inflow (c 1 ) in the model with ξ 3 = 1.2 is smaller than that in the model with ξ 3 = 1.From the angular momentum Equation ( 20), we see that with the decrease of c 1 , the sound speed c 3 also decreases.Therefore, with the same value of α 2 , the sound speed (c 3 ) in the model with ξ 3 = 1.2 is smaller than that in the model with ξ 3 = 1.With the decrease of sound speed, gas pressure gradient force in the radial direction decreases.In order to keep a balance of the forces in the radial direction, the centrifugal force should increase.Therefore, with the same value of α 2 , the value of c 2 in the model with ξ 3 = 1.2 is bigger than that in the model with ξ 3 = 1. .The lower panel plots sound speed (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The viscous coefficient α 1 was set to be 0.1.The ratio of magnetic to gas pressure β was set to be 0.1.ξ 1 is the ratio of outflow radial velocity to inflow radial velocity.ξ 2 is the ratio of outflow rotational velocity to inflow rotational velocity.ξ 3 is the ratio of outflow sound speed to inflow sound speed.We set ξ 1 = ξ 2 = 1.The blue and yellow lines correspond to ξ 3 = 1 and ξ 3 = 1.2, respectively.
The Bernoulli parameter is usually used to judge whether outflow can escape from the black hole gravitational potential to infinity.In the magnetized accretion flow with only a toroidal magnetic field, the Bernoulli parameter is defined as, Here, h = γe/ρ is enthalpy.Both the velocity and gas temperature increase with the increase of the strength of anisotropic pressure (α 2 ).Therefore, the Bernoulli parameter increases with α 2 .The result indicates that the presence of anisotropic pressure makes it easier to generate outflow.

Conclusions
For an extremely low-accretion rate system, the ion's mean-free path can be much larger than its gyroradius, and pressure parallel to the magnetic field is different from that perpendicular to the magnetic field.In this paper, we study the effects of anisotropic pressure on the dynamics of hot accretion flows.We find that when outflow is absent, the radial velocity and sound speed increase with the increase of the strength of anisotropic pressure, the rotational velocity decreases with the increase of the strength of anisotropic pressure.When outflow is present, the radial velocity, rotational velocity, and sound speed increase with the increase of strength of anisotropic pressure.Therefore, in this case, the Bernoulli parameter of the accretion flow increases.This result predicts that when anisotropic pressure is considered, it becomes easier to generate outflow by hot accretion flows.

Figure 1 .
Figure 1.The properties of the accretion flow when outflow is absent (s = −1/2).The top-left panel plots radial infall velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The top-right panel plots rotational velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The lower panel plots sound speed (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The viscous coefficient α 1 was set to be 0.1.The ratio of magnetic to gas pressure β was set to be 0.1.

2 Figure 2 .
Figure 2.The properties of the accretion flow when outflow is present (s = 1/2).The top-left panel plots radial infall velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The top-right panel plots rotational velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The lower panel plots sound speed (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The viscous coefficient α 1 was set to be 0.1.The ratio of magnetic to gas pressure β was set to be 0.1.ξ 1 is the ratio of outflow radial velocity to inflow radial velocity.ξ 2 is the ratio of outflow rotational velocity to inflow rotational velocity.ξ 3 is the ratio of outflow sound speed to inflow sound speed.We set ξ 1 = ξ 3 = 1.The blue and yellow lines correspond to ξ 2 = 1 and ξ 2 = 1.2, respectively.

2 Figure 3 .
Figure 3.The properties of the accretion flow when outflow is present (s = 1/2).The top-left panel plots radial infall velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The top-right panel plots rotational velocity (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The lower panel plots sound speed (in units of Keplerian velocity) as a function of the strength of anisotropic pressure (α 2 ).The viscous coefficient α 1 was set to be 0.1.The ratio of magnetic to gas pressure β was set to be 0.1.ξ 1 is the ratio of outflow radial velocity to inflow radial velocity.ξ 2 is the ratio of outflow rotational velocity to inflow rotational velocity.ξ 3 is the ratio of outflow sound speed to inflow sound speed.We set ξ 1 = ξ 2 = 1.The blue and yellow lines correspond to ξ 3 = 1 and ξ 3 = 1.2, respectively.