基于组件计算的燃耗实验基准题建模分析

组件计算在堆芯核设计中占有重要地位。组件程序的燃耗计算精度对核反应堆堆芯的功率分布、换料寿期及反应性控制设计方面具有重要意义。为了评估用于堆芯燃耗计算的多群常数库的准确性,本文运用DRAGON计算程序建立了燃耗实验计算模型,采用SFCOMPO-2.0燃耗实验基准题提供的乏燃料样品燃耗历史参数及最终核素组分信息,对Takahama-3反应堆、H.B. Robinson-2反应堆及Beznau-1反应堆系列样品进行了建模计算,并将计算结果与SFCOMPO-2.0数据库中的基准实验结果进行了对比和分析。结果表明:多数核素的模拟结果与基准值符合良好,误差在10%以内。同时本文对理论计算值与基准实验值存在差异较大的几种核素进行了相应讨论,并对样品计算结果进行了对比分析。     

PWR堆芯中弥散型可燃毒物的燃耗特性研究

在压水反应堆(PWR)堆芯核设计中,通常采用可燃毒物来补偿反应性和展平功率分布。对于长寿期堆芯设计,可燃毒物的消耗和燃料燃耗的匹配研究更为重要。利用基于蒙特卡罗方法开发的堆芯燃耗计算程序(MOI)对天然元素、人工核素、可溶硼等多种弥散型可燃毒物进行燃耗特性分析。结果表明锕系可燃核素231Pa、240Pu等弥散型可燃毒物可用于长寿期PWR的设计。     

组件燃耗对堆芯裂变中子源参数的影响

在中子输运计算中，准确的裂变中子源是保证计算结果可靠性的前提和基础：随着组件燃耗的加深．^235U等核素不断消耗，^239Pu、^240Pu、^241Pu等核素不断积累，导致每次裂变产生的平均中子数v和裂变释放的平均能量E也随之变化，裂变中子源的归一化因子随燃耗增加而增加。因此，在中子输运计算中，必须考虑组件燃耗对堆芯的裂变中子源参数的影响。     

基于蒙特卡罗程序JMCT模拟计算堆芯物理基准题VERA

JMCT是基于蒙特卡罗方法的中子输运程序，具有几何建模精细、截面数据准确、物理过程真实等特点，计算结果具有更高的精度。通过开发临界硼浓度搜索、输运-热工-燃耗耦合功能模块，使JMCT具备了堆芯物理计算功能。本文利用JMCT程序对CASL项目提出的堆芯物理基准题库VERA进行模拟，获得了k_eff、控制棒价值、反应性系数等启动物理参数以及硼降曲线等堆芯运行参数。JMCT计算结果与蒙特卡罗程序KENO-Ⅵ以及MC21进行了对比，结果符合良好，证明了JMCT具有堆芯物理计算能力，并具有较高的精度。     

燃耗截面基本库对输运燃耗耦合计算的影响分析

基于蒙特卡罗中子输运程序和ORIGEN2点燃耗程序的蒙特卡罗输运燃耗耦合计算方法应用广泛。但现有评价库中子连续截面的核素个数远小于燃耗计算涉及到的核素数量,即通过输运计算得到的燃耗截面不足以完全替代燃耗计算的基本库。采用经过栅元验证的蒙特卡罗燃耗程序MCBMPI,对最新的VERA燃耗计算基准题进行验证计算,对比分析不同的燃耗截面基本库对输运燃耗计算的影响。分析结果表明:1)在实际应用中尽量不要采用典型热中子截面库,会带来较大偏差;2)在燃耗计算核素替换较多的情况下,对该基准题而言,选取典型压水堆基本库还是典型快堆基本库,对结果影响不大,二者keff偏差在8‰以内,燃耗末期235U偏差在4‰以内,135Xe偏差在5‰左右;3)建议选取与研究对象能谱相近的基本库。     

回溯算法在燃耗计算中的应用

采用解析解的线性核素链方法在进行燃耗计算时首先根据分治和递归的策略将燃耗矩阵进行解耦,形成具有马尔科夫特性的线性核素链。然后通过对每一条链的解析计算得到所有相关核素的核密度、活度、衰变热等数据。然而在核素链的构建过程中需对每一个可能的核反应路径进行计算。欲保证计算的精度和效率,需寻求一种既能覆盖所有反应路径、又能根据问题描述和约束条件进行自动搜索的算法。本文通过对各种搜索算法的分析和比较并根据燃耗链构建过程的特点,最终采用回溯算法进行深度优先搜索,在搜索过程中完成燃耗链的构建和计算,从而形成问题相关的、具有高精度的自适应燃耗算法。同时结合燃耗过程和回溯算法的特点进行了解空间和时间复杂性的分析。将所开发的多群点燃耗计算程序与蒙特卡罗输运计算程序MCMG-Ⅱ进行耦合,通过对中国实验快堆首炉堆芯燃耗的计算和分析完成程序的初步验证。     

基于蒙特卡罗程序的AP1000反应堆堆芯布置方案的模拟

AP1000是美国西屋公司研发的大型压水反应堆,采用先进的非能动安全系统。AP1000反应堆有两种堆芯燃料布置方案:D19和Adv。结合两种设计方案的优点提出了一种新的堆芯燃料布置方案。利用MCNP6(Monte Carlo N-particle 6)程序对D19堆芯和新方案堆芯的首循环进行建模,并主要计算了新堆芯的核设计参数随燃耗的变化。结果表明,新堆芯在首循环寿期内满足AP1000的主要核设计准则。通过大规模并行计算表明,带燃耗计算功能的蒙特卡罗程序MCNP6能够在堆芯设计工作中发挥很好的参考作用。     

SRAC程序在中国先进研究堆上的应用研究

SRAC程序是由日本原子能研究机构发布的在Linux系统下运行的堆芯物理计算程序包,包括栅元计算程序、堆芯计算程序和燃耗计算程序.本程序采用基于JENDL或ENDF/B系列的数据库,可处理300多种核素的截面参数.该程序具有计算速度快、堆芯燃料管理方便等优点.本文对SRAC程序在中国先进研究堆(CARR)上的应用进行了初步研究,利用SRAC程序对CARR进行了临界计算和首炉燃耗计算.通过与WIMS-CITATION程序系统计算结果的比较,初步确定了SRAC程序是CARR在线计算工具很好的选择,并为今后的实际应用奠定了良好的技术基础.     

铅基快堆RBEC-M燃耗库制作方法研究及验证

      提出了一套新的方法流程，用来处理和生成燃耗计算所需的数据。利用核数据处理程序NJOY处理评价数据库ENDF-B-Ⅶ.1生成33群的MATXS格式库，再根据具体问题中的材料信息，经截面处理程序MGGC处理得到相关核素的微观、宏观截面，经自编写的处理模块Triso对其进行格式转化、合并，最终得到提供给燃耗计算程序使用的ISOTXS库文件，其中一般核素以微观截面的形式表示，裂变产物以类似宏观截面的伪裂变产物形式表示。对铅冷快堆基准题900 MW RBEC-M进行了计算，采用REBUS-3进行燃耗计算，对比了结果中的有效增殖系数k_eff随燃耗的变化趋势、功率分布以及中子能谱，最终结果与参考报告较为符合，初步验证了这一系列燃耗库制作流程的可行性。      

蒙特卡罗燃耗计算程序MCNTRANS的开发与验证

本文介绍了开发的蒙特卡罗燃耗计算程序MCNTRANS。MCNTRANS的中子学计算参数直接采用MCNP5程序的反应率计算值,燃耗计算方法采用图论算法跟踪燃耗链,同时,对实际燃耗过程进行详细分析以提高计算精度与程序适用性,并使用预估 校正方法以获取较大的燃耗计算步长。程序计算结果通过OECD/NEA与JAERI燃耗基准题实验结果进行验证,并与其他程序的计算结果进行比较。结果表明,MCNTRANS程序在不同燃耗深度下的计算结果和实验值与其他程序的计算值符合较好,部分锕系核素与裂变产物的计算精度更高。     

聚变-裂变混合能源堆一维计算模型燃耗分析

聚变-裂变混合能源堆包括聚变中子源和以天然铀为燃料、水为冷却剂的次临界包层,主要目标是生产电力。利用输运燃耗耦合程序系统MCORGS计算了混合能源堆一维模型的燃耗,给出了中子有效增殖因数k_eff、能量放大倍数M、氚增殖比TBR等物理量随时间的变化。通过分析能谱和重要核素随燃耗时间的变化,说明混合能源堆与核燃料增殖、核废料嬗变混合堆的不同特点。本文给出的结果可作为混合堆中子输运、燃耗分析程序校验的参考数据,为混合堆概念研究提供了基础数据。     

Burnup Calculation Sensitivity and Uncertainty Analysis Method of Fission Yield and Decay Half-life

The calculation model of sensitivity coefficient for decay half-life and fission product yield in burnup calculation was derived based on generalized perturbation theory, which considered the interaction between nuclear concentration and neutron flux. A code was developed to calculate sensitivity and uncertainty of effective neutron multiplication factors and nuclide concentration caused by nuclear data. Covariance matrix of fission yield for a simplified burnup library was generated based on standard deviation data of independent fission yield in evaluated nuclear data library to improve the accuracy of uncertainty quantification. Uncertainties induced by decay half-life and fission yield on infinite neutron multiplication factors and nuclide concentration for TMI-1 pin-cell in the UAM burnup benchmark were quantified based on ENDF/B-Ⅶ.1. The numerical results show that the uncertainty of infinite neutron multiplication factors induced by decay half-lives and fission yields is low, while the uncertainty of concentration of some fission product nuclide is high.     

针对裂变产额和半衰期的燃耗计算灵敏度和不确定度分析方法

基于广义微扰理论推导了裂变产额和半衰期的燃耗灵敏度系数理论模型,该模型考虑了原子核密度和中子通量的相互影响,并开发了燃耗计算中有效增殖因数和原子核密度等响应参数对核数据的灵敏度和不确定度分析程序。基于评价核数据中裂变产物独立产额的标准差数据,产生了针对压缩燃耗数据库的裂变产额协方差矩阵,以提高不确定度的计算精度。基于ENDF/B-Ⅶ.1数据库量化了UAM基准题TMI-1栅元无限增殖因数及重要裂变产物和重核的原子核密度由裂变产额和半衰期引入的不确定度。数值结果表明,对于栅元无限增殖因数,裂变产额和半衰期引入的不确定度很小;对于部分裂变产物的原子核密度,裂变产额和半衰期会引入较大的不确定度。     

CERMET-SNRE堆芯物理计算分析

核火箭发动机功率高、寿命长、比冲大,在执行深空探测和星际航行任务时具有不可替代的优势。小型化是核火箭发动机的一个重要趋势,基于此提出了一种使用钨基金属陶瓷燃料的小型核火箭发动机（CERMET-SNRE）堆芯方案,并采用蒙特卡罗程序（MCNP）进行了精确建模,计算了相关物理参数。计算分析结果表明：CERMET-SNRE堆芯能谱硬,燃耗浅,后备反应性足够,功率分布合理,控制鼓与安全棒价值足够,发射掉落事故下有效增殖因数小于0.98,堆芯方案合理,满足设计要求。     

NAPTH: Neutronics analysis code system for the fusion–fission hybrid reactor with pressure tube type blanket

The fusion–fission hybrid reactor is considered as a potential path to the early application of fusion energy. A new concept with pressure tube type blanket has recently been proposed for a feasible hybrid reactor. In this paper, a code system for the neutronics analysis of the pressure tube type hybrid reactor is developed based on the two-step calculation scheme: the few-group homogeneous constant calculation and the full blanket calculation. The few-group homogeneous constants are calculated using the lattice code DRAGON4. The blanket transport calculation is performed by the multigroup Monte Carlo code. A link procedure for fitting the cross sections is developed between these two steps. An additional procedure is developed to calculate the burnup, power distribution, energy multiplication factor, tritium breeding ratio and neutron multiplication factor. From some numerical results, it is found that the code system NAPTH is reliable and exhibits good calculation efficiency. It can be used for the conceptual design of the pressure tube type hybrid reactor with precise geometry.     

Calculation of fuel burnup and radionuclide inventory in the Syrian miniature neutron source reactor using the GETERA code

Calculations of the fuel burnup and radionuclide inventory in the Syrian miniature neutron source reactor (MNSR) after 10 years (the reactor core expected life) of the reactor operation time are presented in this paper using the GETERA code. The code is used to calculate the fuel group constants and the infinite multiplication factor versus the reactor operating time for 10, 20, and 30 kW operating power levels. The amounts of uranium burntup and plutonium produced in the reactor core, the concentrations and radionuclides of the most important fission products and actinide radionuclides accumulated in the reactor core, and the total radioactivity of the reactor core were calculated using the GETERA code as well. It is found that the GETERA code is better than the WIMSD4 code for the fuel burnup calculation in the MNSR reactor since it is newer, has a bigger library of isotopes, and is more accurate.     

Criticality of the thorium burnup in equilibrium state

The effective neutron multiplication factor (k_eff) as a function of burnup for different volume coolant (CoR) and fuel (FR) to cell ratio is presented. Additionally the Conversion Ratio (CR) of Th-232 to U-233, concentration of U-233, fissile and fission products calculation as a function of burnup are presented. The assembly is a critical reactor which makes volumes of coolant and fuel changes possible. In addition, an analytical model of calculation of k_eff as a function of U-233 and a poison concentration in equilibrium state are presented. One can achieve the criticality of Thorium Breeder Reactor (TBR) for enough high average neutron energy which one can obtain in Fast Breeder Reactor (FBR) only. The maximal value of CR and burnup for case of k_eff ≥ 1 achieves 1.4 and 360 GWd/MTU, correspondently. The calculations were done with a MCNPX 2.7 code using F₂Be, Na and Pb coolants.     

钍基柱状高温气冷堆不同启动燃料特性初步分析

钍是一种可转换材料,将其转换成²³³U能极大提高现有核燃料资源的储量。为实现对钍的合理利用,以模块式柱状高温气冷堆GT-MHR的燃料组件作为研究对象,选取低浓缩铀、武器级钚、核反应堆级钚等作为其启动燃料。利用栅格输运计算程序DRAGON对这3种启动燃料下的钍基柱状燃料组件的寿期初中子能谱、无限增殖系数、燃耗、转换比以及²³³U和²³²Th的含量等参数进行了分析。结果表明,在易裂变物质初装量约为9%时,与低浓缩铀和武器级钚相比,核反应堆级钚作为启动燃料时组件寿期初中子能谱较硬、转换比较高;其燃耗达90 GW•d/tHM;其无限增殖系数在寿期内的波动最小;燃耗为75 GW•d/tHM时组件中²³³U存余量与²³²Th消耗量之比达0.566。     

Effect of Moderator Density Distribution of Annular Flow on Fuel Assembly Neutronic Characteristics in Boiling Water Reactor Cores

The effect of the moderator density distribution of annular flow on the fuel assembly neutronic characteristics in a boiling water nuclear reactor was investigated using the SRAC95 code system. For the investigation, a model of annular flow for fuel assembly calculation was utilized. The results of the assembly calculation with the model (Method 1) and those of the fuel assembly calculation with the uniform void fraction distribution (Method 2) were compared. It was found that Method 2 underestimates the infinite multiplication factor in the fuel assembly including the gadolinia rod (type 1 assembly). This phenomenon is explained by the fact that the capture rate in the thermal energy region in gadolinia fuel is estimated to be smaller when the liquid film of annular flow at the fuel rod surface is considered. A burnup calculation was performed under the condition of a void fraction of 65% and a volumetric fraction of the liquid film in liquid phase of 1. It is found that Method 2 underestimates the infinite multiplication factor in comparison to Method 1 in the early stage of burnup, and that Method 2 becomes to overestimate the factor after a certain degree of burnup. This is because Method 2 overestimates the depletion rate of the gadolinia.