ASTR 370 · Lecture 08

The Solar Wind: Parker’s Transonic Wind Solution

Deriving why the corona cannot simply sit in hydrostatic equilibrium

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.

Reading anchor: OpenStax Astronomy 2e (local extracted reference copy: references/openstax-astronomy-2e-extracted.txt) covers the Sun in Chapter 15 ("The Sun: A Garden-Variety Star"): 15.1 The Structure and Composition of the Sun, 15.2 The Solar Cycle, and 15.3 Solar Activity above the Photosphere. Chapter 15’s section numbering was cross-checked directly against its own embedded figure numbering this session (Figures 15.1 through 15.24 all appear under their matching section headers with no offset), so this citation is verified against the local extracted text rather than extrapolated by analogy from a prior course. OpenStax’s treatment throughout Chapter 15 is descriptive and qualitative; nearly every quantitative derivation in this course (the hydrostatic solar-atmosphere/scale-height argument, the sunspot magnetic-pressure balance, the Parker solar-wind critical-point argument and its numerical transonic solution, flare free-energy budgets, the Chapman-Ferraro magnetopause standoff distance, and the Dessler-Parker-Sckopke ring-current/Dst relation) goes well beyond OpenStax Astronomy 2e’s introductory level; each lecture states explicitly where this course extends beyond the assigned reading.

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).

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.

First question: what here is directly observed or a published/live-verified measurement, and what follows only once a physical law is derived and applied?

Terms for This Lecture

hydrostatic corona (inconsistency)

Definition: The unphysical result that a sufficiently hot static corona retains finite pressure at very large radius instead of matching interstellar conditions.

Why it matters here: 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.

Relationship: Use this term alongside Parker solar wind; the lecture's quantitative model is summarized below.

\[ 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} \]

Boundary or common confusion: 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.

Evidence We Need to Explain

Derivation and Model

Quantitative Tool

\[ 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.

Visual Reasoning

Trace the reasoning chain from raw observation or live-verified/published data to derived physical law to inferred quantity in this lecture’s figure.

  • Which stage is directly observed or a real published measurement?
  • Which stage is the physical law or derivation step?
  • What would change if an assumption in the derivation failed?
Lecture 08: The Solar Wind: Parker’s Transonic Wind Solutionr / r_cv / c_s (transonic solution)Critical point: r=r_c, v=c_sr_c=4.62 R_sun (T_corona=1.5e+06 K, c_s=144 km/s). Numerically solved transonic branch gives v(1 AU)~559 km/s, the right order of magnitude for the real observed slow solar wind (400 km/s).
The Solar Wind: Parker’s Transonic Wind Solution: from observation to inferred quantity.

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 Segment

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

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.

References