Lab 04: CME Transit Time and the Parker Wind's Numerical Solution

Compare naive and drag-based CME transit-time estimates to the real observed value, and numerically solve the Parker transonic wind equation at several radii.

Apparatus and Setup

This lab computes both the naive constant-speed and a simple drag-based transit-time estimate for the real 6 September 2017 CME (Lecture 07), then reproduces Lecture 08's numerical transonic Parker-wind solve at several radii.

Materials and Data

CME LASCO speed: 1571 km/s (level 2). Ambient slow solar wind speed: 400 km/s (level 2, OMNI climatology). Coronal temperature: 1.5e+06 K, mu=0.60. Parker critical radius: r_c=4.62 R_sun (computed in Lecture 08).

Procedure

  1. Compute the naive constant-speed travel time (1 AU/v_CME) and confirm it matches Lecture 07's stated value.
  2. Using the drag-based deceleration functional form from Lecture 07's figure (exponential relaxation toward the ambient wind speed with a 15-hour time constant), numerically integrate or evaluate the CME's distance vs time and estimate the arrival time at 1 AU; compare it to the real observed ~36-hour transit time.
  3. Numerically solve the Parker transonic wind equation (bisection on the physical branch, as in Lecture 08) for u=v/c_s at r/r_c=0.5, 1.0, 2.0, and 47 (1 AU), and tabulate v in km/s at each radius.
  4. Confirm that u=1.0 (v=c_s) exactly at r/r_c=1.0, the defining property of the critical point.

Analysis and Uncertainty

Discuss why the drag-based model's specific functional form (and its 15-hour time constant) is an illustrative, level-3 choice calibrated only to roughly reproduce the real observed transit time, not a first-principles solution of the full drag-force differential equation, which would require knowing the ambient solar wind's density profile in detail.

Deliverables

Assessment Criteria

References and Provenance