尝试看起来成功了,立刻去找最刁钻的反例去攻击它。
他开始动手。
第一次尝试,他把Lundquist数S、电流片长宽比L/a、离子惯性长度比L/d_i、湍流强度σ、边界开放比η,五个量加权相加。
每个量取对数,得到一个组合数,命名为G_1。
他打开工作站,调出已知样本库,用G_1去拟合那些经过无数次观测的经典事件.
太阳耀斑,地球磁层亚暴,托卡马克装置里的撕裂模破裂,以及大学实验室里小型重联装置的实验数据。
结果出来后,G_1在所有样本上都给出了似乎很接近的数值。
但这恰恰是最大的问题:接近,不等于物理上的一致。
太阳耀斑G_1≈2.3。
磁层亚暴G_1≈2.7。
托卡马克G_1≈1.9。
实验室装置G_1≈1.4。
如果G_1这个参数真的准确描述了系统快慢转变的本质,那么所有快速重联系统的G_1应该都大于某个极其明确的临界值G_c,而慢速系统的G_1应该都小于G_c。
从分类上看,太阳耀斑和磁层亚暴绝对属于快速重联,托卡马克撕裂模勉强算中等,实验室装置绝对是慢速。
它们的G_1数值确实呈现出快>中>慢的趋势。
但这有什么用呢?
临界值G_c又在哪里?
如果G_c=2.0,那么实验室装置(1.4)和托卡马克(1.9)都被分到慢速,磁层亚暴(2.7)和太阳耀斑(2.3)被分到快速。
这种划分在统计学上勉强成立,但在物理学上,江临看着这堆数字,感到很不满意。
不是因为数值分布不够漂亮,而是G_1的定义本身,没有任何坚实的物理根基。
为什么是五个量加权相加,权重的分配凭什么定成这个比例,为什么是取对数而不是取指数或者其他非线性函数?
这些他都答不上来。
G_1看起来更像是一个为了凑出结果而事后诸葛亮般拟合出来的数学工具,而不是一个基于第一性原理先验定义出来的物理量。
事后拟合的工具是很危险的。
就像一个量身定制的谎言,能完美拟合你给它的任何已知旧样本,但一旦遇到从未见过的新样本,它就会立刻崩溃,无法做出任何有效的预测。
江临毫不犹豫地把写满G_1推导过程的十几页手稿全部扔到一边。
他决定摒弃这种已有无量纲数加权拼接的方法,从物理过程本身,从最底层的因果关系开始推导。
风机的每一个通道,都有一个特征驱动时间和一个特征响应时间。
驱动时间——指的是这个通道接收到外部施加的扰动,所需的时间尺度。
响应时间——指的是这个通道内部的机制,对这个扰动做出反馈并产生动作所需的时间。
如果驱动时间远大于响应时间(外界变化慢,自己动作快),这个通道就能游刃有余地主导局面。
如果驱动时间远小于响应时间(外界变化极快,自己动作跟不上),这个通道就会严重滞后,甚至被抛弃。
这正是风机稳定运行的判断标准。
风速变化(驱动,几秒钟)远大于偏航电机的动作时间(响应,几百毫秒),所以风机总能稳稳地迎着风。
但如果废土上刮起诡异的微秒级阵风(驱动时间极小),偏航机构根本来不及响应(响应时间太大),风机就会被动挨打,甚至断轴。
江临意识到,如果把磁重联的通道也这样拆解,每一对通道的驱动/响应时间比,都是一个有着清晰物理意义的无量纲数。
他开始计算。
如果用之前列出的六个通道两两组合比对,根据排列组合,会得到整整十五对驱动/响应比。
理论上,这十五个比值构成的高维矩阵,才是真正完整的耦合网络描述。
它不是某种粗暴的加权,而是十五个相互独立的,有血有肉的物理比值。
但问题随之而来。
十五个数,太多了。
人脑根本无法直接处理十五维参数空间里的相变过程。
这在数学上也许很美,但在工程和预测上,毫无实用价值。
江临想,能不能在这个庞大的十五维矩阵里,做一次主成分分析?
能不能提炼出一两个起着决定性作用的关键比值,作为整个系统快慢转变的咽喉代表?
他重新审视这十五个比值的理论推导。
然后发现其中有些比值,在无论是恒星还是实验室的所有样本里,都无限接近于1。
这意味着这两种机制总是同步的,它们之间的耦合并不构成阻碍系统演化的瓶颈,可以忽略。
有些比值在所有样本里都极大或极小,这意味着某个通道始终处于绝对压制或绝对休眠状态,远离快慢相变的临界点,也可以忽略。
经过近一个月的残酷筛查和计算剥离,只有少数几个比值,在不同的样本环境之间,表现出了剧烈而显着的数量级差异。
也就是说,只有这几个比值,才是真正卡住系统脖子,决定其是在缓慢耗散中死亡,还是在狂暴重联中爆发的关键枢纽。
他小心翼翼地把这三个关键比值挑了出来。
第一个:电流片长度L与plasid次级不稳定性临界激发长度L_c的比值Π_1 = L/L_c。
它决定了系统会不会碎裂出磁岛。
第二个:电流片长度L与离子惯性长度d_i的比值Π_2 = L/d_i。
它决定了系统能不能触碰到无碰撞的霍尔物理.
第三个:等离子体的湍流强度σ与某个临界湍流强度σc的比值Π_3 = σ/σ_c。
它决定了背景噪声是推波助澜,还是喧宾夺主。
三个比值,每一项都有坚如磐石的物理意义,每一项都对应着一个关键通道的驱动/响应匹配生存条件。
江临开始把这三个比值,用一种高度凝练的代数拓扑形式组合起来。
Gate Index,G。
门控指标。
为了处理湍流那种非单调影响,他引入了一个归一化响应函数F(σ/σc)。
这个函数在湍流强度适中时取峰值,在过弱和过强两端同时下降。
G的V0.1形式,被他暂时写成——
【plasid项+hall可达项+湍流响应项】
具体代数形式暂时不进入正文,只进入工作日志。
第一个log在物理上相当于一个数学开关,对应着plasid通道是否被强行开启。
第二个log对应着微观的hall通道在当前尺度下是否可达。
第三个带负号的log,则精确地描述了湍流背景是否已经强到了足以压制其他所有单一重联机制的地步。
公式中出现的负号,以及括号里那个看似奇怪的分式结构x/(1-x)并非江临凭空捏造。
那是因为湍流的影响在物理上是非单调的,适度的湍流会增加耗散,帮助快速重联发生。
这种精妙的非单调代数结构,江临是在回忆Lazarian-Vishniac在1999年发表的那篇经典论文某个不起眼的附录里,见过类似的推导形式。
他从不认为自己是凭空创造的神,自己只是在借用巨人的砖块,并修改了水泥的配方。
G参数的定义终于出炉。
大为兴奋的江临强迫自己冷静下来,进入残酷的交叉检验阶段。
他用新的G参数,重新去拟合那些已知样本。
太阳耀斑:G≈4.1。
磁层亚暴:G≈3.8。
托卡马克撕裂模:G≈1.2。
锯齿崩塌:G≈1.5。
实验室小型重联装置:G≈-0.3。
江临看着这些数字,心跳开始加速。
如果设临界值G_c=2.0。
G > 2.0的,毫无争议地全都是宇宙级别的快速重联(太阳耀斑、磁层亚暴,数值高达4左右)。
G落入[1,2]区间的,全都是处于临界状态、具有阵发性特征的事件(托卡马克的破裂现象)。
G < 1甚至是负数的,则是那些死气沉沉的慢速事件(实验室装置)。
这个分类结果,比之前那个胡乱拼凑的G_1清晰锐利了无数倍。
更重要的是,G的定义有深厚的物理根基。
那三个log项,每一个都对应着具体的通道耦合物理方程的渐近展开,绝不是事后为了好看而凑出来的常数。
看着这组漂亮的拟合数据,江临并没有兴奋地跳起来欢呼。
因为他知道这仅仅是第一步,他还没有去主动寻找反例来攻击它。
如果G这个门控指标真的描述了宇宙中等离子体快慢转变的真理,那它必须且一定要满足以下三个苛刻的条件,缺一不可。
第一,必须能准确预测新样本。
把G应用到他之前定义公式时完全没有见过的新实验或观测数据上,它必须仍然能给出正确的物理分类。
不能见光死。
第二,必须对参数变化具备物理敏感性。
如果某个系统的外部参数被人为地从慢速区间调整到了快速区间,G的数值必须在数学上相应地从小变大,且变化趋势必须符合微积分的单调性预期。
第三,必须与已有的经典观测现象兼容。
他开始动手。
第一次尝试,他把Lundquist数S、电流片长宽比L/a、离子惯性长度比L/d_i、湍流强度σ、边界开放比η,五个量加权相加。
每个量取对数,得到一个组合数,命名为G_1。
他打开工作站,调出已知样本库,用G_1去拟合那些经过无数次观测的经典事件.
太阳耀斑,地球磁层亚暴,托卡马克装置里的撕裂模破裂,以及大学实验室里小型重联装置的实验数据。
结果出来后,G_1在所有样本上都给出了似乎很接近的数值。
但这恰恰是最大的问题:接近,不等于物理上的一致。
太阳耀斑G_1≈2.3。
磁层亚暴G_1≈2.7。
托卡马克G_1≈1.9。
实验室装置G_1≈1.4。
如果G_1这个参数真的准确描述了系统快慢转变的本质,那么所有快速重联系统的G_1应该都大于某个极其明确的临界值G_c,而慢速系统的G_1应该都小于G_c。
从分类上看,太阳耀斑和磁层亚暴绝对属于快速重联,托卡马克撕裂模勉强算中等,实验室装置绝对是慢速。
它们的G_1数值确实呈现出快>中>慢的趋势。
但这有什么用呢?
临界值G_c又在哪里?
如果G_c=2.0,那么实验室装置(1.4)和托卡马克(1.9)都被分到慢速,磁层亚暴(2.7)和太阳耀斑(2.3)被分到快速。
这种划分在统计学上勉强成立,但在物理学上,江临看着这堆数字,感到很不满意。
不是因为数值分布不够漂亮,而是G_1的定义本身,没有任何坚实的物理根基。
为什么是五个量加权相加,权重的分配凭什么定成这个比例,为什么是取对数而不是取指数或者其他非线性函数?
这些他都答不上来。
G_1看起来更像是一个为了凑出结果而事后诸葛亮般拟合出来的数学工具,而不是一个基于第一性原理先验定义出来的物理量。
事后拟合的工具是很危险的。
就像一个量身定制的谎言,能完美拟合你给它的任何已知旧样本,但一旦遇到从未见过的新样本,它就会立刻崩溃,无法做出任何有效的预测。
江临毫不犹豫地把写满G_1推导过程的十几页手稿全部扔到一边。
他决定摒弃这种已有无量纲数加权拼接的方法,从物理过程本身,从最底层的因果关系开始推导。
风机的每一个通道,都有一个特征驱动时间和一个特征响应时间。
驱动时间——指的是这个通道接收到外部施加的扰动,所需的时间尺度。
响应时间——指的是这个通道内部的机制,对这个扰动做出反馈并产生动作所需的时间。
如果驱动时间远大于响应时间(外界变化慢,自己动作快),这个通道就能游刃有余地主导局面。
如果驱动时间远小于响应时间(外界变化极快,自己动作跟不上),这个通道就会严重滞后,甚至被抛弃。
这正是风机稳定运行的判断标准。
风速变化(驱动,几秒钟)远大于偏航电机的动作时间(响应,几百毫秒),所以风机总能稳稳地迎着风。
但如果废土上刮起诡异的微秒级阵风(驱动时间极小),偏航机构根本来不及响应(响应时间太大),风机就会被动挨打,甚至断轴。
江临意识到,如果把磁重联的通道也这样拆解,每一对通道的驱动/响应时间比,都是一个有着清晰物理意义的无量纲数。
他开始计算。
如果用之前列出的六个通道两两组合比对,根据排列组合,会得到整整十五对驱动/响应比。
理论上,这十五个比值构成的高维矩阵,才是真正完整的耦合网络描述。
它不是某种粗暴的加权,而是十五个相互独立的,有血有肉的物理比值。
但问题随之而来。
十五个数,太多了。
人脑根本无法直接处理十五维参数空间里的相变过程。
这在数学上也许很美,但在工程和预测上,毫无实用价值。
江临想,能不能在这个庞大的十五维矩阵里,做一次主成分分析?
能不能提炼出一两个起着决定性作用的关键比值,作为整个系统快慢转变的咽喉代表?
他重新审视这十五个比值的理论推导。
然后发现其中有些比值,在无论是恒星还是实验室的所有样本里,都无限接近于1。
这意味着这两种机制总是同步的,它们之间的耦合并不构成阻碍系统演化的瓶颈,可以忽略。
有些比值在所有样本里都极大或极小,这意味着某个通道始终处于绝对压制或绝对休眠状态,远离快慢相变的临界点,也可以忽略。
经过近一个月的残酷筛查和计算剥离,只有少数几个比值,在不同的样本环境之间,表现出了剧烈而显着的数量级差异。
也就是说,只有这几个比值,才是真正卡住系统脖子,决定其是在缓慢耗散中死亡,还是在狂暴重联中爆发的关键枢纽。
他小心翼翼地把这三个关键比值挑了出来。
第一个:电流片长度L与plasid次级不稳定性临界激发长度L_c的比值Π_1 = L/L_c。
它决定了系统会不会碎裂出磁岛。
第二个:电流片长度L与离子惯性长度d_i的比值Π_2 = L/d_i。
它决定了系统能不能触碰到无碰撞的霍尔物理.
第三个:等离子体的湍流强度σ与某个临界湍流强度σc的比值Π_3 = σ/σ_c。
它决定了背景噪声是推波助澜,还是喧宾夺主。
三个比值,每一项都有坚如磐石的物理意义,每一项都对应着一个关键通道的驱动/响应匹配生存条件。
江临开始把这三个比值,用一种高度凝练的代数拓扑形式组合起来。
Gate Index,G。
门控指标。
为了处理湍流那种非单调影响,他引入了一个归一化响应函数F(σ/σc)。
这个函数在湍流强度适中时取峰值,在过弱和过强两端同时下降。
G的V0.1形式,被他暂时写成——
【plasid项+hall可达项+湍流响应项】
具体代数形式暂时不进入正文,只进入工作日志。
第一个log在物理上相当于一个数学开关,对应着plasid通道是否被强行开启。
第二个log对应着微观的hall通道在当前尺度下是否可达。
第三个带负号的log,则精确地描述了湍流背景是否已经强到了足以压制其他所有单一重联机制的地步。
公式中出现的负号,以及括号里那个看似奇怪的分式结构x/(1-x)并非江临凭空捏造。
那是因为湍流的影响在物理上是非单调的,适度的湍流会增加耗散,帮助快速重联发生。
这种精妙的非单调代数结构,江临是在回忆Lazarian-Vishniac在1999年发表的那篇经典论文某个不起眼的附录里,见过类似的推导形式。
他从不认为自己是凭空创造的神,自己只是在借用巨人的砖块,并修改了水泥的配方。
G参数的定义终于出炉。
大为兴奋的江临强迫自己冷静下来,进入残酷的交叉检验阶段。
他用新的G参数,重新去拟合那些已知样本。
太阳耀斑:G≈4.1。
磁层亚暴:G≈3.8。
托卡马克撕裂模:G≈1.2。
锯齿崩塌:G≈1.5。
实验室小型重联装置:G≈-0.3。
江临看着这些数字,心跳开始加速。
如果设临界值G_c=2.0。
G > 2.0的,毫无争议地全都是宇宙级别的快速重联(太阳耀斑、磁层亚暴,数值高达4左右)。
G落入[1,2]区间的,全都是处于临界状态、具有阵发性特征的事件(托卡马克的破裂现象)。
G < 1甚至是负数的,则是那些死气沉沉的慢速事件(实验室装置)。
这个分类结果,比之前那个胡乱拼凑的G_1清晰锐利了无数倍。
更重要的是,G的定义有深厚的物理根基。
那三个log项,每一个都对应着具体的通道耦合物理方程的渐近展开,绝不是事后为了好看而凑出来的常数。
看着这组漂亮的拟合数据,江临并没有兴奋地跳起来欢呼。
因为他知道这仅仅是第一步,他还没有去主动寻找反例来攻击它。
如果G这个门控指标真的描述了宇宙中等离子体快慢转变的真理,那它必须且一定要满足以下三个苛刻的条件,缺一不可。
第一,必须能准确预测新样本。
把G应用到他之前定义公式时完全没有见过的新实验或观测数据上,它必须仍然能给出正确的物理分类。
不能见光死。
第二,必须对参数变化具备物理敏感性。
如果某个系统的外部参数被人为地从慢速区间调整到了快速区间,G的数值必须在数学上相应地从小变大,且变化趋势必须符合微积分的单调性预期。
第三,必须与已有的经典观测现象兼容。