一种有限催化分级隐式数值模拟方法、装置、设备及介质与流程

本发明涉及数值模拟计算领域，特别涉及一种有限催化分级隐式数值模拟方法、装置、设备及介质。

背景技术：

1、在高超声速飞行器再入飞行过程中，如果飞行速度较高（达到马赫10以上），会出现高温气体非平衡效应。飞行器头部激波后的气体温度可高达几千开尔文，高温气体会发生离解、电离、置换、复合等复杂化学反应，产生大量原子、离子气体组分，改变气体性质，对飞行器气动力、热等特性造成显著影响。流场中原子、离子等成分在到达飞行器表面时会发生催化复合反应（即表面催化效应），释放出很大的结合能，不仅影响非平衡流动，而且产生较大组分扩散热流，使飞行器表面气动热环境更加严酷。因此高温气体非平衡效应及表面催化效应的数值模拟技术，一直受到人们的广泛关注。

2、在数值模拟技术中，催化效应常以边界条件的形式出现。完全催化（fcw，(fullycatalytic wall）和完全非催化（ncw，non-catalytic wall）是最容易实现的两种表面催化边界条件，应用十分广泛。但这两者只模拟了表面催化的两种极限状态，要精确模拟催化效应，还须考虑表面材料具体的催化特性，采用有限催化条件（pcw，partially catalyticwall）边界计算模型。

3、高温气体表面有限催化边界的计算模拟相对复杂，需要在非平衡流动模拟过程中，耦合迭代求解表面复杂催化反应体系的质量平衡方程组以计算得到各气体组分的质量分数。在此过程中，由于流动参数与表面反应的相互干扰，常规的显式迭代方法稳定较差，很容易发散；常规的隐式处理方法，需要耦合求解各组分质量平衡方程，并处理与催化反应体系相关联的耦合系数矩阵，实现过程较为复杂，影响计算效率。

技术实现思路

1、有鉴于此，本发明的目的在于提供一种有限催化分级隐式数值模拟方法、装置、设备及介质，能够使得计算各气体组分质量分数的过程更加稳定高效。其具体方案如下：

2、第一方面，本技术公开了一种有限催化分级隐式数值模拟方法，包括：

3、非平衡流动数值模拟过程中，获取表面网格微元气体组分的质量分数的初始值作为当前值，并获取表面法向次层网格微元气体组分的质量分数作为辅助计算值；

4、根据表面催化复合反应体系，基于由低到高的顺序设定各气体组分的优先级；

5、根据所述当前值和表面网格微元材料特性参数计算中间参数；

6、基于所述中间参数、所述当前值和所述辅助计算值，并根据所述优先级分级隐式迭代计算所述表面网格微元气体组分的质量分数的预估值；

7、根据质量分数归一化原理和电中性原则并基于所述预估值，获取所述表面网格微元气体组分的质量分数的新值；

8、判断所述新值是否满足预设收敛条件；

9、若不满足所述预设收敛条件，则将所述新值作为所述当前值，并跳转至所述根据所述当前值和表面网格微元材料特性参数计算中间参数的步骤；

10、若满足所述预设收敛条件，则所述新值即为所述表面网格微元气体组分的质量分数的目标值。

11、可选的，所述根据表面催化复合反应体系，基于由低到高的顺序设定各气体组分的优先级，包括：

12、设置当前优先级n=1，设置所述表面催化复合反应体系为当前体系；

13、将所述当前体系中只在生成物中出现的气体组分的优先级预设为n级，并将所述生成物对应的反应物中气体组分的优先级预设为n+级；所述n+级表示高于n级的优先级；

14、在当前体系中，去除所述生成物的优先级预设为n级的催化复合反应，并将剩余部分重新设置为目标体系；

15、判断所述目标体系是否存在催化复合反应；

16、若存在，则重新设置所述当前优先级n=n+1，并将所述目标体系作为所述当前体系，然后跳转至所述将所述当前体系中只在生成物中出现的气体组分的优先级预设为n级的步骤；

17、若不存在，则对于表面网格微元上所有气体组分，确认优先级预设为a的气体组分的优先级为a，确认优先级预设为a+的气体组分的优先级为a+1，确认没有预设优先级的气体组分的优先级为无优先级。

18、可选的，所述根据所述当前值和表面网格微元材料特性参数计算中间参数，包括：

19、基于所述当前值、表面气体温度和压强计算气体密度；

20、基于所述表面催化复合反应体系中各催化复合反应的催化系数、重粒子反应物的分子量和所述表面气体温度，计算各催化复合反应分别对应的催化反应速率；

21、基于所述当前值计算各所述气体组分的扩散系数。

22、可选的，所述基于所述中间参数、所述当前值和所述辅助计算值，并根据所述优先级分级隐式迭代计算所述表面网格微元气体组分的质量分数的预估值，包括：

23、根据优先级由高到低的顺序选择当前气体组分；在表面催化复合反应体系中，选择所述当前气体组分参与的若干目标催化复合反应；

24、基于所述若干目标催化复合反应的所述催化反应速率，计算当前气体组分催化反应质量生成率；

25、基于所述当前气体组分的所述扩散系数、所述质量生成率和所述辅助计算值，隐式离散所述当前气体组分的表面质量守恒方程，计算得到所述当前气体组分的质量分数的预估值；

26、跳转至所述根据优先级由高到低的顺序选择当前气体组分的步骤，直至得到所有具有优先级的气体组分质量分数的所述预估值。

27、可选的，所述有限催化分级隐式数值模拟方法，还包括：

28、所述表面质量守恒方程中，表面网格微元上所有优先级高于所述当前气体组分的优先级的气体组分的质量分数采用所述预估值；

29、所述表面质量守恒方程中，表面网格微元上所有优先级不高于所述当前气体组分的优先级的气体组分和无优先级的气体组分的质量分数采用所述当前值。

30、可选的，所述根据质量分数归一化原理和电中性原则并基于所述预估值，获取所述表面网格微元气体组分的质量分数的新值，包括：

31、由所述辅助计算值计算无优先级的气体组分的质量分数的预估值；

32、若飞行器表面存在电子组分，则基于电中性原则计算所述电子组分的质量分数的预估值；

33、基于具有优先级的气体组分质量分数的预估值、所述无优先级的气体组分的质量分数的预估值以及可能存在的电子组分质量分数的预估值的质量分数归一化，获取所述表面网格微元气体组分的质量分数的新值。可选的，所述预设收敛条件为最大相对残差不大于预设精度标准；所述最大相对残差为所有所述气体组分根据预设公式计算的预设相对残差中的最大值；所述预设公式为：

34、；

35、；

36、其中，所述为不同所述气体组分的所述预设相对残差，所述表示不同所述气体组分；所述为所述气体组分的所述新值，所述为所述气体组分的所述当前值；所述为所述气体组分的总个数。

37、第二方面，本技术公开了一种有限催化分级隐式数值模拟装置，包括：

38、数值获取模块，用于非平衡流动数值模拟过程中，获取表面网格微元气体组分的质量分数的初始值作为当前值，并获取表面法向次层网格微元气体组分的质量分数作为辅助计算值；

39、优先级设定模块，用于根据表面催化复合反应体系，基于由低到高的顺序设定各气体组分的优先级；

40、参数计算模块，用于根据所述当前值和表面网格微元材料特性参数计算中间参数；

41、预估值计算模块，用于基于所述中间参数、所述当前值和所述辅助计算值，并根据所述优先级分级隐式迭代计算所述表面网格微元气体组分的质量分数的预估值；

42、新值计算模块，用于根据质量分数归一化原理和电中性原则并基于所述预估值，获取所述表面网格微元气体组分的质量分数的新值；

43、判断模块，用于判断所述新值是否满足预设收敛条件；

44、条件不满足模块，用于若不满足所述预设收敛条件，则将所述新值作为所述当前值，并跳转至所述根据所述当前值和表面网格微元材料特性参数计算中间参数的步骤；

45、条件满足模块，用于若满足所述预设收敛条件，所述新值即为表面网格微元气体组分质量分数的目标值。

46、第三方面，本技术公开了一种电子设备，包括处理器和存储器；其中，所述处理器执行所述存储器中保存的计算机程序时实现前述公开的有限催化分级隐式数值模拟方法。

47、第四方面，本技术公开了一种计算机可读存储介质，用于存储计算机程序；其中，所述计算机程序被处理器执行时实现前述公开的有限催化分级隐式数值模拟方法。

48、可见，本技术非平衡流动数值模拟过程中，获取表面网格微元气体组分的质量分数的初始值作为当前值，并获取表面法向次层网格微元气体组分的质量分数作为辅助计算值；根据表面催化复合反应体系，基于由低到高的顺序设定各气体组分的优先级；根据所述当前值和表面网格微元材料特性参数计算中间参数；基于所述中间参数、所述当前值和所述辅助计算值，并根据所述优先级分级隐式迭代计算所述表面网格微元气体组分的质量分数的预估值；根据质量分数归一化原理和电中性原则并基于所述预估值，获取所述表面网格微元气体组分的质量分数的新值；判断所述新值是否满足预设收敛条件；若不满足所述预设收敛条件，则将所述新值作为所述当前值，并跳转至所述根据所述当前值和表面网格微元材料特性参数计算中间参数的步骤；若满足所述预设收敛条件，则所述新值即为所述表面网格微元气体组分的质量分数的目标值。由此可见，本技术利用分级的隐式迭代方法计算质量分数，避免了使用显示迭代方法，也避免了全耦合隐式方法复杂的求解过程，使得计算过程更加稳定高效。