Learning objectives / 学习目标
- 由 steady spherical flow 的质量与动量方程说明 solar wind 和 critical point；
- 推出 Parker-spiral 角度标度，区分 field geometry 与粒子轨道；
- 用 ram/magnetic pressure balance 估算 magnetopause，追踪 reconnection 至 auroral emission 的能量链；
- 区分 accretion 的引力能释放、jet launching/collimation 与 collisionless-shock dissipation；
- 由 induction equation、\(R_m\) 和 turbulent EMF 解释 small-scale/mean-field dynamo 与观测反演。
Prerequisites / 先修知识
- Plasma Fluid and MHD Models 中的 induction equation、frozen flux、magnetic pressure/tension；
- Plasma Waves 与 Plasma Instabilities；
- gravitation、shock conservation laws、synchrotron 和 Faraday rotation 的基本概念。
Scope and conventions / 范围与约定
- 默认 heliocentric inertial frame，\(r\) 向外；地球磁层图中 Sun 在左、solar wind 向右。
- dynamic pressure 取 \(P_{\rm dyn}=\rho v^2\)；若他处使用 \(\rho v^2/2\)，比较前先统一 convention。
- MHD 描述大尺度守恒与 topology；shock ramp、reconnection diffusion region 与 particle injection 需 kinetic theory。
- 天体 jet 可相对论；本课 Newtonian/MHD 关系是标度，不替代 GRMHD。
Notation and modeling contract / 符号与建模契约
| Symbol | Meaning | SI / convention |
|---|---|---|
| \(r,v_r\) | central distance, radial bulk speed | m, \(\mathrm{m\,s^{-1}}\) |
| \(\Omega_\odot\) | solar rotation angular speed | \(\mathrm{rad,s^{-1}}\) |
| \(P_{\rm dyn}\) | momentum-flux pressure | Pa |
| \(M_A=v/v_A\) | Alfvén Mach number | 1 |
| \(r_m\) | magnetopause stand-off distance | m |
| \(\eta_m=1/(\mu_0\sigma)\) | magnetic diffusivity | \(\mathrm{m^2,s^{-1}}\) |
| \(R_m=UL/\eta_m\) | magnetic Reynolds number | 1 |
| \(\boldsymbol{\mathcal E}=\langle\mathbf u'\times\mathbf b'\rangle\) | turbulent electromotive force | \(\mathrm{V,m^{-1}}\) |
1. Solar wind and heliosphere / 太阳风与日球层
Knowledge ID: P13-U04-L02-K01
对 steady spherical wind，质量守恒给出
\[4\pi r^2\rho v_r=\dot M=\mathrm{constant}.\]
isothermal Parker model 还给出
\[\left(v_r-\frac{c_s^2}{v_r}\right)\frac{dv_r}{dr} =\frac{2c_s^2}{r}-\frac{GM_\odot}{r^2}.\]
光滑穿越 \(v_r=c_s\) 的 critical solution 从亚声速过渡到超声速。这是 collective pressure-driven outflow，不是若干独立粒子偶然超过 escape speed。
图 1 区分 solar-wind radial bulk flow 与受星际介质相对流压塑形的非球对称日球边界；nose–tail 方向是示意，真实形状随太阳活动与局域星际条件变化。
图 2 显示 Parker wind 分支穿过 critical point；真实 corona 非等温且有 wave/turbulent heating。
ideal frozen flux 与 solar rotation 产生 Parker spiral：
\[\tan\psi=\frac{\lvert B_\phi\rvert}{B_r} \simeq\frac{\Omega_\odot r}{v_r}, \qquad B_r\propto r^{-2},\quad \lvert B_\phi\rvert\propto r^{-1}.\]
图 3 的曲线是等速稳态近似下的场线，不是单粒子轨道。
2. Magnetospheres and aurorae / 磁层与极光
Knowledge ID: P13-U04-L02-K02
行星磁场与 stellar wind 形成 bow shock、magnetosheath、magnetopause 和 magnetotail。最简 pressure balance 为
\[P_{\rm dyn}+P_{\rm th}+P_{B,\rm sw} \simeq\frac{B_{\rm planet}^2(r_m)}{2\mu_0}+P_{\rm inner}.\]
对 dipole \(B\propto r^{-3}\)，有 \(r_m\propto P_{\rm dyn}^{-1/6}\)：边界对 pressure 只有六次方根响应。
图 4 中 Sun 在左；形状是日风压缩的 schematic。
reconnection 将能量送入 magnetosphere；field-aligned particles 沉降到上层大气，碰撞激发 O、N\(_2\) 等，去激发生成 auroral photons。颜色依赖 species、state、altitude 与 collisional quenching，不是入射电子的“颜色”。
图 5 中红箭头是 guiding-center flux，不是实际 gyro-orbit。
图 6 说明 aurora 是 magnetosphere–ionosphere 耦合系统的可视输出。
3. Accretion, jets, and collisionless shocks / 吸积、喷流与无碰撞激波
Knowledge ID: P13-U04-L02-K03
带角动量 plasma 形成 accretion disk。Maxwell+Reynolds stresses 将角动量向外输送，质量向内；Newtonian energy-release 标度约 \(GM\dot M/(2r_{\rm in})\)。大尺度 magnetic flux 可从 disk/compact-object rotation 取能并 launching/collimating jets；inner disk 需 relativistic treatment。
图 7 在 disk 平面内区分质量向内、应力携带角动量径向向外和 orbital motion；上下红箭头仅表示 bipolar jets。喷流本身也可携带角动量，但不是图中的绿色 disk-stress 箭头；结构未按比例。
space plasma 的 Coulomb mean free path 可远大于 shock ramp，因而 collective fields、waves 和 reflected particles 提供耗散。Rankine–Hugoniot 关系约束上/下游宏观状态，却不解释 microscopic entropy production 和 nonthermal injection。
图 8 中宏观 bulk speed 降低不表示每个粒子单调减速。
4. Dynamo processes and cosmic magnetic fields / 发电机过程与宇宙磁场
Knowledge ID: P13-U04-L02-K04
resistive-MHD induction equation 与 magnetic Reynolds number 为
\[\frac{\partial\mathbf B}{\partial t} =\nabla\times(\mathbf u\times\mathbf B)+\eta_m\nabla^2\mathbf B, \qquad R_m=\frac{UL}{\eta_m}.\]
\(R_m\gg1\) 时 advection/stretching 在大尺度上占优，但局部 reconnection/dissipation 仍可重排 topology。small-scale dynamo 随机 stretch–fold seed fields；mean-field decomposition 产生
\[\boldsymbol{\mathcal E}=\langle\mathbf u'\times\mathbf b'\rangle,\]
与 differential rotation 一起组织 large-scale field。能量来自 rotation/convection/turbulence，不是无中生有。
图 9 是 turbulent dynamo 的 energy-feedback loop；saturation 不表示场停止变化。
宇宙磁场常由 synchrotron intensity/polarization、Faraday rotation 与 inverse-Compton 联合约束。Faraday rotation \(\propto\int n_eB_\parallel dl\)，synchrotron 还依赖 relativistic-electron spectrum，因此从辐射到 \(B\) 是 inverse problem。
图 10 是模型选择图：不同尺度需 kinetic、MHD 或观测 forward model。
Worked examples / 完整例题
Worked Example 1 — Solar-wind dynamic pressure
已知: \(n_p=5.00\ \mathrm{cm^{-3}}\)，\(v=400\ \mathrm{km,s^{-1}}\)；忽略 alpha particles 和 electron mass。
系统/方向: 地球静止，wind 法向入射；\(+x\) 指向地球。
模型: 冷单流体，\(P_{\rm dyn}=\rho v^2\)。
\[n_p=5.00\times10^6\ \mathrm{m^{-3}},\quad \rho=n_pm_p=8.36\times10^{-21}\ \mathrm{kg,m^{-3}},\]
\[P_{\rm dyn}=\rho v^2=1.34\times10^{-9}\ \mathrm{Pa}=1.34\ \mathrm{nPa}.\]
检查: 单位为 Pa，数量级为 nPa；三位有效数字。
Worked Example 2 — Parker angle at 1 AU
已知: \(\Omega_\odot=2.70\times10^{-6}\ \mathrm{rad,s^{-1}}\)，\(r=1.496\times10^{11}\ \mathrm m\)，\(v_r=4.00\times10^5\ \mathrm{m,s^{-1}}\)。
参考系/模型: heliocentric equatorial plane，steady ideal radial wind。
\[\tan\psi=\frac{\Omega_\odot r}{v_r}=1.01, \qquad \psi=45.2^\circ.\]
检查: 比值无量纲；wind 更快时 \(\psi\) 更小、场更径向。
Worked Example 3 — Dipole magnetopause
已知: \(B_0=3.12\times10^{-5}\ \mathrm T\)，\(P_{\rm dyn}=2.00\ \mathrm{nPa}\)。
系统: subsolar stagnation line，忽略 inner pressure、IMF 和 boundary-current geometry factor。
模型: \(B=B_0(R_E/r)^3\)，\(B^2/(2\mu_0)=P_{\rm dyn}\)。
\[\frac{r_m}{R_E}=\left(\frac{B_0^2}{2\mu_0P_{\rm dyn}}\right)^{1/6}=7.60.\]
检查: 六次方根对 pressure 不敏感；实际值的偏差显示理想 dipole balance 的边界，不应隐藏。
Worked Example 4 — Strong-shock compression
已知: \(\gamma=5/3\)，\(v_1=500\ \mathrm{km,s^{-1}}\)，\(M\gg1\)。
参考系: shock rest frame，\(x\) 向 downstream。
模型: ideal nonrelativistic normal hydrodynamic shock。
\[\frac{\rho_2}{\rho_1}=\frac{\gamma+1}{\gamma-1}=4.00, \qquad v_2=v_1\frac{\rho_1}{\rho_2}=125\ \mathrm{km,s^{-1}}.\]
检查: \(\rho v\) 连续；downstream 更密、bulk speed 更低。磁化 oblique shock 还受 \(M_A,\beta\) 和角度影响。
Lesson synthesis / 全课整合
建模流程：声明 object/frame → 比较 \(\lambda_{\rm mfp},\rho_s,d_s,L\) → 算 \(M_s,M_A,\beta,R_m\) → 大尺度用 MHD conservation/topology → shock/reconnection/resonance 切换 kinetic model → 经 radiation/instrument response 生成可观测量 → 带 uncertainty 做 inverse inference。
Formula-and-decision table / 公式与决策表
| Need | Relation | Use only when | Main caution |
|---|---|---|---|
| mass flux | \(4\pi r^2\rho v_r=\dot M\) | steady radial flow | not an acceleration law |
| Parker angle | \(\tan\psi=\Omega r/v_r\) | ideal steady equatorial wind | transients differ |
| dynamic pressure | \(P_{\rm dyn}=\rho v^2\) | stated convention | check factor \(1/2\) |
| stand-off scaling | \(r_m\propto P_{\rm dyn}^{-1/6}\) | dipole pressure balance | geometry/current factors |
| strong-shock ratio | \((\gamma+1)/(\gamma-1)\) | \(M\gg1\), ideal normal shock | kinetic/oblique MHD differs |
| induction | \(\partial_t\mathbf B=\nabla\times(\mathbf u\times\mathbf B)+\eta_m\nabla^2\mathbf B\) | resistive MHD | microphysics hidden |
| \(R_m\) | \(UL/\eta_m\) | characteristic scales defined | large value does not forbid local reconnection |
Model boundaries / 模型边界
- isothermal Parker closure 隐含加热；真实 wind 需 waves、turbulence、conduction 和 multi-fluid effects。
- dipole balance 不描述 reconnection orientation、tail dynamics 和 substorms。
- auroral spectrum 需 atmospheric collisional-radiative model；MHD 单独不足够。
- jump conditions 不解释 collisionless entropy production/nonthermal injection。
- Newtonian disk relation 在 compact-object inner region 失效；jet 可需 GRMHD。
- \(R_m\gg1\) 是大尺度标度，不否定小尺度 nonideal layer。
Misconception clinic / 易错概念诊所
| Misconception | Correction |
|---|---|
| solar wind 是光子风 | 它是 ion-electron plasma bulk flow；radiation pressure 是另一通量 |
| Parker spiral 是粒子螺旋轨道 | 它是 field/flow geometry |
| magnetosphere 完全阻断 wind | shock/reconnection/waves 仍传递质量、动量与能量 |
| aurora 是电子本身发光 | 主要是 atmospheric species 受激发后发光 |
| collisionless shock 没有耗散 | collective fields/waves/scattering 产生熵 |
| dynamo 违反能量守恒 | magnetic energy 来自 rotation/convection/turbulence |
Practice and retrieval / 练习与检索
- [K01] steady wind 远处 \(v_r\) 近似常数，\(\rho(r)\) 如何变？
Answer 1
\(r^2\rho v_r=\mathrm{constant}\)，故 \(\rho\propto r^{-2}\)。
- [K01] 同一 \(r\) 处 wind speed 加倍，Parker angle 怎样变？
Answer 2
\(\tan\psi\) 减半；\(\psi\) 需取 \(\arctan\)，不必精确减半。
- [K02] wind pressure 增为 64 倍，dipole \(r_m\) 变化？
Answer 3
\(64^{-1/6}=1/2\)，故 stand-off distance 减半。
- [K02] 为什么 auroral oxygen 可有不同颜色？
Answer 4
不同 excited transitions 有不同波长，population 还受能量、高度与 quenching 影响。
- [K03] 为什么 disk 角动量须向外输送，质量才能向内？
Answer 5
内移轨道的 specific angular momentum 较小，stress 把多余角动量交给外层。
- [K03] collisionless shock 的“collisionless”指什么？
Answer 6
ramp 远小于 Coulomb mean free path；jump 由 collective electromagnetic processes 维持。
- [K03] \(\gamma=7/5\) 的理想强激波压缩比？
Answer 7
\((\gamma+1)/(\gamma-1)=6\)，仅限 ideal normal hydrodynamic strong shock。
- [K04] \(R_m\ll1\) 和 \(R_m\gg1\) 分别意味什么？
Answer 8
前者 diffusion 占优；后者大尺度 advection/stretching 占优。
- [K04] seed field 由 \(1\ \mathrm{nG}\) 增至 \(1\ \mu\mathrm G\) 需多少 e-folding？
Answer 9
\(N=\ln(10^3)=6.91\)，未计 nonlinear saturation。
- [综合] 为什么 synchrotron brightness 不能唯一给出 \(B\)？
Answer 10
emissivity 同时依赖 relativistic-electron spectrum、\(B_\perp\) 和 path length；需 polarization/Faraday/inverse-Compton 等额外约束。
Exit ticket / 迁移检验
观测到 \(100\ \mathrm{km}\) 宽的 magnetic jump，而 Coulomb mean free path 为 \(10^8\ \mathrm{km}\)。怎样区分宏观 jump 与 ramp 机制？
Exit-ticket benchmark
用 MHD/Rankine–Hugoniot 核查宏观守恒；因 ramp \(\ll\lambda_{\rm mfp}\)，必须检查 waves、reflected particles、cross-shock potential 与 kinetic scales，不能用 binary-collision viscosity。
References and source notes / 参考资料与使用说明
以下来源均在 2026-08-31 实际打开核对；所有 TikZ 依方程原创绘制。
- NASA Science, What Is the Solar Wind? — solar wind、Parker spiral、waves 与 planetary interaction。
- NASA, Aurora Mysteries Unlocked With THEMIS — magnetospheric turbulence、substorms 与 auroral chain。
- NASA Science, Energy Dissipation in Collisionless Shocks — mean-free-path/ramp comparison。
- University of Cambridge Repository, Astrophysical Fluid Dynamics lecture notes — MHD、accretion、jets 与 shocks。
- Max Planck Institute for Radio Astronomy, Cosmic Magnetism — dynamos、galactic fields 与 radio inference。
- NASA Science, Heliosphere — global heliosphere boundary。
Source-to-content map
| Source | Used for | Knowledge ID / Section | Model boundary or convention |
|---|---|---|---|
| NASA Solar Wind | outflow and spiral | K01 | ideal geometry plus observations |
| NASA THEMIS | aurora/substorms | K02 | coupled system |
| NASA Shocks | ramp versus mean free path | K03 | binary collisions insufficient |
| Cambridge AFD | accretion, jets, conservation | K03 | Newtonian/ideal-MHD baseline |
| MPIfR | dynamo and inference | K04 | inverse problem |
| NASA Heliosphere | global boundary | K01 | time-dependent non-spherical bubble |
复习与拓展
- 链接 Black Holes and Accretion 与 Galaxies and AGN，补 relativistic/radiative transfer。
- 链接 Unit 14.3 的 planetary dynamo/Coriolis，比较 terrestrial 与 cosmic dynamos。
- 可扩展：Parker critical integration、diffusive shock acceleration、MRI、\(\alpha\Omega\) dynamo 和 Faraday tomography。