基于三维MC-SN耦合方法的PWR压力容器快中子注量计算基准分析

Fast neutron flux calculation benchmark analysis of PWR pressure vessel based on 3D MC-SN coupled method

蒙特卡罗(MC)-离散纵标(SN)耦合方法是解决同时具有复杂几何和深穿透特点的核装置屏蔽问题的有效方法。本文首次将三维MC-SN耦合方法应用于压水堆屏蔽计算。针对NUREG/CR-6115压水堆基准模型,选取热屏蔽内表面为公共交界面,将其分为几何复杂的MC模拟区和具有深穿透特点的SN模拟区。三维MC程序用于精确描述堆芯到热屏蔽精细模型,并记录穿过热屏蔽内表面的中子径迹信息。接口程序将中子径迹转换为SN计算所需的边界源,提供给三维SN程序进行热屏蔽到压力容器的计算。计算结果包括压力容器内表面、1/4壁厚处及焊缝处快中子注量(E>1.0 MeV)圆周方向分布。三维耦合方法计算结果与基准报告提供的MCNP、DORT结果符合良好,验证了该方法处理圆柱坐标系屏蔽问题的有效性和程序使用的正确性。

The Monte Carlo （MC）-discrete ordinates （SN） coupled method is an efficient approach to solve shielding calculations of nuclear device with complex geometries and deep penetration. The 3D MC-SN coupled method has been used for PWR shielding calculation for the first time. According to characteristics of NUREG/CR-6115 PWR model, the thermal shield is specified as the common surface to link the Monte Carlo complex geometrical model and the deep penetration SN model. 3D Monte Carlo code isemployed to accurately simulate the structure from core to thermal shield. The neutron tracks crossing the thermal shield inner surface are recorded by MC code. The SN boundary source is generated by the interface program and used by the 3D SN code to treat the calculation from thermal shield to pressure vessel. The calculation results include the circular distributions of fast neutron flux at pressure vessel inner wall, pressure vessel T/4 and lower weld locations. The calculation results are performed with comparison to MCNP and DORT solutions of benchmark report and satisfactory agreements are obtained. The validity of the method and the correctness of the programs are proved.