• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.21 no.3
• /
• pp.303-313
• /
• 2023

Spent fuels (SFs) are stored in a storage pool after discharge from nuclear power plants. They can be transferred to for the further processes such as dry storage sites, processing plants, or disposal sites. One of important measures of SF is the burnup. Since the radioactivity of SF is strongly dependent on its burnup, the burnup of SF should be well estimated for the safe management, storage, and final disposal. Published papers about the methodology for the burnup estimation from the known activities of important radioactive sources are somewhat rare. In this study, we analyzed the dependency of the burnup on the important radiation source activities using ORIGEN-ARP, and suggested simple correlations that relate the burnup and the important source activities directly. A burnup estimation equation is suggested for PWR fuels relating burnup with total neutron source intensity (TNSI), initial enrichment, and cooling time. And three burnup estimation equations for major gamma sources, ¹³⁷Cs, ¹³⁴Cs, and ¹⁵⁴Eu are also suggested.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1708-1714
• /
• 2019

The purpose of the critical evaluation of the spent fuel pool (SFP) is to verify that the maximum effective multiplication factor ($K_{eff}$) is less than the critical safety limit at 100% stored condition of the spent fuel with the maximum reactivity. At nuclear power plants, the storage standard of spent fuel, ie, the loading curve, is established to prevent criticality from being generated in SFP. Here, the loading curve refers to a graph showing the minimum discharged burnup versus the initial enrichment of spent fuel. Recently, US NRC proposed the new critical safety assessment guideline (DSS-ISG-2010-01, Revision 0) of PWR SFPs and most of utilities in US is following it. Of course, the licensed criterion of the maximum effective multiplication factor of SFP remains unchanged and it should be less than 0.95 from the 95% probability and the 95% confidence level. However, the new guideline is including the new evaluation methodologies like the application of the axial burnup profile, the validation of depletion and criticality code, and trend analysis. Among the new evaluation methodologies, the most important factor that affects $K_{eff}$ is the axial burnup profile of spent fuel. US NRC recommends to consider the axial burnup profiles presented in NUREG-6801 in criticality analysis. In this paper, criticality effect was evaluated considering three profiles, respectively: i) Axial burnup profiles presented in NUREG-6801. ii) Representative PWR axial burnup profile. iii) Uniform axial burnup profile. As the result, the case applying the axial burnup profiles presented in NUREG-6801 showed the highest $K_{eff}$ among three cases. Therefore, we need to introduce a new methodology because it can be issued if the axial burnup profiles presented in NUREG/CR-6801 are applied to the domestic nuclear power plants without any other consideration.

• Nuclear Engineering and Technology
• /
• v.41 no.8
• /
• pp.1109-1114
• /
• 2009

The evolution of a high burnup structure (HBS) in a light water reactor (LWR) $UO_2$ fuel was simulated using the Potts model. A simulation system for the Potts model was defined as a two-dimensional triangular lattice, for which the stored energy was calculated from both the irradiation damage of the $UO_2$ matrix and the formation of a grain boundary in the newly recrystallized small HBS grains. In the simulation, the evolution probability of the HBS is calculated by the system energy difference between before and after the Monte Carlo simulation step. The simulated local threshold burnup for the HBS formation was 62 MWd/kgU, consistent with the observed threshold burnup range of 60-80 MWd/kgU. The simulation revealed that the HBS was heterogeneously nucleated on the intergranular bubbles in the proximity of the threshold burnup and then additionally on the intragranular bubbles for a burnup above 86 MWd/kgU. In addition, the simulation carried out under a condition of no bubbles indicated that the bubbles played an important role in lowering the threshold burnup for the HBS formation, thereby enabling the HBS to be observed in the burnup range of conventional high burnup fuels.

• Nuclear Engineering and Technology
• /
• v.47 no.5
• /
• pp.596-607
• /
• 2015

In this paper, we describe a computer program, iBEST (inverse Burnup ESTimator), that we developed to accurately estimate the burnup histories of spent nuclear fuels based on sample measurement data. The burnup history parameters include initial uranium enrichment, burnup, cooling time after discharge from reactor, and reactor type. The program uses algebraic equations derived using the simplified burnup chains of major actinides for initial estimations of burnup and uranium enrichment, and it uses the ORIGEN-S code to correct its initial estimations for improved accuracy. In addition, we newly developed a stable bisection method coupled with ORIGEN-S to correct burnup and enrichment values and implemented it in iBEST in order to fully take advantage of the new capabilities of ORIGEN-S for improving accuracy. The iBEST program was tested using several problems for verification and well-known realistic problems with measurement data from spent fuel samples from the Mihama-3 reactor for validation. The test results show that iBEST accurately estimates the burnup history parameters for the test problems and gives an acceptable level of accuracy for the realistic Mihama-3 problems.

• Nuclear Engineering and Technology
• /
• v.47 no.7
• /
• pp.924-933
• /
• 2015

The correlation of the isotopic composition of uranium, plutonium, neodymium, and cesium with the burnup for high burnup pressurized water reactor fuels irradiated in nuclear power reactors has been experimentally investigated. The total burnup was determined by Nd-148 and the fractional $^{235}U$ burnup was determined by U and Pu mass spectrometric methods. The isotopic compositions of U, Pu, Nd, and Cs after their separation from the irradiated fuel samples were measured using thermal ionization mass spectrometry. The contents of these elements in the irradiated fuel were determined through an isotope dilution mass spectrometric method using $^{233}U$, $^{242}Pu$, $^{150}Nd$, and $^{133}Cs$ as spikes. The activity ratios of Cs isotopes in the fuel samples were determined using gamma-ray spectrometry. The content of each element and its isotopic compositions in the irradiated fuel were expressed by their correlation with the total and fractional burnup, burnup parameters, and the isotopic compositions of different elements. The results obtained from the experimental methods were compared with those calculated using the ORIGEN-S code.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.901-910
• /
• 2021

KHNP recently has obtained the approval for the commercialization of the modified 37-element (or 37 M) fuel bundle and now is loading the 37 M fuel bundles in CANDU-6 reactors in KOREA. One of the main issues for approval was the burnup limit. Due to CANDU design and neutronic characteristics, there was no specific burnup restriction of a fuel bundle. The absence of a burnup limit does not mean that a fuel bundle can stay in the CANDU reactor without a time limit. However, the regulator requested traditional design values as well as the burnup limit reflecting the computer code uncertainty. The method for the PWR burnup limit was not applied to the CANDU fuel bundle. Since there was no approved methodology to build the burnup limit with uncertainties, KHNP introduced a Monte-Carlo method coupled with a 95/95 approach to determine the conservative burnup limit from the viewpoint of the centerline temperature, internal pressure, strain measurement deviation. Moreover, to consider the uncertainties of various computing models, a converted power uncertainty was introduced. This paper presents the methodology and puts forward the limits on burnup, evaluated for each of the existing and modified fuel bundles, in consideration of the pressure tube aging condition.

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

• Nuclear Engineering and Technology
• /
• v.48 no.3
• /
• pp.624-634
• /
• 2016

Nuclear criticality safety analyses (NCSAs) considering burnup credit were performed for the GBC-32 cask. The used nuclear fuel assemblies (UNFAs) discharged from Hanbit Nuclear Power Plant Unit 3 Cycle 6 were loaded into the cask. Their axial burnup distributions and average discharge burnups were evaluated using the DeCART and Multi-purpose Analyzer for Static and Transient Effects of Reactors (MASTER) codes, and NCSAs were performed using SCALE 6.1/STandardized Analysis of Reactivity for Burnup Credit using SCALE (STARBUCS) and Monte Carlo N-Particle transport code, version 6 (MCNP 6). The axial burnup distributions were determined for 20 UNFAs with various initial enrichments and burnups, which were applied to the criticality analysis for the cask system. The UNFAs for 20- and 30-year cooling times were assumed to be stored in the cask. The criticality analyses indicated that $k_{eff}$ values for UNFAs with nonuniform axial burnup distributions were larger than those with a uniform distribution, that is, the end effects were positive but much smaller than those with the reference distribution. The axial burnup distributions for 20 UNFAs had shapes that were more symmetrical with a less steep gradient in the upper region than the reference ones of the United States Department of Energy. These differences in the axial burnup distributions resulted in a significant reduction in end effects compared with the reference.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05b
• /
• pp.279-284
• /
• 1998

Applicability of the threshold burnup for rim formation was investigated as a function of temperature by Rest's model. The threshold burnup was the lowest in the intermediate temperature region, while on the other temperature regions the threshold burnup is higher. The rim porosity was predicted by the van der Waals equation based of the rim pore radius of 0.75${\mu}{\textrm}{m}$ and the overpressurization model on rim pores. The calculated centerline temperature is in good agreement with the measured temperature. However, more efforts seem to be necessary for the mechanistic model of the rim effect including rim growth with the fuel burnup.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05c
• /
• pp.265-270
• /
• 1996

The temperature distribution in the pellet was obtained from beginning the general heat conduction equation. The thermal conductivity of pellet used the SIMFUEL data that made clear the effect of burnup on the thermal conductivity degradation. Since the pellet rim acts as the thermal barrier to heat flow. the pellet was subdivided into several rings in which the outer ring was adjusted to play almost the same role as the rim. The local burup in each ring except the outer ring was calculated from the power depression factor based on FASER results. whereas the rim burnup at the outer ring was achieved by the pellet averaged burnup based on the empirical relation. The rim changed to the equivalent Xe film so the predicted temperature shooed the thermal jump across the rim. The observed temperature profiles depended on linear heat generation rate. fuel burnup. and power depression factor. The thermal conductivity degradation modelling can be applied to the fuel performance code to high burnup fuel,