间尺度鸿沟。
高超声速飞行时,流体对流的演化以微秒为单位,而高温下空气分子解离、电离的热化学反应往往在纳秒级完成。
二者跨度可达上千倍,将它们放在同一套方程组内求解,便造就了极强的方程刚性。
传统的数值解法面临两难。
若未匹配化学反应的极速变化,采用极小的时间步长,一次完整仿真的迭代次数会暴涨,即便动用国内顶尖超算,单次工况的演算时长也会拉长至数周,算力消耗达到无法承受的地步。
若是放大时间步长,适配流体流动,化学反应项就会出现严重的数值失稳,残差毫无征兆地猛烈反弹,仿真直接发散崩溃。
过去,行业内的妥协方案大多是简化化学反应组分,删减部分反应路径,牺牲物理精度换取计算稳定。
可一旦简化模型,高温再入超燃内部的热环境预测就会出现巨大偏差,仿真结果失去工程参考价值。
而这个问题更深层的难点是,没有现成的数学边界可供参考。
此前全球范围内,对于多组分热化学耦合刚性系统,不存在严谨解析化的刚性阈值划分。
国内外的项目只能依靠科研人员常年积累的经验,反复试凑调整迭代参数。
研究者能够看见仿真崩溃、残差震荡,却无法从理论上预判在何种参数、何种流场条件下,系统会走向失稳。
相当于医生只能看见病症,却找不到病灶的精确位置,调试全凭盲猜。
激波边界层干扰叠加热化学非平衡效应之后,局部区域的刚性强度会发生非线性突变,一处微小的网格扰动就足以诱发全域求解崩坏,问题的不可预测性进一步加剧。
另外还有一层隐藏的工程壁垒,便是全链路的误差耦合。
热化学模块并非独立运行,它需要和已经完成的置信体系无缝对接。每一处迭代产生的微小偏差,会顺着流体、热场、组分浓度多条路径互相放大。
想要完成误差溯源,就必须理清究竟是网格缺陷、离散格式误差,还是化学反应参数带来的漂移。
传统手段无法区分多重误差的叠加干扰,很容易出现修复一处缺陷,又引发另一处失稳的连锁问题。
所以说这个问题不是算力不够,而是缺少能够驾驭刚性方程的底层数学理论。
国内外很多数学家针对这个问题都进行过研究,直到目前还没有哪个数学家或者是科研机构公开表示解决这个问题。
本章未完,请点击下一页继续阅读!