2.3 MHD如何改写流体方程

2.3 磁场入场：MHD 如何改写流体方程

本节摘要：理想磁流体力学把磁场写成被流体带着走的冻结通量，同时把特征速度从声波扩成一簇磁声波。连续方程里散度为零是定理；离散之后它变成必须看守的约束。约束输运把磁通量放在面上、电场放在棱上，才让喷流、吸积盘和恒星发电机模拟有资格自称 MHD，而不是给速度场涂了一层装饰向量。

本节地图

阅读完本节，你应当能够：

1. 说出理想 MHD 相对纯流体增加了哪两类东西：新的波，以及冻结约束
2. 解释散度误差如何变成假力，为何“只在事后投影”常常来不及
3. 描述约束输运的交错配置在守通量上做了什么几何决定
4. 区分理想 MHD 与带耗散、带霍尔项的扩展：何时必须离开理想冻结
5. 判断一套天体物理 MHD 码最可能在激波、接触面还是低等离子体贝塔区翻车

磁场不是给流体上色

阿尔文在 1940 年代给出的波，让等离子体多了一种沿磁力线的张力传播。理想欧拉方程加上感应方程之后，特征速度变成快磁声波、慢磁声波和阿尔文波。CFL 账单立刻改写：时间步可能被阿尔文速度而不是声速钉死。更重要的是几何。理想条件下磁场被冻在流体里，像染料，又像绷紧的绳子。喷流可以被准直，吸积盘可以靠磁旋转不稳定性输运角动量，分子云可以沿力线塌缩得更慢。这些不是后处理着色能假装出来的。

早期有人把感应方程随便差分在格子中心，跑一段时间后散度噪声长成与物理场同量级。噪声进入洛伦兹力，假力把气体推向不该去的地方，看起来像湍流或喷流扭折。共同体花了很长时间才承认：MHD 的入场券不是“我把磁场数组分配了”，而是“我把面通量守住了”。

纯流体已经有的              MHD 新加上的
  质量动量能量守恒           磁通量守恒与冻结
  声波                       快慢磁声波与阿尔文波
  激波俘获                   旋转间断与磁激波
  密度对比                   低贝塔区的力线张力主导
​

💡 关键直觉：磁场把流体从标量压力的世界推进到有方向的张力世界。方向一旦被数值噪声搞乱，张力就变成假的手。

约束输运：把几何写进格子

Evans 与 Hawley 等人推动的约束输运，做了一个看起来像书呆子的决定：磁通量定义在格子面上，电场定义在棱上，用斯托克斯定理更新面通量。这样一来，进入一个体元的磁通量代数相消，散度在机器精度上保持它开始时的值。你仍可能从初值带来散度，但演化不再养出新的。这是铸币时把面额直接打在面上，而不是事后靠警察去抓假币。

其他维持散度的办法包括：投影法每步解泊松方程把非无散部分剥掉；Powell 源项把误差沿流送出域外；双曲清洁引入额外标量把散度当波传走。投影准但贵，还可能破坏其他守恒。源项和清洁便宜，误差是被转移而不是被消灭。天体物理里长时间、强激波、低贝塔的区域，约束输运往往更受信任，不是因为它更神秘，而是它把定理写成了格子的拓扑。

SPH MHD 更别扭。粒子没有天然的面，冻结和散度都要靠核函数的对称去挤。价格方案、欧拉势、约束的粒子表示，各自留下过坑。成功的 SPH MHD 存在，但谱系更窄，测试清单更长。这也是后来许多喷流和盘的工作更愿意住在网格或运动网格上的原因之一。

图：约束输运把通量放在面上

几何决定比系数决定更硬。电场在棱上绕一圈，面通量才知道该增还是该减。

理想冻结何时必须被打破

理想 MHD 假设电阻为零，磁场完全冻结。恒星内部发电机、吸积盘湍流耗散区、磁重联释放能量的位点，都靠非理想项。电阻、霍尔、双极扩散，把冻结撕开一条缝，让力线可以滑过物质。缝的位置若完全交给数值耗散，重联率就变成格式的函数，而不是物理的函数。这是 MHD 方法史上最尴尬的遗产之一：许多“重联驱动的耀发”其实发生在格子尺度上。

对策分成两类。一类把显式非理想项写进方程，并声明耗散尺度。另一类承认未分辨重联，用亚网格或显式开关。无论哪类，论文都需要报告：重联发生在物理电阻尺度还是网格尺度。只展示电流片漂亮的细丝，不报告这一条，等于拒绝让读者审计张力从哪里被剪断。

低等离子体贝塔区还有另一类翻车：磁压远大于气压，微小的相对误差就能把气体抛飞。高马赫数磁激波则考验黎曼求解器是否承认旋转间断。Athena、PLUTO、Athena++ 一脉把这些测试题变成共同体的入场考试。不会做这些题，就不要声称你模拟了活动星系核的喷流。

⚠️ 常见坑：把理想 MHD 跑出的网格尺度重联，写成“我们发现了快速重联机制”。先把电阻尺度和格子尺度的比报出来。

概念账本：

流体守恒更新
交错：面上磁通，棱上电场
用速度和磁场构造电场
绕棱更新面通量
洛伦兹力从面通量重建后加回动量
诊断：最大散度、磁能、阿尔文步长限制
若出现网格尺度电流片：标明是否物理电阻
​

磁场把流体主线改写成带方向的介质。下一节的辐射则把能量方程拖进更高维的相空间：光子不只是“多一个冷却函数”，它们有方向、有频率，也有自己的墙。

从太阳风到吸积盘：MHD 进入天文的几处码头

磁流体进入天体物理，码头不只有喷流。太阳风与日冕要求阿尔文波和湍动级联；分子云要求磁场支持对抗引力；中子星磁层要求力自由场与电流片；行星磁层要求弓激波与重联。每一处码头对非理想项的需求不同。把理想约束输运码原样搬进日冕，会错过电阻与热传导；把日冕的电阻搬进宇宙学盒子，会在不该耗散的地方剪断力线。

测试题清单因此比纯流体更长：阿尔文波传播、圆偏振波、磁转子、奥苏-泰勒激波、电流片的撕裂。过不了磁转子，说明切向磁场在网格上被错误扩散。过不了奥苏-泰勒，说明黎曼求解器不承认旋转间断。这些题看起来像工程，却直接决定你能不能讨论喷流的内磁旋与外磁旋。

低贝塔加上强剪切时，还有数值各向异性热传导和动量输运被格子方向锁定的问题。有人引入显式各向异性传导，让热沿力线走。这是物理，也是对格子各向异性的对抗。对抗本身必须被分辨率扫描，否则你只是用一种各向异性替换另一种。

磁旋转不稳定性：从线性定理到格子上的湍流

磁旋转不稳定性给吸积盘的角动量输运提供了线性起点：弱磁场、开普勒剪切，模态就能长。线性阶段可以用解析增长率审计格式。非线性阶段进入磁湍流，应力变成统计量，分辨率与耗散决定饱和水平。许多盘模拟的“α 粘性”其实是格子湍流的饱和值。把它写进观测解释之前，必须做分辨率扫描，并声明是否有显式物理耗散。没有扫描的 α，是格式常数冒充流体参数。

喷流的磁通量状态——标准与磁能主导——对边界注入的磁通敏感。内边界若把磁通吃掉或反射，喷流可以假死或假活。结合第 3.4 节的视界讨论：牛顿 MHD 盘已经有这张账单，GRMHD 只是把账单写在弯曲面上。码头不同，审计句式相同：通量从哪进、从哪出、散度养没养假力。

问题：约束输运是不是已经把散度问题彻底结束了？

没有。它阻止演化养出新的散度，不清理初值带来的散度，也不自动提供物理电阻。电流片若落在格子尺度，重联率仍是格式函数。约束输运是入场券，入场之后仍要报告最大散度、磁能分辨率扫描、以及重联尺度与格子尺度之比。低贝塔区还可能因重建误差把气体抛飞。把交错网格写成“从此高枕无忧”，会在喷流和盘的非线性阶段再次交学费。

历史上，太阳风与日冕先把阿尔文波变成必须计算的对象，吸积盘的磁旋转不稳定性再把弱场剪切变成湍流实验室，喷流把磁通量状态变成观测可谈的分类。三条码头共用约束输运这张入场券，却不共用电阻与边界。把日冕的耗散搬进宇宙学盒子，或把理想冻结搬进耀发，都会在错误的尺度上剪断力线。读 MHD 应用时，先认码头，再认格式。磁场改写的是波族和拓扑债，不是在动量方程末尾多写一个数组名。约束输运是入场券：它让演化不再养出新的散度，却不提供物理电阻，也不清理肮脏初值。电流片若贴着格子，重联率仍可能是格式在发言。喷流准直、盘湍流、日冕加热，三条应用河共用这张券，却必须各自交代电阻与边界，不能互相借封面。

本章回顾

• MHD 增加波与冻结，不只增加一个数组
• 散度误差会变成假力：看起来像物理，账在格式上
• 约束输运把通量放在面上：代数相消是入场券
• 投影与清洁是转移误差的办法：有用，但要审计去向
• 理想冻结在重联处必须被打破：否则重联率是格子在发言
• 低贝塔与磁激波是试金石：不会测试题就不要画喷流

下一节进入光子：两流近似、矩方法、蒙特卡罗与射线如何轮流接管恒星大气、再电离和吸积盘。