[发明专利]直接模拟法与预估修正准静态方法的衔接计算方法

说明书

本公开公开了直接模拟法与预估修正准静态方法的衔接计算方法及装置、电子设备和存储介质，涉及数据处理领域，主要技术方案包括：基于直接模拟法计算至预设衔接时刻；在预设衔接时刻，获取求解形状函数所需的第一相关变量；在预设衔接时刻，获取求解点堆方程组需要的第二相关变量；将第一相关变量与第二相关变量输入预设格式文件；获取预设格式文件的中间状态文件，并进行预估修正准静态法计算。与相关技术相比，在预设衔接时刻之前使用直接模拟法进行计算，在预设衔接时刻，使用预估修正准静态方法对第一与第二相关变量进行计算，实现了直接模拟法与预估修正准静态方法的衔接计算，进而实现了对中子场三维几何问题的瞬态计算与高效模拟。

技术领域

本公开涉及数据处理技术领域，尤其涉及一种直接模拟法与预估修正准静态方法的衔接计算方法及装置、电子设备和存储介质。

背景技术

在无源启动计算过程中，初期处于弱中子场时，首先根据随机中子动力学理论计算，随机中子动力学基于直接模拟法求解，理论上可以一直随时间进行输运计算，但直接模拟法的求解较为耗时。

当中子数增多直到进入可用平均中子数描述的强中子场时，可以用计算速度更快的预估修正准静态方法(Predictor-Corrector Quasi-Static，PCQS)描述系统，目前国际上仅研究了点堆的衔接计算，其无法针对中子场三维复杂几何问题进行模拟计算。

在相关技术中还没有一种方法可以针对中子场三维复杂几何问题进行模拟计算。因此，推演一种可以针对中子场三维复杂几何问题进行模拟计算的方法是目前急需解决的问题。

发明内容

本公开提供了一种直接模拟法与预估修正准静态方法的衔接计算方法、装置、电子设备和存储介质。其主要目的在于解决在相关技术中还没有一种方法可以针对中子场三维复杂几何问题进行模拟计算的问题。其主要目的在于推演一种可以针对中子场三维复杂几何问题进行模拟计算的方法且实现对中子场三维几何问题的瞬态计算与高效模拟。

根据本公开的第一方面，提供了一种直接模拟法与预估修正准静态方法的衔接计算方法，其中，包括：

基于直接模拟法计算至预设衔接时刻；

在所述预设衔接时刻，获取求解形状函数所需的第一相关变量，所述第一相关变量包括：缓发中子源与时间吸收源；

在所述预设衔接时刻，获取求解点堆方程组需要的第二相关变量，所述相关变量包括：反应性，有效缓发中子份额与中子代时间；

基于计算机语言脚本，将所述第一相关变量与所述第二相关变量输入预设格式文件，所述预设格式文件与预设预估修正准静态法框架相对应；

基于预设预估修正准静态计算程序，获取所述预设格式文件的中间状态文件，并进行预估修正准静态法计算，所述中间状态文件为每经过一预设时间段，计算得到的所述第一相关变量与所述第二相关变量组成的预设格式文件。

可选的，所述在所述预设衔接时刻，获取求解形状函数所需的第一相关变量之前，所述方法还包括：

构建时空中子动力学方程组，推导得固定源形式的方程，所述时空中子动力学方程组包括中子通量方程与缓发中子先驱核方程，所述固定源为除裂变源外的源项，且由所述时间吸收源和所述缓发中子源组成；

基于所述固定源形式的方程，推导得到形状函数。

可选的，所述在所述预设衔接时刻，获取求解形状函数所需的第一相关变量，包括：

对时间吸收源分布进行计算，得到所述时间吸收源，所述时间吸收源分布为使用直接模拟法计算至预设衔接时刻时，中子库里的中子分布；

基于缓发中子源的公式，对预设衔接时刻的缓发中子源分布进行计算，得到所述缓发中子源。