next up previous
Next: MHD Dynamo Theory Up: Magnetohydrodynamic Fluids Previous: Interplanetary Magnetic Field

Mass and Angular Momentum Loss

Let us consider what impact the solar wind has on solar evolution. The most obvious question is whether the mass loss due to the wind is actually significant. Using typical measured values [i.e., a typical solar wind velocity and particle density at 1 AU of $ 500\,{\rm km}\,{\rm s}^{-1}$ and $ 7\times 10^6\,\,{\rm m}^{-3}$ , respectively (Priest 1984)], the Sun is apparently losing mass at a rate of $ 3\times 10^{-14}\,M_\odot$ per year, where $ M_\odot=2\times 10^{30}\,{\rm kg}$ is the solar mass (Yoder 1995), implying a timescale for significant mass loss of $ 3\times 10^{13}$ years, or some $ 6,000$ times longer than the estimated $ 5\times 10^9$ year age of the Sun (Hansen, Kawaler, and Trimble 2004). Clearly, the mass carried off by the solar wind has a negligible effect on the Sun's evolution. Note, however, that many stars in the universe exhibit significant mass loss via stellar winds. This is particularly the case for late-type stars (Mestel 2012).

Let us now consider the angular momentum carried off by the solar wind. Angular momentum loss is a crucially important topic in astrophysics, because only by losing angular momentum can large, diffuse objects, such as interstellar gas clouds, collapse under the influence of gravity to produce small, compact objects, such as stars and proto-stars (Mestel 2012). Magnetic fields generally play a crucial role in angular momentum loss. This is certainly the case for the solar wind, where the solar magnetic field enforces co-rotation with the Sun out to the Alfvén radius, $ r_A$ . Thus, the angular momentum carried away by a particle of mass $ m$ is $ {\mit\Omega}\,r_A^{\,2}\,
m$ , rather than $ {\mit\Omega}\,R_\odot^{\,2}\,m$ . The angular momentum loss timescale is, therefore, shorter than the mass loss timescale by a factor $ (R_\odot/r_A)^2\simeq 1/2500$ , making the angular momentum loss timescale comparable to the solar lifetime. It is clear that magnetized stellar winds represent a very important vehicle for angular momentum loss in the universe (Mestel 2012). Let us investigate angular momentum loss via stellar winds in more detail.

Under the assumption of spherical symmetry and steady flow, the azimuthal momentum evolution equation for the solar wind, taking into account the influence of the interplanetary magnetic field, is written (Huba 2000a)

$\displaystyle \rho\,[({\bf V}\cdot\nabla){\bf V}]_\phi\equiv \rho\,\frac{u_r}{r...
...({\bf j}\times{\bf B})_\phi \equiv \frac{B_r}{\mu_0\,r}\frac{d(r\,B_\phi)}{dr}.$ (7.90)

The constancy of the mass flux [see Equation (7.64)] and the $ 1/r^{\,2}$ dependence of $ B_r$ [see Equation (7.83)] permit the immediate integration of the previous equation to give

$\displaystyle r\,u_\phi -\frac{r\,B_r\,B_\phi}{\mu_0\,\rho\,u_r} = L,$ (7.91)

where $ L$ is the angular momentum per unit mass carried off by the solar wind. In the presence of an azimuthal wind velocity, the magnetic field and velocity components are related by an expression similar to Equation (7.81):

$\displaystyle \frac{B_r}{B_\phi} = \frac{u_r}{u_\phi - {\mit\Omega}\,r\,\sin\theta}.$ (7.92)

The fundamental physics assumption underlying the previous expression is the absence of an electric field in the frame of reference co-rotating with the Sun. Using Equation (7.92) to eliminate $ B_\phi$ from Equation (7.91), we obtain (in the ecliptic plane, where $ \sin\theta=1$ )

$\displaystyle r\,u_\phi = \frac{L\,M_A^{\,2} - {\mit\Omega}\,r^{\,2}}{M_A^{\,2} - 1},$ (7.93)


$\displaystyle M_A = \sqrt{\frac{u_r^{\,2}}{B_r^{\,2}/\mu_0\,\rho}}$ (7.94)

is the radial Alfvén Mach number. The radial Alfvén Mach number is small near the base of the corona, and about 10 at 1 AU: it passes through unity at the Alfvén radius, $ r_A$ , which is about $ 0.25$ AU from the Sun. The zero denominator on the right-hand side of Equation (7.93) at $ r=r_A$ implies that $ u_\phi$ is finite and continuous only if the numerator is also zero at the Alfvén radius. This condition then determines the angular momentum content of the outflow via

$\displaystyle L = {\mit\Omega}\,r_A^{\,2}.$ (7.95)

Note that the angular momentum carried off by the solar wind is indeed equivalent to that which would be carried off were coronal plasma to co-rotate with the Sun out to the Alfvén radius, and subsequently outflow at constant angular velocity. Of course, the solar wind does not actually rotate rigidly with the Sun in the region $ r<r_A$ : much of the angular momentum in this region is carried in the form of electromagnetic stresses.

It is easily demonstrated that the quantity $ M_A^{\,2}/(u_r\,r^{\,2})$ is a constant (because $ B_r\propto r^{\,-2}$ , and $ r^{\,2}\,\rho\,u_r$ is constant), and can, therefore, be evaluated at $ r=r_A$ to give

$\displaystyle M_A^{\,2} = \frac{u_r\,r^{\,2}}{u_{r\,A}\,r_A^{\,2}},$ (7.96)

where $ u_{r\,A} \equiv u_r(r_A)$ . Equations (7.93), (7.95), and (7.96) can be combined to produce

$\displaystyle u_\phi = \frac{{\mit\Omega}\,r}{u_{r\,A}} \frac{u_{r\,A} - u_r}{1-M_A^{\,2}}.$ (7.97)

In the limit $ r\rightarrow\infty$ , we have $ M_A\gg 1$ , so the previous expression yields

$\displaystyle u_\phi \rightarrow {\mit\Omega}\,r_A\left(\frac{r_A}{r}\right)\left(1-\frac{u_{r\,A}} {u_r}\right)$ (7.98)

at large distances from the Sun. Recall, from Section 7.7, that if the coronal plasma were to simply co-rotate with the Sun out to $ r=r_A$ , and experience no torque beyond this radius, then we would expect

$\displaystyle u_\phi \rightarrow {\mit\Omega}\,r_A\left(\frac{r_A}{r}\right)$ (7.99)

at large distances from the Sun. The difference between the previous two expressions is the factor $ 1-u_{r\,A}/u_r$ , which is a correction for the angular momentum retained by the magnetic field at large $ r$ .

The previous analysis presented was first incorporated into a quantitative coronal expansion model by Weber and Davis (Weber and Davis 1967). The model of Weber and Davis is very complicated. For instance, the solar wind is required to flow smoothly through no less than three critical points. These are associated with the sound speed (as in Parker's original model), the radial Alfvén speed, $ B_r/\!\sqrt{\mu_0\,\rho}$ , (as previously described), and the total Alfvén speed, $ B/\!\sqrt{\mu_0\,\rho}$ . Nevertheless, the simplified analysis outlined in this section captures most of the essential features of the outflow.

next up previous
Next: MHD Dynamo Theory Up: Magnetohydrodynamic Fluids Previous: Interplanetary Magnetic Field
Richard Fitzpatrick 2016-01-23