第95章 林叶的办法
「学长，你们为什么非要执着于在每个网格点上都做特征值分解？对于几百万个网格的流场来说，这计算量哪怕不考虑刚性，本身也是天文数字吧？没有想过其他的方法吗？」
林叶来到了张涛的工位旁，拉过一把椅子坐下，便指着关于本征低维流形的那一页，开门见山地问出了自己的想法。
张涛正盯着屏幕上一堆文献发愁，见到是林叶，便揉了揉眉心，苦笑着叹了口气。
「你以为我们想啊？谁都知道算特征值慢，但是没办法。」
张涛放下滑鼠，转过身来，耐心地解释了起来，「我的林学弟啊，你要明白，这是一个强耦合系统。化学反应的非线性极强，那些快变量和慢变量并不是固定的组分。」
「比如在激波后，氧原子可能是快变量；但在壁面附近，它可能就变成了慢变量。它们是随着流场参数动态变化的。」
「如果我们不实时计算雅可比矩阵的特征值和特征向量，就无法准确地定义出局部相空间里的那个流形方向。一旦定义错了，把快变量当成了慢变量处理，或者反过来，整个降维模型就会失效，算出来的结果就是违背物理定律的。」
张涛顿了顿，双手抱住脑袋，靠在了椅背上，接着说道：「我们之前也尝试过其他方法，比如计算奇异摄动法，但那个算法的叠代过程比算特征值还繁琐。相比之下，本征低维流形虽然笨重，虽然算得慢，但至少在理论上它是最严谨、最有希望保证精度的。就像是走迷宫，这条路虽然远，但我们至少也确定它能通向出口，其他的路————可能全是死胡同。」
林叶听完，沉默了。
确实，站在工程实现的角度，为了求稳而选择笨办法是无可厚非的。
如果数学上无法解耦，那就只能用算力去堆。
但是————真的没有别的路了吗？
林叶并没有就这样被说服。
他的直觉告诉他，这种「把流体力学问题完全转化成线性代数问题」的思路，从根子上就有些不对劲。
当然，他只是一个高中生，即使是发表过论文的高中生，也终究是人微言轻，这样的想法他自然没有表示出来。
工位旁陷入了短暂的寂静。
而张涛也以为林叶已经理解了其中的难处，最后勉励了一句：「好了，别想太多，问题难是肯定的，要是啥问题都能够随便解决，咱们地球人怕不是都几级文明了。」
他拍了拍林叶，随后就转过头，准备继续看文献的时候，却突然听到林叶的声音再次响起，带着一丝不确定的试探：「学长，如果————我们不从代数角度去区分快慢，而是从空间角度去区分呢？」
张涛一愣：「空间角度？你是说分区计算？」
「对，但不完全是普通的分区。」
在刚才那段时间的沉默中，林叶的脑海中却以前所未有的效率运转着。