基于高阶CRAM的燃耗程序开发与验证

本文基于高阶切比雪夫有理近似方法（CRAM）研制了点燃耗程序ICRAM,并内耦合于蒙特卡罗输运程序OpenMC,形成了一套燃耗计算分析程序OPICE。与传统部分分式分解（PFD）形式的CRAM相比,高阶不完全局部分解（IPF）形式的CRAM具有数值稳定性好、计算精度高和步长包容性更好等特点,满足高保真燃耗计算发展的需求。为提高耦合计算精度,OPICE采用了预估-校正和子步法两种耦合策略,支持纯衰变、定通量和定功率3种计算模式。通过OECD/NEA压水堆栅元燃耗基准题和快堆燃耗基准题的验证,程序计算结果与实验值及各参考值吻合良好,初步验证了OPICE的正确性与有效性。     

基于切比雪夫有理逼近方法的蒙特卡罗燃耗计算研究与验证

燃耗计算在反应堆设计、分析研究中起着重要作用。相比于传统点燃耗算法,切比雪夫有理逼近方法(Chebyshev rational approximation method,CRAM)具有计算速度快、精度高的优点。基于超级蒙特卡罗核计算仿真软件系统Super MC(Super Monte Carlo Simulation Program for Nuclear and Radiation Process),采用切比雪夫有理逼近方法和桶排序能量查找方法,进行了蒙特卡罗燃耗计算的初步研究与验证。通过燃料棒燃耗例题以及IAEA-ADS(International Atomic Energy Agency-Accelerator Driven Systems)国际基准题,初步验证了该燃耗计算方法的正确性,且IAEA-ADS基准题测试表明,与统一能量网格方法相比,桶排序能量查找方法在保证了计算效率的同时减少了内存开销。     

基于最佳一致逼近多项式的燃耗计算方法研究

本文研究了一种基于最佳一致逼近多项式（MMPA）的燃耗计算方法求解燃耗方程。相比于切比雪夫有理近似方法（CRAM）和围道积分有理近似方法（QRAM）,MMPA方法只需一次矩阵求逆计算即可求解燃耗方程,且所有计算都是实数运算,具有数值稳定性好、求解效率高等优点。进一步研制了基于MMPA方法的点燃耗程序AMAC,并耦合蒙特卡罗输运程序OpenMC,采用衰变例题、固定辐照例题、OECD/NEA压水堆栅元燃耗基准题和沸水堆组件燃耗基准题进行验证,程序计算结果与实验值及各参考值吻合良好,初步验证了MMPA方法在理论和数值上的正确性和有效性。     

高阶预估-校正燃耗方法在压水堆组件计算中的应用

压水堆燃料组件输运燃耗耦合计算通常采用的是传统的预估-校正(PC)燃耗方法。然而,该方法本身的假设导致其存在一定的计算误差。为进一步提高燃耗计算的精度,本文针对传统的预估-校正燃耗方法的缺陷研究了改进的预估-校正燃耗方法,改进了对核反应率进行修正的高阶预估-校正燃耗方法,并在Bamboo-Lattice程序中进行了程序实现,对该方法进行了验证分析。结果表明:改进的预估-校正燃耗方法和高阶预估-校正燃耗方法在保证计算效率的前提下提高了燃耗计算的精度。     

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

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

Development and validation of the code COUPLE3.0 for the coupled analysis of neutron transport and burnup in ADS

The analysis of the fuel depletion behavior is critical for maintaining the safety of accelerator-driven subcritical systems(ADSs). The code COUPLE2.0 coupling 3-D neutron transport and point burnup calculation was developed by Tsinghua University. A Monte Carlo method is used for the neutron transport analysis, and the burnup calculation is based on a deterministic method. The code can be used for the analysis of targets coupled with a reactor in ADSs. In response to additional ADS analysis requirements at the Institute of Modern Physics at the Chinese Academy of Sciences, the COUPLE3.0 version was developed to include the new functions of(1) a module for the calculation of proton irradiation for the analysis of cumulative behavior using the residual radionuclide operating history,(2) a fixed-flux radiation module for hazard assessment and analysis of the burnable poison, and(3) a module for multi-kernel parallel calculation, which improves the radionuclide replacement for the burnup analysis to balance the precision level and computational efficiency of the program. This paper introduces thevalidation of the COUPLE3.0 code using a fast reactor benchmark and ADS benchmark calculations. Moreover,the proton irradiation module was verified by a comparison with the analytic method of calculating the210 Po accumulation results. The results demonstrate that COUPLE3.0 is suitable for the analysis of neutron transport and the burnup of nuclides for ADSs.     

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

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

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

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

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

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

CRAM在放射性核素存量计算中的应用

在反应堆中,组成材料的稳定核素经受强中子辐照后,会被活化成放射性核素。这些核素及其衰变产物对工作人员的职业辐照剂量具有重要贡献。为了更好地进行人员的辐射防护工作,需要对放射性核素的存量进行精确计算。相对于核素平衡方程的其它求解方法,切比雪夫有理逼近方法(Chebyshev Rational Approximation Method,CRAM)在计算精度和效率方面具有综合性优势。首先介绍了CRAM的基本理论,随后选取典型的例题进行了测试验证。与解析解对比的结果表明,采用CRAM进行中子辐照下的核素活化衰变计算能够取得不错的效果,但是用于核素长期衰变计算可能导致计算错误。针对此问题,将收缩乘方技术与CRAM相结合,取得了正确的计算结果,拓展了CRAM的适用范围。     

聚变-裂变混合能源堆球模型中子学对算研究

利用蒙特卡罗程序和自主开发的蒙特卡罗-燃耗耦合程序MOCouple-s,对北京应用物理与计算数学研究所提出的聚变-裂变混合能源堆球模型进行了对算研究。对初始时刻及各燃耗时刻下的有效增殖因数、能量倍增因子、氚增殖比、中子源强度等堆芯参数进行了比较,结果总体符合较好。对寿期末重要核素的成分进行了详细比较,除个别核素外,偏差很小,表明所采用的计算程序与核参数库一致性良好。对核参数库的选择、铀水体积比等对燃耗计算结果的影响进行敏感性分析,并对外中子源驱动的次临界堆芯的燃耗计算进行详细讨论,提出可行的燃耗计算基准。     

Comparison of depletion algorithms for large systems of nuclides

In this work five algorithms for solving the system of decay and transmutation equations with constant reaction rates encountered in burnup calculations were compared. These are Chebyshev rational approximation method (CRAM), which is a new matrix exponential method, the matrix exponential power series with instant decay and a secular equilibrium approximations for short-lived nuclides, which is used in ORIGEN, and three different variants of transmutation trajectory analysis (TTA), which is also known as the linear chains method. The common feature of these methods is their ability to deal with thousands of nuclides and reactions. Consequently, there is no need to simplify the system of equations and all nuclides can be accounted for explicitly.The methods were compared in single depletion steps using decay and cross-section data taken from the default ORIGEN libraries. Very accurate reference solutions were obtained from a high precision TTA algorithm. The results from CRAM and TTA were found to be very accurate. While ORIGEN was not as accurate, it should still be sufficient for most purposes. All TTA variants are much slower than the other two, which are so fast that their running time should be negligible in most, if not all, applications. The combination of speed and accuracy makes CRAM the clear winner of the comparison.     

HARMONY程序计算中子扩散方程高阶λ本征值问题的基准验证

高阶λ谐波在反应堆堆芯功率重构、换料优化、ADS次临界反应堆物理特性研究等领域有着重要应用价值。为进行高阶λ谐波的计算,本文基于隐式重启动Arnoldi方法(IRAM)编制了可用于一维、二维、三维笛卡尔坐标系中子扩散方程的任意阶λ谐波及本征值计算的HARMONY程序,并进行了基准题的数值验证。结果表明,HARMONY程序能实现高阶λ本征值问题计算,具有较高的精度,为未来基于λ谐波的ADS次临界反应堆物理特性研究奠定了基础。     

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.     

Calculation of radial burnup and nuclides atom density distributions in a VVER-1000 fuel rod using Monte Carlo method

A fundamental knowledge of fuel behavior in different situations is required for safe and economic nuclear power generation. Due to the importance of a fuel rod behavior modelling in high burnup, in this paper, the radial distribution of burnup, fission products, and actinides atom density and their variations by increasing burnup and other factors such as temperature, enrichment and power density are studied in a fuel pellet of a VVER-1000 reactor in an operational cycle using the MCNPX 2.7 Monte Carlo code. A benchmark including a Uranium-Gadolinium (UGD) fuel assembly is used for verification of the developed model in the MCNPX code for radial burnup calculation. A sensitivity study is carried out to investigate the effect of different parameters such as the number of particles per cycle, the number of geometrical radial nodes in the fuel pellet, the number of burnup steps and the selection of different fission-product contents (i.e. those isotopes that are used for particle transport) on the MCNPX model for speed and accuracy compromising. To calculate the radial temperature profiles and to analyze the effect of temperature on the radial burnup distribution and vice versa, the HEATING 7.2 code, which is a general-purpose conduction heat transfer program, and the MCNPX code are applied together. The results show the accuracy and capability of the proposed model in the MCNPX and HEATING codes for radial burnup calculation.     

微观燃耗方程的求解算法比较与性能分析

燃耗方程的求解是燃耗计算的核心。常见的算法包括泰勒方法、Pade方法、子空间方法、切比雪夫有理近似方法和龙格库塔法等。通过数值实验,对每种算法在精度、效率、稳定性方面进行分析比较。结果表明：子空间方法、泰勒方法在计算效率方面具有优势;计算精度及稳定性方面,泰勒方法和Pade方法均占优势。综合考虑,泰勒方法在3个方面均表现突出,可作为燃耗计算的优选算法。     

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

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

Development of core fuel management code system for WWER-type reactors

1 Introduction Over the past decades, although many in-core fuel management code systems for PWRs with square fuel assemblies have been developed, there are only a few codes for the cores with hexagonal assemblies (such as Russian pressurized water type WWER reac- tors). The Tianwan Nuclear Power Station in Jiangsu Province, China, is imported from Russia, which adopts the WWER-1000 reactor, and will be put into operation; therefore, the research of core fuel man- agement for WWER-typ…     

基于蒙特卡罗方法的燃耗计算误差研究

基于蒙特卡罗方法进行燃耗计算时,随着燃耗加深,燃耗的计算误差逐渐增大。本文针对蒙特卡罗方法的燃耗计算误差进行研究,并采取修正措施改善燃耗计算的精度。结果表明：采用无偏差最小方差(UMV)修正可改善统计误差的传递效应,采用密度修正可保证蒙特卡罗输运计算的准确性,在此基础上局部优化燃耗截面库,进一步改善了燃耗计算的精度,为其工程应用奠定了基础。     

三维六角形组件压水堆堆芯燃料管理计算及程序系统研究

介绍所研制的WWER型压水堆堆芯燃料管理计算程序系统TPFAP-H／CSIM-H，六角形组件均匀化计算程序TPFAP-H是在压水堆正方形组件程序TPFAP的基础上，采用穿透概率法与响应矩阵方法相结合计算六角形组件内中子能谱分布，并考虑六角形栅元特点改造开发而成的CSIM-H是以先进六角形节块扩散程序为基础．参照SIMULATE程序功能而研制的物理-热工水力耦合的三维六角形节块PWR堆芯燃料管理程序两者通过接口程序LINK连接起来，可以考虑燃耗，功率、慢化剂密度变化．控制棒、氙等参数的多种反馈效应对IAEA的WWER-1000型Kalinin核电厂基准问题的校算的结果表明，临界硼浓度、功率和燃耗分布等结果与国际各研究机构的结果吻合良好，偏差均在工程要求之内。