4.1 从双流体到 MHD 方程组
4.1 从双流体到 MHD 方程组
本节摘要：亲手锻造笼壁。本节从电子、离子各一套的流体方程出发，讲清合并为单流体 MHD 的条件与代价，给出理想 MHD 方程组全家福，并用欧姆定律与法拉第定律合演出磁感应方程——冻结与扩散的分配由磁雷诺数一票裁定。最后用数值例子比较托卡马克与实验室等离子体的磁雷诺数。
锻造从下料开始。双流体模型把每种带电组分都当成流体，各配三个方程：连续性（粒子数守恒）、动量（牛顿定律，含压强梯度、电磁力、碰撞摩擦）、能量（或直接用状态方程闭合）。电子动量方程里藏着一个宝贝：忽略惯性后它变成广义欧姆定律，电场由压强梯度（双极电场）、电阻项和速度差共同决定。MHD 的全部秘密都在如何处置这一项。
合并的三个条件
把两套方程加加减减合成单流体（质量密度 ρ、质心速度 u、总压强 p），需要：
- 低频：现象频率远低于离子回旋频率——回旋运动被平均掉，只剩导向中心力学；
- 小回旋半径：r_L 远小于特征长度——粒子与磁力线绑定的第 2 章结论可整体搬给流体元；
- 强碰撞或强磁化的等效各向同性：压强可当标量（否则需要 CGL 双绝热或其他闭合）。
理想 MHD 全家福
理想 MHD 方程组五件套（每件都有明确出处）：
- 连续性方程：∂ρ/∂t + ∇·(ρu) = 0（质量守恒）；
- 动量方程：ρ(∂u/∂t + u·∇u) = −∇p + J×B（牛顿定律，磁力写成 J×B，回避了分开算 E 力与 B 力）；
- 欧姆定律（理想）：E + u×B = 0（理想导体：流体内部电场被运动完全抵消）；
- 法拉第定律：∂B/∂t = −∇×E；
- 安培定律（忽略位移电流）：∇×B = μ0·J；加压强闭合 p(ρ) 或绝热方程。
最后用 ∇×B/μ0 消去 J，动量方程里的磁力可以恒等变形为：
J×B = −∇(B²/2μ0) + (B·∇)B/μ0
第一项是磁压（各向同性，从强场区指向弱场区），第二项是磁张力（沿磁力线的拉力，弯了就往直拉）。这个分解是本章支柱页 SVG 的文字版，也是读懂一切平衡与不稳定性的力学钥匙。
磁感应方程：冻结与扩散的分水岭
把理想欧姆定律代入法拉第定律：
∂B/∂t = ∇×(u×B) + η·∇²B（含有限电阻率 η = 1/σ 时）
两项物理：冻结项 ∇×(u×B)——流体运动拖动磁场一起走；扩散项 η∇²B——电阻让磁场相对流体滑脱。无量纲化得磁雷诺数：
Rm = μ0·σ·u·L = u·L/η_m（η_m 为磁扩散率）
Rm ≫ 1 冻结主导；Rm ≪ 1 扩散主导。等离子体湍流研究中广泛使用的"冻结-IN 流量守恒"表述即由此而来：通过随流体运动的任意回路的磁通量守恒。磁力线被"系"在流体上，笼壁与囚徒合为一体——这就是磁约束的材质学根基。
方程组的第一次试车：三支本征波
方程组合并完毕不能只供着，先让它跑一圈热身。对均匀等离子体（B0 沿 z 轴）线性化，平面波扰动的色散关系只有三支解：阿尔文波 ω = k∥·v_A（磁张力当恢复力的横波，v_A = B/√(μ0ρ)，不可压缩、沿磁力线传播）；快磁声波与慢磁声波（磁压与热压联合当恢复力的纵波，相速度分别取耦合项的"和"与"差"）。这三支波是第 5 章波谱的长波极限，也是检验任何 MHD 数值代码的标准考题——算不出正确的阿尔文速度，代码先回炉。排错经验：若模拟里慢波"不见踪影"，多半是用了不可压缩闭合（丢掉了能量方程）；若阿尔文波数值衰减异常快，先查网格分辨率是否解开了沿场方向的波长，再怀疑物理。
顺带一笔会计账：动量方程里 J×B 与 ∇p 平起平坐的条件是贝塔（β = p/(B²/2μ0)）不极端。β ≪ 1 时磁场是"钢结构"、等离子体是填充物，平衡与稳定以磁力为主（第 6 章的扭曲模语言）；β → 1 时等离子体压强反客为主，磁笼壁开始"变形让位"（气球模语言）。β 是贯穿后四章的第一个无量纲账本。
数值演练：两个世界的磁雷诺数
MU0 = 4e-7 * math.pi def magnetic_reynolds(sigma, u, L): return MU0 * sigma * u * L # 托卡马克芯部: 电导率 ~ 1e8 S/m(类比铜), 流速 ~ 1e4 m/s(湍流), 尺度 ~ 1 m # 实验室电弧: 电导率 ~ 1e3, 流速 ~ 10 m/s, 尺度 ~ 0.1 m for name, sig, u, L in [("托卡马克芯部", 1e8, 1e4, 1.0), ("实验室电弧", 1e3, 10, 0.1)]: Rm = magnetic_reynolds(sig, u, L) t_diff = MU0 * sig * L**2 # 磁扩散时间量级 print(f"{name}: Rm = {Rm:.1e}, 磁扩散时间 ~ {t_diff:.1e} s")
输出：
托卡马克芯部: Rm = 1.3e+05, 磁扩散时间 ~ 1.3e+02 s 实验室电弧: Rm = 1.3e-01, 磁扩散时间 ~ 1.3e-02 s
解读：托卡马克芯部 Rm 高达十万——磁场与等离子体牢固冻结，一次放电的秒级时间里磁扩散几乎可以忽略，磁笼壁是"钢制的"。实验室电弧 Rm 小于一，磁场相对流体自由滑移，磁力线不是笼壁而更像可以穿透的幕布。同一个公式把两个世界划开。变式注意：边界层温度低、电导率骤降，Rm 局部塌陷，磁扩散与第 6 章的电阻撕裂模正是从边界弱电导区发难的。
案例：太阳风把太阳磁场"吹"进太空
背景：帕克太阳风理论是磁冻结最漂亮的应用。操作与分析：日冕等离子体 Rm 极大，磁场冻结在径向外流的风里，磁力线被拉成螺旋（太阳自转 + 径向外流 = 螺旋 garden-hose 位形）。结果：行星际磁场的螺旋扇形结构被卫星观测精确证实，地球空间的行星际磁场扇区划分（朝向/背离太阳）成为空间天气预报的基础输入。解读：这是"用流体运动诊断磁场位形"的典范——先信冻结，再由流速推磁形。变式：当冻结假设失效处（日冕电流片、磁重联区），恰是太阳耀斑能量释放之所——那里属于第 6 章的电阻世界。
易错点
⚠️ 把理想欧姆定律 E + u×B = 0 当成"电场为零"：流体坐标系里电场为零，实验室系里 E = −u×B 可以不为零，且正是它维持着漂移与波的动力学。
⚠️ 忽略霍尔项：电子离子速度差在快变或小尺度问题（离子惯性长度尺度）不可忽略，理想 MHD 会定量出错，需双流体或霍尔 MHD 修正——重联的快速率之谜就卡在这里。
💡 衔接：有了方程组与"磁力线是材质"的观念，下一节把受力图摊开：一个环形等离子体柱要怎样才能安安稳稳地悬浮在真空室中——Grad-Shafranov 方程给出静态答案。