• Nuclear Engineering and Technology
• /
• v.53 no.5
• /
• pp.1403-1415
• /
• 2021

The transient capability of a SP₃ based pin-wise core analysis code SPHINCS is developed and verified through the analyses of the C5G7-TD benchmark. Spatial discretization is done by the fine mesh finite difference method (FDM) within the framework of the coarse mesh finite difference (CMFD) formulation. Pin size fine meshes are used in the radial fine mesh kernels. The time derivatives of the odd moments in the time-dependent SP₃ equations are neglected. The pin homogenized group constants and Super Homogenization (SPH) factors generated from the 2D single assembly calculations at the unrodded and rodded conditions are used in the transient calculations via proper interpolation involving the approximate flux weighting method for the cases that involve control rod movement. The simplifications and approximations introduced in SPHINCS are assessed and verified by solving all the problems of C5G7-TD and then by comparing with the results of the direct whole core calculation code nTRACER. It is demonstrated that SPHINCS yields accurate solutions in the transient behaviors of core power and reactivity.

• 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.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.98-109
• /
• 2016

The flow and heat transfer characteristics of a heat exchanger having rectangular pin-fin in the flow direction have been investigated numerically. On the bottom plate, the convective boundary conditions for the hot side was given, and the fins were arranged in a channel-type geometric model using the periodic boundary condition in the span-wise direction. Three-dimensional numerical calculations for the flow and conjugate heat transfer problem were conducted using SIMPLE algorithm and $k-{\varepsilon}$ turbulence model. For the slanted pin-fin models, it was found that the downward cooling flow is generated due to the downward pressure gradient component, which can enhance the heat transfer performance near the bottom surface and the fin stem region. Four different inclined angles were considered in the Reynolds number range of 13,500-55,000. The aero-thermal performance of the slanted pin-fin heat exchangers, such as the volume and area goodness factors, were summarized and compared with the baseline plate-fin type heat exchanger quantitatively.

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.1871-1885
• /
• 2019

The deterministic MOC code STREAM of the Computational Reactor Physics and Experiment (CORE) laboratory of Ulsan National Institute of Science and Technology (UNIST), was initially designed for the calculation of pressurized water reactor two- and three-dimensional assemblies and cores. Since fast reactors play an important role in the generation-IV concept, it was decided that the code should be upgraded for the analysis of fast neutron spectrum reactors. This paper presents a coupled code - TULIP/STREAM, developed for the fast reactor assembly and core calculations. The TULIP code produces self-shielded multi-group cross-sections using a one-dimensional cylindrical model. The generated cross-section library is used in the STREAM code which solves eigenvalue problems for a two-dimensional assembly and a multi-assembly whole reactor core. Multiplication factors and steady-state power distributions were compared with the reference solutions obtained by the continuous energy Monte-Carlo code MCS. With the developed code, a sensitivity study of the number of energy groups, the order of anisotropic P_N scattering, and the multi-group cross-section generation model was performed on the k_eff and power distribution. The 2D core simulation calculations show that the TULIP/STREAM code gives a k_eff error smaller than 200 pcm and the root mean square errors of the pin-wise power distributions within 2%.

• Nuclear Engineering and Technology
• /
• v.53 no.1
• /
• pp.44-60
• /
• 2021

In this paper, we validate the decay heat calculation capability via a two-step method to analyze spent nuclear fuel (SNF) discharged from pressurized water reactors (PWRs). The calculation method is implemented with a lattice code STREAM and a nodal diffusion code RAST-K. One of the features of this method is the direct consideration of three-dimensional (3D) core simulation conditions with the advantage of a short simulation time. Other features include the prediction of the isotope inventory by Lagrange non-linear interpolation and the use of power history correction factors. The validation is performed with 58 decay heat measurements of 48 fuel assemblies (FAs) discharged from five PWRs operated in Sweden and the United States. These realistic benchmarks cover the discharge burnup range up to 51 GWd/MTU, 23.2 years of cooling time, and spanning an initial uranium enrichment range of 2.100-4.005 wt percent. The SNF analysis capability of STREAM is also employed in the code-to-code comparison. Compared to the measurements, the validation results of the FA calculation with RAST-K are within ±4%, and the pin-wise results are within ±4.3%. This paper successfully demonstrates that the developed decay heat calculation method can perform SNF back-end cycle analyses.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3732-3753
• /
• 2023

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, developed based on the Multi-stage Response-function Transport (MRT) methodology, enables real-time simulation of nuclear systems such as reactor cores, spent nuclear fuel pools and casks, and sub-critical facilities. This paper presents the application of a novel fission matrix-based burnup methodology to the well-characterized JSI TRIGA Mark II research reactor. This methodology allows for calculation of nuclear fuel depletion by combination and interpolation of RAPID's burnup dependent fission matrix (FM) coefficients to take into account core changes due to burnup. The methodology is compared to experimentally validated Serpent-2 Monte Carlo depletion calculations. The results show that the burnup methodology for RAPID (bRAPID) implemented into RAPID is capable of accurately calculating the k_eff burnup changes of the reactor core as the average discrepancies throughout the whole burnup interval are 37 pcm. Furthermore, capability of accurately describing 3D fission source distribution changes with burnup is demonstrated by having less than 1% relative discrepancies compared to Serpent-2. Good agreement is observed for axially and pin-wise dependent fuel burnup and nuclear fuel nuclide composition as a function of burnup. It is demonstrated that bRAPID accurately describes burnup in areas with high gradients of neutron flux (e.g. vicinity of control rods). Observed discrepancies for some isotopes are explained by analyzing the neutron spectrum. This paper presents a powerful depletion calculation tool that is capable of characterization of spent nuclear fuel on the fly while the reactor is in operation.