• Energy Research

Abstract This paper presents a performance examination of the Coarse Mesh Finite Difference (CMFD) method for a reduction of inter-cycle correlations in Monte Carlo (MC) active cycle simulations. Sensitivity tests using the UNIST in-house MC code ‘MCS’ reveal the characteristics of CMFD effectiveness depending on problem sizes, dimensions , and cross-section dependency on neutron energy (multi-group vs. continuous). Also, by applying this method to the three-dimensional BEAVRS whole core problem, it is successfully demonstrated that the CMFD can reduce the inter-cycle correlations of practical, large size, and high dominance ratio problems, and that the figure of merit of efficiency for a practical reactor problem increases by a factor of six.

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. Keywords: Monte Carlo, Core simulation, Zero power physics test, Validation, MCS

The theoretical aspects behind the reactor depletion capability of the Monte Carlo code MCS developed at the Ulsan National Institute of Science and Technology (UNIST) and practical results of this depletion feature for a Material-Testing Reactor (MTR) with plate-type fuel are described in this paper. A verification of MCS results is first performed against MCNP6 to confirm the suitability of MCS for the criticality and depletion analysis of the MTR. Then, the dependence of the effective neutron multiplication factor to the number of axial and radial depletion cells adopted in the fuel plates is performed with MCS in order to determine the minimum spatial segmentation of the fuel plates. Monte Carlo depletion results with 37,800 depletion cells are provided by MCS within acceptable calculation time and memory usage. The results show that at least 7 axial meshes per fuel plate are required to reach the same precision as the reference calculation whereas no significant differences are observed when modeling 1 or 10 radial meshes per fuel plate. This study demonstrates that MCS can address the need for Monte Carlo codes capable of providing reference solutions to complex reactor depletion problems with refined meshes for fuel management and research reactor applications.

Abstract This paper presents the verification and validation of the radiation source term calculation capability implemented in the reactor analysis code STREAM for pressurized water reactor (PWR) spent nuclear fuel (SNF) analysis. Activity, decay heat, neutron and gamma source spectra of irradiated PWR fuel assemblies are calculated with STREAM and ORIGEN and the results obtained are compared. The verification is performed with Westinghouse 15 × 15, 17 × 17 and Combustion Engineering 16 × 16 designs with enrichment of 3–5 wt.% 235U, burnup of 30–60 GWd/t and cooling times of 0.3–1000 years. The radiation source term comparison between STREAM and ORIGEN indicates that the implementation of the source term equations is correct. The validation is done by comparing STREAM calculation results against 91 decay heat calorimetric measurements for 52 PWR fuel assemblies performed at the Swedish central storage facility for spent nuclear fuel, CLAB, and United States Hanford Engineering Development Laboratory (HEDL) and General Electric Morris Operations. The measured fuel assemblies have enrichment of 2–4 wt.% 235U, burnup of 28–51 GWd/t and cooling times of 2–27 years. The comparison between calculated and measured decay heat shows good agreement. For most of the measurements analyzed, the calculated decay heat is within the range of the measurement uncertainty. The average, overall C/E over 71 CLAB decay heat measurements is 1.000 ± 0.017, and 0.997 ± 0.013 for 20 US measurements. Overall, this work demonstrates that STREAM can be reliably used to predict PWR SNF radiation source terms.

Abstract This work investigates the depletion capability implemented in lattice physics code STREAM for the prediction of pressurized water reactor (PWR) uranium dioxide (UO2) spent nuclear fuel (SNF) isotopic inventory. The validation of this capability is performed by comparison of STREAM calculation results to measured SNF assay data obtained from PWRs Takahama-3, Calvert Cliffs and GKN II. The depletion analysis is conducted with the ENDF/B-VII.0 library and uses a pin cell model of the fuel rods from which the fuel samples were taken. The Chebyshev Rational Approximation Method (CRAM) is used to solve the depletion equation with about 1300–1600 isotopes in the depletion chain. 16 actinides and 23 fission products are analyzed in 14 spent UO2 fuel samples. The actinides are isotopes of uranium, neptunium, plutonium, americium and curium. The fission products nuclides include isotopes of cesium, neodymium, europium, samarium as well as 106Ru, 144Ce, 155Gd, 99Tc, 90Sr, 109Ag, and 103Rh. The sensitivity of some of the nuclides to the details of the power history and the adjustment of the fuel sample burnup is discussed. The impact of using ENDF/B-VII.0 library instead of ENDF/B-VI.8 is also discussed. Most of the nuclides analyzed are well predicted within ±7% of the experiment for actinides and fission products. STREAM depletion results are also compared to the codes SWAT, HELIOS and SCALE results based on publicly available information in literature, to check the performance of STREAM relative to other codes for the prediction of SNF isotopic inventory. The comparison to other code systems shows that the implementation in STREAM is of comparable accuracy. Overall, this paper demonstrates that the depletion capability in STREAM can be reliably applied to predict the isotopic inventory of PWR UO2 SNF for burnup ranging from 14 to 54 GWd/t and initial enrichment ranging from 3.0 to 4.1 wt% 235U.

Abstract A depletion chain simplification method is applied to UNIST lattice code STREAM (Steady state and Transient REactor Analysis code with Method of Characteristics) in this paper to alleviate the computational burden of depletion calculation associated with a large depletion matrix. A simplified burnup matrix (burnup chain) of 464 nuclides and 10,638 transitions is thus created from a large detailed burnup matrix containing 3,837 nuclides and 43,416 transitions from ENDF/B-VII.0 nuclear decay data. The simplified burnup matrix is optimized for the purpose of calculating effective neutron multiplication factors (keff) and power distributions with reduced computation time and memory usage. The nuclide selection method relies on the Generalized Perturbation Theory (GPT): by exploiting the adjoint function of nuclide number densities as derived with GPT, a set of nuclides which must be included in the simplified chain (nuclides with important contribution to reactivity) is determined. Numerical verification of the simplified burnup chain is conducted using the deterministic neutron transport analysis code STREAM, developed to perform whole light water reactor (LWR) core calculations with the direct transport analysis method and the two-step method. The simplified burnup chain is tested on the 16 LWR fuel assembly depletion problems from the Virtual Environment for Reactor Application (VERA) benchmarks, including burnable poison (BP) pin cells commonly used in nuclear reactor design such as gadolinia, Pyrex, AIC, B4C and IFBA, and one OPR-1000 fuel assembly with gadolinium bearing fuel depletion problem. For burnup up to 80 MWd/kg and 235U enrichment ranging from 2.1 w/o to 4.6 w/o, the calculations with the simplified burnup chain predict the keff variations within 10 pcm and identical power distributions, with a speed-up factor from 20 to 40 in depletion calculation and a reduction factor by 3–4 of the total simulation time compared to the ones with the detailed burnup chain.

Abstract This paper extends the methodology for treating the resonance self-shielding effect of doubly heterogeneous (DH) tristructural-isotropic fuel particles to take the coating layers of fuel grains into account. Due to the two-region assumption of fuel kernel and graphite matrix in the original derivation of Williams’ method, which was developed to replace the time consuming SCALE continuous energy spectrum calculations with more efficient equivalence theory and intermediate resonance approximation, the method can be applied only to the fuel kernels without coating layers. In order to overcome this limitation, Ji’s analytical expression of the Dancoff factor for multi-layer spheres has been incorporated into the DH method and the applicability has been successfully extended to the kernels with coating layers. It is demonstrated that the new method can predict the DH effects within 200 pcm order of errors against the reference Monte Carlo calculations for wide range of block type gas cooled reactor designs with coating layers.

This paper presents the verification and validation (V&V) of the STREAM/RAST-K 2.0 code system for a pressurized water reactor (PWR) analysis. A lattice physics code STREAM and a nodal diffusion code RAST-K 2.0 have been developed by a computational reactor physics and experiment laboratory (CORE) of Ulsan National Institute of Science and Technology (UNIST) for an accurate two-step PWR analysis. The calculation modules of each code were already verified against various benchmark problems, whereas this paper focuses on the V&V of linked code system. Three PWR type reactor cores, OPR-1000, three-loop Westinghouse reactor core, and APR-1400, are selected as V&V target plants. This code system, for verification, is compared against the conventional code systems used for the calculations in nuclear design reports (NDRs) and validated against measured plant data. Compared parameters are as follows: critical boron concentration (CBC), axial shape index (ASI), assembly-wise power distribution, burnup distribution and peaking factors. STREAM/RAST-K 2.0 shows the RMS error of critical boron concentration within 20 ppm, and the RMS error of assembly power within 1.34% for all the cycles of all reactors. Keywords: Verification and validation, PWR core, Two-step approach, STREAM, RAST-K 2.0