基于蒙特卡洛程序的快堆少群截面制作方法研究

肖鹏; 罗琦; 夏榜样; 张广春; 姚栋; 周亚婧; 方超; 李天涯

doi:10.13832/j.jnpe.2023.S2.0017

基于蒙特卡洛程序的快堆少群截面制作方法研究

doi: 10.13832/j.jnpe.2023.S2.0017

1.
中国核动力研究设计院核反应堆系统设计技术重点实验室，成都，610213
2.
中国核工业集团有限公司，北京，100045
3.
哈尔滨工业大学能源科学与工程学院，哈尔滨，150001

详细信息

作者简介:
肖　鹏（1988—），男，高级工程师，现主要从事反应堆物理设计与数值计算方法研究，Email：who545@163.com

中图分类号: TP29
计量
- 文章访问数: 75
- HTML全文浏览量: 13
- PDF下载量: 29
- 被引次数: 0
出版历程
- 收稿日期: 2023-07-11
- 修回日期: 2023-08-03
- 刊出日期: 2023-12-30

Research on Few Group Cross-Section Generation Method For Fast Reactor Based on Monte Carlo Code

1.
Science and Technology on Reactor System Design Technology Laboratory, Nuclear Power Institute of China, Chengdu, 610213, China
2.
China National Nuclear Corporation, Beijing, 100045, China
3.
School of Energy Science and Engineering, Harbin Institute of Technology, Harbin, 150001, China

摘要

摘要: 为了提高快堆中子学分析精度，解决传统确定论方法制作快堆少群截面面临的复杂共振问题，基于组件均匀化-堆芯输运计算两步法，从各向异性散射、能群结构、泄漏修正等方面，开展了采用蒙特卡洛程序制作堆芯输运计算所需的组件均匀化少群截面的方法研究。数值结果表明：采用蒙特卡洛程序制作快堆少群截面，可以很好地处理中等质量核素的共振效应；采用1阶散射截面或者输运修正后的少群截面可以有效考虑各向异性散射的影响；少群截面的能群数目减小会导致堆芯有效增殖因子（k_eff）减小；泄漏修正会导致堆芯k_eff增大。采用蒙特卡洛程序制作快堆少群截面，堆芯输运得到的k_eff偏差小于100pcm（1pcm=10⁻⁵），功率分布偏差小于2%。
- 少群截面生成 /
- 蒙特卡洛 /
- 堆芯输运 /
- 快堆
Abstract: In order to improve the accuracy of neutronics analysis for fast reactors and solve the complex resonance problem faced by traditional deterministic method in making few group cross-sections of fast reactors, this paper, based on the two-step method of assembly homogenization-core transport calculation, studies the method of assembly homogenized few group cross-section by using the Monte Carlo code from the aspects of anisotropic scattering, energy group structure and leakage correction. The numerical results show that the resonance effect of medium-mass nuclides can be well treated in the fast reactor few group cross-section generation based on the Monte Carlo code; the effects of anisotropic scattering can be effectively considered with 1-order scattering cross-section or transport modified few group section. Core effective multiplication factor ( k_eff ) may decrease due to the decrease in the number of energy groups of few group cross-section; core k_eff may increase due to the leakage correction. When Monte Carlo code is used in few group cross-sections generation of fast reactor, the k_eff deviation of core transport is less than 100 pcm (1pcm=10⁻⁵), and the power distribution deviation is less than 2%.
- Few group cross-section generation /
- Monte Carlo /
- Core transport /
- Fast reactor

HTML全文

图 1 MET-1000堆芯径向示意图

Figure 1. Schematic Diagram of Radial Core Layout of MET-1000　　　　　

下载: 全尺寸图片幻灯片

图 2 二维六边形精细燃料组件模型

Figure 2. 2D Hexagonal Fine Fuel Assembly Model

下载: 全尺寸图片幻灯片

图 3 二维超栅元精细模型

Figure 3. 2D Super-cell Fine Models

下载: 全尺寸图片幻灯片

图 4 二维径向反射层模型

Figure 4. 2D Model for Radial Reflector

下载: 全尺寸图片幻灯片

图 5 堆芯径向功率偏差分布

Figure 5. Core Radial Power Difference Distribution

下载: 全尺寸图片幻灯片

表 1 MET-1000堆芯主要参数

Table 1. Main Parameters of MET-1000

参数名	参数值
燃料	(U,Pu)Zr
活性区高度/cm	85.82
组件中心距/cm	16.2471
组件壁厚度/cm	0.3966
包壳外半径/cm	0.3857
包壳内半径/cm	0.3236
燃料棒半径/cm	0.3236
结构材料	HT-9钢
冷却剂温度/K	705
结构材料温度/K	705
燃料温度/K	807

下载: 导出CSV

表 2 粒子数对堆芯k_eff的影响

Table 2. Influence of Particle Number on Core k_eff

算例编号	每代粒子数/万	活跃代数/代	k_eff	计算偏差/pcm
1	50	1000	1.03530	−19
2	50	1500	1.03584	33
3	200	2000	1.03653	99
参考解			1.03550
1pcm=10⁻⁵

下载: 导出CSV

表 3 各向异性散射对堆芯k_eff的影响

Table 3. Influence of Anisotropic Scattering on Core k_eff

算例	k_eff	计算偏差/pcm
参考解	1.0355
MOCS 0阶	1.07744	3970
MOCS 1阶	1.02351	−1165
MOCS输运修正	1.02653	−870
Serpent 0阶	1.06844	3132
Serpent 1阶	1.03549	−1
Serpent输运修正	1.03653	99

下载: 导出CSV

表 4 不同能群下堆芯k_eff计算结果对比

Table 4. Comparison of Calculation Results of Core k_eff for Different Energy Groups

算例	k_eff	计算偏差/pcm
参考解	1.03550
111g-24g	1.03652	98
47g-24g	1.03653	99
24g-24g	1.03653	99
111g-9g	1.03590	39
47g-9g	1.03589	38
24g-9g	1.03589	38

下载: 导出CSV

表 5 泄漏效应对堆芯k_eff的影响

Table 5. Influence of Leakage Effect on Core k_eff

算例	k_eff	计算偏差/pcm	泄漏修正影响/pcm
参考解	1.03550
111g-24g	1.03652	98	218
111g-24g-B1	1.03878	316	218
47g-24g	1.03653	99	161
47g-24g-B1	1.03820	260	161
111g-9g	1.03590	39	936
111g-9g-B1	1.04565	975	936
47g-9g	1.03589	38	875
47g-9g-B1	1.04500	913	875
24g-9g	1.03589	38	701
24g-9g-B1	1.04318	739	701

下载: 导出CSV

参考文献(6)

[1]	FRIDMAN E, SHWAGERAUS E. Modeling of SFR cores with Serpent–DYN3D codes sequence[J]. Annals of Nuclear Energy, 2013, 53: 354-363. doi: 10.1016/j.anucene.2012.08.006
[2]	杜夏楠,吴宏春,郑友琦. 蒙特卡罗方法在快堆组件参数计算中的应用[J]. 核动力工程,2014, 35(S2): 67-70. doi: 10.13832/j.jnpe.2014.S2.0067
[3]	Nuclear Energy Agency. Benchmark for neutronic analysis of sodium-cooled fast reactor cores with various fuel types and core sizes: NEA/NSC/R(2015)9[R]. Paris: OECD/NEA, 2016.
[4]	LEPPÄNEN J, PUSA M, FRIDMAN E, Overview of methodology for spatial homogenization in the Serpent 2 Monte Carlo code[J]. Annals of Nuclear Energy, 2016, 96: 126-136.
[5]	ZHANG B, WANG L J, LOU L, et al. Development and verification of lead-bismuth cooled fast reactor calculation code system Mosasaur[J]. Frontiers in Energy Research, 2023, 10: 1055405. doi: 10.3389/fenrg.2022.1055405
[6]	RUGGIERI J M, TOMMASI J, LEBRAT J F, et al. ERANOS 2.1: international code system for GEN IV fast reactor analysis[C]//Proceedings of the 2006 International Congress on Advances in Nuclear Power Plants (ICAPP’06). Reno, Nevada: American Nuclear Society, 2006: 2432-2439.