(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211249477.4 (22)申请日 2022.10.12 (71)申请人 国家超级计算天津中心 地址 300457 天津市滨 海新区经济技 术开 发区信环西路19号5号楼5102 申请人 先进计算与关键软件 （信创） 海河实 验室 (72)发明人 杨灿群　仲彦旭　王伟　黄颖杰　 夏梓峻　李健增　郑伟龙　卢海林　 (74)专利代理 机构 北京开阳星知识产权代理有 限公司 1 1710 专利代理师 吴崇 (51)Int.Cl. G06F 30/23(2020.01) G06F 30/20(2020.01) (54)发明名称 冲击波仿真方法、 装置、 设备及 介质 (57)摘要 本公开实施例涉及一种冲击波仿真方法、 装 置、 设备及介质， 其中该方法包括： 确定当前时间 步计算区域内的冲击波网格； 当计算区域的边界 区域中包括冲击波网格， 根据第一物理参数以及 第二物理参数， 构建计算区域的边界对应的边界 黎曼问题； 确定边界黎曼问题的求解结果， 根据 求解结果确定计算区域外的虚拟网格对应的第 一虚拟物理参数； 根据虚拟网格的第一虚拟物理 参数和计算区域内网格当前时间步的中间物理 参数， 确定预设流体控制方程的离散求解结果， 根据离散求解结果确定下一时间步的冲击波仿 真结果。 本公开实施例有效处理冲击波边界传播 问题， 抑制冲击波在计算区域边界出现的与真实 情况不符的非物理反射， 提高了仿真结果的准确 性。 权利要求书2页 说明书16页 附图9页 CN 115329646 A 2022.11.11 CN 115329646 A 1.一种冲击波仿真方法， 其特 征在于， 包括： 确定当前时间步计算区域内的冲击波网格； 当所述计算区域的边界区域中包括所述冲击波 网格， 根据第 一物理参数以及第 二物理 参数， 构建所述计算区域的边界对应的边界黎曼问题； 其中， 第一物理参数为所述计算区域 之外无穷远处网格的物理参数， 所述第二物理参数为所述计算区域内与其边界相 邻的网格 当前时间步的物理参数； 确定所述边界黎曼问题的求解结果， 根据 所述求解结果确定所述计算 区域外的虚拟网 格对应的第一虚拟物理参数； 根据所述虚拟网格的所述第一虚拟物理参数和所述计算区域内网格当前时间步的中 间物理参数， 确定预设流体控制方程的离散求解结果， 根据所述离散求解结果确定下一时 间步的冲击波仿真结果。 2.根据权利要求1所述的方法， 其特征在于， 所述确定当前时间步计算 区域内的冲击波 网格， 包括： 将所述计算区域中每 个网格依次确定为待处 理网格； 根据所述待处理网格以及所述待处理网格的前一网格当前时间步的中间物理参数， 确 定所述待处 理网格的逼近物理参数； 若所述待处理网格当前时间步的中间物理参数和所述逼近物理参数的差值绝对值大 于等于预设阈值， 则确定所述待处理网格为冲击波网格； 其中， 所述预设阈值为多分辨率参 数和网格尺寸 参数的乘积。 3.根据权利要求1所述的方法， 其特征在于， 所述根据第一物理参数以及第二物理参 数， 构建所述计算区域的边界对应的边界黎曼问题， 包括： 将所述第一物理参数作为所述 边界黎曼问题中第一侧对应的第一 黎曼参数值； 将所述第二物理参数作为所述 边界黎曼问题中第二侧对应的第二 黎曼参数值； 其中， 若所述第一侧为所述边界黎曼问题中的左侧， 则所述第二侧为所述边界黎曼问 题中的右侧； 若所述第一侧 为所述边界黎曼问题中的右侧， 则所述第二侧 为所述边界黎曼 问题中的左侧。 4.根据权利要求1所述的方法， 其特征在于， 所述确定所述边界黎曼问题的求解结果， 包括： 采用黎曼问题求解器对所述边界黎曼问题进行求解， 获取所述求解结果； 其中所述黎 曼问题求 解器包括： 兰金 ‑雨贡纽关系式， 以及接触间断速度与压力连续关系。 5.根据权利要求1所述的方法， 其特征在于， 所述根据所述求解结果确定所述计算区域 外的虚拟网格对应的第一虚拟物理参数， 包括： 提取所述求解结果中第 一物理参数对应的计算物 理参数， 将所述计算物理参数确定为 所述虚拟网格对应的所述第一虚拟物理参数。 6.根据权利要求1所述的方法， 其特征在于， 所述根据所述虚拟网格的所述第 一虚拟物 理参数和所述计算区域内网格当前时间步的中间物理参数， 确定预设流体控制方程的离散 求解结果， 根据所述离 散求解结果确定下一时间步的冲击波仿真结果， 包括： 将所述计算区域内相互接触的多物质进行解耦； 基于解耦后的多物质、 所述中间物理参数、 所述第 一虚拟物理参数， 对所述预设流体控权　利　要　求　书 1/2 页 2 CN 115329646 A 2制方程进行时间维度和空间维度的求 解， 获得所述预设流体控制方程对应的离 散方程； 基于所述离 散方程确定所述计算区域内网格下一时间步的中间物理参数。 7.根据权利要求1所述的方法， 其特 征在于， 所述方法还 包括： 当所述计算区域的边界区域中不包括所述冲击波网格， 将所述第 二物理参数确定为所 述虚拟网格对应的第二虚拟物理参数； 根据所述第 二虚拟物 理参数和所述计算区域内网格当前时间步的中间物理参数， 确定 预设流体控制方程的离散求解结果， 根据所述离散求解结果确定下一时间步的冲击波仿真 结果。 8.一种冲击波仿真装置， 其特 征在于， 包括： 第一确定模块， 用于确定当前时间步计算区域内的冲击波网格； 构建模块， 用于当所述计算区域的边界区域中包括所述冲击波网格， 根据第一物理参 数以及第二物理参数， 构建所述计算区域的边界对应的边界黎曼问题； 其中， 第一物理参数 为所述计算区域之外无穷远处网格的物理参数， 所述第二物理参数为所述计算区域内与其 边界相邻的网格当前时间步的物理参数； 求解模块， 用于确定所述边界黎曼问题的求解结果， 根据所述求解结果确定所述计算 区域外的虚拟网格对应的第一虚拟物理参数； 第二确定模块， 用于根据所述虚拟网格的所述第 一虚拟物理参数和所述计算 区域内网 格当前时间步的中间物理参数， 确定预设流体控制方程的离散求解结果， 根据所述离散求 解结果确定下一时间步的冲击波仿真结果。 9.一种电子设备， 其特 征在于， 所述电子设备包括： 处理器； 用于存储所述处 理器可执行指令的存 储器； 所述处理器， 用于从所述存储器中读取所述可执行指令， 并执行所述指令以实现上述 权利要求1 ‑7中任一所述的冲击波仿真方法。 10.一种计算机可读存储介质， 其特征在于， 所述存储介质存储有计算机程序， 所述计 算机程序用于执 行上述权利要求1 ‑7中任一所述的冲击波仿真方法。权　利　要　求　书 2/2 页 3 CN 115329646 A 3