Stellar Wind -- Parker Wind Theory

Stellar winds are observed around various type of stars. Early type (massive) stars have large luminosities; the photon absorbed by the bound-bound transition transfers its outward momentum to the gas. This line-driven mechanism seems to work around the early type stars. On the other hand, acceleration mechanism of less massive stars are thought to be related to the coronal activity or dust driven mechanism (dusts absorb the emission and obtain outward momentum from the emission).

Here, we will see the identical mechanism in $\S$2.8.1 and $\S$2.8.2 works to accelerate the wind from a star. Consider a steady state and ignore $\partial /\partial t=0$. The continuity equation (2.1) gives

\begin{displaymath}
r^2 \rho u ={\rm const,}
\end{displaymath} (2.94)

where we used
\begin{displaymath}
{\rm div}\rho {\bf v}=\frac{1}{r^2}\frac{\partial }{\partial r}\left(r^2 \rho u \right).
\end{displaymath} (2.95)

This leads to
\begin{displaymath}
\frac{2}{r}+\frac{1}{r}\frac{\partial \rho}{\partial r}+\frac{1}{u}\frac{\partial u}{\partial r}=0.
\end{displaymath} (2.96)

The equation of motion (2.2) is as follows:
\begin{displaymath}
u \frac{\partial u}{\partial r} = - {{c_s^2 } \over \rho }\frac{\partial \rho }{\partial r} - {{GM_*} \over {r^2 }},
\end{displaymath} (2.97)

where we used $g=-GM_*/r^2$ ($M_*$ means the mass of the central star). From these two equations (2.96) and (2.97) we obtain
\begin{displaymath}
\left( {{\cal M}^2 - 1} \right){1 \over {u}}\frac{\partial u...
...al r} = {2 \over r} - {{GM_*} \over {c_s^2 }}{1 \over {r^2 }},
\end{displaymath} (2.98)

where ${\cal M}\equiv u/c_s$ represents the Mach number of the radial velocity. Take notice that this is similar to equations (2.89) and (2.93). That is, the fact that the rhs of equation (2.98) is positive corresponds to increasing the cross-section $d S/d x>0$. On the contrary, when the rhs is negative, the fluid acts as the cross-section $\S$ is decreasing.

Figure 2.7: Right-hand side of equation (2.98) is plotted against the distance from the center $r$.
\begin{figure}\centering\leavevmode
\epsfxsize =.45\columnwidth \epsfbox{eps/wind-rhs.ps}\end{figure}

For simplicity, we assume the gas is isothermal. The rhs of equation (2.98) varies shown in Figure 2.7. Therefore, near to the star, the flow acts as the cross-section of nozzle is decreasing and far from the star it does as the cross-section is increasing. This is the same situation that the gas flows through the Laval nozzle.

Using a normalized distance $x\equiv r/(GM_*/2c_{is}^2)$, equation (2.98) becomes

\begin{displaymath}
\left( {M^2 - 1} \right){1 \over {\cal M}}{{\partial {\cal M}} \over {\partial x}} = {2 \over x} - {2 \over {x^2 }}.
\end{displaymath} (2.99)

This is rewritten as
\begin{displaymath}
{d \over {dx}}\left( {{{{\cal M}^2 } \over 2} - \log {\cal M} - 2\log x - {2 \over x}} \right) = 0,
\end{displaymath} (2.100)

we obtain the solution of equation (2.99) as
\begin{displaymath}
{\cal M}^2 - 2\log {\cal M} = 4\log x + {4 \over x} + C.
\end{displaymath} (2.101)

This gives how the Mach number ${\cal M}$ varies changing $x$. To explore this, we define two functions:
\begin{displaymath}
f({\cal M}) = 2{\cal M}^2 - 2\log {\cal M}
\end{displaymath} (2.102)

is a function only depending on ${\cal M}$ and
\begin{displaymath}
g(x) = 4\log x + {4 \over x} + C
\end{displaymath} (2.103)

is a function only depending on $x$.

Figure 2.8: $f({\cal M})$ (left) and $g(x)$ for $C=-3$ (right)
\begin{figure}\centering\leavevmode
\epsfxsize =.45\columnwidth \epsfbox{eps/f_M.ps} \hfil
\epsfxsize =.45\columnwidth \epsfbox{eps/g_x.ps}\end{figure}

Since the minima of $f({\cal M})$ and $g(x)$ are respectively $1$ and $C+4$, the permitted region in $(x,{\cal M})$ changes for the values of $C$.

  1. If $C=-3$, for all values of $x>0$ there exist ${\cal M}$ which satisfies $f({\cal M})=g(x)$. This corresponds to the two curves which pass through a critical X-point of $(x,{\cal M})=(1,1)$ in Figure 2.9.

  2. If $C < -3$, the minimum of $g(x)$ is smaller than that of $f({\cal M})$. In this case, for $x$ where $g(x)<1=\min(f({\cal M}))$, there is no solution. Thus, $f({\cal M})=g(x)$ has solutions for $x<x_1$ and $x>x_2$, where $x_1 < x_2$, and $g(x_1)=g(x_2)=1$. This corresponds to the curves running perpendicularly in Figure 2.9.

  3. If $C > -3$, the minimum of $f({\cal M})$ is smaller than that of $g(x)$. In this case, $f({\cal M})=g(x)$ has solutions for ${\cal M}<{\cal M}_1$ and ${\cal M}>{\cal M}_2$, where ${\cal M}_1 < {\cal M}_2$, and $f({\cal M}_1)=f({\cal M}_2)=C+4$. This corresponds to curves running horizontally in Figure 2.9.

Figure 2.9: ${\cal M}$ vs $x$ (the $x$-axis is the normalized distance from the center $x\equiv r/(GM_*/2c_{is}^2)$ and the $y$-axis is the Mach number ${\cal M}$.)
\begin{figure}{\centering\leavevmode
\epsfxsize =.9\columnwidth \epsfbox{eps/M-x.ps}}\vspace*{5mm}\end{figure}

Out of the two solutions of $C=-3$, a wind solution is one with increasing ${\cal M}$ while departing from the star. This shows us the outflow speed is slow near the star but it is accelerated and a supersonic wind blows after passing the critical point. Since the equations are unchanged even if we replace $u$ with $-u$, the above solution is valid for an accreting flow $u < 0$. Considering such a flow, the solution running from the lower-right corner to the upper-left corner represents the accretion flow, in which the inflow velocity is accelerated reaching the star and finally accretes on the star surface with a supersonic speed.

Kohji Tomisaka 2012-10-03