Validation of UNIST Monte Carlo Code MCS using VERA Progression Problems

This paper presents the validation of UNIST in-house Monte Carlo code MCS used for the high-fidelity simulation of commercial pressurized water reactors (PWRs). Its focus is on the accurate, spatially detailed neutronic analyses of startup physics tests for the initial core of the Watts Bar Nuclear 1 reactor, which is a vital step in evaluating core phenomena in an operating nuclear power reactor. The MCS solutions for the Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) core physics benchmark progression problems 1 to 5 were verified with KENO-VI and Serpent 2 solutions for geometries ranging from a single-pin cell to a full core. MCS was also validated by comparing with results of reactor zero-power physics tests in a full-core simulation. MCS exhibits an excellent consistency against the measured data with a bias of +/- 3 pcm at the initial criticality whole-core problem. Furthermore, MCS solutions for rod worth are consistent with measured data, and reasonable agreement is obtained for the isothermal temperature coefficient and soluble boron worth. This favorable comparison with measured parameters exhibited by MCS continues to broaden its validation basis. These results provide confidence in MCS's capability in high-fidelity calculations for practical PWR cores. (C) 2019 Korean Nuclear Society, Published by Elsevier Korea LLC

Physics Study of Canada Deuterium Uranium Lattice with Coolant Void Reactivity Analysis

This study presents a coolant void reactivity analysis of Canada Deuterium Uranium (CANDU)-6 and Advanced Canada Deuterium Uranium Reactor-700 (ACR-700) fuel lattices using a Monte Carlo code. The reactivity changes when the coolant was voided were assessed in terms of the contributions of four factors and spectrum shifts. In the case of single bundle coolant voiding, the contribution of each of the four factors in the ACR-700 lattice is large in magnitude with opposite signs, and their summation becomes a negative reactivity effect in contrast to that of the CANDU-6 lattice. Unlike the coolant voiding in a single fuel bundle, the 2 ?? 2 checkerboard coolant voiding in the ACR-700 lattice shows a positive reactivity effect. The neutron current between the no-void and voided bundles, and the four factors of each bundle were analyzed to figure out the mechanism of the positive coolant void reactivity of the checkerboard voiding case. Through a sensitivity study of fuel enrichment, type of burnable absorber, and moderator to fuel volume ratio, a design strategy for the CANDU reactor was suggested in order to achieve a negative coolant void reactivity even for the checkerboard voiding case.ope

Comparative Neutronics Analysis of DIMPLE S06 Criticality Benchmark with Contemporary Reactor Core Analysis Computer Code Systems

A high-leakage core has been known to be a challenging problem not only for a two-step homogenization approach but also for a direct heterogeneous approach. In this paper the DIMPLE S06 core, which is a small high-leakage core, has been analyzed by a direct heterogeneous modeling approach and by a two-step homogenization modeling approach, using contemporary code systems developed for reactor core analysis. The focus of this work is a comprehensive comparative analysis of the conventional approaches and codes with a small core design, DIMPLE S06 critical experiment. The calculation procedure for the two approaches is explicitly presented in this paper. Comprehensive comparative analysis is performed by neutronics parameters: multiplication factor and assembly power distribution. Comparison of two-group homogenized cross sections from each lattice physics codes shows that the generated transport cross section has significant difference according to the transport approximation to treat anisotropic scattering effect. The necessity of the ADF to correct the discontinuity at the assembly interfaces is clearly presented by the flux distributions and the result of two-step approach. Finally, the two approaches show consistent results for all codes, while the comparison with the reference generated by MCNP shows significant error except for another Monte Carlo code, SERPENT2open0

Incorporation of anisotropic scattering into the method of characteristics

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 signif-icantly 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.(c) 2022 Korean Nuclear Society, Published by Elsevier Korea LLC. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

Neutronic simulation of the CEFR experiments with the nodal diffusion code system RAST-F

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 l e ft 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. (C) 2022 Korean Nuclear Society, Published by Elsevier Korea LLC

Monte Carlo Analysis of the Accelerator-Driven System at Kyoto University Research Reactor Institute

An accelerator-driven system consists of a subcritical reactor and a controllable external neutron source. The reactor in an accelerator-driven system can sustain fission reactions in a subcritical state using an external neutron source, which is an intrinsic safety feature of the system. The system can provide efficient transmutations of nuclear wastes such as minor actinides and long-lived fission products and generate electricity. Recently at Kyoto University Research Reactor Institute (KURRI; Kyoto, Japan), a series of reactor physics experiments was conducted with the Kyoto University Critical Assembly and a Cockcroft-Walton type accelerator, which generates the external neutron source by deuterium-tritium reactions. In this paper, neutronic analyses of a series of experiments have been re-estimated by using the latest Monte Carlo code and nuclear data libraries. This feasibility study is presented through the comparison of Monte Carlo simulation results with measurements.clos

Development of nodal diffusion code RAST-V for Vodo-Vodyanoi Energetichesky reactor analysis

This paper presents the development of a nodal diffusion code, RAST-V, and its verification and validation for VVER (vodo-vodyanoi energetichesky reactor) analysis. A VVER analytic solver has been imple-mented in an in-house nodal diffusion code, RAST-K. The new RAST-K version, RAST-V, uses the triangle -based polynomial expansion nodal method. The RAST-K code provides stand-alone and two-step computation modes for steady-state and transient calculations. An in-house lattice code (STREAM) with updated features for VVER analysis is also utilized in the two-step method for cross-section gen-eration. To assess the calculation capability of the formulated analysis module, various verification and validation studies have been performed with Rostov-II, and X2 multicycles, Novovoronezh-4, and the Atomic Energy Research benchmarks. In comparing the multicycle operation, rod worth, and integrated temperature coefficients, RAST-V is found to agree with measurements with high accuracy which RMS differences of each cycle are within +/- 47 ppm in multicycle operations, and +/- 81 pcm of the rod worth of the X2 reactor. Transient calculations were also performed considering two different rod ejection sce-narios. The accuracy of RAST-V was observed to be comparable to that of conventional nodal diffusion codes (DYN3D, BIPR8, and PARCS).(c) 2022 Korean Nuclear Society, Published by Elsevier Korea LLC. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

Validation of nuclide depletion capabilities in Monte Carlo code MCS

In this work, the depletion capability implemented in Monte Carlo code MCS is investigated to predict the isotopic compositions of spent nuclear fuel (SNF). By comparison of MCS calculation results to post irradiation examination (PIE) data obtained from one pressurized water reactor (PWR), the validation of this capability is conducted. The depletion analysis is performed with the ENDF/B-VII.1 library and a fuel assembly model. The transmutation equation is solved by the Chebyshev Rational Approximation Method (CRAM) with a depletion chain of 3820 isotopes. 18 actinides and 19 fission products are analyzed in 14 SNF samples. The effect of statistical uncertainties on the calculated number densities is discussed. On average, most of the actinides and fission products analyzed are predicted within +/- 6% of the experiment. MCS depletion results are also compared to other depletion codes based on publicly reported information in literature. The code-to-code analysis shows comparable accuracy. Overall, it is demonstrated that the depletion capability in MCS can be reliably applied in the prediction of SNF isotopic inventory. (c) 2020 Korean Nuclear Society, Published by Elsevier Korea LLC. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

Use of Monte Carlo code MCS for multigroup cross section generation for fast reactor analysis

Multigroup cross section (MG XS) generation by the UNIST in-house Monte Carlo (MC) code MCS for fast reactor analysis using nodal diffusion codes is reported. The feasibility of the approach is quantified for two sodium fast reactors (SFRs) specified in the OECD/NEA SFR benchmark: a 1000 MWth metal-fueled SFR (MET-1000) and a 3600 MWth oxide-fueled SFR (MOX-360 0). The accuracy of a few-group XSs generated by MCS is verified using another MC code, Serpent 2. The neutronic steady-state whole-core problem is analyzed using MCS/RAST-K with a 24-group XS set. Various core parameters of interest (core keff, power profiles, and reactivity feedback coefficients) are obtained using both MCS/RAST-K and MCS. A code-to-code comparison indicates excellent agreement between the nodal diffusion solution and sto-chastic solution; the error in the core keff is less than 110 pcm, the root-mean-square error of the power profiles is within 1.0%, and the error of the reactivity feedback coefficients is within three standard deviations. Furthermore, using the super-homogenization-corrected XSs improves the prediction accu-racy of the control rod worth and power profiles with all rods in. Therefore, the results demonstrate that employing the MCS MG XSs for the nodal diffusion code is feasible for high-fidelity analyses of fast reactors. (C) 2021 Korean Nuclear Society, Published by Elsevier Korea LLC

Interpretation of two SINBAD photon-leakage benchmarks with nuclear library ENDF/B-VIII.0 and Monte Carlo code MCS

A review of the documentation and an interpretation of the NEA-1517/74 and NEA-1517/80 shielding benchmarks (measurements of photon leakage flux from a hollow sphere with a central 14 MeV neutron source) from the SINBAD database with the Monte Carlo code MCS and the most up-to-date ENDF/B-VIII.0 neutron data library are conducted. The two analyzed benchmarks describe satisfactorily the energy resolution of the photon detector and the geometry of the spherical samples with inner beam tube, tritium target and cooling water circuit, but lack information regarding the detector geometry and the distances of shields and collimators relatively to the neutron source and the detector. Calculations are therefore conducted for a sphere model only. A preliminary verification of MCS neutron-photon calculations against MCNP6.2 is first conducted, then the impact of modelling the inner beam tube, tritium target and cooling water circuit is assessed. Finally, a comparison of calculated results with the libraries ENDF/B-VII.1 and ENDF/B-VIII.0 against the measurements is conducted and shows reasonable agreement. The MCS and MCNP inputs used for the interpretation are available as supplementary material of this article. (C) 2019 Korean Nuclear Society, Published by Elsevier Korea LLC