第328章 这是什么情况（2/2）

韩柏岩开口的时候嗓子有些乾涩，哪怕是肖宿，也不能隨隨便便就否认他们研究院辛苦的成果，总得说个一二三吧。

肖宿听出他的认真，看著他的眼睛，郑重的解释道：

“这个计算模型当然能求出解。

或者换句话说，任何一个不完备的离散化格式，只要网格划得够密、叠代步数给得够多，总能给出一个数值上的结果。

但是能求出解和能用是两回事。”

他顿了顿，用手指点了点被圈起来的那部分。

“你们想用第二步叠代来修正离散误差，这个想法本身是有道理的。

但问题是，修正项本身在连续极限下是不收敛的，截断误差里藏了一个非零的拓扑阻碍项。

这个阻碍项在网格加密的时候不会自动消失，反而会在边界条件变化的时候被放大。

你们搭的这座桥，不仅没有弥合两个框架的裂缝，反而正好把力量传导到了这个隱藏的裂缝上，把它给撕开了。

换句话说，你们这个方法，算出正確答案的概率比算出错误答案的概率还小。

你们说结果跟实验数据对得上，只是因为碰巧阻碍项在你们选的那组参数下没有触发而已，但碰巧对了一次不代表这个方法是对的。

换一组边界条件，换一个来流马赫数，你再去跑，结果还能对上吗？”

韩柏岩没说话，但他脸上的表情已经给出了答案。

他们確实跑过好几组不同的参数，有些组收敛得很好，有些组直接发散，还有几组跑出来的结果明显是错的，他们当时以为是初值给得不好，但是反覆调了几次初值也没找到规律。

原来，不是初值不好，而是方法本身就不可靠。

那边，肖宿又仔细翻了翻韩柏岩他们提供的数据，在脑海中计算了一下，说道：

“在实验的时候，你们的飞行高度大概在四十到六十公里之间，来流马赫数应该在十五到二十之间，对吧。

取的原始壁面温度边界条件是辐射平衡假设下的等温壁，壁面发射率取的零点八五左右。”

这话一出，韩柏岩浑身一僵，后背一阵发凉。

这些数据他们根本没写在材料里。

他们给肖宿的材料里只包含了模型结构和部分经过处理的输出结果，原始输入参数和具体工况数据出於安全考虑全部摘掉了。

可肖宿竟然仅凭这套模型，就精准还原了他们所有未公开的实验细节，说出来的每一个数字还都跟实际情况分毫不差！

这是什么情况。

肖宿可不管他们心里怎么想的，继续推道：

“然后，你们把壁面沿轴向切成了若干个控制体，每个控制体上对嵌入曲率做局部线性化，再用有限体积法在每一个控制体上分別求解热流密度。

接著，你们把相邻控制体之间的热流密度差值作为一个额外的源项，重新代回到嵌入曲率的连续表示里做了第二步叠代。

你们想用这一步在离散解和连续表示之间搭一座桥，弥补两种框架的兼容性问题，但是恰好就在这一步把拓扑阻碍项激活了。”