第327章 不是很难(2/2)
他们把法旗方向曲率分解中的嵌入曲率部分单独提出来做了积分。
肖宿越看,眉头皱的越紧。
他们求最优解,竟然是在实空间里沿著激波面做线积分,而不是在模式空间里沿著嵌入曲率的最速下降方向做路径积分。
这就好比你要爬一座山,明明有索道可以直接上去,你偏要从山脚下绕盘山路走,花的功夫是別人的几十倍,最后得到的结果还未必对。
连积分路径都能选错。
肖宿抿著唇,一言不发地翻到了第三部分。
航天科工自己做的这部分问题也很大。
他们竟然把嵌入曲率分解和传统的有限体积法硬拼在了一起。
要知道,嵌入曲率分解是基於辛几何的连续谱表示,有限体积法是基於离散网格的代数求解,这两个框架的数学基础根本不兼容的。
他们显然也认识到了这一点,在后面又做了一系列修正,试图用相邻控制体之间的热流密度差值作为额外源项,重新代回连续框架里做第二步叠代,想靠这个闭环把两个框架的误差给磨平。
在工程上,这是一种相当有魄力的尝试。
可问题是,这两套框架的数学根基从根本上就是拧著的,你想在它们之间搭桥,就像想把两根不同材质的绳子编成一根,每一股力都较著劲,最后谁先断,就看这股劲憋在哪个看不见的地方了。
肖宿的表情越看越凝重。
他嘴唇紧紧抿著,翻页的速度越来越快。
看到最后一页的时候,他手指压在某一行公式上停了好一会儿。
这个热流密度梯度的离散化格式,公式本身推得倒是真长,从第一行到最后一行写了將近一页纸。
但肖宿一眼就看出来,这个离散化格式在连续极限下根本不可能收敛到原来的嵌入曲率表示,它的截断误差项里藏著一个非零的拓扑阻碍。
这个阻碍项在网格加密的时候不会消失,反而会在某些边界条件下被放大。
换句话说,他们搭的那座桥,正好一头杵在了这个隱藏的拓扑阻碍上。
桥没建成,反倒把它给激活了。
林瑾和韩柏岩一直盯著肖宿的表情。
看到他眉头越皱越紧、翻页越来越快,韩柏岩的心臟就跟著越跳越快。
肖宿脸上很少有明显情绪,之前二把手开会的时候他都是波澜不惊的样子,现在眉头紧锁,嘴唇紧抿,显然不是看到了什么令他满意的东西。
“肖教授,”林瑾忍不住问了一句,“这些问题很难吗?”
肖宿抬眼看了他们一眼,抿了抿唇。
“不是很难。”