如遇到章节错误,请关闭浏览器的阅读/畅读/小说模式并且关闭广告屏蔽过滤功能,稍后尝试刷新。
他们————没有看错吧?
这些东西,都是这个高中生搞出来的?!
除了李建国和王海峰,剩下三个团队的人都是满脸的不可思议。
想想他们最开始还满是对林叶的不相信,此时只感觉脸有些疼。
至於李建国和王海峰,此时也沉浸在一阵震撼当中。
这个年轻人,可真是————嚇了他们一跳。
这个思路,他们此前竟然从来没有想到过。
但两位教授的震撼,持续时间也很短暂。
不多时,作为计算专家的王海峰教授,就直接看出来了这个方案的致命问题。
他坐直了身体,双手交叉放在了桌上。
“想法不错,林同学。”王海峰开口了,语气虽然客气,但问题直指核心,“但是,你怎么处理这两个区域的交界面?显式算法和隱式算法的通量计算方式完全不同,如果在界面处直接拼接,必然会產生虚假的数值反射和振盪。这对於高超声速流动来说,任何一点虚假波动都会被激波放大,导致计算发散。”
“王教授说得对。”林叶早有准备,他自信地切换到下一页ppt,“硬拼接肯定不行。所以,我准备在两个区域之间,引入一个缓衝过渡层。”
“在这个过渡层里,不做那些非黑即白的切换,而是构造一个基於sigmoid函数的权重因子α,然后在过渡层內求解一个混合通量,大概长这样——”
【f—hybr id=af—implicit+(1—a)f—eplicit】
“同时,”林叶指著公式中的一项修正项,“为了保证守恆性,我利用了梅林变换技巧,从数学上推导了一个通量修正项。理论上,这能確保在过渡层內的质量和能量守恆误差控制在二阶精度以內。”
梅林变换!
看著林叶给出的这一步,下面的张涛也露出微微感慨的神情。
他和林叶前几天用到的方法,终究还是发挥了作用。
而王海峰看著那个精妙的权重函数构造,眼中的怀疑逐渐消退,取而代之的是一丝惊讶。
这个处理手法极其老练,完全不像是个新手。
“理论上————似乎是自洽的。”他沉吟著,脑海中开始推导后面的步骤,思考著具体的可行性。
“等等,我也有个问题。”一直沉默的李建国教授推了推眼镜,他是做物理模型的,更关心物理真实性,“你这个平衡区的判定,如果太宽泛,会不会漏掉一些虽然缓慢但很关键的中间反应?比如一氧化氮的生成反应比较慢,如果你把它当成平衡態处理了,热流预测可能会偏低。”
“这就是这个判据中渐近项的作用。”林叶解释道,“李教授,我引入的判据不是固定的閾值,而是与流场的对流时间尺度耦合的。如果流速很快,即使反应慢,判据也会收紧,强制將其划入刚性区计算。”
“从渐近分析的角度来看,只要判据的截断误差小於我们的离散误差,物理精度就是可以保证的。”
林叶並没有拿出具体的数据,因为他还没跑过,但他用严密的数学逻辑和物理图像,构建了一个看起来无懈可击的理论框架。
会议室里陷入了短暂的沉默。
王海峰教授盯著那页ppt看了足足半分钟,然后转过头,看向周文渊,语气中带著一丝认真:“老周,这路子有点野,但逻辑上是通的。比那个算不动的特徵值分解要有操作性得多。”
李建国教授也微微頷首:“只要那个渐近判据能如他所说的那样工作,物理上確实没问题。这实际上是用空间换时间,很聪明的做法。”
周文渊看著两位被说服的老搭档,心中一块大石落地,哈哈大笑了起来,他站起身,环视全场:“还有人有问题吗?”
三个团队的人都是你看看我,我看看你。
开玩笑,两位教授都没问题了,他们还能有什么问题呢?
“那,既然各位都觉得没问题了,那我们就试一试。海峰,这几天可能要辛苦你们团队了,把这个算法植入到求解器里去验证一下。”
“没问题!”王海峰大手一挥,“只要不再让我算那什么让人头大的特徵值就行了,写几行新代码算什么。今晚我就安排人手,咱们爭取三天內跑出个结果!”
>
本章未完,点击下一页继续阅读。