Evaluation of Shutdown Gamma-ray Dose Rates around the Duct Penetration by Three-Dimensional Monte Carlo Decay Gamma-ray Transport Calculation with Variance Reduction Method

For the evaluation of gamma-ray dose rates around the duct penetrations after shutdown of nuclear fusion reactor, the calculation method is proposed with an application of the Monte Carlo neutron and decay gamma-ray transport calculation. For the radioisotope production rates during operation, the Monte Carlo calculation is conducted by the modification of the nuclear data library replacing a prompt gamma-ray spectrum with a decay gamma-ray spectrum. By multiplying each correction factor, which is ratio of the actual activation level after shutdown to the production rate during operation, with each decay gamma-ray flux due to each radioisotope, the decay gamma-ray dose rate is evaluated. In order to improve the statistical error, a variance reduction method is proposed by the application of the weight window importance technique and the specification of the decay gamma-ray generation location. We identify the cell producing the decay gamma-ray which can contribute the decay gamma-ray flux in evaluation locations, and forcibly terminate the gamma-ray transport calculation in the cells except for the identified cells. In order to validate the effectiveness of the method, shielding calculation for actual ITER (International Thermonuclear Experimental Reactor) configuration is performed, and small statistical errors below criteria are obtained. The effectiveness of the proposed method for ITER design analysis is demonstrated.