2.2 流体SPH与网格分叉

2.2 流体：SPH 与网格方法如何分道扬镳

本节摘要：可压缩天体物理流体在 1970 年代分成两条几乎不来往的河。光滑粒子流体动力学把气体做成带核函数的拉格朗日粒子；欧拉网格把气体做成格子上的守恒量，并用限制器俘获激波。密度对比、接触间断、开尔文-亥姆霍兹不稳定性成为试金石。运动网格是两条河都暴露缺点之后的改道，不是一开始就有的中庸之道。

核心问题

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

1. 说出 SPH 在 1977 年前后被发明时，究竟要解决网格做不到的哪件事
2. 解释守恒形式与黎曼求解器如何让欧拉方法在激波上更老实
3. 用开尔文-亥姆霍兹和接触间断，说明两条路线各自的经典失败模式
4. 区分固定网格、自适应加密、运动网格：谁在动，谁在加密
5. 判断星系形成论文里选择 SPH 或网格，常常是在继承哪一棵代码谱系

1977 年为什么会有人把气体做成粒子

网格流体在工程和武器计算里已经有人工粘性和激波俘获。天文却面临一种网格极不喜欢的对象：高度真空里突然出现几个坍缩核，密度对比可以拉开许多个数量级，边界还在运动。固定笛卡尔格子要么把真空填满浪费，要么在核上太粗。自适应加密当时还没有长成后来的 AMR 工业。Lucy、Gingold 与 Monaghan 把流体写成带质量的粒子，用核函数在邻近粒子上插值密度和压力梯度。粒子跟着流动走，真空不占格子，坍缩自然把粒子聚集到核里。

这是拉格朗日刀对天体物理密度对比的一次准确还击。星形成、星系并合、潮汐尾，SPH 看起来像为它们定制的。代价立刻出现：邻近关系不再由格子拓扑免费提供，必须用核半径去搜邻居；激波要靠人工粘性；接触间断两侧的粒子可以互不渗透，形成表面张力一样的伪力；剪切会被粘性煮熟。网格阵营后来常把这些写成 SPH 的原罪。公平的说法是：SPH 用粒子换来了自适应，就把拓扑和间断处理变成了要单独发明的东西。

网格当时难的                 SPH 用粒子换来的              立刻欠下的
  真空与运动边界             粒子自动去稠密处              邻居搜索
  巨大密度对比               分辨率跟随质量                接触间断伪表面张力
  自由表面                   不需要贴体网格                激波靠人工粘性
​

💡 关键直觉：SPH 不是“不会写网格的人的流体”，它是密度对比先于激波精度时的那条河。

欧拉河：通量、限制器，以及 AMR 何时才长成

欧拉方法把未知量放在格子里，更新写成界面通量。守恒形式保证质量、动量、能量在整块区域上只通过边界进出。Godunov 思路把每个界面看成局部黎曼问题，激波和稀疏波按特征来，熵增落在间断上而不是均匀洒在光滑区。高阶重建必须配合限制器，否则会在间断两侧振荡出负压。这套纪律让喷流、超新星冲击波、吸积激波看起来比早期 SPH 更像教科书里的间断。

固定网格仍怕密度对比。Berger 与 Oliger 的自适应网格细化把细格子只铺在梯度大的地方，粗格子覆盖空洞。FLASH、Enzo、RAMSES 把 AMR 做成天体物理里可继承的基础设施。加密准则变成一种物理口味：按密度加密，星形成核被看见；按激波加密，喷流被看见；按电离锋加密，再电离被看见。错误的准则会让你在漂亮的地方浪费格子，在真正决定反馈的地方仍然粗糙。

两条河的经典对质发生在开尔文-亥姆霍兹不稳定性与星系盘的湍流级联。SPH 若粘性过大，混合被抑制，看起来像“气体不容易搅”。网格若耗散或对齐误差处理不当，切向对流会抹掉接触面。共同体为此打过多年架，有些架是格式的，有些架其实是亚网格冷却和星形成配方不同。把配方差异写成“SPH 与网格谁对”，是方法史上常见的张冠李戴。

图：SPH 与网格对照

不是谁淘汰谁。左边跟着质量走，右边看守通量。当代星系码往往住在两者之间的焊接带。

运动网格：兼得的尝试，以及配方比格式更决定胜负

Springel 的运动网格把格子顶点大致随流，体积变化承担大部分膨胀与剪切，界面仍用黎曼通量处理激波。盘可以少受固定网格的切向耗散，喷流又不必完全回到 SPH 的粘性开关。AREPO 以及后续的网格重构方法，让“SPH 还是网格”这句老口号开始过时。过时不等于消失：实现复杂度、负载均衡、网格畸变后的重连，都是新账单。有人把运动网格写成终极和解，忽略了它仍然要选重建阶数、仍然要亚网格冷却、仍然会在未分辨的星形成区求助配方。

到了星系形成，格式之争常常被星形成与反馈配方盖过。同样的 SPH 或同样的网格，改超新星能量如何注入、动量如何耦合、辐射压如何降阶，星系的质量函数和盘厚度可以换一副面孔。方法史的教训是：先把流体测试题上的差异说清楚，再把配方差异单独记账。混在一张图里争论“谁更像观测”，读者分不清是河道赢了，还是铸币面额表赢了。

⚠️ 常见坑：用开尔文-亥姆霍兹测试否定整条 SPH 谱系，或用一次漂亮的并合图否定整条网格谱系。测试题只否决特定粘性和限制器，否决不了一棵活了四十年的树。

概念上，一次流体更新可以记成：

若 SPH：
    搜邻居，估密度与压力梯度
    开粘性开关于压缩区
    对称更新动量与能量
若网格：
    重建界面左右状态
    限制器防止新极值
    黎曼通量更新格子平均值
若运动网格：
    先移顶点，再算残余通量
无论哪条河：
    检查负密度、熵、角动量
​

流体主线在此分成三岔，磁场将在下一节把第四种纪律——散度为零——叠上来。没有这门纪律，喷流和吸积盘里的“磁”只是装饰用的向量着色。

测试题如何变成两条河的公共法庭

开尔文-亥姆霍兹、瑞利-泰勒、塞多夫、流体静力学大气、开普勒盘，构成流体天体物理的入场考试。SPH 在 KH 上的早期失败被反复引用，后来的开关粘性与人工热传导改善了混合，却不能让核函数变成黎曼求解器。网格在开普勒盘上的切向耗散同样被反复引用，运动网格与高阶重建改善了角动量，却不能让固定笛卡尔格子自动跟着圆轨道走。法庭是公共的，辩护词必须针对具体格式版本，不能针对“SPH 姓”或“网格姓”。

星形成区把考试升级。密度对比、冷却、自引力同时出现，任何一种刀都会在某一项上失分。有人用混合：致密核用粒子，外流用网格，接口再做一次翻译。混合增加了第 2.5 节那种漏油点。是否值得，取决于你的科学问题更怕接触间断还是更怕真空浪费。没有与问题无关的正确刀。

共同体对比项目曾经让多家 SPH 与网格码跑同一星系并合或同一盘。结果通常是：大尺度形态相似，核区星形成率和激波加热不同。相似说明动力学骨架稳健，不同说明条款与粘性仍在发言。读这些对比，应把“相似”和“不同”分开记，不要只截取对自己谱系有利的一半。

自引力流体：刀还在，重力把考试改写了

一旦流体感受到自引力，杰恩斯不稳定性成为新的试金石。解析增长率已知，格式必须在线性阶段贴住它，再谈非线性碎裂。SPH 的核半径与网格的格子尺度都会给出一个最小可塌缩质量。小于它的“星”是分辨率物体。星形成阈值若设在这个最小质量附近，阈值与刀的选择缠死，跨码对比会变成跨分辨率对比。

开普勒盘加自引力后出现旋臂与碎裂判据，对人工粘性和限制器极度敏感。粘性过大，旋臂长不出来；限制器过猛，角动量外运被格式接管，看起来像物理上的吸积。读盘论文，应要求无自引力对照和线性增长率对照。缺这两项，碎裂图只是刀在发言。运动网格在盘上的优势正是把这张发言压低，但它仍要过同一法庭，不能凭谱系豁免。

问题：今天还要不要在 SPH 和网格之间站队？

短答案：站队已经过时，审计没有过时。运动网格与混合方法把河道接在一起，但每一条河留下的失败模式仍在。你要问的不是作者姓哪条河，而是：接触间断、剪切、密度对比、自引力线性增长率这四张卷，这份实现交了哪几张。四张都交了，站队无意义；只交了并合封面图，站队也救不了。跨项目对比时，先统一测试题版本，再谈星系质量函数。否则你在对比粘性开关和加密准则，不在对比宇宙。分叉史留下的不是姓氏，是每条河必须继续偿还的债。

重点提炼

• SPH 为密度对比而发明：粒子跟随质量，欠下邻居、粘性与接触面的债
• 欧拉守恒形式为间断而发明：通量与限制器让激波落在正确位置
• AMR 把固定网格从浪费里救出来：加密准则本身就是物理口味
• KH 与接触间断是试金石：失败模式不同，不要互相写成总判决
• 运动网格是改道不是中庸：兼得的是耗散与俘获，不是免去配方
• 星系形成里配方常比格式更决定终态：先分开记账

下一节让磁场进入同一组流体：冻结定理如何从连续陈述变成格子上必须守住的面通量。