公开(公告)号：CN115329689A

公开(公告)日：2022-11-11

申请号：CN202210846089.8

申请日：2022-07-05

Applicant: 北京航空航天大学

Inventor： 韩明泽 ,  姜振华 ,  闫超

IPC: G06F30/28 ,  G06F30/22 ,  G06F30/15 ,  G06F17/11 ,  G06F111/10 ,  G06F113/08 ,  G06F119/14

Abstract: 本发明公开一种基于伪非定常时间推进的复杂湍流流动高效率计算方法，包括：步骤1，建立非定常湍流模拟的控制方程；步骤2，确定求解区域，在区域内部划分网格单元，对所有的网格单元赋流场初值；步骤3，利用当前tn时刻网格单元变量，计算所有网格单元的无粘通量和粘性通量，并求和，作为当前tn时刻的网格单元流场变量时间导数值；步骤4，利用当前tn时刻的流场变量的时间导数值求解下一步tn+1时刻的流场变量值；步骤5，判断下一步tn+1时刻是否大于或等于计算终止tend时刻，如果满足，则转到步骤6；否则令n＝n+1，返回步骤2；步骤6，输出最终结果。本发明方法可有效保证矩阵对角占优性质，保证数值计算稳定性，从而提高计算效率，减小计算周期和计算成本。

公开(公告)号：CN115169259A

公开(公告)日：2022-10-11

申请号：CN202210782333.9

申请日：2022-07-05

Applicant: 北京航空航天大学

Inventor： 姜振华 ,  韩明泽 ,  闫超

IPC: G06F30/28 ,  G06F30/22 ,  G06F30/15 ,  G06F17/13 ,  G06F17/11 ,  G06F17/10 ,  G06F111/10 ,  G06F113/08 ,  G06F119/14

Abstract: 本发明公开一种高分辨率激波稳定的高超声速流动预测方法，主要包括：流动变量重构步骤：在MUSCL插值框架下，引入双曲正切函数近似，通过相邻单元的插值将数值近似解在单元里线性分布，得到网格单元界面两侧守恒变量；控制方程求解步骤：用重构步中得到的网格单元界面两侧守恒变量，计算界面两侧对流通量，利用网格单元中心点的守恒变量计算网格单元界面的特征速度，该方法能够有效减弱激波异常现象，进一步增强高超声速流动计算稳定性，并且不会对预测精度产生较大影响。本发明方法可以更好的实现激波稳定计算；可以有效提高对间断的分辨率；能够改善高超声速流动计算稳定性，尤其是能更好的实现各类激波干扰复杂流动的激波稳定计算。