(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 20221081713 0.9 (22)申请日 2022.07.12 (71)申请人 清华大学 地址 100084 北京市海淀区清华园 (72)发明人 王侃　贾从龙　李泽光　王武　 (74)专利代理 机构 北京清亦华知识产权代理事 务所(普通 合伙) 11201 专利代理师 孟洋 (51)Int.Cl. G06F 30/20(2020.01) G06F 119/14(2020.01) (54)发明名称 直接模拟法与预估修正准静态方法的衔接 计算方法 (57)摘要 本公开公开了直接模拟法与预估修正准静 态方法的衔接计算方法及装置、 电子设备和存储 介质， 涉及数据处理领域， 主要 技术方案包括： 基 于直接模拟法计算至预设衔接时刻； 在预设衔接 时刻， 获取求解形状函数所需的第一相关变量； 在预设衔接时刻， 获取求解点堆方程组需要的第 二相关变量； 将第一相关变量与第二相关变量输 入预设格式文件； 获取预设格式文件的中间状态 文件， 并进行预估修正准静态法计算。 与相关技 术相比， 在预设衔接时刻之前使用直接模拟法进 行计算， 在预设衔接时刻， 使用预估修正准静态 方法对第一与第二相关变量进行计算， 实现了直 接模拟法与预估修正准静态方法的衔接计算， 进 而实现了对中子场三维几何问题的瞬态计算与 高效模拟。 权利要求书2页 说明书11页 附图4页 CN 115374598 A 2022.11.22 CN 115374598 A 1.一种直接模拟法与预估修 正准静态方法的衔接计算方法， 其特 征在于， 包括： 基于直接模拟法计算至预设衔接时刻； 在所述预设衔接时刻， 获取求解形状函数所需的第一相关变量， 所述第一相关变量包 括： 缓发中子源与时间吸 收源； 在所述预设衔接时刻， 获取求解点堆方程组需要的第 二相关变量， 所述相关变量包括： 反应性， 有效缓发中子份额与中子代时间； 基于计算机语言脚本， 将所述第一相关变量与所述第二相关变量输入预设格式文件， 所述预设格式文件与预设预估修 正准静态法框架相对应； 基于预设预估修正准静态计算程序， 获取所述预设格式文件的中间状态文件， 并进行 预估修正准静态法计算， 所述中间状态文件为每经过一预设时间段， 计算得到的所述第一 相关变量与所述第二相关变量组成的预设格式文件。 2.根据权利要求1所述的直接模拟法与预估修正准静态方法的衔接计算方法， 其特征 在于， 所述在所述预设衔接时刻， 获取求解形状函数所需的第一相关变量之前， 所述方法还 包括： 构建时空中子动力学方程组， 推导得固定源形式的方程， 所述时空中子动力学方程组 包括中子通量方程与缓发中子先驱核方程， 所述固定源为除裂变源外的源项， 且由所述时 间吸收源和所述缓发中子源组成； 基于所述固定源形式的方程， 推导得到形状函数。 3.根据权利要求1所述的直接模拟法与预估修正准静态方法的衔接计算方法， 其特征 在于， 所述在所述预设衔接时刻， 获取求解形状函数所需的第一相关变量， 包括： 对时间吸 收源分布进行计算， 得到所述时间吸收源， 所述时间吸收源分布为使用直接模拟法计算至 预设衔接时刻时， 中子库里的中子分布； 基于缓发中子源的公式， 对预设衔接时刻的缓发中子源分布进行计算， 得到所述缓发 中子源。 4.根据权利要求1所述的直接模拟法与预估修正准静态方法的衔接计算方法， 其特征 在于， 所述在所述预设衔接时刻， 获取求解点堆方程组需要的第二相关变量之前， 所述方法 还包括： 基于时空中子动力学方程组， 引入临界伴随通量加权函数， 推导出点堆方程组， 所述点 堆方程组为幅度函数随时间变化的描述。 5.根据权利要求1所述的直接模拟法与预估修正准静态方法的衔接计算方法， 其特征 在于， 所述在所述预设衔接时刻， 获取求 解点堆方程组需要的第二相关变量， 包括： 基于临界计算， 得到有效增殖因子keff， 使用所述有效增殖因子keff计算， 得到反应 性； 基于反复裂 变几率法计算， 得到有效缓发中子份额与中子代时间。 6.一种直接模拟法与预估修 正准静态方法的衔接计算装置， 其特 征在于， 包括： 第一计算单 元， 用于基于直接模拟法计算至预设衔接时刻； 第一获取单元， 用于在所述预设衔接时刻， 获取求解形状函数所需的第 一相关变量， 所 述第一相关变量包括： 缓发中子源与时间吸 收源； 第二获取单元， 用于在所述预设衔接时刻， 获取求解点堆方程组需要的第 二相关变量，权　利　要　求　书 1/2 页 2 CN 115374598 A 2所述相关变量包括： 反应性， 有效缓发中子份额与中子代时间； 输入单元， 用于基于计算机语言脚本， 将所述第一相关变量与所述第二相关变量输入 预设格式文件， 所述预设格式文件与预设预估修 正准静态法框架相对应； 第二计算单元， 用于基于预设预估修正准静态计算程序， 获取所述预设格式文件的中 间状态文件， 并进 行预估修正准静态法计算， 所述中间状态文件为每经过一预设时间段， 计 算得到的所述第一相关变量与所述第二相关变量组成的预设格式文件。 7.根据权利要求6所述的直接模拟法与预估修正准静态方法的衔接计算装置， 其特征 在于， 所述装置还 包括： 第一推导单元， 用于构建时空中子动力学方程组， 推导得固定源形式的方程， 所述时空 中子动力学方程组包括中子通量方程与缓发中子先驱核方程， 所述固定源为除裂变源外的 源项， 且由所述时间吸 收源和所述缓发中子源组成； 第二推导单 元， 用于基于所述固定源形式的方程， 推导得到形状函数。 8.根据权利要求6所述的直接模拟法与预估修正准静态方法的衔接计算装置， 其特征 在于， 所述第一获取 单元包括： 第一计算模块， 用于对时间吸收源分布进行计算， 得到所述时间吸收源， 所述时间吸收 源分布为使用直接模拟法计算至预设衔接时刻时， 中子库里的中子分布； 第二计算模块， 用于基于缓发中子源的公式， 对预设衔接时刻的缓发中子源分布进行 计算， 得到所述缓发中子源。 9.根据权利要求6所述的直接模拟法与预估修正准静态方法的衔接计算装置， 其特征 在于， 所述装置还 包括： 第三推导单元， 用于基于时空中子动力学方程组， 引入临界伴随通量加权函数， 推导出 点堆方程组， 所述 点堆方程组为幅度函数随时间变化的描述。 10.根据权利要求6所述的直接模拟法与 预估修正准静态方法的衔接计算装置， 其特征 在于， 所述第二获取 单元包括： 第一计算模块， 用于基于临界计算， 得到有效增殖因子keff， 使用 所述有效增殖因子 keff计算， 得到反应性； 第二计算模块， 用于基于反复裂 变几率法计算， 得到有效缓发中子份额与中子代时间。 11.一种电子设备， 其特 征在于， 包括： 至少一个处 理器； 以及 与所述至少一个处 理器通信连接的存 储器； 其中， 所述存储器存储有可被所述至少一个处理器执行的指令， 所述指令被所述至少一个处 理器执行， 以使所述至少一个处 理器能够执 行权利要求1 ‑5中任一项所述的方法。 12.一种存储有计算机指令的非瞬时计算机可读存储介质， 其特征在于， 所述计算机指 令用于使所述计算机执 行根据权利要求1 ‑5中任一项所述的方法。 13.一种计算机程序产品， 其特征在于， 包括计算机程序， 所述计算机程序在被处理器 执行时实现根据权利要求1 ‑5中任一项所述的方法。权　利　要　求　书 2/2 页 3 CN 115374598 A 3