• Nuclear Engineering and Technology
• /
• v.52 no.12
• /
• pp.2733-2742
• /
• 2020

This paper presents an assessment of applicability of the multigroup cross sections generated with Monte Carlo tools to the fast reactor analysis based on transport calculations. 33-group cross section sets were generated for simple one- (1-D) and two-dimensional (2-D) sodium-cooled fast reactor problems using the SERPENT code and applied to deterministic steady-state and depletion calculations. Relative to the reference continuous-energy SERPENT results, with the transport corrected P₀ scattering cross section, the k-eff value was overestimated by 506 and 588 pcm for 1-D and 2-D problems, respectively, since anisotropic scattering is important in fast reactors. When the scattering order was increased to P₅, the 1-D and 2-D problem errors were increased to 577 and 643 pcm, respectively. A sensitivity and uncertainty analysis with the PERSENT code indicated that these large k-eff errors cannot be attributed to the statistical uncertainties of cross sections and they are likely due to the approximate anisotropic scattering matrices determined by scalar flux weighting. The anisotropic scattering cross sections were alternatively generated using the MC²-3 code and merged with the SERPENT cross sections. The mixed cross section set consistently reduced the errors in k-eff, assembly powers, and nuclide densities. For example, in the 2-D calculation with P3 scattering order, the k-eff error was reduced from 634 pcm to -223 pcm. The maximum error in assembly power was reduced from 2.8% to 0.8% and the RMS error was reduced from 1.4% to 0.4%. The maximum error in the nuclide densities at the end of 12-month depletion that occurred in ²³⁷Np was reduced from 3.4% to 1.5%. The errors of the other nuclides are also reduced consistently, for example, from 1.1% to 0.1% for ²³⁵U, from 2.2% to 0.7% for ²³⁸Pu, and from 1.6% to 0.2% for ²⁴¹Pu. These results indicate that the scalar flux weighted anisotropic scattering cross sections of SERPENT may not be adequate for application to fast reactors where anisotropic scattering is important.

• Nuclear Engineering and Technology
• /
• v.50 no.3
• /
• pp.340-355
• /
• 2018

The methods and performance of a fast reactor multigroup cross section (XS) generation code EXUS-F are described that is capable of directly processing Evaluated Nuclear Data File format nuclear data files. RECONR of NJOY is used to generate pointwise XS data, and Doppler broadening is incorporated by the Gauss-Hermite quadrature method. The self-shielding effect is incorporated in the ultrafine group XSs in the resolved and unresolved resonance ranges. Functions to generate scattering transfer matrices and fission spectrum matrices are realized. The extended transport approximation is used in zero-dimensional calculations, whereas the collision probability method and the method of characteristics are used for one-dimensional cylindrical geometry and two-dimensional hexagonal geometry problems, respectively. Verification calculations are performed first for various homogeneous mixtures and cylindrical problems. It is confirmed that the spectrum calculations and the corresponding multigroup XS generations are performed adequately in that the reactivity errors are less than 50 pcm with the McCARD Monte Carlo solutions. The nTRACER core calculations are performed with the EXUS-F-generated 47 group XSs for the two-dimensional Advanced Burner Reactor 1000 benchmark problem. The reactivity error of 160 pcm and the root mean square error of the pin powers of 0.7% indicate that EXUF-F generates properly the broad-group XSs.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2002.10a
• /
• pp.49.1-49
• /
• 2002

• Proceedings of the Korean Nuclear Society Conference
• /
• 2005.05a
• /
• pp.117-118
• /
• 2005

• Nuclear Engineering and Technology
• /
• v.55 no.7
• /
• pp.2474-2482
• /
• 2023

The development of fast reactors with complex designs and operation status requires more accurate and effective simulation. The Monte-Carlo method can generate multi-group cross-sections in arbitrary geometry without approximation on resonances treatment and leads to good results in combination with diffusion codes. However, in previous studies, the coupling of Monte-Carlo generated multi-group cross-sections (MC-MGXS) and transport solvers has shown relatively large biases in fast reactor problems. In this paper, the main contribution to the biases is proved to be the neglect of the angle-dependence of the total cross-sections. The flux-moment homogenization technique (MHT) is proposed to take into account this dependence. In this method, the angular dependence is attributed to the transfer cross-sections, keeping an independent form for the total sections. For the MET-1000 benchmark, the multi-group transport simulation results with MC-MGXS generated with MHT are improved by 700 pcm and an additional 120 pcm with higher order scattering. The factors that cause the residual bias are discussed. The core power distribution bias is also significantly reduced when MHT is used. It proves that the MCMGXS with MHT can be applicable with transport solvers in fast reactor analysis.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1287-1300
• /
• 2017

Based on the observation that ignoring the angle dependency of multigroup resonance cross sections within a fuel pellet would result in nontrivial underestimation of the spatial self-shielding of flux, a parametrized spectral superhomogenization (SPH) factor library (PSSL) method is developed as a practical means of resolving the problem. Region-wise spectral SPH factors are calculated by the normal and transport corrected SPH iterations after ultrafine group slowing down calculations over various light water reactor pin-cell configurations. The parametrization is done with fuel temperature, U-238 number density, fuel radius, moderator source represented by ${\Sigma}_{mod}V_{mod}$, and the number density ratio of resonance nuclides to that of U-238 in a form of resonance interference correction factors. The parametrization is successful in that the root mean square errors of the interpolated SPH factors over the fuel regions of various pin-cells are within 0.1%. The improvement in reactivity error of the PSSL method is shown to be superior to that by the original SPH method in that the reactivity bias of -200 pcm to -300 pcm vanishes almost completely. It is demonstrated that the environment effect takes only about 4% in the reactivity improvement so that the pin-cell based PSSL method is effective in the assembly problems.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1209-1217
• /
• 2019

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.

• Nuclear Engineering and Technology
• /
• v.55 no.3
• /
• pp.1084-1096
• /
• 2023

On-the-fly(OTF) generation of Spectral Superhomogenization(SSPH) factors is analyzed in the multigroup(MG) whole core calculation employing pin-wise continuous energy(CE) slowing-down solutions. The motivation for the work is to avoid the huge computing time required for the generation of a parametrized SSPH factor library(PSSL) which is used to resolve the angular dependency of MG resonance cross sections, and also to exploit the advantage of flexible choice of a MG structure by using CE slowing-down solutions. Two pin-wise CE slowing-down methods, the equivalent Dancoff cell method and the shadowing effect correction method, are evaluated with the OTF SSPH method. The effectiveness of the OTF SSPH method is examined for various simplified and realistic core problems with various MG structures. It is demonstrated that the computing time overhead of this method is negligible whereas the solution accuracy is considerably enhanced.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1240-1249
• /
• 2017

An essential component of the neutron transport solver is the resonance self-shielding calculation used to determine equivalence cross sections. The neutron transport code, MPACT, is currently using the subgroup self-shielding method, in which the method of characteristics (MOC) is used to solve purely absorbing fixed-source problems. Recent efforts incorporating multigroup kernels to the MOC solvers in MPACT have reduced runtime by roughly $2{\times}$. Applying the same concepts for self-shielding and developing a novel lumped parameter approach to MOC, substantial improvements have also been made to the self-shielding computational efficiency without sacrificing any accuracy. These new multigroup and lumped parameter capabilities have been demonstrated on two test cases: (1) a single lattice with quarter symmetry known as VERA (Virtual Environment for Reactor Applications) Progression Problem 2a and (2) a two-dimensional quarter-core slice known as Problem 5a-2D. From these cases, self-shielding computational time was reduced by roughly $3-4{\times}$, with a corresponding 15-20% increase in overall memory burden. An azimuthal angle sensitivity study also shows that only half as many angles are needed, yielding an additional speedup of $2{\times}$. In total, the improvements yield roughly a $7-8{\times}$ speedup. Given these performance benefits, these approaches have been adopted as the default in MPACT.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1204-1212
• /
• 2024

To accurately calculate the heating distribution of the fast reactor, a neutron-photon library in MATXS format named Knight-B7.1-1968n × 94γ was processed based on the ENDF/B-VII.1 library for ultrafine groups. The neutron cross-section processing code MGGC2.0 was used to generate few-group neutron cross sections in ISOTXS format. Additionally, the self-developed photon cross-section processing code NGAMMA was utilized to generate photon libraries for neutron-photon coupled heating calculations, including photo-atom cross sections for the ISOTXS format, prompt photon production cross sections, and kinetic energy release in materials (KERMA) factors for neutrons and photons, and the self-shielding effect from the capture and fission cross sections of neutron to photon have been taken into account when the photon source generated by neutron is calculated. The interface code GSORCAL was developed to generate the photon source distribution and interface with the DIF3D code to calculate the neutron-photon coupling heating distribution of the fast reactor core. The neutron-photon coupled heating calculation route was verified using the ZPPR-9 benchmark and the RBEC-M benchmark, and the results of the coupled heating calculations were analyzed in comparison with those obtained from the Monte Carlo code MCNP. The calculations show that the library was accurately processed, and the results of the fast reactor neutron-photon coupled heating calculations agree well with those obtained from MCNP.