第十三年三月十七日清晨，废土的季风刮过石屋的屋脊，发出沉闷的呜咽。
江临在巡田结束，吃过早餐之后，做了一件酝酿已久的事。
把工作站里所有的二手教材，所有的综述讲义，所有的中文翻译版，所有被前人整出来的精神流食，全部框选。
移到了一个名叫【digested】的归档文件夹里。
然后挂载了一块新硬盘。
硬盘里只有一样东西。
【original_papers】
这是他在现实世界里，从能找到的所有学术网站，学位论文库，影印本档案，扫描期刊里搬运过来的原始论文集。
按年份排序，从一九四〇年代一直到二零二一年。
七千四百二十六篇。
每一篇都是它最初被写出来的样子。
没有教材的简化，没有综述的事后归纳，没有教学逻辑的美化，没有后人为了把它讲清楚而强行抹平的犹豫。
然后将石屋的东墙擦干抹净，写下一行字。
教材是别人替你准备的精神流食，原稿才是先驱者笔尖的犹豫。
写完，他在桌前坐下，像是一个即将踏入未开发原始森林的探险者，点开第一个文件。
Sweet，p.A.(1958).the Neutral point theory of Solar Flares.
中性点理论的太阳耀斑。
一九五八年，伦敦皇家天文学会的一次会议论文集。
扫描件的画质有些糙，英文字母带着打字机特有的毛边。
作为所有磁重联教材的开篇，是基础中的基础，江临之前在教材上看过Sweet-parker模型不下十遍。
对于那个标志性的标度律，他倒背如流。
他能闭着眼睛在画出那条又长又薄，像面条一样的电流片，标出流入区的磁场和流出区的等离子体喷流。
但他未曾知道Sweet本人当年是怎么把它写下来的。
第一页，背景介绍，引入问题。
第二页，磁流体力学方程的初步化简。
第三页，质量守恒与能量守恒的几何假设。
读到第四页时，江临停了下来。
教材上有一句话，几乎所有的磁重联教材都会写。
【在Sweet-parker模型中，电流片的几何被简化为一个长宽比远大于一的薄层。】
简化。
多么轻描淡写的两个字。
任何读过教材的学生都会理所当然地认为，这是作者为了让复杂的偏微分方程在数学上可解，而做出的合理近似。
但Sweet在原稿里，用了整整两页的篇幅来解释这件事。
文字里透着一种面对未知时的局促和坦诚。
他写，他非常清楚这个几何形状在真实的太阳大气里不一定成立。
他知道太阳耀斑爆发的时间尺度只有几十分钟，而这个模型给出的预测时间可能长达几个月甚至几年。
他知道这个模型一旦发表，绝对会被同行批评。
但是他必须这么写。
因为在一九五八年那个没有超级计算机，没有高级数值算法的年代，如果不这么写，那些非线性的mhd方程组根本就解不出来。
在他能找到的所有可解的简化里，这一个是他在数学能力范围内，对真实物理图像破坏最小的那一个。
扫描件底部，Sweet在原稿脚注里写下这样一句话。
【这是一个出发点，而不是一个答案。】
教材从来不会保留这种脚注。
它只会保留那些看似坚不可摧的方程。
所以教材会让所有后来者以为，那个长宽比远大于一的薄层，就是磁重联本来的样子，是某种被确立的真理。
但Sweet知道它不是。
江临于是在东墙上，精神流食那句话下面添了一句。
【Sweet的电流片不是真理，是他能找到的最不糟糕的简化。】
写完，他感觉胸口郁结了很久的一块石头松动了一点。
回到桌前，继续读。
第六页，Sweet终于给出了那个着名的标度律。
教材上，这个公式通常是被放在一个带着灰色阴影的方框里展示的，旁边配上粗体字，仿佛它就是磁重联的固有性质。
但Sweet在原稿里，是这样写下它的。
【在上述简化下，重联速率不可能比1/√S更快。这意味着，对于太阳大气中典型的Lundquist数值，重联将慢到无法解释观测到的耀斑时间尺度。】
紧接着，Sweet又写。
【作者不认为这是磁重联的最终答案。作者只是希望，通过给出这个最简单的模型，能让后续工作有一个明确的批评对象。】
江临轻笑了一下。
教材里写Sweet-parker模型的时候，从来不会保留这种谦卑。
只会用一种高高在上的全知视角，把它当作一个被petschek、hall、plasmoid等后续理论超越的过时模型来批判。
但Sweet本人在写下它的当天，就知道它会被超越。
他甚至在殷切地希望它被超越。
从而主动把自己变成一个靶子，插在未知的荒原上，告诉后来者。
朝我开枪，然后踩着我的尸体往前走。
江临翻到下一篇。
parker，E. N.(1957). Sweets mechanism for merging magnetic fields in conducting fluids.
这是物理学界的巨擘，尤金·帕克。
他在Sweet提出模型的同一年，独立从另一个角度做了相同的推导。
后来的学术界图省事，把这两篇合并称为Sweet-parker模型。
但江临一行行啃完parker的原稿后，后背渗出了一层冷汗。
后世教材为了教学，把parker原稿中更宽的磁雷诺数讨论压缩成了 Sweet-parker 标度律。
其实，parker和Sweet的出发点完全不一样。
Sweet是在解一个具体的几何模型，而parker的野心大得多。
他想从导电流体的磁雷诺数这个更基础的无量纲量出发，推导出在什么物理条件下磁场会被等离子体冻结，而在什么条件下磁场会发生自由扩散。
Sweet-parker标度律，仅仅是parker那长达十几页推导过程中的一个副产品。
parker在原稿里花了大量篇幅去讨论磁雷诺数本身。
甚至在纸面上给出了几个完全不同的物理图像猜想。
每一个图像，都对应着一种磁场演化的极限情形。
Sweet的这套长面条电流片，只是parker备选库里的其中一种可能。
但傲慢的教材大笔一挥，抹掉了parker原本想要建立的那个探讨磁场拓扑本质的图像，只保留了Sweet-parker电流片。
江临再次走向墙壁，再次添上一句。
【教材给你一个公式，原稿给你一整个被遗弃的备选方案库。】
那天傍晚，江临在工作站里建立了一个新的顶层文件夹。
【original_Insights_Lost_In_textbooks（那些迷失在教材中的原初洞见）】
他把Sweet和parker的论文复制了进去，并附上了长达三千字的个人批注。
而这，只是一个开始。
第十三年的整个春天和夏天，江临都在做同一件事。
大量精读，扫读，复现，索引。
他没有按时间顺序顺流而下，而是像个拿着金属探测器的工兵，在重联研究史的雷区里反复跳跃。
读完Sweet和parker，他直接跳到了一九六四年的petschek。
petschek，h.E.(1964). magnetic field annihilation. NASA Special publication Sp-50.
在所有的教科书里，这都是磁重联研究史上最漂亮的一篇论文。
petschek是个天才。
他敏锐地察觉到了Sweet-parker模型中长宽比过大导致物质排不出去的物理瓶颈。
于是大笔一挥，提出了一个全新的几何拓扑。
他不再用那条又长又薄，像面条一样的电流片。
而是引入了一个x型的重联点，并在四周加上了四道向外扩张的慢模激波。
在这个极其优雅的新几何下，大部分等离子体不需要挤进狭窄的中心扩散区，而是直接穿过慢模激波被加热和加速。
磁能释放的速度，理论上可以逼近局地阿尔芬速度。
这是理论上的最快可能。
很好地解决了耀斑爆发的快速时间尺度问题。
教材上写到petschek模型的时候，字里行间总是洋溢着一种英雄史诗般的赞美。
【一九六四年，petschek提出了快磁重联的几何，彻底解决了Sweet-parker模型的慢速困境，开创了快重联的新纪元。】
江临以前也是这么信的。
直到他读到了petschek原稿的最后一节。
视线扫过结论部分时，他整个人猛地在椅子上坐直了。
petschek在原稿里，白纸黑字地写了这样一段话。
【上述模型的成立，完全依赖于在中心扩散区存在某种局域的反常电阻机制。本文无法提供这种机制的微观物理推导。作者坦率承认，这一假设是本模型最为脆弱的部分。如果未来的工作不能在微观动力学尺度上证明这种反常电阻机制的存在，那么本模型只能被视为一种几何拓扑上的可能性，而非一种物理上的必然。】
江临把这段话反反复复读了三遍。
每读一遍，教材的公信力就在心里崩塌一层。
教材永远不会引用这一段，它只会不遗余力地展示petschek那个漂亮绝顶的x型几何，画出那四条威风凛凛的慢模激波。
但petschek自己，在模型诞生之初，就已经把刀架在了自己的脖子上。
他清楚地知道，那四条华丽的激波，那个名叫反常电阻的微观机制，是建立在一个没有被证明的沙盒之上的。
后来的几十年物理史，简直就像是针对petschek这句预言的精准打击。
超级计算机发展起来后，无数的数值模拟反复表明。
在简单的均匀电阻mhd条件下，petschek的x型结构根本站不住脚。
只要稍微跑长一点时间，那个漂亮的x就会坍缩，退化回Sweet-parker那条又长又无聊的面条。
如果不人为地在代码中心区域塞进去一个局域的反常电阻项，petschek的快速重联几何就像是用纸牌搭的城堡，一阵微风就能让它垮塌。
而教材把这段惨烈的历史轻描淡写地写成了，petschek模型在后来的高精度数值模拟中遇到了一些挑战。
挑战？
那更像是物理条件在拒绝这种几何的无条件成立。
petschek本人，在一九六四年就已经预见到了这一切。
他在原稿里明明白白地把自己的软肋暴露给了世界。
江临第三次起身，在东墙的空白处画了一个巨大的x。
在x的左边，写【几何成立】。
在x的右边，写【微观机制未证】。
在x的下方，写【petschek知道，petschek写了，教材没写。】
写完，他热血难凉，走出石屋。
风机二号的尾舵在低风中轻轻摆动。
看上去就像是风机在追着风转。
如果有人现在问他，风机为什么能追着风转？
如果用最简单的教材语言来回答，那就是一句话。
【风机的尾舵面积产生了空气动力学力矩，使得机舱迎风。】
这句话对吗？
对。
它完美地描述了宏观现象。
但这句话有什么用吗？
没用。
因为真正决定风机能不能追风，追得快不快，会不会在湍流中被撕裂的，根本不是那个宏观的尾舵面积。
而是底下那些极其繁琐复杂，毫不起眼的参数的耦合：阻尼弹簧的刚度，偏航轴承的摩擦系数，叶轮旋转时的陀螺推力偏置，控制系统的阵风频率响应……
只要其中一个微观参数崩溃，那个宏观上看起来天经地义的追风动作，立刻就会变成一场灾难。
江临坐在沙地上，脑子里突然闪过petschek论文里的那个x型几何。
下一刻，他冲回石屋。
站在那个大大的x的下方，添上了一句话。
【漂亮的宏观几何背后，永远站着一个没被证明的微观机制。】
读完Sweet、parker、petschek的原稿之后，江临意识到一个严重的问题。
自己一直在纸面上做评判，从来没有真正弄脏过自己的手。
理论物理的尽头，如果不能用实验去验证，就只能靠数值模拟去逼近。
教材上说petschek模型在简单的均匀电阻mhd模拟中退化成了Sweet-parker，这只是一个结论。
江临决定摸一遍这背后的泥。
他要亲手把petschek的模型放进代码里，看着它崩溃。
他选了数值模拟历史上非常经典的一篇论文。
biskamp，d.(1986).magnetic reconnection via current sheets.
这篇论文的核心工作，就是用直接数值模拟研究二维电阻mhd中电流片的演化，并给出了petschek结构无法长期维持的实锤。
江临在工作站里新建了一个项目。
【Replicate_biskamp_1986】
他要用c++，从零开始，手敲每一行求解器代码。
二维笛卡尔网格，均匀电阻假设，反平行磁场初始条件，空间离散采用中心差分，时间推进采用半隐式的预测-校正法。
他试图严格按照论文里的描述去设置Lundquist数、初始扰动幅度、边界条件和网格尺度。
但真正开始写代码时，他立刻发现，老论文里的文字描述，根本不足以完整复现一个数值实验。
论文里写，我们在边界上采用了自由流出条件。
但自由流出是一阶外推还是二阶外推？
有没有做特征波分解处理无反射边界？
如果是简单的零梯度外推，在非线性演化后期绝对会产生数值反射回波，污染中心扩散区。
论文里写，我们引入了一个小幅度的磁通量扰动。
但扰动函数是什么？
是高斯型还是正弦型？
截断范围怎么选？
扰动的高频分量有没有在初始化时做平滑处理？
……
这些细节，在论文里统统没有。
或者说，当年受限于篇幅，作者不可能把几万行代码的细节逐一写进文章里。
但对于复现者来说，每一个缺失的细节，都是一个致命的变量。
江临没有把这些不确定性糊弄过去了事。
他在项目的REAdmE文件的第一行写——
# 非严格复现。尽可能逼近。所有未明实现项，将单独列出参数空间进行扫描。
第一次正式运行，工作站硬盘灯连续亮了十七天。
第十七天下午，终端上的时间步循环终于停下。
江临颤抖着手，用python的matplotlib库调出了最后一个时间步的磁场流函数等值线图。
屏幕上出现了一条缓慢演化的电流片。
他调出时间序列动画。
在模拟的早期，也就是刚加入初始扰动的时候，系统确实在一瞬间出现过petschek那种华丽的x型几何。
那是一个令人心跳加速的瞬间。
但随着时间的推移，中心扩散区无法将堆积的物质及时排出。
那个x的夹角越来越小，重联点被慢慢拉长，变薄。
最后，在动画的结尾，那个x完全坍缩，最后退回那条熟悉得令人不舒服的长电流片。
和biskamp论文里的截图几乎一模一样。
很漂亮。
江临却看得脊背发凉。
太顺了。
他决定做压力测试。
回到配置文件，把Lundquist数提高一个量级。
但这一次，他做了一个违背常规的操作，故意没有同步加密网格。
他想知道，当前这套网格的承受边界到底在哪里。
系统是真的因为物理规律演化成这样，还是因为数值格式把它限制成了这样？
代码重新编译，运行。
这一次，程序跑到第八天深夜，工作站发出一声尖锐的报错蜂鸣。
ERRoR: Floating point exception (core dumped).
跑崩了。
江临立刻调出崩溃前几个时间步的快照。
不是因为更高的Lundquist数更容易算。
而是因为电阻减小后，物理上的电流片扩散区厚度变得极其薄，被压到了接近甚至小于单个网格尺度的地步。
当物理尺度小于网格尺度时，代码里的数值结构就强行接管了物理演化。
江临在屏幕上看到了可怕的一幕。
x型几何刚刚尝试形成，电流片的中央就开始出现非物理的高频振荡。
红蓝相间的色块像病毒一样在网格间蔓延，几个时间步之后，整幅精密的磁场图散成了一片毫无意义的马赛克。
也就是所谓的数值爆炸。
江临感觉这片马赛克在嘲笑他。
于是关掉可视化窗口，一头扎进几万行的c++代码里。
开始排查
第一天，他怀疑是数值黏性不够，压不住高频振荡，于是引入了人工耗散项。
结果振荡没了，但重联速率被人工耗散接管，变成了假数据。
第二天，他怀疑是网格尺度太大，局部做了AmR。
结果交界面上出现了虚假的反射波。
第三天，他怀疑是边界反射，把流出边界改成了完全匹配层。
第四天到第二十天……
他像个在黑暗中乱撞的困兽，在这三件事之间反复来回。
每修掉一个问题，另一个更诡异的数值伪影就跳出来。
有时候是质量不守恒了，有时候是磁场的散度不再为零，导致产生了单极磁荷。
第二十一天的凌晨，江临满眼血丝地看着满屏幕被注释掉的测试代码，突然感到一阵无力。
打开工作日志，敲下了几行字。
【biskamp的图不是错的。】
【但它不能脱离那组参数、边界和数值实现，被当作无条件成立的图像。它是在那特定的参数，特定的网格分辨率，特定的边界条件和特定的隐式数值耗散共同妥协下，勉强维持出来的一个结果。】
【论文里展示的，是被提纯后打扮得漂漂亮亮的图像。】
【而我现在踩到的，才是这张图像诞生前，那片充满未显式选择、数值妥协和实现细节的泥地。】
写完这一段，江临扫了一眼显示器下方，三年前写下的一张便签。
【数值图像绝对不是物理证据，它只是一台没有感情的机器，负责把你脑子里那些愚蠢的假设画出来而已。】
这张便签历经三年的风尘，已经泛黄，胶条的边缘打着卷。
江临伸出手指，用力把它压牢。
然后在它的正下方，又贴了一张崭新的便签。
【Replicate_biskamp_1986：暂缓。】
【原因：不是因为biskamp的论文错了，而是因为我还没有足够的能力，去审讯这类复杂的数值结果。】
【下一步：全面补习高阶数值偏微分方程。】
回头看着那长达二十天的排查记录，江临终于意识到了一个更深的问题。
他这二十天失败的，不是排查本身。
而是他根本不知道自己应该排查什么。
他的每一次修改，都是凭直觉和经验去盲猜的。
数值黏性不够？
他凭感觉知道。
网格不够密？
他凭以前犯过的错知道。
边界反射？
他凭听说过的案例知道。
但这些，只是他听说过的最浅薄的几种陷阱。
mhd方程组是一个高度非线性的强耦合系统。
他根本不知道在这庞大的参数空间里，还有多少种他连名字都没听说过的陷阱，正安安静静地潜伏在他的代码里，伪造着看似合理的物理现象。
以前学过的《数值分析》，他自学泰勒展开推导差分格式，计算截断误差，用冯·诺依曼方法分析一下线性方程的稳定性。
那个版本的知识，只能停留在【知道哪些格式能跑通】的幼儿园阶段。
要真正写出一篇站得住脚的数值模拟论文，或者去甄别别人的模拟是不是垃圾，需要的是【知道每一种格式，在什么样的非线性参数下，会以什么样的方式对你撒谎】的能力。
这一层硬核能力，他没有。
一念及此，他果断关掉所有关于磁重联的论文窗口，新建了一个文件夹。
【Numerical_pdE_Advanced（高阶数值偏微分方程）】
从硬盘调出几本“砖头”。
LeVeque的《Finite Volume methods for hyperbolic problems》。
trefethen的《Spectral methods in mAtLAb》。
hairer和wanner的《Solving ordinary differential Equations II》。
Strikwerda的《Finite difference Schemes and partial differential Equations》。
四本神作。
加起来一千八百多页的致密数学。
江临把【Replicate_biskamp_1986】整个文件夹打了一个压缩包，贴上一个红色的标签。
【pending: 等我把这堆数学啃完再回来收拾你。】
然后，他点开了LeVeque的第一页。
第十四年的整个秋天和冬天，江临完全沉浸在这四本书里。
每天平均学习六到八小时，有时候只能推进十页。
第一遍，拿着草稿纸，把书上的每一个推导亲自推一遍。
第二遍，把书上介绍的每一个数值格式，无论是Godunov方法、Roe格式还是高阶wENo格式，全部自己写一段一维的c++代码，跑出对比图。
第三遍，也是最痛苦的一遍。
构造反例。
他故意设置极端的初始条件，比如强激波或接触间断，去把这些看似完美的格式搞崩，然后分析它们是怎么崩溃的。
简直就像是一个拆弹专家，在学习各种炸弹的起爆原理。
在这期间，他终于明白自己以前对稳定性这个词的理解有多么肤浅。
以前学的冯·诺依曼稳定性分析，只能告诉你一个线性方程在无限大均匀网格下的表现。
但现实世界里的mhd是非线性的。
他学会了用频谱分析去看待误差，明白了什么是色散关系误差，明白了什么是群速度误差，明白了保结构算法为什么比单纯的精度更重要。
每一种数学工具，都是一台照妖镜。
每一种都能识别出他之前在泥地里瞎撞时完全忽略的失败模式。
第十五年的春天，废土的冰雪开始消融，风力发电机的效率达到了峰值。
江临终于合上了最后一本书。
他没有像个毛头小子一样立刻去解冻biskamp的复现项目。
而是打开终端，在工作站里郑重其事地写下一份文档。
【数值健康检查清单V1】
这是他用一年半的时间，用一千八百页的数学底蕴，熬出来的一份审讯大纲。
一共十二项必查项。
他在开头写道。
【对于任何数值模拟结果，以下十二项未经自证，结果一律不予采信，视为伪造物理。】
第一项：线性稳定性与 cFL 余量。
第二项：守恒量检查：质量、动量、能量。
第三项：??b=0约束误差。
第四项：网格分辨率与扩散区厚度比例。
第五项：时间步长收敛测试。
第六项：边界反射测试。
第七项：初始扰动频谱检查。
第八项：人工黏性/数值耗散敏感性。
第九项：网格加密收敛。
第十项：改变Lundquist数后的拓扑稳健性。
第十一项：与解析标度律的量级比对。
第十二项：可视化结果不得单独作为证据。
这份清单后来被他抄写出来，贴在工作站显示器旁边的石子墙上。
每次跑任何新的模拟之前，他必须像飞行员起飞前一样，逐项打勾。
在跑biskamp复现崩溃之前，他从来不知道世界上还需要这样一份清单。
以为只要把公式敲进电脑，电脑就会给他真理。
现在他知道了。
墙上这份薄薄的清单，比那张被顶级期刊发表的，漂漂亮亮的biskamp电流片图，更接近真实的物理法则。
第十五年的夏天，江临的目光转向了磁重联领域在二十一世纪初最大的突破。
plasmoid Instability（等离子体团/磁岛不稳定性）。
他调出了两篇里程碑式的文献： Loureiro，N.F.，Schekochihin，A.A.，& cowley，S.c.(2007).Instability of current sheets and formation of plasmoid chains.
bhattacharjee，A.，huang，Y.m.， Yang，h.，& Rogers，b.(2009). Fast reconnection in high-Lundquist-number plasmas due to the plasmoid instability.
这两篇论文给出了一条逃离Sweet-parker慢结构的重要路线。
它们的核心结论暴烈得就像是一条泛滥的大河。
当一条Sweet-parker电流片被拉得足够长、长宽比突破某一个临界值时，它就不再是那条安静死板的面条了。
它会自己失稳。
从内部断裂，碎成一串像珍珠项链一样的磁岛。
而每个大磁岛之间，会形成新的更短的微型电流片。
这些微型电流片如果还足够长，会继续碎裂出更小的磁岛。
这是一种分形几何般的联级断裂。
因为每个微型电流片都很短，所以物质排出极快，整体的重联速率直接跳出了慢速深渊，变成了一个与Lundquist数几乎无关的快重联。
在相应二维电阻mhd参数区间内，整体重联率开始呈现弱依赖甚至近似无关的趋势。
极其优美。
江临满怀期待地开始读Loureiro2007年的原稿。
前十六页，物理图像清晰，推导动机明确。
但当他翻到第十七页，看到作者用来证明不稳定性增长率的数学推导时，江临的手再次停住了。
论文动用了一个极其硬核的解析数学工具。
渐近匹配展开。
作者把整个系统强行切成了两半。
内层：电流片中心的极薄区域，这里耗散项主导，采用奇异摄动展开，方程保留了高阶导数。
外层：电流片外部的宏观区域，这里理想mhd主导，采用正则展开。
然后，最恐怖的一步来了。
要求内层解的外部极限，必须严丝合缝地等于外层解的内部极限。
两层之间通过一套极其复杂的匹配条件衔接，最终解出本征值。
江临的眉头拧成了一个疙瘩。
他以前学的《数学物理方法》当然讲过奇异摄动和渐近匹配。
但他清楚地记得，那是宝宝巴士版本。
视频里的教授在黑板上画一个含有极小参数?的常微分方程，解一个最规则的边界层问题。
但Loureiro用的这套版本，简直是重武器级别的。
内层是一个非线性的mhd偏微分方程组的线性化扰动方程，外层是一个已经被驱动的动态电流片几何背景，两层之间的边界根本不是固定的，而是随着物理参数动态耦合的。
那个匹配条件本身，就需要证明在这个渐近极限下是良态的。
江临从第十七页开始，读得越来越慢。
读到第二十页的时候，他颓然地放下了笔。
他读不懂。
这不是指他看不懂某个具体的代数变形。
如果只看纸面，论文里的每一行偏导数展开，每一项系数合并，他都能验算得一字不差。
但他没办法判断那一整套浩大的渐近匹配框架，在物理上到底为什么成立。
他无法一眼看透为什么外层磁场的泰勒展开只取到一阶，而内层的速度场却必须保留二阶项？
为什么在匹配区，某个特定的流函数分量一定会衰减？
他没有【看穿】这个工具的上帝视角。
他只有【机械验算】的奴隶本能。
机械验算，只能保证他这道算术题没算错。
无法保证他真正理解了这个工具的灵魂。
如果稍微改变一下物理背景，换一个方程，他绝不可能自己独立建立起这样一套匹配框架。
江临没有强迫自己继续往下读。
假装看懂是对自己最大的背叛。
他默默地关掉Loureiro的论文pdF。
在工作站里，新建文件夹。
【渐近方法深潜】
然后在硬盘里点开了一部史诗级砖头书。
bender，c.m.， & orszag，S.A.(1978). Advanced mathematical methods for Scientists and Engineers.
一共七百九十三页。
江临打开第一章，从最简单的wKb近似开始，重新走了一遍他以前自以为学过的所有内容。
很快，他出了一身冷汗。
因为他发现自己当年真的只是记住了考点，并没有真正理解。
比如wKb近似。
当年他只会套公式求相函数。
现在，跟着bender和orszag的魔鬼推导，他才真正看懂wKb为什么在转折点附近会失效甚至发散。
原来是因为在那里，波的局部波长趋于无限大，渐近假设崩溃。
所以才必须引入Airy函数，在转折点附近建立一个过渡层，把两边的wKb解缝合起来。
比如多尺度展开，为什么不能简单地把快慢尺度叠加？
是因为非线性项会产生长期项，导致级数在长时间后爆炸。
必须通过引入慢时间变量来主动消除长期项，这其实是一种高明的物理隔离。
最让他震撼的，是渐近级数与普通收敛级数的根本区别。
收敛级数是项数越多，误差越小，但前几项可能离真值很远。
而渐近级数却是，虽然项数趋于无穷时它会发散到无穷大，但在保留前两三项时，它却能给出极其精确的局部近似。
渐近方法，根本不是追求严密数学收敛的乖宝宝，而是在悬崖边缘跳舞，用野蛮但精准的直觉去逼近真理的刺客工具。
这些东西，他都学过。
但学过和理解，中间隔着十万八千里。
第十五年的整个秋天和冬天，江临像个被困在迷宫里的老鼠，日夜啃噬着这块名叫bender&orszag的砖头。
读了和mhd强相关的核心章节，奇异摄动，边界层理论，匹配渐近展开。
然后他发现地基不稳，又回过头去死磕主导平衡分析，渐近级数敛散性。
甚至为了理解奇点附近的奇异性，他去研究了斯托克斯现象。
最后，顺手把应用部分也刷了。
延迟微分方程，积分方程的Laplace方法，非线性振荡的Lindstedt-poincaré摄动。
七百九十三页的厚度，几乎每一页平均有一到两个极其烧脑的课后习题。
江临不敢放过任何一个。
第十六年的春天，当废土的第一场春雨敲打玻璃时，江临关掉了这本书。
走到墙边，拿起笔，写下一行他用半条命换来的认知。
渐近匹配展开绝不是单纯的代数数学技巧，它是物理直觉的严格数学化翻译。
如果你的物理直觉错了，怎么匹配都是一团乱麻。
如果你的直觉对了，匹配过程只是把你的直觉翻译成一种可被同行验证的形式语言。
写完这三句话，他感觉眼眶微微发热。
立即回到桌前，重新双击点开了Loureiro 2007年的论文pdF。
再次从第十七页开始读。
这一次，他不再是那个在数字里打转的验算工，感觉完全变了。
当他看到Loureiro构建内外层匹配边界时，他会在脑海中先停下来，问自己。
如果我是作者，在这里，内层的某一个物理量，在外层应该按什么规律衰减？
是代数衰减还是指数衰减？”
他先用物理直觉给出答案，然后再去看Loureiro的方程。
完全吻合。
他发现，Loureiro的每一步看似复杂的代数匹配背后，都对应着一个精准的物理直觉。
因为直觉对了，所以内层展开的高阶项能和外层的极限完美咬合。
因为咬合对了，所以最后解出来的特征值方程，给出了正确的plasmoid不稳定性增长率阈值。
江临一口气读完了整篇论文。
通体舒泰，像是在三伏天喝了一桶冰水。
紧接着，他趁热打铁，又去读了bhattacharjee 2009的非线性演化论文。
这一次，他毫无阻碍。
虽然bhattacharjee用的非线性渐近分析更加繁杂，但核心的直觉和工具与Loureiro是一脉相承的。
他的视线在公式间游刃有余地穿梭，甚至能看出作者在某一步做截断时的巧妙考量。
当翻到论文最后一页的结论时，一个念头，像一颗种子一样，在他脑海深处破土而出。
他模糊地感觉到自己已经掌握了这套名叫渐近匹配的重武器。
也许，未来有一天，他可以用这套同样的工具，去推导一个和当前文献稍微不一样的物理问题。
比如，加入某种特定的废土环境背景场？
这种创造的冲动，是他这十几年钻在故纸堆里前所未有的。
不过江临并没有急于求成，立即去尝试。
他打开手边的笔记本，郑重地记下一行字。
【plasmoid类磁岛衍生问题，渐近工具已彻底具备。暂不发散，作为核心弹药，储备待用。】