Superhomogenization Method Research Articles

FSAR: Full neutron spectrum code for advanced reactor simulation

The full neutron spectrum code for advanced reactor simulation named FSAR has recently been developed at Nuclear Power Institute of China in order to meet the requirements of advanced reactor with large neutron energy range. Based on the two-step calculation scheme, FSAR consists of two-dimensional lattice calculation code and core calculation code. In two-dimensional lattice calculation, the subgroup method with ultrafine energy groups was implemented in the two-dimensional resonance self-shielding by incorporating the MOC to get the accurate self-shielded cross-section. For better consideration of the strong space coupling in different geometric sizes due to different mean free path of neutrons, the super-homogenization method and the leakage model were applied. In the core calculation, the discrete-ordinate method and micro burnup calculation method were used to simulate the core neutron transport and depletion. Preliminary calculation results showed that for the problems with wide spectrum, the self-shielded cross-sections have a good agreement with the Monte Carlo solution. The results shown in this paper indicate that FSAR has good performances of the cross-section generation in full neutron spectrum problem simulation.

DEVELOPMENT OF THE MPACT 69-GROUP LIBRARY FOR MAGNOX REACTOR ANALYSIS USING CASL VERA

The Consortium for Advanced Simulation of Light Water Reactors (CASL) has developed the CASL toolset, Virtual Environment for Reactor Analysis (VERA), for pressurized water reactor (PWR) analysis. Recently the CASL VERA was improved for Magnox reactor analysis, which required the development of a new cross section library and new geometrical and thermal feedback capabilities for graphite-moderated Magnox reactors. The MPACT neutronics module of the CASL core simulator is a 3D whole core transport code, which requires a new cross section library with a different energy group structure due to the different neutronic characteristics of Magnox compared with PWR. A new 69-group structure was developed based on the MPACT 51-group structure to have more thermal energy groups and to be a subset of the SCALE 252-group structure. The ENDF/B-VII.1 MPACT 69-group library was developed for Magnox reactor analysis using the SCALE/AMPX and VERA-XSTools for which a super-homogenization method was applied, and transport cross sections were generated for graphite using a neutron leakage conservation method. Benchmark results show that new MPACT 69-group library works reasonably well for Magnox reactor analysis.

On the equivalence of reaction rate in energy collapsing of fast reactor code SARAX

Scattering resonance of medium mass nuclides leads complex spectrum in the fast reactor, which requires thousands of energy groups in the spectrum calculation. When the broad-group cross sections are collapsed, reaction rate cannot be completely conserved. To eliminate the error from energy collapsing, the Super-homogenization method in energy collapsing (ESPH) was employed in the fast reactor code SARAX. An ESPH factor was derived based on the ESPH-corrected SN transport equation. By applying the factor in problems with reflective boundary condition, both the effective multiplication factor and reaction rate were conserved. The fixed-source iteration was used to ensure the stability of ESPH iteration. However, in the energy collapsing process of SARAX, the vacuum boundary condition was adopted, which was necessary for fast reactors with strong heterogeneity. To further reduce the error caused by leakage, an additional conservation factor was proposed to correct the neutron current in energy collapsing. To evaluate the performance of ESPH with conservation factor, numerical benchmarks of fast reactors were calculated. The results of broad-group calculation agreed well with the direct full-core Monte-Carlo calculation, including the effective multiplication factor, radial power distribution, total control rod worth and sodium void worth.

Generalized Equivalence Theory Used with Spatially Linear Sources in the Method of Characteristics for Neutron Transport

A recent hybrid stochastic-deterministic calculation scheme using Monte Carlo–tallied group cross sections in a deterministic solver uses the best of both worlds for accurate and fast reactor agnostic transport simulations. However, neglecting the angular dependence of group cross sections induces large self-shielding errors in resonance groups, causing a large reactivity bias up to 300 pcm in light water reactors. To recover this error, we introduce a two-scale assembly transport calculation scheme: cross sections are tallied at the assembly level, while equivalence parameters are computed in a two-dimensional (2-D) pin cell system. We validate a novel equivalence method based on jump conditions on angular fluxes by comparing to the well-established superhomogenization method for 2-D and three-dimensional (3-D) linear source method of characteristics calculations. Test cases include 2-D and 3-D assemblies of two different enrichments with homogeneous and discretized cross-section discretizations. The linear source approximation enables using coarse source-region discretization for these hot zero-power problems. Both equivalence techniques perform similarly, recover the reactivity bias, and achieve near preservation of reaction rates, supporting this multiscale approach to equivalence.

Development of the multigroup cross section library for the CASL neutronics simulator MPACT: Method and procedure

The MPACT neutronics module of the Consortium for Advanced Simulation of Light Water Reactors (CASL) core simulator is a 3D whole core transport code being developed for the CASL toolset known as the Virtual Environment for Reactor Analysis (VERA). Key characteristics of the MPACT code include a subgroup method for resonance self-shielding and a whole-core transport solver with a 2D/1D synthesis method. The AMPX and SCALE code packages have been significantly improved to include intermediate resonance parameters, as well as new modules to enhance the accuracy of resonance data and scattering matrices for the Bondarenko approach and transport cross sections of hydrogen. In addition, CASL XSTools have been developed to generate the MPACT multigroup cross section library. A simple super-homogenization method for 238U has been developed to resolve reaction rate discrepancy issues caused by angle-dependent total cross sections, poor scattering matrices, and a poor resonance interference model at coarse group structure. Due to the limitation of computing capacity that can occur even in high performance computing, a new 51-group structure has been developed for efficient simulation that is applicable for both pressurized water reactor (PWR) and boiling water reactor (BWR) simulations. A new ENDF/B-VII.1 MPACT 51-group library has been developed using the AMPX/SCALE and XSTools code packages, and it has been verified and validated by performing various benchmark calculations.

One-step Monte Carlo global homogenization based on RMC code

Abstract Due to the limitation of the computers, the conventional homogenization method is based on many assumptions and approximations, and some tough problems such as energy spectrum and boundary condition are faced. To deal with those problems, the Monte Carlo global homogenization is adopted. The Reactor Monte Carlo code RMC is used to study the global homogenization method, and the one-step global homogenization method is proposed. The superimposed mesh geometry is also used to divide the physical models, leading to better geometric flexibility. A set of multigroup homogenization cross sections is online generated for each mesh under the real neutron energy spectrum and boundary condition, the cross sections are adjusted by the superhomogenization method, and no leakage correction is required. During the process of superhomogenization, the author-developed reactor core program NLSP3 is used for global calculation, so the global flux distribution and equivalent homogenization cross sections could be solved simultaneously. Meanwhile, the calculated homogenization cross section could accurately reconstruct the non-homogenization flux distribution and could also be used for fine calculation. This one-step global homogenization method was tested by a PWR assembly and a small reactor model, and the results show the validity.

Utilization of Regionwise Even-Parity Discontinuity Factor to Reduce Discretization Error of MOC

To reduce angular and spatial discretization error of the method of characteristics with a coarse calculation condition, the regionwise even-parity discontinuity factor (EPDF) for transport calculations is evaluated through an iterative procedure using only the regionwise scalar flux, i.e., without the odd-parity angular flux, the partial current, or the net current at the region boundary. The regionwise EPDF is evaluated in a single-assembly geometry with the reflective boundary condition. The evaluated EPDF is applied to a 2 × 2 colorset assembly and core configurations, and the performance is compared to that of the conventional superhomogenization (SPH) method. The calculation results indicate that (1) no convergence issue is observed during the iteration process to estimate the EPDF, (2) the performance of the regionwise EPDF is better than that of the conventional SPH method, and (3) the normalization of the EPDF is necessary to incorporate different surface scalar flux levels among different types of fuel assemblies.

Application of the SPH method to account for the angular dependence of multigroup resonant cross sections in thermal reactor calculations

The issue of the angular dependence of the multigroup total cross sections in the multigroup transport equation solution was recognized a long time ago. Only recently has its impact on resonant self-shielding been identified. A viable solution for this problem is the superhomogenization (SPH) corrections. In the present work, two applications of the SPH method are presented, one applying to the APOLLO3® deterministic calculations with the multigroup cross sections condensed from the TRIPOLI-4® Monte Carlo simulation results and the other applying to a subgroup method based on the fine-structure equation. The numerical results show the effectiveness of the SPH method in dealing with this difficulty in both situations. This work also demonstrates that the subgroup method based on the fine-structure equation plus the SPH corrections is a practical, effective and accurate resonant self-shielding solution for the PWR fuel pin-cell and fuel assembly calculations.

A Newton solution for the Superhomogenization method: The PJFNK-SPH

This work presents two novel topics regarding the Superhomogenization method: 1) the formalism for the implementation of the method with the linear Boltzmann Transport Equation, and 2) a Newton algorithm for the solution of the nonlinear problem that arises from the method. These new ideas have been implemented in a continuous finite element discretization in the MAMMOTH reactor physics application. The traditional solution strategy for this nonlinear problem uses a Picard, fixed-point iterative process whereas the new implementation relies on MOOSE’s Preconditioned Jacobian-Free Newton Krylov method to allow for a direct solution. The PJFNK-SPH can converge problems that were either intractable or very difficult to converge with the traditional iterative approach, including geometries with reflectors and vacuum boundary conditions. This is partly due to the underlying Scalable Nonlinear Equations Solvers in PETSc, which are integral to MOOSE and offer Newton damping, line search and trust region methods. The PJFNK-SPH has been implemented and tested for various discretizations of the transport equation included in the Rattlesnake transport solver. Speedups of five times for diffusion and ten to fifteen times for transport were obtained when compared to the traditional Picard approach. The three test problems cover a wide range of applications including a standard Pressurized Water Reactor lattice with control rods, a Transient Reactor Test facility control rod supercell and a prototype fast-thermal reactor. The reference solutions and initial cross sections were obtained from the Serpent 2 Monte Carlo code. The SPH-corrected cross sections yield eigenvalues that are near exact, relative to reference solutions, for reflected geometries, even with reflector regions. In geometries with vacuum boundary conditions the accuracy is problem dependent and solutions can be within a few to a few hundred pcm. The root-mean-square error in the power distribution is below 0.8% of the reference Monte Carlo. There is little benefit from SPH-corrected transport in typical scoping calculations, but for more detailed analyses it can yield superior convergence of the solution in some of the test problems. This PJFNK-SPH approach is currently being used in the modeling of the Transient Test Reactor at Idaho National Laboratory, where full reactor core SPH-corrected cross sections are employed to reduce the homogenization errors in transient or multi-physics calculations. This base implementation of the PJFNK-SPH provides an extremely robust solver and a springboard to further improve the Superhomogenization method in order to better preserve neutron currents, one of the primary deficiencies of the method.

Super homogenization method applied in diffusion calculation for pebble-bed fluoride-salt-cooled high-temperature reac

Super homogenization method applied in diffusion calculation for pebble-bed fluoride-salt-cooled high-temperature reac

Reduction of MOC discretization errors through a minimization of source ratio variances

ABSTRACTA new technique to reduce discretization errors for ray tracing in the method of characteristics (MOC) is proposed focusing on depletion calculations of single and multi-assembly geometries. In order to efficiently carry out depletion calculations, a calculation scheme using the superhomogenization (SPH) method can be used. However, the discretization errors are caused by changes of neutron sources and total cross sections according to a depletion. This fact means that improvement of accuracy cannot be expected by the calculation scheme with the SPH method when changes of the above parameters are significant. In order to mitigate this problem, a new approach is developed. In the new approach, the discretization errors are reduced by minimizing a variance of a certain parameter which is composed of a ratio of neutron source to total cross section. The verification results suggest that accuracy is degraded by the SPH method as expected especially in a geometry where neutron sources and total cross sections are drastically changing through a depletion. On the other hand, the new approach gives more accurate results compared to the conventional MOC in all calculation cases. Consequently, improvement of calculation efficiency by the new approach is confirmed.

Equivalence treatment in homogenization via Monte Carlo method

The multi-group constants generated via continuous energy Monte Carlo method do not satisfy the equivalence between heterogeneous and homogeneous calculation. To meet the requirement of equivalence, the general equivalence theory (GET) and the superhomogenization method (SPH) are applied to the Monte Carlo based group constants. Both of the methods improved accuracy. Inheriting the essences of GET and SPH, the super equivalence method (SPE) is developed in this paper to overcome the limitations of the two conventional equivalence techniques. The three methods have been implemented into the Monte Carlo based homogenization code MCMC. Test calculations were carried out in a simplified PWR reactor core and the C5G7 MOX fuel benchmark. The study has shown that the equivalence treatment is necessary in Monte Carlo homogenization calculations and the SPE method can meet the requirement.

05/00876 Entropy generation due to non-Newtonian flow around a spherical particle at low Reynolds number

In this paper, performances of spatial homogenization methods for fuel or non-fuel cells are compared in slab geometry in order to facilitate pin-by-pin core calculations. Since the spatial homogenization methods were mainly developed for fuel assemblies, systematic study of their performance for the cell-level homogenization has not been carried out. Importance of cell-level homogenization is recently increasing since the pin-by-pin mesh core calculation in actual three-dimensional geometry, which is less approximate approach than current advanced nodal method, is getting feasible. Four homogenization methods were investigated in this paper; the flux-volume weighting, the generalized equivalence theory, the superhomogenization (SPH) method and the nonlinear iteration method. The last one, the nonlinear iteration method, was tested as the homogenization method for the first time. The calculations were carried out in simplified colorset assembly configurations of PWR, which are simulated by slab geometries, and homogenization performances were evaluated through comparison with the reference cell-heterogeneous calculations. The calculation results revealed that the generalized equivalence theory showed best performance. Though the nonlinear iteration method can significantly reduce homogenization error, its performance was not as good as that of the generalized equivalence theory. Through comparison of the results obtained by the generalized equivalence theory and the superhomogenization method, important byproduct was obtained; deficiency of the current superhomogenization method, which could be improved by incorporating the “cell-level discontinuity factor between assemblies”, was clarified.

Cell homogenization methods for pin-by-pin core calculations tested in slab geometry

In this paper, performances of spatial homogenization methods for fuel or non-fuel cells are compared in slab geometry in order to facilitate pin-by-pin core calculations. Since the spatial homogenization methods were mainly developed for fuel assemblies, systematic study of their performance for the cell-level homogenization has not been carried out. Importance of cell-level homogenization is recently increasing since the pin-by-pin mesh core calculation in actual three-dimensional geometry, which is less approximate approach than current advanced nodal method, is getting feasible. Four homogenization methods were investigated in this paper; the flux-volume weighting, the generalized equivalence theory, the superhomogenization (SPH) method and the nonlinear iteration method. The last one, the nonlinear iteration method, was tested as the homogenization method for the first time. The calculations were carried out in simplified colorset assembly configurations of PWR, which are simulated by slab geometries, and homogenization performances were evaluated through comparison with the reference cell-heterogeneous calculations. The calculation results revealed that the generalized equivalence theory showed best performance. Though the nonlinear iteration method can significantly reduce homogenization error, its performance was not as good as that of the generalized equivalence theory. Through comparison of the results obtained by the generalized equivalence theory and the superhomogenization method, important byproduct was obtained; deficiency of the current superhomogenization method, which could be improved by incorporating the “cell-level discontinuity factor between assemblies”, was clarified.