船用反应堆屏蔽设计的可视化与快速计算功能开发

船用反应堆的屏蔽设计问题直接关系到核能能否安全的用作舰船的动力系统。MCNP在船用反应堆的屏蔽计算中应用十分广泛,但其输入程序的编写及输出结果的整理较为繁琐,为了使用户更加简便的编写MCNP输入文件,直观的分析输出结果,本文开发了针对MCNP输入与输出文件的可视化软件。此外,在船用反应堆的屏蔽设计过程中需要MCNP进行大量屏蔽计算,所耗时间过长,为了实现在一定误差范围内的快速计算功能,本文采用BP神经网络模拟学习MCNP的计算过程,仅需给出指定的输入变量即可预测屏蔽计算输出结果,解决了MCNP计算耗时过长问题,提高了屏蔽设计优化效率。     

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

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

MCMG蒙特卡罗多群-连续截面耦合中子输运计算

本文针对多群蒙特卡罗计算省时但共振自屏处理存在缺陷,以及连续截面蒙特卡罗输运计算精度高但计算费时的问题,发展了一种多群-连续截面耦合计算方法。该方法在自主研发的三维中子-光子耦合输运蒙特卡罗程序MCMG中得到应用,通过多个模型的计算验证了方法的有效性。MCMG耦合计算取得了与连续点截面MCNP程序一致的结果,其计算速度较MCNP的提高了1倍左右。     

三维多群P5中子输运蒙特卡罗程序MCMG及检验

本工作介绍了自主开发研制的三维多群P₅中子输运蒙特卡罗程序MCMG及从ENDF/B-Ⅶ库制作的47群P₅中子截面库G47B7P5N。MCMG程序发展了针对物质的碰撞机制,适合ANISN格式和非标准ANISN格式的多群中子截面。程序计算了6个临界基准模型和2个外源问题,模拟取得了与实验和连续截面MCNP程序一致的结果,计算速度较MCNP程序提高3倍以上。     

基于BP神经网络的船用反应堆屏蔽设计快速计算功能研究

船用反应堆的屏蔽设计问题直接关系到核能能否安全的用作舰船的动力系统,而在屏蔽设计问题中屏蔽计算是十分重要的环节。为了实现在一定误差范围内实现快速计算功能,采用BP神经网络模拟学习MCNP的计算过程,仅需给出指定的输入变量即可预测屏蔽计算输出结果,解决了MCNP计算耗时过长问题,提高了屏蔽设计优化效率。     

压水堆主冷却剂16N的比活度计算

对压水堆主冷却剂中16N的快中子活化生成和衰变平衡机理进行了分析,建立了基于堆芯流道求和的16N比活度计算模型,此模型以堆芯各流道内流速及平均中子注量率为输入数据,计算反应堆运行时主冷却剂中16N的稳定浓度。以秦山二期核电厂为研究对象,在对堆芯进行精确MCNP建模的基础上,对堆芯中子注量率分布进行了MCNP模拟计算。并将模拟计算数据代入16N比活度计算模型,对主冷却剂中16N比活度进行了综合计算,计算结果与工程实用参考值吻合。     

基于蒙特卡罗-离散纵标双向耦合方法的快中子注量基准分析

蒙特卡罗（MC）-离散纵标（S_N）双向耦合方法是解决大型复杂核装置屏蔽问题的有效方法。本文针对三维MC-S_N双向耦合方法在大型压水堆核电站屏蔽计算中的应用,进行了基准验证分析。基于美国核管会（NRC）发布的NUREG/CR-6115压水堆基准模型,采用自主开发的三维MC-S_N双向耦合屏蔽计算分析方法,利用MCNP4C精确计算堆芯到热屏蔽精细模型以及位于压力容器内部计算区域的精确模型,三维S N 程序TORT用于进行热屏蔽到第2下降区外表面间的计算。通过自主研发的接口程序实现MC粒子概率分布与S_N角通量密度间的相互转换,实现MC和S_N 双向耦合计算。三维MC-S_N双向耦合方法计算结果与基准报告提供的MCNP、DORT结果符合良好,初步验证了该方法解决大型复杂核装置屏蔽问题的可行性。     

MCAM在ITER装置TRIPOLI三维中子学建模中的应用

目前,国际热核聚变实验反应堆ITER装置仅有针对MCNP程序的三维中子学基准模型(ITERA-lite4),因此无法使用TRIPOLI程序对ITER装置进行中子学计算分析.本文利用蒙特卡罗计算自动建模软件MCAM 5.1创建ITER装置的三维中子学TRIPOLI模型,并通过TRIPOLI程序对其进行中子学计算分析.计算...     

EPR堆芯严重事故下安全壳内的辐射屏蔽计算

MCNP程序可以从粒子输运、扩散方程的角度来模拟计算堆芯在严重事故下安全壳内的辐射剂量水平。文章以EPR堆芯为例,采用MCNP 5程序及其核数据库CCC-710建立了精确的三维蒙特卡罗模型,在此基础上对EPR严重事故下安全壳内的辐射剂量率进行了计算分析,为判断堆芯情况和制定应急防护行动提供了数据参考。     

三维多群中子输运蒙特卡罗程序MCMG-Ⅱ基准检验

三维多群P₃中子输运蒙特卡罗程序MCMG通过版本更新和功能扩充，能够配备各种多群微观、宏观截面库，截面输入文件进一步简化，版本从Ⅰ发展到Ⅱ，参数更新到ENDF/B-Ⅶ库。程序发展了针对物质的碰撞机制，具有并行计算功能。对于12个临界基准问题和1个外源问题，MCMG-Ⅱ计算取得了与实验和连续截面MCNP 5程序一致的结果，计算速度较MCNP-5提高了3~6倍。     

Research on Coupling of Neutronics and Thermal-hydraulics for Fuel Assembly of Thorium Molten Salt Reactor Moderatedby Zirconium Hydride Rod

Based on Monte Carlo particle transport code MCNP and self-developed sub-channel thermal-hydraulic code SubTH, a code system MCNP-SubTH coupling neutronics with thermal-hydraulics was developed, which was suitable for steady state analysis for a thorium molten salt reactor moderated by zirconium hydride rod (ZrH-MSR). It solved the difficulties in the neutronics and thermal-hydraulics coupling code due to different mesh types, and has a general validity. MCNP-SubTH exchanged data between MCNP and SubTH by an external coupling. The power density field obtained from MCNP was provided as a SubTH solution file to give a user-specified source term, and then the density and temperature field from SubTH was updated and as a new MCNP input file by MCNP-SubTH to realize iterative calculation. The accuracy of MCNP-SubTH was verified by each relatively independent module. MCNP-SubTH application in the fuel assembly of ZrH-MSR was studied, and its validity was verified.     

棒状氢化锆慢化钍基熔盐堆燃料组件稳态核热耦合程序开发

基于蒙特卡罗粒子输运程序MCNP与自主开发的子通道热工水力学程序SubTH，开发了棒状氢化锆慢化钍基熔盐堆燃料组件稳态核热耦合程序MCNP-SubTH，解决核热耦合程序因网格类型不同难以耦合的问题，程序具有普适性。MCNP-SubTH通过外耦合的方式进行MCNP和SubTH之间的数据交换，将MCNP计算得到的功率场加载到SubTH的求解文件中，然后将SubTH计算得到的密度和温度场更新到MCNP的输入卡中，实现程序迭代计算。分模块验证了MCNP-SubTH的准确性，并用MCNP-SubTH对棒状氢化锆慢化钍基熔盐堆燃料组件进行了稳态核热耦合计算，验证了核热耦合方法的有效性。     

The comparison of MCNP perturbation technique with MCNP difference method in critical calculation

For a nuclear fission system, we calculated Δk_eff, which arise from system material composition changes, by two different approaches, the MCNP perturbation technique and the MCNP difference method. For every material composition change, we made four different runs, each run with different cycles or each cycle generating different neutrons, then we compared the two Δk_eff that are obtained by two different approaches. As a material composition change in any particular cell of the nuclear fission system is small compared to the material compositions in the whole nuclear fission system, in other words, this composition change can be treated as a small perturbation, the Δk_eff results obtained from the MCNP perturbation technique are much quicker, much more efficient and reliable than the results from the MCNP difference method.When a material composition change in any particular cell of the nuclear fission system is significant compared to the material compositions in the whole nuclear fission system, both the MCNP perturbation technique and the MCNP difference method can give satisfactory results. But for the run with the same cycles and each cycle generating the same neutrons, the results obtained from the MCNP perturbation technique are systemically less than the results obtained from the MCNP difference method. To further confirm our calculation results from the MCNP4C, we run the exact same MCNP4C input file in MCNP5, the calculation results from MCNP5 are the same as the calculation results from MCNP4C.We need caution when using the MCNP perturbation technique to calculate the Δk_eff as the material composition change is large compared to the material compositions in the whole nuclear fission system, even though the material composition changes of any particular cell of the fission system still meet the criteria of MCNP perturbation technique.     

中国评价核衰变数据库管理程序系统

简要介绍了中国评价核衰变数据库的组织结构及其库的管理程序系统与具体应用。     

Performance of the improved version of the Monte Carlo A3MCNP code in the design of cask shielding

Abstract

A³MCNP (automatic adjoint accelerated MCNP) is a revised version of the MCNP Monte Carlo code that automatically prepares variance reduction parameters for the CADIS (consistent adjoint driven importance sampling) methodology. Using a deterministic 'importance' (or adjoint) function, CADIS performs source and radiation particle transport biasing within the weight-window technique. The current version of A³MCNP uses the 3-D Sn transport TORT code to determine a 3-D importance function distribution. Based on simulation of several real-life problems, it is demonstrated that A³MCNP provides precise calculation results with a remarkably short computation time by using the proper and objective variance reduction parameters. However, since the first version of A³MCNP provided only a point source configuration option for large-scale shielding problems, such as spent fuel transport casks, a large amount of memory may be necessary to store enough points to properly represent the source. Hence, we have developed an improved version of A³MCNP (referred to as A³MCNPV) which has a volumetric source configuration option. This paper describes the successful use of A³MCNPV in solving the cask neutron and gamma-ray shielding problem.     

快堆用蒙特卡罗多群数据库制作方法研究

基于离散角方法，开发了蒙特卡罗多群数据库生成程序MGXSMC，该程序可以实现从输入文件读取截面数据或者从指定格式的截面库中读取截面，产生可供蒙特卡罗程序MCNP或RMC计算的数据库，并且可自动生成相应的索引文件列表。采用二维两群不带反射层的国际原子能机构（IAEA）压水堆（PWR）基准题和铅基快堆（RBEC-M）基准题对MGXSMC程序加工产生的核数据进行验证，计算结果表明，采用P5阶近似多群截面与连续点截面计算的有效增殖系数（keff）结果相差24 pcm（1pcm=10-5），而采用P0阶近似多群截面与连续点截面计算的keff结果相差较大。由此说明蒙特卡罗多群数据库的制作方法和所开发的程序是正确的，同时，中子各向异性散射对铅基快堆计算结果影响较大，故制作蒙特卡罗多群数据库时应加入中子散射角数据。      

快堆用蒙特卡罗多群数据库制作方法研究

Based on the discrete angle method, a Monte Carlo multi-group cross section generation program MGXSMC was developed. This program can read the cross section data from an input file or read the cross section from a library in a specified format to generate the multi-group cross section for MCNP or RMC. The corresponding index file list can be automatically generated. The two-dimensional two-group IAEA pressurized water reactor (PWR) benchmark and lead-based fast reactor (RBEC-M) benchmark were used to verify the cross section library generated by the MGXSMC program. The calculation results show that the difference between the calculated result of the P5 order approximate multigroup section and the continuous point cross section is 24 pcm (1pcm = 10-5), and the difference of the keff result calculated by the P0 order approximate multigroup section and the continuous point section is large. This shows that the method and the program developed for the Monte Carlo Group Section Library are correct. At the same time, the neutron anisotropic scattering has a large impact on the calculation results of the lead-based fast reactor. Therefore, when the Monte Carlo Group Section library is produced, the neutron scattering angle data should be added.