Shutdown dose rate analysis with CAD geometry, Cartesian/tetrahedral mesh, and advanced variance reduction

Abstract In fusion energy systems (FES) high-energy neutrons born from burning plasma activate system components to form radionuclides. The biological dose rate that results from photons emitted by these radionuclides after shutdown—the shutdown dose rate (SDR)—must be quantified for maintenance planning. This can be done using the Rigorous Two-Step (R2S) method, which involves separate neutron and photon transport calculations, coupled by a nuclear inventory analysis code. The geometric complexity and highly attenuating configuration of FES motivates the use of CAD geometry and advanced variance reduction for this analysis. An R2S workflow has been created with the new capability of performing SDR analysis directly from CAD geometry with Cartesian or tetrahedral meshes and with biased photon source sampling, enabling the use of the Consistent Adjoint Driven Importance Sampling (CADIS) variance reduction technique. This workflow has been validated with the Frascati Neutron Generator (FNG)-ITER SDR benchmark using both Cartesian and tetrahedral meshes and both unbiased and biased photon source sampling. All results are within 20.4% of experimental values, which constitutes satisfactory agreement. Photon transport using CADIS is demonstrated to yield speedups as high as 8.5·105 for problems using the FNG geometry.

[1]  J. Sanz,et al.  Development of the R2SUNED Code System for Shutdown Dose Rate Calculations , 2016, IEEE Transactions on Nuclear Science.

[2]  D. Henderson,et al.  Validation of the Alara Activation Code , 1998 .

[3]  Timothy J. Tautges,et al.  MOAB : a mesh-oriented database. , 2004 .

[4]  Jonathan Cole Smith,et al.  An Analysis of the Alias Method for Discrete Random-Variate Generation , 2005, INFORMS J. Comput..

[5]  Philip R. Page,et al.  ENDF/B-VII.0: Next Generation Evaluated Nuclear Data Library for Nuclear Science and Technology , 2006 .

[6]  T. L. Edwards,et al.  CUBIT mesh generation environment. Volume 1: Users manual , 1994 .

[7]  Kevin T. Clarno,et al.  Denovo: A New Three-Dimensional Parallel Discrete Ordinates Code in SCALE , 2010 .

[8]  Ulrich Fischer,et al.  Rigorous mcnp based shutdown dose rate calculations: computational scheme, verification calculations and application to ITER , 2002 .

[9]  Arkady Serikov,et al.  Verification and validation of the R2Smesh approach for the calculation of high resolution shutdown dose rate distributions , 2012 .

[10]  Timothy J. Tautges,et al.  ACCELERATION TECHNIQUES FOR DIRECT USE OF CAD-BASED GEOMETRIES IN MONTE CARLO RADIATION TRANSPORT , 2009 .

[11]  M.E. Sawan,et al.  Application of CAD-neutronics coupling to geometrically complex fusion systems , 2009, 2009 23rd IEEE/NPSS Symposium on Fusion Engineering.

[12]  Andrew A. Davis,et al.  Comparison of global variance reduction techniques for Monte Carlo radiation transport simulations of ITER , 2011 .

[13]  Alireza Haghighat,et al.  MONTE CARLO VARIANCE REDUCTION WITH DETERMINISTIC IMPORTANCE FUNCTIONS , 2003 .

[14]  Paolo Cignoni,et al.  Generating Random Points in a Tetrahedron , 2000, J. Graphics, GPU, & Game Tools.

[15]  A. Turner,et al.  MCR2S Unstructured Mesh Capabilities for use in Shutdown Dose Rate Analysis , 2015 .

[16]  Andrew A. Davis,et al.  Benchmarking the MCR2S system for high-resolution activation dose analysis in ITER , 2010 .

[17]  Seth R. Johnson,et al.  ADVANTG An Automated Variance Reduction Parameter Generator , 2013 .

[18]  Paul P. H. Wilson,et al.  ITER Neutronics Modeling Using Hybrid Monte Carlo/Deterministic and CAD-Based Monte Carlo Methods , 2011 .

[19]  Douglas E. Peplow,et al.  FW-CADIS Method for Global and Regional Variance Reduction of Monte Carlo Radiation Transport Calculations , 2014 .