第72章 模型的尸骨（2/2）

江临感觉这片马赛克在嘲笑他。

於是关掉可视化窗口，一头扎进几万行的c++代码里。

开始排查

第一天，他怀疑是数值黏性不够，压不住高频振盪，於是引入了人工耗散项。

结果振盪没了，但重联速率被人工耗散接管，变成了假数据。

第二天，他怀疑是网格尺度太大，局部做了ar。

结果交界面上出现了虚假的反射波。

第三天，他怀疑是边界反射，把流出边界改成了完全匹配层。

第四天到第二十天……

他像个在黑暗中乱撞的困兽，在这三件事之间反覆来回。

每修掉一个问题，另一个更诡异的数值偽影就跳出来。

有时候是质量不守恆了，有时候是磁场的散度不再为零，导致產生了单极磁荷。

第二十一天的凌晨，江临满眼血丝地看著满屏幕被注释掉的测试代码，突然感到一阵无力。

打开工作日誌，敲下了几行字。

【biskap的图不是错的。】

【但它不能脱离那组参数、边界和数值实现，被当作无条件成立的图像。它是在那特定的参数，特定的网格解析度，特定的边界条件和特定的隱式数值耗散共同妥协下，勉强维持出来的一个结果。】

【论文里展示的，是被提纯后打扮得漂漂亮亮的图像。】

【而我现在踩到的，才是这张图像诞生前，那片充满未显式选择、数值妥协和实现细节的泥地。】

写完这一段，江临扫了一眼显示器下方，三年前写下的一张便签。

【数值图像绝对不是物理证据，它只是一台没有感情的机器，负责把你脑子里那些愚蠢的假设画出来而已。】

这张便签歷经三年的风尘，已经泛黄，胶条的边缘打著卷。

江临伸出手指，用力把它压牢。

然后在它的正下方，又贴了一张崭新的便签。

【replicate_biskap_1986：暂缓。】

【原因：不是因为biskap的论文错了，而是因为我还没有足够的能力，去审讯这类复杂的数值结果。】

【下一步：全面补习高阶数值偏微分方程。】

回头看著那长达二十天的排查记录，江临终於意识到了一个更深的问题。

他这二十天失败的，不是排查本身。

而是他根本不知道自己应该排查什么。

他的每一次修改，都是凭直觉和经验去盲猜的。

数值黏性不够

他凭感觉知道。

网格不够密

他凭以前犯过的错知道。

边界反射

他凭听说过的案例知道。

但这些，只是他听说过的最浅薄的几种陷阱。

hd方程组是一个高度非线性的强耦合系统。

他根本不知道在这庞大的参数空间里，还有多少种他连名字都没听说过的陷阱，正安安静静地潜伏在他的代码里，偽造著看似合理的物理现象。

以前学过的《数值分析》，他自学泰勒展开推导差分格式，计算截断误差，用冯诺依曼方法分析一下线性方程的稳定性。

那个版本的知识，只能停留在【知道哪些格式能跑通】的幼儿园阶段。

要真正写出一篇站得住脚的数值模擬论文，或者去甄別別人的模擬是不是垃圾，需要的是【知道每一种格式，在什么样的非线性参数下，会以什么样的方式对你撒谎】的能力。

这一层硬核能力，他没有。

一念及此，他果断关掉所有关於磁重联的论文窗口，新建了一个文件夹。

【nurical_pde_advanced（高阶数值偏微分方程）】

从硬碟调出几本“砖头”。

leveque的《fite vo thods for hyperbolic probles》。

trefethen的《spectral thods atb》。

hairer和wanner的《solvg ordary differential equations ii》。

strikwerda的《fite difference sches and partial differential equations》。

四本神作。

加起来一千八百多页的致密数学。

江临把【replicate_biskap_1986】整个文件夹打了一个压缩包，贴上一个红色的標籤。

【pendg: 等我把这堆数学啃完再回来收拾你。】

然后，他点开了leveque的第一页。

第十四年的整个秋天和冬天，江临完全沉浸在这四本书里。

每天平均学习六到八小时，有时候只能推进十页。

第一遍，拿著草稿纸，把书上的每一个推导亲自推一遍。

第二遍，把书上介绍的每一个数值格式，无论是godunov方法、roe格式还是高阶weno格式，全部自己写一段一维的c++代码，跑出对比图。

第三遍，也是最痛苦的一遍。

构造反例。

他故意设置极端的初始条件，比如强激波或接触间断，去把这些看似完美的格式搞崩，然后分析它们是怎么崩溃的。

简直就像是一个拆弹专家，在学习各种炸弹的起爆原理。

在这期间，他终於明白自己以前对稳定性这个词的理解有多么肤浅。

以前学的冯诺依曼稳定性分析，只能告诉你一个线性方程在无限大均匀网格下的表现。

但现实世界里的hd是非线性的。

他学会了用频谱分析去看待误差，明白了什么是色散关係误差，明白了什么是群速度误差，明白了保结构算法为什么比单纯的精度更重要。

每一种数学工具，都是一台照妖镜。

每一种都能识別出他之前在泥地里瞎撞时完全忽略的失败模式。

第十五年的春天，废土的冰雪开始消融，风力发电机的效率达到了峰值。

江临终於合上了最后一本书。

他没有像个毛头小子一样立刻去解冻biskap的復现项目。

而是打开终端，在工作站里郑重其事地写下一份文档。

【数值健康检查清单v1】

这是他用一年半的时间，用一千八百页的数学底蕴，熬出来的一份审讯大纲。

一共十二项必查项。

他在开头写道。

【对於任何数值模擬结果，以下十二项未经自证，结果一律不予採信，视为偽造物理。】

第一项：线性稳定性与 cfl 余量。

第二项：守恆量检查：质量、动量、能量。

第三项：??b=0约束误差。

第四项：网格解析度与扩散区厚度比例。

第五项：时间步长收敛测试。

第六项：边界反射测试。

第七项：初始扰动频谱检查。

第八项：人工黏性/数值耗散敏感性。

第九项：网格加密收敛。

第十项：改变ndquist数后的拓扑稳健性。

第十一项：与解析標度律的量级比对。

第十二项：可视化结果不得单独作为证据。

这份清单后来被他抄写出来，贴在工作站显示器旁边的石子墙上。

每次跑任何新的模擬之前，他必须像飞行员起飞前一样，逐项打勾。

在跑biskap復现崩溃之前，他从来不知道世界上还需要这样一份清单。

以为只要把公式敲进电脑，电脑就会给他真理。

现在他知道了。

墙上这份薄薄的清单，比那张被顶级期刊发表的，漂漂亮亮的biskap电流片图，更接近真实的物理法则。

第十五年的夏天，江临的目光转向了磁重联领域在二十一世纪初最大的突破。

psoid stability（等离子体团/磁岛不稳定性）。

他调出了两篇里程碑式的文献： loureiro，n.f.，schekochih，a.a.， wley，s.c..stability of current sheets and foration of psoid chas.

bhattacharjee，a.，huang，y..， yang，h.， rogers，b.. fast reion high-ndquist-nuber psas due to the psoid stability.

这两篇论文给出了一条逃离sweet-parker慢结构的重要路线。

它们的核心结论暴烈得就像是一条泛滥的大河。

当一条sweet-parker电流片被拉得足够长、长宽比突破某一个临界值时，它就不再是那条安静死板的麵条了。

它会自己失稳。

从內部断裂，碎成一串像珍珠项炼一样的磁岛。

而每个大磁岛之间，会形成新的更短的微型电流片。

这些微型电流片如果还足够长，会继续碎裂出更小的磁岛。

这是一种分形几何般的联级断裂。

因为每个微型电流片都很短，所以物质排出极快，整体的重联速率直接跳出了慢速深渊，变成了一个与ndquist数几乎无关的快重联。

在相应二维电阻hd参数区间內，整体重联率开始呈现弱依赖甚至近似无关的趋势。

极其优美。

江临满怀期待地开始读loureiro2007年的原稿。

前十六页，物理图像清晰，推导动机明確。

但当他翻到第十七页，看到作者用来证明不稳定性增长率的数学推导时，江临的手再次停住了。

论文动用了一个极其硬核的解析数学工具。

渐近匹配展开。

作者把整个系统强行切成了两半。

內层：电流片中心的极薄区域，这里耗散项主导，採用奇异摄动展开，方程保留了高阶导数。

外层：电流片外部的宏观区域，这里理想hd主导，採用正则展开。

然后，最恐怖的一步来了。

要求內层解的外部极限，必须严丝合缝地等於外层解的內部极限。

两层之间通过一套极其复杂的匹配条件衔接，最终解出本徵值。

江临的眉头拧成了一个疙瘩。

他以前学的《数学物理方法》当然讲过奇异摄动和渐近匹配。

但他清楚地记得，那是宝宝巴士版本。

视频里的教授在黑板上画一个含有极小参数?的常微分方程，解一个最规则的边界层问题。

但loureiro用的这套版本，简直是重武器级別的。

內层是一个非线性的hd偏微分方程组的线性化扰动方程，外层是一个已经被驱动的动態电流片几何背景，两层之间的边界根本不是固定的，而是隨著物理参数动態耦合的。

那个匹配条件本身，就需要证明在这个渐近极限下是良態的。

江临从第十七页开始，读得越来越慢。

读到第二十页的时候，他颓然地放下了笔。

他读不懂。

这不是指他看不懂某个具体的代数变形。

如果只看纸面，论文里的每一行偏导数展开，每一项係数合併，他都能验算得一字不差。

但他没办法判断那一整套浩大的渐近匹配框架，在物理上到底为什么成立。

他无法一眼看透为什么外层磁场的泰勒展开只取到一阶，而內层的速度场却必须保留二阶项

为什么在匹配区，某个特定的流函数分量一定会衰减

他没有【看穿】这个工具的上帝视角。

他只有【机械验算】的奴隶本能。

机械验算，只能保证他这道算术题没算错。

无法保证他真正理解了这个工具的灵魂。

如果稍微改变一下物理背景，换一个方程，他绝不可能自己独立建立起这样一套匹配框架。

江临没有强迫自己继续往下读。

假装看懂是对自己最大的背叛。

他默默地关掉loureiro的论文pdf。

在工作站里，新建文件夹。

【渐近方法深潜】

然后在硬碟里点开了一部史诗级砖头书。

bender，c..， ，s.a.. advaheatical thods for stists and engeers.

一共七百九十三页。

江临打开第一章，从最简单的wkb近似开始，重新走了一遍他以前自以为学过的所有內容。

很快，他出了一身冷汗。

因为他发现自己当年真的只是记住了考点，並没有真正理解。

比如wkb近似。

当年他只会套公式求相函数。

现在，跟著bender和的魔鬼推导，他才真正看懂wkb为什么在转折点附近会失效甚至发散。

原来是因为在那里，波的局部波长趋於无限大，渐近假设崩溃。

所以才必须引入airy函数，在转折点附近建立一个过渡层，把两边的wkb解缝合起来。

比如多尺度展开，为什么不能简单地把快慢尺度叠加

是因为非线性项会產生长期项，导致级数在长时间后爆炸。

必须通过引入慢时间变量来主动消除长期项，这其实是一种高明的物理隔离。

最让他震撼的，是渐近级数与普通收敛级数的根本区別。

收敛级数是项数越多，误差越小，但前几项可能离真值很远。

而渐近级数却是，虽然项数趋於无穷时它会发散到无穷大，但在保留前两三项时，它却能给出极其精確的局部近似。

渐近方法，根本不是追求严密数学收敛的乖宝宝，而是在悬崖边缘跳舞，用野蛮但精准的直觉去逼近真理的刺客工具。

这些东西，他都学过。

但学过和理解，中间隔著十万八千里。

第十五年的整个秋天和冬天，江临像个被困在迷宫里的老鼠，日夜啃噬著这块名叫bender的砖头。

读了和hd强相关的核心章节，奇异摄动，边界层理论，匹配渐近展开。

然后他发现地基不稳，又回过头去死磕主导平衡分析，渐近级数敛散性。

甚至为了理解奇点附近的奇异性，他去研究了斯托克斯现象。

最后，顺手把应用部分也刷了。

延迟微分方程，积分方程的pce方法，非线性振盪的ldstedt-pocaré摄动。

七百九十三页的厚度，几乎每一页平均有一到两个极其烧脑的课后习题。

江临不敢放过任何一个。

第十六年的春天，当废土的第一场春雨敲打玻璃时，江临关掉了这本书。

走到墙边，拿起笔，写下一行他用半条命换来的认知。

渐近匹配展开绝不是单纯的代数数学技巧，它是物理直觉的严格数学化翻译。

如果你的物理直觉错了，怎么匹配都是一团乱麻。

如果你的直觉对了，匹配过程只是把你的直觉翻译成一种可被同行验证的形式语言。

写完这三句话，他感觉眼眶微微发热。

立即回到桌前，重新双击点开了loureiro 2007年的论文pdf。

再次从第十七页开始读。

这一次，他不再是那个在数字里打转的验算工，感觉完全变了。

当他看到loureiro构建內外层匹配边界时，他会在脑海中先停下来，问自己。

如果我是作者，在这里，內层的某一个物理量，在外层应该按什么规律衰减

是代数衰减还是指数衰减”

他先用物理直觉给出答案，然后再去看loureiro的方程。

完全吻合。

他发现，loureiro的每一步看似复杂的代数匹配背后，都对应著一个精准的物理直觉。

因为直觉对了，所以內层展开的高阶项能和外层的极限完美咬合。

因为咬合对了，所以最后解出来的特徵值方程，给出了正確的psoid不稳定性增长率閾值。

江临一口气读完了整篇论文。

通体舒泰，像是在三伏天喝了一桶冰水。

紧接著，他趁热打铁，又去读了bhattacharjee 2009的非线性演化论文。

这一次，他毫无阻碍。

虽然bhattacharjee用的非线性渐近分析更加繁杂，但核心的直觉和工具与loureiro是一脉相承的。

他的视线在公式间游刃有余地穿梭，甚至能看出作者在某一步做截断时的巧妙考量。

当翻到论文最后一页的结论时，一个念头，像一颗种子一样，在他脑海深处破土而出。

他模糊地感觉到自己已经掌握了这套名叫渐近匹配的重武器。

也许，未来有一天，他可以用这套同样的工具，去推导一个和当前文献稍微不一样的物理问题。

比如，加入某种特定的废土环境背景场

这种创造的衝动，是他这十几年钻在故纸堆里前所未有的。

不过江临並没有急於求成，立即去尝试。

他打开手边的笔记本，郑重地记下一行字。

【psoid类磁岛衍生问题，渐近工具已彻底具备。暂不发散，作为核心弹药，储备待用。】