Stellar Winds Star v w Stellar Winds Geoffrey V. Bicknell
1 Characteristics of stellar winds Solar wind Velocity at earth s orbit: Density: Temperature: Speed of sound: v 400 km/s n 10 7 m 3 c s T 10 5 K 50 km/s (1) (2) (3) (4) Astrophysical Gas Dynamics: Stellar Winds 2/66
Mass flux (spherically symmetric wind): M 4nmv r r 2 310 14 solar masses /yr (5) Other stars Red giants: M O&B type stars: 10 11 M sun /yr M 10 7 10 6 M sun /yr Protostars: M 10 4 M sun /yr Astrophysical Gas Dynamics: Stellar Winds 3/66
2 Why winds? Yohkoh Soft X-ray Telescope (SXT) full-field images from the Hiraiso Solar Terrestrial Research Center / CRL White-light Mk. 4 coronameter images http://umbra.nascom.nasa.gov/images/ Astrophysical Gas Dynamics: Stellar Winds 4/66
The Sun's outer atmosphere as it appears in ultraviolet light emitted by ionized oxygen flowing away from the Sun to form the solar wind (region outside black circle), and the disk of the Sun in light emitted by ionized iron at temperatures near two million degrees Celsius (region inside circle). This composite image taken by two instruments (UVCS, outer region and EIT, inner region) aboard the SOHO spacecraft shows dark areas called coronal holes at the poles and across the disk of the Sun where the highest speed solar wind originates. UVCS has discovered that the oxygen atoms flowing out of these regions have extremely high energies corresponding to temperatures of over 200 million degrees Celsius and accelerate to supersonic outflow velocities within 1.5 solar radii of the solar surface. The structure of the corona is controlled by the Sun's magnetic field which forms the bright active regions and the ray-like structures originating in the coronal holes. The composite image allows one to trace these structures from the base of the corona to millions of kilometers above the solar surface. (http://sohowww.nascom.nasa.gov/ gallery/uvcs/) Astrophysical Gas Dynamics: Stellar Winds 5/66
2.1 Hydrostatic atmospheres Analysis of the X-ray and radio emission from the solar corona indicates a temperature 10 6 K. Gas at this temperature cannot be held in by the gravitational field of the Sun.To show this, we first consider hydrostatic atmospheres and show that an hydrostatic atmosphere is not feasible. The momentum equations are: v i t + v i v j x j 1p -- xi x i Gravitational potential (6) (7) Astrophysical Gas Dynamics: Stellar Winds 6/66
Hydrostatic solutions v i 0 1p -- xi x i (8) Spherical symmetry 1 -- dp ----- dr d ----- dr (9) Astrophysical Gas Dynamics: Stellar Winds 7/66
Suppose we consider isothermal solutions with the pressure defined by the following equation: Boltzmanns constant 1.3810 J K 1 Temperature 10 6 K kt p --------- m Atomic mass unit Molecular weight 0.62 23 27 1.6610 Kg (10) Astrophysical Gas Dynamics: Stellar Winds 8/66
The potential exterior to a spherical star is Newtons gravitational constant 6.6710 11 SI units GM -------- r Mass of sun 30 2.010 kg (11) Astrophysical Gas Dynamics: Stellar Winds 9/66
For an isothermal atmosphere, the hydrostatic equation becomes: Integrating kt ------- 1 m --d ----- dr ln d ----- dr m ------- + Constant kt (12) Astrophysical Gas Dynamics: Stellar Winds 10/66
Suppose the density at a point r 0 within the corona is 0. Then the constant of integration is Constant ln 0 ln----- 0 m + ------- kt 0 m ------- kt 0 GMm ----------------- 1 kt -- r 1 ---- r 0 (13) ----- 0 exp GMm ----------------- 1 kt -- r 1 ---- r 0 Astrophysical Gas Dynamics: Stellar Winds 11/66
Since the pressure is proportional to the density, the pressure in this atmosphere is given by: p ----- p 0 GMm ----------------- 1 kt -- 1 exp r ---- r 0 exp GMm ----------------- 1 ktr 0 r 0 ---- r (14) The pressure at r the pressure cooker model The pressure described by the above equation decreases as r. However, it does not decrease to zero. The asymptotic value is P ------ P 0 exp GMm ----------------- ktr 0 (15) Astrophysical Gas Dynamics: Stellar Winds 12/66
The interpretation of this pressure is that this is what is required in addition to the gravitational field to confine the atmosphere. In the case of the sun, this confining pressure has to be provided by the interstellar medium. Astrophysical Gas Dynamics: Stellar Winds 13/66
2.2 Some numbers Excerpt from the Handbook of Astronomy and Astrophysics, available from: http://adsabs.harvard.edu/ books/hsaa/idx.html Astrophysical Gas Dynamics: Stellar Winds 14/66
For indicative calculations let us insert some numbers into the above equations, with the view in mind of determining what pressure we would need to hold the sun s corona in. The following parameters relate to the corona at a solar radius. 0.62 G 6.6710 6 m 1.6610 27 kg 11 T 1.510 K Solar radius 6.9610 m Electron density 1.5510 cm 3 1.5510 m 3 8 SI 8 14 (16) Astrophysical Gas Dynamics: Stellar Winds 15/66
For order of magnitude purposes, it is not necessary to distinguish between electron density and total number density. However, it is useful to go though the exercise, linking electron density to total number density and mass density. Number density, electron density and mass density in a fully ionised plasma Suppose that the plasma consists of ionised Hydrogen, and fully ionised Helium and the associated electrons. Then, Electron density n e n H + 2n He n H 1 + 2 n He --------- n H (17) Astrophysical Gas Dynamics: Stellar Winds 16/66
The total number density is: n n H + n He + n e 2n H n H 2 + + 3n He 3 n He --------- n H (18) The mass density is given by: n H m + 4n He m n H m1 + 4 n He --------- n H (19) Astrophysical Gas Dynamics: Stellar Winds 17/66
If we adopt a solar abundance of Helium, then n He --------- 0.085 n e 1.17n H n 2.26n H n H 1.34n H m 0.59m (20) Usually, based on more accurate composition and allowing for metals, i.e. elements heavier than Helium, we adopt and 0.62 Also from the above numbers: n 0.62nm 1.9n e (21) (22) Astrophysical Gas Dynamics: Stellar Winds 18/66
Required pressure for the pressure-cooker model p ------ p 0 exp GMm ----------------- ktr 0 (23) At a solar radius r 0 6.9610 m, we have 8 GMm ----------------- 7.2 ktr 0 p ------ 7.510 p 0 4 (24) The pressure at the base of the corona, is 3 p 0 nkt 1.9n e kt 8.110 N m 2 (25) Astrophysical Gas Dynamics: Stellar Winds 19/66
Therefore, the ISM pressure required by the pressure cooker model is: 4 3 p 7.510 8.110 N m 2 6.110 N m 2 (26) The ISM is a multi-phase medium, with all phases in pressure equilibrium. Take the warm phase, for example: 6 n 1 cm 3 10 6 m 3 T 10 4 K 13 p ISM 1.410 N m 2 (27) The ISM pressure fails by a factor of about 410 to contain the solar atmosphere! Hence, the solar atmosphere flows out in a wind. 7 Astrophysical Gas Dynamics: Stellar Winds 20/66
2.3 Preliminary estimate of mass flux In order to gain some idea of the mass flux we might expect, let us suppose that the solar wind flows out spherically from a solar radius at the escape velocity from the Sun. Escape velocity 2GM ------------ 12 / 5 6.210 m s 1 r 0 620 km s 1 (28) Astrophysical Gas Dynamics: Stellar Winds 21/66
Mass flux 4 r2 0 0 V esc 8 4 6.9610 2 1.2n e m 6.210 12 1.210 kg s 1 11 2.010 solar masses per year (29) This estimate is out by about three orders of magnitude mainly because our estimate of the density of the outflowing wind is in error. It is interesting to see how the correct theory makes allowance for this. 5 Astrophysical Gas Dynamics: Stellar Winds 22/66
3 Analysis of spherically symmetric winds 3.1 Fundamental equations Euler equations v i t Spherical symmetry + v i v j x j v r v r r 1p -- xi (30) See Landau & Lifshitz Fluid Mechanics for expressions for the gas dynamics equations in cylindrical and spehrical coordinates. 1p -- r x i r Astrophysical Gas Dynamics: Stellar Winds 23/66
For a spherical star of mass Equation of state r M GM -------- r GM -------- r 2 (31) Adiabatic: p Polytropic: p Ks (32) In the latter, is not necessarily c p c v. 5 3 Heating (in an expanding flow). C Astrophysical Gas Dynamics: Stellar Winds 24/66
Isothermal equation of state (33) Source of heat for the solar wind: Dissipation by waves in solar wind generated in photosphere/chromosphere. Speed of sound p kt ------- T constant m a s 2 dp d C 1 (34) where the derivative is no longer at constant entropy. We also define the Isothermal sound speed kt ------- m (35) Astrophysical Gas Dynamics: Stellar Winds 25/66
Radial momentum equation Pressure gradient: dp dr dp d ddr 2d a s dr (36) so that dv r a2 s d v r ----- dr dr GM -------- r 2 (37) Astrophysical Gas Dynamics: Stellar Winds 26/66
Mass flux t + x i vi r 2 v r 1 d ---- r 2 v dr r 0 r 2 Constant M ----- 4 (38) M From now on v. v r v i n i ds Sphere 4r 2 v r (39) Astrophysical Gas Dynamics: Stellar Winds 27/66
v r Also note that this equation gives us a way of estimating the density given the velocity. Solving for : M -------------- 4vr 2 (40) Astrophysical Gas Dynamics: Stellar Winds 28/66
3.2 Sonic point An important feature of all winds is the existence of a sonic point where the flow makes a transition from subsonic to supersonic flow. To show the existence of this we consider both the mass flux equation and the momentum equation. Mass flux 4vR 2 M (41) Astrophysical Gas Dynamics: Stellar Winds 29/66
Taking logs and differentiating: ln4 + ln + lnv + 2lnR lnm 1d 1dv 2 -- + -- + -- 0 dr v dr r (42) 1d -- dr 1dv 2 -- + -- v dr r Astrophysical Gas Dynamics: Stellar Winds 30/66
Substitute d 1 dr into momentum equation: dv v dr a2 s dv ----- v dr + 2 a s 2 ----- r GM -------- r 2 v 2 a2 dv s dr v 2 a s 2 ----- r GM -------- r 2 (43) dv dr v 2 a s 2 GM ----- -------- r r 2 --------------------------------- v 2 a2 s Astrophysical Gas Dynamics: Stellar Winds 31/66
The equation for v has a critical point where the numerator and denominator of the right hand side of this equation are both zero. That is, where: v 2 a s 2 ----- r a s GM -------- r c r 2 kt c ----------- m GM -------- 2a s 2 GMm ----------------- 2kT c (44) Astrophysical Gas Dynamics: Stellar Winds 32/66
e.g. the Sun ( 1): r c 4.810 9 v c T 170 km/s-------- 12 / 10 6 T -------- 1 T 10 6 metres 6.9R sun -------- 1 10 6 (45) Solar radius: R sun 6.9610 m. 8 Astrophysical Gas Dynamics: Stellar Winds 33/66
The actual measured coronal temperature is approximately 6 210 K. This implies that: 9 r c 2.010 m v c 340 km s 1 3.5 solar radii (46) 3.3 A better estimate of mass flux Let us suppose that the corona of the Sun is approximately determined by the hydrostatic solution out to the sonic point. Then ----- 0 exp GMm ----------------- 1 ktr 0 r ---- 0 r (47) Astrophysical Gas Dynamics: Stellar Winds 34/66
We have GMm ----------------- 7.2 ktr 0 r ---- 3.5 r 0 ----- 2.110 0 3 (48) Hence our mass flux estimate becomes: M 8 4 3.5 6.9610 2 2.110 1.2 1.5510 m3.410 13 We are now only about an order of magnitude out! 14 2.610 solar masses per year 5 3 (49) Astrophysical Gas Dynamics: Stellar Winds 35/66
4 Critical point analysis 4.1 Splitting equations If we want a wind solution which accelerates to supersonic then we need to negotiate the critical point: dv dr v 2a s 2 GM -------- -------- r r 2 --------------------------------- v 2 a2 s (50) Astrophysical Gas Dynamics: Stellar Winds 36/66
Split this equation into 2 by introducing a parameter the equations as: dr ----- r du 2 v 2 a2 s f 1 r v dv ----- v2a du 2 s r GM f 2 r v and write (51) These equations have no potential infinities anywhere and are much better behaved numerically and easier to analyse mathematically. u Astrophysical Gas Dynamics: Stellar Winds 37/66
4.2 Critical point This is now defined by: dr ----- du dv ----- 0 du Expansion in neighbourhood of critical point: dr du r r c + r' v v c + v' d dv r' du du d v' du (52) (53) Astrophysical Gas Dynamics: Stellar Winds 38/66
where d du r' v' f 1 r f 2 r f 1 v f 2 v f 1 r 2 v 2 a2 s f 2 v2a 2 s r GM r' v' (54) (55) Astrophysical Gas Dynamics: Stellar Winds 39/66
The functions f 1 and f 2 involve the function a2 s and in order to evaluate the partial derivatives of these functions, we need to evaluate the partial derivatives of Bernoulli s equation. a s 2. This involves evaluating 4.3 Bernoulli s equation: Since then p C 1 -- dp ----- C dr 2d ----- dr (56) (57) Astrophysical Gas Dynamics: Stellar Winds 40/66
and a s 2 dp ----- C d 1 (58) Astrophysical Gas Dynamics: Stellar Winds 41/66
Hence 1 --v 2 2 dv v dr C 2 d dr d 1 --v dr2 2 C 2 d dr 1 --v 2 2 C ---------- 1 1 C + ---------- 1 1 GM -------- constant r a2 s 1constant GM -------- r 2 GM -------- r 2 GM + -------- + constant r + GM -------- r 1 --v 2 2 (59) Astrophysical Gas Dynamics: Stellar Winds 42/66
The last equation implies: Since as r 2 1 GM -------- r 2 as 2 1v v f 1 r 2 v 2 a2 s f 2 v2a 2 s r GM (60) (61) then, at the critical point: f 1 2rv 2 a2 r s r 2 as 2 1GM r (62) f 1 v as 2 r 2 2v v r2 + 1 v (63) Astrophysical Gas Dynamics: Stellar Winds 43/66
f 2 v 2a2 r s + 2r as 2 v GM -------- 2 1 GM -------- r r r f 2 3 2 r v GM -------- r (64) f 2 2a v 2 s r GM + 2vr as 2 2 1a 2 v s r f 2 1GM v (65) Astrophysical Gas Dynamics: Stellar Winds 44/66
f 1 r f 2 r f 1 v f 2 v 1GM + 1v c r2 c GM -------- v r c 3 2 1GM c (66) and d du r' v' 1GM + 1v c r 2 GM -------- v r c 3 2 1GM c r' v' (67) Astrophysical Gas Dynamics: Stellar Winds 45/66
The solution of these equations requires the eigenvalues and eigenvectors of the matrix. Denoting the eigenvalues by : 1GM + 1v c r2 c GM -------- v c 3 2 + 1GM r c 0 (68) This simplifies to: 2 GM 2 5 3 ------------------ 2 GM 5 -------------- 3 12 / 2 (69) Astrophysical Gas Dynamics: Stellar Winds 46/66
Eigenvectors: 1GM + 1vr2 c X Y u 1 0 u 2 0 (70) 1GMu 1 + 1v c r 2 c u 2 0 Take u 2 1 u 1 + 1 v c r c 2 ---------- GM ------------------------------------------------------ 5 3 -------------- 12 / 1 2 (71) Astrophysical Gas Dynamics: Stellar Winds 47/66
The general solution in the neighbourhood of the critical point is: r' v' A 1 u 1 exp 1 u + A 2 u 2 exp 2 u (72) where u 1 and u 2 are the two independent eigenvectors. Slopes of lines through critical point and topology of the critical point: Astrophysical Gas Dynamics: Stellar Winds 48/66
V u 1 Critical Point u 2 u 2 r dv ----- dr u 2 ----- u 1 2v c ------- r c 5 3 -------------- 12 / 1 2 --------------------------------------------------- + 1 (73) Astrophysical Gas Dynamics: Stellar Winds 49/66
5 Scaling of wind equations for numerical solutions When we have analytical solutons of differential (or other) equations, the dependence on physical parameters is evident. However, when we construct numerical solutions, it is best to scale the equations in order to derive solutions depending upon the smallest number of physical parameters. Numerical solutions of the wind equations were given in section 5. The physical equations should be scaled before numerical solutions can be determined. One reason for this scaling is that the numerical solutions can then be presented using the smallest number of parameters. Various parameters (specifically the velocity, radius and density at the critical point) then enter as scaling parameters. Astrophysical Gas Dynamics: Stellar Winds 50/66
The wind equations in physical units are: dr ----- du r 2 v 2 a2 s dv ----- du v2a 2 s r GM (74) For numerical purposes we scale these by the values of velocity and radius at the critical point. Hence r r c r v v c v a s v c a s c (75) NB The primed variables here are not the same as the perturbations that were used to study solutions in the neighbourhood of the critical point. Astrophysical Gas Dynamics: Stellar Winds 51/66
5.1 Scaling of the sound speed The solutions that were shown in this lecture are all constrained to have the same mass flux. Hence, M ----- vr 4 2 c v c r2 c ----- c v ---- 1 r ---- 2 v c r c (76) i.e. v 1 r 2 (77) Also, a2 s C 1 a2 c C 1 c a s 2 ---- a c ----- 1 c (78) Astrophysical Gas Dynamics: Stellar Winds 52/66
i.e. a s 2 1 (79) Astrophysical Gas Dynamics: Stellar Winds 53/66
5.2 Normalisation of the differential equations dr r c ----- r du 2 c v c r 2 v 2 a s 2 1 dr' ---------- ----- r r c v2 du 2 v 2 a s 2 c dv' v c ------ r du c v 3 c v2a s 2 r GM ---------- r c v2 c (80) 2r c v 3 c va s 2 r 1 1 ---------- dv ------- 2a r c v2 du s 2 r 1 c Astrophysical Gas Dynamics: Stellar Winds 54/66
The last equation follows from the condition at the critical point GM ---------- 2 r c v2 c (81) We now simply make the transformation of the parameter u: du r c v c 2 du (82) and the normalised equations become: dr ------- r du 2 v 2 a s 2 dv ------- 2a du s 2 r 1 (83) Astrophysical Gas Dynamics: Stellar Winds 55/66
with v 1r 2 a s 2 1 (84) These can be easily solved with a numerical differential equation solver (I simply use the NAG library). To start solutions near the critical point, the perturbation equations, developed in this lecture, need to be used. Astrophysical Gas Dynamics: Stellar Winds 56/66
6 Numerical solution of equations V/V c 3 2 1 Velocity vs radius for γ 1.2 This plot shows the numerical integration of the two curves through the critical point and the numerical solution of other curves not passing through the critical point. 0 0 1 2 3 4 5 r/r c Blue curve: wind Red curve: accretion Astrophysical Gas Dynamics: Stellar Winds 57/66
1e+03 1e+02 Density for the two critical curves 1e+01 Density vs radius for γ1.2 Blue: wind ρ/ ρ c 1e+00 Red: accretion 1e-01 1e-02 1e-01 1e+00 1e+01 r/r c Astrophysical Gas Dynamics: Stellar Winds 58/66
7 Mass flux We can estimate the mass flux from the sun with this simple model. Estimate parameters at the critical point assuming a temperature of T c 10 6 K there. The mass flux is: M 4 c v c r c 2 (85) where the velocity and critical radius are given by v c kt -------- 1.27 10 m 2 km/s Astrophysical Gas Dynamics: Stellar Winds 59/66 (86)
and r c GM GMm -------- ----------------- 4.14 10 2a2 2kT 9 m 5.96 R sun s c (87) Using our numerical solution, we estimate that the density at the solar radius, i.e. at r --------- is 230. 5.96 c r c Hence, we can write the mass flux as M 4 c v c r2 c M c 4 0 ----- vc r2 c 0 (88) Astrophysical Gas Dynamics: Stellar Winds 60/66
Earlier we used an electron number density at the base of the corona of 8 n e 1.5510 cm 3 so that the density and 16 0 1.9n e m 3.010 gm cm 3 M 4 3.010 16 7 11 ------------------------ 1.2710 4.110 2 230 13 M 3.610 gm s 1 13 5.710 M sun yr 1 (89) (90) (91) Astrophysical Gas Dynamics: Stellar Winds 61/66
This is about a factor of 20 larger than the observed value. In order to obtain better agreement with observation, theory needs to take into account the effect of magnetic fields and coronal holes. Astrophysical Gas Dynamics: Stellar Winds 62/66
8 Asymptotic wind solutions 8.1 Deductions from Bernoulli s equation As r the density in the wind goes to zero and therefore since 1 --v 2 2 C + ---------- 1 GM -------- r constant 1 --v 2 a2 s GM + ---------- -------- 1 r constant (92) Astrophysical Gas Dynamics: Stellar Winds 63/66
We can use Bernoulli s equation to estimate the asymptotic velocity as follows. The constant on the right hand side of the above equation can be estimated from the conditions at the critical point. constant a c 2 1 --v 2 GM 2 c + ---------- -------- v2 1 1 c -- 2 5 3 ------------------ v2 2 1 c r c As r, a2 GM s 0 and -------- 0 so that r + 1 ---------- 2 1 (93) v2 5 3 ---------------v 2 1 c v 5 3 ---------------v 1 c (94) Astrophysical Gas Dynamics: Stellar Winds 64/66
We have been using 1.2, in which case v 2.65v c 2.65 127 km s 1 320 km s 1 (95) Thhis desn t compare too badly with the observed value of 400 km s 1 measured at the Earth. Reasons for the discrepancy between theory and observation Magnetic effects are important in the initial phase of the solar wind The effective outside the solar region is greater than 1.2 8.2 Density in the asymptotic region Since M 4vr 2 Astrophysical Gas Dynamics: Stellar Winds 65/66
then asymptotically M ------------------ 4v r 2 (96) This expression for the density is frequently used for stellar winds well outside the sonic point where the wind can be considered to have achieved its terminal velocity,. As we have sen v with the sun, the sonic point is reasonably close to the star. Astrophysical Gas Dynamics: Stellar Winds 66/66