This study introduces a transient analysis module developed for RAST-V and validates it using the Kalinin-3 benchmark problem. For the benchmark analysis, RAST-V standalone and STREAM/RAST-V ...calculations were performed. STREAM supplies the few-group constants and RAST-V conducts a 3D core simulation utilizing few-group cross-sectional data. To improve accuracy, the main solver was developed based on the advanced semi-analytic nodal method. To evaluate the computational capability of the transient analysis module in RAST-V, Kalinin-3 benchmark is employed. Kalinin-3 represents a coolant transient benchmark that offers experimental data during the deactivation of the Main Circulation Pumps. Consequently, the transient calculations reflected the changes in the reactor flow rate. Benchmark comprising steady-state and transient calculations. During the steady state, the STREAM/RAST-V combination demonstrated a 30 ppm root mean square difference from 0 to 128.50 EFPD. For the transient calculations, STREAM/RAST-V showed power differences within ±7 % over a range of 0–300 s. Axial offset differences were within ±3 %, and the RMS difference in radial power ranged within 2.596 % at both 0 and 300 s. Overall, this study effectively demonstrated the newly developed transient solver in RAST-V and validated it using the Kalinin-3 benchmark problem.
This paper introduces the creation of a module for Uncertainty Quantification (UQ) specifically designed for VVER analysis through the implementation of the STREAM/RAST-V two-step approach. The aim ...was to expand the range of use by developing a UQ module tailored for analyzing VVER. This research presents two innovative computational functionalities: (1) development of a library for the pin-based pointwise energy slowing down method (PSM), and (2) extension of the analysis area to study hexagonal-geometry fuel assemblies. The proposed UQ scheme was evaluated through verification using UAM benchmark, and comparative analysis between codes using SCALE 6.2.2 for. STREAM provides an accuracy comparable to that of SCALE 6.2.2. Additionally, a PSM covariance library was utilized in the calculations, achieving 0.7941% and 0.7907% accuracies in the hot full power and hot zero power calculations, respectively. To assess the UQ sequences in the two-step method, the STREAM/RAST-V calculation scheme was verified using the STREAM lattice code. To conclude, this study furnishes comprehensive insights into the development of the UQ module within the two-step method for VVER analysis, and it validates its performance through utilization of the UAM benchmark.
This paper presents a high-fidelity simulation of the Organization for Economic Co-operation and Development (OECD) Nuclear Energy Agency (NEA) 3D whole-core Watts Bar benchmark using the UNIST ...in-house STREAM3D (Steady State and Transient Reactor Analysis code with Method of Characteristics) neutronic code. The benchmark encompasses various whole-core exercises, including single physics problems, multi-physics simulations, and depletion problems. When comparing parameters during the zero-power physics tests, including ITC, DBW, CRW, and criticality tests, STREAM3D results indicate a strong agreement with the measured data and KENO-VI. The comparison with the MC21/CTF code in 3D HFP BOC condition demonstrated strong agreement, with only a 0.42% difference in the normalized radial power distribution, a 0.38 K difference in the RMS of the assembly coolant exit temperature, and a mere 4 ppm difference in CBC.
The China Experimental Fast Reactor (CEFR) is a small, sodium-cooled fast reactor with 20 MW(e) of power. Start-up tests of the CEFR were performed from 2010 to 2011. The China Institute of Atomic ...Energy made some of the neutronics start-up-test data available to the International Atomic Energy Agency (IAEA) as part of an international neutronics benchmarking exercise by distributing the experimental data to interested organizations from the member states of the IAEA. This benchmarking aims to validate and verify the physical models and neutronics simulation codes with the help of the recorded experimental data. The six start-up tests include evaluating criticality, control-rod worth, reactivity effects, and neutron spectral characteristics. As part of this coordinated research, the IAEA performed neutronics calculations using the Monte Carlo codes Serpent 2 and OpenMC, which can minimize modeling assumptions and produce reference solutions for code verification. Both codes model a three-dimensional heterogeneous core with an ENDF/B-VII.1 cross-section library. This study presents the calculation results with a well-estimated criticality and a reasonably good estimation of reactivities. The description and analysis of the core modeling assumptions, challenges in modeling a dense SFR core, results of the first phase of this project, and comparative analysis with measurements are presented.
STREAM, developed by the Computational Reactor Physics and Experiment laboratory (CORE) of the Ulsan National Institute of Science and Technology (UNIST), is a deterministic neutron- and ...photon-transport code primarily designed for light water reactor (LWR) analysis. Initially, the photon module in STREAM did not account for fluorescence and bremsstrahlung photons. This article presents recent developments regarding the integration of atomic relaxation and bremsstrahlung models into the existing photon module, thus allowing for the transport of secondary photons. The photon flux and photon heating computed with the newly incorporated models is compared to results obtained with the Monte Carlo code MCS. The incorporation of secondary photons has substantially improved the accuracy of photon flux calculations, particularly in scenarios involving strong gamma emitters. However, it is essential to note that despite the consideration of secondary photon sources, there is no noticeable improvement in the photon heating for LWR problems when compared to the photon heating obtained with the previous version of STREAM.
A general solution strategy of the modified power iteration method for calculating higher eigenmodes has been developed and applied in continuous energy Monte Carlo simulation. The new approach ...adopts four features: 1) the eigen decomposition of transfer matrix, 2) weight cancellation for higher modes, 3) population control with higher mode weights, and 4) stabilization technique of statistical fluctuations using multi-cycle accumulations. The numerical tests of neutron transport eigenvalue problems successfully demonstrate that the new strategy can significantly accelerate the fission source convergence with stable convergence behavior while obtaining multiple higher eigenmodes at the same time. The advantages of the new strategy can be summarized as 1) the replacement of the cumbersome solution step of high order polynomial equations required by Booth's original method with the simple matrix eigen decomposition, 2) faster fission source convergence in inactive cycles, 3) more stable behaviors in both inactive and active cycles, and 4) smaller variances in active cycles. Advantages 3 and 4 can be attributed to the lower sensitivity of the new strategy to statistical fluctuations due to the multi-cycle accumulations. The application of the modified power method to continuous energy Monte Carlo simulation and the higher eigenmodes up to 4th order are reported for the first time in this paper.
•Modified power method is applied to continuous energy Monte Carlo simulation.•Transfer matrix is introduced to generalize the modified power method.•All mode based population control is applied to get the higher eigenmodes.•Statistic fluctuation can be greatly reduced using accumulated tally results.•Fission source convergence is accelerated with higher mode solutions.
This article presents the convergence analysis of the coarse mesh finite difference (CMFD) method applied to two-group (2-G) three-dimensional (3D) neutron diffusion problem. Two CMFD algorithms are ...examined: one-node (1-N) CMFD and two-node (2-N) CMFD. Two test problems are used for the study of the convergence behavior: a model problem of homogeneous 2-G 3D eigenvalue problem and the NEACRP LWR transient benchmark problem. The convergence rates of the 1-N and 2-N CMFD algorithms are numerically measured in terms of the convergence of current correction factors (CCFs). The numerical test results are presented as well as the comparison with the previous analytical study. Overall, 1-N CMFD with the CCF relaxation shows a comparable performance to 2-N CMFD for the realistic 3D rod ejection transients.
In this study, we incorporate an anisotropic scattering scheme involving spherical harmonics into the method of characteristics (MOC). The neutron transport solution in a light water reactor can be ...significantly improved because of the impact of an anisotropic scattering source with the MOC flat source approximation. Several problems are selected to verify the proposed scheme and investigate its effects and accuracy. The MOC anisotropic scattering source is based on the expansion of spherical harmonics with Legendre polynomial functions. The angular flux, scattering source, and cross section are expanded in terms of the surface spherical harmonics. Later, the polynomial is expanded to achieve the odd and even parity of the source components. Ultimately, the MOC angular and scalar fluxes are calculated from a combination of two sources. This paper presents various numerical examples that represent the hot and cold conditions of a reactor core with boron concentration, burnable absorbers, and control rod materials, with and without a reflector or baffle. Moreover, a small critical core problem is considered which involves significant neutron leakage at room temperature. We demonstrate that an anisotropic scattering source significantly improves solution accuracy for the small core high-leakage problem, as well as for practical large core analyses.
CEFR is a small core-size sodium-cooled fast reactor (SFR) using high enrichment fuel with stainless-steel reflectors, which brings a significant challenge to the deterministic methodologies due to ...the strong spectral effect. The neutronic simulation of the start-up experiments conducted at the CEFR have been performed with a deterministic code system RAST-F, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver. The RAST-F results were compared against the measurement data. Moreover, the characteristic of neutron spectrum in the fuel rings, and adjacent reflectors was evaluated using different models for generation of accurate nuclear libraries. The numerical solution of RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of keff and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental values within two-sigma of experimental uncertainty. The good agreement of these results indicating that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code system.
This study incorporates a high-fidelity transient analysis solver based on multigroup CMFD in the MOC code STREAM. Transport modeling with heterogeneous geometries of the reactor core increases ...computational cost in terms of memory and time, whereas the multigroup CMFD reduces the computational cost. The reactor condition does not change at every time step, which is a vital point for the utilization of CMFD. CMFD correction factors are updated from the transport solution whenever the reactor core condition changes, and the simulation continues until the end. The transport solution is adjusted once CMFD achieves the solution. The flux-weighted method is used for rod decusping to update the partially inserted control rod cell material, which maintains the solution's stability. A smaller time-step size is needed to obtain an accurate solution, which increases the computational cost. The adaptive step-size control algorithm is robust for controlling the time step size. This algorithm is based on local errors and has the potential capability to accept or reject the solution. Several numerical problems are selected to analyze the performance and numerical accuracy of parallel computing, rod decusping, and adaptive time step control. Lastly, a typical pressurized LWR was chosen to study the rod-ejection accident.