征低维流形ild思路，试图通过特征值分解来区分快慢变量。但这涉及到了繁重的矩阵运算，时间复杂度高达0（n3），而且在刚性极强时极不稳定。」

林叶切换ppt，屏幕上出现了一张高超声速流场的示意图，激波层和边界层被标红。

「我的想法是，跳出代数思维，回归流体力学的物理本质。化学反应虽然快，但剧烈的非平衡过程只发生在激波后和壁面附近的反应层里，在流场其余90的区域，反应其实处于局部平衡态。」

「因此，我构造了一个基于组分梯度和局部达姆科勒数的渐近判据。」

林叶在白板上刷刷写下了一行不等式，那是他这两天根据自己的思路，最终优化出来的判据公式。

「利用这个判据，我们可以动态地将网格标记为刚性区和平衡区。在刚性区，我们继续使用隐式算法；而在平衡区，我们直接使用显式算法或代数平衡关系。」

说到这里，他的声音顿了顿，转头看向了在场的众人。

而此时此刻，大概除了已经和他了解过思路的周文渊和张涛之外，其他所有人都已经是一副看呆了的模样，看着他，以及ppt和他在白板上面写下的那行不等式。

他们————没有看错吧？

这些东西，都是这个高中生搞出来的？！

除了李建国和王海峰，剩下三个团队的人都是满脸的不可思议。

想想他们最开始还满是对林叶的不相信，此时只感觉脸有些疼。

至于李建国和王海峰，此时也沉浸在一阵震撼当中。

这个年轻人，可真是————吓了他们一跳。

这个思路，他们此前竟然从来没有想到过。

但两位教授的震撼，持续时间也很短暂。

不多时，作为计算专家的王海峰教授，就直接看出来了这个方案的致命问题。

他坐直了身体，双手交叉放在了桌上。

「想法不错，林同学。」王海峰开口了，语气虽然客气，但问题直指核心，「但是，你怎么处理这两个区域的交界面？显式算法和隐式算法的通量计算方式完全不同，如果在界面处直接拼接，必然会产生虚假的数值反射和振荡。这对于高超声速流动来说，任何一点虚假波动都会被激波放大，导致计算发散。」

「王教授说得对。」林叶早有准备，他自信地切换到下一页ppt，「硬拼接肯定不行。所以，我准备在两个区域之间，引入一个缓冲过渡层。」

「在这个过渡层里，不做那些非黑即白的切换，而是构造一个基于sigoid函数的权重因子α，然后在过渡层内求解一个混合通量，大概长这样——」

【f—hybr id=af—iplicit+（1—a）f—eplicit】

「同时，」林叶指着公式中的一项修正项，「为了保证守恒性，我利用了梅林变换技巧，从数学上推导了一个通量修正项。理论上，这能确保在过渡层内的质量和能量守恒误差控制在二阶精度以内。」

梅林变换！

看着林叶给出的这一步，下面的张涛也露出微微感慨的神情。

他和林叶前几天用到的方法，终究还是发挥了作用。

章节内容不完整，请退出阅读模式查看完整内容！