Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT

A consistent “2D/1D” neutron transport method is derived from the 3D Boltzmann transport equation, to calculate fuel-pin-resolved neutron fluxes for realistic full-core Pressurized Water Reactor (PWR) problems. The 2D/1D method employs the Method of Characteristics to discretize the radial variables and a lower order transport solution to discretize the axial variable. This paper describes the theory of the 2D/1D method and its implementation in the MPACT code, which has become the whole-core deterministic neutron transport solver for the Consortium for Advanced Simulations of Light Water Reactors (CASL) core simulator VERA-CS. Several applications have been performed on both leadership-class and industry-class computing clusters. Results are presented for whole-core solutions of the Watts Bar Nuclear Power Station Unit 1 and compared to both continuous-energy Monte Carlo results and plant data.

[1]  Han Gyu Joo,et al.  Solution of the C5G7MOX benchmark three-dimensional extension problems by the DeCART direct whole core calculation code , 2006 .

[2]  K. S. Smith.  Nodal method storage reduction by nonlinear iteration , 1983 .

[3]  Han Gyu Joo,et al.  Practical numerical reactor employing direct whole core neutron transport and subchannel thermal/hydraulic solvers , 2013 .

[4]  Mathieu Hursin.  Full Core, Heterogeneous, Time Dependent Neutron Transport Calculations with the 3D Code DeCART , 2010 .

[5]  N. Z. Cho,et al.  Fusion of method of characteristics and nodal method for 3-D whole-core transport calculation , 2002 .

[6]  Benoit Forget,et al.  Improved diffusion coefficients generated from Monte Carlo codes , 2013 .

[7]  Blake W. Kelley.  An Investigation of 2D/1D Approximations to the 3D Boltzmann Transport Equation. , 2015 .

[8]  Tara M. Pandya,et al.  Implementation, capabilities, and benchmarking of Shift, a massively parallel Monte Carlo radiation transport code , 2016, J. Comput. Phys..

[9]  Han Gyu Joo,et al.  Methods and performance of a three-dimensional whole-core transport code DeCART , 2004 .

[10]  Joseph W. Michels,et al.  The Pennsylvania State University Kaminaljuyu Project , 1973 .

[11]  T. M. Sutton,et al.  Diffusion theory methods for spatial kinetics calculations , 1994 .

[12]  A. Yamamoto,et al.  Derivation of Optimum Polar Angle Quadrature Set for the Method of Characteristics Based on Approximation Error for the Bickley Function , 2007 .

[13]  Brendan Kochunas,et al.  Overview of development and design of MPACT: Michigan parallel characteristics transport code , 2013 .

[14]  Ryan G. McClarren.  Theoretical Aspects of the Simplified Pn Equations , 2010 .

[15]  K. Smith,et al.  ASSEMBLY HOMOGENIZATION TECHNIQUES FOR LIGHT WATER REACTOR ANALYSIS , 1986 .

[16]  Thomas J. Downar,et al.  High-Fidelity Light Water Reactor Analysis with the Numerical Nuclear Reactor , 2007 .

[17]  R. Roy,et al.  DEVELOPMENT AND PARALLELIZATION OF THE THREE-DIMENSIONAL CHARACTERISTICS SOLVER MCI OF DRAGON , 2002 .

[18]  Brendan Kochunas,et al.  Assessment of the 2D MOC solver in MPACT: Michigan parallel characteristics transport code , 2013 .

[19]  Edward W. Larsen,et al.  A consistent 2D/1D approximation to the 3D neutron transport equation , 2015 .

[20]  R. Roy,et al.  A new characteristics algorithm for 3D transport calculations , 2003 .

[21]  Dorothea Wiarda,et al.  Development of a New 47-Group Library for the CASL Neutronics Simulators , 2015 .

[22]  Brendan Matthew Kochunas,et al.  A Hybrid Parallel Algorithm for the 3-D Method of Characteristics Solution of the Boltzmann Transport Equation on High Performance Compute Clusters. , 2013 .

[23]  Bradley T Rearden,et al.  Monte Carlo Criticality Methods and Analysis Capabilities in SCALE , 2011 .

[24]  M. Avramova,et al.  CTF Theory Manual , 2016 .

[25]  Simone Santandrea,et al.  APOLLO2 YEAR 2010 , 2010 .

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