Lecture 08: The Solar Wind: Parker’s Transonic Wind Solution

ASTR 370 Solar Physics and Space Weather

Context and Why This Matters

Lecture 07’s CME kinematics assumed an ambient solar wind for the CME to decelerate against; this lecture derives, from first principles, why the corona cannot be static at all and must instead flow outward as a supersonic wind -- one of the foundational results of heliophysics, first derived by Eugene Parker (the Parker Solar Probe’s namesake, Unit 5).

Learning Goals

  1. Show that a static, hot corona in simple hydrostatic equilibrium is physically inconsistent with a vanishing pressure at infinity.
  2. Derive the isothermal Parker wind’s critical-point condition and numerically solve for the transonic wind-speed profile.
  3. Compare the numerically solved solar wind speed to real, published solar wind observations at 1 AU.

Opening Phenomenon

A naive hydrostatic corona (using Lecture 01’s scale-height framework at the corona’s real temperature, T=1.5e+06 K) predicts a finite, non-zero gas pressure even at an infinite distance from the Sun -- an unphysical result, since the corona must eventually meet the near-vacuum of interstellar space, and this contradiction is precisely what led Eugene Parker in 1958 to propose that the corona instead flows outward continuously as a wind, reaching supersonic speeds well before 1 AU.

Terms Developed in Context

hydrostatic corona (inconsistency)

The unphysical result that a sufficiently hot static corona retains finite pressure at very large radius instead of matching interstellar conditions. In this lecture, the term is needed to interpret the following evidence or model claim: A naive hydrostatic corona (using Lecture 01’s scale-height framework at the corona’s real temperature, T=1.5e+06 K) predicts a finite, non-zero gas pressure even at an infinite distance from the Sun -- an unphysical result, since the corona must eventually meet the near-vacuum of interstellar space, and this contradiction is precisely what led Eugene Parker in 1958 to propose that the corona instead flows outward continuously as a wind, reaching supersonic speeds well before 1 AU.

Its relationship to Parker solar wind is part of the lecture's reasoning chain rather than a synonym: For a static, isothermal corona, hydrostatic equilibrium (Lecture 01) integrates to a pressure that falls off but asymptotes to a finite, non-zero value as r goes to infinity -- physically impossible if the corona must match onto a near-vacuum interstellar medium, forcing the conclusion that the corona cannot be static. The working relation \( u^2 - \ln(u^2) = 4\ln\!\left(\frac{r}{r_c}\right) + \frac{4 r_c}{r} - 3, \qquad r_c = \frac{GM_\odot}{2c_s^2},\ u \equiv \frac{v}{c_s} \) makes the relevant observable or model dependence explicit.

Limitation and usage: Assuming the Parker wind equation has a simple algebraic closed-form solution for v(r) at all radii. Only the critical point itself (r=r_c, v=c_s) has a simple closed form; the full transonic profile requires numerically solving a transcendental equation (as this lecture’s figure does), and naively rearranging the equation algebraically for v(r) is not possible in general. This caution sets the boundary on what can be inferred when using hydrostatic corona (inconsistency) in an explanation.

Parker solar wind

The steady hydrodynamic outflow obtained when a hot corona passes smoothly from subsonic to supersonic speed. In this lecture, the term is needed to interpret the following evidence or model claim: Parker Solar Probe (Unit 5) has directly sampled the solar wind acceleration region within a few tens of solar radii, providing direct in-situ tests of the transonic wind theory this lecture derives.

Its relationship to critical point is part of the lecture's reasoning chain rather than a synonym: A naive hydrostatic corona (using Lecture 01’s scale-height framework at the corona’s real temperature, T=1.5e+06 K) predicts a finite, non-zero gas pressure even at an infinite distance from the Sun -- an unphysical result, since the corona must eventually meet the near-vacuum of interstellar space, and this contradiction is precisely what led Eugene Parker in 1958 to propose that the corona instead flows outward continuously as a wind, reaching supersonic speeds well before 1 AU. The working relation \( u^2 - \ln(u^2) = 4\ln\!\left(\frac{r}{r_c}\right) + \frac{4 r_c}{r} - 3, \qquad r_c = \frac{GM_\odot}{2c_s^2},\ u \equiv \frac{v}{c_s} \) makes the relevant observable or model dependence explicit.

Limitation and usage: Assuming the Parker wind equation has a simple algebraic closed-form solution for v(r) at all radii. Only the critical point itself (r=r_c, v=c_s) has a simple closed form; the full transonic profile requires numerically solving a transcendental equation (as this lecture’s figure does), and naively rearranging the equation algebraically for v(r) is not possible in general. This caution sets the boundary on what can be inferred when using Parker solar wind in an explanation.

critical point

The radius where the Parker wind speed equals the sound speed and the differential equation becomes singular unless its numerator also vanishes. In this lecture, the term is needed to interpret the following evidence or model claim: Parker (1958) instead solved the steady-state, spherically symmetric momentum equation for an isothermal outflow, finding that physically acceptable solutions exist only along specific branches of a family of curves, with the unique physically realized solution (matching low speed near the Sun and continuing to accelerate outward, becoming supersonic) passing through a critical point at r_c=GM_sun/(2 c_s^2), where the flow speed equals the isothermal sound speed c_s.

Its relationship to transonic solution is part of the lecture's reasoning chain rather than a synonym: Parker’s transonic wind solution -- derived here by confronting hydrostatic equilibrium’s failure at infinity and numerically solving the resulting critical-point equation -- explains why the corona continuously flows outward as a supersonic wind, the medium Unit 6 will show couples the Sun’s activity directly to Earth’s magnetosphere. The working relation \( u^2 - \ln(u^2) = 4\ln\!\left(\frac{r}{r_c}\right) + \frac{4 r_c}{r} - 3, \qquad r_c = \frac{GM_\odot}{2c_s^2},\ u \equiv \frac{v}{c_s} \) makes the relevant observable or model dependence explicit.

Limitation and usage: Assuming the Parker wind equation has a simple algebraic closed-form solution for v(r) at all radii. Only the critical point itself (r=r_c, v=c_s) has a simple closed form; the full transonic profile requires numerically solving a transcendental equation (as this lecture’s figure does), and naively rearranging the equation algebraically for v(r) is not possible in general. This caution sets the boundary on what can be inferred when using critical point in an explanation.

transonic solution

The unique physically continuous wind branch that crosses the critical point from subsonic to supersonic flow. In this lecture, the term is needed to interpret the following evidence or model claim: The full nonlinear equation for this transonic solution, u^2 - ln(u^2) = 4 ln(r/r_c) + 4 r_c/r - 3 (u=v/c_s), has no closed-form algebraic solution for v(r) and must be solved numerically (root-finding on the correct physical branch) at each radius -- exactly the numerical approach this lecture’s worked example and figure use, rather than an algebraic shortcut.

Its relationship to isothermal sound speed is part of the lecture's reasoning chain rather than a synonym: Parker’s transonic wind solution -- derived here by confronting hydrostatic equilibrium’s failure at infinity and numerically solving the resulting critical-point equation -- explains why the corona continuously flows outward as a supersonic wind, the medium Unit 6 will show couples the Sun’s activity directly to Earth’s magnetosphere. The working relation \( u^2 - \ln(u^2) = 4\ln\!\left(\frac{r}{r_c}\right) + \frac{4 r_c}{r} - 3, \qquad r_c = \frac{GM_\odot}{2c_s^2},\ u \equiv \frac{v}{c_s} \) makes the relevant observable or model dependence explicit.

Limitation and usage: Assuming the Parker wind equation has a simple algebraic closed-form solution for v(r) at all radii. Only the critical point itself (r=r_c, v=c_s) has a simple closed form; the full transonic profile requires numerically solving a transcendental equation (as this lecture’s figure does), and naively rearranging the equation algebraically for v(r) is not possible in general. This caution sets the boundary on what can be inferred when using transonic solution in an explanation.

isothermal sound speed

The characteristic pressure-wave speed c_s equal to the square root of kT divided by mean-particle mass for an isothermal gas. In this lecture, the term is needed to interpret the following evidence or model claim: Parker (1958) instead solved the steady-state, spherically symmetric momentum equation for an isothermal outflow, finding that physically acceptable solutions exist only along specific branches of a family of curves, with the unique physically realized solution (matching low speed near the Sun and continuing to accelerate outward, becoming supersonic) passing through a critical point at r_c=GM_sun/(2 c_s^2), where the flow speed equals the isothermal sound speed c_s.

Its relationship to terminal wind speed is part of the lecture's reasoning chain rather than a synonym: A naive hydrostatic corona (using Lecture 01’s scale-height framework at the corona’s real temperature, T=1.5e+06 K) predicts a finite, non-zero gas pressure even at an infinite distance from the Sun -- an unphysical result, since the corona must eventually meet the near-vacuum of interstellar space, and this contradiction is precisely what led Eugene Parker in 1958 to propose that the corona instead flows outward continuously as a wind, reaching supersonic speeds well before 1 AU. The working relation \( u^2 - \ln(u^2) = 4\ln\!\left(\frac{r}{r_c}\right) + \frac{4 r_c}{r} - 3, \qquad r_c = \frac{GM_\odot}{2c_s^2},\ u \equiv \frac{v}{c_s} \) makes the relevant observable or model dependence explicit.

Limitation and usage: Assuming the Parker wind equation has a simple algebraic closed-form solution for v(r) at all radii. Only the critical point itself (r=r_c, v=c_s) has a simple closed form; the full transonic profile requires numerically solving a transcendental equation (as this lecture’s figure does), and naively rearranging the equation algebraically for v(r) is not possible in general. This caution sets the boundary on what can be inferred when using isothermal sound speed in an explanation.

terminal wind speed

The asymptotic or far-heliosphere bulk speed approached by a solar-wind stream, not a universal constant for all streams. In this lecture, the term is needed to interpret the following evidence or model claim: A naive hydrostatic corona (using Lecture 01’s scale-height framework at the corona’s real temperature, T=1.5e+06 K) predicts a finite, non-zero gas pressure even at an infinite distance from the Sun -- an unphysical result, since the corona must eventually meet the near-vacuum of interstellar space, and this contradiction is precisely what led Eugene Parker in 1958 to propose that the corona instead flows outward continuously as a wind, reaching supersonic speeds well before 1 AU.

Its relationship to hydrostatic corona (inconsistency) is part of the lecture's reasoning chain rather than a synonym: The real, measured solar wind speed at 1 AU is bimodal: a "slow" wind (~400 km/s, originating from streamer-belt regions near the boundary between open and closed magnetic field) and a "fast" wind (~700-800 km/s, originating from coronal holes, regions of open magnetic field), both real, standard published climatological ranges (OMNI database). The working relation \( u^2 - \ln(u^2) = 4\ln\!\left(\frac{r}{r_c}\right) + \frac{4 r_c}{r} - 3, \qquad r_c = \frac{GM_\odot}{2c_s^2},\ u \equiv \frac{v}{c_s} \) makes the relevant observable or model dependence explicit.

Limitation and usage: Assuming the Parker wind equation has a simple algebraic closed-form solution for v(r) at all radii. Only the critical point itself (r=r_c, v=c_s) has a simple closed form; the full transonic profile requires numerically solving a transcendental equation (as this lecture’s figure does), and naively rearranging the equation algebraically for v(r) is not possible in general. This caution sets the boundary on what can be inferred when using terminal wind speed in an explanation.

Evidence

Derivation and Model

Working Equation

\[ u^2 - \ln(u^2) = 4\ln\!\left(\frac{r}{r_c}\right) + \frac{4 r_c}{r} - 3, \qquad r_c = \frac{GM_\odot}{2c_s^2},\ u \equiv \frac{v}{c_s} \]

Worked Example

  1. For the corona’s real characteristic temperature (T=1.5e+06 K, fully ionized, mu=0.60), the isothermal sound speed is c_s=144 km/s, giving a critical radius r_c=3.22x10^9 m = 4.62 R_sun -- well within the corona, meaning the wind should already be supersonic before reaching a few solar radii.
  2. Numerically solving the transonic-branch equation at r=1 AU (r/r_c=46.5) gives a predicted wind speed of v(1 AU)~559 km/s -- the right order of magnitude compared to the real, published slow solar wind speed at 1 AU (400 km/s), though this simplified isothermal, single-temperature model does not reproduce the real wind’s observed bimodal (slow/fast) structure, which requires a more realistic, non-isothermal, magnetically structured corona -- an honest, disclosed limitation, not a claimed precise match.

Common Misconception

Assuming the Parker wind equation has a simple algebraic closed-form solution for v(r) at all radii. Only the critical point itself (r=r_c, v=c_s) has a simple closed form; the full transonic profile requires numerically solving a transcendental equation (as this lecture’s figure does), and naively rearranging the equation algebraically for v(r) is not possible in general.

Active-Learning Guidance

Using the numerically solved transonic branch, explain qualitatively why choosing the wrong root of the transcendental equation (the subsonic-everywhere "breeze" solution rather than the physically realized transonic solution) would predict a solar wind that never becomes supersonic, contradicting the direct spacecraft evidence cited above.

Lab Connection

Lab 04’s second half reproduces this lecture’s numerical transonic-wind solve at several radii and compares the resulting 1 AU wind speed to real, published OMNI-database climatological values.

Synthesis Questions

References