• 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.52 no.11
• /
• pp.2431-2441
• /
• 2020

Currently, the pin-by-pin homogenized solvers are a very active research field as they can, unlike the nodal codes, directly predict the local power, while requiring significantly less computational resources than the heterogeneous transport codes. This paper presents a recently developed pin-by-pin diffusion/SP3 solver Tortin, its spatial discretization method and the reflector treatment. Regarding the spatial discretization, it was observed that the finite difference method applied on pin-cell size mesh does not properly capture the big flux change between MOX and uranium fuel, while the nodal expansion method is more accurate but too slow. If the finite difference method is used with a finer mesh in the outer two pin rows of the fuel assembly, it increases the required computation time by only 50%, but decreases the pin power errors below 1% with respect to lattice code reference solutions. The paper further describes the coupling of Tortin with a microscopic depletion solver. Several verification tests show that the SP3 pin-by-pin solver can reproduce the heterogeneous transport solvers results with very good accuracy, even for fuel cycle depletion of very heterogeneous core employing MOX fuel or inserted control rods, while being two orders of magnitude faster.