A Neutron Transport Characteristics Method for 3D Axially Extruded Geometries Coupled with a Fine Group Self-Shielding Environment

Abstract In this paper we describe some recent developments in the Method of Characteristics (MOC) for three-dimensional (3D) extruded geometries in the nuclear reactor analysis code APOLLO3®. We discuss the parallel strategies implemented for the transport sweep of the MOC solver in the OpenMP framework, and introduce the 3D version of the DPN operator that is customarily used in APOLLO2 to accelerate MOC convergence. In order to provide good physical results, we have also coupled the MOC with the self-shielding environment of APOLLO3®. We describe, in particular, the coupling techniques necessary to implement a full subgroup cross-section self-shielding method and a specialized version of the Tone self-shielding technique. In this framework, we use part of the tracking method used for the 3D calculation to provide the two-dimensional Collision Probability Method (CPM) coefficients necessary to produce the self-shielding calculations. We will show some important computational speedups also in the CPM of APOLLO3® with respect to the APOLLO2 CPM equivalent implementation, including the parallelization issue. Finally, we will compare our approach toward a Monte Carlo calculation of a fast breeder reactor hexagonal assembly representing a fertile-fissile interface.

[1]  A. Hébert.  High-Order Diamond Differencing along Finite Characteristics , 2011 .

[2]  Richard Sanchez,et al.  PROSPECTS IN DETERMINISTIC THREE-DIMENSIONAL WHOLE-CORE TRANSPORT CALCULATIONS , 2012 .

[3]  D. Sciannandrone.  Acceleration and higher order schemes of a characteristic solver for the solution of the neutron transport equation in 3D axial geometries , 2015 .

[4]  Fausto Malvagi,et al.  TRIPOLI-4®, CEA, EDF and AREVA Reference Monte Carlo Code , 2015, ICS 2014.

[5]  Hugo J. Hahn.  ORGANISATION FOR ECONOMIC CO-OPERATION AND DEVELOPMENT, NUCLEAR ENERGY AGENCY , 1983 .

[6]  Mohamed Dahmani,et al.  Parallel Solver Based on the Three-Dimensional Characteristics Method: Design and Performance Analysis , 2005 .

[7]  Alain Hébert,et al.  Applied reactor physics , 2009 .

[8]  K. Kim,et al.  Unstructured partial- and net-current based coarse mesh finite difference acceleration applied to the extended step characteristics method in NEWT , 2011 .

[9]  L. Cao,et al.  A 2D/1D coupling neutron transport method based on the matrix MOC and NEM methods , 2013 .

[10]  K. Kim,et al.  Axial SPN and Radial MOC Coupled Whole Core Transport Calculation , 2007 .

[12]  N. Cho,et al.  Whole-Core Transport Solutions with 2-D/1-D Fusion Kernel via p-CMFD Acceleration and p-CMFD Embedding of Nonoverlapping Local/Global Iterations , 2015 .

[13]  A. Hébert,et al.  Tracking algorithms for multi-hexagonal assemblies (2D and 3D) , 2014 .

[14]  S. Kosaka,et al.  Transport Theory Calculation for a Heterogeneous Multi-Assembly Problem by Characteristics Method with Direct Neutron Path Linking Technique , 2000 .

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

[17]  Eugene P. Wigner,et al.  The Physical Theory of Neutron Chain Reactors , 1958 .

[18]  Simone Santandrea,et al.  Analysis and improvements of the DPN acceleration technique for the method of characteristics in unstructured meshes , 2005 .

[19]  Jasmina L. Vujic,et al.  ANEMONA — a neutron transport code for general geometry reactor assemblies based on the method of characteristics and R-function solid modeler , 2001 .

[20]  J. Taylor,et al.  A METHOD OF CHARACTERISTICS SOLUTION TO THE OECD/NEA 3D NEUTRON TRANSPORT BENCHMARK PROBLEM , 2007 .

[21]  S. Santandrea,et al.  A Linear Surface Characteristics Approximation for Neutron Transport in Unstructured Meshes , 2008 .

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

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

[24]  Liangzhi Cao,et al.  A new three-dimensional method of characteristics for the neutron transport calculation , 2011 .

[25]  Andrew R. Siegel,et al.  Parallel performance results for the OpenMOC neutron transport code on multicore platforms , 2016, Int. J. High Perform. Comput. Appl..

[26]  M. Williams.  Correction of multigroup cross sections for resolved resonance interference in mixed absorbers , 1983 .

[27]  Gérald Rimpault,et al.  THE ERANOS CODE AND DATA SYSTEM FOR FAST REACTOR NEUTRONIC ANALYSES , 2002 .

[28]  R. Wiersma,et al.  The Collision Probability Method in Today’s Computer Environment , 2016 .

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

[30]  Richard Sanchez,et al.  Approximate Solutions of the Two-Dimensional Integral Transport Equation by Collision Probability Methods , 1977 .

[31]  Nam Zin Cho,et al.  2D/1D fusion method solutions of the three-dimensional transport OECD benchmark problem C5G7 MOX , 2006 .

[32]  Igor Zmijarevic,et al.  New Numerical Solution with the Method of Short Characteristics for 2-D Heterogeneous Cartesian Cells in the APOLLO2 Code: Numerical Analysis and Tests , 2009 .

[33]  Norman J. McCormick,et al.  A Review of Neutron Transport Approximations , 1982 .

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

[35]  Brendan Kochunas,et al.  Parallel 3-D method of characteristics in MPACT , 2013 .

[36]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[37]  François Févotte,et al.  Vectorization of a 2D–1D Iterative Algorithm for the 3D Neutron Transport Problem in Prismatic Geometries , 2014, ICS 2014.

[38]  Pietro Mosca.  Conception et développement d'un mailleur énergétique adaptatif pour la génération des bibliothèques multigroupes des codes de transport , 2009 .

[39]  Makoto Ishikawa.  Consistency evaluation of JUPITER experiment and analysis for large FBR cores , 1996 .

[40]  Nader Satvat,et al.  MOCUM: A two-dimensional method of characteristics code based on constructive solid geometry and unstructured meshing for general geometries , 2012 .

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

[42]  S. Santandrea,et al.  Optimized tracking strategies for step MOC calculations in extruded 3D axial geometries , 2016 .

[43]  Development of Axially Simplified Method of Characteristics in Three-Dimensional Geometry , 2008 .

[44]  Simone Santandrea,et al.  An Integral Multidomain DPN Operator as Acceleration Tool for the Method of Characteristics in Unstructured Meshes , 2007 .

[45]  F. Fevotte,et al.  MICADO: Parallel implementation of a 2D-1D iterative algorithm for the 3D neutron transport problem in prismatic geometries , 2013 .

[46]  Edward W. Larsen,et al.  2D/1D approximations to the 3D neutron transport equation. I: Theory , 2013 .

[47]  Edward W. Larsen,et al.  Fast iterative methods for discrete-ordinates particle transport calculations , 2002 .

[48]  R. Lenain,et al.  APOLLO3: a common project of CEA, AREVA and EDF for the development of a new deterministic multi-purpose code for core physics analysis , 2009 .

[49]  Tatsuzo Tone,et al.  A Numerical Study of Heterogeneity Effects in Fast Reactor Critical Assemblies , 1975 .

[50]  M. J. Abbate,et al.  Methods of Steady-State Reactor Physics in Nuclear Design , 1983 .

[51]  Robert Roy,et al.  Solving Three-Dimensional Large-Scale Neutron Transport Problems Using Hybrid Shared-Distributed Parallelism and Characteristics Method , 2007 .

[52]  Shanjie Xiao.  Hardware accelerated high performance neutron transport computation based on AGENT methodology , 2009 .

[53]  William R. Martin,et al.  New approach to synthetic acceleration of transport calculations , 1988 .

[54]  Toshikazu Takeda,et al.  Simple and Efficient Parallelization Method for MOC Calculation , 2010 .

[55]  Jean-François Sauvage,et al.  Sodium-cooled fast reactors: the ASTRID plant project , 2011 .

[56]  Bradley T Rearden,et al.  Methods and issues for the combined use of integral experiments and covariance data: Results of a NEA international collaborative study , 2014 .

[57]  Toshikazu Takeda,et al.  Verification of 3D Heterogeneous Core Transport Calculation Utilizing Non-linear Iteration Technique , 2004 .

[58]  Benoit Forget,et al.  The OpenMOC method of characteristics neutral particle transport code , 2014 .

[59]  Simone Santandrea,et al.  Treatment of Boundary Conditions in Trajectory-Based Deterministic Transport Methods , 2002 .

[60]  G. Chiba,et al.  A Combined Method to Evaluate the Resonance Self Shielding Effect in Power Fast Reactor Fuel Assembly Calculation , 2003 .