2.5 多物理耦合的拼接史

2.5 多物理耦合：从模块拼接到一体求解

本节摘要：真实天体很少只服从一种方程。重力拉、流体推、冷却漏能量、辐射改电离、磁场给张力、核反应改成分。把它们联立成单一巨型隐式系统，至今仍非常少见。历史选择是算子分裂：每个子步只跑一个模块，用对方刚留下的状态当源。模块因此能活，误差也因此在接口安家。理解耦合，就是审计这些接口，而不是收集更多物理名词。

本节导读

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

1. 用算子分裂说明“包含多种物理”如何在时间上被切开
2. 指出重力与流体、冷却与流体、辐射与流体三类接口最常见的漏油方式
3. 解释刚度来自时标分离，而不是来自代码行数
4. 区分单向耦合、双向耦合、亚网格源项三种深度
5. 读耦合流程图时，能问出：哪个子步用的是过时的对方状态

为什么几乎总是先切开再粘

联立求解意味着一个时间步里，重力势、流体通量、辐射矩、感应方程互相看见彼此的最新值。矩阵巨大，非线性，还可能破坏各模块已经调好的守恒技巧。星团纯 N 体可以不切；一旦气体出现，冷却时标在稠密核里可以短到让流体 CFL 变成笑话。再叠辐射，光时标更短。共同体选择了工程师更熟悉的办法：把一步拆成重力、流体、源项、辐射等子步，像工厂流水线。Godunov 分裂、Strang 分裂，决定源项在流体前后各做一次还是只做一次。偶数阶分裂能把某些交叉误差压低，不能消灭“用过时状态”这件事本身。

切开的好处是谱系可以独立进化。树码组改开角，不必重写黎曼求解器；辐射组改 M1 闭合，不必动约束输运。开源码能活二十年，多半因为接口比完美联立更允许人活。代价是物理上本该同时发生的事，被排成了队列。队列顺序会改变激波加热与冷却的竞速，从而改变能否形成稠密核。

本该同时发生的                 分裂后的队列
  压缩加热与辐射冷却           先流体压缩，再冷却，或反过来
  电离与复合                   先更新辐射，再用旧密度复合
  磁张力与气体压力             流体与感应交错半步
  自引力与坍缩                 势场用的是上一流体态
​

💡 关键直觉：耦合图上的箭头方向，就是误差允许流动的方向。没有箭头的“物理”，只是注释。

三类最爱漏油的接口

重力与流体。粒子或网格把质量交给泊松求解器，势再以加速度或引力源项回到动量。若质量在流体步里被通量改写，势却要等下一步才更新，自由落体在细格子上会相位滞后。个别时间步让问题更尖：有的区域走了好几步，势还是老的。宇宙学码用中点势或预测质量来减轻，但不能假装重力与流体在同一瞬间握手。

冷却与流体。冷却率对温度和密度极度敏感，显式冷却可以把内能扣成负数。历史对策是对冷却取隐式或半隐式，甚至限制单步最多扣掉若干成内能。限制本身是亚网格：它在说“未分辨的微观过程不允许温度跳进我的格式不能处理的谷”。星形成配方往往就接在这个谷上——密度高、温度低到某阈值，就变成恒星粒子。于是冷却接口的数值限制，直接参与了恒星的诞生。把星形成写成纯物理阈值，而不提冷却步的限制器，是接口漏油被镀成理论。

辐射与流体。加热、电离、辐射压必须写回气体。若辐射用的是上一步的密度，电离锋会滞后或超前一个格子。光子守恒若只在辐射模块内成立，流体平流把电离物质带走时又可能破坏预算。双向耦合意味着辐射看见新密度，流体看见新辐射力，通常要子循环。单向耦合——辐射照一次，流体自己走很久——只适合光源几乎不变的场面。

刚度、亚网格源项，以及“一体求解”仍有多远

刚度来自时标分离。核反应、冷却、化学反应网可以把常微分方程变得极硬，流体自己还在按 CFL 走。把硬源项丢进显式流体，是在邀请爆炸。方法史的回答是：硬的部分交给隐式常微分方程求解器，软的部分留给双曲格式。化学网络与核网络因此常常像挂在流体身上的另一台小机器。小机器的误差可以主导元素丰度，却几乎不出现在你看的密度图上。

亚网格源项是更彻底的承认：有些过程永远不会被格子看见，超新星、H II 区、黑洞吸积盘的内区。它们以质量、动量、能量、光子的源项形式注入。注入方式——热注入、动能注入、延迟冷却、随机位置——每一项都改变星系尺度的结果。第 5 章会专门讲亚网格如何成为当代拐点。此处只需把账记在耦合史上：源项是接口，不是装饰性物理清单。

真正的一体求解并非没有人做。隐式辐射流体、全相对论的 GRMHD 辐射，在致密天体和少数恒星问题上推进过。代价是线性代数、内存、以及难以把宇宙学盒子装进去。当代大型星系形成模拟仍然是分裂的联邦：N 体重力、流体、冷却库、星形成模块、黑洞模块、偶尔加上辐射。联邦能跑完哈勃时间，一体求解还在特区里。

⚠️ 常见坑：在方法表里列出六种物理，就声称“自洽耦合”。先看时间步流程图里有没有回流迭代。没有回流，就是单向通知。

概念账本：

一个大步：
    更新重力势
    流体双曲子步
    隐式冷却与化学
    辐射隐式或稳态
    磁通约束更新
    亚网格源项注入
    若紧耦合：重复辐射与流体直到残差下降
诊断：各模块守恒、接口处负密度、锋面位置相对解析极限
​

第 2 章的零件到此收齐。第 3 章把它们装进对象：恒星会优先要隐式结构和核网；星系要联邦；宇宙学要先长暗物质骨架再挂气体；黑洞要 GRMHD；行星要尘埃与 N 体并存。对象选择零件，不是零件平等地服务所有对象。

联邦制的治理：谁有权改接口

大型开源码把模块分成基础设施与物理包。改黎曼求解器要过测试题库；改星形成阈值往往只需改参数文件。这种治理让科学用户能活，也让接口误差被参数化得过于容易。一个延迟冷却时间可以在不碰任何通量公式的情况下改写星系风。审稿若只盯流体格式阶数，会漏掉真正移动结果的那一条款。

紧耦合迭代出现在辐射流体和化学反应里：同一大步内反复交换密度与辐射场，直到残差下降。迭代次数成为隐式的分辨率——迭代不够，锋面仍用过时状态。把最大迭代次数写进方法节，和写 CFL 一样必要。有的运行在稀薄区很快收敛，在电离锋上达到上限仍未收敛，日志里应能看见。看不见日志的可视化，无法证明握手完成。

还有一种政治性的耦合：两套码松耦合，例如一维恒星演化定期把剖面交给三维对流盒，再把混合效率送回。松耦合承认时标分离，也承认反馈可能过时。它比假装一体求解更诚实，只要送回的量被当作模型而不是当作测量。

算子顺序作为隐藏参数：冷却与激波的竞速实验

把同一初值跑两遍，一遍先流体后冷却，一遍先冷却后流体，核区密度可以给出可观测的差别。这不是新物理，是分裂顺序。Strang 对称分裂能减轻但不能消灭对强冷却的敏感。诚实的方法节应报告顺序，并至少做一次对调实验。若对调后科学结论翻转，接口比模块更决定胜负，应把计算资源用于紧耦合，而不是用于更高阶的无源流体。

化学网络的隐式求解器有自己的时间步与容差。容差过松，丰度振荡；过紧，化学比流体还贵。网络大小——从氢氦到包含金属线——会改冷却曲线，从而改星形成。把网络版本钉死，和把流体格式钉死属于同一纪律。联邦治理若允许物理包随意升级网络而不跑回归，产额表会在沉默中漂移。

问题：列出六种物理是否等于自洽耦合？

不等于。看时间步流程图有没有回流。单向通知——辐射照一次、流体走很久——只适合光源几乎不变的场面。冷却与激波竞速时，算子顺序本身是隐藏参数，应对调实验。化学网络的容差、最大迭代次数、延迟冷却时间，都能在不碰黎曼求解器的情况下改写星系风。方法表里的物理名词是菜单；箭头、容差和顺序才是治理。没有回流就不要写自洽。

联邦制能活二十年，是因为接口允许模块单独进化；也正因为如此，延迟冷却时间这类参数可以在不碰通量公式的情况下改写风。审稿若只盯格式阶数，会把真正移动结果的治理漏掉。把顺序、容差、最大迭代写进方法节，和写 CFL 一样基本。对调实验便宜，却能揭示接口是否比模块更决定胜负。

温故知新

• 算子分裂让模块能活：也让误差住在接口
• 箭头方向就是误差方向：清单上的物理没有箭头就不算耦合
• 冷却限制器会参与造星：数值谷可以冒充物理阈值
• 刚度要把硬源项交给隐式 ODE：不要丢进显式通量
• 亚网格源项是联邦的宪法：注入方式属于方法，不只属于“模型”
• 一体求解仍在特区：盒子宇宙学靠分裂联邦跑完时间

下一章按对象重走一遍历史：从一维恒星演化到三维爆发，从无碰撞晕到重子星系，从粒子网格宇宙到 GRMHD 视界。