Advanced Burnup Calculation Code System in a Subcritical State with Continuous-Energy Monte Carlo Code for Fusion-Fission Hybrid Reactor

The fusion-fission (FF) hybrid reactor is a promising energy source that is thought to act as a bridge between the existing fission reactor and the genuine fusion reactor in the future. The burnup calculation system that aims at precise burnup calculations of a subcritical system was developed for the detailed design of the FF hybrid reactor, and the system consists of MCNP, ORIGEN, and postprocess codes. In the present study, the calculation system was substantially modified to improve the calculation accuracy and at the same time the calculation speed as well. The reaction rate estimation can be carried out accurately with the present system that uses track-length (TL) data in the continuous-energy treatment. As for the speed-up of the reaction rate calculation, a new TL data bunching scheme was developed so that only necessary TL data are used as long as the accuracy of the point-wise nuclear data is conserved. With the present system, an example analysis result for our proposed FF hybrid reactor is described, showing that the computation time could really be saved with the same accuracy as before.     

Validation of IRBURN calculation code system through burnup benchmark analysis

Assessment of the reactor fuel composition during the irradiation time, fuel management and criticality safety analysis require the utilization of a validated burnup calculation code system. In this work a newly developed burnup calculation code system, IRBURN, is introduced for the estimation and analysis of the fuel burnup in LWR reactors. IRBURN provides the full capabilities of the Monte Carlo neutron and photon transport code MCNP4C as well as the versatile code for calculating the buildup and decay of nuclides in nuclear materials, ORIGEN2.1, along with other data processing and linking subroutines. This code has the capability of using different depletion calculation schemes.     

基于蒙特卡罗方法的三维燃耗计算研究

采用通过编写连接MCNP程序和ORIGEN2程序的接口处理程序的方法进行快中子系统的燃耗计算。由MCNP、ORIGEN2、接口处理程序和截面文件组成的软件系统可用于燃料或堆芯非均匀布置快中子系统的燃料同位素成分和燃耗反应性损失计算,在燃耗反应性损失计算中采用了伪裂变产物的方法。介绍程序系统的研制情况,并给出用该软件系统计算中国实验快堆首炉堆芯和OECD/NEAMOX燃料快堆基准题的燃耗计算结果。     

基于MCNP和ORIGEN的熔盐快堆燃耗分析计算

熔盐堆是6种第4代先进核能系统中唯一使用液态燃料设计的反应堆型,其堆芯一回路中循环流动的熔盐既是燃料,也是冷却剂。这一特征在省去燃料元件加工制造步骤的同时,也使得熔盐堆能进行在线处理和在线添料的操作。因此,传统固态反应堆燃耗分析程序不再适用于熔盐堆。本文以熔盐快堆(MSFR)为分析对象,基于物理分析程序MCORE(MCNP+ORIGEN),将上述熔盐堆特点考虑进去,开发出能进行熔盐堆燃耗分析的MCORE-MS。初步分析表明,233 U-started模式下,熔盐在线处理可有效降低堆芯熔盐中裂变产物的含量,提高中子经济性。MSFR运行过程中能够一直保持负的温度反应性系数。     

ITER中国液态锂铅实验包层模块中子学分析

基于ITER装置全模型,借助于MCNP自动建模程序MCAM,将TBM模块插入该模型的赤道窗口,使用MCNP/4C和FENDL1.0数据库,对DLL和SLL两个典型子模块设计进行三维中子学计算和分析,给出TBM模块核热功率密度分布以及氚增殖能力.     

堆用蒙卡程序燃耗计算功能开发

堆用蒙卡程序(RMC)是由清华大学工程物理系REAL实验室自主开发的用于反应堆物理分析的中子输运蒙卡程序,本文主要介绍其燃耗计算功能的开发与验证。RMC的燃耗计算功能具有的特点:内部耦合ORIGEN,相比于外耦合方式,更加灵活和高效;使用基于能谱的单群截面统计方法,可在保证精度的前提下,显著提高计算效率;采取预估修正和中点近似等多种燃耗步策略,减小大燃耗步长时的计算误差。通过计算压水堆栅元、沸水堆组件、快堆等一系列基准题和算例,验证了RMC燃耗计算的正确性和速度优势。     

IAEA加速器驱动系统(ADS)中子学基准问题的计算

为了验证加速器驱动洁净核能系统研究拟采用的程序系统,对IAEA加速器驱动系统(ADS)中子学第一阶段基准问题进行了校算。其中,散裂中子源的计算采用LAHET程序;中子输运计算采用MCNP程序,核素的燃耗计算采用0RIGEN2。计算结果与IAEA的ADS研究协调项目组(ADSCRP)成员俄罗斯和瑞典的结果吻合较好。     

球床高温堆平衡态燃耗计算程序的开发

基于MCNP5和ORIGEN2耦合方法,开发了平衡态下球床高温堆的燃耗计算程序PBRE,用于堆的性能价值分析。为节省蒙特卡罗计算时间,对迭代收敛的方法进行优化,使之可在10个迭代步内收敛。使用PBRE对清华大学HTR-10进行建模计算,得到的平均卸料燃耗深度与文献报道值一致,表明PBRE程序适用于球床堆平衡态的燃耗分析。     

CERMET-SNRE堆芯物理计算分析

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

基于切比雪夫有理近似方法的点燃耗程序开发

基于切比雪夫有理近似方法（CRAM）开发了点燃耗求解程序。程序采用2套燃耗数据库，精细燃耗数据库和简化燃耗数据库，并将点燃耗程序与输运系统耦合。计算了定注量率辐照问题和衰变问题，以及JAEA轻水堆基准题，计算结果与国际知名程序对比。结果表明，程序在定注量率辐照问题和衰变问题的计算上，核子密度精度与ORIGEN2相当，单栅元和组件计算结果与HELIOS1.11以及参考解吻合良好。      

反应堆临界-燃耗耦合蒙特卡罗计算

基于连续点截面MCNP程序 ,研制了三维多群P3 中子输运蒙特卡罗程序MCMG ,并与栅元均匀化程序WIMS耦合 ,实现了临界 燃耗耦合计算。采用WIMS产生的 69群共振、自屏宏观中子截面和BUGLE 80u47群微观中子截面 ,分别计算了简单反应堆和临界实验堆问题 ,计算结果与其它输运方法的计算结果和试验结果一致。在相同计算精度下 ,MCMG的计算时间较MCNP的计算时间少     

基于CRAM的压水堆燃耗程序开发及验证

本文基于Cinder90燃耗数据库开发了燃耗求解程序MCRAM，并耦合MCNP程序对重要的锕系核素和裂变产物核素的反应截面进行了修正。以OECD/NEA乏燃料成分基准数据库中的Takahama-3压水堆燃料组件为基准题，对MCRAM程序的计算结果进行了验证，并与其他程序的计算结果进行了比较。结果表明，MCRAM程序对重要裂变产物和主要锕系核素的计算结果相对偏差小于5%，计算精度与ORIGEN2程序的相当。与此同时，同一例题的计算效率MCRAM较之MCNTRANS程序提高了近200倍。     

Time dependent burn-up and fission products inventory calculations in the discharged fuel of the Syrian MNSR

The Syrian Miniature Neutron Source Reactor (MNSR), a 30 kW, 89.8% HEU fueled (U-Al), went critical in March, 1996. By operating the reactor at nominal power for 2.5 h/day, the estimated core life is 10 years. This paper presents the results of fuel burn-up and depletion analysis of the MNSR fuel lattice using the ORIGEN 2 code. A one-group cross-section data base for the ORIGEN 2 computer code was developed for the Syrian MNSR research reactor. The ORIGEN 2 predicted burn-up dependent actinide compositions of MNSR spent fuel using the newly developed data base show a good agreement with the published results in the literature. In addition, the burn-up characteristics of MNSR spent fuel was analyzed with the new data base. Finally, to study the effect of burn-up on the reactivity, the microscopic cross-sections of the fission products calculated by the WlMS code (using the number densities of fission products generated by the ORIGEN 2 code as a function of burn-up time), were used as an input for the CITATION code calculations. The results contained in this paper could be used in performing criticality safety analysis and shielding calculations for the design of a spent fuel storage cask for the MNSR core.     

Validation of a new continuous Monte Carlo burnup code using a Mox fuel assembly

The reactivity of nuclear fuel decreases with irradiation (or burnup) due to the transformation of heavy nuclides and the formation of fission products. Burnup credit studies aim at accounting for fuel irradiation in criticality studies of the nuclear fuel cycle (transport, storage, etc…). The principal objective of this study is to evaluate the potential capabilities of a newly developed burnup code called “BUCAL1”. BUCAL1 differs in comparison with other burnup codes as it does not use the calculated neutron flux as input to other computer codes to generate the nuclide inventory for the next time step. Instead, BUCAL1 directly uses the neutron reaction tally information generated by MCNP for each nuclide of interest to determine the new nuclides inventory. This allows the full capabilities of MCNP to be incorporated into the calculation and a more accurate and robust analysis to be performed. Validation of BUCAL1 was processed by code-to-code comparisons using predictions of several codes from the NEA/OCED. Infinite multiplication factors (k_∞) and important fission product and actinide concentrations were compared for a MOX core benchmark exercise. Results of calculations are analysed and discussed.     

Development and validation of depletion code system IMPC Burnup for ADS

Depletion calculation is important for studying the transmutation efficiency of minor actinides and longlife fission products in accelerator-driven subcritical reactor system(ADS). Herein the Python language is used to develop a burnup code system called IMPC-Burnup by coupling FLUKA, OpenMC, and ORIGEN2. The program is preliminarily verified by OECD-NEA pin cell and IAEAADS benchmarking by comparison with experimental values and calculated results from other studies. Moreover,the physics design scheme of the CIADS subcritical core is utilized to test the feasibility of IMPC-Burnup program in the burnup calculation of ADS system. Reference results are given by the COUPLE3.0 program. The results of IMPC-Burnup show good agreement with those of COUPLE3.0. In addition, since the upper limit of the neutron transport energy for OpenMC is 20 MeV, neutrons with energies greater than 20 MeV in the CIADS subcritical core cannot be transported; thus, an equivalent flux method has been proposed to consider neutrons above 20 MeV in the OpenMC transport calculation. The results are compared to those that do not include neutrons greater than 20 MeV. The conclusion is that the accuracy of the actinide nuclide mass in the burnup calculation is improved when the equivalent flux method is used. Therefore, the IMPC-Burnup code is suitable for burnup analysis of the ADS system.     

Verification of Source Term Estimation Method against Measured Data for Decommissioning Waste from a CANDU Reactor

The method for the establishment of an equilibrium core model proposed in the previous paper and the source term calculation method proposed in this paper for the characterization of decommissioning waste were verified by comparing the nuclide inventory estimated by MCNP/ORIGEN2 simulations with the measured nuclide inventory according to a chemical assay in an irradiated pressure tube discharged from Wolsong Unit 1 in 1994. At first, the time-average pseudoequilibrium full-core model of Wolsong Unit 1 was developed on the basis of the previously proposed modeling method for the activation of in-core and ex-core structural components. Then, the application level of the neutron flux and cross section in the radionuclide buildup calculation were compromised. Fourteen major actinides and fission products were considered to represent the irradiated fuel condition, and a geometry simplification was also introduced in the burned full-core model for MCNP simulation. The assumption of a constant neutron flux and capture cross section as a function of the irradiation time was applied in the radionuclide buildup calculation in ORIGEN2. As a result, the values estimated from the analysis system agreed with the measured data within a difference range of 30%. Therefore, it was found that the MCNP/ORIGEN system and source term characterization method proposed can be viable to estimate the source terms of the decommissioning waste from a CANDU reactor.     

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栅元无限增殖因数及重要裂变产物和重核的原子核密度由裂变产额和半衰期引入的不确定度。数值结果表明，对于栅元无限增殖因数，裂变产额和半衰期引入的不确定度很小；对于部分裂变产物的原子核密度，裂变产额和半衰期会引入较大的不确定度。     

SPRR-300反应堆混凝土屏蔽层内中子注量率分布研究

采用MCNP程序与ANISN程序结合的计算方案获取了SPRR-300反应堆混凝土屏蔽层内的中子注量率分布情况，同时采用固体核径迹探测器测量了混凝土屏蔽层外低水平中子注量率，两者吻合较好，说明了计算结果的可信性。上述结果为反应堆退役工作提供了放射性源项的计算依据。     

Effect of fission products yield library on the calculation of the nuclide composition of the fuel in sodium-cooled fast reactors

The calculation of the composition of irradiated fuel for different degrees of burnup is a basic problem in the analysis of nuclear-radiological safety of objects holding spent fuel assemblies. The yield of fission products is one of the important initial indicators in burnup calculations. Methods for compiling libraries of fission products yield on the basis of the ENDF/B up-to-date evaluated nuclear data files are described. The nuclide composition of uranium oxide and uranium-plutonium-zirconium metal fuel in sodium-cooled fast reactors is analyzed by means of high-precision calculations performed with different fission product yields libraries using different computer codes MONTEBURNS–MCNP5–ORIGEN2 and the results are presented.