MC21 v.6.0 – A continuous-energy Monte Carlo particle transport code with integrated reactor feedback capabilities

Abstract MC21 is a continuous-energy Monte Carlo radiation transport code for the calculation of the steady-state spatial distributions of reaction rates in three-dimensional models. The code supports neutron and photon transport in fixed source problems, as well as iterated-fission-source (eigenvalue) neutron transport problems. MC21 has been designed and optimized to support large-scale problems in reactor physics, shielding, and criticality analysis applications. The code also supports many in-line reactor feedback effects, including depletion, thermal feedback, xenon feedback, eigenvalue search, and neutron and photon heating. MC21 uses continuous-energy neutron/nucleus interaction physics over the range from 10−5 eV to 20 MeV. The code treats all common neutron scattering mechanisms, including fast-range elastic and non-elastic scattering, and thermal- and epithermal-range scattering from molecules and crystalline materials. For photon transport, MC21 uses continuous-energy interaction physics over the energy range from 1 keV to 100 GeV. The code treats all common photon interaction mechanisms, including Compton scattering, pair production, and photoelectric interactions. All of the nuclear data required by MC21 is provided by the NDEX system of codes, which extracts and processes data from EPDL-, ENDF-, and ACE-formatted source files. For geometry representation, MC21 employs a flexible constructive solid geometry system that allows users to create spatial cells from first- and second-order surfaces. The system also allows models to be built up as hierarchical collections of previously defined spatial cells, with interior detail provided by grids and template overlays. Results are collected by a generalized tally capability which allows users to edit integral flux and reaction rate information. Results can be collected over the entire problem or within specific regions of interest through the use of phase filters that control which particles are allowed to score each tally. The tally system has been optimized to maintain a high level of efficiency, even as the number of edit regions becomes very large.

[1]  David P. Griesheimer,et al.  Reactivity Effects of Spatial Homogenization of Thermal Feedback Regions in Monte Carlo Reactor Calculations , 2012 .

[2]  Paul K. Romano,et al.  On-the-Fly Monte Carlo Dominance Ratio Calculation Using the Noise Propagation Matrix (Selected Papers of the Joint International Conference of Supercomputing in Nuclear Applications and Monte Carlo : SNA + MC 2010) , 2011 .

[3]  Brian R. Nease,et al.  Time Series Analysis of Monte Carlo Fission Sources—III: Coarse-Mesh Projection , 2007 .

[4]  D. P. Griesheimer,et al.  A generalized framework for in-line energy deposition during steady-state Monte Carlo radiation transport , 2013 .

[5]  E. D. Blakeman,et al.  Advanced variance reduction strategies for optimizing mesh tallies in MAVRIC , 2007 .

[6]  Benoit Forget,et al.  Benchmark for evaluation and validation of reactor simulations (BEAVRS) , 2013 .

[7]  Philip F. Rose,et al.  ENDF-6 Formats Manual , 1997 .

[8]  T. M. Sutton,et al.  A method for reducing the largest relative errors in Monte Carlo iterated-fission-source calculations , 2013 .

[9]  David L. Millman,et al.  Geometric Templates for Improved Tracking Performance in Monte Carlo Codes , 2014, ICS 2014.

[10]  R. A. Forster,et al.  Estimation and Interpretation of keff Confidence Intervals in MCNP , 1995 .

[11]  R. J. Brissenden,et al.  Biases in the estimation of Keff and its error by Monte Carlo methods , 1986 .

[12]  Martin J. Berger,et al.  Bremsstrahlung spectra from electron interactions with screened atomic nuclei and orbital electrons , 1985 .

[13]  J. H. Hubbell,et al.  EPDL97: the evaluated photo data library `97 version , 1997 .

[14]  David C. Carpenter A COMPARISON OF CONSTANT POWER DEPLETION ALGORITHMS , 2009 .

[15]  T. M. Sutton,et al.  Comparison of Some Monte Carlo Models for Bound Hydrogen Scattering , 2009 .

[16]  Jerome Spanier Two Pairs of Families of Estimators for Transport Problems , 1966 .

[17]  G. D. Byrne,et al.  VODE: a variable-coefficient ODE solver , 1989 .

[18]  W. Heitler,et al.  The quantum theory of radiation , 1936 .

[19]  D. J. Kelly,et al.  MC21 analysis of the MIT PWR benchmark: Hot zero power results , 2013 .

[20]  W. C. Dickinson,et al.  CALCULATION OF FORWARD BREMSSTRAHLUNG SPECTRA FROM THICK TARGETS. , 1968 .

[21]  Alireza Haghighat,et al.  Automated variance reduction of Monte Carlo shielding calculations using the discrete ordinates adjoint function , 1998 .

[22]  L. Levitt,et al.  The Probability Table Method for Treating Unresolved Neutron Resonances in Monte Carlo Calculations , 1972 .

[23]  T Donovan,et al.  Geometric Representations in the Developmental Monte Carlo Transport Code MC21 , 2006 .

[24]  D. C. Carpenter,et al.  The MC21 Monte Carlo Transport Code , 2007 .

[25]  Stephen M. Bowman,et al.  Validation of the scale system for PWR spent fuel isotopic composition analyses , 1995 .

[26]  Ely M. Gelbard,et al.  Computation of standard deviations in Eigenvalue calculations , 1990 .

[27]  E. D. Blakeman,et al.  Forward-weighted CADIS method for global variance reduction , 2007 .

[28]  M. J. Berger,et al.  Density effect for the ionization loss of charged particles in various substances , 1984 .

[29]  Forrest B. Brown,et al.  Stationarity Modeling and Informatics-Based Diagnostics in Monte Carlo Criticality Calculations , 2005 .

[30]  W. W. Scott Angular distribution of thick-target bremsstrahlung which includes multiple electron scatterings , 1967 .

[31]  D. J. Kelly,et al.  MC21 analysis of the nuclear energy agency Monte Carlo performance benchmark problem , 2012 .

[32]  Bojan. Petrovic,et al.  MONTE CARLO PERFORMANCE BENCHMARK FOR DETAILED POWER DENSITY CALCULATION IN A FULL SIZE REACTOR CORE Benchmark specifications Revision 1 , 2010 .

[33]  R.E. MacFarlane,et al.  The NJOY Nu-clear Data Processing System Version 91 , 1994 .