一种含气泡介质的中子输运计算方法技术

本发明专利技术公开了一种含气泡介质的中子输运计算方法，首先分别计算核反应堆活性区内各个网格中气泡区域和非气泡区域的宏观中子反应截面；根据各网格内气泡尺寸，计算气泡区域和非气泡区域间的中子碰撞概率；利用中子碰撞概率计算网格内气泡区域和非气泡区域的中子通量密度；最后，通过通量

【技术实现步骤摘要】
一种含气泡介质的中子输运计算方法


[0001]本专利技术涉及核反应堆物理分析计算
，具体涉及一种含气泡介质的中子输运计算方法。

技术介绍

[0002]在反应堆正常运行过程中，堆芯活性区不会出现气泡，因此在采用两步法进行物理分析时，组件计算只计算包含固体（如燃料或其他结构材料）和液体（如冷却剂）的均匀化截面。然而，在某些情况下，如反应性严重事故工况，由于裂变反应剧增、反应堆温度超限，堆芯内会产生气泡，气泡的产生改变了局部的材料布置特性，对反应堆临界特性产生较大影响，因此必须考虑这一现象。
[0003]目前，在堆芯输运计算时，若堆芯活性区产生气泡，则根据气泡体积改变网格内核素的平均核子密度，即将气泡均匀打混至液体中去。这一做法相当于体积权重，忽略了气泡区域和非气泡区域的中子通量的差别，无法保证网格内的中子反应率守恒，在中子输运计算中引入了误差。

技术实现思路

[0004]针对现有技术存在的问题，本专利技术的目的在于提供一种含气泡介质的中子输运计算方法，首先分别计算核反应堆活性区网格内气泡区域和非气泡区域的宏观中子反应截面，根据各网格内气泡尺寸计算气泡区域和非气泡区域间的中子碰撞概率，利用中子碰撞概率计算网格内气泡区域和非气泡区域的中子通量密度，通过通量
‑
体积权重获得网格内的均匀化宏观中子反应截面，从而考虑了气泡区域和非气泡区域的通量分布的区别，提供了准确的均匀化中子反应截面，提高了后续中子输运计算的精度。
[0005]为了实现以上目的，本专利技术采取如下的技术方案予以实施：一种含气泡介质的中子输运计算方法，首先分别计算核反应堆活性区网格内气泡区域和非气泡区域的宏观中子反应截面，根据各网格内气泡尺寸计算气泡区域和非气泡区域间的中子碰撞概率，利用中子碰撞概率计算网格内气泡区域和非气泡区域的中子通量密度，通过通量
‑
体积权重获得网格内的均匀化宏观中子反应截面，从而考虑了气泡区域和非气泡区域的通量分布的区别，提供了准确的均匀化中子反应截面，提高了后续中子输运计算的精度，包括如下步骤：步骤1：对于每个计算网格，分别利用公式（1）和公式（2）计算气泡区域和非气泡区域的宏观中子反应截面：
ꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
ꢀꢀꢀꢀꢀ
(2)
式中：—— 气泡区域的宏观中子反应截面；I —— 气泡区域的核素总数；i —— 气泡区域的第i种核素； —— 气泡区域第i种核素的核子密度；—— 气泡区域第i种核素的微观中子反应截面； —— 非气泡区域的宏观中子反应截面；J —— 非气泡区域的核素总数；j —— 气泡区域的第j种核素；—— 非气泡区域第j种核素的核子密度；—— 非气泡区域第j种核素的微观中子反应截面；步骤2：根据各网格内气泡尺寸，计算气泡区域和非气泡区域间的中子碰撞概率；根据各网格内气泡尺寸，计算出气泡的中子穿透概率和中子逃脱概率为：
ꢀꢀꢀꢀꢀꢀꢀ
(3) (4)式中： —— 气泡的中子穿透概率，即中子进入气泡后能飞出气泡的概率；—— 气泡的中子逃脱概率，即气泡内的中子未经碰撞逃脱气泡的概率；—— 气泡区域的宏观中子输运截面；—— 气泡的平均半径；计算气泡区域和非气泡区域间的中子碰撞概率：
ꢀꢀꢀꢀꢀꢀ
(5)
ꢀꢀ
(6)
ꢀꢀꢀꢀꢀ
(7)
ꢀꢀꢀꢀꢀ
(8)式中： —— 气泡区域到非气泡区域的中子碰撞概率； —— 非气泡区域到非气泡区域的中子碰撞概率； —— 气泡区域到气泡区域的中子碰撞概率； —— 非气泡区域到气泡区域的中子碰撞概率； —— 非气泡区域的宏观中子输运截面； —— 非气泡区域的平均弦长；步骤3：利用中子碰撞概率计算网格内气泡区域和非气泡区域的中子通量密度；利用中子碰撞概率，建立以下气泡区域和非气泡区域的中子平衡方程：
ꢀꢀꢀꢀ
(9)
ꢀꢀꢀ
(10)式中： —— 气泡区域的中子通量密度； —— 气泡区域的宏观中子散射截面； —— 非气泡区域的中子通量密度； —— 非气泡区域的宏观中子散射截面； —— 核反应堆的有效中子增殖因子；方程（9）和方程（10）包含的未知量为和，两个方程两个未知数，因此通过求解获得气泡区域的中子通量密度和非气泡区域的中子通量密度；步骤4：通过通量
‑
体积权重对网格中气泡区域和非气泡区域的宏观中子反应截面加权平均获得网格内的均匀化宏观中子反应截面，从而进行后续中子输运计算，获得各网格内的中子通量密度；网格内的均匀化宏观中子反应截面计算式为：
ꢀꢀꢀꢀ
(11)
式中： —— 网格的均匀化宏观中子反应截面； —— 气泡区域的体积； —— 非气泡区域的体积；获得网格内的均匀化宏观中子反应截面后，便采用中子输运计算方法求解中子输运方程，最终获得所有网格内的中子通量密度。
[0006]优选的，获得网格内的均匀化宏观中子反应截面后，利用离散纵标计算方法求解中子输运方程，获得所有网格内的中子通量密度。
[0007]与现有技术相比，本专利技术有如下优点：本专利技术通过计算网格内气泡区域和非气泡区域间的中子碰撞概率，获得了网格内气泡区域和非气泡区域的中子通量密度，进而可以通过体积
‑
通量权重获得网格均匀化截面，获得了保证反应率守恒的高精度的网格均匀化截面，提高了堆芯中子输运计算的精度。
附图说明
[0008]图1为一种含气泡介质的中子输运计算方法总体流程图。
[0009]图2为某一网格内气泡区域和非气泡区域示意图。
具体实施方式
[0010]下面结合附图和具体实施方式对本专利技术做进一步详细说明：如图1所示，本专利技术一种含气泡介质的中子输运计算方法，包括如下步骤：步骤1：对于核反应堆活性区的计算网格，分别利用公式（1）和公式（2）计算气泡区域和非气泡区域的宏观中子反应截面。例如，脉冲堆在大脉冲功率运行工况下，会产生局部高温，导致水中产生气泡，如图2所示，分别计算图2中气泡区域和非气泡区域的宏观中子反应截面。
[0011]ꢀꢀꢀꢀ
(1)
ꢀꢀꢀꢀꢀ
(2)式中： —— 气泡区域的宏观中子反应截面；I —— 气泡区域的核素总数；i —— 气泡区域的第i种核素； —— 气泡区域第i种核素的核子密度； —— 气泡区域第i种核素的微观中子反应截面； —— 非气泡区域的宏观中子反应截面；
J —— 非气泡区域的核素总数；j —— 气泡区域的第j种核素； —— 非气泡区域第j种核素的核子密度； —— 非气泡区域第j种核素的微观中子反应截面；步骤2：中子碰撞概率的计算与气泡尺寸密切相关，因此需要根据各网格内气泡尺寸，计算气泡区域和非气泡区域间的中子碰撞概率；根据各网格内气泡尺寸，计算出气泡的中子穿透概率和中子逃脱概率为：
ꢀꢀꢀ
(3)
ꢀꢀꢀꢀ
(4)式中：—— 气泡的中子穿透概率，即中子进入气泡后能飞出气泡的概率；—— 气泡的中子逃脱概率，即气泡内的中子未经碰撞逃脱本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种含气泡介质的中子输运计算方法，其特征在于：首先分别计算核反应堆活性区网格内气泡区域和非气泡区域的宏观中子反应截面，根据各网格内气泡尺寸计算气泡区域和非气泡区域间的中子碰撞概率，利用中子碰撞概率计算网格内气泡区域和非气泡区域的中子通量密度，通过通量
‑
体积权重获得网格内的均匀化宏观中子反应截面，包括如下步骤：步骤1：对于每个计算网格，分别利用公式（1）和公式（2）计算气泡区域和非气泡区域的宏观中子反应截面：
ꢀꢀꢀ
(1)
ꢀꢀꢀꢀꢀꢀꢀ
(2)式中： —— 气泡区域的宏观中子反应截面；I —— 气泡区域的核素总数；i —— 气泡区域的第i种核素； —— 气泡区域第i种核素的核子密度； —— 气泡区域第i种核素的微观中子反应截面；—— 非气泡区域的宏观中子反应截面；J —— 非气泡区域的核素总数；j —— 气泡区域的第j种核素； —— 非气泡区域第j种核素的核子密度； —— 非气泡区域第j种核素的微观中子反应截面；步骤2：根据各网格内气泡尺寸，计算气泡区域和非气泡区域间的中子碰撞概率；根据各网格内气泡尺寸，计算出气泡的中子穿透概率和中子逃脱概率为：
ꢀꢀꢀꢀꢀ
(3)
ꢀꢀꢀꢀꢀ
(4)式中：—— 气泡的中子穿透概率，即中子进入气泡后能飞出气泡的概率；—— 气泡的中子逃脱概率，即气泡内的中子未经碰撞逃脱气泡的概率；—— 气泡区域的宏观中子输运截面； —— 气泡的平均半径；
计算气泡区域和非气泡区域间的中子碰撞概率：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
ꢀꢀꢀꢀ<...

【专利技术属性】
技术研发人员：王永平，郑友琦，杜夏楠，吴宏春，曹良志，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：