第95章 林叶的办法

  第95章 林叶的办法

  “学长，你们为什么非要执著於在每个网格点上都做特徵值分解？对於几百万个网格的流场来说，这计算量哪怕不考虑刚性，本身也是天文数字吧？没有想过其他的方法吗？”

  林叶来到了张涛的工位旁，拉过一把椅子坐下，便指著关於本徵低维流形的那一页，开门见山地问出了自己的想法。

  张涛正盯著屏幕上一堆文献发愁，见到是林叶，便揉了揉眉心，苦笑著嘆了口气。

  “你以为我们想啊？谁都知道算特徵值慢，但是没办法。”

  张涛放下滑鼠，转过身来，耐心地解释了起来，“我的林学弟啊，你要明白，这是一个强耦合系统。化学反应的非线性极强，那些快变量和慢变量並不是固定的组分。”

  “比如在激波后，氧原子可能是快变量；但在壁面附近，它可能就变成了慢变量。它们是隨著流场参数动態变化的。”

  “如果我们不实时计算雅可比矩阵的特徵值和特徵向量，就无法准確地定义出局部相空间里的那个流形方向。一旦定义错了，把快变量当成了慢变量处理，或者反过来，整个降维模型就会失效，算出来的结果就是违背物理定律的。”

  张涛顿了顿，双手抱住脑袋，靠在了椅背上，接著说道：“我们之前也尝试过其他方法，比如计算奇异摄动法，但那个算法的叠代过程比算特徵值还繁琐。相比之下，本徵低维流形虽然笨重，虽然算得慢，但至少在理论上它是最严谨、最有希望保证精度的。就像是走迷宫，这条路虽然远，但我们至少也確定它能通向出口，其他的路————可能全是死胡同。”

  林叶听完，沉默了。

  確实，站在工程实现的角度，为了求稳而选择笨办法是无可厚非的。

  如果数学上无法解耦，那就只能用算力去堆。

  但是————真的没有別的路了吗？

  林叶並没有就这样被说服。

  他的直觉告诉他，这种“把流体力学问题完全转化成线性代数问题”的思路，从根子上就有些不对劲。