边界条件呢?”他重新戴上眼镜,盯着白板上的?p,“在一个带无滑移边界(no-slip boundary)的弯曲流形上,设定u=0在02上。这个勒贝格积分还能不能把爆破点死死卡在时间 T?”风向变了。费弗曼已经放弃了找茬,开始顺着这个新模型的逻辑,去试探它到底能抗住多极端的物理推演。“当然可以。”林允宁在数位板上快速写下一行边界转换,“通过引入边界层的庞加菜不等式,边界上的涡量生成会被限制在一个紧致子集内。“积分阈值不会改变,依然是唯一的奇点界限。”“那要是流体的初始雷诺数(Reynolds number)奔着无穷大去了呢?”费弗曼重新拿起了笔,“极限状况下的无粘欧拉方程(Euler equations),你这拓扑判据还能起作用么?”“欧拉方程没了VAu耗散项,爆发会提前,但积分发散的本质没变。”林允宁回道,“T会提前,但临界点依旧存在。”话音刚落,右侧那个黑白头像上的红色静音斜杠突然消失。一阵带着高频底噪的呼吸声挤进了会议频道。应用流体力学专家阿瑟·雷诺兹终于出声了:“林先生......”这低沉的中年男声里透着明显的紧绷,“你的意思是,我根本不需要费劲去跑全网格的偏微分方程……………“只要用传感器阵列实时抓取几个关键节点的动能积分,我就能直接锁定系统彻底崩溃的确切时间 T?”“对。”林允宁敲击键盘切回主屏幕,“三次方积分一旦发散,系统就会不可逆转地崩溃。”屏幕右侧标着“Arthur_Reynolds”的黑白头像随之消失,雷诺兹直接切了屏幕共享。林允宁的电脑中央弹出一个界面极简的内部工程播放器,去掉了所有标识,只剩下左上角跳动的绿色时间戳和右边密集的流场遥测参数。“这是03年波音和NASA联合搞的高超音速风洞吹除测试记录。”伴随着底噪,雷诺兹语速极快,“内部编号HX-77,测试目标是马赫数8.5的边界层转捩极限。“这批数据压在机房七年,四个超算团队跑了上千万小时的数值解,也没扯清楚事故的具体物理诱因。”鼠标指针滑到播放键上。“注意音量。”雷诺兹点下播放。音箱纸盆猛地一震,高分贝的低频啸叫冲出屏幕。马赫数8.5的气流被硬生生挤进狭窄风洞,发出刺耳的撕裂声。画面里,一个楔形钛合金模型死死固定在支架上。激波从锐角边缘切过,摩擦高温把周围空气烧出了暗红色的等离子体辉光,驻点温度逼近两千开尔文。左上角的绿色时间戳快速跳动。T = 3.8125。模型表面平稳的红光毫无预兆地开始紊乱,钛合金翼根处绽开一条黑色细线。尖锐的金属断裂声盖过了气流啸叫,失去气动外形引导的高压气流瞬间失控成湍流漩涡。T = 4.120s。测试台在画面中轰然解体,两亿美元的硬件化作高速飞溅的金属破片直接砸碎了摄像头。屏幕切成雪花,啸叫声骤停。扬声器里骤然安静,只剩林允宁电脑底壳风扇的嗡嗡声。“残骸重构显示,彻底解体发生在T=4.125,肉眼可见的气动弹性发散起点是T=3.815。”雷诺兹拖动进度条,把画面卡在金属撕裂前的一帧,顺手拉出一张红蓝相间的时序图,“超算给出的结论是,T=3.8s之前,流体一直处于可控的莱维-霍普夫平滑期。“林先生,能不能用你的拓扑判据来算,临界点到底在哪一秒?”林允宁根本没管那张爆炸截图。“把 T = 0到 T = 3.0之间,模型根部附面层内的原始速度场张量矩阵发给我。”没多久,一个 2.4GB的.csv压缩包从加密通道传了过来。林允宁把它拖进本地计算集群终端,开始敲击代码,将离散的传感器阵列读数重新映射到连续的拓扑流形上:构建网格、导入 u(x,t)矢量场、计算三次方,最后在指定的有界积分域内执行勒贝格积分。import numpynpfro import simps代码写得很快,在敲完最后一行后,他按下回车。CPU利用率瞬间顶满100%。庞大的算力强行咀嚼着七年前的旧数据,终端窗口里的数字成排滚落。四十五秒后,进程结束。屏幕底部拋出了一行红色结果。林允宁拿起电子笔,在那张时序图里画了一条红线。“T = 1.452s。”他语调平淡,随口报出结果,“出问题的是这里。”扬声器那头立刻传来一声不以为然的冷哼。“不可能。”雷诺兹声音高了八度,“T=1.45s时边界层压力波动才0.02%。“跨音速震颤阶段出点局部毛刺再正常不过,当初工程师都认定那是传感器的热背景噪声。整个钛合金结构在那会儿连个划痕都没有。”“从欧几里得几何的视角看,它确实完整。”林允宁把终端的计算过程投到共享屏幕上,“但在拓扑流形上,这结构已经崩了。”他圈出T = 1.452s处那个微小的参数波动。“看这点的勒贝格积分:f_Q|u(x,t)|^3 dx。就在这一微秒内,积分突破紧致子集上限走向发散。“这说明什么?说明涡管根本没在单纯震动,其拓扑结构已经彻底断裂。”画面卡在那个看似毫无异样的瞬间。“能量级联通道发生不可逆的拓扑破缺,系统越过奇异性边界。”林允宁把笔搁回桌上,“所以T=1
本章未完，请点击下一页继续阅读》》