The solar photosphere radiates into space at a modest effective temperature of (5{,}778 ext{ K}). Yet just a few thousand kilometers above the visible surface, the tenuous solar corona suddenly leaps to temperatures exceeding two to three million Kelvin. In direct violation of naive thermodynamics, heat flows from a cold reservoir to a hot one. In 1958, Eugene Parker resolved how this scorching atmosphere continuously blows out into interplanetary space as the Solar Wind. In this treatise, we derive Parker’s transonic expansion and explore solutions to the coronal heating paradox.
1. Parker’s Hydrodynamic Solar Wind Model
Consider a steady, spherically symmetric, isothermal corona of temperature (T_0). The plasma consists of ionized hydrogen with mass density (
ho(r) = m_p n(r)) and isothermal sound speed (c_s quiv sqrt{2 k_B T_0 / m_p} = ext{const}). The equation of state is (P =
ho c_s^2).
The steady-state mass continuity equation is:
[
rac{d}{dr}left( 4pi r^2
ho v
ight) = 0 quad implies quad r^2
ho v = ext{const} quad implies quad rac{1}{
ho}rac{d
ho}{dr} = -rac{2}{r} – rac{1}{v}rac{dv}{dr}
]
The radial momentum conservation equation under solar gravity is:
[
v rac{dv}{dr} = -rac{1}{
ho}rac{dP}{dr} – rac{G M_odot}{r^2} = -c_s^2 rac{1}{
ho}rac{d
ho}{dr} – rac{G M_odot}{r^2}
]
Substituting the density gradient from mass continuity:
[
v rac{dv}{dr} = -c_s^2 left( -rac{2}{r} – rac{1}{v}rac{dv}{dr}
ight) – rac{G M_odot}{r^2}
]
Regrouping velocity derivatives yields the celebrated Parker Solar Wind Equation:
[
left( v – rac{c_s^2}{v}
ight) rac{dv}{dr} = rac{2 c_s^2}{r} left( 1 – rac{G M_odot}{2 c_s^2 r}
ight) = rac{2 c_s^2}{r} left( 1 – rac{r_c}{r}
ight)
]
where the Sonic Critical Radius (r_c) is defined by:
[
r_c quiv rac{G M_odot}{2 c_s^2}
]
2. Transonic de Laval Nozzle Topology
The Parker equation is structurally identical to the flow of compressible gas through a de Laval nozzle, where gravity and spherical geometric divergence combine to act as an effective converging-diverging throat:
- For (r < r_c), the right-hand side is negative ((1 – r_c/r 0)), the left-hand side must also be negative, requiring subsonic flow: (v < c_s).
- For (r > r_c), the right-hand side is positive ((1 – r_c/r > 0)). Sustained acceleration requires (v > c_s), yielding supersonic expansion.
- At the critical radius (r = r_c), the flow passes smoothly through the sound barrier: (v(r_c) = c_s).
Integrating the Parker equation produces the transcendental energy integral:
[
rac{v^2}{c_s^2} – lnleft( rac{v^2}{c_s^2}
ight) = 4 lnleft( rac{r}{r_c}
ight) + 4 rac{r_c}{r} – 3
]
Among the family of mathematical trajectories, only the unique transonic solution starts subsonic near the solar base ((v ll c_s)) and expands asymptotically into a supersonic wind with finite terminal pressure matching the interstellar medium, exactly as confirmed by the Mariner 2 and Parker Solar Probe missions.
3. Resolving the Coronal Heating Paradox
Why does the corona reach millions of Kelvin in the first place? Energy cannot be transported by thermal conduction or photons. Instead, mechanical convective motions in the photosphere inject non-thermal energy into magnetic fields via two primary mechanisms:
- Wave Heating (AC Mechanism): Convective granulation shakes photospheric footpoints, launching high-frequency shear Alfvén waves upward into coronal loops. As these waves propagate along magnetic flux tubes, nonlinear wave-wave interactions cascade energy to kinetic scales, where ion-cyclotron resonance and Landau damping dissipate the waves into thermal plasma motion.
- Magnetic Braiding & Nanoflares (DC Mechanism): Proposed by Eugene Parker (1988), random footpoint motions entangle and twist coronal magnetic field lines into complex topologies. When stress exceeds critical thresholds, the tangled field lines undergo microscopic magnetic reconnection events—Nanoflares of energy (sim 10^{24} ext{ ergs})—which impulsively ignite the ambient coronal plasma.
