Calculation of Higher-Order Fluxes in Symmetric Cores—II: Implementation

Abstract In neutron chain systems with material symmetries, various k-eigenvalues of the neutron balance equation beyond the dominant one may be degenerated. As shown in a companion paper, the power iteration method can be used to compute higher eigenfunctions in symmetric systems, provided that the global problem is partitioned into symmetry class—related lower-sized problems with appropriate boundary conditions. Those boundary conditions have been implemented in the diffusion solver of the ERANOS code system in rectangular geometry and within the framework of a discontinuous Galerkin spatial approximation of the multigroup discrete ordinates transport equation in the SNATCH solver. Numerical results in homogeneous geometry are provided for verification purposes, as well as the first eigenfunctions of the Takeda benchmarks. Finally, the transport effect on the first flux harmonics for an industrial-sized reactor ZPPR-18 is discussed.

[1]  J. Kuttler,et al.  EIGENVALUES OF THE LAPLACIAN ON REGULAR POLYGONS AND POLYGONS RESULTING FROM THEIR DISECTION , 1999 .

[2]  A. Sargeni,et al.  A computational technique for evaluating eigenfunctions of symmetrical nuclear systems , 1994 .

[3]  W. Martin,et al.  Theory and applications of the fission matrix method for continuous-energy Monte Carlo , 2014 .

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

[5]  T. Booth.  Power Iteration Method for the Several Largest Eigenvalues and Eigenfunctions , 2006 .

[6]  D. Fournier,et al.  Reactivity Perturbation Formulation for a Discontinuous Galerkin–Based Transport Solver and Its Use with Adaptive Mesh Refinement , 2011 .

[7]  Gumersindo Verdú,et al.  3D λ-modes of the neutron-diffusion equation , 1994 .

[8]  R. S. Modak,et al.  Use of a New Boundary Condition in Computational Neutron Transport , 2009 .

[9]  Damien Fournier,et al.  High-order discrete ordinate transport in hexagonal geometry: A new capability in ERANOS , 2010 .

[10]  Gérald Rimpault,et al.  Calculation of Higher-Order Fluxes in Symmetric Cores—I: Theory , 2016 .

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

[12]  Bruno Fontaine,et al.  Pre-conceptual design study of ASTRID core , 2012 .

[13]  Toshikazu Takeda,et al.  3-D Neutron Transport Benchmarks , 1991 .

[14]  Toshikazu Takeda,et al.  Analysis of First-Harmonic Eigenvalue Separation Experiments on KUCA Coupled-Core , 1998 .

[15]  T. Takeda,et al.  Evaluation of Eigenvalue Separation by the Monte Carlo Method , 2002 .