第73章 文献的战场(1/2)
第十六年的夏天,江临终於觉得,自己有资格回去面对那片泥地了。
於是打开了尘封已久的【replicate_biskamp_1986】项目。
带著墙上那份严苛的【数值健康检查清单v1】,以及从bender&orszag那里学来的对物理多尺度的敬畏,他把c++代码全部推翻重写。
这一次,他不再盲目地调大参数。
特意把电流片的长宽比,精准地拉到了loureiro理论推导出的那个不稳定性临界閾值之上。
同时,利用高阶数值格式压制了数值耗散,保证网格解析度能够完美解析內层扩散区。
程序跑了整整二十三天。
江临没有像上次那样焦虑,他甚至有心情去检修风机。
第二十三天的傍晚,结果出来了。
江临调出动画。
屏幕上,那条原本像麵条一样的sweet-parker电流片,在演化的中期,中央突然鼓出了一个小小的黑点。
反平行磁场的零点,磁岛的雏形。
江临屏住呼吸,没有动。
紧接著,第二个小黑点出现,第三个, 第五个。
短短十几个时间步之后,整条绷紧的电流片就像一根承受了极限压力的玻璃丝,內部的应力网络全面崩溃。
从內部砰地碎裂开来,变成了一长串大小不一、相吞噬的磁岛结构。
完美的非线性plasmoid不稳定性再现!
江临终於亲眼在自己的代码里看到了原来sweet-parker的长麵条,並不是永远只能绝望地慢下去。
当它慢到极致,被拉伸到物理极限的时候,它会选择玉石俱焚式的碎裂。
整个庞大的僵死结构,在演化中不断细分为无数个微观的快重联区,再不断重组,最终爆发出惊人的能量释放率。
绝境之中,物理规律自己找到了出口。
江临没有被胜利冲昏头脑,立刻转过身,看著墙上的十二项【健康检查清单】。
开始严格执行全套扰动测试。
把网格解析度整体提高一倍,磁岛链依然出现,位置微调,但拓扑不变。
证明收敛。
把lundquist数再提高一个量级,磁岛变得更密集,碎裂更剧烈,符合理论预测。
改变初始扰动的频谱分布,虽然早期的演化路径有差异,但最终的非线性状態殊途同归。
每一项严苛的拷问之后,磁岛串作为一个整体的物理现象,宛如一块经过千锤百炼的钢铁,稳稳地立在那里。
江临打开项目日誌,用严谨的语调写下结案陈词。
【此数值结果在网格加密及参数扰动后仍稳定存在,通过健康清单十二项检验。】
【但请注意:结论仅限於当前的二维mhd参数区间,不可轻易外推至废土天空的高维真实尺度。】
【归档名:mr-plm-04(magnetic reconnection - plasmoid - 04)。】
然后,他走到石屋的北墙。
这面墙是他用来绘製研究史地图的地方。
拿起笔,在plasmoid理论的区块下方,画了一个重重的勾,並打上了一个標记。
【有限可信】
江临看著它,长长地吐出了一口浊气。
第十七年和第十八年,江临的探索触角继续向微观层面延伸,开始啃无碰撞重联和hall效应。
这是单流体mhd理论开始全面崩溃的地带。
当等离子体变得极其稀薄,比如在地球磁层或是废土的高层大气,粒子之间的碰撞几乎消失。
这时候,把等离子体看作一种均匀导电流体的假设就不成立了。
电子和离子因为质量相差悬殊,开始分道扬鑣。
gem magnetic reconnection challenge的论文是绕不过去的一座丰碑。
这是一次罕见的学术界集体大阅兵。
九个顶尖的国际研究组,带著各自引以为傲的九种截然不同的数值代码,去算同一个设定好的重联问题。
论文的结论极其刺眼,就像是在单流体mhd的脸上狠狠扇了一巴掌。
不管你用什么天花乱坠的数值方法,只要你的物理方程里包含了hall效应,所有代码算出的重联速率,最终都会惊人地收敛到一个相近的快重联值。
反之,只要你固执地只用单流体mhd,哪怕你的代码写出花来,结果毫无例外都会悲惨地滑回sweet-parker那条绝望的慢路上。
如果是以前看教材,江临只会看到教科书上总结的一句乾巴巴的话。
【研究表明,hall效应在小尺度扩散区上变得异常重要,是触发快重联的关键机制。】
但现在的江临,已经是对这种教科书体產生免疫抗体的老兵。
他根本不看正文的讚美,直接翻到附录。
像个苛刻的检察官,拿著放大镜去对比那九个研究组设定的网格参数,边界条件,初始扰动类型和数值耗散格式。
看了一周后,他明白了。
gem挑战赛的结果当然极其伟大,它是歷史的转折点。
但它不是终极真理。
它真正证明的是,在那一组被九个组共同约定的特殊的初始扰动和边界条件的笼子里,hall效应是统治者。
如果你换一种边界,或者换一种三维的湍流环境,hall效应还能不能这么立竿见影?
没人敢打包票。
江临走到北墙的研究史地图前,给gem挑战赛画了一个圈。
【强证据,但限定条件非常明確,不可隨意泛化外推。】
这已经是江临能给出的极高评价了。
真正让江临陷入长时间沉默的,是另一篇论文。
burch,j.l.,et al.(2016). electron-scale measurements of magnetic reconnection in space.
发表在《science》上。
这是人类观测史上的奇蹟。
mms卫星编队。
四颗造价高昂的卫星,排成一个紧密的四面体阵列,卫星之间的间距只有区区几公里。
这在广袤的太空中简直就是贴脸。
它们像四根精细的银针,一头扎进了地球的磁层顶。
然后,它们直接在太空中,抓到了那个传说中只存在於纸面和超算里的电子扩散区。
那是一个只有几十公里尺度的微观物理结构。
为了在这转瞬即逝的太空中看它一眼,人类燃烧了价值十亿美元级別的经费,动用了最顶尖的宇航工程技术,才换来这张数据图。
夜里。
江临坐在石屋的桌前,盯住论文里那张带有四面体坐標的相空间分布图。
看了很久,很久。
他的观测点a和b,加上一条移动观测线,组成的最大观测基线不过是可怜的两点五公里。
低照度相机,民用磁力计,风速计,自製时间基准……
这些东西能记录废土天空,能告诉他红带不是错觉,却永远不可能把他送进高空几十公里处那个精微的电子扩散区。
他想在废土上用观测去验证hall电场?
根本不可能。
他想捕捉像mms那样快速重联的实时微观粒子分布?
痴人说梦。
他想像四颗探针那样扎进磁层顶?
这辈子都別想了。
航天工程、探测器阵列和原位观测的技术壁垒,横亘在他和那一层真相之间。
不过江临倒没有像热血漫画男主那样感到屈辱或不甘。
毕竟不是第一次认识边界的残酷性了。
他起身开门,走出石屋。
已经是冬天的深处。废土的天空呈现出病態的灰蓝色。
北方低空那条常年隱隱绰绰的淡红色亮带,今晚並没有出现,似乎被厚重的辐射云层遮挡了。
风机二號在塔尖上缓慢而执著地转动著。
江临双手插在衝锋衣的口袋里,在塔下站了足足半个小时,直到睫毛上结出了白霜。
回到石屋后,他翻开笔记本,在最新的一页写下。
【mms卫星编队,花了我几辈子也拿不到的钱,去看到了一个我用肉眼凡胎永远也看不到的微观世界。】
【我能做的,是从他们已经在山顶確认存在的那个结构出发,往回推。如果那么强烈的电子扩散区真的在头顶的高空发生,它向下辐射的过程中,会不会在我旧相机里,在这2.5公里长的基线上,在这条夜空的红带里,留下某种山脚下的宏观投影?】
【我看不清山顶的雪。】
【但如果山脚下发生了雪崩,我在泥地里测绘雪崩的残骸分布,也是一种测绘。用宏观推演微观,这是我唯一能走的路。】
写完这几段话,他走到北墙前,在burch 2016那篇闪耀的论文旁边,打了一个前所未有的特殊標记。
然后在旁边写下。
【不可触及的高点】
並在下方標註了一行小字。
【这是我用来在黑夜里定向的星,虽然我知道,那绝不是我这具躯壳能到达的地。】
心如止水。
第十九年的春天,废土迎来了一个罕见的平静期,没有沙尘暴,也没有酸雨。
但江临迎来了他研究生涯中的第三次大停顿。
这一次的停顿,比前两次都要绝望。
前两次,无论是c++代码崩溃,还是渐近匹配的代数推导看不穿,好歹都在物理和工程的范畴內。
只要花时间,总能找到抓手。
但这次,他卡在了一篇关於plasmoid不稳定性数学严格性证明的理论文献前。
论文试图从纯数学的偏微分方程分析角度,证明在某些极限定理下,磁岛的破裂是不可避免的。
然后江临发现自己成了一个文盲。
这篇论文里的每一行不等式变形,拆开看他似乎都认识。
但他完全读不懂。
完全无法理解那些符號背后真正承担的,属於现代数学体系的恐怖重量。
文章里充斥著这样的词汇。
sobolev范数估计, 弱解的全局存在性,先验能量不等式,空间正则性的提升……
这些词,江临在自学高等数学和泛函分析时隱约看过。
只勉强能在草稿纸上默写出內积空间或完备性的定义。
可当这些乾瘪的定义,被硬生生地塞进一个极其复杂的非线性mhd方程组里,用来作为砖块,推动一条长达三十页的证明链往前走时,江临第一次清楚地感觉到,他站在一堵高耸入云的光滑铁墙门外。
他能凭藉强大的代数功底,去机械验算论文里的某一个积分不等式是否放缩正確。
但他完全不能判断,作者为什么必须从左边这个空间跳到右边那个空间?
为什么这里要用一个弱收敛,那里又要用一个紧性定理?
这条逻辑的桥,为什么必须这么修?
机械验算,只能证明他这步算术题没算错。
它完全不能证明他理解了这套庞大而精密的高阶数学工具。
那天晚上,江临坐在屏幕前,看著那篇满是积分號的论文,自嘲地摇了摇头。
果断关掉了文献阅读器。
在工作站的根目录下,新建文件夹。
【sobolev_energy_methods】
然后熟练地调出三本在纯数学界大名鼎鼎,足以让无数研究生熬禿头的神书。
evans的《partial differential equations》。
brezis的《functional analysis, sobolev spaces and partial differential equations》。
temam的《navier-stokes equations: theory and numerical analysis》。
三本电子书的图標在桌面上排开,像三块表面长满青苔的拦路石,压在了他通往真理的门槛上。
江临在纸笔记本上,翻开新的一页,重重写下。
【这是我离独立在数学上严格证明一个磁重联领域定理,距离最近的一次。】
【但我可悲地发现,我缺少了一整套用来修桥的现代数学工具。】
【別废话了,先去学怎么造工具。】
接下来的整整一年多时间里,废土石屋里的节奏,变得单调到了近乎自虐的地步。
像一个孤独的数学信徒,一头扎进evans和brezis构建的那个抽象的拓扑世界里。
evans把他按进偏微分方程的泥地里反覆摩擦。
在那里,江临的世界观被重塑了。
以前,方程有解就是算出一个带著x和t的表达式,比如u(x,t)=sin(x?ct)。
物理直觉告诉他,波就是这么传的。
但在evans的世界里,对於非线性方程,有没有解根本不是显式写出一个公式。
那是一条由无数个极其脆弱的逻辑拼凑起来的窄桥。
第一步,先构造一系列平滑的近似解。
第二步,利用极度精妙的先验估计,比如让某项乘以方程本身,然后分部积分,证明这些近似解的某种能量是被控制住的。
第三步,利用这种能量的有界性,引出弱收敛。
第四步,再利用紧性定理,在弱收敛的序列中抽出一个强收敛的子列。
最后一步,让极限量取极限,证明这个影子般的极限,竟然真的满足原方程。
这桥上的每一步,只要有一个不等式放缩漏掉了一丁点误差,整座桥就会轰然坍塌。
没有任何糊弄过去的空间。
江临第一次真正地理解了,为什么一个看似毫无物理意义的数学能量不等式,竟然能拥有如此巨大的权力。
能决定一整族近似解,到底有没有资格活到极限那个彼岸里去。
而卡他卡得最久的,是brezis里的sobolev嵌入定理。
以前,它只是课本上几条排版难看的定理。
现在,它在江临眼里,变成了一整张等级森严的废土边界地图。
“在这个世界里,你想拥有某种特权,就必须用另一种资源来交换。”
江临在手写板上画著图,喃喃自语。
空间维度不同,可微性指数k不同,可积指数p不同,你所能嵌入的空间性质就截然不同。
哪怕你只差了一丁点指標,越过那条线,数学就会毫不留情地甩给你一个发散的反例,告诉你此路不通。
江临是个极其固执的人。
他用了整整两个月的时间,把所有高维低维情况下的常用嵌入定理,全部推导一遍,並列成了一张巨大的表格。
在这个表格里,每一种嵌入定理的后面,他都硬生生地附上了一个自己构造的,或是书上找来的反例。
那张写满反例的巨大表格,后来被他用胶带贴在了石屋的北墙上。
不是为了展示学问。
是为了每天抬头提醒自己一句话。
在非线性方程面前,正则性从来不是你想要就能有的愿望。
它是你必须拿命去换的条件。
brezis补齐了他泛函的底盘。
弱收敛、紧性、对偶空间、自反空间……
这些名词,以前在他脑子里就像是一堆毫无生气的標本名字。
读完brezis后,它们活了过来,变成了他手里锋利的解剖刀。
江临终於能看懂,前人证明文章里那些抽象的词,在整条逻辑链里到底在承受哪一段重量。
强收敛,那是真金白银的极限。
但强收敛太难拿到了。
所以退而求其次,弱收敛只给你一个模糊的影子。
而紧性定理,则是那个变魔术的手,负责在某些特定的低阶空间里,把那些模糊的影子,重新凝聚成可以触摸的实体。
能不能把这种重新凝聚的实体带回原非线性方程里?
那就得看你一开始第一步的先验能量估计有没有做对,有没有漏水。
一环扣一环,美不胜收。
到了temam的纳维-斯托克斯方程。
这已经是流体力学的圣经了。
纳维-斯托克斯方程虽然不是包含了磁场的mhd,但它们拥有极其相似的数学骨架。
不可压条件。
討厌的非线性对流项。
本章未完,点击下一页继续阅读。