第100章 彻底解决刚性方程的办法 首席设计师
题。」
林叶一怔:「您请说。」
周文渊竖起一根手指:「第一,bch公式里涉及到了矩阵指数的计算,也就是e(at)。对于一个小矩阵来说,这很好算。但我们的流场有几百万个网格,每个网格有几十个组分,想要直接计算它的矩阵指数,其计算量是天文数字,比直接解线性方程组还要慢得多。」
「第二,就是李括号的计算代价。」周文渊竖起第二根手指:「你要修正误差,就要计算[f,s]。」
「在数学上,这等于计算ab和ba而已,没有区别。但是在物理上,这意味着你要计算流体雅可比矩阵和化学反应雅可比矩阵的乘积。对于复杂的化学反应机理,这个雅可比矩阵极其稠密且复杂,每算一次李括号,对于计算机来说都是一种极大的负担,而如果高阶项太多,这个修正过程本身就会成为新的计算瓶颈。」
说到这里,周文渊深深地看了林叶一眼,总结道:「所以,这个方法在理论上是完美的,但在工程实现上,往往因为计算代价过大而并不被接受,这也是它至今没有在高超声速cfd领域普及的原因。」
「你想做这个,就必须找到一种巧妙的方法,既能利用李代数的结构优势,又能绕过这些昂贵的矩阵运算。」
「这其中的难度,我想即使不用我多说,你应该也能够明白。」
周文渊说完,静静地看着林叶,似乎在等待他的反应。
这个年轻人是会知难而退?还是————
林叶听完这番分析,眼中的光芒非但没有黯淡,反而变得更加炽热。
矩阵指数难算?李括号代价高?
随着脑海中闪现过好几种之前学习过的方法和知识,他的心中顿时就隐隐有了一些新的想法。
他的脸上露出笑容。
「谢谢周老师专门帮我分析了,虽然您说的这两个问题可能确实存在,不过我还是挺想试一试的。」
周文渊的眼睛一瞪。
不是,你还真想试啊?
他算是看出来了,林叶现在是跟这个问题杠上了。
虽然从某种程度上来说,对于问题的坚持,是一种很不错的习惯。
但是他们现在讨论的这个问题,可是刚性方程问题,他们整个课题组都被这个问题给难住了。
早在之前的一场组会上,他们也讨论过能否搞出一种数学理论,将这个问题给解决。
但是,当时得出的就结论就是,很难,非常难。
他们三位子课题负责人,以及陈景明院士都表示过,如非必要,就不要将时间浪费在这个方法上面了,到时候未必能够搞出什么成果,还是老实走他们之前那种降维思路就行,稳扎稳打,至少搞出来的东西是看得到的。
再加上他们每位教授平时都有其他事情要忙,当然就更没空把太多时间放在这种事情上面了。
就比如周文渊,如果他手上只有两个项目的话,他可能还会尝试一下,但是他现在手上可是有六个项目,那就只能敬谢不敏了。
而林叶虽然不像是他们那样事务繁忙,但对于他也想要尝试这个方向,周文渊也是不支持的。
毕竟,像是林叶这样有天赋的学生,完全没必要在这
章节内容不完整,请退出阅读模式查看完整内容!