Issue |
A&A
Volume 654, October 2021
|
|
---|---|---|
Article Number | A63 | |
Number of page(s) | 11 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/202141047 | |
Published online | 11 October 2021 |
Outflows in the presence of cosmic rays and waves with cooling
1
Graduate Institute of Astronomy, National Central University, Zhongli Dist.,
Taoyuan City
32001,
Taiwan, ROC
e-mail: bilal@gm.astro.ncu.edu.tw; cmko@astro.ncu.edu.tw
2
Department of physics, National Central University, Zhongli Dist.,
Taoyuan City
32001,
Taiwan, ROC
3
I. E. Tamm Theoretical Physics Division of P. N. Lebedev Institute of Physics,
Leninskii pr. 53,
119991
Moscow,
Russia
4
Moscow Institute of Physics and Technology (State University),
9, Institutsky lane,
Dolgoprudny
141707,
Russia
Received:
10
April
2021
Accepted:
17
August
2021
Context. Plasma outflow from a gravitational potential well with cosmic rays and self-excited Alfvén waves with cooling and wave damping is studied in the hydrodynamics regime.
Aims. We study outflows in the presence of cosmic ray and Alfvén waves including the effect of cooling and wave damping. We seek physically allowable steady-state subsonic-supersonic transonic solutions.
Methods. We adopted a multi-fluid hydrodynamical model for the cosmic ray plasma system. Thermal plasma, cosmic rays, and self-excited Alfvén waves are treated as fluids. Interactions such as cosmic-ray streaming instability, cooling, and wave damping were fully taken into account. We considered one-dimensional geometry and explored steady-state solutions. The model is reduced to a set of ordinary differential equations, which we solved for subsonic-supersonic transonic solutions with given boundary conditions at the base of the gravitational potential well.
Results. We find that physically allowable subsonic-supersonic transonic solutions exist for a wide range of parameters. We studied the three-fluid system (considering only forward-propagating Alfvén waves) in detail. We examined the cases with and without cosmic ray diffusion separately. Comparisons of solutions with and without cooling and with and without wave damping for the same set of boundary conditions (on density, pressures of thermal gas, cosmic rays and waves) are presented. We also present the interesting case of a four-fluid system (both forward- and backward-propagating Alfvén waves are included), highlighting the intriguing relation between different components.
Key words: hydrodynamics / waves / plasmas / cosmic rays / ISM: magnetic fields / stars: winds, outflows
© ESO 2021
1 Introduction
Cosmic ray energy density and the energy density of other components of the interstellar medium (e.g., different phases of gas, magnetic field) are on the same order of magnitude (e.g., Ferrière 2001; Cox 2005). The equipartition of energy between the different components indicates there are significant interactions among these components of the interstellar medium. Cosmic rays can play a dynamical role in the structure and evolution of the interstellar as well as intergalactic medium. For instances, from the hydrodynamic perspective, cosmic rays can affect instabilities, such as, Parker instability, Jeans instability, magnetorotional instability (e.g., Parker 1969, 1966; Hanasz & Lesch 2000, 2003; Ryu et al. 2003; Kuwabara et al. 2004; Kuwabara & Ko 2006, 2015, 2020; Ko & Lo 2009; Lo et al. 2011; Heintz & Zweibel 2018; Heintz et al. 2020), and they can modify structures and outflows (e.g., Ko et al. 1991; Yang et al. 2012; Girichidis et al. 2016; Dorfi & Breitschwerdt 2012; Recchia et al. 2016; Ruszkowski et al. 2017; Mao & Ostriker 2018; Farber et al. 2018; Holguin et al. 2019; Dorfi et al. 2019; Yu et al. 2020; Recchia 2020; Ramzan et al. 2020).
The cosmic-ray driven wind mechanism was first discussed by Ipavich (1975). Subsequently, Breitschwerdt et al. (1991, 1993), and Zirakashvili et al. (1996), laid down the proper formulation of cosmic ray propagation with diffusion and wave damping for cosmic-ray driven winds. Later studies (e.g., Everett et al. 2008; Recchia et al. 2016; Wiener et al. 2017; Ruszkowski et al. 2017; Farber et al. 2018) showed that the continuous gas outflows from galaxies can affect their evolution directly. For more details, we refer to the works of Zhang (2018) and Recchia (2020).
Steady-state thermal wind transonic solutions with cooling for a wide range of mass-loading factors, energy-loading factors, galaxy mass, and galaxy radius was studied by Bustard et al. (2016). These authors showed that radiative losses or cooling can enormously affect the wind solutions. The cooling of the gas is important for finding smooth transition from the subsonic to supersonic branch (transonic solutions). Bustard et al. (2016) presented the cases of thermal driven wind with cooling and provided feasible transonic solutions. In this article, we investigate the effect of cooling on outflows under the dynamical effect of cosmic rays and waves. We explore different possibilities of subsonic-supersonic transonic solutions when cosmic ray diffusion and wave damping mechanism are considered. We note that the terms “subsonic” and “supersonic” in this article do not refer to the propagation speed of a disturbance in the medium, but rather the speed that is determined by the critical point condition of the system. We call this speed the “effective sound speed”.
The paper is organized as follows. In Sect. 2, we present the four-fluid model of the cosmic-ray-plasma system in flux-tube formation. The four fluids are thermal plasma, cosmic ray, forward- and backward-propagating Alfvén waves. A simpler system, namely, a system with only a forward-propagating wave (i.e., a three-fluid system or a one-wave system) is studied in Sect. 4 and outflows with and without cosmic ray diffusion are discussed in detail. Section 5 describes an example of the four-fluid system. In Sect. 6, we present our remarks and a summary.
2 Hydrodynamical model for a cosmic ray-plasma system
Cosmic rays interact with magnetised plasma via the embedded magnetic field, in particular, magnetic irregularities or waves. As a result, cosmic rays advect with and get diffused through the plasma. As cosmic rays stream through the plasma, they excite hydromagnetic waves via streaming instabilities. When the energy density of cosmic rays (and waves) are on the same order as that of the thermal plasma, their feedback onto the thermal plasma cannot be neglected. Energy exchange is facilitated by various pressure gradients and wave damping. To study the dynamics of the cosmic ray plasma system, we adopted the four-fluid hydrodynamical model (or two-wave system) put forward by Ko (1992) (see also Zweibel 2017). In this model, the thermal plasma, cosmic rays, and two Alfvén waves (one propagating forward along the magnetic field and one propagating backward) are treated as fluids. In the flux-tube formulation (i.e., considering the dynamics along the magnetic field only) and for the steady state, the governing equations are (Ramzan et al. 2020): (1) (2) (3) (4) (5) (6)
where the energy fluxes are (7) (8) (9)
Equations (1) and (2) come from the divergent free magnetic field and mass continuity equation, ψB and ψm are called the magnetic flux and the mass flow rate, respectively; Δ(ξ) is the cross-sectional area of the magnetic flux tube and ξ is the coordinate along the magnetic field. In Eq. (3), ρ and U are the density and flow velocity of the plasma, Pg, Pc, and are the pressures of the thermal plasma, cosmic rays, and waves (± denote forward- and backward-propagating Alfvén waves). Also, ge is the external gravitational field. Equation (4) describes the influence of gas by cooling, Γ, and by heating via wave damping, . The terms ∓ e±VAdPc∕dξ in Eqs. (5) and (6) are wave excitation by cosmic ray streaming instability. In Eq. (8), κ is the diffusion coefficient of cosmic rays and Dc = κdEc∕dξ is the cosmic ray diffusive flux, Pc∕τ represents stochastic acceleration. Also, is the Alfvén speed. We adopted the following working model for e±, κ and 1∕τ (e.g., Skilling 1975; Ko 1992):
where ν± are the collision frequencies of cosmic rays by forward- and backward-propagating waves. We take in the calculations in Sects. 3 and 4. The smaller the η, the stronger the coupling strength.
3 Wind equation in constant flux-tube geometry
In the following sections, we focus on outflows (especially transonic outflows) under the influence of gas cooling, wave damping, and cosmic ray diffusion as well. To avoid complications arising from a divergent flux-tube, we restrict our discussions to constant flux-tube geometry only, that is, we take Δ to be a constant. To study outflows, it is useful to express the momentum equation (Eq. (3)) in the form of a wind equation (cf. the classic stellar wind in Parker 1958), (11)
where is the Alfvén Mach number and Meff = U∕aeff is the “effective Mach number.” dPc∕dξ appears explicitly only if cosmic ray diffusion coefficient is non-zero. The case without cosmic ray diffusion is discussed in Sect. 4.1. Here the “effective sound speed” is: (12)
We note that the “effective sound speed” is not necessarily the speed of the propagation of a disturbance in the medium. It is determined by the “critical point” of the system, as seen in the term on the left-hand side of Eq. (11); more details are given in Appendix A.
Equations for the gas, cosmic rays, and waves are: (14) (15) (16)
In the following discussion, we take the galactic wind as our fiducial example of outflows. As a working model, we adopted non-linear Landau damping as the wave damping mechanism, , and bremsstrahlung radiation for cooling mechanism, . The corresponding proportionality constants are denoted by δ± and δc in the following calculations. The gravitational field is (cf. a Kuzmin-like field near the center of the disk), where a is the characteristic scale of the field.
4 Transonic outflows in a three-fluid system
In this section, we consider a simpler system, namely, a three-fluid or one-wave system. In this system, only the forward-propagating wave is considered, namely, e+ = 1, e− =0, 1∕τ = 0. We seek physically allowable solutions from a given set of boundary conditions at the base of the potential well. Physically allowable solution means the outflow quantities are single-value with respect to ξ, and the velocity and pressures are non-negative and finite. In particular, we are interested in subsonic-supersonic transonic solutions. Here, we discuss two different cases: (1) systems without cosmic ray diffusion and (2) systems with cosmic-ray diffusion. For subsonic or supersonic solutions, these two cases are not expected be too different (Ramzan et al. 2020). However, for transonic solutions, they can be very different because the effective sound speeds for the two cases are very different (cf. Eqs. (12) and (18)).
4.1 Three-fluid outflow without cosmic ray diffusion
In this subsection, we study three-fluid outflows by neglecting cosmic ray diffusion (i.e., cosmic ray is strongly coupled to thermal gas). Setting e+ = 1, e− = 0, 1∕τ = 0 and κ = 0 in Eq. (16), the cosmic-ray pressure can be expressed analytically as Eq. (A.1) (or see Ramzan et al. 2020) and the three-fluid wind equation becomes (cf. Eq. (11)): (17)
where the “effective sound speed” is modified to (cf. Eq. (12)): (18)
and . We explore outflow against gravity, thus, ge is negative (and approaches zero at large distances away from the gravitational potential well). For convenience, we call the right hand side of the wind equation (for this subsection it is Eq. (17)) the driving term. If there is no cooling and no wave damping, then the driving term is always negative. Thus, in the subsonic regime (Meff < 1), d U∕d ξ > 0 and the outflow speed is monotonically increasing with ξ. We may call this accelerating outflow. There are two possibilities: (1) the accelerating outflow is always subsonic and approaches an asymptotic value at large distances or (2) the flow speed of the solution approaches the effective sound speed at finite distance, then dU∕dξ →∞ and the solution becomes unphysical. By the same token, in the supersonic regime (Meff > 1), d U∕d ξ < 0, and the outflow is a decelerating outflow. Here, there are also two possibilities: (1) a decelerating supersonic outflow that approaches a finite velocity at large distances or (2) the solution approaches the effective sound speed and becomes unphysical.
If cooling is included (but still ignoring wave damping), then the driving term can be positive or negative. Thus, a new possibility arises. It is then possible to have a solution that as its flow speed approaches the effective sound speed and driving term also approaches zero simultaneously. dU∕dξ can be finite at this point (the sonic point or critical point), and the solution can pass through the critical point. A transonic solution is born. A more detailed analysis of this possibility can be found in Appendix A.1.
When wave damping is included, its contribution to the driving term can be positive or negative. For sub-Alfvénic flow the contribution is always negative if γg > 5∕4, while for super-Alfvénic flow the contribution is always negative if γg > 3∕2 and always positive if γg < 5∕4. In this work we only consider cases with γg = 5∕3. We seek physically allowable solutions to the wind equation and pressure equations (Eqs. (14)–(17)) subjected to a set of boundary conditions of at the base of the gravitational potential well. However, not every set of boundary conditions gives us a physically allowable solution. Physically allowable solutions are divided into subsonic solutions, supersonic solutions, and transonic solutions. For the purposes of this study, we are interested in subsonic-supersonic transonic solutions only. At the outset we chose a set at the base of the potential well and then adjusted U to obtain the transonic solution.
In Fig. 1, we show a typical result of three-fluid outflows with cooling. We present different cases side by side. The middle column is the case without wave damping and the right column is the case with wave damping. The boundary conditions at the base for both cases are the same, but Ub are differentin order to obtain a transonic solution for the two cases. As shown in the figure, the velocity andpressure profiles of the two cases are quite similar. We also included a pure thermal outflow for comparison (the left column). The boundary condition for the pure thermal outflow has the same (Pgb, ρb, Bb) as the three-fluid cases, while Ub is adjusted to get the transonic solution. The velocity of the pure thermal outflow is significantly smaller than the three-fluid cases (black solid line in the upper row of Fig. 1 is flow speed). This is expected as the corresponding effective sound speeds are different (thin brown line in the upper row of Fig. 1 is effective sound speed). The effective sound speed of the three-fluid cases has extra contributions from waves and cosmic rays (see Eq. (18)). Nevertheless, the thermal gas pressureprofiles of the three cases are shown to be rather similar (see the lower row of Fig. 1).
Fig. 1 Transonic solutions of a three-fluid system without cosmic ray diffusion. The transonic solution is subsonic-supersonic transition. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1. Left column:pure thermal outflow with cooling (δc = 0.1). Middle column: three-fluid outflow with cooling (δc = 0.1) but without wave damping (δ+ = 0). Right column: three-fluid outflow with cooling (δc = 0.1) and wave damping (δ+ = 0.1). Upper row: black line and thin brown lines show the flow speed and the effective sound speed of the system. For the three-fluid cases, Alfvén speed is also shown (thin gray lines). Lower row: blue, orange, and green lines show the pressures of the thermal gas, cosmic ray, and forward-propagating wave, respectively. The boundary conditions (Pgb, ρb, Bb) are the same for all three cases, and are the same for the two three-fluid cases. In this example, Pgb = 1, Pcb = 0.5, , ρb = 1, . The boundary condition Ub is adjusted to obtain the transonic solutions. |
4.2 Three-fluid outflow with cosmic ray diffusion
An important feature of cosmic ray transport is the diffusion of cosmic rays through the thermal plasma. The diffusion coefficient indicates the coupling of cosmic rays with the plasma; the larger the diffusion coefficient, the weaker the coupling. We modeled the diffusioncoefficient as inversely proportional to the wave pressure (thus, as the wave pressure tends to zero, cosmic rays tend to be decoupled from the plasma). The wind equation is Eq. (11) with only the forward-propagating wave, and e+ =1, e− =0, 1∕τ = 0. We note that cosmic ray pressure gradient appears explicitly in the driving term of the wind equation (the right hand side of Eq. (11)). In comparison with the cases without cosmic ray diffusion (Sect. 4.1), the presence of cosmic ray diffusion shifts the narrative of the existence of transonic solution in two ways: the cosmic ray contribution to the effective sound speed disappears (see Eq. (12)), and the role of cosmic rays in the driving term. In fact, transonic solution can exist in three-fluid system with diffusion even without cooling and wave damping (see Appendix A.2). However, we have opted to focus on the role of cooling on transonic solution in the present study.
We explore solutions for the wind equation (Eq. (11)), together with the corresponding pressure equations, namely, Eqs. (14)–(16) that are subject to boundary conditions at the base of the gravitational potential well. In addition to , we need onemore boundary condition for the gradient of cosmic-ray pressure (or cosmic-ray diffusive flux Db) at the base. In Fig. 2, we compare the case with both cooling and wave damping are “OFF” to cases where wave damping is “ON” or cooling is “ON,” or both are “ON.” In this example, wave damping is “OFF” and “ON,” referring to δ+ = 0 and 0.1 and cooling is “OFF” and “ON,” referring to δc = 0 and 0.1. All cases have the same boundary conditions except Dcb. For those for which cooling is “ON,” we adjust Dcb to obtain the transonic solution (subsonic-supersonic transition), while for those for which cooling is “OFF,” we again adjust Dcb to obtain the physically allowable solution (i.e., subsonic solution that extends to large distances).
In Fig. 2, the upper row shows the outflow speed (black), effective sound speed (thin brown) and Alfvén speed (thin gray). The lower row shows the thermal pressure (blue), cosmic ray pressure (orange) and wave pressure (green). The case where both wave damping and cooling are “OFF” is compared to the case that only wave damping is “ON” in the left column, as well as to the case that only cooling is “ON” in the middle column, and to the case that both wave damping and cooling are “ON” in the right column.
In the left column, all the cases are without cooling and they all belong to the subsonic outflow. The dotted line is the case without wave damping and cooling. The dashed, dot-dashed, and solid lines belong to wave damping case and all with the same δ+ = 0.1, but with slightly different Dcb. Here, Dcb decreases from dashed to dot-dashed to solid and the solid line has the minimum Dcb beyond which there is no physically allowable solution in this example. The solid line is barely different from the dotted line. The dashed and dot-dashed lines have a rapid increase in U and Pc (and a rapid decrease in Pg) somewhere along the flow where becomes diminish and cosmic rays decouple from the thermal gas. These are typical quasi-thermal outflows discussed in Ramzan et al. (2020). Appendix B gives a summary of different flow profiles in three-fluid system (with forward-propagating wave).
In the middle and the right columns, cooling is included and the solutions are subsonic-supersonic transonic flows (solid lines). The profiles carry the flavor of quasi-thermal outflows. It appears that the transition (sonic point) takes place at a further distance for the case without wave damping (the middle column). However, the profiles at large distances for both cooling cases are quite similar in this example.
We are also interested in the effect of the cosmic-ray coupling strength on the results with cooling (in which the physically allowable solutions are transonic solutions). The cosmic ray diffusion coefficient in the three-fluid case is . The coupling is stronger for smaller η. In Fig. 2, η = 1. For the cases with cooling in Fig. 2, we work out the corresponding solutions for η = 5, 2, 1, and 0.5 (by adjusting Dcb), and the result are shown in Fig. 3. In the figure, the left column is the case with cooling “ON” and wave damping “OFF” (δc = 0.1, δ+ = 0) and the right column with both cooling and wave damping “ON” (δc = 0.1, δ+ = 0.1). The rapid increase in the flow speed and cosmic ray pressure (and decrease in the thermal pressure) become steeper when η is smaller. The solutions will be unphysical for even smaller η. Moreover, we note that for η = 0 (strong coupling limit), the system belongs to Sect. 4.1. The corresponding solution with the same boundary conditions as those in Fig. 2 is unphysical.
Fig. 2 Comparison of transonic and subsonic solutions of a three-fluid system with cosmic ray diffusion. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, η = 1. Cases with wave damping and/or cooling are compared to the cases without both wave damping and cooling. All cases have the same boundary conditions at the base except Dcb. The diffusive flux (or cosmic ray pressure gradient) at the base is adjusted to get physically allowable solutions (subsonic solutions for the cases without cooling and transonic solution for the cases with cooling). In this example, Ub = 1, Pgb = 2.5, Pcb = 1, , ρb = 2∕3, . Dotted line displays the case in which both wave damping and cooling are “OFF” (δc = 0, δ+ = 0). Left column: dashed, dot-dashed, and solid lines are the cases in which only wave damping is “ON” (δc = 0, δ+ = 0.1). The three solutions correspond to slightly different Dcb (see text). Middle column: solid line is the case with only cooling is “ON” (δc = 0.1, δ+ = 0). Right column: solid line is the case in which both cooling and wave damping are “ON” (δc = 0.1, δ+ = 0.1). Upper row: black, thin brown, and grey lines show the flow speed, the effective sound speed and the Alfvén speed, respectively. Lower row: blue, orange, and green lines show the pressures of the thermal gas, cosmic ray, and forward-propagating wave, respectively. |
Fig. 3 Comparison of results with different cosmic ray coupling strength. The cosmic ray diffusion coefficient is defined as . Here, we compare solutions for different η of the cases with cooling in Fig. 2. Left column: cooling is “ON” and wave damping is “OFF” (δc = 0.1, δ+ = 0). Right column: both cooling and wave damping are “ON” (δc = 0.1, δ+ = 0.1). Four different cases are shown, (1) η = 5 (dotted line), (2) η = 2 (dashed line), (3) η = 1 (dot-dashed line), and (4) η = 0.5 (solid line). All other parameters are the same as Fig. 2, except that Dcb is adjusted to produce transonic solutions. |
5 Example of a four-fluid outflow
The four-fluid model is the most comprehensive hydrodynamic model for cosmic ray-plasma system. It involves both forward- and backward-propagating waves. It is far more complicated than the three-fluid model discussed inSect. 4. When both waves exist, stochastic acceleration cannot be neglected (1∕τ > 0). In addition, it is possible for the cosmic ray streaming instability (represented by the cosmic ray pressure gradient term in Eq. (15)) to have an opposite effect on the two waves. Here, we present an interesting example to illustrate the intriguing relation between different factors. In Fig. 4, we compare the case with cooling to the case without cooling. Both cases have the same wave damping and the same boundary conditions at the base of the gravitational potential well. By adjusting the cosmic ray diffusive flux Dcb at the base, wesolve the wind equation (Eq. (11)) together with the corresponding pressure equations (Eqs. (14)–(16)) for transonic solution for the case with cooling and for subsonic solution that can extend to large distances (i.e., physically allowable solution) for the case without cooling.
The case with cooling has a number of distinct features, as shown in Fig. 4. The velocity profile (black solid line in the left panel of Fig. 4) has an obvious hump. Non-monotonicity is not observed in cases without cooling. The rapid decrease of velocity is related to the (second) increase of cosmic ray pressure (orange solid line in the right panel) and the rise in the pressure of the backward-propagating wave (purple solid line in the right panel). Moreover, the backward-propagating wave survives at large distances because wave damping approaches to zero. Cooling causes the temperature (T ∝Pg∕ρ) to decrease. At this point, we recall the working model for cooling () and wave damping (). For details, see the end of Sect. 2. When the temperature becomes very low, both cooling and wave damping becomes ineffective.
It is worth pointing out that in this particular example, the wave pressures of the case with cooling are practically zero between ξ = 10 and 30. Thus, cosmic ray diffusion coefficient becomes very large and cosmic ray is decoupled from the plasma. However, when the backward-propagating wave rises again (around ξ = 40), cosmic rays are recoupled to the plasma again.
Another interesting result of this particular example is the flow speed at large distances of the subsonic solution is actually larger than the transonic solution (black dashed and solid lines in the left panel of Fig. 4).
Fig. 4 Example of transonic and subsonic solutions of a four-fluid system with both forward- and backward-propagating waves. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, η = 1. Two cases are compared, (1) dashed line: cooling is “OFF” (δc = 0) and (2) solid line: cooling is “ON” (δc = 0.1). Wave damping is “ON” for both cases (δ± = 1). Left panel: black, thin brown, and grey lines show the flow speed, the effective sound speed, and the Alfvén speed, respectively. Right panel: blue, orange, green, and purple lines show the pressures of thermal gas, cosmic ray, forward-propagating wave, and backward-propagating wave, respectively. The boundary conditions are the same for the two cases. In this example, Ub = 2∕3, Pgb = 1, Pcb = 1, , , ρb = 1, . The cosmic ray pressure gradient or diffusive flux at the base is adjusted to get the transonic solution if cooling is “ON” and get a physically allowable subsonic solution if the cooling is “OFF.” |
6 Summary and discussion
In this article, we study the role of cosmic rays and waves on outflows with cooling. Our fiducial outflow example is that of galactic winds. For simplicity, we only considered steady-state outflows in a constant flux tube (i.e., one-dimensional flow). For a pure thermal outflow in subsonic regime, cooling causes the outflow speed to decrease towards small values and zero at finite distances. The density becomes very large at finite distances and it would be difficult to match the density, for instance, in the intergalacticmedium. However, we find that if cosmic rays and waves are present and the cosmic ray diffusion can be neglected, then still it is possible to have an outflow to large distances even in the subsonic regime. Let us examine the case of three-fluid system with no cosmic ray diffusion and no wave damping in further detail. In this case, Pc and can be expressed in terms of U explicitly (see Eqs. (A.1) and (A.2)). To facilitate discussion, let us write down the equations governing the flow speed, thermal pressure, and temperature as: (19) (20) (21)
where Meff = U∕aeff, aeff is given by Eq. (18), , , Γ† = (γg − 1)Γ∕ρU > 0, and (22)
We note that Meff = Mg for pure thermal outflow. In the supersonic regime, βT > βP > 1, while in the subsonic regime βT < βP < 1 (βP can be negative, and βT < 0 if βP < 1∕γg).
In the subsonic regime ():
(a1): if βTΓ† < βPΓ† < Γ† < , then dU∕dξ > 0, dPg∕dξ < 0, dT∕dξ < 0
(a2): if βTΓ† < βPΓ† < < Γ†, then dU∕dξ < 0, dPg∕dξ < 0, dT∕dξ < 0
(a3): if βTΓ† < < βPΓ† < Γ†, then dU∕dξ < 0, dPg∕dξ > 0, dT∕dξ < 0
(a4): if < βTΓ† < βPΓ† < Γ†, then dU∕dξ < 0, dPg∕dξ > 0, dT∕dξ > 0
In the supersonic regime ():
(b1): If Γ† < βPΓ† < βTΓ† < , then dU∕dξ < 0, dPg∕dξ > 0, dT∕dξ > 0
(b2): If Γ† < βPΓ† < < βTΓ†, then dU∕dξ < 0, dPg∕dξ > 0, dT∕dξ < 0
(b3): If Γ† < < βPΓ† < βTΓ†, then dU∕dξ < 0, dPg∕dξ < 0, dT∕dξ < 0
(b4): if < Γ† < βPΓ† < βTΓ†, then dU∕dξ > 0, dPg∕dξ < 0, dT∕dξ < 0
We note that a system can exist in a range (a4) of the subsonic regime only for a small extent of ξ (if at all). In range (a4), as ξ increases, U decreases, and Pg increases, which implies that βP decreases. When βP is less than 1∕γg, the value of βT is negative and the system will leave the current range (a4) and enter range (a3).
In the subsonic regime when cooling is present, the flow speed decreases as ξ increases if the system is in range (a2), (a3) or (a4). To prevent the flow from stalling at finite distance, the cooling has to be turned off before this happens. The system will then stay in range (a1), and the flow speed will gradually increase towards a finite value as ge ≈ 0 at large distances. Suppose the cooling function vanishes as temperature goes to a small value (e.g., the cooling function approaches zero as T approaches zero). When the flow is in the range of (a2) or (a3) of the subsonic regime, it has a tendency to become stalled. In these ranges, the temperature also decreases. When the temperature decreases rapidly enough, the cooling can be turned off before the speed of the flow reaches zero. It is then possible to have a physically allowable solution that extends to large distances. We then consider what the role of cosmic rays and waves is on the flow in the subsonic regime. We thus denote and for the pure thermal case and and for the three-fluid case (see Eq. (22) for the definition of βP and βT). If U and Pg are the same in the two cases, then and . When βT is smaller, the temperature decreases more rapidly. Hence, the presence of cosmic rays and waves helps the temperature to drop faster and turns off the cooling. As a result, it is easier to have a subsonic flow when cosmic rays and waves are included.
Figure 5 illustrates the effect of cosmic rays and waves on outflows, in particular, subsonic flows. Radiative cooling by bremsstrahlung () is adopted here. The left column shows the pure thermal case. The middle and right columns show the three-fluid case (without cosmic ray diffusion and wave damping). The figure shows four types of solutions, subsonic, subsonic-supersonic, supersonic-subsonic, and supersonic. All the given solutions have the same boundary conditions (Pgb, ρb, Bb) for both cases and, in addition, the same in the case of three-fluid systems. Also, Ub is adjusted to get the two transonic solutions. Ub for the subsonic (supersonic) solution is slightly smaller (larger) than the one for the subsonic-supersonic (supersonic-subsonic) transonic solution. As shown in the top-left panel of the figure for the pure thermal outflow, the subsonic (dot-dashed line) and the supersonic-subsonic (long dashed line) solutions are not capable of being extended to large distances and the flow velocity goes to zero at finite distance. The top-middle panel shows that with the help of cosmic rays and waves, the temperature drops to small value and the cooling is off before the flow velocity goes to zero. The flow is then able to enter range (a1) and approaches a finite flow speed at large distances (see middle column of Fig. 5. The sudden change in U (and Pc, ) in the subsonic region occurs when the temperature drops to zero and the cooling function also becomes zero. In this particular case, the flow goes from range (a2) to range (a1).
The orange dot in the top-left and top-middle panels of Fig. 4 is the sonic point of the subsonic-supersonic transonic solution, i.e., the location where the flow changes from subsonic to supersonic. At this point βT = βP = 1, and the flow moves from range (a1) in the subsonic regime to range (b4) in the supersonic regime. The brown dot is the sonic point of the supersonic-subsonic transonic solution, and the flow moves from range (b1) in the supersonic regime to range (a4) in the subsonic regime (and leaves range (a4) soon). The apparent crossing of the two transonic solutions in the figure is only a projection effect.
We conclude that in the case of three-fluid system without cosmic ray diffusion and wave damping, a subsonic outflow is possible provided that the amount of cosmic ray and wave is large enough. We caution that if cosmic ray diffusion is allowed, the line of reasoning above does not apply and the existence of subsonic flow will be in question (we have tried a range of parameters for physically allowable subsonic solutions but without success).
As discussed above and shown in Fig. 5, a three-fluid subsonic outflow is possible. However, the flow speed at large distances is small, that is, the density is large (e.g., the example in Fig. 5 has a density at large distances larger than at the base of the potential well). We are thus more interested in subsonic-supersonic transonic solutions.
In summary, in this article, we present a study of the steady-state subsonic-supersonic transonic outflow of cosmic-ray plasma systems with cooling. We only considered one-dimensional flow. Our fiducial model is based on the galactic wind. We choose bremsstrahlung radiation as working model for cooling and non-linear Landau damping for wave damping.
We examined the three-fluid model or one-wave model (which comprises thermal plasma, cosmic ray, and self-excited forward-propagating Alfvén wave) without cosmic ray diffusion in detail (Sect. 4.1). Here are set at the base of the gravitational potential well. Here, the value of U at the base is adjusted to obtain a transonic solution. With the same pressures, density, and magnetic field at the base, the profiles of the flow speed and pressures of the cases with and without wave damping are rather similar. However, even with the same thermal pressure and density at the base, the flow speed of the pure thermal transonic outflow is significantly lower than the three-fluid flow, which is expected, as the effective sound speeds are different for both cases (see Fig. 1).
When cosmic ray diffusion is taken into consideration in the three-fluid system, the behavior exhibited by the solutions can be quite different, particularly when the cosmic-ray diffusion coefficient depends on waves (Sect. 4.2). In this case, are set at the base and the cosmic ray pressure gradient or diffusive flux (Dcb) is adjusted to get transonic solutions. We compare the case of without both cooling and wave damping to cases of either cooling or wave damping, or both. In this example, cases with either cooling or wave damping show “rapid change” in velocity andpressure profiles. This is a characteristic feature of quasi-thermal outflows (Ramzan et al. 2020; see Appendix B). When cooling is considered, subsonic-supersonic transition is the only physical solution. In the presence of wave damping, the “sonic point” seems to be pulled closer to the base (compare the middle and the right columns of Fig. 2). When the coupling between cosmic rays and the thermal gas increases (η decreases), the “rapid change” in the profiles becomes “steeper” (see Fig. 3). In this example, there will be no physical solution for a small-enough value for η.
In Sect. 5, we gave an example on the most comprehensive cosmic ray-MHD system, the four-fluid model or two-wave model (where both forward- and backward-propagating waves are included). Due to complicated interactions between different components, the profiles of the flow and pressures may show some interesting features. For instance, in this particular example, when cooling is “ON” (wave damping is also “ON”) the flow speed has a hump that is not observed in cases without cooling, Moreover, along the flow cosmic ray is first coupled to the plasma, subsequently decoupled and then recoupled again at the further down the flow (indicated by the amount of waves along the flow, shown as green and purple solid lines in the right panel of Fig. 4). In addition, in this particular example, the flow speed at large distances of the subsonic solution canbe larger than that of the subsonic-supersonic transonic solution (left panel of Fig. 4).
Fig. 5 Comparison of solutions of a three-fluid system with cooling but without cosmic ray diffusion and wave damping. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, δc = 0.1, δ+ = 0. Four solutions are shown: (1) solid line: subsonic-supersonic transonic solution; (2) Long-dashed line: supersonic-subsonic transonic solution; (3) dot-dashed line: Subsonic solution; and (4) dashed line: supersonic solution. The three-fluid case is shown in the middle and right columns, with boundary conditions . For comparison, a corresponding pure thermal case is shown in the left column with boundary conditions (Pgb, ρb, Bb)=(1, 1, 0.25). In the top-left and top-middle panels, the orange and brown dots are the sonic points of the subsonic-supersonic and supersonic-subsonic transonic solutions, respectively. |
Acknowledgements
C.M.K. and B.R. are supported in part by the Taiwan Ministry of Science and Technology grant MOST 109-2112-M-008-005 and MOST 110-2112-M-008-005.
Appendix A Critical behaviors of three-fluid systems
In this appendix, we analyse the critical behaviour (or fixed points) of two three-fluid systems. This is important for the study of transonic solutions (i.e., solutions passing through fixed points or critical points).
Appendix A.1 Three-fluid system with cooling but without cosmic ray diffusion and wave damping
We considered a three-fluid system (one-wave system with forward-propagating wave) without cosmic ray diffusion and wave damping. The cosmic ray pressure and wave pressure can be expressed in terms of flow speed by solving Equations (15) and (16) (see Ramzan et al. (2020)): (A.1) (A.2)
where , . The governing set of equations can be written as a three-dimensional autonomous system, (A.3)
where Meff= U∕aeff, is given by Equation (18). The fixed point or critical point of the system is determined by setting the right hand side of the Equation (A.3) equal to zero. This can be achieved by: (A.7)
which gives a critical line in the three-dimensional (ξ,U,Pg) space, (A.8) (A.9) (A.10)
where μ is the parametric variable of the critical line and Γs = Γs(ρs, Pgs). For illustration, we take and show the result in Figure A.1. The critical line is the long-dashed pink line. Linear analysis show that the points on the critical line are saddle-points except for those close to Pg = 0. In the figure, we show transonic solutions from three points on the critical line. For each point, there are two transonic solutions: one from subsonic to supersonic (solid line) and the other from supersonic to subsonic (dashed line).
Appendix A.2 Three-fluid system with cosmic ray diffusion but without cooling and wave damping
We also considered a three-fluid system (with forward-propagating wave) in which cosmic ray diffusion is included but cooling and wave damping are neglected. The diffusion coefficient is inversely proportional to the wave pressure. In addition to magnetic flux and mass flux, the total energy flux and the wave-action integral are constant in this system: (A.11) (A.12)
, , and Fg, Fc and are given by Equations (7), (8), and (9). Here Ψe is the external gravitational potential and: (A.13)
Solving Equations (A.11) and (A.12) for and , (A.14) (A.15)
The governing set of equations can be written as a set of three-dimensional autonomous systems: (A.16)
where (A.17) (A.18) (A.19) (A.20)
where Meff= U∕aeff, is given by Equation (12) with . The fixed point or critical point of the system is determined by setting the right hand side of Equation (A.16) equal to zero. This can be achieved via: (A.21)
which gives a critical line in three dimensional (ξ,U,Pc) space, (A.22) (A.23) (A.24)
where Ψes = Ψe(ξes), ges = ge(ξes), and ρs, Pgs and are given by Equations (A.13) and (A.14); also, μ is the parametric variable of the critical line.
Fig. A.1 Critical line and three sets of transonic solutions of a three-fluid system without cosmic ray diffusion. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, , , δc = 0.1, Ac = 0.5, . Magenta long-dashed line is the critical line. Solid lines: subsonic-supersonic transitions. Dashed lines: supersonic-subsonic transitions. |
Appendix B Summary of three-fluid outflows
In this appendix, we discuss some general properties of the three-fluid system with a forward-propagating wave and cosmic-ray diffusion. In the realm of physically allowable solutions starting in the subsonic region at the base of potential, the profiles of U, Pc, and Pg can be categorized into two regimes (Ramzan et al. 2020): regime A (called cosmic ray accompanied outflows) and regime B (quasi-thermal outflows). In regime A, the value of U increases (Pc and Pg decreasing) gradually to a constant value at large distances (e.g., Figure 1). Cosmic rays are coupled to the thermal gas at all distances (i.e., remains finite all the way). In regime B, values of U and Pc have a general trend of increasing towards a finite value at large distances (Pg decreases). The profile of U and Pc has a characteristic rapid increase or a “jump” at some distance from the base (Pg has a rapid decrease or a “drop,” e.g., Figure 2). Also, decreases monotonically and has a “drop” at the same location and approaches zero afterward. Cosmic rays are decoupled from the thermal gas after the “jump” and Pc becomes constant afterward. Wave damping and cooling are not necessary conditions for the existence of the two regimes, although they may affect parameters defining the regimes. In regime A, the set of boundary conditions at the base that gives a physically allowable solution is “isolated” in the sense that any slight change in one (and only one) parameter will push the system into the unphysical domain (however, if more than one parameter can be changed, then it is possible to get into some other physically allowable domain). In contrast, in regime B, a small change in one (and only one) parameter of a physically allowable set may also give physical solutions. This can be understood as follows.
To simplify the argument, let us consider a basic case (system without wave damping and cooling). In this case, only the first two terms on the right hand side of Equation (11) remain (with e+ = 1 and e− = 0), and the second and third terms on the right hand side of Equation (16) vanish. Equation (16) can be viewed as an equation for the cosmic ray diffusive flux Dc = κdEc∕dξ (). When Meff < 1 (“subsonic” region), Equation (11) and Equation (16) form a positive feedback loop in the following sense. If both dU∕dξ and Dc are negative (positive), then both will be more negative (positive) downstream. There are two possible routes, route (a) and route (b), for avoiding divergence and ending up with a physically allowable solution.
Let the solution start at the base with Dc < 0 and dU∕dξ > 0. Dc then increases while dU∕dξ decreases. For one particular set of boundary values, Dc becomes less and less negative and approaches to zero at large distances (and dU∕dξ approaches zero from the positive side). This is route (a).
If a solution starts with the same boundary values as route (a) except that Dc is (slightly) smaller, then dU∕dξ will turn negative before Dc goes to zero. The positive feedback loop activates and the solution will go to the unphysical domain (U and Pc become negative). However, if a solution starts with the same boundary values as route (a) except that Dc is larger, the situation will be different.
Supposing the solution starts from Dc < 0 and dU∕dξ > 0, Dc increases and becomes positive at some finite ξ. The positive feedback loop kicks in and both Dc and dU∕dξ increase rapidly. A positive Dc, namely, a positive cosmic-ray pressure gradient, will cause to diminish via the cosmic-ray streaming instability. Cosmic rays will be decoupled from the thermal gas and the outflow behaves like a thermal outflow. If the velocity and thermal pressure at that location belong to the physically allowable region of a pure thermal outflow, then a physically allowable solution is obtained: this is route (b). The set of parameters that allows for route (b) belonging to regime B (quasi-thermal). If the velocity and thermal pressure after the “rapid change” do not belonging to any physically allowable region of a pure thermal outflow, then the set of parameters only allows for route (a) and it belongs to regime A mentioned above.
Once again, the characteristics of a quasi-thermal outflow is a “jump” in velocity and cosmic-ray pressure (and a “drop” in thermal pressure) at the location where diminishes towards zero. A small change in Dc at the base will change the location of the “jump.” Smaller Dc (more negative) pushes the “jump” to larger distances. We point out that there is a minimum Dc such that the “jump” occurs at “infinity” (e.g., the left column of Figure 2). This particular solution coincides with route (a). The solutions become unphysical for value of Dc that are smaller than the minimum.
References
- Breitschwerdt,D., McKenzie, J. F., & Völk, H. J. 1991, A&A, 245, 79 [NASA ADS] [Google Scholar]
- Breitschwerdt,D., McKenzie, J. F., & Völk, H. J. 1993, A&A, 269, 54 [NASA ADS] [Google Scholar]
- Bustard, C., Zweibel, E. G., & D’Onghia, E. 2016, ApJ, 29, 819 [Google Scholar]
- Chevalier, R. A., & Clegg, A. W. 1985, Nature, 317, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Cox, D. P. 2005, ARA&A, 43, 1 [Google Scholar]
- Dorfi, E. A., & Breitschwerdt, D. 2012, A&A, 540, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dorfi, E. A., Steiner, D., Ragossnig, F., & Breitschwerdt, D. 2019, A&A, 630, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Everett, J. E., Zweibel, E. G., Benjamin, R. A., et al. 2008, ApJ, 674, 258 [NASA ADS] [CrossRef] [Google Scholar]
- Farber, R., Ruszkowski, M., Yang, H.-Y. K., & Zweibel, E. G. 2018, ApJ, 856, 112 [NASA ADS] [CrossRef] [Google Scholar]
- Ferrière, K. M. 2001, Rev. Mod. Phys. 73, 1031 [NASA ADS] [CrossRef] [Google Scholar]
- Girichidis, P., Naab, T., Walch, S., et al. 2016, ApJ, 816, L19 [NASA ADS] [CrossRef] [Google Scholar]
- Hanasz, M., & Lesch, H. 2000, ApJ, 543, 235 [NASA ADS] [CrossRef] [Google Scholar]
- Hanasz, M., & Lesch, H. 2003, A&A, 412, 331 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Heintz, E., & Zweibel, E. G. 2018, ApJ, 860, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Heintz, E., Bustard, C., & Zweibel, E. G. 2020, ApJ, 891, 157 [NASA ADS] [CrossRef] [Google Scholar]
- Holguin, F., Ruszkowski, M., Lazarian, A., Farber, R., & Yang, H.-Y. K. 2019, MNRAS, 490, 1271 [NASA ADS] [CrossRef] [Google Scholar]
- Ipavich, F. M. 1975, ApJ, 196, 107 [NASA ADS] [CrossRef] [Google Scholar]
- Ko, C. M. 1992, A&A, 259, 377 [NASA ADS] [Google Scholar]
- Ko, C. M., & Lo, Y. Y. 2009, ApJ, 691, 2 [Google Scholar]
- Ko, C. M., Dougherty, M. K., & McKenzie, J. F. 1991, A&A, 241, 62 [NASA ADS] [Google Scholar]
- Kulsrud, R. M. 2005, Plasma Physics for Astrophysics (Princeton, NJ: Princeton University Press) [CrossRef] [Google Scholar]
- Kuwabara, T., & Ko, C. M. 2006, ApJ, 636, 290 [NASA ADS] [CrossRef] [Google Scholar]
- Kuwabara, T., & Ko, C. M. 2015, ApJ, 798, 79 [Google Scholar]
- Kuwabara, T., & Ko, C. M. 2020, ApJ, 899, 72 [NASA ADS] [CrossRef] [Google Scholar]
- Kuwabara, T., Nakamura, K., & Ko, C. M. 2004, ApJ, 607, 828 [NASA ADS] [CrossRef] [Google Scholar]
- Li, M., Ostriker, J. P., Cen, R., Bryan, G. L., & Naab, T. 2015, ApJ, 814, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Lo, Y. Y, Ko, C. M., & Wang, C. Y. 2011, Comput. Phys. Commun., 182, 177 [NASA ADS] [CrossRef] [Google Scholar]
- Mao, S. A., & Ostriker, E. C. 2018, ApJ, 854, 89 [NASA ADS] [CrossRef] [Google Scholar]
- Parker, E. N. 1958, ApJ, 128, 664 [Google Scholar]
- Parker, E. N. 1963, Interplanetary Dynamical Processes (New York: Interscience Publication) [Google Scholar]
- Parker, E. N. 1966, ApJ, 145, 811 [NASA ADS] [CrossRef] [Google Scholar]
- Parker, E. N. 1969, Space Sci. Rev., 9, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Ramzan, B., Ko, C. M., & Chernyshov, D. O. 2020, ApJ, 117, 905 [Google Scholar]
- Recchia, S. 2020 Int. J Mod. Phys. D 29, 07 [Google Scholar]
- Recchia, S., Blasi, P., & Morlino, G., 2016, MNRAS, 462, 4227 [NASA ADS] [CrossRef] [Google Scholar]
- Ruszkowski, M., Yang, H.-Y. K., & Zweibel, E. 2017, ApJ, 834, 208 [NASA ADS] [CrossRef] [Google Scholar]
- Ryu, D., Jongsoo, K., Seung S. H., & Jones, T. W. 2003, ApJ, 589, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Skilling, J. 1975, MNRAS, 172, 557 [NASA ADS] [CrossRef] [Google Scholar]
- Wiener, J., Pfrommer, C., & Oh, S. P. 2017, MNRAS, 467, 646 [NASA ADS] [Google Scholar]
- Yang, H.-Y. K., Ruszkowski, M., Ricker, P. M., Zweibel, E., & Lee, D. 2012, ApJ, 761, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Yu, B. P. B., Owen, E. R., Wu, K., & Ferreras, I. 2020, MNRAS, 492, 3179 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, D. 2018, Galaxies, 6, 114 [NASA ADS] [CrossRef] [Google Scholar]
- Zirakashvili, V. N., Breitschwerdt, D., Ptuskin, V. S., & Völk, H. J. 1996, A&A, 311, 113 [NASA ADS] [Google Scholar]
- Zweibel, E. G. 2017, Phys. Plasmas, 24, 055402 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Transonic solutions of a three-fluid system without cosmic ray diffusion. The transonic solution is subsonic-supersonic transition. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1. Left column:pure thermal outflow with cooling (δc = 0.1). Middle column: three-fluid outflow with cooling (δc = 0.1) but without wave damping (δ+ = 0). Right column: three-fluid outflow with cooling (δc = 0.1) and wave damping (δ+ = 0.1). Upper row: black line and thin brown lines show the flow speed and the effective sound speed of the system. For the three-fluid cases, Alfvén speed is also shown (thin gray lines). Lower row: blue, orange, and green lines show the pressures of the thermal gas, cosmic ray, and forward-propagating wave, respectively. The boundary conditions (Pgb, ρb, Bb) are the same for all three cases, and are the same for the two three-fluid cases. In this example, Pgb = 1, Pcb = 0.5, , ρb = 1, . The boundary condition Ub is adjusted to obtain the transonic solutions. |
|
In the text |
Fig. 2 Comparison of transonic and subsonic solutions of a three-fluid system with cosmic ray diffusion. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, η = 1. Cases with wave damping and/or cooling are compared to the cases without both wave damping and cooling. All cases have the same boundary conditions at the base except Dcb. The diffusive flux (or cosmic ray pressure gradient) at the base is adjusted to get physically allowable solutions (subsonic solutions for the cases without cooling and transonic solution for the cases with cooling). In this example, Ub = 1, Pgb = 2.5, Pcb = 1, , ρb = 2∕3, . Dotted line displays the case in which both wave damping and cooling are “OFF” (δc = 0, δ+ = 0). Left column: dashed, dot-dashed, and solid lines are the cases in which only wave damping is “ON” (δc = 0, δ+ = 0.1). The three solutions correspond to slightly different Dcb (see text). Middle column: solid line is the case with only cooling is “ON” (δc = 0.1, δ+ = 0). Right column: solid line is the case in which both cooling and wave damping are “ON” (δc = 0.1, δ+ = 0.1). Upper row: black, thin brown, and grey lines show the flow speed, the effective sound speed and the Alfvén speed, respectively. Lower row: blue, orange, and green lines show the pressures of the thermal gas, cosmic ray, and forward-propagating wave, respectively. |
|
In the text |
Fig. 3 Comparison of results with different cosmic ray coupling strength. The cosmic ray diffusion coefficient is defined as . Here, we compare solutions for different η of the cases with cooling in Fig. 2. Left column: cooling is “ON” and wave damping is “OFF” (δc = 0.1, δ+ = 0). Right column: both cooling and wave damping are “ON” (δc = 0.1, δ+ = 0.1). Four different cases are shown, (1) η = 5 (dotted line), (2) η = 2 (dashed line), (3) η = 1 (dot-dashed line), and (4) η = 0.5 (solid line). All other parameters are the same as Fig. 2, except that Dcb is adjusted to produce transonic solutions. |
|
In the text |
Fig. 4 Example of transonic and subsonic solutions of a four-fluid system with both forward- and backward-propagating waves. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, η = 1. Two cases are compared, (1) dashed line: cooling is “OFF” (δc = 0) and (2) solid line: cooling is “ON” (δc = 0.1). Wave damping is “ON” for both cases (δ± = 1). Left panel: black, thin brown, and grey lines show the flow speed, the effective sound speed, and the Alfvén speed, respectively. Right panel: blue, orange, green, and purple lines show the pressures of thermal gas, cosmic ray, forward-propagating wave, and backward-propagating wave, respectively. The boundary conditions are the same for the two cases. In this example, Ub = 2∕3, Pgb = 1, Pcb = 1, , , ρb = 1, . The cosmic ray pressure gradient or diffusive flux at the base is adjusted to get the transonic solution if cooling is “ON” and get a physically allowable subsonic solution if the cooling is “OFF.” |
|
In the text |
Fig. 5 Comparison of solutions of a three-fluid system with cooling but without cosmic ray diffusion and wave damping. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, δc = 0.1, δ+ = 0. Four solutions are shown: (1) solid line: subsonic-supersonic transonic solution; (2) Long-dashed line: supersonic-subsonic transonic solution; (3) dot-dashed line: Subsonic solution; and (4) dashed line: supersonic solution. The three-fluid case is shown in the middle and right columns, with boundary conditions . For comparison, a corresponding pure thermal case is shown in the left column with boundary conditions (Pgb, ρb, Bb)=(1, 1, 0.25). In the top-left and top-middle panels, the orange and brown dots are the sonic points of the subsonic-supersonic and supersonic-subsonic transonic solutions, respectively. |
|
In the text |
Fig. A.1 Critical line and three sets of transonic solutions of a three-fluid system without cosmic ray diffusion. In this figure, γg = 5∕3, γc = 4∕3, GM = 1, a = 1, , , δc = 0.1, Ac = 0.5, . Magenta long-dashed line is the critical line. Solid lines: subsonic-supersonic transitions. Dashed lines: supersonic-subsonic transitions. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.