公开(公告)号：CN111584104A

公开(公告)日：2020-08-25

申请号：CN202010382151.3

申请日：2020-05-08

Applicant: 西安交通大学

Inventor： 贺清明 ,  黄展鹏 ,  郑琪 ,  曹良志 ,  吴宏春

IPC: G21C17/00 ,  G21C17/108

Abstract: 一种基于多次碰撞源-蒙特卡罗耦合获取反应堆外探测器响应的方法，将反应堆区域分为堆内A区与堆外B区，根据多次碰撞源技术，碰撞源产生的中子通量密度为碰撞中子通量密度；将B区中堆外探测器处中子通量密度按中子来源分为两部分，第一部分为A区内第0→n次碰撞中子通量密度总和第二部分为A区内大于n次碰撞中子通量密度总和与B区内碰撞中子通量密度之和φM(r,E)；用确定论方法计算 和第n+1次碰撞源 用蒙特卡罗方法从 出发计算φM(r,E)；将两部分中子通量密度相加得到目标探测器处的中子通量密度；本发明利用多次碰撞源技术与蒙特卡罗方法在A、B区域中计算的优点，提升堆外探测器响应计算的统计效率与精度，为堆外探测器在实际反应堆探测时的响应提供了参考，提高反应堆安全性。

公开(公告)号：CN107122545B

公开(公告)日：2020-08-04

申请号：CN201710286988.6

申请日：2017-04-27

Applicant: 西安交通大学

Inventor： 刘宙宇 ,  王博 ,  曹良志 ,  吴宏春

IPC: G06F30/10 ,  G06F30/20 ,  G06F119/14 ,  G06F111/10

Abstract: 一种精确计算核反应堆时空中子分布的方法，包括如下步骤：1、使用一步法直接求解稳态中子输运方程获得稳态状态下的中子通量密度分布以及先驱核浓度；2、在0.25ms内求解时空中子输运方程，通过0.25时刻的中子通量密度与中子通量密度形状函数获得0.25ms的幅值，此幅值作为点堆计算的初始值；3、在大步长内使用一步法计算时空中子输运方程，获得预估解，在极小步长内求解点堆方程，用点堆的幅值去校正中子通量密度；本发明通过改进现有的预估校正准静态方法，在不显著增加计算时间的前提下，能够精确捕捉核反应堆在极短时间内中子的迅速变化，彻底克服传统预估校正准静态方法的缺陷，使其在高保真时空中子动力学计算中真正发挥作用。

公开(公告)号：CN107092781B

公开(公告)日：2020-03-24

申请号：CN201710213114.8

申请日：2017-04-05

Applicant: 西安交通大学

Inventor： 祖铁军 ,  李继奎 ,  曹良志

IPC: G06F17/10

Abstract: 一种获得核反应堆高保真共振群常数的计算方法，基于ESSM方法的高效率迭代计算，获得目标共振核素在目标能群背景截面，根据非均匀共振积分表拟合得到的背景截面和栅元尺寸关系公式，得到等效栅元尺寸，保持燃料、包壳等几何、材料、温度不变，建立等效2‑D栅元；求解相关超细群中子慢化方程，从而获得高保真共振群常数；本发明相比非确定论方法能显著提高计算效率，提高经济性；同时相比基于直接求解问题相关慢化方程，该方法可以显著减小碰撞概率计算时间，提高计算效率，同时规避存储量巨大的缺陷；而相比基于共振积分表的共振计算方法，2‑D栅元中子慢化方程的求解可以精确处理共振干涉、空间自屏和温度分布效应，获得每个燃料内精细分布的高保真共振群常数。

公开(公告)号：CN110826205A

公开(公告)日：2020-02-21

申请号：CN201911028486.9

申请日：2019-10-28

Applicant: 西安交通大学

Inventor： 李云召 ,  王思成 ,  曹良志 ,  吴宏春

IPC: G06F30/20

Abstract: 本发明公开了一种适用于压水堆Pin-by-pin计算的能群结构优化方法，1、选取特征组件集合，2、从特征组件集合中首先选定一个目标组件，进行全反射边界条件下的单组件非均匀计算，3、将其余组件与选定的组件组合，进行考虑该目标组件环境效应的多组件计算，4、求得的各栅元归一化能谱求比值，5、依据规则进行少群能群结构能群边界的划分，6、重复步骤2到步骤5的计算，得到一系列少群能群结构划分方案，7、评价各个能群结构划分边界，得到一套优化后的能群结构优化方案；通过本发明优化方法得到一套能群数少于七群且能够实际应用于压水堆Pin-by-pin计算的少群能群结构；在保证计算精度的前提下，达到减少计算存储、提高计算效率的目标。

公开(公告)号：CN110579790A

公开(公告)日：2019-12-17

申请号：CN201910982838.8

申请日：2019-10-16

Applicant: 西安交通大学

Inventor： 张清民 ,  桑耀东 ,  邓邦杰 ,  曹良志 ,  李云召

IPC: G01T7/00

Abstract: 一种自给能中子探测器静电效应的评估方法，步骤如下：(1)由蒙特卡罗程序模拟得到零时刻绝缘层电荷沉积速率Vq,dep,0(r,t)，并且此时电荷流失速率Vq,out,0(r,t)为零；(3)由Vq,dep,0(r,t)和Vq,out,0(r,t)得到零时刻电荷净沉积速率Vq,0(r,t)；(4)设定时间步长为Δt，将前一个时间点的电荷沉积速率、电荷流失速率作为当前时间点的初值，得到当前时间绝缘层电荷分布Q(r,t)；(5)根据Q(r,t)以及绝缘层体积获得绝缘层电荷密度ρ0(r,t)，并拟合得到其解析式ρ(r,t)；(6)将ρ(r,t)代入泊松方程并求解，得到绝缘层电势与电场；(7)将电场加入蒙特卡罗模拟过程中得到当前时间点新的电荷净沉积速率，并重复步骤(4)-(7)，直到得到稳定的电荷净沉积速率Vq(r,t)；(8)若Vq(r,t)满足收敛要求则输出电场与电势，否则重复上述步骤。

公开(公告)号：CN107103189B

公开(公告)日：2019-08-09

申请号：CN201710240986.3

申请日：2017-04-13

Applicant: 西安交通大学

Inventor： 李云召 ,  张斌 ,  吴宏春 ,  曹良志

IPC: G06F17/10

Abstract: 一种进行反应堆临界曲率搜索的方法，包括如下步骤：1、统计栅元信息，对非均匀组件进行高阶输运计算，得到组件中各个栅元的均匀化截面以及中子通量密度分布；2、求解泄漏修正模型进行堆芯临界曲率搜索计算，根据堆芯临界曲率，求解堆芯特征值，判断该特征值是否为真实堆芯特征值，若是则临界曲率搜索结束，若不是则更新堆芯临界曲率并再次进行泄露修正模型的计算，直至算出的特征值为真实堆芯特征值并得到渐近能谱；3、对栅格计算得到的中子通量密度乘上渐近能谱和栅格计算中无限中子能谱的比值并结束泄漏修正的临界曲率搜索；本发明方法基于中子平衡方程进行推导，能够快速找到对应临界堆芯下的临界曲率。

公开(公告)号：CN108733903B

公开(公告)日：2019-06-11

申请号：CN201810429333.4

申请日：2018-05-08

Applicant: 西安交通大学

Inventor： 李云召 ,  郑琪 ,  吴宏春 ,  曹良志

IPC: G06F17/50

Abstract: 一种按中子能量进行个性化处理的中子输运数值模拟方法，通过建立相同的平源区结构和能群结构，在快中子能量段和热中子能量段使用特征线方法扫描计算，在共振能段使用蒙特卡罗模拟，不同能量分段之间通过散射源和裂变源项互相耦合，根据收敛情况再进行裂变源迭代，以得到中子输运数值模拟需要的特征值和通量分布；本发明能利用特征线方法处理非共振能量段，保持在这些能量段中较快的计算速度；能采用蒙特卡罗方法直接抽样模拟特征线方法计算不够准确的共振能量段；在缩短计算时间的同时，还能精确处理共振能量段，提高计算精度。

公开(公告)号：CN109493924A

公开(公告)日：2019-03-19

申请号：CN201811347675.8

申请日：2018-11-13

Applicant: 西安交通大学

Inventor： 祖铁军 ,  尹文 ,  曹良志 ,  吴宏春 ,  贺清明 ,  刘宙宇

IPC: G16C20/40

Abstract: 一种获取FCM燃料有效多群截面的方法，在共振能量段，基于超细群方法，求解包含TRISO颗粒和基质材料的一维球模型得到超细群缺陷因子，通过超细群缺陷因子修正所有核素的超细群截面即可将颗粒和基质均匀化，利用丹可夫因子等价得到FCM燃料的各个燃料棒的等价一维棒模型，对一维棒模型求解超细群慢化方程从而得到共振能量段的有效自屏截面，在热能段，利用穿透概率和碰撞概率等价得到多群缺陷因子，通过多群缺陷因子修正所有核素的多群截面将燃料和基质均匀化并得到热能段有效多群截面；本发明能有效地处理FCM燃料的双重非均匀性效应，得到精确的有效自屏截面。

公开(公告)号：CN106126928B

公开(公告)日：2018-07-03

申请号：CN201610471991.0

申请日：2016-06-24

Applicant: 西安交通大学

Inventor： 曹良志 ,  李志峰 ,  吴宏春 ,  祖铁军

IPC: G06F19/00

Abstract: 获得固态及液态氟盐热中子散射数据库的方法及数据库，该方法步骤为：1、对于固态氟盐晶体及液态氟盐溶液，使用基于第一性量子力学原理的晶格动力学程序计算出其声子谱分布；此外，对于液态氟盐溶液，使用基于分子动力学的统计物理方法统计出各氟盐离子的扩散系数；2、将步骤1所获得的声子谱分布、扩散系数提供给核数据加工程序NJOY，获得固态及液体氟盐的热中子散射数据；3、根据不同的应用目的，配套相应的核数据库索引文件，索引文件与步骤2获得的热中子散射数据一起构成完整的热中子散射数据库；本发明还公开了提供该方法获得的数据库；填补了国际上固态氟盐晶体及液体氟盐溶液多群输运及子群共振计算用热中子散射数据库制作领域的空白。

公开(公告)号：CN107423542A

公开(公告)日：2017-12-01

申请号：CN201710240987.8

申请日：2017-04-13

Applicant: 西安交通大学

Inventor： 李云召 ,  张斌 ,  吴宏春 ,  曹良志

IPC: G06F19/00

Abstract: 一种适用于逐棒计算的非均匀泄漏修正方法，包括如下步骤：1、统计栅元信息，首先对非均匀组件进行高阶输运计算，然后计算得到组件中各个栅元的均匀化截面以及平均中子通量密度、各栅元表面的净中子流密度和中子通量密度；2、根据等效均匀化原理为非均匀泄漏修正的求解方法产生不连续因子；采用纵坐标离散法对非均匀泄漏修正模型进行求解；3、利用非均匀泄漏修正产生的临界能谱修正中子通量密度分布，重新归并得到均匀化截面，并产生相应的扩散系数，并对均匀化截面和扩散系数乘上栅元均匀化组件与非均匀组件的通量比值；本发明可以利用非均匀泄漏修正模型求得各个栅元各自的能谱修正因子以及其泄漏系数。