CN107273582B - 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 - Google Patents
一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 Download PDFInfo
- Publication number
- CN107273582B CN107273582B CN201710368718.XA CN201710368718A CN107273582B CN 107273582 B CN107273582 B CN 107273582B CN 201710368718 A CN201710368718 A CN 201710368718A CN 107273582 B CN107273582 B CN 107273582B
- Authority
- CN
- China
- Prior art keywords
- burnup
- neutron
- group
- nuclide
- reactor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
一种用于快中子反应堆中子输运燃耗耦合分析的计算方法,1、将堆芯径向结构剖分为三角形网格,建立三棱柱网格,以组件轴向分段为单位划分燃耗区,将堆内循环划分为若干燃耗步,对各燃耗步执行如下流程;2、计算各燃耗区少群宏观截面,采用基于三棱柱网格的中子输运计算方法进行堆芯中子输运计算;3、计算燃耗步初时各燃耗区的燃耗矩阵,采用切比雪夫有理近似方法求解燃耗方程;4、根据燃耗步末时各燃耗区核子密度向量,进行堆芯中子输运计算;5、对燃耗步初和末时的燃耗矩阵求平均得到各燃耗区的平均燃耗矩阵,对各燃耗区重新从燃耗步初进行燃耗计算;6、重复步骤4和5,直至相邻两次计算得到的燃耗步末时各燃耗区的核子密度向量收敛。
Description
技术领域
本发明涉及核反应堆物理计算和设计分析领域,是一种用于快中子反应堆中子输运燃耗耦合分析的计算方法。
背景技术
随着核能的不断发展和应用需求的不断扩大,新型快中子反应堆设计方案被不断地提出,堆芯内部不再是单一的、规则的正方形或六边形组件排布,快中子反应堆的设计分析对堆芯中子输运燃耗耦合分析方法的精度和效率均提出新的挑战。现有的快中子反应堆中子输运燃耗耦合分析方法主要基于六边形节块的中子输运或扩散计算方法,并采用线性链解析燃耗算法或矩阵指数燃耗算法用于堆芯燃耗计算,中子输运计算和燃耗计算的耦合方式为子步法或预估校正法。
现有的快中子反应堆中子输运燃耗耦合计算方法主要存在以下两方面问题:第一,只适用于堆芯规则几何结构,无法处理具有不规则堆芯布置的新型快中子反应堆的设计分析;第二,以子步法或预估校正法对中子输运计算和燃耗计算进行耦合,在燃耗步长较大时,计算精度有所损失。
因此,亟需研究一种快中子反应堆中子输运燃耗耦合分析的计算方法,模拟具有规则及不规则堆芯布置的快中子反应堆堆内燃耗循环过程,并可在采用较大的燃耗步长时保证计算精度,提高计算效率。
发明内容
为了解决上述快中子反应堆中子输运燃耗耦合计算方法存在的问题,本发明提出了一种用于快中子反应堆中子输运燃耗耦合分析的计算方法,可用于具有规则及不规则堆芯布置的快中子反应堆堆内燃耗循环过程的模拟,并可在采用较大的燃耗步长时保证计算精度,提高计算效率。
为了达到上述目的,本发明的技术方案如下:
一种用于快中子反应堆中子输运燃耗耦合分析的计算方法,包括如下步骤:
步骤1:将快中子反应堆的径向几何结构剖分为三角形网格,通过轴向分层,建立快中子反应堆的三棱柱空间网格,并以组件轴向分段为单位确定快中子反应堆的燃耗区网格;将快中子反应堆堆内燃耗循环过程按时间划分为若干个燃耗步,对每个燃耗步依次执行如下的步骤2至步骤6中子输运燃耗耦合分析流程;
步骤2:根据当前燃耗步初时堆芯各燃耗区内各燃耗核素的核子密度,结合各燃耗核素的少群微观截面,计算得到堆芯中子输运计算需要的各燃耗区少群宏观截面如公式(1)-(3)所示;
式中:
ν∑f,g——第g群宏观中子产生截面;
νσf,g,i——燃耗核素i的第g群微观中子产生截面;
Ni——燃耗核素i的核子密度;
式中:
∑s,g′→g——燃耗区从第g′群到第g群的宏观中子散射截面;
σs,g′→g,i——燃耗核素i的从第g′群到第g群的微观中子散射截面;
Ni——燃耗核素i的核子密度;
式中:
∑t,g——燃耗区第g群宏观总截面;
σa,g,i——燃耗核素i的第g群中子吸收截面;
Ni——燃耗核素i的核子密度;
∑s,g→g′——燃耗区从第g群到第g′群的宏观中子散射截面;
采用基于三棱柱空间网格的中子输运计算方法对快中子反应堆进行中子输运计算,如公式(4)所示将堆芯功率归一化至额定功率,计算得到堆芯各燃耗区内的中子通量密度分布;
式中:
P——堆芯额定功率;
f——堆芯功率归一化因子;
Vk——燃耗区k的体积;
k——平均每次裂变在堆芯内沉积的能量;
∑f,g,k——燃耗区k的第g群宏观裂变截面;
φg,k——燃耗区k的第g群中子通量密度;
步骤3:对堆芯各燃耗区分别进行如下本步骤内的燃耗计算,如公式(5)所示,由当前燃耗区内各燃耗核素的微观截面及中子通量密度,计算得到燃耗计算需要的描述各燃耗核素间相互转化关系的燃耗矩阵;如公式(6)所示,采用切比雪夫有理近似方法以计算矩阵指数的方式求解燃耗方程,得到当前燃耗步末时当前燃耗区内各燃耗核素的核子密度;
式中:
Aij——燃耗矩阵A的第(i,j)个元素;
φg——第g群中子通量密度;
γij——核素j经衰变反应产生i的份额;
λi——核素i的衰变常数/s-1;
λj——核素j的衰变常数/s-1;
式中:
N(t)——燃耗步末时的核子密度向量;
A——燃耗矩阵;
t——燃耗时间步长;
N(0)——燃耗步初时的核子密度向量;
a0——函数在趋于无穷时的极限值;
Re——复数取实部的函数;
k——切比雪夫有理近似式的展开阶数;
aj——在极点θj的留数;
I——单位矩阵;
步骤4:根据当前燃耗步末时堆芯各燃耗区内各燃耗核素的核子密度,重新计算得到堆芯各燃耗区少群宏观截面,并进行基于三棱柱空间网格的中子输运计算,计算得到燃耗步末时堆芯各燃耗区内的中子通量密度分布;
步骤5:由燃耗步末时各燃耗区内各燃耗核素的微观截面及中子通量密度建立燃耗矩阵,并由公式(7)对燃耗步初和末时的燃耗矩阵求平均得到各燃耗区的平均燃耗矩阵,利用平均燃耗矩阵重新对各燃耗区从燃耗步初进行燃耗计算至燃耗步末;
式中:
tp——燃耗步初的时刻;
tp+1——燃耗步末的时刻;
步骤6:重复执行步骤4和5,直至相邻两次计算得到的燃耗步末时各燃耗区内各燃耗核素的核子密度收敛为止;收敛的判断准则为相邻两次计算得到的燃耗步末时各燃耗区核子密度向量的相对误差2范数e均小于用户声明的收敛准则值,如公式(8)所示;
e=||N(q)-N(q-1)||2 公式(8)
e——相邻两次计算得到的燃耗步末时核子密度向量的相对误差2范数;
N(q)——第q次计算得到的燃耗步末时的核子密度向量;
N(q-1)——第q-1次计算得到的燃耗步末时的核子密度向量;
与传统的基于六边形结构网格的中子输运或扩散燃耗耦合计算方法相比,本发明有如下突出优点:
1.本发明在堆芯中子输运计算时采用三棱柱空间网格对快中子反应堆中出现的结构和非结构几何进行逼近,具有良好的几何适应性;
2.通过将堆芯中子输运计算与燃耗计算过程进行迭代求解,在选择较大燃耗步长时保证计算结果的精度不受损失,提高计算效率。
附图说明
图1快中子反应堆三棱柱网格剖分示意图。
图2快中子反应堆中子输运燃耗耦合分析流程图。
具体实施方式
本发明基于三棱柱空间网格的中子输运计算方法以及切比雪夫有理近似的矩阵指数燃耗算法,提出了一套基于三棱柱空间网格的中子输运燃耗耦合计算方法,可用于具有规则及不规则堆芯布置的快中子反应堆堆内燃耗循环过程的模拟,并可在采用较大的燃耗步长时保证计算精度,提高计算效率。
本发明包括以下方面:
1)将快中子反应堆的堆芯径向结构剖分为三角形网格,并通过轴向分层,建立快中子反应堆的三棱柱空间网格;
2)基于三棱柱空间网格的中子输运计算方法,计算得到快中子反应堆各燃耗区的中子通量密度分布;
3)对各燃耗区,根据燃耗步初和末时的燃耗矩阵计算平均燃耗矩阵,并由切比雪夫有理近似的矩阵指数燃耗算法进行燃耗计算,得到燃耗步末时的核子密度向量;
4)对堆芯中子输运计算与燃耗计算过程进行迭代求解,直到各燃耗区在燃耗步末时的核子密度向量收敛为止。
如图2所示,具体实施方式如下:
步骤1:如图1所示,使用针对二维几何结构的三角形网格剖分工具,如通用有限元分析软件ANSYS的几何预处理工具,将快中子反应堆的径向几何结构剖分为三角形网格,通过轴向分层,建立快中子反应堆的三棱柱空间网格,并以组件轴向分段为单位确定快中子反应堆的燃耗区网格;将快中子反应堆堆内燃耗循环过程按时间划分为若干个燃耗步,对每个燃耗步执行如下的中子输运燃耗耦合分析流程;
步骤2:根据该燃耗步初时堆芯各燃耗区内各燃耗核素的核子密度,结合各燃耗核素的少群微观截面,计算得到堆芯中子输运计算需要的各燃耗区少群宏观截面如公式(1)-(3)所示;
式中:
ν∑f,g——第g群宏观中子产生截面;
νσf,g,i——燃耗核素i的第g群微观中子产生截面;
Ni——燃耗核素i的核子密度;
式中:
∑s,g′→g——燃耗区从第g′群到第g群的宏观中子散射截面;
σs,g′→g,i——燃耗核素i的从第g′群到第g群的微观中子散射截面;
Ni——燃耗核素i的核子密度;
式中:
∑t,g——燃耗区第g群宏观总截面;
σa,g,i——燃耗核素i的第g群中子吸收截面;
Ni——燃耗核素i的核子密度;
∑s,g→g′——燃耗区从第g群到第g′群的宏观中子散射截面;
采用基于三棱柱空间网格的中子输运计算方法,如基于横向积分的离散纵标节块输运计算方法,对快中子反应堆进行堆芯中子输运计算,如公式(4)所示将堆芯功率归一化至额定功率,计算得到堆芯各燃耗区内的中子通量密度分布;
式中:
P——堆芯额定功率;
f——堆芯功率归一化因子;
Vk——燃耗区k的体积;
k——平均每次裂变在堆芯内沉积的能量;
∑f,g,k——燃耗区k的第g群宏观裂变截面;
φg,k——燃耗区k的第g群中子通量密度;
步骤3:对堆芯各燃耗区分别进行如下的燃耗计算,如公式(5)所示,由该燃耗区内各燃耗核素的微观截面及中子通量密度,计算得到燃耗计算需要的描述各燃耗核素间相互转化关系的燃耗矩阵;如公式(6)所示,采用切比雪夫有理近似方法以计算矩阵指数的方式求解燃耗方程,得到该燃耗步末时该燃耗区内各燃耗核素的核子密度;
式中:
Aij——燃耗矩阵A的第(i,j)个元素;
φg——第g群中子通量密度;
γij——核素j经衰变反应产生i的份额;
λi——核素i的衰变常数/s-1;
λj——核素j的衰变常数/s-1;
式中:
N(t)——燃耗步末时的核子密度向量;
A——燃耗矩阵;
t——燃耗时间步长;
N(0)——燃耗步初时的核子密度向量;
α0——函数在趋于无穷时的极限值;
Re——复数取实部的函数;
k——切比雪夫有理近似式的展开阶数;
αj——在极点θj的留数;
I——单位矩阵;
步骤4:根据该燃耗步末时堆芯各燃耗区内各燃耗核素的核子密度,重新计算得到堆芯各燃耗区少群宏观截面,并进行基于三棱柱空间网格的中子输运计算,计算得到燃耗步末时堆芯各燃耗区内的中子通量密度分布;
步骤5:由燃耗步末时各燃耗区内各燃耗核素的微观截面及中子通量密度建立燃耗矩阵,并由公式(7)对燃耗步初和末时的燃耗矩阵求平均得到各燃耗区的平均燃耗矩阵,利用平均燃耗矩阵重新对各燃耗区从燃耗步初进行燃耗计算至燃耗步末;
式中:
tp——燃耗步初的时刻;
tp+1——燃耗步末的时刻;
步骤6:重复执行步骤4和5,直至相邻两次计算得到的燃耗步末时各燃耗区内各燃耗核素的核子密度收敛为止。收敛的判断准则为相邻两次计算得到的燃耗步末时各燃耗区核子密度向量的相对误差2范数e均小于用户声明的收敛准则值,如公式(8)所示。
e=||N(q)-N(q-1)||2 公式(8)
e——相邻两次计算得到的燃耗步末时核子密度向量的相对误差2范数;
N(q)——第q次计算得到的燃耗步末时的核子密度向量;
N(q-1)——第q-1次计算得到的燃耗步末时的核子密度向量。
Claims (1)
1.一种用于快中子反应堆中子输运燃耗耦合分析的计算方法,其特征在于:包括如下步骤:
步骤1:将快中子反应堆的径向几何结构剖分为三角形网格,通过轴向分层,建立快中子反应堆的三棱柱空间网格,并以组件轴向分段为单位确定快中子反应堆的燃耗区网格;将快中子反应堆堆内燃耗循环过程按时间划分为若干个燃耗步,对每个燃耗步依次执行如下的步骤2至步骤6中子输运燃耗耦合分析流程;
步骤2:根据当前燃耗步初时堆芯各燃耗区内各燃耗核素的核子密度,结合各燃耗核素的少群微观截面,计算得到堆芯中子输运计算需要的各燃耗区少群宏观截面如公式(1)-(3)所示;
式中:
νΣf,g——第g群宏观中子产生截面;
νσf,g,i——燃耗核素i的第g群微观中子产生截面;
Ni——燃耗核素i的核子密度;
式中:
Σs,g′→g——燃耗区从第g′群到第g群的宏观中子散射截面;
σs,g′→g,i——燃耗核素i的从第g′群到第g群的微观中子散射截面;
Ni——燃耗核素i的核子密度;
式中:
Σt,g——燃耗区第g群宏观总截面;
σa,g,i——燃耗核素i的第g群中子吸收截面;
Ni——燃耗核素i的核子密度;
Σs,g→g′——燃耗区从第g群到第g′群的宏观中子散射截面;
采用基于三棱柱空间网格的中子输运计算方法对快中子反应堆进行中子输运计算,如公式(4)所示将堆芯功率归一化至额定功率,计算得到堆芯各燃耗区内的中子通量密度分布;
式中:
P——堆芯额定功率;
f——堆芯功率归一化因子;
Vk——燃耗区k的体积;
κ——平均每次裂变在堆芯内沉积的能量;
Σf,g,k——燃耗区k的第g群宏观裂变截面;
φg,k——燃耗区k的第g群中子通量密度;
步骤3:对堆芯各燃耗区分别进行如下本步骤内的燃耗计算,如公式(5)所示,由当前燃耗区内各燃耗核素的微观截面及中子通量密度,计算得到燃耗计算需要的描述各燃耗核素间相互转化关系的燃耗矩阵;如公式(6)所示,采用切比雪夫有理近似方法以计算矩阵指数的方式求解燃耗方程,得到当前燃耗步末时当前燃耗区内各燃耗核素的核子密度;
式中:
Aij——燃耗矩阵A的第(i,j)个元素;
φg——第g群中子通量密度;
γij——核素j经衰变反应产生i的份额;
λi——核素i的衰变常数,单位为每秒;
λj——核素j的衰变常数,单位为每秒;
式中:
N(t)——燃耗步末时的核子密度向量;
A——燃耗矩阵;
t——燃耗时间步长;
N(0)——燃耗步初时的核子密度向量;
α0——函数在趋于无穷时的极限值;
Re——复数取实部的函数;
k——切比雪夫有理近似式的展开阶数;
αj——在极点θj的留数;
I——单位矩阵;
步骤4:根据当前燃耗步末时堆芯各燃耗区内各燃耗核素的核子密度,重新计算得到堆芯各燃耗区少群宏观截面,并进行基于三棱柱空间网格的中子输运计算,计算得到燃耗步末时堆芯各燃耗区内的中子通量密度分布;
步骤5:由燃耗步末时各燃耗区内各燃耗核素的微观截面及中子通量密度建立燃耗矩阵,并由公式(7)对燃耗步初和末时的燃耗矩阵求平均得到各燃耗区的平均燃耗矩阵,利用平均燃耗矩阵重新对各燃耗区从燃耗步初进行燃耗计算至燃耗步末;
式中:
tp——燃耗步初的时刻;
tp+1——燃耗步末的时刻;
步骤6:重复执行步骤4和5,直至相邻两次计算得到的燃耗步末时各燃耗区内各燃耗核素的核子密度收敛为止;收敛的判断准则为相邻两次计算得到的燃耗步末时各燃耗区核子密度向量的相对误差2范数e均小于用户声明的收敛准则值,如公式(8)所示;
e=||N(q)-N(q-1)||2 公式(8)
e——相邻两次计算得到的燃耗步末时核子密度向量的相对误差2范数;
N(q)——第q次计算得到的燃耗步末时的核子密度向量;
N(q-1)——第q-1次计算得到的燃耗步末时的核子密度向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710368718.XA CN107273582B (zh) | 2017-05-23 | 2017-05-23 | 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710368718.XA CN107273582B (zh) | 2017-05-23 | 2017-05-23 | 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107273582A CN107273582A (zh) | 2017-10-20 |
CN107273582B true CN107273582B (zh) | 2020-03-17 |
Family
ID=60065630
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710368718.XA Active CN107273582B (zh) | 2017-05-23 | 2017-05-23 | 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107273582B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108664714B (zh) * | 2018-04-26 | 2022-03-29 | 华南理工大学 | 克雷洛夫子空间加速求解燃耗方程的数值计算方法 |
CN108846190B (zh) * | 2018-06-05 | 2022-04-12 | 哈尔滨工程大学 | 一种压水堆燃料组件的核热耦合仿真方法 |
CN109086506B (zh) * | 2018-07-20 | 2020-08-28 | 西安交通大学 | 一种适用于液态燃料熔盐堆的燃耗分析计算方法 |
FR3085788B1 (fr) * | 2018-09-12 | 2020-11-27 | Framatome Sa | Procede de protection d'un reacteur nucleaire et reacteur nucleaire correspondant |
CN109522510B (zh) * | 2018-11-19 | 2020-05-22 | 西安交通大学 | 一种熔盐堆氚输运特性耦合计算方法 |
CN110717254A (zh) * | 2019-09-21 | 2020-01-21 | 哈尔滨工程大学 | 一种快速预测数字反应堆中子通量分布的混合预处理方法 |
CN111950177B (zh) * | 2020-07-22 | 2024-02-09 | 核工业西南物理研究院 | 一种固态产氚包层的多物理场耦合中子学自动优化方法 |
CN112632771A (zh) * | 2020-12-21 | 2021-04-09 | 中国核动力研究设计院 | 一种求解燃耗计算响应灵敏度的方法及设备 |
CN113312791B (zh) * | 2021-06-17 | 2022-02-22 | 中国核动力研究设计院 | 一种基于sp3方程的反应堆中子噪声分析方法及系统 |
CN113536580B (zh) * | 2021-07-22 | 2023-06-27 | 中国核动力研究设计院 | 确定核反应堆考验回路功率和中子通量密度的方法和系统 |
CN113504993B (zh) * | 2021-07-23 | 2023-06-06 | 中国核动力研究设计院 | 基于有效共振截面的控制棒尖齿效应处理方法及系统 |
CN114003856B (zh) * | 2021-11-02 | 2022-08-05 | 西安交通大学 | 一种核热推进反应堆停堆状态外环境辐射场计算方法 |
CN114510677B (zh) * | 2022-01-19 | 2024-06-11 | 西北核技术研究所 | 基于间断有限元的中子输运方程处理方法、计算机程序产品 |
CN114913936B (zh) * | 2022-07-18 | 2022-10-25 | 西安交通大学 | 一种针对铀钚混合氧化物燃料的多物理燃料性能分析方法 |
CN115130330B (zh) * | 2022-08-30 | 2022-12-09 | 西安交通大学 | 一种含气泡介质的中子输运计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020468A (zh) * | 2012-12-26 | 2013-04-03 | 中山大学 | 核反应堆核热耦合计算方法 |
CN103150424A (zh) * | 2013-02-05 | 2013-06-12 | 西安交通大学 | 一种获取反应堆堆芯三维中子通量密度精细分布的方法 |
CN103294899A (zh) * | 2013-05-10 | 2013-09-11 | 西安交通大学 | 一种计算小型实验反应堆堆芯中子通量分布的方法 |
CN104021278A (zh) * | 2014-05-16 | 2014-09-03 | 中国核动力研究设计院 | 一种反应堆内非燃料可燃毒物燃耗的计算方法 |
CN105426659A (zh) * | 2015-10-30 | 2016-03-23 | 西安交通大学 | 获取不同燃耗下有效增殖因子对截面的灵敏度系数的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120027150A1 (en) * | 2010-08-02 | 2012-02-02 | Ryoichi Wada | High flux fast neutron generator |
-
2017
- 2017-05-23 CN CN201710368718.XA patent/CN107273582B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103020468A (zh) * | 2012-12-26 | 2013-04-03 | 中山大学 | 核反应堆核热耦合计算方法 |
CN103150424A (zh) * | 2013-02-05 | 2013-06-12 | 西安交通大学 | 一种获取反应堆堆芯三维中子通量密度精细分布的方法 |
CN103294899A (zh) * | 2013-05-10 | 2013-09-11 | 西安交通大学 | 一种计算小型实验反应堆堆芯中子通量分布的方法 |
CN104021278A (zh) * | 2014-05-16 | 2014-09-03 | 中国核动力研究设计院 | 一种反应堆内非燃料可燃毒物燃耗的计算方法 |
CN105426659A (zh) * | 2015-10-30 | 2016-03-23 | 西安交通大学 | 获取不同燃耗下有效增殖因子对截面的灵敏度系数的方法 |
Non-Patent Citations (3)
Title |
---|
Youqi Zheng 等.Application of the wavelet expansion method in spatial-angular discretization of the neutron transport equation.《Annals of Nuclear Energy》.2012, * |
周生诚 等.加速器驱动次临界堆芯中子学程序开发.《原子能科学技术》.2013, * |
张延云 等.基于蒙特卡罗方法与指数欧拉法耦合的反应堆三维燃耗计算程序的研发.《核科学与工程》.2011, * |
Also Published As
Publication number | Publication date |
---|---|
CN107273582A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107273582B (zh) | 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 | |
CN111414722B (zh) | 一种核反应堆堆芯物理与热工耦合的模拟方法 | |
CN105653869B (zh) | 一种超临界水堆堆芯三维瞬态性能分析方法 | |
Liu et al. | BEAVRS full core burnup calculation in hot full power condition by RMC code | |
Ivanov et al. | High fidelity simulation of conventional and innovative LWR with the coupled Monte-Carlo thermal-hydraulic system MCNP-SUBCHANFLOW | |
Kochunas et al. | VERA core simulator methodology for PWR cycle depletion | |
CN112906271A (zh) | 一种反应堆瞬态物理热工全耦合精细数值模拟方法及系统 | |
CN115565617B (zh) | 基于能谱环境的超钚同位素生产效率快速评估方法 | |
CN112906272B (zh) | 一种反应堆稳态物理热工全耦合精细数值模拟方法及系统 | |
CN107122564B (zh) | 一种压水堆pin-by-pin计算少群常数参数化方法 | |
Attom et al. | Neutronic analysis of thorium MOX fuel blocks with different driver fuels in advanced block-type HTRs | |
Leppänen | SERPENT Monte Carlo reactor physics code | |
Shi et al. | Improvement of sensitivity and uncertainty analysis capabilities of generalized response in Monte Carlo code RMC | |
Fiorina et al. | Creation of an OpenFOAM fuel performance class based on FRED and integration into the GeN-foam multi-physics code | |
Gougar et al. | Multiscale analysis of pebble bed reactors | |
Barragán-Martínez et al. | MCNPX and HELIOS-2 comparison for the neutronics calculations of a Supercritical Water Reactor HPLWR | |
Liu et al. | Reaction rate tally and depletion calculation with on-the-fly temperature treatment | |
Hakim et al. | Source Convergence in Monte-Carlo Criticality Simulation of CANDU-6 Reactor | |
Bao | Development of the model for the multi-physics analysis of Molten Salt Reactor Experiment using GeN-Foam code | |
Labarile et al. | Comparison of SERPENT and SCALE methodology for LWRs transport calculations and additionally uncertainty analysis for cross-section perturbation with SAMPLER module | |
Liu et al. | Burnup-dependent core neutronics analysis of plate-type research reactor using deterministic and stochastic methods | |
Jiménez et al. | Comparative analysis of neutronics/thermal-hydraulics multi-scale coupling for LWR analysis | |
Sohail et al. | Improved multi-group cross sections with resonance treatment for future applications involving minor actinides | |
Zhang et al. | Adjoint-based global variance reduction approach for reactor analysis problems | |
Kępisty et al. | Assesment of advanced step models for steady state Monte Carlo burnup calculations in application to prismatic HTGR |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |